Patterns of strain in the macaque tibia during functional activity

Mechanical and morphometric approaches to body mass estimation in rhesus macaques: A test of skeletal variables

OBJECTIVES

Estimation of body mass from skeletal metrics can reveal important insights into the paleobiology of archeological or fossil remains. The standard approach constructs predictive equations from postcrania, but studies have questioned the reliability of traditional measures. Here, we examine several skeletal features to assess their accuracy in predicting body mass.

MATERIALS AND METHODS

Antemortem mass measurements were compared with common skeletal dimensions from the same animals postmortem, using 115 rhesus macaques (male: n = 43; female: n = 72). Individuals were divided into training (n = 58) and test samples (n = 57) to build and assess Ordinary Least Squares or multivariate regressions by residual sum of squares (RSS) and AIC weights. A leave-one-out approach was implemented to formulate the best fit multivariate models, which were compared against a univariate and a previously published catarrhine body-mass estimation model.

RESULTS

Femur circumference represented the best univariate model. The best model overall was composed of four variables (femur, tibia and fibula circumference and humerus length). By RSS and AICw, models built from rhesus macaque data (RSS = 26.91, AIC = -20.66) better predicted body mass than did the catarrhine model (RSS = 65.47, AIC = 20.24).

CONCLUSION

Body mass in rhesus macaques is best predicted by a 4-variable equation composed of humerus length and hind limb midshaft circumferences. Comparison of models built from the macaque versus the catarrhine data highlight the importance of taxonomic specificity in predicting body mass. This paper provides a valuable dataset of combined somatic and skeletal data in a primate, which can be used to build body mass equations for fragmentary fossil evidence.

Cortical and Trabecular Bone Modeling and Implications for Bone Functional Adaptation in the Mammalian Tibia

Bone modeling involves the addition of bone material through osteoblast-mediated deposition or the removal of bone material via osteoclast-mediated resorption in response to perceived changes in loads by osteocytes. This process is characterized by the independent occurrence of deposition and resorption, which can take place simultaneously at different locations within the bone due to variations in stress levels across its different regions. The principle of bone functional adaptation states that cortical and trabecular bone tissues will respond to mechanical stimuli by adjusting (i.e., bone modeling) their morphology and architecture to mechanically improve their mechanical function in line with the habitual in vivo loading direction. This principle is relevant to various research areas, such as the development of improved orthopedic implants, preventative medicine for osteopenic elderly patients, and the investigation of locomotion behavior in extinct species. In the present review, the mammalian tibia is used as an example to explore cortical and trabecular bone modeling and to examine its implications for the functional adaptation of bones. Following a short introduction and an exposition on characteristics of mechanical stimuli that influence bone modeling, a detailed critical appraisal of the literature on cortical and trabecular bone modeling and bone functional adaptation is given. By synthesizing key findings from studies involving small mammals (rodents), large mammals, and humans, it is shown that examining both cortical and trabecular bone structures is essential for understanding bone functional adaptation. A combined approach can provide a more comprehensive understanding of this significant physiological phenomenon, as each structure contributes uniquely to the phenomenon.

Changes in limb bone diaphyseal structure in chimpanzees during development

OBJECTIVES

This study tests if femoral and humeral cross-sectional geometry (CSG) and cross-sectional properties (CSPs) in an ontogenetic series of wild-caught chimpanzees (Pan troglodytes ssp.) reflect locomotor behavior during development. The goal is to clarify the relationship between limb bone structure and locomotor behavior during ontogeny in Pan.

MATERIALS AND METHODS

The latex cast method was used to reconstruct cross sections at the midshaft femur and mid-distal humerus. Second moments of area (SMAs) (Ix, Iy, Imax, Imin), which are proportional to bending rigidity about a specified axis, and the polar SMA (J), which is proportional to average bending rigidity, were calculated at section locations. Cross-sectional shape (CSS) was assessed from Ix/Iy and Imax/Imin ratios. Juvenile and adult subsamples were compared.

RESULTS

Juveniles and adults have significantly greater femoral J compared to humeral J. Mean interlimb proportions of J are not significantly different between the groups. There is an overall decreasing trend in diaphyseal circularity between the juvenile phase of development and adulthood, although significant differences are only found in the humerus.

DISCUSSION

Juvenile chimpanzee locomotion includes forelimb- and hindlimb-biased behaviors. Juveniles and adults preferentially load their hindlimbs relative to their forelimbs. This may indicate similar locomotor behavior, although other explanations including a diversity of hindlimb-biased locomotor behaviors in juveniles cannot be ruled out. Different ontogenetic trends in forelimb and hindlimb CSS are consistent with limb bone CSG reflecting functional adaptation, albeit the complex nature of bone functional adaptation requires cautious interpretations of skeletal functional morphology from biomechanical analyses.

Strain-mode-specific mechanical testing and the interpretation of bone adaptation in the deer calcaneus

The artiodactyl (deer and sheep) calcaneus is a model that helps in understanding how many bones achieve anatomical optimization and functional adaptation. We consider how the dorsal and plantar cortices of these bones are optimized in quasi-isolation (the conventional view) versus in the context of load sharing along the calcaneal shaft by "tension members" (the plantar ligament and superficial digital flexor tendon). This load-sharing concept replaces the conventional view, as we have argued in a recent publication that employs an advanced analytical model of habitual loading and fracture risk factors of the deer calcaneus. Like deer and sheep calcanei, many mammalian limb bones also experience prevalent bending, which seems problematic because the bone is weaker and less fatigue-resistant in tension than compression. To understand how bones adapt to bending loads and counteract deleterious consequences of tension, it is important to examine both strain-mode-specific (S-M-S) testing (compression testing of bone habitually loaded in compression; tension testing of bone habitually loaded in tension) and non-S-M-S testing. Mechanical testing was performed on individually machined specimens from the dorsal "compression cortex" and plantar "tension cortex" of adult deer calcanei and were independently tested to failure in one of these two strain modes. We hypothesized that the mechanical properties of each cortex region would be optimized for its habitual strain mode when these regions are considered independently. Consistent with this hypothesis, energy absorption parameters were approximately three times greater in S-M-S compression testing in the dorsal/compression cortex when compared to non-S-M-S tension testing of the dorsal cortex. However, inconsistent with this hypothesis, S-M-S tension testing of the plantar/tension cortex did not show greater energy absorption compared to non-S-M-S compression testing of the plantar cortex. When compared to the dorsal cortex, the plantar cortex only had a higher elastic modulus (in S-M-S testing of both regions). Therefore, the greater strength and capacity for energy absorption of the dorsal cortex might "protect" the weaker plantar cortex during functional loading. However, this conventional interpretation (i.e., considering adaptation of each cortex in isolation) is rejected when critically considering the load-sharing influences of the ligament and tendon that course along the plantar cortex. This important finding/interpretation has general implications for a better understanding of how other similarly loaded bones achieve anatomical optimization and functional adaptation.

