Free body analysis, beam mechanics, and finite element modeling of the mandible of Alligator mississippiensis

Revisiting Old Questions With New Methods: The Effect of Embryonic Motility on Skull Development in the Domestic Chick

Muscle loading is known to influence skeletal morphology. Therefore, modification of the biomechanical environment is expected to cause coordinated morphological changes to the bony and cartilaginous tissues. Understanding how this musculoskeletal coordination contributes to morphological variation has relevance to health sciences, developmental biology, and evolutionary biology. To investigate how muscle loading influences skeletal morphology, we replicate a classic in ovo embryology experiment in the domestic chick (Gallus gallus domesticus) while harnessing modern methodologies that allow us to quantify skeletal anatomy more precisely and in situ. We induced rigid muscle paralysis in developing chicks mid-incubation, then compared the morphology of the cranium and mandible between immobilized and untreated embryos using microcomputed tomography and landmark-based geometric morphometric methods. Like earlier studies, we found predictable differences in the size and shape of the cranium and mandible in paralyzed chicks. These differences were concentrated in areas known to experience high strains during feeding, including the jaw joint and jaw muscle attachment sites. These results highlight specific areas of the skull that appear to be mechanosensitive and suggest muscles that could produce the biomechanical stimuli necessary for normal hatchling morphology. Interestingly, these same areas correspond to areas that show the greatest disparity and fastest evolutionary rates across the avian diversity, which suggests that the musculoskeletal integration observed during development extends to macroevolutionary scales. Thus, selection and evolutionary changes to muscle physiology and architecture could generate large and predictable changes to skull morphology. Building upon previous work, the adoption of modern imaging and morphometric techniques allows richer characterization of musculoskeletal integration that empowers researchers to understand how tissue-to-tissue interactions contribute to overall phenotypic variation.

Morphological evolution and functional consequences of giantism in tyrannosauroid dinosaurs

Tyrannosauroids are a clade of theropod dinosaur taxa that varied greatly in their body size distribution. We investigated the feeding performance of six tyrannosaur genera of variable body size and skull morphology. We used 3D finite element analysis to test whether skull shape becomes more or less resistant to feeding-induced forces. Cranial and mandibular models were scaled by adult Tyrannosaurus's surface area to analyze the influence of shape on skull function. It was found that Tyrannosaurus experienced higher absolute stresses compared to small-bodied relatives. When surface area values were equalized across genera to account for the effect of size and test efficiency of skull shape, smaller individuals experience notably greater stresses than larger relatives due to the robust cranial osteology characterized in the allometry of tyrannosaurids. These results may indicate that the wide crania of tyrannosaurids convey a functional advantage that basal tyrannosauroids, juvenile tyrannosauroids, and alioramins lacked.

A new biomechanical model of the mammal jaw based on load path analysis

The primary function of the tetrapod jaw is to transmit jaw muscle forces to bite points. The routes of force transfer in the jaw have never been studied but can be quantified using load paths - the shortest, stiffest routes from regions of force application to support constraints. Here, we use load path analysis to map force transfer from muscle attachments to bite point and jaw joint, and to evaluate how different configurations of trabecular and cortical bone affect load paths. We created three models of the mandible of the Virginia opossum, Didelphis virginiana, each with a cortical bone shell, but with different material properties for the internal spaces: (1) a cortical-trabecular model, in which the interior space is modeled with bulk properties of trabecular bone; (2) a cortical-hollow model, in which trabeculae and mandibular canal are modeled as hollow; and (3) a solid-cortical model, in which the interior is modeled as cortical bone. The models were compared with published in vivo bite force and bone strain data, and the load paths calculated for each model. The trabecular model, which is preferred because it most closely approximates the actual morphology, was best validated by in vivo data. In all three models, the load path was confined to cortical bone, although its route within the cortex varied depending on the material properties of the inner model. Our analysis shows that most of the force is transferred through the cortical, rather than trabecular bone, and highlights the potential of load path analysis for understanding form-function relationships in the skeleton.

Comparative cranial biomechanics reveal that Late Cretaceous tyrannosaurids exerted relatively greater bite force than in early-diverging tyrannosauroids

Tyrannosaurus has been an exemplar organism in feeding biomechanical analyses. An adult Tyrannosaurus could exert a bone-splintering bite force, through expanded jaw muscles and a robust skull and teeth. While feeding function of adult Tyrannosaurus has been thoroughly studied, such analyses have yet to expand to other tyrannosauroids, especially early-diverging tyrannosauroids (Dilong, Proceratosaurus, and Yutyrannus). In our analysis, we broadly assessed the cranial and feeding performance of tyrannosauroids at varying body sizes. Our sample size included small (Proceratosaurus and Dilong), medium-sized (Teratophoneus), and large (Tarbosaurus, Daspletosaurus, Gorgosaurus, and Yutyrannus) tyrannosauroids, and incorporation of tyrannosaurines at different ontogenetic stages (small juvenile Tarbosaurus, Raptorex, and mid-sized juvenile Tyrannosaurus). We used jaw muscle force calculations and finite element analysis to comprehend the cranial performance of our tyrannosauroids. Scaled subtemporal fenestrae areas and calculated jaw muscle forces show that broad-skulled tyrannosaurines (Tyrannosaurus, Daspletosaurus, juvenile Tyrannosaurus, and Raptorex) exhibited higher jaw muscle forces than other similarly sized tyrannosauroids (Gorgosaurus, Yutyrannus, and Proceratosaurus). The large proceratosaurid Yutyrannus exhibited lower cranial stress than most adult tyrannosaurids. This suggests that cranial structural adaptations of large tyrannosaurids maintained adequate safety factors at greater bite force, but their robust crania did not notably decrease bone stress. Similarly, juvenile tyrannosaurines experienced greater cranial stress than similarly-sized earlier tyrannosauroids, consistent with greater adductor muscle forces in the juveniles, and with crania no more robust than in their small adult predecessors. As adult tyrannosauroid body size increased, so too did relative jaw muscle forces manifested even in juveniles of giant adults.

Assessing the adductor musculature and jaw mechanics of Proterochampsa nodosa (Archosauriformes: Proterochampsidae) through finite element analysis

Proterochampsids are a group of South American nonarchosaurian archosauromorphs whose general morphology has been historically likened to that of the extant Crocodylia, which purportedly exhibited similar habits by convergence. Taxa from the genus Proterochampsa, for example, show platyrostral skulls with dorsally faced orbits and external nares and elongated snouts that might indicate a feeding habit similar to that of crocodilians. Nonetheless, some aspects of their craniomandibular anatomy are distinct. Proterochampsa has comparatively larger skull temporal fenestrae, and a unique morphology of the mandibular adductor chamber, with a remarkably large surangular shelf and a fainter retroarticular region in the mandible. In light of this, we conducted biomechanical tests on a 3-dimensional model of Proterochampsa nodosa including the first Finite Element Analysis for proterochampsians and compared it with models of the extant crocodylians Tomistoma schlegelii and Alligator mississippiensis. Our analyses suggested that, despite the differences in adductor chamber, Proterochampsa was able to perform bite forces comparable to those modeled for Alligator and significantly higher than Tomistoma. However, the morphology of the surangular shelf and the adductor chamber of Proterochampsa renders it more prone to accumulate stresses resulting from muscle contraction, when compared with both analogs. The elongated lower jaw of Proterochampsa, like that of Tomistoma, is more susceptible to bending, when compared with Alligator. As a result, we suggest that Proterochampsa might employ anteriorly directed bites only when handling small and soft-bodied prey. In addition, Proterochampsa exemplifies the diversity of arrangements that the adductor musculature adopted in different diverging archosauromorph groups.

