Modeling masticatory muscle force in finite element analysis: Sensitivity analysis using principal coordinates analysis

Bite hard: Linking cranial loading mechanics to ecological differences in gnawing behavior in caviomorph rodents

The mammalian skull is very malleable and has notably radiated into highly diverse morphologies, fulfilling a broad range of functional needs. Although gnawing is relatively common in mammals, this behavior and its associated morphology are diagnostic features for rodents. These animals possess a very versatile and highly mechanically advantageous masticatory apparatus, which, for instance, allowed caviomorph rodents to colonize South America during the Mid-Eocene and successfully radiate in over 200 extant species throughout most continental niches. Previous work has shown that differences in bite force within caviomorphs could be better explained by changes in muscle development than in mechanical advantages (i.e., in cranial overall morphology). Considering the strong bites they apply, it is interesting to assess how the reaction forces upon the incisors (compression) and the powerful adductor musculature pulling (tension) mechanically affect the cranium, especially between species with different ecologies (e.g., chisel-tooth digging). Thus, we ran finite element analyses upon crania of the subterranean Talas' tuco-tuco Ctenomys talarum, the semi-fossorial common degu Octodon degus, and the saxicolous long-tailed chinchilla Chinchilla lanigera to simulate: (A) in vivo biting in all species, and (B) rescaled muscle forces in non-ctenomyid rodents to match those of the tuco-tuco. Results show that the stress patterns correlate with the mechanical demands of distinctive ecologies, on in vivo-based simulations, with the subterranean tuco-tuco being the most stressed species. In contrast, when standardizing all three species (rescaled models), non-ctenomyid models exhibited a several-fold increase in stress, in both magnitude and affected areas. Detailed observations evidenced that this increase in stress was higher in lateral sections of the snout and, mainly, the zygomatic arch; between approximately 2.5-3.5 times in the common degu and 4.0-5.0 times in the long-tailed chinchilla. Yet, neither species, module, nor simulation condition presented load factor levels that would imply structural failure by strong, incidental biting. Our results let us conclude that caviomorphs have a high baseline for mechanical strength of the cranium because of the inheritance of a very robust "rodent" model, while interspecific differences are associated with particular masticatory habits and the concomitant level of development of the adductor musculature. Especially, the masseteric and zygomaticomandibular muscles contribute to >80% of the bite force, and therefore, their contraction is responsible for the highest strains upon their origin sites, that is, the zygomatic arch and the snout. Thus, the robust crania of the subterranean and highly aggressive tuco-tucos allow them to withstand much stronger forces than degus or chinchillas, such as the ones produced by their hypertrophied jaw adductor muscles or imparted by the soil reaction.

Functional adaptation of the infant craniofacial system to mechanical loadings arising from masticatory forces

The morphology and biomechanics of infant crania undergo significant changes between the pre- and post-weaning phases due to increasing loading of the masticatory system. The aims of this study were to characterize the changes in muscle forces, bite forces and the pattern of mechanical strain and stress arising from the aforementioned forces across crania in the first 48 months of life using imaging and finite element methods. A total of 51 head computed tomography scans of normal individuals were collected and analysed from a larger database of 217 individuals. The estimated mean muscle forces of temporalis, masseter and medial pterygoid increase from 30.9 to 87.0 N, 25.6 to 69.6 N and 23.1 to 58.9 N, respectively (0-48 months). Maximum bite force increases from 90.5 to 184.2 N (3-48 months). There is a change in the pattern of strain and stress from the calvaria to the face during postnatal development. Overall, this study highlights the changes in the mechanics of the craniofacial system during normal development. It further raises questions as to how and what level of changes in the mechanical forces during the development can alter the morphology of the craniofacial system.

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'.

Dietary niche partitioning in Early Jurassic ichthyosaurs from Strawberry Bank