Breed and loading history influence in vivo skeletal strain patterns in pre-pubertal female chickens

The influence of loading history on in vivo strains within a given specie remains poorly understood, and although in vivo strains have been measured at the hindlimb bones of various species, strains engendered during modes of activity other than locomotion are lacking, particularly in non-human species. For commercial egg-laying chickens specifically, there is an interest in understanding their bones' mechanical behaviour, particularly during youth, to develop early interventions to prevent the high incidence of osteoporosis in this population. We measured in vivo mechanical strains at the tibiotarsus midshaft during steady activities (ground, uphill, downhill locomotion) and non-steady activities (perching, jumping, aerial transition landing) in 48 pre-pubescent female (egg-laying) chickens from two breeds that were reared in three different housing systems, allowing varying amounts and types of physical activity. Mechanical strain patterns differed between breeds, and were dependent on the activity performed. Mechanical strains were also affected by rearing environment: chickens that were restricted from performing dynamic load bearing activity due to caged-housing generally exhibited higher mechanical strain levels during steady, but not non-steady activities, compared to chickens with prior dynamic load-bearing activity experience. Among chickens with prior experience of dynamic load bearing activity, those reared in housing systems that allowed more frequent physical activity did not exhibit lower mechanical strains. In all groups, the tibiotarsus was subjected to a loading environment consisting of a combination of axial compression, bending, and torsion, with torsion being the predominant source of strain. Aerial transition landing produced the highest strain levels with unusual strain patterns compared to other activities, suggesting it may produce the strongest anabolic response. These results exemplify how different breeds within a given specie adapt to maintain different patterns of mechanical strains, and how benefits of physical activity in terms of resistance to strain are activity-type dependent and do not necessarily increase with increased physical activity. These findings directly inform controlled loading experiments aimed at studying the bone mechanoresponse in young female chickens and can also be associated to measures of bone morphology and material properties to understand how these features influence bone mechanical properties in vivo.

Knuckle-walking in Sahelanthropus? Locomotor inferences from the ulnae of fossil hominins and other hominoids

Because the ulna supports and transmits forces during movement, its morphology can signal aspects of functional adaptation. To test whether, like extant apes, some hominins habitually recruit the forelimb in locomotion, we separate the ulna shaft and ulna proximal complex for independent shape analyses via elliptical Fourier methods to identify functional signals. We examine the relative influence of locomotion, taxonomy, and body mass on ulna contours in Homo sapiens (n = 22), five species of extant apes (n = 33), two Miocene apes (Hispanopithecus and Danuvius), and 17 fossil hominin specimens including Sahelanthropus, Ardipithecus, Australopithecus, Paranthropus, and early Homo. Ulna proximal complex contours correlate with body mass but not locomotor patterns, while ulna shafts significantly correlate with locomotion. African apes' ulna shafts are more robust and curved than Asian apes and are unlike other terrestrial mammals (including other primates), curving ventrally rather than dorsally. Because this distinctive curvature is absent in orangutans and hylobatids, it is likely a function of powerful flexors engaged in wrist and hand stabilization during knuckle-walking, and not an adaptation to climbing or suspensory behavior. The OH 36 (purported Paranthropus boisei) and TM 266 (assigned to Sahelanthropus tchadensis) fossils differ from other hominins by falling within the knuckle-walking morphospace, and thus appear to show forelimb morphology consistent with terrestrial locomotion. Discriminant function analysis classifies both OH 36 and TM 266 with Pan and Gorilla with high posterior probability. Along with its associated femur, the TM 266 ulna shaft contours and its deep, keeled trochlear notch comprise a suite of traits signaling African ape-like quadrupedalism. While implications for the phylogenetic position and hominin status of S. tchadensis remain equivocal, this study supports the growing body of evidence indicating that S. tchadensis was not an obligate biped, but instead represents a late Miocene hominid with knuckle-walking adaptations.

Does a suspensory lifestyle result in increased tensile strength?: Organ level material properties of sloth limb bones

