In vivo bone strain and finite-element modeling of the craniofacial haft in catarrhine primates

Biomechanics in anthropology

Biomechanics is the set of tools that explain organismal movement and mechanical behavior and links the organism to the physicality of the world. As such, biomechanics can relate behaviors and culture to the physicality of the organism. Scale is critical to biomechanical analyses, as the constitutive equations that matter differ depending on the scale of the question. Within anthropology, biomechanics has had a wide range of applications, from understanding how we and other primates evolved to understanding the effects of technologies, such as the atlatl, and the relationship between identity, society, culture, and medical interventions, such as prosthetics. Like any other model, there is great utility in biomechanical models, but models should be used primarily for hypothesis testing and not data generation except in the rare case where models can be robustly validated. The application of biomechanics within anthropology has been extensive, and holds great potential for the future.

The relative size of the calcaneal tuber reflects heel strike plantigrady in African apes and humans

OBJECTIVES

The positional repertoire of the human-chimpanzee last common ancestor is critical for reconstructing the evolution of bipedalism. African apes and humans share a heel strike plantigrade foot posture associated with terrestriality. Previous research has established that modern humans have a relatively large and intrinsically robust calcaneal tuber equipped to withstand heel strike forces associated with bipedal walking and running. However, it is unclear whether African apes have a relatively larger calcaneal tuber than non-heel-striking primates, and how this trait might have evolved among anthropoids. Here, I test the hypothesis that heel-striking primates have a relatively larger calcaneal tuber than non-heel-striking primates.

METHODS

The comparative sample includes 331 individuals and 53 taxa representing hominoids, cercopithecoids, and platyrrhines. Evolutionary modeling was used to test for the effect of foot posture on the relative size of the calcaneal tuber in a phylogenetic framework that accounts for adaptation and inertia. Bayesian evolutionary modeling was used to identify selective regime shifts in the relative size of the calcaneal tuber among anthropoids.

RESULTS

The best fitting evolutionary model was a Brownian motion model with regime-dependent trends characterized by relatively large calcaneal tubers among African apes and humans. Evolutionary modeling provided support for an evolutionary shift toward a larger calcaneal tuber at the base of the African ape and human clade.

CONCLUSIONS

The results of this study support the view that African apes and humans share derived traits related to heel strike plantigrady, which implies that humans evolved from a semi-terrestrial quadrupedal ancestor.

Dynamic finite element modelling of the macaque mandible during a complete mastication gape cycle

Three-dimensional finite element models (FEMs) are powerful tools for studying the mechanical behaviour of the feeding system. Using validated, static FEMs we have previously shown that in rhesus macaques the largest food-related differences in strain magnitudes during unilateral postcanine chewing extend from the lingual symphysis to the endocondylar ridge of the balancing-side ramus. However, static FEMs only model a single time point during the gape cycle and probably do not fully capture the mechanical behaviour of the jaw during mastication. Bone strain patterns and moments applied to the mandible are known to vary during the gape cycle owing to variation in the activation peaks of the jaw-elevator muscles, suggesting that dynamic models are superior to static ones in studying feeding biomechanics. To test this hypothesis, we built dynamic FEMs of a complete gape cycle using muscle force data from in vivo experiments to elucidate the impact of relative timing of muscle force on mandible biomechanics. Results show that loading and strain regimes vary across the chewing cycle in subtly different ways for different foods, something which was not apparent in static FEMs. These results indicate that dynamic three-dimensional FEMs are more informative than static three-dimensional FEMs in capturing the mechanical behaviour of the jaw during feeding by reflecting the asymmetry in jaw-adductor muscle activations during a gape cycle. This article is part of the theme issue 'Food processing and nutritional assimilation in animals'.

Sagittal suture strain in capuchin monkeys (Sapajus and Cebus) during feeding

OBJECTIVES

Morphological variation in cranial sutures is used to infer aspects of primate feeding behavior, including diet, but strain regimes across sutures are not well documented. Our aim is to test hypotheses about sagittal suture morphology, strain regime, feeding behavior, and muscle activity relationships in robust Sapajus and gracile Cebus capuchin primates.

MATERIALS AND METHODS

Morphometrics of sinuosity in three regions of the sagittal suture were compared among museum specimens of Sapajus and Cebus, as well as in robust and gracile lab specimens. In vivo strains and bilateral electromyographic (EMG) activity were recorded from these regions in the temporalis muscles of capuchin primates while they fed on mechanically-varying foods.

RESULTS

Sapajus and the anterior suture region exhibited greater sinuosity than Cebus and posterior regions. In vivo data reveal minor differences in strain regime between robust and gracile phenotypes but show higher strain magnitudes in the middle suture region and higher tensile strains anteriorly. After gage location, feeding behavior has the most consistent and strongest impact on strain regime in the sagittal suture. Strain in the anterior suture has a high tension to compression ratio compared to the posterior region, especially during forceful biting in the robust Sapajus-like individual.

DISCUSSION

Sagittal suture complexity in robust capuchins likely reflects feeding behaviors associated with mechanically challenging foods. Sutural strain regimes in other anthropoid primates may also be affected by activity in feeding muscles.

Assessment of the mechanical role of cranial sutures in the mammalian skull: Computational biomechanical modelling of the rat skull

Cranial sutures are fibrocellular joints between the skull bones that are progressively replaced with bone throughout ontogeny, facilitating growth and cranial shape change. This transition from soft tissue to bone is reflected in the biomechanical properties of the craniofacial complex. However, the mechanical significance of cranial sutures has only been explored at a few localised areas within the mammalian skull, and as such our understanding of suture function in overall skull biomechanics is still limited. Here, we sought to determine how the overall strain environment is affected by the complex network of cranial sutures in the mammal skull. We combined two computational biomechanical methods, multibody dynamics analysis and finite element analysis, to simulate biting in a rat skull and compared models with and without cranial sutures. Our results show that including complex sutures in the rat model does not substantially change overall strain gradients across the cranium, particularly strain magnitudes in the bones overlying the brain. However, local variations in strain magnitudes and patterns can be observed in areas close to the sutures. These results show that, during feeding, sutures may be more important in some regions than others. Sutures should therefore be included in models that require accurate local strain magnitudes and patterns of cranial strain, particularly if models are developed for analysis of specific regions, such as the temporomandibular joint or zygomatic arch. Our results suggest that, for mammalian skulls, cranial sutures might be more important for allowing brain expansion during growth than redistributing biting loads across the cranium in adults.