Jurassic ichthyosaurs dominated upper trophic levels of marine ecosystems. Many species coexisted alongside each another, and it is uncertain whether they competed for the same array of food or divided dietary resources, each specializing in different kinds of prey. Here, we test whether feeding differences existed between species, applying finite element analysis to ichthyosaurs for the first time. We examine two juvenile ichthyosaur specimens, referred to Hauffiopteryx typicus and Stenopterygius triscissus, from the Strawberry Bank Lagerstätte, a shallow marine environment from the Early Jurassic of southern England (Toarcian, ~183 Ma). Snout and cranial robusticity differ between the species, with S. triscissus having a more robust snout and cranium and specializing in slow biting of hard prey, and H. typicus with its slender snout specializing in fast, but weaker bites on fast-moving, but soft prey. The two species did not differ in muscle forces, but stress distributions varied in the nasal area, reflecting differences when biting at different points along the tooth row: the more robustly snouted Stenopterygius resisted increases or shifts in stress distribution when the bite point was shifted from the posterior to the mid-point of the tooth row, but the slender-snouted Hauffiopteryx showed shifts and increases in stress distributions between these two bite points. The differences in cranial morphology, dentition and inferred stresses between the two species suggest adaptations for dietary niche partitioning.

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.

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.

Population affinity and variation of sexual dimorphism in three-dimensional facial forms: comparisons between Turkish and Japanese populations

Examining the extent to which sex differences in three-dimensional (3D) facial soft tissue configurations are similar across diverse populations could suggest the source of the indirect evolutionary benefits of facial sexual dimorphism traits. To explore this idea, we selected two geographically distinct populations. Three-dimensional model faces were derived from 272 Turkish and Japanese men and women; their facial morphologies were evaluated using landmark and surface-based analyses. We found four common facial features related to sexual dimorphism. Both Turkish and Japanese females had a shorter lower face height, a flatter forehead, greater sagittal cheek protrusion in the infraorbital region but less prominence of the cheek in the parotid-masseteric region, and an antero-posteriorly smaller nose when compared with their male counterparts. The results indicated the possible phylogenetic contribution of the masticatory organ function and morphogenesis on sexual dimorphism of the human face in addition to previously reported biological and psychological characteristics, including sexual maturity, reproductive potential, mating success, general health, immune response, age, and personality.

Back to the bones: do muscle area assessment techniques predict functional evolution across a macroevolutionary radiation?

Measures of attachment or accommodation area on the skeleton are a popular means of rapidly generating estimates of muscle proportions and functional performance for use in large-scale macroevolutionary studies. Herein, we provide the first evaluation of the accuracy of these muscle area assessment (MAA) techniques for estimating muscle proportions, force outputs and bone loading in a comparative macroevolutionary context using the rodent masticatory system as a case study. We find that MAA approaches perform poorly, yielding large absolute errors in muscle properties, bite force and particularly bone stress. Perhaps more fundamentally, these methods regularly fail to correctly capture many qualitative differences between rodent morphotypes, particularly in stress patterns in finite-element models. Our findings cast doubts on the validity of these approaches as means to provide input data for biomechanical models applied to understand functional transitions in the fossil record, and perhaps even in taxon-rich statistical models that examine broad-scale macroevolutionary patterns. We suggest that future work should go back to the bones to test if correlations between attachment area and muscle size within homologous muscles across a large number of species yield strong predictive relationships that could be used to deliver more accurate predictions for macroevolutionary and functional studies.

Experimental validation of finite element simulation of a new custom-designed fixation plate to treat mandibular angle fracture

BACKGROUND

The objective of the study was to validate biomechanical characteristics of a 3D-printed, novel-designated fixation plate for treating mandibular angle fracture, and compare it with two commonly used fixation plates by finite element (FE) simulations and experimental testing.

METHODS

A 3D virtual mandible was created from a patient's CT images as the master model. A custom-designed plate and two commonly used fixation plates were reconstructed onto the master model for FE simulations. Modeling of angle fracture, simulation of muscles of mastication, and defining of boundary conditions were integrated into the theoretical model. Strain levels during different loading conditions were analyzed using a finite element method (FEM). For mechanical test design, samples of the virtual mandible with angle fracture and the custom-designed fixation plates were printed using selective laser sintering (SLS) and selective laser melting (SLM) printing methods. Experimental data were collected from a testing platform with attached strain gauges to the mandible and the plates at different 10 locations during mechanical tests. Simulation of muscle forces and temporomandibular joint conditions were built into the physical models to improve the accuracy of clinical conditions. The experimental vs the theoretical data collected at the 10 locations were compared, and the correlation coefficient was calculated.