Regional variation of the cortical and trabecular bone material properties in the rabbit skull

The material properties of some bones are known to vary with anatomical location, orientation and position within the bone (e.g., cortical and trabecular bone). Details of the heterogeneity and anisotropy of bone is an important consideration for biomechanical studies that apply techniques such as finite element analysis, as the outcomes will be influenced by the choice of material properties used. Datasets detailing the regional variation of material properties in the bones of the skull are sparse, leaving many finite element analyses of skulls no choice but to employ homogeneous, isotropic material properties, often using data from a different species to the one under investigation. Due to the growing significance of investigating the cranial biomechanics of the rabbit in basic science and clinical research, this study used nanoindentation to measure the elastic modulus of cortical and trabecular bone throughout the skull. The elastic moduli of cortical bone measured in the mediolateral and ventrodorsal direction were found to decrease posteriorly through the skull, while it was evenly distributed when measured in the anteroposterior direction. Furthermore, statistical tests showed that the variation of elastic moduli between separate regions (anterior, middle and posterior) of the skull were significantly different in cortical bone, but was not in trabecular bone. Elastic moduli measured in different orthotropic planes were also significantly different, with the moduli measured in the mediolateral direction consistently lower than that measured in either the anteroposterior or ventrodorsal direction. These findings demonstrate the significance of regional and directional variation in cortical bone elastic modulus, and therefore material properties in finite element models of the skull, particularly those of the rabbit, should consider the heterogeneous and orthotropic properties of skull bone when possible.

Gape drives regional variation in temporalis architectural dynamics in tufted capuchins

Dynamic changes in jaw movements and bite forces depend on muscle architectural and neural factors that have rarely been compared within the same muscle. Here we investigate how regional muscle architecture dynamics-fascicle rotation, shortening, lengthening and architectural gear ratio (AGR)-vary during chewing across a functionally heterogeneous muscle. We evaluate whether timing in architecture dynamics relates to gape, food material properties and/or muscle activation. We also examine whether static estimates of temporalis fibre architecture track variation in dynamic architecture. Fascicle-level architecture dynamics were measured in three regions of the superficial temporalis of three adult tufted capuchins (Sapajus apella) using biplanar videoradiography and the XROMM workflow. Architecture dynamics data were paired with regional fine-wire electromyography data from four adult tufted capuchins. Gape accounted for most architectural change across the temporalis, but architectural dynamics varied between regions. Mechanically challenging foods were associated with lower AGRs in the anterior region. The timing of most dynamic architectural changes did not vary between regions and differed from regional variation in static architecture. Collectively these findings suggest that, when modelling temporalis muscle force production in extant and fossil primates, it is important to account for the effects of gape, regionalization and food material properties. This article is part of the theme issue 'Food processing and nutritional assimilation in animals'.

A three-dimensional geometric morphometric analysis of the morphological transformation of Caiman lower jaw during post-hatching ontogeny

Shape ontogenetic changes of the lower jaw in crocodylians are poorly understood. In order to answer some questions related to the inter- and intraspecific morphological variation of the mandible of two extant Caiman species, we performed a three-dimensional geometric morphometric approach. For this purpose, we used landmarks and semilandmarks on two ontogenetic mandibular series of 48 and 15 post-hatching specimens of C. yacare and C. latirostris, respectively. We have also examined the relationship between these anatomical transformations and ontogenetic shifts in diet. We performed a principal component analysis (PCA) for the two species, and regression and partial least squares (PLS) analyses for each species, separately. As a result, species were segregated along the PC1 with specimens of C. yacare showing more gracile mandibles, and specimens of C. latirostris more robust ones. The PC2 and regression analyses showed an age gradient and represented ontogenetic shape changes. Adult caiman mandibles are higher and wider than juvenile ones, and shape changes are more conspicuous in C. latirostris. The PLS analyses showed a significant relationship between shape and diet. Morphological changes of the PLS1 of block-1 match with those of the regression analysis for both species. We have detected morphological transformations in areas where the musculature in charge of mandibular movements is attached. Common morphological changes occurring during ontogeny seem to reflect the same mechanical properties required for crushing and killing in both species, driven by an ontogenetic shift in the diet from invertebrates to vertebrates. Additionally, interspecific differences were also found to be correlated to ontogenetic changes in diet and could be related to dissimilar feeding mechanical requirements (e.g., stiffness and toughness of the item consumed), and to different habitat preferences. Robust mandibles would be more suitable for shallow and fully vegetated environments, as it can be seen in C. latirostris, whereas slender jaws seem to be more suitable for more aquatic species such as C. yacare.

Humeri under external load: Mechanical implications of differing bone curvature in American otter and honey badger

The mechanical properties of limb long bones are impacted by bone shape and especially curvature, which is therefore likely to be of adaptive value. We use finite element analysis to compare the mechanical properties of humeri of the closely related American otter and honey badger under external loads, and to analyze the significance of bone curvature. We simulate the effects generated by loads applied in directions that differ relative to the humeral longitudinal axes, and then compare the stress characteristics with a series of humerus-inspired abstracted curved structures with increasing ratio (C/R) of eccentricity C to radius of cross section R. The humeri of the two species differ in bone curvature, with C/R of 0.6201 and 0.8752, respectively. Our analysis shows that the peak and mean stress values found within the sampling line of bone models reach a minimum when the directions of loads are 105 ± 5°, and the humerus of the American otter always experienced lower stress values than those of the honey badger in the sampling line. An analysis of stress distribution in abstract curved structures showed the greatest reduction in stress when the direction of external load was equal or greater than 95°. This suggests that the variability of the direction of external loads is an important determinant of bone curvature, and should be accounted for when assessing load carrying capacity. This study provides a basis for biomechanics research and yields insight into the form-function relationship of nature's structural elements within limbs. It potentially contributes to the design of biomimetic robots while also highlighting the functional significance of humeral bone curvature in mammals.

Multiple pathways to herbivory underpinned deep divergences in ornithischian evolution

The extent to which evolution is deterministic is a key question in biology,^{1,2,3,4,5,6,7,8,9} with intensive debate on how adaptation^{6,10,11,12,13} and constraints^14,15,16 might canalize solutions to ecological challenges.^4,5,6 Alternatively, unique adaptations^1,9,17 and phylogenetic contingency^1,3,18 may render evolution fundamentally unpredictable.³ Information from the fossil record is critical to this debate,^1,2,11 but performance data for extinct taxa are limited.⁷ This knowledge gap is significant, as general morphology may be a poor predictor of biomechanical performance.^17,19,20 High-fiber herbivory originated multiple times within ornithischian dinosaurs,²¹ making them an ideal clade for investigating evolutionary responses to similar ecological pressures.²² However, previous biomechanical modeling studies on ornithischian crania^17,23,24,25 have not compared early-diverging taxa spanning independent acquisitions of herbivory. Here, we perform finite-element analysis on the skull of five early-diverging members of the major ornithischian clades to characterize morphofunctional pathways to herbivory. Results reveal limited functional convergence among ornithischian clades, with each instead achieving comparable performance, in terms of reconstructed patterns and magnitudes of functionally induced stress, through different adaptations of the feeding apparatus. Thyreophorans compensated for plesiomorphic low performance through increased absolute size, heterodontosaurids expanded jaw adductor muscle volume, ornithopods increased jaw system efficiency, and ceratopsians combined these approaches. These distinct solutions to the challenges of herbivory within Ornithischia underpinned the success of this diverse clade. Furthermore, the resolution of multiple solutions to equivalent problems within a single clade through macroevolutionary time demonstrates that phenotypic evolution is not necessarily predictable, instead arising from the interplay of adaptation, innovation, contingency, and constraints.^{1,2,3,7,8,9,18}.