Does the model reflect the system? When two-dimensional biomechanics is not 'good enough'

Models are mathematical representations of systems, processes or phenomena. In biomechanics, finite-element modelling (FEM) can be a powerful tool, allowing biologists to test form-function relationships in silico, replacing or extending results of in vivo experimentation. Although modelling simplifications and assumptions are necessary, as a minimum modelling requirement the results of the simplified model must reflect the biomechanics of the modelled system. In cases where the three-dimensional mechanics of a structure are important determinants of its performance, simplified two-dimensional modelling approaches are likely to produce inaccurate results. The vertebrate mandible is one among many three-dimensional anatomical structures routinely modelled using two-dimensional FE analysis. We thus compare the stress regimes of our published three-dimensional model of the chimpanzee mandible with a published two-dimensional model of the chimpanzee mandible and identify several fundamental differences. We then present a series of two-dimensional and three-dimensional FE modelling experiments that demonstrate how three key modelling parameters, (i) dimensionality, (ii) symmetric geometry, and (iii) constraints, affect deformation and strain regimes of the models. Our results confirm that, in the case of the primate mandible (at least), two-dimensional FEM fails to meet this minimum modelling requirement and should not be used to draw functional, ecological or evolutionary conclusions.

Mechanical compensation in the evolution of the early hominin feeding apparatus

Australopiths, a group of hominins from the Plio-Pleistocene of Africa, are characterized by derived traits in their crania hypothesized to strengthen the facial skeleton against feeding loads and increase the efficiency of bite force production. The crania of robust australopiths are further thought to be stronger and more efficient than those of gracile australopiths. Results of prior mechanical analyses have been broadly consistent with this hypothesis, but here we show that the predictions of the hypothesis with respect to mechanical strength are not met: some gracile australopith crania are as strong as that of a robust australopith, and the strength of gracile australopith crania overlaps substantially with that of chimpanzee crania. We hypothesize that the evolution of cranial traits that increased the efficiency of bite force production in australopiths may have simultaneously weakened the face, leading to the compensatory evolution of additional traits that reinforced the facial skeleton. The evolution of facial form in early hominins can therefore be thought of as an interplay between the need to increase the efficiency of bite force production and the need to maintain the structural integrity of the face.

More Challenging Diets Sustain Feeding Performance: Applications Toward the Captive Rearing of Wildlife

The rescue and rehabilitation of young fauna is of substantial importance to conservation. However, it has been suggested that incongruous diets offered in captive environments may alter craniofacial morphology and hinder the success of reintroduced animals. Despite these claims, to what extent dietary variation throughout ontogeny impacts intrapopulation cranial biomechanics has not yet been tested. Here, finite element models were generated from the adult crania of 40 rats (n = 10 per group) that were reared on 4 different diet regimes and stress magnitudes compared during incisor bite simulations. The diets consisted of (1) exclusively hard pellets from weaning, (2) exclusively soft ground pellet meal from weaning, (3) a juvenile switch from pellets to meal, and (4) a juvenile switch from meal to pellets. We hypothesized that a diet of exclusively soft meal would result in the weakest adult skulls, represented by significantly greater stress magnitudes at the muzzle, palate, and zygomatic arch. Our hypothesis was supported at the muzzle and palate, indicating that a diet limited to soft food inhibits bone deposition throughout ontogeny. This finding presents a strong case for a more variable and challenging diet during development. However, rather than the "soft" diet group resulting in the weakest zygomatic arch as predicted, this region instead showed the highest stress among rats that switched as juveniles from hard pellets to soft meal. We attribute this to a potential reduction in number and activity of osteoblasts, as demonstrated in studies of sudden and prolonged disuse of bone. A shift to softer foods in captivity, during rehabilitation after injury in the wild for example, can therefore be detrimental to healthy development of the skull in some growing animals, potentially increasing the risk of injury and impacting the ability to access full ranges of wild foods upon release. We suggest captive diet plans consider not just nutritional requirements but also food mechanical properties when rearing wildlife to adulthood for reintroduction.

Fracture Fixation Technique and Chewing Side Impact Jaw Mechanics in Mandible Fracture Repair

Lower jaw (mandible) fractures significantly impact patient health and well‐being due to pain and difficulty eating, but the best technique for repairing the most common subtype—angle fractures—and rehabilitating mastication is unknown. Our study is the first to use realistic in silico simulation of chewing to quantify the effects of Champy and biplanar techniques of angle fracture fixation. We show that more rigid, biplanar fixation results in lower strain magnitudes in the miniplates, the bone around the screws, and in the fracture zone, and that the mandibular strain regime approximates the unfractured condition. Importantly, the strain regime in the fracture zone is affected by chewing laterality, suggesting that both fixation type and the patient's post‐fixation masticatory pattern—ipsi‐ or contralateral to the fracture— impact the bone healing environment. Our study calls for further investigation of the impact of fixation technique on chewing behavior. Research that combines in vivo and in silico approaches can link jaw mechanics to bone healing and yield more definitive recommendations for fixation, hardware, and postoperative rehabilitation to improve outcomes. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

The influence of jaw-muscle fibre-type phenotypes on estimating maximum muscle and bite forces in primates

Numerous anthropological studies have been aimed at estimating jaw-adductor muscle forces, which, in turn, are used to estimate bite force. While primate jaw adductors show considerable intra- and intermuscular heterogeneity in fibre types, studies generally model jaw-muscle forces by treating the jaw adductors as either homogeneously slow or homogeneously fast muscles. Here, we provide a novel extension of such studies by integrating fibre architecture, fibre types and fibre-specific tensions to estimate maximum muscle forces in the masseter and temporalis of five anthropoid primates: Sapajus apella (N = 3), Cercocebus atys (N = 4), Macaca fascicularis (N = 3), Gorilla gorilla (N = 1) and Pan troglodytes (N = 2). We calculated maximum muscle forces by proportionally adjusting muscle physiological cross-sectional areas by their fibre types and associated specific tensions. Our results show that the jaw adductors of our sample ubiquitously express MHC α-cardiac, which has low specific tension, and hybrid fibres. We find that treating the jaw adductors as either homogeneously slow or fast muscles potentially overestimates average maximum muscle forces by as much as approximately 44%. Including fibre types and their specific tensions is thus likely to improve jaw-muscle and bite force estimates in primates.