RESULTS

The results show that use of the novel-designated fixation plate has significant mechanical advantages compared to the two commonly used fixation plates. The results of measured strains at each location show a very high correlation between the physical model and the virtual mandible of their biomechanical behaviors under simulated occlusal loading conditions when treating angle fracture of the mandible.

CONCLUSIONS

Based on the results from our study, we validate the accuracy of our computational model which allows us to use it for future clinical applications under more sophisticated biomechanical simulations and testing.

Bite force and its relationship to jaw shape in domestic dogs

Previous studies based on two-dimensional methods have suggested that the great morphological variability of cranial shape in domestic dogs has impacted bite performance. Here, we used a three-dimensional biomechanical model based on dissection data to estimate the bite force of 47 dogs of various breeds at several bite points and gape angles. In vivo bite force for three Belgian shepherd dogs was used to validate our model. We then used three-dimensional geometric morphometrics to investigate the drivers of bite force variation and to describe the relationships between the overall shape of the jaws and bite force. The model output shows that bite force is rather variable in dogs and that dogs bite harder on the molar teeth and at lower gape angles. Half of the bite force is determined by the temporal muscle. Bite force also increased with size, and brachycephalic dogs showed higher bite forces for their size than mesocephalic dogs. We obtained significant covariation between the shape of the upper or lower jaw and absolute or residual bite force. Our results demonstrate that domestication has not resulted in a disruption of the functional links in the jaw system in dogs and that mandible shape is a good predictor of bite force.

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.

Principal coordinate analysis assisted chromatographic analysis of bacterial cell wall collection: A robust classification approach

In the present work, Principal coordinate analysis (PCoA) is introduced to develop a robust model to classify the chromatographic data sets of peptidoglycan sample. PcoA captures the heterogeneity present in the data sets by using the dissimilarity matrix as input. Thus, in principle, it can even capture the subtle differences in the bacterial peptidoglycan composition and can provide a more robust and fast approach for classifying the bacterial collection and identifying the novel cell wall targets for further biological and clinical studies. The utility of the proposed approach is successfully demonstrated by analysing the two different kind of bacterial collections. The first set comprised of peptidoglycan sample belonging to different subclasses of Alphaproteobacteria. Whereas, the second set that is relatively more intricate for the chemometric analysis consist of different wild type Vibrio Cholerae and its mutants having subtle differences in their peptidoglycan composition. The present work clearly proposes a useful approach that can classify the chromatographic data sets of chromatographic peptidoglycan samples having subtle differences. Furthermore, present work clearly suggest that PCoA can be a method of choice in any data analysis workflow.

Predicting calvarial growth in normal and craniosynostotic mice using a computational approach

During postnatal calvarial growth the brain grows gradually and the overlying bones and sutures accommodate that growth until the later juvenile stages. The whole process is coordinated through a complex series of biological, chemical and perhaps mechanical signals between various elements of the craniofacial system. The aim of this study was to investigate to what extent a computational model can accurately predict the calvarial growth in wild-type (WT) and mutant type (MT) Fgfr2C342Y/+ mice displaying bicoronal suture fusion. A series of morphological studies were carried out to quantify the calvarial growth at P3, P10 and P20 in both mouse types. MicroCT images of a P3 specimen were used to develop a finite element model of skull growth to predict the calvarial shape of WT and MT mice at P10. Sensitivity tests were performed and the results compared with ex vivo P10 data. Although the models were sensitive to the choice of input parameters, they predicted the overall skull growth in the WT and MT mice. The models also captured the difference between the ex vivoWT and MT mice. This modelling approach has the potential to be translated to human skull growth and to enhance our understanding of the different reconstruction methods used to manage clinically the different forms of craniosynostosis, and in the long term possibly reduce the number of re-operations in children displaying this condition and thereby enhance their quality of life.

Merging cranial histology and 3D-computational biomechanics: a review of the feeding ecology of a Late Triassic temnospondyl amphibian

