Palatal Biomechanics and Its Significance for Cranial Kinesis in Tyrannosaurus rex

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.

Cranial functional specialisation for strength precedes morphological evolution in Oviraptorosauria

Oviraptorosaurians were a theropod dinosaur group that reached high diversity in the Late Cretaceous. Within oviraptorosaurians, the later diverging oviraptorids evolved distinctive crania which were extensively pneumatised, short and tall, and had a robust toothless beak, interpreted as providing a powerful bite for their herbivorous to omnivorous diet. The present study explores the ability of oviraptorid crania to resist large mechanical stresses compared with other theropods and where this adaptation originated within oviraptorosaurians. Digital 3D cranial models were constructed for the earliest diverging oviraptorosaurian, Incisivosaurus gauthieri, and three oviraptorids, Citipati osmolskae, Conchoraptor gracilis, and Khaan mckennai. Finite element analyses indicate oviraptorosaurian crania were stronger than those of other herbivorous theropods (Erlikosaurus and Ornithomimus) and were more comparable to the large, carnivorous Allosaurus. The cranial biomechanics of Incisivosaurus align with oviraptorids, indicating an early establishment of distinctive strengthened cranial biomechanics in Oviraptorosauria, even before the highly modified oviraptorid cranial morphology. Bite modelling, using estimated muscle forces, suggests oviraptorid crania may have functioned closer to structural safety limits. Low mechanical stresses around the beaks of oviraptorids suggest a convergently evolved, functionally distinct rhamphotheca, serving as a cropping/feeding tool rather than for stress reduction, when compared with other herbivorous theropods.

Functional morphology of the Triassic apex predator Saurosuchus galilei (Pseudosuchia: Loricata) and convergence with a post-Triassic theropod dinosaur

Pseudosuchian archosaurs, reptiles more closely related to crocodylians than to birds, exhibited high morphological diversity during the Triassic and are thus associated with hypotheses of high ecological diversity during this time. One example involves basal loricatans which are non-crocodylomorph pseudosuchians traditionally known as "rauisuchians." Their large size (5-8+ m long) and morphological similarities to post-Triassic theropod dinosaurs, including dorsoventrally deep skulls and serrated dentitions, suggest basal loricatans were apex predators. However, this hypothesis does not consider functional behaviors that can influence more refined roles of predators in their environment, for example, degree of carcass utilization. Here, we apply finite element analysis to a juvenile but three-dimensionally well-preserved cranium of the basal loricatan Saurosuchus galilei to investigate its functional morphology and to compare with stress distributions from the theropod Allosaurus fragilis to assess degrees of functional convergence between Triassic and post-Triassic carnivores. We find similar stress distributions and magnitudes between the two study taxa under the same functional simulations, indicating that Saurosuchus had a somewhat strong skull and thus exhibited some degree of functional convergence with theropods. However, Saurosuchus also had a weak bite for an animal of its size (1015-1885 N) that is broadly equivalent to the bite force of modern gharials (Gavialis gangeticus). We infer that Saurosuchus potentially avoided tooth-bone interactions and consumed the softer parts of carcasses, unlike theropods and other basal loricatans. This deduced feeding mode for Saurosuchus increases the known functional diversity of basal loricatans and highlights functional differences between Triassic and post-Triassic apex predators.

The dual function of prokinesis in the feeding and locomotor systems of parrots

Prokinesis, a mode of avian cranial kinesis involving motion between the neurocranium and upper beak, has long been investigated in biomechanical analyses of avian feeding and drinking. However, the modern avian beak is also used for non-feeding functions. Here, we investigate the dual function of prokinesis in the feeding and locomotor systems of the rosy-faced lovebird (Agapornis roseicollis). Lovebirds and other parrots utilize their beak both during feeding and as a third limb during vertical climbing. Thus, we experimentally measured both force-generating potential and movement of the rosy-faced lovebird mandible and maxilla (via prokinetic flexion of the craniofacial hinge) during tripedal climbing and mandibular/maxillary adduction. We found that whereas the maxilla is primarily responsible for generating force during locomotion, the mandible is primarily responsible for generating force during forceful jaw adduction, hinting at a remarkable capacity to alter prokinetic function with differing neuromuscular control. The ability of the prokinetic apparatus to perform functions with competing optimality criteria via modulation of motor control illustrates the functional plasticity of the avian cranial kinesis and sheds new light on the adaptive significance of cranial mobility.

Orbital soft tissues, bones, and allometry: Implications for the size and position of crocodylian eyes