The cranial biomechanics and feeding performance of Homo floresiensis

Homo floresiensis is a small-bodied hominin from Flores, Indonesia, that exhibits plesiomorphic dentognathic features, including large premolars and a robust mandible, aspects of which have been considered australopith-like. However, relative to australopith species, H. floresiensis exhibits reduced molar size and a cranium with diminutive midfacial dimensions similar to those of later Homo, suggesting a reduction in the frequency of forceful biting behaviours. Our study uses finite-element analysis to examine the feeding biomechanics of the H. floresiensis cranium. We simulate premolar (P³) and molar (M²) biting in a finite-element model (FEM) of the H. floresiensis holotype cranium (LB1) and compare the mechanical results with FEMs of chimpanzees, modern humans and a sample of australopiths (MH1, Sts 5, OH5). With few exceptions, strain magnitudes in LB1 resemble elevated levels observed in modern Homo. Our analysis of LB1 suggests that H. floresiensis could produce bite forces with high mechanical efficiency, but was subject to tensile jaw joint reaction forces during molar biting, which perhaps constrained maximum postcanine bite force production. The inferred feeding biomechanics of H. floresiensis closely resemble modern humans, suggesting that this pattern may have been present in the last common ancestor of Homo sapiens and H. floresiensis.

Anthroengineering: an independent interdisciplinary field

In recent decades, funding agencies, institutes and professional bodies have recognized the profound benefits of transdisciplinarity in tackling targeted research questions. However, once questions are answered, the previously abundant support often dissolves. As such, the long-term benefits of these transdisciplinary approaches are never fully achieved. Over the last several decades, the integration of anthropology and engineering through inter- and multidisciplinary work has led to advances in fields such as design, human evolution and medical technologies. The lack of formal recognition, however, of this transdisciplinary approach as a unique entity rather than a useful tool or a subfield makes it difficult for researchers to establish laboratories, secure permanent jobs, fund long-term research programmes and train students in this approach. To facilitate the growth and development and witness the long-term benefits of this approach, we propose the integration of anthropology and engineering be recognized as a new, independent field known as anthroengineering. We present a working definition for anthroengineering and examples of how anthroengineering has been used. We discuss the necessity of recognizing anthroengineering as a unique field and explore potential novel applications. Finally, we discuss the future of anthroengineering, highlighting avenues for moving the field forward.

Comparative biomechanics of the Pan and Macaca mandibles during mastication: finite element modelling of loading, deformation and strain regimes

The mechanical behaviour of the mandibles of Pan and Macaca during mastication was compared using finite element modelling. Muscle forces were calculated using species-specific measures of physiological cross-sectional area and scaled using electromyographic estimates of muscle recruitment in Macaca. Loading regimes were compared using moments acting on the mandible and strain regimes were qualitatively compared using maps of principal, shear and axial strains. The enlarged and more vertically oriented temporalis and superficial masseter muscles of Pan result in larger sagittal and transverse bending moments on both working and balancing sides, and larger anteroposterior twisting moments on the working side. The mandible of Pan experiences higher principal strain magnitudes in the ramus and mandibular prominence, higher transverse shear strains in the top of the symphyseal region and working-side corpus, and a predominance of sagittal bending-related strains in the balancing-side mandible. This study lays the foundation for a broader comparative study of Hominidae mandibular mechanics in extant and fossil hominids using finite element modelling. Pan's larger and more vertical masseter and temporalis may make it a more suitable model for hominid mandibular biomechanics than Macaca.

Comparative cranial biomechanics in two lizard species: impact of variation in cranial design

Cranial morphology in lepidosaurs is highly disparate and characterised by the frequent loss or reduction of bony elements. In varanids and geckos, the loss of the postorbital bar is associated with changes in skull shape, but the mechanical principles underlying this variation remain poorly understood. Here, we sought to determine how the overall cranial architecture and the presence of the postorbital bar relate to the loading and deformation of the cranial bones during biting in lepidosaurs. Using computer-based simulation techniques, we compared cranial biomechanics in the varanid Varanus niloticus and the teiid Salvator merianae, two large, active foragers. The overall strain magnitude and distribution across the cranium were similar in the two species, despite lower strain gradients in V. niloticus. In S. merianae, the postorbital bar is important for resistance of the cranium to feeding loads. The postorbital ligament, which in varanids partially replaces the postorbital bar, does not affect bone strain. Our results suggest that the reduction of the postorbital bar impaired neither biting performance nor the structural resistance of the cranium to feeding loads in V. niloticus. Differences in bone strain between the two species might reflect demands imposed by feeding and non-feeding functions on cranial shape. Beyond variation in cranial bone strain related to species-specific morphological differences, our results reveal that similar mechanical behaviour is shared by lizards with distinct cranial shapes. Contrary to the situation in mammals, the morphology of the circumorbital region, calvaria and palate appears to be important for withstanding high feeding loads in these lizards.

Summary:In vivo measurements and computer-based simulations of the cranial mechanics of two large lizards indicate that similar mechanical behaviour is shared by lizards with distinct cranial architecture, and show the importance of the postorbital bar in resisting the feeding loads.

Ingestive behaviors in bearded capuchins (Sapajus libidinosus)

The biomechanical and adaptive significance of variation in craniodental and mandibular morphology in fossil hominins is not always clear, at least in part because of a poor understanding of how different feeding behaviors impact feeding system design (form–function relationships). While laboratory studies suggest that ingestive behaviors produce variable loading, stress, and strain regimes in the cranium and mandible, understanding the relative importance of these behaviors for feeding system design requires data on their use in wild populations. Here we assess the frequencies and durations of manual, ingestive, and masticatory behaviors from more than 1400 observations of feeding behaviors video-recorded in a wild population of bearded capuchins (Sapajus libidinosus) at Fazenda Boa Vista in Piauí, Brazil. Our results suggest that ingestive behaviors in wild Sapajus libidinosus were used for a range of food material properties and typically performed using the anterior dentition. Coupled with previous laboratory work indicating that ingestive behaviors are associated with higher mandibular strain magnitudes than mastication, these results suggest that ingestive behaviors may play an important role in craniodental and mandibular design in capuchins and may be reflected in robust adaptations in fossil hominins.