Finite Element Analysis (FEA) is a useful method for understanding form and function. However, modelling of fossil taxa invariably involves assumptions as a result of preservation-induced loss of information in the fossil record. To test the validity of predictions from FEA, given such assumptions, these results could be compared to independent lines of evidence for cranial mechanics. In the present study a new concept of using bone microstructure to predict stress distribution in the skull during feeding is put forward and a correlation between bone microstructure and results of computational biomechanics (FEA) is carried out. The bony framework is a product of biological optimisation; bone structure is created to meet local mechanical conditions. To test how well results from FEA correlate to cranial mechanics predicted from bone structure, the well-known temnospondyl Metoposaurus krasiejowensis was used as a model. A crucial issue to Temnospondyli is their feeding mode: did they suction feed or employ direct biting, or both? Metoposaurids have previously been characterised either as active hunters or passive bottom dwellers. In order to test the correlation between results from FEA and bone microstructure, two skulls of Metoposaurus were used, one modelled under FE analyses, while for the second one 17 dermal bone microstructure were analysed. Thus, for the first time, results predicting cranial mechanical behaviour using both methods are merged to understand the feeding strategy of Metoposaurus. Metoposaurus appears to have been an aquatic animal that exhibited a generalist feeding behaviour. This taxon may have used two foraging techniques in hunting; mainly bilateral biting and, to a lesser extent, lateral strikes. However, bone microstructure suggests that lateral biting was more frequent than suggested by Finite Element Analysis (FEA). One of the potential factors that determined its mode of life may have been water levels. During optimum water conditions, metoposaurids may have been more active ambush predators that were capable of lateral strikes of the head. The dry season required a less active mode of life when bilateral biting is particularly efficient. This, combined with their characteristically anteriorly positioned orbits, was optimal for ambush strategy. This ability to use alternative modes of food acquisition, independent of environmental conditions, might hold the key in explaining the very common occurrence of metoposaurids during the Late Triassic.

Dynamic Musculoskeletal Functional Morphology: Integrating diceCT and XROMM

The tradeoff between force and velocity in skeletal muscle is a fundamental constraint on vertebrate musculoskeletal design (form:function relationships). Understanding how and why different lineages address this biomechanical problem is an important goal of vertebrate musculoskeletal functional morphology. Our ability to answer questions about the different solutions to this tradeoff has been significantly improved by recent advances in techniques for quantifying musculoskeletal morphology and movement. Herein, we have three objectives: (1) review the morphological and physiological parameters that affect muscle function and how these parameters interact; (2) discuss the necessity of integrating morphological and physiological lines of evidence to understand muscle function and the new, high resolution imaging technologies that do so; and (3) present a method that integrates high spatiotemporal resolution motion capture (XROMM, including its corollary fluoromicrometry), high resolution soft tissue imaging (diceCT), and electromyography to study musculoskeletal dynamics in vivo. The method is demonstrated using a case study of in vivo primate hyolingual biomechanics during chewing and swallowing. A sensitivity analysis demonstrates that small deviations in reconstructed hyoid muscle attachment site location introduce an average error of 13.2% to in vivo muscle kinematics. The observed hyoid and muscle kinematics suggest that hyoid elevation is produced by multiple muscles and that fascicle rotation and tendon strain decouple fascicle strain from hyoid movement and whole muscle length. Lastly, we highlight current limitations of these techniques, some of which will likely soon be overcome through methodological improvements, and some of which are inherent. Anat Rec, 301:378-406, 2018. © 2018 Wiley Periodicals, Inc.

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.

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.

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.

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.

Jaw morphology and fighting forces in stag beetles

The jaws of different species of stag beetles show a large variety of shapes and sizes. The male jaws are used as weapons in fights, and they may exert a very forceful bite in some species. We investigated in 16 species whether and how their forcefulness is reflected in their jaw morphology. We found a large range of maximal muscle forces (1.8N-33N; factor 18). Species investing in large bite muscles, also have disproportionately large jaw volumes. They use this additional jaw volume to elongate their jaws, increasing their winning chances in battles. The fact that this also decreases the mechanical advantage, is largely compensated by elongated in-levers. As a result, high muscle forces are correlated with elevated bite forces (0.27N-7.6N; factor 28). Despite the large difference in forcefulness, all investigated species experience similar Von Mises stresses in their jaws while biting (29MPa–114MPa; factor 4.0; calculated with Finite Element simulations). Hence, stag beetles have successfully adapted their jaw anatomy according to their bite force in fights.