Although the visual system of crocodylians has attracted interest regarding optical parameters and retinal anatomy, fundamental questions remain about the allometry of the eyeball and whether such scaling is the same across all crown groups of crocodylians. In addition, anatomy and identities of adnexal soft tissues that interact with the visual system are not well understood in many cases. We used contrast-enhancing iodine stain and high-resolution micro-computed tomography to assess the anatomy of orbital soft tissues, including extraocular muscles and glands, in crocodylians. We also used regression analysis to estimate the allometric relationship between the bony orbit and eyeball across Alligator mississippiensis and Crocodylus niloticus for the first time. Results revealed tight, negatively allometric relationships between the bony orbit and eyeball. Notably, the eyes of C. niloticus were larger for a given orbit size than the eyes of A. mississippiensis, although the slope of the relationship was no different between these two crown crocodylian groups. Among the findings from our anatomical study, new details were uncovered about the homologies of muscles of the abducens complex. In particular, M. rectus lateralis of crocodylians is revealed to have a more complex form than previously appreciated, being adhered to the tendon of the nictitating membrane, which may be apomorphic for Crocodylia. Our calculation of the orbit-eyeball allometric relationship and study of the adnexal soft tissues of the crocodylian visual system, in combination with previous work by other teams in other crown saurian clades, is a critical, formerly missing, piece in the Extant Phylogenetic Bracket for restoring the visual apparatus of extinct crocodyliforms and other archosauriform groups.

Turtle Shell Kinesis Underscores Constraints and Opportunities in the Evolution of the Vertebrate Musculoskeletal System

Species groups that feature traits with a low number of potentially variable (evolvable) character states are more likely to repeatedly evolve similar phenotypes, that is, convergence. To evaluate this phenomenon, this present paper addresses anatomical alterations in turtles that convergently evolved shell kinesis, for example, the movement of shell bones to better shield the head and extremities. Kinesis constitutes a major departure from the evolutionarily conserved shell of modern turtles, yet it has arisen independently at least 8 times. The hallmark signature of kinesis is the presence of shell bone articulations or "hinges," which arise via similar skeletal remodeling processes in species that do not share a recent common ancestor. Still, the internal biomechanical components that power kinesis may differ in such distantly related species. Complex diarthrodial joints and modified muscle connections expand the functional boundaries of the limb girdles and neck in a lineage-specific manner. Some lineages even exhibit mobility of thoracic and sacral vertebrae to facilitate shell closure. Depending on historical contingency and structural correlation, a myriad of anatomical alterations has yielded similar functional outcomes, that is, many-to-one mapping, during the convergent evolution of shell kinesis. The various iterations of this intricate phenotype illustrate the potential for the vertebrate musculoskeletal system to undergo evolutionary change, even when constraints are imposed by the development and structural complexity of a shelled body plan. Based on observations in turtles and comparisons to other vertebrates, a hypothetical framework that implicates functional interactions in the origination of novel musculoskeletal traits is presented.

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

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.

One step further in biomechanical models in palaeontology: a nonlinear finite element analysis review

Finite element analysis (FEA) is no longer a new technique in the fields of palaeontology, anthropology, and evolutionary biology. It is nowadays a well-established technique within the virtual functional-morphology toolkit. However, almost all the works published in these fields have only applied the most basic FEA tools i.e., linear materials in static structural problems. Linear and static approximations are commonly used because they are computationally less expensive, and the error associated with these assumptions can be accepted. Nonetheless, nonlinearities are natural to be used in biomechanical models especially when modelling soft tissues, establish contacts between separated bones or the inclusion of buckling results. The aim of this review is to, firstly, highlight the usefulness of non-linearities and secondly, showcase these FEA tool to researchers that work in functional morphology and biomechanics, as non-linearities can improve their FEA models by widening the possible applications and topics that currently are not used in palaeontology and anthropology.

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.

Estimating bite force in extinct dinosaurs using phylogenetically predicted physiological cross-sectional areas of jaw adductor muscles

I present a Bayesian phylogenetic predictive modelling (PPM) framework that allows the prediction of muscle parameters (physiological cross-sectional area, A _Phys) in extinct archosaurs from skull width (W _Sk) and phylogeny. This approach is robust to phylogenetic uncertainty and highly versatile given its ability to base predictions on simple, readily available predictor variables. The PPM presented here has high prediction accuracy (up to 95%), with downstream biomechanical modelling yielding bite force estimates that are in line with previous estimates based on muscle parameters from reconstructed muscles. This approach does not replace muscle reconstructions but one that provides a powerful means to predict A _Phys from skull geometry and phylogeny to the same level of accuracy as that measured from reconstructed muscles in species for which soft tissue data are unavailable or difficult to obtain.

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.