Growth constraints set an upper limit to theropod dinosaur body size

Despite nearly 200 years of scientific collecting and study, none of the extinct, bipedal, predatory, theropod dinosaurs have been reliably shown to exceed 12 m in length. Using digital 3D models of theropods with lengths spanning 80 cm to 12 m, their body masses were found to scale to the 3.5 power of body lengths. The lateral area of the pelvis and the cross-sectional area of the tail base of these animals corresponds to the cross-sectional areas of key muscle groups important for balance and locomotion, and both scale to the 2.4 power of body length. Body accelerations in the lateral and forward directions are, using F = ma, given by dividing muscle area (force proxy) by body mass. Plotting these acceleration estimates against body length shows them to decrease exponentially. The largest theropods with body lengths of 10-12 m have less than 10% of the acceleration capacity of the smaller forms. The distinct lack of fossil remains of theropods demonstrably longer than 12 m suggests that the theropod body plan had an upper size limit based on a minimum acceleration threshold. Rotational inertia of the theropod body was found to be proportional to body length raised to the 5.5 power, and with increasing length, the capacity for agility would rapidly diminish. The tight relationship between theropod pelvic area and body length allows for the estimation of body lengths of specimens lacking complete axial skeletons, and this is done for four, large, well-preserved pelves.

The ecological diversification and evolution of Teleosauroidea (Crocodylomorpha, Thalattosuchia), with insights into their mandibular biomechanics

Throughout the Jurassic, a plethora of marine reptiles dominated ocean waters, including ichthyosaurs, plesiosaurs and thalattosuchian crocodylomorphs. These Jurassic ecosystems were characterized by high niche partitioning and spatial variation in dietary ecology. However, while the ecological diversity of many marine reptile lineages is well known, the overall ecological diversification of Teleosauroidea (one of the two major groups within thalattosuchian crocodylomorphs) has never been explored. Teleosauroids were previously deemed to have a morphologically conservative body plan; however, they were in actuality morphofunctionally more diverse than previously thought. Here we investigate the ecology and feeding specializations of teleosauroids, using morphological and functional cranio-dental characteristics. We assembled the most comprehensive dataset to date of teleosauroid taxa (approximately 20 species) and ran a series of principal component analyses (PC) to categorize them into various feeding ecomorphotypes based on 17 dental characteristics (38 specimens) and 16 functionally significant mandibular characters (18 specimens). The results were examined in conjunction with a comprehensive thalattosuchian phylogeny (153 taxa and 502 characters) to evaluate macroevolutionary patterns and significant ecological shifts. Machimosaurids display a well-developed ecological shift from: (1) slender, pointed tooth apices and an elongate gracile mandible; to (2) more robust, pointed teeth with a slightly deeper mandible; and finally, (3) rounded teeth and a deep-set, shortened mandible with enlarged musculature. Overall, there is limited mandibular functional variability in teleosaurids and machimosaurids, despite differing cranial morphologies and habitat preferences in certain taxa. This suggests a narrow feeding ecological divide between teleosaurids and machimosaurids. Resource partitioning was primarily related to snout and skull length as well as habitat; only twice did teleosauroids manage to make a major evolutionary leap to feed distinctly differently, with only the derived machimosaurines successfully radiating into new feeding ecologies.

The efficacy of computed tomography scanning versus surface scanning in 3D finite element analysis

Finite element analysis (FEA) is a commonly used application in biomechanical studies of both extant and fossil taxa to assess stress and strain in solid structures such as bone. FEA can be performed on 3D structures that are generated using various methods, including computed tomography (CT) scans and surface scans. While previous palaeobiological studies have used both CT scanned models and surface scanned models, little research has evaluated to what degree FE results may vary when CT scans and surface scans of the same object are compared. Surface scans do not preserve the internal geometries of 3D structures, which are typically preserved in CT scans. Here, we created 3D models from CT scans and surface scans of the same specimens (crania and mandibles of a Nile crocodile, a green sea turtle, and a monitor lizard) and performed FEA under identical loading parameters. It was found that once surface scanned models are solidified, they output stress and strain distributions and model deformations comparable to their CT scanned counterparts, though differing by notable stress and strain magnitudes in some cases, depending on morphology of the specimen and the degree of reconstruction applied. Despite similarities in overall mechanical behaviour, surface scanned models can differ in exterior shape compared to CT scanned models due to inaccuracies that can occur during scanning and reconstruction, resulting in local differences in stress distribution. Solid-fill surface scanned models generally output lower stresses compared to CT scanned models due to their compact interiors, which must be accounted for in studies that use both types of scans.

Functional and ecomorphological evolution of orbit shape in mesozoic archosaurs is driven by body size and diet

The orbit is one of several skull openings in the archosauromorph skull. Intuitively, it could be assumed that orbit shape would closely approximate the shape and size of the eyeball resulting in a predominantly circular morphology. However, a quantification of orbit shape across Archosauromorpha using a geometric morphometric approach demonstrates a large morphological diversity despite the fact that the majority of species retained a circular orbit. This morphological diversity is nearly exclusively driven by large (skull length > 1000 mm) and carnivorous species in all studied archosauromorph groups, but particularly prominently in theropod dinosaurs. While circular orbit shapes are retained in most herbivores and smaller species, as well as in juveniles and early ontogenetic stages, large carnivores adopted elliptical and keyhole-shaped orbits. Biomechanical modelling using finite element analysis reveals that these morphologies are beneficial in mitigating and dissipating feeding-induced stresses without additional reinforcement of the bony structure of the skull.

From fibre to function: are we accurately representing muscle architecture and performance?