The material composition of vertebrate connective tissue is highly conserved across taxa. Existing data suggest that the compressive and tensile strength of limb bones are very similar despite marked variation in limb posture and locomotor patterns. However, the material properties of limb bone tissue from suspensory taxa have not been formally evaluated. Sloths are nearly obligatory in their use of below-branch suspensory locomotion and posture, thus placing their limb bones and associated soft tissue structures under routine tensile loading. It is possible that sloth limb bones are modified for enhanced tensile strength, perhaps at the expense of compressive strength. Fore- and hindlimb bones of two-toed (Choloepus hoffmanni) and three-toed (Bradypus variegatus) sloths were tested in compression and bending to evaluate this hypothesis. Strength and elastic (Young's) modulus were similarly lower in sloth limb bones during both compression and bending, as compared to pronograde taxa. Ratios of peak bending strength to compressive strength additionally were elevated (sloths: 1.4-1.7; upright taxa: 0.6-1.2) for sloth limb bones. Overall, the material properties measured from the limb bones of tree sloths support our hypothesis of predicted function in a tensile limb system. Future studies should aim to directly test bones in tension to confirm indications of elevated axial tensile strength. Nevertheless, the results herein expand understanding of functional adaptation in mammalian tissue for a range of locomotor/postural behaviors that were previously unexplored.

Efficacy of shear strain gradients as an osteogenic stimulus

The question of which mechanical variables are responsible for inducing osteogenic activity is unresolved despite extensive experimental and theoretical investigation. Candidate variables include strain magnitude, loading frequency, the interaction of magnitude and frequency (strain rate), and strain gradients. An additional challenge is discerning the coordination of periosteal and endosteal expansion during growth, and whether this coordination (or lack thereof) is fully dependent or partially independent of the local mechanical environment. In this study, under the assumption that calculated stresses correspond to relative strain magnitudes, we specify alternative growth algorithms of bone cross-sectional size and geometry to explore skeletal growth under alternative scenarios of osteogenic activity that are tracking 1) an attractor stress, 2) local stress magnitude or 3) steepness of stress gradients. These developmental simulations are initiated from two initial geometries (symmetrical and asymmetrical ellipses) under a time-varying torsional load whose magnitude is proportional to body size growth in a model primate. In addition, we model endosteal expansion under three conditions hypothesized in the literature, in which endosteal expansion is 1) independent of the mechanical milieu, 2) completely dependent on the mechanical milieu, and 3) a "hybrid" model in which intrinsic biological (independent) growth is operative early but gives way to mechanically-sensitive (dependent) growth at later ages. Three variables were recorded over each growth simulation: the safety factor (ratio of yield stress to actual stress), an efficiency ratio (invested bone area per unit of stress), and proximity to an isostress condition (an optimal design criterion in which stress is invariant throughout the structure). The attractor stress algorithm produces the most "adapted" bones in terms of mechanical competence and economy of material. Localized osteogenic activity that is guided in direct proportion to stress magnitude produces competent bones but with variable adult geometries depending on conditions of endosteal expansion. Stress gradients also produce functional but relatively inefficient bones, with widely variable safety factors during growth and heterogeneous stress fields. If, in fact, the osteocyte network monitors strain gradients to generate osteogenic signals, the resulting morphology is competent but falls well short of an optimal mechanical solution.

Ontogenetic and morphological variation in primate long bones reflects signals of size and behavior

OBJECTIVES

Many primates change their locomotor behavior as they mature from infancy to adulthood. Here we investigate how long bone cross-sectional geometry in Pan, Gorilla, Pongo, Hylobatidae, and Macaca varies in shape and form over ontogeny, including whether specific diaphyseal cross sections exhibit signals of periosteal adaptation or canalization.

MATERIALS AND METHODS

Diaphyseal cross sections were analyzed in an ontogenetic series across infant, juvenile, and adult subgroups. Three-dimensional laser-scanned long bone models were sectioned at midshaft (50% of biomechanical length) and distally (20%) along the humerus and femur. Traditional axis ratios acted as indices of cross-sectional circularity, while geometric morphometric techniques were used to study cross-sectional allometry and ontogenetic trajectory.

RESULTS

The humeral midshaft is a strong indicator of posture and locomotor profile in the sample across development, while the mid-femur appears more reflective of shifts in size. By comparison, the distal diaphyses of both limb elements are more ontogenetically constrained, where periosteal shape is largely static across development relative to size, irrespective of a given taxon's behavior or ecology.

DISCUSSION

Primate limb shape is not only highly variable between taxa over development, but at discrete humeral and femoral diaphyseal locations. Overall, periosteal shape of the humeral and femoral midshaft cross sections closely reflects ontogenetic transitions in behavior and size, respectively, while distal shape in both bones appears more genetically constrained across intraspecific development, regardless of posture or size. These findings support prior research on tradeoffs between function and safety along the limbs.

Haversian remodeling corresponds to load frequency but not strain magnitude in the macaque (Macaca fascicularis) skeleton

One way bone adapts to its mechanical environment is by Haversian remodeling, a repair process in which existing bone is resorbed and replaced by new bone. Haversian remodeling forms interconnected, cylindrical structures called secondary osteons. The amount of remodeling that occurs is related to the nature of mechanical loading and accrual of microdamage, but it is uncertain whether habitual loads of high magnitude versus high frequency result in more remodeling. The answer to this question is important if remodeling is to be a tool for inferring loading environments, and thus behavior, in past populations. Here, secondary osteon population density (OPD), osteon cross-sectional area (On.Ar), and percent Haversian bone (%HAV) were compared among mid-diaphysis femora, tibia, fibulae, and mid-level ribs of five adult crab-eating macaques (Macaca fascicularis). Ribs experience relatively low strains but have a high daily loading frequency (~33 times per minute). Limb bones are loaded for fewer cycles per day, but the femur and tibia have high load magnitudes due to gravitational forces. Strain magnitudes in the fibula are a fraction of those in the femur and tibia. Analyses of variance demonstrated significant differences in OPD (P = 0.010) and On.Ar (P < 0.001) among the bones. Pairwise t-tests revealed greater OPD but lower On.Ar in the rib than all other bones. The high rib OPD suggests that Haversian remodeling is more responsive to load frequency than strain magnitude. The fact that osteons are smaller in ribs than any other bone may be an effect of remodeling in comparatively narrow cortices.

Computed tomography shows high fracture prevalence among physically active forager-horticulturalists with high fertility

Modern humans have more fragile skeletons than other hominins, which may result from physical inactivity. Here, we test whether reproductive effort also compromises bone strength, by measuring using computed tomography thoracic vertebral bone mineral density (BMD) and fracture prevalence among physically active Tsimane forager-horticulturalists. Earlier onset of reproduction and shorter interbirth intervals are associated with reduced BMD for women. Tsimane BMD is lower versus Americans, but only for women, contrary to simple predictions relying on inactivity to explain skeletal fragility. Minimal BMD differences exist between Tsimane and American men, suggesting that systemic factors other than fertility (e.g. diet) do not easily explain Tsimane women's lower BMD. Tsimane fracture prevalence is also higher versus Americans. Lower BMD increases Tsimane fracture risk, but only for women, suggesting a role of weak bone in women's fracture etiology. Our results highlight the role of sex-specific mechanisms underlying skeletal fragility that operate long before menopause.

Selection for longer limbs in mice increases bone stiffness and brittleness, but does not alter bending strength