Nonlinear dynamic analysis of the pure "buckling" mechanism during blow-out trauma of the human orbit

Considering the interplay between orbital bones and intraorbital soft tissues, commonly accepted patterns of the blow-out type of trauma within the human orbit require more thorough investigation to assess the minimal health-threatening impact value. Two different three-dimensional finite element method (FEM) models of the human orbital region were developed to simulate the pure “buckling” mechanism of orbital wall fracture in two variants: the model of orbital bone elements and the model of orbital bone, orbit and intraorbital tissue elements. The mechanical properties of the so-defined numerical skull fragment were applied to the model according to the unique laboratory tensile stress tests performed on small and fragile specimens of orbital bones as well as using the data available in the literature. The nonlinear transient analysis of the contact problem between bodies that differ substantially in terms of the Young’s modulus was carried out to investigate the interaction of different bodies within an instant injury. Potential damage areas were found within the lower orbital wall as well as the destructive load values for both FEM skull models (7,660 N and 8,520 N). Moreover, numerical simulations were validated by comparing them with computed tomography scans of real injuries.

Biomechanics of the mandible of Macaca mulatta during the power stroke of mastication: Loading, deformation, and strain regimes and the impact of food type

Mandible morphology has yet to yield definitive information on primate diet, probably because of poor understanding of mandibular loading and strain regimes, and overreliance on simple beam models of mandibular mechanics. We used a finite element model of a macaque mandible to test hypotheses about mandibular loading and strain regimes and relate variation in muscle activity during chewing on different foods to variation in strain regimes. The balancing-side corpus is loaded primarily by sagittal shear forces and sagittal bending moments. On the working side, sagittal bending moments, anteroposterior twisting moments, and lateral transverse bending moments all reach similar maxima below the bite point; sagittal shear is the dominant loading regime behind the bite point; and the corpus is twisted such that the mandibular base is inverted. In the symphyseal region, the predominant loading regimes are lateral transverse bending and negative twisting about a mediolateral axis. Compared with grape and dried fruit chewing, nut chewing is associated with larger sagittal and transverse bending moments acting on balancing- and working-side mandibles, larger sagittal shear on the working side, and larger twisting moments about vertical and transverse axes in the symphyseal region. Nut chewing is also associated with higher minimum principal strain magnitudes in the balancing-side posterior ramus; higher sagittal shear strain magnitudes in the working-side buccal alveolar process and the balancing-side oblique line, recessus mandibulae, and endocondylar ridge; and higher transverse shear strains in the symphyseal region, the balancing-side medial prominence, and the balancing-side endocondylar ridge. The largest food-related differences in maximum principal and transverse shear strain magnitudes are in the transverse tori and in the balancing-side medial prominence, extramolar sulcus, oblique line, and endocondylar ridge. Food effects on the strain regime are most salient in areas not traditionally investigated, suggesting that studies seeking dietary effects on mandible morphology might be looking in the wrong places.

The use of extruded finite-element models as a novel alternative to tomography-based models: a case study using early mammal jaws

Finite-element (FE) analysis has been used in palaeobiology to assess the mechanical performance of the jaw. It uses two types of models: tomography-based three-dimensional (3D) models (very accurate, not always accessible) and two-dimensional (2D) models (quick and easy to build, good for broad-scale studies, cannot obtain absolute stress and strain values). Here, we introduce extruded FE models, which provide fairly accurate mechanical performance results, while remaining low-cost, quick and easy to build. These are simplified 3D models built from lateral outlines of a relatively flat jaw and extruded to its average width. There are two types: extruded (flat mediolaterally) and enhanced extruded (accounts for width differences in the ascending ramus). Here, we compare mechanical performance values resulting from four types of FE models (i.e. tomography-based 3D, extruded, enhanced extruded and 2D) in Morganucodon and Kuehneotherium. In terms of absolute values, both types of extruded model perform well in comparison to the tomography-based 3D models, but enhanced extruded models perform better. In terms of overall patterns, all models produce similar results. Extruded FE models constitute a viable alternative to the use of tomography-based 3D models, particularly in relatively flat bones.

Warped finite element models predict whole shell failure in turtle shells

Finite element (FE) models have become increasingly popular in comparative biomechanical studies, with researchers continually developing methods such as 'warping' preexisting models to facilitate analyses. However, few studies have investigated how well FE models can predict biologically crucial whole-structure performance or whether 'warped' models can provide useful information about the mechanical behavior of actual specimens. This study addresses both of these issues through a validation of warped FE models of turtle shells. FE models for 40 turtle specimens were built using 3D landmark coordinates and thin-plate spline interpolations to warp preexisting turtle shell models. Each actual turtle specimen was loaded to failure, and the load at failure and mode of fracture were then compared with the behavior predicted by the models. Overall, the models performed very well, despite the fact that many simplifying assumptions were made for analysis. Regressions of observed on predicted loads were significant for the dataset as a whole, as well as in separate analyses within two turtle species, and the direction of fracture was generally consistent with the patterns of stresses observed in the models. This was true even when size (an important factor in determining strength) was removed from analyses - the models were also able to predict which shells would be relatively stronger or weaker. Although some residual variation remains unexplained, this study supports the idea that warped FE models run with simplifying assumptions at least can provide useful information for comparative biomechanical studies.

Analyzing the sclerocarpy adaptations of the Pitheciidae mandible

Primates are interpreted to be ancestrally adapted to frugivory, although some modern groups show clear adaptations to other diets. Among them, pitheciids stand out for specifically predating seeds. This dietary specialization is known as sclerocarpy and refers to the extraction of seeds from surrounding hard tissues using the anterior dentition followed by the mastication of seeds by the molars. It has been proposed that Callicebus-Pithecia-Chiropotes-Cacajao represent a morphocline of increasingly specialized anatomical traits for sclerocarpic foraging. This study addresses whether there is a sclerocarpic specialization gradient in the mandibular morphology of pitheciids. Finite element analysis (FEA) was used to simulate two biting scenarios and the obtained stress values were compared between different pitheciids. Geometric morphometrics (GM) were used to display the morphological variation of this group. No support was found for the morphocline hypothesis from a biomechanical viewpoint since all pitheciins showed similar stress values and on average Chiropotes rather than Cacajao exhibited the strongest mandible. From a morphological perspective, it was found that there is indeed relative "robusticity" continuum in the pitheciid mandible for some aspects of shape as expected for the morphocline hypothesis, but this gradient could be related to other factors rather than sclerocarpic specialization. The present results are expected to contribute to a better insight regarding the ecomorphological relationship between mandibular morphology and mechanical performance among pitheciids.