The size and arrangement of fibres play a determinate role in the kinetic and energetic performance of muscles. Extrapolations between fibre architecture and performance underpin our understanding of how muscles function and how they are adapted to power specific motions within and across species. Here we provide a synopsis of how this 'fibre to function' paradigm has been applied to understand muscle design, performance and adaptation in animals. Our review highlights the widespread application of the fibre to function paradigm across a diverse breadth of biological disciplines but also reveals a potential and highly prevalent limitation running through past studies. Specifically, we find that quantification of muscle architectural properties is almost universally based on an extremely small number of fibre measurements. Despite the volume of research into muscle properties, across a diverse breadth of research disciplines, the fundamental assumption that a small proportion of fibre measurements can accurately represent the architectural properties of a muscle has never been quantitatively tested. Subsequently, we use a combination of medical imaging, statistical analysis, and physics-based computer simulation to address this issue for the first time. By combining diffusion tensor imaging (DTI) and deterministic fibre tractography we generated a large number of fibre measurements (>3000) rapidly for individual human lower limb muscles. Through statistical subsampling simulations of these measurements, we demonstrate that analysing a small number of fibres (n < 25) typically used in previous studies may lead to extremely large errors in the characterisation of overall muscle architectural properties such as mean fibre length and physiological cross-sectional area. Through dynamic musculoskeletal simulations of human walking and jumping, we demonstrate that recovered errors in fibre architecture characterisation have significant implications for quantitative predictions of in-vivo dynamics and muscle fibre function within a species. Furthermore, by applying data-subsampling simulations to comparisons of muscle function in humans and chimpanzees, we demonstrate that error magnitudes significantly impact both qualitative and quantitative assessment of muscle specialisation, potentially generating highly erroneous conclusions about the absolute and relative adaption of muscles across species and evolutionary transitions. Our findings have profound implications for how a broad diversity of research fields quantify muscle architecture and interpret muscle function.

A toolbox for the retrodeformation and muscle reconstruction of fossil specimens in Blender

Accurate muscle reconstructions can offer new information on the anatomy of fossil organisms and are also important for biomechanical analysis (multibody dynamics and finite-element analysis (FEA)). For the sake of simplicity, muscles are often modelled as point-to-point strands or frustra (cut-off cones) in biomechanical models. However, there are cases in which it is useful to model the muscle morphology in three dimensions, to better examine the effects of muscle shape and size. This is especially important for fossil analyses, where muscle force is estimated from the reconstructed muscle morphology (rather than based on data collected in vivo). The two main aims of this paper are as follows. First, we created a new interactive tool in the free open access software Blender to enable interactive three-dimensional modelling of muscles. This approach can be applied to both palaeontological and human biomechanics research to generate muscle force magnitudes and lines of action for FEA. Second, we provide a guide on how to use existing Blender tools to reconstruct distorted or incomplete specimens. This guide is aimed at palaeontologists but can also be used by anatomists working with damaged specimens or to test functional implication of hypothetical morphologies.

A toolbox for the retrodeformation and muscle reconstruction of fossil specimens in Blender

Accurate muscle reconstructions can offer new information on the anatomy of fossil organisms and are also important for biomechanical analysis (multibody dynamics and finite-element analysis (FEA)). For the sake of simplicity, muscles are often modelled as point-to-point strands or frustra (cut-off cones) in biomechanical models. However, there are cases in which it is useful to model the muscle morphology in three dimensions, to better examine the effects of muscle shape and size. This is especially important for fossil analyses, where muscle force is estimated from the reconstructed muscle morphology (rather than based on data collected in vivo). The two main aims of this paper are as follows. First, we created a new interactive tool in the free open access software Blender to enable interactive three-dimensional modelling of muscles. This approach can be applied to both palaeontological and human biomechanics research to generate muscle force magnitudes and lines of action for FEA. Second, we provide a guide on how to use existing Blender tools to reconstruct distorted or incomplete specimens. This guide is aimed at palaeontologists but can also be used by anatomists working with damaged specimens or to test functional implication of hypothetical morphologies.

A toolbox for the retrodeformation and muscle reconstruction of fossil specimens in Blender

Accurate muscle reconstructions can offer new information on the anatomy of fossil organisms and are also important for biomechanical analysis (multibody dynamics and finite-element analysis (FEA)). For the sake of simplicity, muscles are often modelled as point-to-point strands or frustra (cut-off cones) in biomechanical models. However, there are cases in which it is useful to model the muscle morphology in three dimensions, to better examine the effects of muscle shape and size. This is especially important for fossil analyses, where muscle force is estimated from the reconstructed muscle morphology (rather than based on data collected in vivo). The two main aims of this paper are as follows. First, we created a new interactive tool in the free open access software Blender to enable interactive three-dimensional modelling of muscles. This approach can be applied to both palaeontological and human biomechanics research to generate muscle force magnitudes and lines of action for FEA. Second, we provide a guide on how to use existing Blender tools to reconstruct distorted or incomplete specimens. This guide is aimed at palaeontologists but can also be used by anatomists working with damaged specimens or to test functional implication of hypothetical morphologies.

New frontiers in imaging, anatomy, and mechanics of crocodylian jaw muscles

New imaging and biomechanical approaches have heralded a renaissance in our understanding of crocodylian anatomy. Here, we review a series of approaches in the preparation, imaging, and functional analysis of the jaw muscles of crocodylians. Iodine-contrast microCT approaches are enabling new insights into the anatomy of muscles, nerves, and other soft tissues of embryonic as well as adult specimens of alligators. These imaging data and other muscle modeling methods offer increased accuracy of muscle sizes and attachments without destructive methods like dissection. 3D modeling approaches and imaging data together now enable us to see and reconstruct 3D muscle architecture which then allows us to estimate 3D muscle resultants, but also measurements of pennation in ways not seen before. These methods have already revealed new information on the ontogeny, diversity, and function of jaw muscles and the heads of alligators and other crocodylians. Such approaches will lead to enhanced and accurate analyses of form, function, and evolution of crocodylians, their fossil ancestors and vertebrates in general.

The effects of skull flattening on suchian jaw muscle evolution

Jaw muscles are key features of the vertebrate feeding apparatus. The jaw musculature is housed in the skull whose morphology reflects a compromise between multiple functions, including feeding, housing sensory structures, and defense, and the skull constrains jaw muscle geometry. Thus, jaw muscle anatomy may be suboptimally oriented for the production of bite force. Crocodylians are a group of vertebrates that generate the highest bite forces ever measured with a flat skull suited to their aquatic ambush predatory style. However, basal members of the crocodylian line (e.g., Prestosuchus) were terrestrial predators with plesiomorphically tall skulls, and thus the origin of modern crocodylians involved a substantial reorganization of the feeding apparatus and its jaw muscles. Here, we reconstruct jaw muscles across a phylogenetic range of crocodylians and fossil suchians to investigate the impact of skull flattening on muscle anatomy. We used imaging data to create 3D models of extant and fossil suchians that demonstrate the evolution of the crocodylian skull, using osteological correlates to reconstruct muscle attachment sites. We found that jaw muscle anatomy in early fossil suchians reflected the ancestral archosaur condition but experienced progressive shifts in the lineage leading to Metasuchia. In early fossil suchians, musculus adductor mandibulae posterior and musculus pterygoideus (mPT) were of comparable size, but by Metasuchia, the jaw musculature is dominated by mPT. As predicted, we found that taxa with flatter skulls have less efficient muscle orientations for the production of high bite force. This study highlights the diversity and evolution of jaw muscles in one of the great transformations in vertebrate evolution.

Cranial muscle reconstructions quantify adaptation for high bite forces in Oviraptorosauria