The ability of a bone to withstand loads depends on its structural and material properties. These tend to differ among species with different modes of locomotion, reflecting their unique loading patterns. The evolution of derived limb morphologies, such as the long limbs associated with jumping, may compromise overall bone strength. We evaluated bone mechanical properties in the Longshanks mouse, which was selectively bred for increased tibia length relative to body mass. We combined analyses of 3D shape and cross-sectional geometry of the tibia, with mechanical testing and bone composition assays, to compare bone strength, elastic properties and mineral composition in Longshanks mice and randomly bred controls. Our data show that, despite being more slender, cortical geometry and predicted bending strength of the Longshanks tibia were similar to controls. In whole bone bending tests, measures of bone bending strength were similar across groups; however, Longshanks tibiae were significantly more rigid, more brittle, and required less than half the energy to fracture. Tissue-level elastic properties were also altered in Longshanks mice, but the bones did not differ from the control in water content, ash content or density. These results indicate that while Longshanks bones are as strong as control tibiae, selection for increased tibia length has altered its elastic properties, possibly through changes in organic bony matrix composition. We conclude that selection for certain limb morphologies, and/or selection for rapid skeletal growth, can lead to tissue-level changes that can increase the risk of skeletal fracture, which in turn may favor the correlated evolution of compensatory mechanisms to mitigate increased fracture risk, such as delayed skeletal maturity.

Differences in the Cortical Structure of the Whole Fibula and Tibia Between Long-Distance Runners and Untrained Controls. Toward a Wider Conception of the Biomechanical Regulation of Cortical Bone Structure

The cortical structure of human fibula varies widely throughout the bone suggesting a more selective adaptation to different mechanical environments with respect to the adjacent tibia. To test this hypothesis, serial-pQCT scans of the dominant fibulae and tibiae of 15/15 men/women chronically trained in long-distance running were compared with those of 15/15 untrained controls. When compared to controls, the fibulae of trained individuals had similar (distally) or lower (proximally) cortical area, similar moments of inertia (MI) for anterior-posterior bending (xMI) and lower for lateral bending (yMI) with a lower "shape-index" (yMI/xMI ratio) throughout, and higher resistance to buckling distally. These group differences were more evident in men and independent of group differences in bone mass. These results contrast with those observed in the tibia, where, as expected, structural indicators of bone strength were greater in trained than untrained individuals. Proximally, the larger lateral flexibility of runners' fibulae could improve the ability to store energy, and thereby contribute to fast-running optimization. Distally, the greater lateral fibular flexibility could reduce bending strength. The latter appears to have been compensated by a higher buckling strength. Assuming that these differences could be ascribed to training effects, this suggests that usage-derived strains in some bones may modify their relative structural resistance to different kinds of deformation in different regions, not only regarding strength, but also concerning other physiological roles of the skeleton.

Cancellous bone and theropod dinosaur locomotion. Part I-an examination of cancellous bone architecture in the hindlimb bones of theropods

This paper is the first of a three-part series that investigates the architecture of cancellous ('spongy') bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and has previously been used to infer locomotor biomechanics in extinct tetrapod vertebrates, especially primates. Despite great promise, cancellous bone architecture has remained little utilized for investigating locomotion in many other extinct vertebrate groups, such as dinosaurs. Documentation and quantification of architectural patterns across a whole bone, and across multiple bones, can provide much information on cancellous bone architectural patterns and variation across species. Additionally, this also lends itself to analysis of the musculoskeletal biomechanical factors involved in a direct, mechanistic fashion. On this premise, computed tomographic and image analysis techniques were used to describe and analyse the three-dimensional architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs for the first time. A comprehensive survey across many extant and extinct species is produced, identifying several patterns of similarity and contrast between groups. For instance, more stemward non-avian theropods (e.g. ceratosaurs and tyrannosaurids) exhibit cancellous bone architectures more comparable to that present in humans, whereas species more closely related to birds (e.g. paravians) exhibit architectural patterns bearing greater similarity to those of extant birds. Many of the observed patterns may be linked to particular aspects of locomotor biomechanics, such as the degree of hip or knee flexion during stance and gait. A further important observation is the abundance of markedly oblique trabeculae in the diaphyses of the femur and tibia of birds, which in large species produces spiralling patterns along the endosteal surface. Not only do these observations provide new insight into theropod anatomy and behaviour, they also provide the foundation for mechanistic testing of locomotor hypotheses via musculoskeletal biomechanical modelling.

Do forelimb shape and peak forces co-vary in strepsirrhines?

OBJECTIVES

In this study, we explore whether ground reaction forces recorded during horizontal walking co-vary with the shape of the long bones of the forelimb in strepsirrhines. To do so, we quantify (1) the shape of the shaft and articular surfaces of each long bone of the forelimb, (2) the peak vertical, mediolateral, and horizontal ground reaction forces applied by the forelimb during arboreal locomotion, and (3) the relationship between the shape of the forelimb and peak forces.

MATERIALS AND METHODS

Geometric morphometric approaches were used to quantify the shape of the bones. Kinetic data were collected during horizontal arboreal walking in eight species of strepsirrhines that show variation in habitual substrate use and morphology of the forelimb. These data were then used to explore the links between locomotor behavior, morphology, and mechanics using co-variation analyses in a phylogenetic framework.

RESULTS

Our results show significant differences between slow quadrupedal climbers (lorises), vertical clinger and leapers (sifaka), and active arboreal quadrupeds (ring-tailed lemur, ruffed lemur) in both ground reaction forces and the shape of the long bones of the forelimb, with the propulsive and medially directed peak forces having the highest impact on the shape of the humerus. Co-variation between long bone shape and ground reaction forces was detected in both the humerus and ulna even when accounting for differences in body mass.

DISCUSSION

These results demonstrate the importance of considering limb-loading beyond just peak vertical force, or substrate reaction force. A re-evaluation of osseous morphology and functional interpretations is necessary in light of these findings.

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.

Collagen Fiber Orientation in Primate Long Bones

Studies of variation in orientation of collagen fibers within bone have lead to the proposition that these are preferentially aligned to accommodate different kinds of load, with tension best resisted by fibers aligned longitudinally relative to the load, and compression best resisted by transversely aligned fibers. However, previous studies have often neglected to consider the effect of developmental processes, including constraints on collagen fiber orientation (CFO), particularly in primary bone. Here we use circularly polarized light microscopy to examine patterns of CFO in cross-sections from the midshaft femur, humerus, tibia, radius, and ulna in a range of living primate taxa with varied body sizes, phylogenetic relationships and positional behaviors. We find that a preponderance of longitudinally oriented collagen is characteristic of both periosteal primary and intracortically remodeled bone. Where variation does occur among groups, it is not simply understood via interpretations of mechanical loads, although prioritized adaptations to tension and/or shear are considered. While there is some suggestion that CFO may correlate with body size, this relationship is neither consistent nor easily explicable through consideration of size-related changes in mechanical adaptation. The results of our study indicate that there is no clear relationship between CFO and phylogenetic status. One of the principle factors accounting for the range of variation that does exist is primary tissue type, where slower depositing bone is more likely to comprise a larger proportion of oblique to transverse collagen fibers. Anat Rec, 300:1189-1207, 2017. © 2017 Wiley Periodicals, Inc.