Validation experiments on finite element models of an ostrich (Struthio camelus) cranium

The first finite element (FE) validation of a complete avian cranium was performed on an extant palaeognath, the ostrich (Struthio camelus). Ex-vivo strains were collected from the cranial bone and rhamphotheca. These experimental strains were then compared to convergence tested, specimen-specific finite element (FE) models. The FE models contained segmented cortical and trabecular bone, sutures and the keratinous rhamphotheca as identified from micro-CT scan data. Each of these individual materials was assigned isotropic material properties either from the literature or from nanoindentation, and the FE models compared to the ex-vivo results. The FE models generally replicate the location of peak strains and reflect the correct mode of deformation in the rostral region. The models are too stiff in regions of experimentally recorded high strain and too elastic in regions of low experimentally recorded low strain. The mode of deformation in the low strain neurocranial region is not replicated by the FE models, and although the models replicate strain orientations to within 10° in some regions, in most regions the correlation is not strong. Cranial sutures, as has previously been found in other taxa, are important for modifying both strain magnitude and strain patterns across the entire skull, but especially between opposing the sutural junctions. Experimentally, we find that the strains on the surface of the rhamphotheca are much lower than those found on nearby bone. The FE models produce much higher principal strains despite similar strain ratios across the entirety of the rhamphotheca. This study emphasises the importance of attempting to validate FE models, modelling sutures and rhamphothecae in birds, and shows that whilst location of peak strain and patterns of deformation can be modelled, replicating experimental data in digital models of avian crania remains problematic.

The impact of simplifications on the performance of a finite element model of a Macaca fascicularis cranium

In recent years finite element analysis (FEA) has emerged as a useful tool for the analysis of skeletal form-function relationships. While this approach has obvious appeal for the study of fossil specimens, such material is often fragmentary with disrupted internal architecture and can contain matrix that leads to errors in accurate segmentation. Here we examine the effects of varying the detail of segmentation and material properties of teeth on the performance of a finite element model of a Macaca fascicularis cranium within a comparative functional framework. Cranial deformations were compared using strain maps to assess differences in strain contours and Procrustes size and shape analyses, from geometric morphometrics, were employed to compare large scale deformations. We show that a macaque model subjected to biting can be made solid, and teeth altered in material properties, with minimal impact on large scale modes of deformation. The models clustered tightly by bite point rather than by modeling simplification approach, and fell out as being distinct from another species. However localized fluctuations in predicted strain magnitudes were recorded with different modeling approaches, particularly over the alveolar region. This study indicates that, while any model simplification should be undertaken with care and attention to its effects, future applications of FEA to fossils with unknown internal architecture may produce reliable results with regard to general modes of deformation, even when detail of internal bone architecture cannot be reliably modeled.

Comparative finite element analysis of the cranial performance of four herbivorous marsupials

Marsupial herbivores exhibit a wide variety of skull shapes and sizes to exploit different ecological niches. Several studies on teeth, dentaries, and jaw adductor muscles indicate that marsupial herbivores exhibit different specializations for grazing and browsing. No studies, however, have examined the skulls of marsupial herbivores to determine the relationship between stress and strain, and the evolution of skull shape. The relationship between skull morphology, biomechanical performance, and diet was tested by applying the finite element method to the skulls of four marsupial herbivores: the common wombat (Vombatus ursinus), koala (Phascolarctos cinereus), swamp wallaby (Wallabia bicolor), and red kangaroo (Macropus rufus). It was hypothesized that grazers, requiring stronger skulls to process tougher food, would have higher biomechanical performance than browsers. This was true when comparing the koala and wallaby (browsers) to the wombat (a grazer). The cranial model of the wombat resulted in low stress and high mechanical efficiency in relation to a robust skull capable of generating high bite forces. However, the kangaroo, also a grazer, has evolved a very different strategy to process tough food. The cranium is much more gracile and has higher stress and lower mechanical efficiency, but they adopt a different method of processing food by having a curved tooth row to concentrate force in a smaller area and molar progression to remove worn teeth from the tooth row. Therefore, the position of the bite is crucial for the structural performance of the kangaroo skull, while it is not for the wombat which process food along the entire tooth row. In accordance with previous studies, the results from this study show the mammalian skull is optimized to resist forces generated during feeding. However, other factors, including the lifestyle of the animal and its environment, also affect selection for skull morphology to meet multiple functional demands.