Structure-function covariation with nonfeeding ecological variables influences evolution of feeding specialization in Carnivora

Skull shape convergence is pervasive among vertebrates. Although this is frequently inferred to indicate similar functional underpinnings, neither the specific structure-function linkages nor the selective environments in which the supposed functional adaptations arose are commonly identified and tested. We demonstrate that nonfeeding factors relating to sexual maturity and precipitation-related arboreality also can generate structure-function relationships in the skulls of carnivorans (dogs, cats, seals, and relatives) through covariation with masticatory performance. We estimated measures of masticatory performance related to ecological variables that covary with cranial shape in the mammalian order Carnivora, integrating geometric morphometrics and finite element analyses. Even after accounting for phylogenetic autocorrelation, cranial shapes are significantly correlated to both feeding and nonfeeding ecological variables, and covariation with both variable types generated significant masticatory performance gradients. This suggests that mechanisms of obligate shape covariation with nonfeeding variables can produce performance changes resembling those arising from feeding adaptations in Carnivora.

Biting mechanics and niche separation in a specialized clade of primate seed predators

We analyzed feeding biomechanics in pitheciine monkeys (Pithecia, Chiropotes, Cacajao), a clade that specializes on hard-husked unripe fruit (sclerocarpy) and resistant seeds (seed predation). We tested the hypothesis that pitheciine crania are well-suited to generate and withstand forceful canine and molar biting, with the prediction that they generate bite forces more efficiently and better resist masticatory strains than the closely-related Callicebus, which does not specialize on unripe fruits and/or seeds. We also tested the hypothesis that Callicebus-Pithecia-Chiropotes-Cacajao represent a morphocline of increasing sclerocarpic specialization with respect to biting leverage and craniofacial strength, consistent with anterior dental morphology. We found that pitheciines have higher biting leverage than Callicebus and are generally more resistant to masticatory strain. However, Cacajao was found to experience high strain magnitudes in some facial regions. We therefore found limited support for the morphocline hypothesis, at least with respect to the mechanical performance metrics examined here. Biting leverage in Cacajao was nearly identical (or slightly less than) in Chiropotes and strain magnitudes during canine biting were more likely to follow a Cacajao-Chiropotes-Pithecia trend of increasing strength, in contrast to the proposed morphocline. These results could indicate that bite force efficiency and derived anterior teeth were selected for in pitheciines at the expense of increased strain magnitudes. However, our results for Cacajao potentially reflect reduced feeding competition offered by allopatry with other pitheciines, which allows Cacajao species to choose from a wider variety of fruits at various stages of ripeness, leading to reduction in the selection for robust facial features. We also found that feeding biomechanics in sympatric Pithecia and Chiropotes are consistent with data on food structural properties and observations of dietary niche separation, with the former being well-suited for the regular molar crushing of hard seeds and the latter better adapted for breaching hard fruits.

Bite force and cranial bone strain in four species of lizards

In vivo bone strain data provide direct evidence of strain patterns in the cranium during biting. Compared to mammals, in vivo bone strains in lizard skulls are poorly documented. This paper presents strain data from the skulls of Anolis equestris, Gekko gecko, Iguana iguana and Salvator merianae during transducer biting. Analysis of variance was used to investigate effects of bite force, bite point, diet, cranial morphology and cranial kinesis on strain magnitudes. Within individuals the most consistent determinants of variance in bone strain magnitudes are gage location and bite point, with the importance of bite force varying between individuals. Inter-site variance in strain magnitudes—strain gradient—is present in all individuals, and varies with bite point. Between individuals within species, variance in strain magnitude is driven primarily by variation in bite force, not gage location or bite point, suggesting that inter-individual variation in patterns of strain magnitude is minimal. Between species, variation in strain magnitudes is significantly impacted by bite force and species membership, as well as by interactions between gage location, species, and bite point. Independent of bite force, species differences in cranial strain magnitudes may reflect selection for different cranial morphology in relation to feeding function, but what these performance criteria are is not clear. The relatively low strain magnitudes in Iguana and Uromastyx compared to other lizards may be related to their herbivorous diet. Cranial kinesis and the presence or absence of postorbital and supratemporal bars are not important determinants of inter-specific variation in strain magnitudes.

The Biomechanics of Bony Facial "Buttresses" in South African Australopiths: An Experimental Study Using Finite Element Analysis

Australopiths exhibit a number of derived facial features that are thought to strengthen the face against high and/or repetitive loads associated with a diet that included mechanically challenging foods. Here, we use finite element analysis (FEA) to test hypotheses related to the purported strengthening role of the zygomatic root and "anterior pillar" in australopiths. We modified our previously constructed models of Sts 5 (Australopithecus africanus) and MH1 (A. sediba) to differ in the morphology of the zygomatic root, including changes to both the shape and positioning of the zygomatic root complex, in addition to creating variants of Sts 5 lacking anterior pillars. We found that both an expanded zygomatic root and the presence of "anterior pillars" reinforce the face against feeding loads. We also found that strain orientations are most compatible with the hypothesis that the pillar evolved to resist loads associated with premolar loading, and that this morphology has an ancillary effect of strengthening the face during all loading regimes. These results provide support for the functional hypotheses. However, we found that an anteriorly positioned zygomatic root increases strain magnitudes even in models with an inflated/reinforced root complex. These results suggest that an anteriorly placed zygomatic root complex evolved to enhance the efficiency of bite force production while facial reinforcement features, such as the anterior pillar and the expanded zygomatic root, may have been selected for in part to compensate for the weakening effect of this facial configuration. Anat Rec, 300:171-195, 2017. © 2016 Wiley Periodicals, Inc.

In vivo bone strain and finite element modeling of a rhesus macaque mandible during mastication