Horticultural activity predicts later localized limb status in a contemporary pre-industrial population

OBJECTIVES

Modern humans may have gracile skeletons due to low physical activity levels and mechanical loading. Tests using pre-historic skeletons are limited by the inability to assess behavior directly, while modern industrialized societies possess few socio-ecological features typical of human evolutionary history. Among Tsimane forager-horticulturalists, we test whether greater activity levels and, thus, increased loading earlier in life are associated with greater later-life bone status and diminished age-related bone loss.

MATERIALS AND METHODS

We used quantitative ultrasonography to assess radial and tibial status among adults aged 20+ years (mean ± SD age = 49 ± 15; 52% female). We conducted systematic behavioral observations to assess earlier-life activity patterns (mean time lag between behavioural observation and ultrasound = 12 years). For a subset of participants, physical activity was again measured later in life, via accelerometry, to determine whether earlier-life time use is associated with later-life activity levels. Anthropometric and demographic data were collected during medical exams.

RESULTS

Structural decline with age is reduced for the tibia (female: -0.25 SDs/decade; male: 0.05 SDs/decade) versus radius (female: -0.56 SDs/decade; male: -0.20 SDs/decade), which is expected if greater loading mitigates bone loss. Time allocation to horticulture, but not hunting, positively predicts later-life radial status (β_Horticulture  = 0.48, p = 0.01), whereas tibial status is not significantly predicted by subsistence or sedentary leisure participation.

DISCUSSION

Patterns of activity- and age-related change in bone status indicate localized osteogenic responses to loading, and are generally consistent with the logic of bone functional adaptation. Nonmechanical factors related to subsistence lifestyle moderate the association between activity patterns and bone structure.

Utility of osteon circularity for determining species and interpreting load history in primates and nonprimates

OBJECTIVES

Histomorphological analyses of bones are used to estimate an individual's chronological age, interpret a bone's load history, and differentiate species. Among various histomorphological characteristics that can influence mechanical properties of cortical bone, secondary osteon (Haversian system) population density and predominant collagen fiber orientation are particularly important. Cross-sectional shape characteristics of secondary osteons (On.Cr = osteon circularity, On.El = osteon ellipticality) are considered helpful in these contexts, but more robust proof is needed. We sought to determine if variations in osteon shape characteristics are sufficient for accurately differentiating species, load-complexity categories, and regional habitual strain-mode distributions (e.g., tension vs. compression regions).

MATERIALS AND METHODS

Circularly polarized light images were obtained from 100-micron transverse sections from diaphyses of adult deer calcanei; sheep calcanei, radii, and tibiae; equine calcanei, radii, and third metacarpals (MC3s); chimpanzee femora; and human femora and fibulae. Osteon cross-sectional area (On.Ar), On.Cr, and On.El were quantified indiscriminately and in the contexts of load-complexity and regional strain-mode distributions.

RESULTS

On.Cr and On.El, when examined independently in terms of all data, or mean (nested) data, for each bone, exceeded 80% accuracy in the inter-species comparisons only with respect to distinguishing humans from nonhumans. Correct classification among the nonhuman species was <70%. When On.Cr and On.El were coupled together and with On.Ar in discriminant function analyses (nested and unnested data) there were high misclassifications in all but human vs. nonhuman comparisons.

DISCUSSION

Frequent misclassifications in nonhuman comparisons might reflect influences of habitual load complexity and/or strain-mode distributions, or other factors not accounted for by these two considerations.

In vitro bone strain distributions in a sample of primate pelves

The pelvis is a critical link in the hindlimb locomotor system and has a central role in resisting loads associated with locomotion, but our understanding of its structural biomechanics is quite limited. Empirical data on how the pelvis responds to the loads it encounters are important for understanding pelvic adaptation to locomotion, and for testing hypotheses regarding how the pelvis is adapted to its mechanical demands. This paper presents in vitro strain gauge data on a sample of monkey and ape cadaveric specimens (Macaca, Papio, Ateles, Hylobates), and assesses strain magnitudes and distributions through the bones of the pelvis: the ilium, ischium and pubis. Pelves were individually mounted in a materials testing system, loads were applied across three hindlimb angular positions, and strains were recorded from 18 locations on the pelvic girdle. Peak principal strains range from 2000 to 3000 με, similar to peak strains recorded from other mammals in vivo. Although previous work has suggested that the bones of the pelvis may act as bent beams, this study suggests that there are likely additional loading regimes superimposed on bending. Specifically, these data suggest that the ilium is loaded in axial compression and torsion, the ischium in torsion, the pubic rami in mediolateral bending, and the pubic symphysis is loaded in a combination of compression and torsion. Compressive strains dominate the pelves of all species representatives. Shear strains change with limb position; hip flexion at 45° induces smaller shear strains than mid-stance (90°) or hip extension (105°). The pelvic girdle is a complex structure that does not lend itself easily to modeling, but finite element analyses may prove useful to generate and refine hypotheses of pelvic biomechanics.

Exercise-induced bone formation is poorly linked to local strain magnitude in the sheep tibia

Functional interpretations of limb bone structure frequently assume that diaphyses adjust their shape by adding bone primarily across the plane in which they are habitually loaded in order to minimize loading-induced strains. Here, to test this hypothesis, we characterize the in vivo strain environment of the sheep tibial midshaft during treadmill exercise and examine whether this activity promotes bone formation disproportionately in the direction of loading in diaphyseal regions that experience the highest strains. It is shown that during treadmill exercise, sheep tibiae were bent in an anteroposterior direction, generating maximal tensile and compressive strains on the anterior and posterior shaft surfaces, respectively. Exercise led to significantly increased periosteal bone formation; however, rather than being biased toward areas of maximal strains across the anteroposterior axis, exercise-related osteogenesis occurred primarily around the medial half of the shaft circumference, in both high and low strain regions. Overall, the results of this study demonstrate that loading-induced bone growth is not closely linked to local strain magnitude in every instance. Therefore, caution is necessary when bone shaft shape is used to infer functional loading history in the absence of in vivo data on how bones are loaded and how they actually respond to loading.