Skull Sutures and Cranial Mechanics in the Permian Reptile Captorhinus aguti and the Evolution of the Temporal Region in Early Amniotes

While most early limbed vertebrates possessed a fully-roofed dermatocranium in their temporal skull region, temporal fenestrae and excavations evolved independently at least twice in the earliest amniotes, with several different variations in shape and position of the openings. Yet, the specific drivers behind this evolution have been only barely understood. It has been mostly explained by adaptations of the feeding apparatus as a response to new functional demands in the terrestrial realm, including a rearrangement of the jaw musculature as well as changes in strain distribution. Temporal fenestrae have been retained in most extant amniotes but have also been lost again, notably in turtles. However, even turtles do not represent an optimal analog for the condition in the ancestral amniote, highlighting the necessity to examine Paleozoic fossil material. Here, we describe in detail the sutures in the dermatocranium of the Permian reptile Captorhinus aguti (Amniota, Captorhinidae) to illustrate bone integrity in an early non-fenestrated amniote skull. We reconstruct the jaw adductor musculature and discuss its relation to intracranial articulations and bone flexibility within the temporal region. Lastly, we examine whether the reconstructed cranial mechanics in C. aguti could be treated as a model for the ancestor of fenestrated amniotes. We show that C. aguti likely exhibited a reduced loading in the areas at the intersection of jugal, squamosal, and postorbital, as well as at the contact between parietal and postorbital. We argue that these “weak” areas are prone for the development of temporal openings and may be treated as the possible precursors for infratemporal and supratemporal fenestrae in early amniotes. These findings provide a good basis for future studies on other non-fenestrated taxa close to the amniote base, for example diadectomorphs or other non-diapsid reptiles.

Diet of Mesozoic toothed birds (Longipterygidae) inferred from quantitative analysis of extant avian diet proxies

Background

Birds are key indicator species in extant ecosystems, and thus we would expect extinct birds to provide insights into the nature of ancient ecosystems. However, many aspects of extinct bird ecology, particularly their diet, remain obscure. One group of particular interest is the bizarre toothed and long-snouted longipterygid birds. Longipterygidae is the most well-understood family of enantiornithine birds, the dominant birds of the Cretaceous period. However, as with most Mesozoic birds, their diet remains entirely speculative.

Results

To improve our understanding of longipterygids, we investigated four proxies in extant birds to determine diagnostic traits for birds with a given diet: body mass, claw morphometrics, jaw mechanical advantage, and jaw strength via finite element analysis. Body mass of birds tended to correspond to the size of their main food source, with both carnivores and herbivores splitting into two subsets by mass: invertivores or vertivores for carnivores, and granivores + nectarivores or folivores + frugivores for herbivores. Using claw morphometrics, we successfully distinguished ground birds, non-raptorial perching birds, and raptorial birds from one another. We were unable to replicate past results isolating subtypes of raptorial behaviour. Mechanical advantage was able to distinguish herbivorous diets with particularly high values of functional indices, and so is useful for identifying these specific diets in fossil taxa, but overall did a poor job of reflecting diet. Finite element analysis effectively separated birds with hard and/or tough diets from those eating foods which are neither, though could not distinguish hard and tough diets from one another. We reconstructed each of these proxies in longipterygids as well, and after synthesising the four lines of evidence, we find all members of the family but Shengjingornis (whose diet remains inconclusive) most likely to be invertivores or generalist feeders, with raptorial behaviour likely in Longipteryx and Rapaxavis.

Conclusions

This study provides a 20% increase in quantitatively supported fossil bird diets, triples the number of diets reconstructed in enantiornithine species, and serves as an important first step in quantitatively investigating the origins of the trophic diversity of living birds. These findings are consistent with past hypotheses that Mesozoic birds occupied low trophic levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01294-3.

2D and 3D visualizations of archosaur jaw muscle mechanics, ontogeny and phylogeny using ternary diagrams and 3D modeling

Comparing patterns of performance and kinematics across behavior, development and phylogeny is crucial to understand the evolution of complex musculoskeletal systems such as the feeding apparatus. However, conveying 3D spatial data of muscle orientation throughout a feeding cycle, ontogenetic pathway or phylogenetic lineage is essential to understanding the function and evolution of the skull in vertebrates. Here, we detail the use of ternary plots for displaying and comparing the 3D orientation of muscle data. First, we illustrate changes in 3D jaw muscle resultants during jaw closing taxa the American alligator (Alligator mississippiensis). Second, we show changes in 3D muscle resultants of jaw muscles across an ontogenetic series of alligators. Third, we compare 3D resultants of jaw muscles of avian-line dinosaurs, including extant (Struthio camelus, Gallus gallus, Psittacus erithacus) and extinct (Tyrannosaurus rex) species to outline the reorganization of jaw muscles that occurred along the line to modern birds. Finally, we compare 3D resultants of jaw muscles of the hard-biting species in our sample (A. mississippiensis, T. rex, P. erithacus) to illustrate how disparate jaw muscle resultants are employed in convergent behaviors in archosaurs. Our findings show that these visualizations of 3D components of jaw muscles are immensely helpful towards identifying patterns of cranial performance, growth and diversity. These tools will prove useful for testing other hypotheses in functional morphology, comparative biomechanics, ecomorphology and organismal evolution.

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.

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

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

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