Finite element analysis (FEA) is a commonly used tool in musculoskeletal biomechanics and vertebrate paleontology. The accuracy and precision of finite element models (FEMs) are reliant on accurate data on bone geometry, muscle forces, boundary conditions and tissue material properties. Simplified modeling assumptions, due to lack of in vivo experimental data on material properties and muscle activation patterns, may introduce analytical errors in analyses where quantitative accuracy is critical for obtaining rigorous results. A subject-specific FEM of a rhesus macaque mandible was constructed, loaded and validated using in vivo data from the same animal. In developing the model, we assessed the impact on model behavior of variation in (i) material properties of the mandibular trabecular bone tissue and teeth; (ii) constraints at the temporomandibular joint and bite point; and (iii) the timing of the muscle activity used to estimate the external forces acting on the model. The best match between the FEA simulation and the in vivo experimental data resulted from modeling the trabecular tissue with an isotropic and homogeneous Young's modulus and Poisson's value of 10GPa and 0.3, respectively; constraining translations along X,Y, Z axes in the chewing (left) side temporomandibular joint, the premolars and the m₁; constraining the balancing (right) side temporomandibular joint in the anterior-posterior and superior-inferior axes, and using the muscle force estimated at time of maximum strain magnitude in the lower lateral gauge. The relative strain magnitudes in this model were similar to those recorded in vivo for all strain locations. More detailed analyses of mandibular strain patterns during the power stroke at different times in the chewing cycle are needed.

The jaw is a second-class lever in Pedetes capensis (Rodentia: Pedetidae)

The mammalian jaw is often modelled as a third-class lever for the purposes of biomechanical analyses, owing to the position of the resultant muscle force between the jaw joint and the teeth. However, it has been proposed that in some rodents the jaws operate as a second-class lever during distal molar bites, owing to the rostral position of the masticatory musculature. In particular, the infraorbital portion of the zygomatico-mandibularis (IOZM) has been suggested to be of major importance in converting the masticatory system from a third-class to a second-class lever. The presence of the IOZM is diagnostic of the hystricomorph rodents, and is particularly well-developed in Pedetes capensis, the South African springhare. In this study, finite element analysis (FEA) was used to assess the lever mechanics of the springhare masticatory system, and to determine the function of the IOZM. An FE model of the skull of P. capensis was constructed and loaded with all masticatory muscles, and then solved for biting at each tooth in turn. Further load cases were created in which each masticatory muscle was removed in turn. The analyses showed that the mechanical advantage of the springhare jaws was above one at all molar bites and very close to one during the premolar bite. Removing the IOZM or masseter caused a drop in mechanical advantage at all bites, but affected strain patterns and cranial deformation very little. Removing the ZM had only a small effect on mechanical advantage, but produced a substantial reduction in strain and deformation across the skull. It was concluded that the masticatory system of P. capensis acts as a second class lever during bites along almost the entire cheek tooth row. The IOZM is clearly a major contributor to this effect, but the masseter also has a part to play. The benefit of the IOZM is that it adds force without substantially contributing to strain or deformation of the skull. This may help explain why the hystricomorphous morphology has evolved multiple times independently within Rodentia.

Elastic Properties of Chimpanzee Craniofacial Cortical Bone

Relatively few assessments of cranial biomechanics formally take into account variation in the material properties of cranial cortical bone. Our aim was to characterize the elastic properties of chimpanzee craniofacial cortical bone and compare these to the elastic properties of dentate human craniofacial cortical bone. From seven cranial regions, 27 cylindrical samples were harvested from each of five chimpanzee crania. Assuming orthotropy, axes of maximum stiffness in the plane of the cortical plate were derived using modified equations of Hooke's law in a Mathcad program. Consistent orientations among individuals were observed in the zygomatic arch and alveolus. The density of cortical bone showed significant regional variation (P < 0.001). The elastic moduli demonstrated significant differences between sites, and a distinct pattern where E₃  > E₂  > E₁ . Shear moduli were significantly different among regions (P < 0.001). The pattern by which chimpanzee cranial cortical bone varies in elastic properties resembled that seen in humans, perhaps suggesting that the elastic properties of craniofacial bone in fossil hominins can be estimated with at least some degree of confidence. Anat Rec, 299:1718-1733, 2016. © 2016 Wiley Periodicals, Inc.

Review of In Vivo Bone Strain Studies and Finite Element Models of the Zygomatic Complex in Humans and Nonhuman Primates: Implications for Clinical Research and Practice

The craniofacial skeleton is often described in the clinical literature as being comprised of vertical bony pillars, which transmit forces from the toothrow to the neurocranium as axial compressive stresses, reinforced transversely by buttresses. Here, we review the literature on bony microarchitecture, in vivo bone strain, and finite-element modeling of the facial skeleton of humans and nonhuman primates to address questions regarding the structural and functional existence of facial pillars and buttresses. Available bone material properties data do not support the existence of pillars and buttresses in humans or Sapajus apella. Deformation regimes in the zygomatic complex emphasize bending and shear, therefore conceptualizing the zygomatic complex of humans or nonhuman primates as a pillar obscures its patterns of stress, strain, and deformation. Human fossil relatives and chimpanzees exhibit strain regimes corroborating the existence of a canine-frontal pillar, but the notion of a zygomatic pillar has no support. The emerging consensus on patterns of strain and deformation in finite element models (FEMs) of the human facial skeleton corroborates hypotheses in the clinical literature regarding zygomatic complex function, and provide new insights into patterns of failure of titanium and resorbable plates in experimental studies. It is suggested that the "pillar and buttress" model of human craniofacial skeleton function be replaced with FEMs that more accurately and precisely represent in vivo function, and which can serve as the basis for future research into implants used in restoration of occlusal function and fracture repair. Anat Rec, 299:1753-1778, 2016. © 2016 Wiley Periodicals, Inc.

Cell Mechanics of Craniosynostosis

Craniosynostosis is the premature fusion of the calvarial sutures that is associated with a number of physical and intellectual disabilities spanning from pediatric to adult years. Over the past two decades, techniques in molecular genetics and more recently, advances in high-throughput DNA sequencing have been used to examine the underlying pathogenesis of this disease. To date, mutations in 57 genes have been identified as causing craniosynostosis and the number of newly discovered genes is growing rapidly as a result of the advances in genomic technologies. While contributions from both genetic and environmental factors in this disease are increasingly apparent, there remains a gap in knowledge that bridges the clinical characteristics and genetic markers of craniosynostosis with their signaling pathways and mechanotransduction processes. By linking genotype to phenotype, outlining the role of cell mechanics may further uncover the specific mechanotransduction pathways underlying craniosynostosis. Here, we present a brief overview of the recent findings in craniofacial genetics and cell mechanics, discussing how this information together with animal models is advancing our understanding of craniofacial development.