Reconstruction of muscle fascicle architecture from iodine-enhanced microCT images: A combined texture mapping and streamline approach

Skeletal muscle models are used to investigate motion and force generation in both biological and bioengineering research. Yet, they often lack a realistic representation of the muscle's internal architecture which is primarily composed of muscle fibre bundles, known as fascicles. Recently, it has been shown that fascicles can be resolved with micro-computed tomography (µCT) following staining of the muscle tissue with iodine potassium iodide (I2KI). Here, we present the reconstruction of the fascicular spatial arrangement and geometry of the superficial masseter muscle of a dog based on a combination of pattern recognition and streamline computation. A cadaveric head of a dog was incubated in I2KI and µCT-scanned. Following segmentation of the masseter muscle a statistical pattern recognition algorithm was applied to create a vector field of fascicle directions. Streamlines were then used to transform the vector field into a realistic muscle fascicle representation. The lengths of the reconstructed fascicles and the pennation angles in two planes (frontal and sagittal) were extracted and compared against a tracked fascicle field obtained through cadaver dissection. Both fascicle lengths and angles were found to vary substantially within the muscle confirming the complex and heterogeneous nature of skeletal muscle described by previous studies. While there were significant differences in the pennation angle between the experimentally derived and µCT-reconstructed data, there was congruence in the fascicle lengths. We conclude that the presented approach allows for embedding realistic fascicle information into finite element models of skeletal muscles to better understand the functioning of the musculoskeletal system.

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.

Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull

Background

Which factors influence the distribution patterns of morphological diversity among clades? The adaptive radiation model predicts that a clade entering new ecological niche will experience high rates of evolution early in its history, followed by a gradual slowing. Here we measure disparity and rates of evolution in Carnivora, specifically focusing on the terrestrial-aquatic transition in Pinnipedia. We analyze fissiped (mostly terrestrial, arboreal, and semi-arboreal, but also including the semi-aquatic otter) and pinniped (secondarily aquatic) carnivorans as a case study of an extreme ecological transition. We used 3D geometric morphometrics to quantify cranial shape in 151 carnivoran specimens (64 fissiped, 87 pinniped) and five exceptionally-preserved fossil pinnipeds, including the stem-pinniped Enaliarctos emlongi. Range-based and variance-based disparity measures were compared between pinnipeds and fissipeds. To distinguish between evolutionary modes, a Brownian motion model was compared to selective regime shifts associated with the terrestrial-aquatic transition and at the base of Pinnipedia. Further, evolutionary patterns were estimated on individual branches using both Ornstein-Uhlenbeck and Independent Evolution models, to examine the origin of pinniped diversity.

Results

Pinnipeds exhibit greater cranial disparity than fissipeds, even though they are less taxonomically diverse and, as a clade nested within fissipeds, phylogenetically younger. Despite this, there is no increase in the rate of morphological evolution at the base of Pinnipedia, as would be predicted by an adaptive radiation model, and a Brownian motion model of evolution is supported. Instead basal pinnipeds populated new areas of morphospace via low to moderate rates of evolution in new directions, followed by later bursts within the crown-group, potentially associated with ecological diversification within the marine realm.

Conclusion

The transition to an aquatic habitat in carnivorans resulted in a shift in cranial morphology without an increase in rate in the stem lineage, contra to the adaptive radiation model. Instead these data suggest a release from evolutionary constraint model, followed by aquatic diversifications within crown families.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0285-5) contains supplementary material, which is available to authorized users.

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.

In vivo cranial bone strain and bite force in the agamid lizard Uromastyx geyri