Oviraptorosaurians are an unusual and probably herbivorous group of theropod dinosaurs that evolved pneumatised crania with robust, toothless jaws, apparently adapted for producing a strong bite. Using 3D retrodeformed skull models of oviraptorid oviraptorosaurians Citipati, Khaan, and Conchoraptor, along with the earliest diverging oviraptorosaurian, Incisivosaurus, we digitally reconstruct jaw adductor musculature and estimate bite force to investigate cranial function in each species. We model muscle length change during jaw opening to constrain optimal and maximum gape angles. Results demonstrate oviraptorids were capable of much stronger bite forces than herbivorous theropods among Ornithomimosauria and Therizinosauria, relative to body mass and absolutely. Increased bite forces in oviraptorid oviraptorosaurians compared to the earliest diverging oviraptorosaurian result from expanded muscular space and different cranial geometry, not changes in muscular arrangement. Estimated optimal and maximum possible gapes are much smaller than published estimates for carnivorous theropods, being more similar to the herbivorous therizinosaurian theropod Erlikosaurus and modern birds. Restrictive gape and high bite force may represent adaptation towards exploiting tough vegetation, suggesting cranial function and dietary habits differed between oviraptorids and other herbivorous theropods. Differences in the relative strength of jaw adductor muscles between co-occurring oviraptorids may be a factor in niche partitioning, alongside body size.

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.

Cranial functional morphology of the pseudosuchian Effigia and implications for its ecological role in the Triassic

Pseudosuchians, archosaurian reptiles more closely related to crocodylians than to birds, exhibited high morphological diversity during the Triassic with numerous examples of morphological convergence described between Triassic pseudosuchians and post-Triassic dinosaurs. One example is the shuvosaurid Effigia okeeffeae which exhibits an "ostrich-like" bauplan comprising a gracile skeleton with edentulous jaws and large orbits, similar to ornithomimid dinosaurs and extant palaeognaths. This bauplan is regarded as an adaptation for herbivory, but this hypothesis assumes morphological convergence confers functional convergence, and has received little explicit testing. Here, we restore the skull morphology of Effigia, perform myological reconstructions, and apply finite element analysis to quantitatively investigate skull function. We also perform finite element analysis on the crania of the ornithomimid dinosaur Ornithomimus edmontonicus, the extant palaeognath Struthio camelus and the extant pseudosuchian Alligator mississippiensis to assess the degree of functional convergence with a taxon that exhibit "ostrich-like" bauplans and its closest extant relatives. We find that Effigia possesses a mosaic of mechanically strong and weak features, including a weak mandible that likely restricted feeding to the anterior portion of the jaws. We find limited functional convergence with Ornithomimus and Struthio and limited evidence of phylogenetic constraints with extant pseudosuchians. We infer that Effigia was a specialist herbivore that likely fed on softer plant material, a niche unique among the study taxa and potentially among contemporaneous Triassic herbivores. This study increases the known functional diversity of pseudosuchians and highlights that superficial morphological similarity between unrelated taxa does not always imply functional and ecological convergence.

Modeling tooth enamel in FEA comparisons of skulls: Comparing common simplifications with biologically realistic models

Palaeontologists often use finite element analyses, in which forces propagate through objects with specific material properties, to investigate feeding biomechanics. Teeth are usually modeled with uniform properties (all bone or all enamel). In reality, most teeth are composed of pulp, dentine, and enamel. We tested how simplified teeth compare to more realistic models using mandible models of three reptiles. For each, we created models representing enamel thicknesses found in extant taxa, as well as simplified models (bone, dentine or enamel). Our results suggest that general comparisons of stress distribution among distantly related taxa do not require representation of dental tissues, as there was no noticeable effect on heatmap representations of stress. However, we find that representation of dental tissues impacts bite force estimates, although magnitude of these effects may differ depending on constraints. Thus, as others have shown, the detail necessary in a biomechanical model relates to the questions being examined.

The diet of early birds based on modern and fossil evidence and a new framework for its reconstruction

Birds are some of the most diverse organisms on Earth, with species inhabiting a wide variety of niches across every major biome. As such, birds are vital to our understanding of modern ecosystems. Unfortunately, our understanding of the evolutionary history of modern ecosystems is hampered by knowledge gaps in the origin of modern bird diversity and ecosystem ecology. A crucial part of addressing these shortcomings is improving our understanding of the earliest birds, the non-avian avialans (i.e. non-crown birds), particularly of their diet. The diet of non-avian avialans has been a matter of debate, in large part because of the ambiguous qualitative approaches that have been used to reconstruct it. Here we review methods for determining diet in modern and fossil avians (i.e. crown birds) as well as non-avian theropods, and comment on their usefulness when applied to non-avian avialans. We use this to propose a set of comparable, quantitative approaches to ascertain fossil bird diet and on this basis provide a consensus of what we currently know about fossil bird diet. While no single approach can precisely predict diet in birds, each can exclude some diets and narrow the dietary possibilities. We recommend combining (i) dental microwear, (ii) landmark-based muscular reconstruction, (iii) stable isotope geochemistry, (iv) body mass estimations, (v) traditional and/or geometric morphometric analysis, (vi) lever modelling, and (vii) finite element analysis to reconstruct fossil bird diet accurately. Our review provides specific methodologies to implement each approach and discusses complications future researchers should keep in mind. We note that current forms of assessment of dental mesowear, skull traditional morphometrics, geometric morphometrics, and certain stable isotope systems have yet to be proven effective at discerning fossil bird diet. On this basis we report the current state of knowledge of non-avian avialan diet which remains very incomplete. The ancestral dietary condition in non-avian avialans remains unclear due to scarce data and contradictory evidence in Archaeopteryx. Among early non-avian pygostylians, Confuciusornis has finite element analysis and mechanical advantage evidence pointing to herbivory, whilst Sapeornis only has mechanical advantage evidence indicating granivory, agreeing with fossilised ingested material known for this taxon. The enantiornithine ornithothoracine Shenqiornis has mechanical advantage and pedal morphometric evidence pointing to carnivory. In the hongshanornithid ornithuromorph Hongshanornis only mechanical advantage evidence indicates granivory, but this agrees with evidence of gastrolith ingestion in this taxon. Mechanical advantage and ingested fish support carnivory in the songlingornithid ornithuromorph Yanornis. Due to the sparsity of robust dietary assignments, no clear trends in non-avian avialan dietary evolution have yet emerged. Dietary diversity seems to increase through time, but this is a preservational bias associated with a predominance of data from the Early Cretaceous Jehol Lagerstätte. With this new framework and our synthesis of the current knowledge of non-avian avialan diet, we expect dietary knowledge and evolutionary trends to become much clearer in the coming years, especially as fossils from other locations and climates are found. This will allow for a deeper and more robust understanding of the role birds played in Mesozoic ecosystems and how this developed into their pivotal role in modern ecosystems.

Unexpected bite-force conservatism as a stable performance foundation across mesoeucrocodylian historical diversity

Effective interpretation of historical selective regimes requires comprehensive in vivo performance evaluations and well-constrained ecomorphological proxies. The feeding apparatus is a frequent target of such evolutionary studies due to a direct relationship between feeding and survivorship, and the durability of craniodental elements in the fossil record. Among vertebrates, behaviors such as bite force have been central to evaluation of clade dynamics; yet, in the absence of detailed performance studies, such evaluations can misidentify potential selective factors and their roles. Here, we combine the results of a total-clade performance study with fossil-inclusive, phylogenetically informed methods to assess bite-force proxies throughout mesoeucrocodylian evolution. Although bite-force shifts were previously thought to respond to changing rostrodental selective regimes, we find body-size dependent conservation of performance proxies throughout the history of the clade, indicating stabilizing selection for bite-force potential. Such stasis reveals that mesoeucrocodylians with dietary ecologies as disparate as herbivory and hypercarnivory maintain similar bite-force-to-body-size relationships, a pattern which contrasts the precept that vertebrate bite forces should vary most strongly by diet. Furthermore, it may signal that bite-force conservation supported mesoeucrocodylian craniodental disparity by providing a stable performance foundation for the exploration of novel ecomorphospace.