Human feeding biomechanics: performance, variation, and functional constraints

The evolution of the modern human (Homo sapiens) cranium is characterized by a reduction in the size of the feeding system, including reductions in the size of the facial skeleton, postcanine teeth, and the muscles involved in biting and chewing. The conventional view hypothesizes that gracilization of the human feeding system is related to a shift toward eating foods that were less mechanically challenging to consume and/or foods that were processed using tools before being ingested. This hypothesis predicts that human feeding systems should not be well-configured to produce forceful bites and that the cranium should be structurally weak. An alternate hypothesis, based on the observation that humans have mechanically efficient jaw adductors, states that the modern human face is adapted to generate and withstand high biting forces. We used finite element analysis (FEA) to test two opposing mechanical hypotheses: that compared to our closest living relative, chimpanzees (Pan troglodytes), the modern human craniofacial skeleton is (1) less well configured, or (2) better configured to generate and withstand high magnitude bite forces. We considered intraspecific variation in our examination of human feeding biomechanics by examining a sample of geographically diverse crania that differed notably in shape. We found that our biomechanical models of human crania had broadly similar mechanical behavior despite their shape variation and were, on average, less structurally stiff than the crania of chimpanzees during unilateral biting when loaded with physiologically-scaled muscle loads. Our results also show that modern humans are efficient producers of bite force, consistent with previous analyses. However, highly tensile reaction forces were generated at the working (biting) side jaw joint during unilateral molar bites in which the chewing muscles were recruited with bilateral symmetry. In life, such a configuration would have increased the risk of joint dislocation and constrained the maximum recruitment levels of the masticatory muscles on the balancing (non-biting) side of the head. Our results do not necessarily conflict with the hypothesis that anterior tooth (incisors, canines, premolars) biting could have been selectively important in humans, although the reduced size of the premolars in humans has been shown to increase the risk of tooth crown fracture. We interpret our results to suggest that human craniofacial evolution was probably not driven by selection for high magnitude unilateral biting, and that increased masticatory muscle efficiency in humans is likely to be a secondary byproduct of selection for some function unrelated to forceful biting behaviors. These results are consistent with the hypothesis that a shift to softer foods and/or the innovation of pre-oral food processing techniques relaxed selective pressures maintaining craniofacial features that favor forceful biting and chewing behaviors, leading to the characteristically small and gracile faces of modern humans.

The Effect of Varying Jaw-elevator Muscle Forces on a Finite Element Model of a Human Cranium

Finite element analyses simulating masticatory system loading are increasingly undertaken in primates, hominin fossils and modern humans. Simplifications of models and loadcases are often required given the limits of data and technology. One such area of uncertainty concerns the forces applied to cranial models and their sensitivity to variations in these forces. We assessed the effect of varying force magnitudes among jaw-elevator muscles applied to a finite element model of a human cranium. The model was loaded to simulate incisor and molar bites using different combinations of muscle forces. Symmetric, asymmetric, homogeneous, and heterogeneous muscle activations were simulated by scaling maximal forces. The effects were compared with respect to strain distribution (i.e., modes of deformation) and magnitudes; bite forces and temporomandibular joint (TMJ) reaction forces. Predicted modes of deformation, strain magnitudes and bite forces were directly proportional to total applied muscle force and relatively insensitive to the degree of heterogeneity of muscle activation. However, TMJ reaction forces and mandibular fossa strains decrease and increase on the balancing and working sides according to the degree of asymmetry of loading. These results indicate that when modes, rather than magnitudes, of facial deformation are of interest, errors in applied muscle forces have limited effects. However the degree of asymmetric loading does impact on TMJ reaction forces and mandibular fossa strains. These findings are of particular interest in relation to studies of skeletal and fossil material, where muscle data are not available and estimation of muscle forces from skeletal proxies is prone to error. Anat Rec, 299:828-839, 2016. © 2016 Wiley Periodicals, Inc.

Dental materials for cleft palate repair

Numerous bone and soft tissue grafting techniques are followed to repair cleft of lip and palate (CLP) defects. In addition to the gold standard surgical interventions involving the use of autogenous grafts, various allogenic and xenogenic graft materials are available for bone regeneration. In an attempt to discover minimally invasive and cost effective treatments for cleft repair, an exceptional growth in synthetic biomedical graft materials have occurred. This study gives an overview of the use of dental materials to repair cleft of lip and palate (CLP). The eligibility criteria for this review were case studies, clinical trials and retrospective studies on the use of various types of dental materials in surgical repair of cleft palate defects. Any data available on the surgical interventions to repair alveolar or palatal cleft, with natural or synthetic graft materials was included in this review. Those datasets with long term clinical follow-up results were referred to as particularly relevant. The results provide encouraging evidence in favor of dental and other related biomedical materials to fill the gaps in clefts of lip and palate. The review presents the various bones and soft tissue replacement strategies currently used, tested or explored for the repair of cleft defects. There was little available data on the use of synthetic materials in cleft repair which was a limitation of this study. In conclusion although clinical trials on the use of synthetic materials are currently underway the uses of autologous implants are the preferred treatment methods to date.

Validity and sensitivity of a human cranial finite element model: implications for comparative studies of biting performance

Finite element analysis (FEA) is a modelling technique increasingly used in anatomical studies investigating skeletal form and function. In the case of the cranium this approach has been applied to both living and fossil taxa to (for example) investigate how form relates to function or infer diet or behaviour. However, FE models of complex musculoskeletal structures always rely on simplified representations because it is impossible completely to image and represent every detail of skeletal morphology, variations in material properties and the complexities of loading at all spatial and temporal scales. The effects of necessary simplifications merit investigation. To this end, this study focuses on one aspect, model geometry, which is particularly pertinent to fossil material where taphonomic processes often destroy the finer details of anatomy or in models built from clinical CTs where the resolution is limited and anatomical details are lost. We manipulated the details of a finite element (FE) model of an adult human male cranium and examined the impact on model performance. First, using digital speckle interferometry, we directly measured strains from the infraorbital region and frontal process of the maxilla of the physical cranium under simplified loading conditions, simulating incisor biting. These measured strains were then compared with predicted values from FE models with simplified geometries that included modifications to model resolution, and how cancellous bone and the thin bones of the circum-nasal and maxillary regions were represented. Distributions of regions of relatively high and low principal strains and principal strain vector magnitudes and directions, predicted by the most detailed FE model, are generally similar to those achieved in vitro. Representing cancellous bone as solid cortical bone lowers strain magnitudes substantially but the mode of deformation of the FE model is relatively constant. In contrast, omitting thin plates of bone in the circum-nasal region affects both mode and magnitude of deformation. Our findings provide a useful frame of reference with regard to the effects of simplifications on the performance of FE models of the cranium and call for caution in the interpretation and comparison of FEA results.