Torsion and antero-posterior bending in the in vivo human tibia loading regimes during walking and running

Bending, in addition to compression, is recognized to be a common loading pattern in long bones in animals. However, due to the technical difficulty of measuring bone deformation in humans, our current understanding of bone loading patterns in humans is very limited. In the present study, we hypothesized that bending and torsion are important loading regimes in the human tibia. In vivo tibia segment deformation in humans was assessed during walking and running utilizing a novel optical approach. Results suggest that the proximal tibia primarily bends to the posterior (bending angle: 0.15°-1.30°) and medial aspect (bending angle: 0.38°-0.90°) and that it twists externally (torsion angle: 0.67°-1.66°) in relation to the distal tibia during the stance phase of overground walking at a speed between 2.5 and 6.1 km/h. Peak posterior bending and peak torsion occurred during the first and second half of stance phase, respectively. The peak-to-peak antero-posterior (AP) bending angles increased linearly with vertical ground reaction force and speed. Similarly, peak-to-peak torsion angles increased with the vertical free moment in four of the five test subjects and with the speed in three of the test subjects. There was no correlation between peak-to-peak medio-lateral (ML) bending angles and ground reaction force or speed. On the treadmill, peak-to-peak AP bending angles increased with walking and running speed, but peak-to-peak torsion angles and peak-to-peak ML bending angles remained constant during walking. Peak-to-peak AP bending angle during treadmill running was speed-dependent and larger than that observed during walking. In contrast, peak-to-peak tibia torsion angle was smaller during treadmill running than during walking. To conclude, bending and torsion of substantial magnitude were observed in the human tibia during walking and running. A systematic distribution of peak amplitude was found during the first and second parts of the stance phase.

Geometric properties and comparative biomechanics of Homo floresiensis mandibles

The hypodigm of Homo floresiensis from the cave of Liang Bua on Flores Island in the archipelago of Indonesia includes two mandibles (LB1/2 and LB6/1). The morphology of their symphyses and corpora has been described as sharing similarities with both australopiths and early Homo despite their Late Pleistocene age. Although detailed morphological comparisons of these mandibles with those of modern and fossil hominin taxa have been made, a functional analysis in the context of masticatory biomechanics has yet to be performed. Utilizing data on cortical bone geometry from computed tomography scans, we compare the mechanical attributes of the LB1 and LB6 mandibles with samples of modern Homo, Pan, Pongo, and Gorilla, as well as fossil samples of Paranthropus robustus, Australopithecus africanus and South African early Homo. Structural stiffness measures were derived from the geometric data to provide relative measures of mandibular corpus strength under hypothesized masticatory loading regimes. These mechanical variables were evaluated relative to bone area, mandibular length and estimates of body size to assess their functional affinities and to test the hypothesis that the Liang Bua mandibles can be described as scaled-down variants of either early hominins or modern humans. Relative to modern hominoids, the H. floresiensis material appears to be relatively strong in terms of rigidity in torsion and transverse bending, but is relatively weak under parasagittal bending. Thus, they are 'robust' relative to modern humans (and comparable with australopiths) under some loads but not others. Neither LB1 nor LB6 can be described simply as 'miniaturized' versions of modern human jaws since mandible length is more or less equivalent in Homo sapiens and H. floresiensis. The mechanical attributes of the Liang Bua mandibles are consistent with previous inferences that masticatory loads were reduced relative to australopiths but remained elevated relative to modern Homo.

Cortical bone thickness can adapt locally to muscular loading while changing with age

Mechanical loading of muscle action is concentrated at muscle attachment sites; thus there may be a potential for site-specific variation in cortical bone thickness. Humeri from an early 20th-century Finnish (Helsinki) and two medieval English (Newcastle, Blackgate and York, Barbican) populations were subjected to pQCT scanning to calculate site-specific cross-sectional cortical bone area (CA) for four locations and to measure cortical thickness at muscle attachment sites and non-attachment sites. We found that CA at 80% of humerus length was significantly reduced compared to more distal cross-sections, which can be due to reduced stresses at the proximal shaft. The principal direction of loading at 80% humerus length was towards mediolateral plane, likely due to fixing the humerus close to the torso. At 35% the main direction of loading was towards anteroposterior plane, reflecting elbow flexing forces. The principal direction of loading varied between populations, sides and sexes at 50% humerus length due to preference between elbow and shoulder joint; thus this location might be useful when trying to infer differences in activity. These changes are likely due to overall shaft adaptation to forces acting at the humerus. In addition, we found a potential for site-specific variation in cortical thickness; cortical bone at muscle attachment sites was significantly thicker compared to non-attachment sites. Lastly, CA at 35% of humerus length and cortical thickness at non-attachment sites decreased with age. These results underline the importance of muscle loading for bone mass preservation as well as indicate that a site-specific variation of bone mass is possible.

In vivo bone strain and finite element modeling of the mandible of Alligator mississippiensis

Forces experienced during feeding are thought to strongly influence the morphology of the vertebrate mandible; in vivo strain data are the most direct evidence for deformation of the mandible induced by these loading regimes. Although many studies have documented bone strains in the mammalian mandible, no information is available on strain magnitudes, orientations or patterns in the sauropsid lower jaw during feeding. Furthermore, strain gage experiments record the mechanical response of bone at a few locations, not across the entire mandible. In this paper, we present bone strain data recorded at various sites on the lower jaw of Alligator mississippiensis during in vivo feeding experiments. These data are used to understand how changes in loading regime associated with changes in bite location are related to changes in strain regime on the working and balancing sides of the mandible. Our results suggest that the working side mandible is bent dorsoventrally and twisted about its long-axis during biting, and the balancing side experiences primarily dorsoventral bending. Strain orientations are more variable on the working side than on the balancing side with changes in bite point and between experiments; the balancing side exhibits higher strain magnitudes. In the second part of this paper, we use principal strain orientations and magnitudes recorded in vivo to evaluate a finite element model of the alligator mandible. Our comparison demonstrates that strain orientations and mandibular deformation predicted by the model closely match in vivo results; however, absolute strain magnitudes are lower in the finite element model.