Biomechanical performance of the cranio-mandibular complex of the small notosuchian Araripesuchus gomesii (Notosuchia, Uruguaysuchidae)

Notosuchia is a clade of crocodyliforms that was highly successful and diverse in the Cretaceous of Gondwana. Araripesuchus gomesii is a small notosuchian from the Early Cretaceous of Brazil that belongs to Uruguaysuchidae, one of the subgroups of notosuchians that first radiated, during the Aptian-Albian. Here we present a finite element analysis of A. gomesii based on a model reconstructed from CT scans and performed using published bone properties for crocodiles. The adductor musculature and their respective attachment areas were reconstructed based on Extant Phylogenetic Bracket. Different functional scenarios were tested applying an estimated 158 N bite force: unilateral bite, bilateral bite, pullback, head-shake, and head-twist. The results obtained were compared with those of Alligator mississippiensis, one of its closest living relatives. In the different simulations, the skull and lower jaws of Araripesuchus suffers more stress in the head-shake movement, followed by the unilateral and pullback bites with stress focalized in the premaxillary region. In contrast, the head-twist is the one with smaller stress values. Araripesuchus possess an oreinirostral skull that may provide greater overall resistance in the different scenarios on average, unlike Alligator that has a platyrostral skull with less resistance to dorsoventral mechanical loads. Previous hypotheses that considered A. gomesii as omnivorous coupled with our results, its small size, and likely limited bite force, suggest this taxon probably fed on small prey and other trophic items that could catch and handle entirely with its mouth, such as insects and small vertebrates.

Biomechanics of juvenile tyrannosaurid mandibles and their implications for bite force: Evolutionary biology

The tyrannosaurids are among the most well-studied dinosaurs described by science, and analysis of their feeding biomechanics allows for comparison between established tyrannosaurid genera and across ontogeny. 3D finite element analysis (FEA) was used to model and quantify the mechanical properties of the mandibles (lower jaws) of three tyrannosaurine tyrannosaurids of different sizes. To increase evolutionary scope and context for 3D tyrannosaurine results, a broader sample of validated 2D mandible FEA enabled comparisons between ontogenetic stages of Tyrannosaurus rex and other large theropods. It was found that mandibles of small juvenile and large subadult tyrannosaurs experienced lower stress overall because muscle forces were relatively lower, but experienced greater simulated stresses at decreasing sizes when specimen muscle force is normalized. The strain on post-dentary ligaments decreases stress and strain in the posterior region of the dentary and where teeth impacted food. Tension from the lateral insertion of the looping m. ventral pterygoid muscle increases compressive stress on the angular but may decrease anterior bending stress on the mandible. Low mid-mandible bending stresses are congruent with ultra-robust teeth and high anterior bite force in adult T. rex. Mandible strength increases with size through ontogeny in T. rex and phylogenetically among other tyrannosaurids, in addition to that tyrannosaurid mandibles exceed the mandible strength of other theropods at equivalent ramus length. These results may indicate separate predatory strategies used by juvenile and mature tyrannosaurids; juvenile tyrannosaurids lacked the bone-crunching bite of adult specimens and hunted smaller prey, while adult tyrannosaurids fed on larger prey.

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.

Quantitatively assessing mekosuchine crocodile locomotion by geometric morphometric and finite element analysis of the forelimb

Morphological shifts observed in the fossil record of a lineage potentially indicate concomitant shifts in ecology of that lineage. Mekosuchine crocodiles of Cenozoic Australia display departures from the typical eusuchian body-plan both in the cranium and postcranium. Previous qualitative studies have suggested that these crocodiles had a more terrestrial habitus than extant crocodylians, yet the capacity of mekosuchine locomotion remains to be tested. Limb bone shape, such as diaphyseal cross-section and curvature, has been related to habitual use and locomotory function across a wide variety of taxa. Available specimens of mekosuchine limbs, primarily humeri, are distinctly columnar compared with those of extant crocodylians. Here we apply a quantitative approach to biomechanics in mekosuchine taxa using both geomorphic morphometric and finite element methods to measure bone shape and estimate locomotory stresses in a comparative context. Our results show mekosuchines appear to diverge from extant semi-aquatic saltwater and freshwater crocodiles in cross-sectional geometry of the diaphysis and generate different structural stresses between models that simulate sprawling and high-walk gaits. The extant crocodylians display generally rounded cross-sectional diaphyseal outlines, which may provide preliminary indication of resistance to torsional loads that predominate during sprawling gait, whereas mekosuchine humeri appear to vary between a series of elliptical outlines. Mekosuchine structural stresses are comparatively lower than those of the extant crocodylians and reduce under high-walk gait in some instances. This appears to be a function of bending moments induced by differing configurations of diaphyseal curvature. Additionally, the neutral axis of structural stresses is differently oriented in mekosuchines. This suggests a shift in the focus of biomechanical optimisation, from torsional to axial loadings. Our results lend quantitative support to the terrestrial habitus hypothesis in so far as they suggest that mekosuchine humeri occupied a different morphospace than that associated with the semi-aquatic habit. The exact adaptational trajectory of mekosuchines, however, remains to be fully quantified. Novel forms appear to emerge among mekosuchines during the late Cenozoic. Their adaptational function is considered here; possible applications include navigation of uneven terrain and burrowing.

A unique predator in a unique ecosystem: modelling the apex predator within a Late Cretaceous crocodyliform-dominated fauna from Brazil

Theropod dinosaurs were relatively scarce in the Late Cretaceous ecosystems of southeast Brazil. Instead, hypercarnivorous crocodyliforms known as baurusuchids were abundant and probably occupied the ecological role of apex predators. Baurusuchids exhibited a series of morphological adaptations hypothesized to be associated with this ecological role, but quantitative biomechanical analyses of their morphology have so far been lacking. Here, we employ a biomechanical modelling approach, applying finite element analysis (FEA) to models of the skull and mandibles of a baurusuchid specimen. This allows us to characterize the craniomandibular apparatus of baurusuchids, as well as to compare the functional morphology of the group with that of other archosaurian carnivores, such as theropods and crocodylians. Our results support the ecological role of baurusuchids as specialized apex predators in the continental Late Cretaceous ecosystems of South America. With a relatively weak bite force (~600 N), the predation strategies of baurusuchids likely relied on other morphological specializations, such as ziphodont dentition and strong cervical musculature. Comparative assessments of the stress distribution and magnitude of scaled models of other predators (the theropod Allosaurus fragilis and the living crocodylian Alligator mississippiensis) consistently show different responses to loadings under the same functional scenarios, suggesting distinct predatory behaviors for these animals. The unique selective pressures in the arid to semi-arid Late Cretaceous ecosystems of southeast Brazil, which were dominated by crocodyliforms, possibly drove the emergence and evolution of the biomechanical features seen in baurusuchids, which are distinct from those previously reported for other predatory taxa.