The relationship between cranial structure, biomechanical performance and ecological diversity in varanoid lizards

Skull structure is intimately associated with feeding ability in vertebrates, both in terms of specific performance measures and general ecological characteristics. This study quantitatively assessed variation in the shape of the cranium and mandible in varanoid lizards, and its relationship to structural performance (von Mises strain) and interspecific differences in feeding ecology. Geometric morphometric and linear morphometric analyses were used to evaluate morphological differences, and finite element analysis was used to quantify variation in structural performance (strain during simulated biting, shaking and pulling). This data was then integrated with ecological classes compiled from relevant scientific literature on each species in order to establish structure-function relationships. Finite element modelling results showed that variation in cranial morphology resulted in large differences in the magnitudes and locations of strain in biting, shaking and pulling load cases. Gracile species such as Varanus salvadorii displayed high strain levels during shaking, especially in the areas between the orbits. All models exhibit less strain during pull back loading compared to shake loading, even though a larger force was applied (pull =30N, shake = 20N). Relationships were identified between the morphology, performance, and ecology. Species that did not feed on hard prey clustered in the gracile region of cranial morphospace and exhibited significantly higher levels of strain during biting (P = 0.0106). Species that fed on large prey clustered in the elongate area of mandible morphospace. This relationship differs from those that have been identified in other taxonomic groups such as crocodiles and mammals. This difference may be due to a combination of the open 'space-frame' structure of the varanoid lizard skull, and the 'pull back' behaviour that some species use for processing large prey.

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.

Descriptive anatomy and three-dimensional reconstruction of the skull of the early tetrapod Acanthostega gunnari Jarvik, 1952

The early tetrapod Acanthostega gunnari is an iconic fossil taxon exhibiting skeletal morphology reflecting the transition of vertebrates from water onto land. Computed tomography data of two Acanthostega skulls was segmented using visualization software to digitally separate bone from matrix and individual bones of the skull from each other. A revised description of cranial and lower jaw anatomy in this taxon based on CT data includes new details of sutural morphology, the previously undescribed quadrate and articular bones, and the mandibular symphysis. Sutural morphology is used to infer loading regime in the skull during feeding, and suggests Acanthostega used its anterior jaws to initially seize prey while smaller posterior teeth were used to restrain struggling prey during ingestion. Novel methods were used to repair and retrodeform the skull, resulting in a three-dimensional digital reconstruction that features a longer postorbital region and more strongly hooked anterior lower jaw than previous attempts while supporting the presence of a midline gap between the nasals and median rostrals.

What does feeding system morphology tell us about feeding?

Feeding is the set of behaviors whereby organisms acquire and process the energy required for survival and reproduction. Thus, feeding system morphology is presumably subject to selection to maintain or improve feeding performance. Relationships among feeding system morphology, feeding behavior, and diet not only explain the morphological diversity of extant primates, but can also be used to reconstruct feeding behavior and diet in fossil taxa. Dental morphology has long been known to reflect aspects of feeding behavior and diet but strong relationships of craniomandibular morphology to feeding behavior and diet have yet to be defined.

The feeding biomechanics and dietary ecology of Paranthropus boisei

The African Plio-Pleistocene hominins known as australopiths evolved derived craniodental features frequently interpreted as adaptations for feeding on either hard, or compliant/tough foods. Among australopiths, Paranthropus boisei is the most robust form, exhibiting traits traditionally hypothesized to produce high bite forces efficiently and strengthen the face against feeding stresses. However, recent mechanical analyses imply that P. boisei may not have been an efficient producer of bite force and that robust morphology in primates is not necessarily strong. Here we use an engineering method, finite element analysis, to show that the facial skeleton of P. boisei is structurally strong, exhibits a strain pattern different from that in chimpanzees (Pan troglodytes) and Australopithecus africanus, and efficiently produces high bite force. It has been suggested that P. boisei consumed a diet of compliant/tough foods like grass blades and sedge pith. However, the blunt occlusal topography of this and other species suggests that australopiths are adapted to consume hard foods, perhaps including grass and sedge seeds. A consideration of evolutionary trends in morphology relating to feeding mechanics suggests that food processing behaviors in gracile australopiths evidently were disrupted by environmental change, perhaps contributing to the eventual evolution of Homo and Paranthropus.

Biomechanical implications of intraspecific shape variation in chimpanzee crania: moving toward an integration of geometric morphometrics and finite element analysis

In a broad range of evolutionary studies, an understanding of intraspecific variation is needed in order to contextualize and interpret the meaning of variation between species. However, mechanical analyses of primate crania using experimental or modeling methods typically encounter logistical constraints that force them to rely on data gathered from only one or a few individuals. This results in a lack of knowledge concerning the mechanical significance of intraspecific shape variation that limits our ability to infer the significance of interspecific differences. This study uses geometric morphometric methods (GM) and finite element analysis (FEA) to examine the biomechanical implications of shape variation in chimpanzee crania, thereby providing a comparative context in which to interpret shape-related mechanical variation between hominin species. Six finite element models (FEMs) of chimpanzee crania were constructed from CT scans following shape-space Principal Component Analysis (PCA) of a matrix of 709 Procrustes coordinates (digitized onto 21 specimens) to identify the individuals at the extremes of the first three principal components. The FEMs were assigned the material properties of bone and were loaded and constrained to simulate maximal bites on the P(3) and M(2) . Resulting strains indicate that intraspecific cranial variation in morphology is associated with quantitatively high levels of variation in strain magnitudes, but qualitatively little variation in the distribution of strain concentrations. Thus, interspecific comparisons should include considerations of the spatial patterning of strains rather than focus only on their magnitudes.