Evolutionary biomechanics: hard tissues and soft evidence?

Biomechanical modelling is a powerful tool for quantifying the evolution of functional performance in extinct animals to understand key anatomical innovations and selective pressures driving major evolutionary radiations. However, the fossil record is composed predominantly of hard parts, forcing palaeontologists to reconstruct soft tissue properties in such models. Rarely are these reconstruction approaches validated on extant animals, despite soft tissue properties being highly determinant of functional performance. The extent to which soft tissue reconstructions and biomechanical models accurately predict quantitative or even qualitative patterns in macroevolutionary studies is therefore unknown. Here, we modelled the masticatory system in extant rodents to objectively test the ability of current muscle reconstruction methods to correctly identify quantitative and qualitative differences between macroevolutionary morphotypes. Baseline models generated using measured soft tissue properties yielded differences in muscle proportions, bite force, and bone stress expected between extant sciuromorph, myomorph, and hystricomorph rodents. However, predictions from models generated using reconstruction methods typically used in fossil studies varied widely from high levels of quantitative accuracy to a failure to correctly capture even relative differences between macroevolutionary morphotypes. Our novel experiment emphasizes that correctly reconstructing even qualitative differences between taxa in a macroevolutionary radiation is challenging using current methods. Future studies of fossil taxa should incorporate systematic assessments of reconstruction error into their hypothesis testing and, moreover, seek to expand primary datasets on muscle properties in extant taxa to better inform soft tissue reconstructions in macroevolutionary studies.

Bite force data suggests relationship between acrodont tooth implantation and strong bite force

Extant and extinct reptiles exhibit numerous combinations of tooth implantation and attachment. Tooth implantation ranges from those possessing roots and lying within a socket (thecodonty), to teeth lying against the lingual wall of the jawbone (pleurodonty), to teeth without roots or sockets that are attached to the apex of the marginal jawbones (acrodonty). Attachment may be ligamentous (gomphosis) or via fusion (ankylosis). Generally speaking, adaptative reasonings are proposed as an underlying driver for evolutionary changes in some forms of tooth implantation and attachment. However, a substantiated adaptive hypothesis is lacking for the state of acrodont ankylosis that is seen in several lineages of Lepidosauria, a clade that is plesiomorphically pleurodont. The convergent evolution of acrodont ankylosis in several clades of lepidosaurs suggests a selective pressure shaped the evolution of the trait. We hypothesize that acrodont ankylosis as seen in Acrodonta and Sphenodon punctatus, is an adaptation either resulting from or allowing for a stronger bite force. We analyzed bite force data gathered from the literature to show that those taxa possessing acrodont dentition possess a stronger bite force on average than those taxa with pleurodont dentition. Dietary specialists with pleurodont dentition may also possess relatively high bite forces, though body size may also play a role in their ability to bite hard. Furthermore, our results have implications for the evolution of acrodont ankylosis and potential behaviors related to strong bite force that influenced the evolution of acrodonty within Acrodonta and Rhynchocephalia.

The craniomandibular anatomy of the early archosauriform Euparkeria capensis and the dawn of the archosaur skull

Archosauria (birds, crocodilians and their extinct relatives) form a major part of terrestrial ecosystems today, with over 10 000 living species, and came to dominate the land for most of the Mesozoic (over 150 Myr) after radiating following the Permian-Triassic extinction. The archosaur skull has been essential to this diversification, itself diversified into myriad forms. The archosauriform Euparkeria capensis from the Middle Triassic (Anisian) of South Africa has been of great interest since its initial description in 1913, because its anatomy shed light on the origins and early evolution of crown Archosauria and potentially approached that of the archosaur common ancestor. Euparkeria has been widely used as an outgroup in phylogenetic analyses and when investigating patterns of trait evolution among archosaurs. Although described monographically in 1965, subsequent years have seen great advances in the understanding of early archosaurs and in imaging techniques. Here, the cranium and mandible of Euparkeria are fully redescribed and documented using all fossil material and computed tomographic data. Details previously unclear are fully described, including vomerine dentition, the epiptergoid, number of premaxillary teeth and palatal arrangement. A new diagnosis and cranial and braincase reconstruction is provided, and an anatomical network analysis is performed on the skull of Euparkeria and compared with other amniotes. The modular composition of the cranium suggests a flexible skull well adapted to feeding on agile food, but with a clear tendency towards more carnivorous behaviour, placing the taxon at the interface between ancestral diapsid and crown archosaur ecomorphology, corresponding to increases in brain size, visual sensitivity, upright locomotion and metabolism around this point in archosauriform evolution. The skull of Euparkeria epitomizes a major evolutionary transition, and places crown archosaur morphology in an evolutionary context.