Measuring strain using digital image correlation of second harmonic generation images

The micromechanical environment of bone is crucial to understanding both bone fracture and mechanobiological responses of osteocytes, yet few techniques exist that are capable of measuring strains on the micrometer scale. A method for measuring micrometer level strains has been developed based on digital image correlation (DIC) of second harmonic generation microscopy (SHGM) images. Bovine tibias milled into thin sections were imaged using SHGM under loads of 0 and 15 MPa. Strains were measured using DIC and compared to applied strain values. First and second principal strains decreased in magnitude as the analysis region area increased from 1750 µm(2) to 60,920 µm(2), converging to 1.23 ± 0.74 and -0.745 ± 0.9816 times the applied strain respectively. A representative sample histogram revealed regions of pure tensile and compressive strain, and that strains were highly heterogeneous ranging from 8410 to -8840 microstrain for an applied 2870 microstrain. Comparison with applied strain measures suggested that analysis sizes of 1750 µm(2) and greater were measuring strains on the tissue scale, and higher resolution is required for collagen fibrillar strains. Regions of low SHGM intensity ("dark" regions) were seen which are believed to be lacunar and perilacunar regions of low collagen density. However, no significant differences in strain magnitude were present in dark regions versus regions of high signal intensity. The proposed technique is effective for strains on the size order of bone microarchitecture, and would be useful for studies into the mechanical microenvironment during loading. The technique also has potential for in vivo studies in small animal models.

A comparison of proximal humeral cancellous bone of great apes and humans

The shoulder is the most mobile joint in the primate body, and is involved in both locomotor and manipulative activities. The presumed functional sensibility of trabecular bone can offer a way of decoding the activities to which the forelimbs of fossil primates were subjected. We examine the proximal humeral trabecular architecture in a relatively closely related group of similarly sized hominids (Pongo pygmaeus, Pan troglodytes, and Homo sapiens), in order to evaluate the effect of diverging habitual motion behaviors of the shoulder complex in a coherent phylogenetic group. In order to characterize and compare the humeral trabecular architectures of the three species, we imaged a large sample by high-resolution computed tomography (HrCT) and quantified their trabecular architectures by standard bone 3D morphometric parameters. Univariate statistical analysis was performed, showing significant differences among the species. However, univariate statistics could not highlight the structural particularity in the cancellous bone of each species. A principal component analysis also showed clear separation of the three taxa and enabled a structural characterization of the humeral trabecular bone of each species. We conclude that the differences in the architectural setup of the three hominids likely reflect multiple differences in their habitual activity patterns of their shoulder joint, although individual structural features are difficult to relate to specific loading conditions.

Locomotor loading mechanics in the hindlimbs of tegu lizards (Tupinambis merianae): comparative and evolutionary implications

Skeletal elements are usually able to withstand several times their usual load before they yield, and this ratio is known as the bone's safety factor. Limited studies on amphibians and non-avian reptiles have shown that they have much higher limb bone safety factors than birds and mammals. It has been hypothesized that this difference is related to the difference in posture between upright birds and mammals and sprawling ectotherms; however, limb bone loading data from a wider range of sprawling species are needed in order to determine whether the higher safety factors seen in amphibians and non-avian reptiles are ancestral or derived conditions. Tegus (family Teiidae) are an ideal lineage with which to expand sampling of limb bone loading mechanics for sprawling taxa, particularly for lizards, because they are from a different clade than previously sampled iguanas and exhibit different foraging and locomotor habits (actively foraging carnivore versus burst-activity herbivore). We evaluated the mechanics of locomotor loading for the femur of the Argentine black and white tegu (Tupinambus merianae) using three-dimensional measurements of the ground reaction force and hindlimb kinematics, in vivo bone strains and femoral mechanical properties. Peak bending stresses experienced by the femur were low (tensile: 10.4 ± 1.1 MPa; compressive: -17.4 ± 0.9 MPa) and comparable to those in other reptiles, with moderate shear stresses and strains also present. Analyses of peak femoral stresses and strains led to estimated safety factor ranges of 8.8-18.6 in bending and 7.8-17.5 in torsion, both substantially higher than typical for birds and mammals but similar to other sprawling tetrapods. These results broaden the range of reptilian and amphibian taxa in which high femoral safety factors have been evaluated and further indicate a trend for the independent evolution of lower limb bone safety factors in endothermic taxa.

Loading mechanics of the femur in tiger salamanders (Ambystoma tigrinum) during terrestrial locomotion

Salamanders are often used as representatives of the basal tetrapod body plan in functional studies, but little is known about the loads experienced by their limb bones during locomotion. Although salamanders' slow walking speeds might lead to low locomotor forces and limb bone stresses similar to those of non-avian reptiles, their highly sprawled posture combined with relatively small limb bones could produce elevated limb bone stresses closer to those of avian and mammalian species. This study evaluates the loads on the femur of the tiger salamander (Ambystoma tigrinum) during terrestrial locomotion using three-dimensional measurements of the ground reaction force (GRF) and hindlimb kinematics, as well as anatomical measurements of the femur and hindlimb muscles. At peak stress (29.8 ± 2.0% stance), the net GRF magnitude averaged 0.42 body weights and was directed nearly vertically for the middle 20-40% of the contact interval, essentially perpendicular to the femur. Although torsional shear stresses were significant (4.1 ± 0.3 MPa), bending stresses experienced by the femur were low compared with other vertebrate lineages (tensile: 14.9 ± 0.8 MPa; compressive: -18.9 ± 1.0 MPa), and mechanical property tests indicated yield strengths that were fairly standard for tetrapods (157.1 ± 3.7 MPa). Femoral bending safety factors (10.5) were considerably higher than values typical for birds and mammals, and closer to the elevated values calculated for reptilian species. These results suggest that high limb bone safety factors may have an ancient evolutionary history, though the underlying cause of high safety factors (e.g. low limb bone loads, high bone strength or a combination of the two) may vary among lineages.

In vivo strains in the femur of the Virginia opossum (Didelphis virginiana) during terrestrial locomotion: testing hypotheses of evolutionary shifts in mammalian bone loading and design