In vivo bone strain data are the most direct evidence of deformation and strain regimes in the vertebrate cranium during feeding and can provide important insights into skull morphology. Strain data have been collected during feeding across a wide range of mammals; in contrast, in vivo cranial bone strain data have been collected from few sauropsid taxa. Here we present bone strain data recorded from the jugal of the herbivorous agamid lizard Uromastyx geyri along with simultaneously recorded bite force. Principal and shear strain magnitudes in Uromastyx geyri were lower than cranial bone strains recorded in Alligator mississippiensis, but higher than those reported from herbivorous mammals. Our results suggest that variations in principal strain orientations in the facial skeleton are largely due to differences in feeding behavior and bite location, whereas food type has little impact on strain orientations. Furthermore, mean principal strain orientations differ between male and female Uromastyx during feeding, potentially because of sexual dimorphism in skull morphology.

Beware the black box: investigating the sensitivity of FEA simulations to modelling factors in comparative biomechanics

Finite element analysis (FEA) is a computational technique of growing popularity in the field of comparative biomechanics, and is an easily accessible platform for form-function analyses of biological structures. However, its rapid evolution in recent years from a novel approach to common practice demands some scrutiny in regards to the validity of results and the appropriateness of assumptions inherent in setting up simulations. Both validation and sensitivity analyses remain unexplored in many comparative analyses, and assumptions considered to be ‘reasonable’ are often assumed to have little influence on the results and their interpretation.

Here we report an extensive sensitivity analysis where high resolution finite element (FE) models of mandibles from seven species of crocodile were analysed under loads typical for comparative analysis: biting, shaking, and twisting. Simulations explored the effect on both the absolute response and the interspecies pattern of results to variations in commonly used input parameters. Our sensitivity analysis focuses on assumptions relating to the selection of material properties (heterogeneous or homogeneous), scaling (standardising volume, surface area, or length), tooth position (front, mid, or back tooth engagement), and linear load case (type of loading for each feeding type).

Our findings show that in a comparative context, FE models are far less sensitive to the selection of material property values and scaling to either volume or surface area than they are to those assumptions relating to the functional aspects of the simulation, such as tooth position and linear load case. Results show a complex interaction between simulation assumptions, depending on the combination of assumptions and the overall shape of each specimen. Keeping assumptions consistent between models in an analysis does not ensure that results can be generalised beyond the specific set of assumptions used. Logically, different comparative datasets would also be sensitive to identical simulation assumptions; hence, modelling assumptions should undergo rigorous selection. The accuracy of input data is paramount, and simulations should focus on taking biological context into account. Ideally, validation of simulations should be addressed; however, where validation is impossible or unfeasible, sensitivity analyses should be performed to identify which assumptions have the greatest influence upon the results.

Bite of the cats: relationships between functional integration and mechanical performance as revealed by mandible geometry

Cat-like carnivorous mammals represent a relatively homogeneous group of species whose morphology appears constrained by exclusive adaptations for meat eating. We present the most comprehensive data set of extant and extinct cat-like species to test for evolutionary transformations in size, shape and mechanical performance, that is, von Mises stress and surface traction, of the mandible. Size and shape were both quantified by means of geometric morphometrics, whereas mechanical performance was assessed applying finite element models to 2D geometry of the mandible. Additionally, we present the first almost complete composite phylogeny of cat-like carnivorans for which well-preserved mandibles are known, including representatives of 35 extant and 59 extinct species of Felidae, Nimravidae, and Barbourofelidae. This phylogeny was used to test morphological differentiation, allometry, and covariation of mandible parts within and among clades. After taking phylogeny into account, we found that both allometry and mechanical variables exhibit a significant impact on mandible shape. We also tested whether mechanical performance was linked to morphological integration. Mechanical stress at the coronoid process is higher in sabertoothed cats than in any other clade. This is strongly related to the high degree of covariation within modules of sabertooths mandibles. We found significant correlation between integration at the clade level and per-clade averaged stress values, on both original data and by partialling out interclade allometry from shapes when calculating integration. This suggests a strong interaction between natural selection and the evolution of developmental and functional modules at the clade level.

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.