A high-resolution growth series of Tyrannosaurus rex obtained from multiple lines of evidence

Background

During the growth of complex multicellular organisms, chronological age, size and morphology change together in a hierarchical and coordinated pattern. Among extinct species, the growth of Tyrannosaurus rex has received repeated attention through quantitative analyses of relative maturity and chronological age. Its growth series shows an extreme transformation from shallow skulls in juveniles to deep skulls in adults along with a reduction in tooth count, and its growth curve shows that T. rex had a high growth rate in contrast to its closest relatives. However, separately, these sets of data provide an incomplete picture of the congruence between age, size, and relative maturity in this exemplar species. The goal of this work is to analyze these data sets together using cladistic analysis to produce a single hypothesis of growth that includes all of the relevant data.

Methods

The three axes of growth were analyzed together using cladistic analysis, based on a data set of 1,850 morphological characters and 44 specimens. The analysis was run in TNT v.1.5 under a New Technology search followed by a Traditional search. Correlation tests were run in IBM SPSS Statistics v. 24.0.0.0.

Results

An initial analysis that included all of the specimens recovered 50 multiple most parsimonious ontograms a series of analyses identified 13 wildcard specimens. An analysis run without the wildcard specimens recovered a single most parsimonious tree (i.e., ontogram) of 3,053 steps. The ontogram is composed of 21 growth stages, and all but the first and third are supported by unambiguously optimized synontomorphies. T. rex ontogeny can be divided into five discrete growth categories that are diagnosed by chronological age, morphology, and, in part, size (uninformative among adults). The topology shows that the transition from shallow to deep skull shape occurred between 13 and 15 years of age, and the size of the immediate relatives of T. rex was exceeded between its 15th and 18th years. Although size and maturity are congruent among juveniles and subadults, congruence is not seen among adults; for example, one of the least mature adults (RSM 2523.8) is also the largest and most massive example of the species. The extreme number of changes at the transition between juveniles and subadults shows that the ontogeny of T. rex exhibits secondary metamorphosis, analogous to the abrupt ontogenetic changes that are seen at sexual maturity among teleosts. These results provide a point of comparison for testing the congruence between maturity and chronological age, size, and mass, as well as integrating previous work on functional morphology into a rigorous ontogenetic framework. Comparison of the growth series of T. rex with those of outgroup taxa clarifies the ontogenetic trends that were inherited from the common ancestor of Archosauriformes.

The roles of joint tissues and jaw muscles in palatal biomechanics of the savannah monitor (Varanus exanthematicus) and their significance for cranial kinesis

Numerous vertebrates exhibit cranial kinesis, or movement between bones of the skull and mandible other than at the jaw joint. Many kinetic species possess a particular suite of features to accomplish this movement, including flexible cranial joints and protractor musculature. Whereas the musculoskeletal anatomy of these kinetic systems is well understood, how these joints are biomechanically loaded, how different soft tissues affect joint loading and kinetic capacity, and how the protractor musculature loads the skull remain poorly understood. Here, we present a finite element model of the savannah monitor, Varanus exanthematicus, a modestly kinetic lizard, to better elucidate the roles of soft tissue in mobile joints and protractor musculature in cranial loading. We describe the 3D resultants of jaw muscles and the histology of palatobasal, otic and jaw joints. We tested the effects of joint tissue type, bite point and muscle load to evaluate the biomechanical role of muscles on the palate and braincase. We found that the jaw muscles have significant mediolateral components that can impart stability across palatocranial joints. Articular tissues affect the magnitude of strains experienced around the palatobasal and otic joints. Without protractor muscle loading, the palate, quadrate and braincase experience higher strains, suggesting this muscle helps insulate the braincase and palatoquadrate from high loads. We found that the cross-sectional properties of the bones of V. exanthematicus are well suited for performing under torsional loads. These findings suggest that torsional loading regimes may have played a more important role in the evolution of cranial kinesis in lepidosaurs than previously appreciated.