Terrestrial locomotion can impose substantial loads on vertebrate limbs. Previous studies have shown that limb bones from cursorial species of eutherian mammals experience high bending loads with minimal torsion, whereas the limb bones of non-avian reptiles (and amphibians) exhibit considerable torsion in addition to bending. It has been hypothesized that these differences in loading regime are related to the difference in limb posture between upright mammals and sprawling reptiles, and that the loading patterns observed in non-avian reptiles may be ancestral for tetrapod vertebrates. To evaluate whether non-cursorial mammals show loading patterns more similar to those of sprawling lineages, we measured in vivo strains in the femur during terrestrial locomotion of the Virginia opossum (Didelphis virginiana), a marsupial that uses more crouched limb posture than most mammals from which bone strains have been recorded, and which belongs to a clade phylogenetically between reptiles and the eutherian mammals studied previously. The presence of substantial torsion in the femur of opossums, similar to non-avian reptiles, would suggest that this loading regime likely reflects an ancestral condition for tetrapod limb bone design. Strain recordings indicate the presence of both bending and appreciable torsion (shear strain: 419.1 ± 212.8 με) in the opossum femur, with planar strain analyses showing neutral axis orientations that placed the lateral aspect of the femur in tension at the time of peak strains. Such mediolateral bending was unexpected for a mammal running with near-parasagittal limb kinematics. Shear strains were similar in magnitude to peak compressive axial strains, with opossum femora experiencing similar bending loads but higher levels of torsion compared with most previously studied mammals. Analyses of peak femoral strains led to estimated safety factor ranges of 5.1-7.2 in bending and 5.5-7.3 in torsion, somewhat higher than typical mammalian values for bending, but approaching typical reptilian values for shear. Loading patterns of opossum limb bones therefore appear intermediate in some respects between those of eutherian mammals and non-avian reptiles, providing further support for hypotheses that high torsion and elevated limb bone safety factors may represent persistent ancestral conditions in the evolution of tetrapod limb bone loading and design.

Metacarpal trabecular architecture variation in the chimpanzee (Pan troglodytes): Evidence for locomotion and tool-use?

Trabecular architecture was assessed by 3D micro-computed tomography from spherical volumes of interest located within the head and base of metacarpals (MC) 1, 2, and 5 from n = 19 adult common chimpanzees. Two subspecies, West African Pan troglodytes verus from the Taï Forest, Côte d'Ivoire (n = 12) and Central African P. t. troglodytes from Cameroon (n = 7), were studied. For the combined sample, the metacarpal head is distinguished by greater bone volume fraction across all metacarpals, though the MC 1 is distinctive in having thicker, more plate-like trabeculae. The architecture in the MC 2 and MC 5 can be related to strains associated with terrestrial knuckle-walking. In particular, the relatively robust MC 5 head architecture may result from functional loading incurred during braking and use of a palm-in hand posture. Examining differences between samples, we found that the Cameroon chimpanzees possess a more robust architecture across all metacarpals in the form of greater bone volume fraction, higher connectivity, and somewhat more plate-like structure. These differences are not explicable in terms of population distinctions in body size or daily travel distance, but possibly reflect a combination of more terrestrial knuckle-walking in the Cameroon sample and more diverse hand postures and precision handling required of nut-cracking in West African chimpanzees.

Characterizing gait induced normal strains in a murine tibia cortical bone defect model

The critical role that mechanical stimuli serve in mediating bone repair is recognized but incompletely understood. Further, previous attempts to understand this role have utilized application of externally applied mechanical loads to study the tissue's response. In this project, we have therefore endeavored to capitalize on bone's own consistently diverse loading environment to develop a novel model that would enable assessment of the influence of physiologically engendered mechanical stimuli on cortical defect repair. We used an inverse dynamics approach with finite element analysis (FEA) to first quantify normal strain distributions generated in mouse tibia during locomotion. The strain environment of the tibia, as previously reported for other long bones, was found to arise primarily due to bending and was consistent in orientation through the stance phase of gait. Based on these data, we identified three regions within a transverse cross-section of the mid-diaphysis as uniform locations of either peak tension, peak compression, or the neutral axis of bending (i.e. minimal strain magnitude). We then used FEA to quantify the altered strain environment that would be produced by a 0.6mm diameter cylindrical cortical bone defect at each diaphyseal site and, in an in situ study confirmed our ability to accurately place defects at the desired diaphyseal locations. The resulting model will enable the exploration of cortical bone healing within the context of physiologically engendered mechanical strain.

Dynamic pressure patterns in the hands of olive baboons (Papio anubis) during terrestrial locomotion: implications for cercopithecoid primate hand morphology

Habitually terrestrial monkeys adopt digitigrade hand postures at slow speeds to increase effective forelimb length and reduce distal limb joint moments. As these primates move faster, however, their hands transition to a more palmigrade posture, which is likely associated with the inability of wrist and hand joints to resist higher ground reaction forces (GRF) associated with faster speeds. Transitioning to a palmigrade posture may serve to distribute GRFs over a larger surface area (i.e., increased palmar contact), ultimately reducing stresses in fragile hand bones. To test this hypothesis, dynamic palmar pressure data were collected on two adult baboons (Papio anubis) walking, running, and galloping across a runway integrated with a dynamic pressure mat (20 steps each; speed range: 0.46-4.0 m/s). Peak GRF, contact area, peak pressure, and pressure-time integral were quantified in two regions of the hand: fingers and palms (including metacarpal heads). At slower speeds with lower GRFs, the baboons use digitigrade postures resulting in small palmar contact area (largely across the metacarpal heads). At faster speeds with higher GRFs, they used less digitigrade hand postures resulting in increased palmar contact area. Finger contact area did not change across speeds. Despite higher GRFs at faster speeds, metacarpal pressure was moderated across speeds due to increased palmar contact area as animals transitioned from digitigrady to palmigrady. In contrast, the pressure in the fingers increased with faster speeds. Results indicate that the transition from digitigrady to palmigrady distributes increased forces over a larger palmar surface area. Such dynamic changes in palmar pressure likely moderate strain in the gracile bones of the hand, a structure that is integral not only for locomotion, but also feeding and social behaviors in primates.