Palatal Biomechanics and Its Significance for Cranial Kinesis in Tyrannosaurus rex

The extinct nonavian dinosaur Tyrannosaurus rex, considered one of the hardest biting animals ever, is often hypothesized to have exhibited cranial kinesis, or, mobility of cranial joints relative to the braincase. Cranial kinesis in T. rex is a biomechanical paradox in that forcefully biting tetrapods usually possess rigid skulls instead of skulls with movable joints. We tested the biomechanical performance of a tyrannosaur skull using a series of static positions mimicking possible excursions of the palate to evaluate Postural Kinetic Competency in Tyrannosaurus. A functional extant phylogenetic bracket was employed using taxa, which exhibit measurable palatal excursions: Psittacus erithacus (fore-aft movement) and Gekko gecko (mediolateral movement). Static finite element models of Psittacus, Gekko, and Tyrannosaurus were constructed and tested with different palatal postures using anatomically informed material properties, loaded with muscle forces derived from dissection, phylogenetic bracketing, and a sensitivity analysis of muscle architecture and tested in orthal biting simulations using element strain as a proxy for model performance. Extant species models showed lower strains in naturally occurring postures compared to alternatives. We found that fore-aft and neutral models of Tyrannosaurus experienced lower overall strains than mediolaterally shifted models. Protractor muscles dampened palatal strains, while occipital constraints increased strains about palatocranial joints compared to jaw joint constraints. These loading behaviors suggest that even small excursions can strain elements beyond structural failure. Thus, these postural tests of kinesis, along with the robusticity of other cranial features, suggest that the skull of Tyrannosaurus was functionally akinetic. Anat Rec, 303:999-1017, 2020. © 2019 Wiley Periodicals, Inc.

Digital dissection of the head of the rock dove (Columba livia) using contrast-enhanced computed tomography

The rock dove (or common pigeon), Columba livia, is an important model organism in biological studies, including research focusing on head muscle anatomy, feeding kinematics, and cranial kinesis. However, no integrated computer-based biomechanical model of the pigeon head has yet been attempted. As an initial step towards achieving this goal, we present the first three-dimensional digital dissection of the pigeon head based on a contrast-enhanced computed tomographic dataset achieved using iodine potassium iodide as a staining agent. Our datasets enable us to visualize the skeletal and muscular anatomy, brain and cranial nerves, and major sense organs of the pigeon, including very small and fragile features, as well as maintaining the three-dimensional topology of anatomical structures. This work updates and supplements earlier anatomical work on this widely used laboratory organism. We resolve several key points of disagreement arising from previous descriptions of pigeon anatomy, including the precise arrangement of the external adductor muscles and their relationship to the posterior adductor. Examination of the eye muscles highlights differences between avian taxa and shows that pigeon eye muscles are more similar to those of a tinamou than they are to those of a house sparrow. Furthermore, we present our three-dimensional data as publicly accessible files for further research and education purposes. Digital dissection permits exceptional visualisation and will be a valuable resource for further investigations into the head anatomy of other bird species, as well as efforts to reconstruct soft tissues in fossil archosaurs.

Anatomy and Ontogeny of the Mandibular Symphysis in Alligator mississippiensis

Crocodylians evolved some of the most characteristic skulls of the animal kingdom with specializations for semiaquatic and ambush lifestyles, resulting in a feeding apparatus capable of tolerating high biomechanical loads and bite forces and a head with a derived sense of trigeminal-nerve-mediated touch. The mandibular symphysis accommodates these specializations being both at the end of a biomechanical lever and an antenna for sensation. Little is known about the anatomy of the crocodylian mandibular symphysis, hampering our understanding of form, function, and evolution of the joint in extant and extinct lineages. We explore mandibular symphysis anatomy of an ontogenetic series of Alligator mississippiensis using imaging, histology, and whole mount methods. Complex sutural ligaments emanating about a midline-fused Meckel's cartilage bridge the symphysis. These tissues organize during days 37-42 of in ovo development. However, interdigitations do not manifest until after hatching. These soft tissues leave a hub and spoke-like bony morphology of the symphyseal plate, which never fuses. Interdigitation morphology varies within the symphysis suggesting differential loading about the joint. Neurovascular canals extend throughout the mandibles to alveoli, integument, and bone adjacent to the symphysis. These features suggest the Alligator mandibular symphysis offers compliance in an otherwise rigid skull. We hypothesize a fused Meckel's cartilage offers stiffness in hatchling mandibles prior to the development of organized sutural ligaments and mineralized bone while offering a scaffold for somatic growth. The porosity of the dentaries due to neurovascular tissues likely allows transmission of sensory and proprioceptive information from the surroundings and the loaded symphysis. Anat Rec, 302:1696-1708, 2019. © 2019 American Association for Anatomy.

The importance of muscle architecture in biomechanical reconstructions of extinct animals: a case study using Tyrannosaurus rex

Functional reconstructions of extinct animals represent a crucial step towards understanding palaeocological interactions, selective pressures and macroevolutionary patterns in the fossil record. In recent years, computational approaches have revolutionised the field of 'evolutionary biomechanics' and have, in general, resulted in convergence of quantitative estimates of performance on increasingly narrow ranges for well studied taxa. Studies of body mass and locomotor performance of Tyrannosaurus rex - arguably the most intensively studied extinct animal - typify this pattern, with numerous independent studies predicting similar body masses and maximum locomotor speeds for this animal. In stark contrast to this trend, recent estimates of maximum bite force in T. rex vary considerably (> 50%) despite use of similar quantitative methodologies. Herein we demonstrate that the mechanistic causes of these disparate predictions are indicative of important and underappreciated limiting factors in biomechanical reconstructions of extinct organisms. Detailed comparison of previous models of T. rex bite force reveals that estimations of muscle fibre lengths and architecture are the principal source of disagreement between studies, and therefore that these parameters represents the greatest source of uncertainty in these reconstructions, and potentially therefore extinct animals generally. To address the issue of fibre length and architecture estimation in extinct animals we present data tabulated from the literature of muscle architecture from over 1100 muscles measured in extant terrestrial animals. Application of this dataset in a reanalysis of T. rex bite force emphasises the need for more data on jaw musculature from living carnivorous animals, alongside increased sophistication of modelling approaches. In the latter respect we predict that implementing limits on skeletal loading into musculoskeletal models will narrow predictions for T. rex bite force by excluding higher-end estimates.

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.

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.

Ontogeny of bite force in a validated biomechanical model of the American alligator

Three-dimensional computational modeling offers tools with which to investigate forces experienced by the skull encountered during feeding and other behaviors. American alligators (Alligator mississippiensis) generate some of the highest measured bite forces among extant tetrapods. A concomitant increase in bite force accompanies ontogenetic increases in body mass, which has been linked with dietary changes as animals increase in size. Because the flattened skull of crocodylians has substantial mediolaterally-oriented muscles, they are an excellent model taxon in which to explore the role of mediolateral force components experienced by the feeding apparatus. Many previous modeling studies of archosaur cranial function focused on planar analysis, ignoring the mediolateral aspects of cranial forces. Here we use three-dimensionally accurate anatomical data to resolve 3D muscle forces. Using dissection, imaging, and computational techniques, we developed lever and finite element models of an ontogenetic series of alligators to test the effects of size and shape on cranial loading and compared estimated bite forces to those previously measured in vivo in Alligator mississippiensis. We found that modeled forces matched in vivo data well for intermediately sized individuals, and somewhat overestimated force in smaller specimens and underestimated force in larger specimens, suggesting that ontogenetically static muscular parameters and bony attachment sites alone cannot account for all the variation in bite force. Adding aponeurotic muscle attachments would likely improve force predictions, but such data are challenging to model and integrate into analyses of extant taxa and are generally unpreserved in fossils. We conclude that anatomically accurate modeling of muscles can be coupled with finite element and lever analyses to produce reliable, reasonably accurate estimate bite forces and thus both skeletal and joint loading, with known sources of error, which can be applied to extinct taxa.

The Predentary Bone and Its Significance in the Evolution of Feeding Mechanisms in Ornithischian Dinosaurs

The characteristic predentary bone in ornithischian dinosaurs is a unique, unpaired element located at the midline of the mandibular symphysis. Although traditionally thought to only be a plant "nipping" bone, the true functional significance of this bone among feeding mechanisms of ornithischian dinosaurs is poorly known. Recent studies of a select few ornithischian genera have suggested rotation of the mandibular corpora around their long axes relative to their midline joint articulation with the predentary bone. This study aims to re-evaluate these hypotheses as well as provide in-depth qualitative comparative descriptions of predentary bone morphology in ornithischian genera throughout all subclades, including heterodontosaurids, thyreophorans, ornithopods, and marginocephalians. Descriptions evaluate overall shape of the predentary, its articular surfaces contacting the rostral ends of the dentaries, and the morphology of the rostral extent of the dentaries and their midline symphysis. Functionally relevant morphologies in each predentary morphotype are accentuated for further speculation of feeding mechanisms. Three predentary morphotypes are described throughout ornithischian subclades and each plays a unique role in feeding adaptations. Most notably, the predentary likely evolved as a midline axial point of the mandibular symphysis for simultaneous variable movement or rotation of the mandibular corpora in many, but not all, taxa. This simultaneous movement of the hemimandibles would have aided in feeding on both sides of the jaw at once. The function of the predentary as well as other jaw adaptations is discussed for genera throughout all subclades, focusing on both general shape and joint morphology. Anat Rec, 299:1358-1388, 2016. © 2016 Wiley Periodicals, Inc.

Comparative 3D analyses and palaeoecology of giant early amphibians (Temnospondyli: Stereospondyli)

Macroevolutionary, palaeoecological and biomechanical analyses in deep time offer the possibility to decipher the structural constraints, ecomorphological patterns and evolutionary history of extinct groups. Here, 3D comparative biomechanical analyses of the extinct giant early amphibian group of stereospondyls together with living lissamphibians and crocodiles, shows that: i) stereospondyls had peculiar palaeoecological niches with proper bites and stress patterns very different than those of giant salamanders and crocodiles; ii) their extinction may be correlated with the appearance of neosuchians, which display morphofunctional innovations. Stereospondyls weathered the end-Permian mass extinction, re-radiated, acquired gigantic sizes and dominated (semi) aquatic ecosystems during the Triassic. Because these ecosystems are today occupied by crocodilians, and stereospondyls are extinct amphibians, their palaeobiology is a matter of an intensive debate: stereospondyls were a priori compared with putative living analogous such as giant salamanders and/or crocodilians and our new results try to close this debate.

Estimating cranial musculoskeletal constraints in theropod dinosaurs

Many inferences on the biology, behaviour and ecology of extinct vertebrates are based on the reconstruction of the musculature and rely considerably on its accuracy. Although the advent of digital reconstruction techniques has facilitated the creation and testing of musculoskeletal hypotheses in recent years, muscle strain capabilities have rarely been considered. Here, a digital modelling approach using the freely available visualization and animation software Blender is applied to estimate cranial muscle length changes and optimal and maximal possible gape in different theropod dinosaurs. Models of living archosaur taxa (Alligator mississippiensis, Buteo buteo) were used in an extant phylogenetically bracketed framework to validate the method. Results of this study demonstrate that Tyrannosaurus rex, Allosaurus fragilis and Erlikosaurus andrewsi show distinct differences in the recruitment of the jaw adductor musculature and resulting gape, confirming previous dietary and ecological assumptions. While the carnivorous taxa T. rex and Allo. fragilis were capable of a wide gape and sustained muscle force, the herbivorous therizinosaurian E. andrewsi was constrained to small gape angles.

The sensitivity of biological finite element models to the resolution of surface geometry: a case study of crocodilian crania

The reliability of finite element analysis (FEA) in biomechanical investigations depends upon understanding the influence of model assumptions. In producing finite element models, surface mesh resolution is influenced by the resolution of input geometry, and influences the resolution of the ensuing solid mesh used for numerical analysis. Despite a large number of studies incorporating sensitivity studies of the effects of solid mesh resolution there has not yet been any investigation into the effect of surface mesh resolution upon results in a comparative context. Here we use a dataset of crocodile crania to examine the effects of surface resolution on FEA results in a comparative context. Seven high-resolution surface meshes were each down-sampled to varying degrees while keeping the resulting number of solid elements constant. These models were then subjected to bite and shake load cases using finite element analysis. The results show that incremental decreases in surface resolution can result in fluctuations in strain magnitudes, but that it is possible to obtain stable results using lower resolution surface in a comparative FEA study. As surface mesh resolution links input geometry with the resulting solid mesh, the implication of these results is that low resolution input geometry and solid meshes may provide valid results in a comparative context.

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.

Bone-breaking bite force of Basilosaurus isis (Mammalia, Cetacea) from the late Eocene of Egypt estimated by finite element analysis

Bite marks suggest that the late Eocence archaeocete whale Basilosaurus isis (Birket Qarun Formation, Egypt) fed upon juveniles of the contemporary basilosaurid Dorudon atrox. Finite element analysis (FEA) of a nearly complete adult cranium of B. isis enables estimates of its bite force and tests the animal's capabilities for crushing bone. Two loadcases reflect different biting scenarios: 1) an intitial closing phase, with all adductors active and a full condylar reaction force; and 2) a shearing phase, with the posterior temporalis active and minimized condylar force. The latter is considered probable when the jaws were nearly closed because the preserved jaws do not articulate as the molariform teeth come into occulusion. Reaction forces with all muscles active indicate that B. isis maintained relatively greater bite force anteriorly than seen in large crocodilians, and exerted a maximum bite force of at least 16,400 N at its upper P3. Under the shearing scenario with minimized condylar forces, tooth reaction forces could exceed 20,000 N despite lower magnitudes of muscle force. These bite forces at the teeth are consistent with bone indentations on Dorudon crania, reatract-and-shear hypotheses of Basilosaurus bite function, and seizure of prey by anterior teeth as proposed for other archaeocetes. The whale's bite forces match those estimated for pliosaurus when skull lengths are equalized, suggesting similar tradeoffs of bite function and hydrodynamics. Reaction forces in B. isis were lower than maxima estimated for large crocodylians and carnivorous dinosaurs. However, comparison of force estimates from FEA and regression data indicate that B. isis exerted the largest bite forces yet estimated for any mammal, and greater force than expected from its skull width. Cephalic feeding biomechanics of Basilosaurus isis are thus consistent with habitual predation.