Phylogenetically structured variance in felid bite force: the role of phylogeny in the evolution of biting performance

Evolutionary patterns of cat-like carnivorans unveil drivers of the sabertooth morphology

The sabertooth morphology stands as a classic case of convergence, manifesting recurrently across various vertebrate groups, prominently within two carnivorans clades: felids and nimravids. Nonetheless, the evolutionary mechanisms driving these recurring phenotypes remain insufficiently understood, lacking a robust phylogenetic and spatiotemporal framework. We reconstruct the tempo and mode of craniomandibular evolution of Felidae and Nimravidae and evaluate the strength of the dichotomy between conical and saber-toothed species, as well as within saber-toothed morphotypes. To do so, we investigate morphological variation, convergence, phenotypic integration, and evolutionary rates, employing a comprehensive dataset of nearly 200 3D models encompassing mandibles and crania from both extinct and extant feline-like carnivorans, spanning their entire evolutionary timeline. Our results reject the hypothesis of a distinctive sabertooth morphology, revealing instead a continuous spectrum of feline-like phenotypes in both the cranium and mandible, with sporadic instances of unequivocal convergence. Disparity peaked at the end of the Miocene and is usually higher in clades containing taxa with extreme sabertoothed adaptations. We show that taxa with saberteeth exhibit a lower degree of craniomandibular integration, allowing to exhibit a greater range of phenotypes. Those same groups usually show a burst of morphological evolutionary rate at the beginning of their evolutionary history. Consequently, we propose that a reduced degree of integration coupled with rapid evolutionary rates emerge as key components in the development of a sabertooth morphology in multiple clades.

A bite force database of 654 insect species

Bite force is a decisive performance trait in animals because it plays a role for numerous life history components such as food consumption, inter- and intraspecific interactions, and reproductive success. Bite force has been studied across a wide range of vertebrate species, but only for 32 species of insects, the most speciose animal lineage. Here we present the insect bite force database with bite force measurements for 654 insect species covering 476 genera, 111 families, and 13 orders with body lengths ranging from 3.76 to 180.12 mm. In total we recorded 1906 bite force series from 1290 specimens, and, in addition, present basal head, body, and wing metrics. As such, the database will facilitate a wide range of studies on the characteristics, predictors, and macroevolution of bite force in the largest clade of the animal kingdom and may serve as a basis to further our understanding of macroevolutionary processes in relation to bite force across all biting metazoans.

Sagittal crest morphology decoupled from relative bite performance in Pleistocene tapirs (Perissodactyla: Tapiridae)

Bite force is often associated with specific morphological features, such as sagittal crests. The presence of a pronounced sagittal crest in some tapirs (Perissodactyla: Tapiridae) was recently shown to be negatively correlated with hard-object feeding, in contrast with similar cranial structures in carnivorans. The aim of this study was to investigate bite forces and sagittal crest heights across a wide range of modern and extinct tapirs and apply a comparative investigation to establish whether these features are correlated across a broad phylogenetic scope. We examined a sample of 71 specimens representing 15 tapir species (5 extant, 10 extinct) using the dry-skull method, linear measurements of cranial features, phylogenetic reconstruction, and comparative analyses. Tapirs were found to exhibit variation in bite force and sagittal crest height across their phylogeny and between different biogeographical realms, with high-crested morphologies occurring mostly in Neotropical species. The highest bite forces within tapirs appear to be driven by estimates for the masseter-pterygoid muscle complex, rather than predicted forces for the temporalis muscle. Our results demonstrate that relative sagittal crest height is poorly correlated with relative cranial bite force, suggesting high force application is not a driver for pronounced sagittal crests in this sample. The divergent biomechanical capabilities of different contemporaneous tapirids may have allowed multiple species to occupy overlapping territories and partition resources to avoid excess competition. Bite forces in tapirs peak in Pleistocene species, independent of body size, suggesting possible dietary shifts as a potential result of climatic changes during this epoch.

Plant-eating carnivores: Multispecies analysis on factors influencing the frequency of plant occurrence in obligate carnivores

Plant-eating behavior is one of the greatest mysteries in obligate carnivores. Despite unsuitable morphological and physiological traits for plant consumption, the presence of plants in scat or stomach contents has been reported in various carnivorous species. However, researchers' interpretations of this subject are varied, and knowledge about it is scarce, without any multispecies studies. This study assessed the extent of variation in the frequency of plant occurrence in scat and stomach contents, as well as its relationship with various factors in 24 felid species using data from 213 published articles. Since the frequency of plant occurrence has not always been reported, we created two-part models and estimated parameters in a Bayesian framework. We found a significant negative relationship between the frequency of plant occurrence and body mass. This may be because plant-eating behavior reduces the energy loss caused by parasites and increases the efficiency of energy intake, which has a greater importance in smaller animals that have relatively high metabolic rates. This exploratory study highlights the importance of considering plant consumption in dietary studies on carnivorous species to understand the adaptive significance of this behavior and the relationship between obligate carnivores and plants.

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

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

Radular force performance of stylommatophoran gastropods (Mollusca) with distinct body masses

The forces exerted by the animal's food processing structures can be important parameters when studying trophic specializations to specific food spectra. Even though molluscs represent the second largest animal phylum, exhibiting an incredible biodiversity accompanied by the establishment of distinct ecological niches including the foraging on a variety of ingesta types, only few studies focused on the biomechanical performance of their feeding organs. To lay a keystone for future research in this direction, we investigated the in vivo forces exerted by the molluscan food gathering and processing structure, the radula, for five stylommatophoran species (Gastropoda). The chosen species and individuals have a similar radular morphology and motion, but as they represent different body mass classes, we were enabled to relate the forces to body mass. Radular forces were measured along two axes using force transducers which allowed us to correlate forces with the distinct phases of radular motion. A radular force quotient, AFQ = mean Absolute Force/bodymass^0.67, of 4.3 could be determined which can be used further for the prediction of forces generated in Gastropoda. Additionally, some specimens were dissected and the radular musculature mass as well as the radular mass and dimensions were documented. Our results depict the positive correlation between body mass, radular musculature mass, and exerted force. Additionally, it was clearly observed that the radular motion phases, exerting the highest forces during feeding, changed with regard to the ingesta size: all smaller gastropods rather approached the food by a horizontal, sawing-like radular motion leading to the consumption of rather small food particles, whereas larger gastropods rather pulled the ingesta in vertical direction by radula and jaw resulting in the tearing of larger pieces.

Phenotypic integration in feliform carnivores: Covariation patterns and disparity in hypercarnivores versus generalists

The skeleton is a complex arrangement of anatomical structures that covary to various degrees depending on both intrinsic and extrinsic factors. Among the Feliformia, many species are characterized by predator lifestyles providing a unique opportunity to investigate the impact of highly specialized hypercarnivorous diet on phenotypic integration and shape diversity. To do so, we compared the shape of the skull, mandible, humerus, and femur of species in relation to their feeding strategies (hypercarnivorous vs. generalist species) and prey preference (predators of small vs. large prey) using three-dimensional geometric morphometric techniques. Our results highlight different degrees of morphological integration in the Feliformia depending on the functional implication of the anatomical structure, with an overall higher covariation of structures in hypercarnivorous species. The skull and the forelimb are not integrated in generalist species, whereas they are integrated in hypercarnivores. These results can potentially be explained by the different feeding strategies of these species. Contrary to our expectations, hypercarnivores display a higher disparity for the skull than generalist species. This is probably due to the fact that a specialization toward high-meat diet could be achieved through various phenotypes. Finally, humeri and femora display shape variations depending on relative prey size preference. Large species feeding on large prey tend to have robust long bones due to higher biomechanical constraints.

Body Size and Bite Force of Stray and Feral Cats-Are Bigger or Older Cats Taking the Largest or More Difficult-to-Handle Prey?

As carnivorans rely heavily on their head and jaws for prey capture and handling, skull morphology and bite force can therefore reflect their ability to take larger or more difficult-to-handle prey. For 568 feral and stray cats (Felis catus), we recorded their demographics (sex and age), source location (feral or stray) and morphological measures (body mass, body condition); we estimated potential bite force from skull measurements for n = 268 of these cats, and quantified diet composition from stomach contents for n = 358. We compared skull measurements to estimate their bite force and determine how it varied with sex, age, body mass, body condition. Body mass had the strongest influence of bite force. In our sample, males were 36.2% heavier and had 20.0% greater estimated bite force (206.2 ± 44.7 Newtons, n = 168) than females (171.9 ± 29.3 Newtons, n = 120). However, cat age was the strongest predictor of the size of prey that they had taken, with older cats taking larger prey. The predictive power of this relationship was poor though (r² < 0.038, p < 0.003), because even small cats ate large prey and some of the largest cats ate small prey, such as invertebrates. Cats are opportunistic, generalist carnivores taking a broad range of prey. Their ability to handle larger prey increases as the cats grow, increasing their jaw strength, and improving their hunting skills, but even the smallest cats in our sample had tackled and consumed large and potentially 'dangerous' prey that would likely have put up a defence.

Dry versus wet and gross: Comparisons between the dry skull method and gross dissection in estimations of jaw muscle cross-sectional area and bite forces in sea otters

Bite force is a measure of feeding performance used to elucidate links between animal morphology, ecology, and fitness. Obtaining live individuals for in vivo bite-force measurements or freshly deceased specimens for bite force modeling is challenging for many species. Thomason's dry skull method for mammals relies solely on osteological specimens and, therefore, presents an advantageous approach that enables researchers to estimate and compare bite forces across extant and even extinct species. However, how accurately the dry skull method estimates physiological cross-sectional area (PCSA) of the jaw adductor muscles and theoretical bite force has rarely been tested. Here, we use an ontogenetic series of southern sea otters (Enhydra lutris nereis) to test the hypothesis that skeletomuscular traits estimated from the dry skull method accurately predicts test traits derived from dissection-based biomechanical modeling. Although variables from these two methods exhibited strong positive relationships across ontogeny, we found that the dry skull method overestimates PCSA of the masseter and underestimates PCSA of the temporalis. Jaw adductor in-levers for both jaw muscles and overall bite force are overestimated. Surprisingly, we reveal that sexual dimorphism in craniomandibular shape affects temporalis PCSA estimations; the dry skull method predicted female temporalis PCSA well but underestimates male temporalis PCSA across ontogeny. These results highlight the importance of accounting for sexual dimorphism and other intraspecific variation when using the dry skull method. Together, we found the dry skull method provides an underestimation of bite force over ontogeny and that the underlying anatomical components driving bite force may be misrepresented.

An assessment of the role of the falx cerebri and tentorium cerebelli in the cranium of the cat (Felis silvestris catus)

The falx cerebri and the tentorium cerebelli are two projections of the dura mater in the cranial cavity which ossify to varying degrees in some mammalian species. The idea that the ossification of these structures may be necessary to support the loads arising during feeding has been proposed and dismissed in the past, but never tested quantitatively. To address this, a biomechanical model of a domestic cat (Felis silvestris catus) skull was created and the material properties of the falx and tentorium were varied for a series of loading regimes incorporating the main masticatory and neck muscles during biting. Under these loading conditions, ossification of the falx cerebri does not have a significant impact on the stress in the cranial bones. In the case of the tentorium, however, a localized increase in stress was observed in the parietal and temporal bones, including the tympanic bulla, when a non-ossified tentorium was modelled. These effects were consistent across the different analyses, irrespective of loading regime. The results suggest that ossification of the tentorium cerebelli may play a minor role during feeding activities by decreasing the stress in the back of the skull.

Directional selection in the evolution of elongated upper canines in clouded leopards and sabre-toothed cats

Extremely developed or specialized traits such as the elongated upper canines of extinct sabre-toothed cats are often not analogous to those of any extant species, which limits our understanding of their evolutionary cause. However, an extant species may have undergone directional selection for a similar extreme phenotype. Among living felids, the clouded leopard, Neofelis nebulosa, has exceptionally long upper canines for its body size. We hypothesized that directional selection generated the elongated upper canines of clouded leopards in a manner similar to the process in extinct sabre-toothed cats. To test this, we developed an approach that compared the effect of directional selection among lineages in a phylogeny using a simulation of trait evolution and approximate Bayesian computation. This approach was applied to analyse the evolution of upper canine length in the Felidae phylogeny. Our analyses consistently showed directional selection favouring longer upper canines in the clouded leopard lineage and a lineage leading to the sabre-toothed cat with the longest upper canines, Smilodon. Most of our analyses detected an effect of directional selection for longer upper canines in the lineage leading to another sabre-toothed cat, Homotherium, although this selection may have occurred exclusively in the primitive species. In all the analyses, the clouded leopard and Smilodon lineages showed comparable directional selection. This implies that clouded leopards share a selection advantage with sabre-toothed cats in having elongated upper canines.

An integrative method for testing form-function linkages and reconstructed evolutionary pathways of masticatory specialization

Morphology serves as a ubiquitous proxy in macroevolutionary studies to identify potential adaptive processes and patterns. Inferences of functional significance of phenotypes or their evolution are overwhelmingly based on data from living taxa. Yet, correspondence between form and function has been tested in only a few model species, and those linkages are highly complex. The lack of explicit methodologies to integrate form and function analyses within a deep-time and phylogenetic context weakens inferences of adaptive morphological evolution, by invoking but not testing form-function linkages. Here, we provide a novel approach to test mechanical properties at reconstructed ancestral nodes/taxa and the strength and direction of evolutionary pathways in feeding biomechanics, in a case study of carnivorous mammals. Using biomechanical profile comparisons that provide functional signals for the separation of feeding morphologies, we demonstrate, using experimental optimization criteria on estimation of strength and direction of functional changes on a phylogeny, that convergence in mechanical properties and degree of evolutionary optimization can be decoupled. This integrative approach is broadly applicable to other clades, by using quantitative data and model-based tests to evaluate interpretations of function from morphology and functional explanations for observed macroevolutionary pathways.

Accuracy of mandibular force profiles for bite force estimation and feeding behavior reconstruction in extant and extinct carnivorans

Mandibular force profiles apply the principles of beam theory to identify mandibular biomechanical properties that reflect the bite force and feeding strategies of extant and extinct predators. While this method uses external dimensions of the mandibular corpus to determine its biomechanical properties, more accurate results could potentially be obtained by quantifying its internal cortical bone distribution. To test this possibility, mandibular force profiles were calculated using both external mandibular dimensions (‘solid mandible model’) and quantification of internal bone distribution of the mandibular corpus obtained from CT scans (‘hollow mandible model’) for five carnivorans (Canis lupus, Crocuta crocuta, Panthera leo, Neofelis nebulosa, and the extinct Canis dirus). Comparison reveals that the solid model slightly overestimates mandibular biomechanical properties, but the pattern of change in biomechanical properties along the mandible remains the same. As such, feeding behavior reconstructions are consistent between the two models and are not improved by computed tomography. Bite force estimates produced by the two models are similar, except for Crocuta where the solid model underestimates bite force by 10%-14%. This discrepancy is due to the more solid nature of the Crocuta mandible relative to other carnivorans. Therefore, computed tomography improves bite force estimation accuracy for taxa with thicker mandibular corpora, but not significantly so otherwise. Bite force estimates derived from mandibular force profiles are far closer to empirically-measured bite force than those inferred from jaw musculature dimension. Consequently, bite force estimates derived from this method can be used to calibrate finite-element analysis models.

Constructional morphology within the head of hammerhead sharks (sphyrnidae)

The study of functional trade-offs is important if a structure, such as the cranium, serves multiple biological roles, and is, therefore, shaped by multiple selective pressures. The sphyrnid cephalofoil presents an excellent model for investigating potential trade-offs among sensory, neural, and feeding structures. In this study, hammerhead shark species were chosen to represent differences in head form through phylogeny. A combination of surface-based geometric morphometrics, computed tomography (CT) volumetric analysis, and phylogenetic analyses were utilized to investigate potential trade-offs within the head. Hammerhead sharks display a diversity of cranial morphologies where the position of the eyes and nares vary among species, with only minor changes in shape, position, and volume of the feeding apparatus through phylogeny. The basal winghead shark, Eusphyra blochii, has small anteriorly positioned eyes. Through phylogeny, the relative size and position of the eyes change, such that derived species have larger, more medially positioned eyes. The lateral position of the external nares is highly variable, showing no phylogenetic trend. Mouth size and position are conserved, remaining relatively unchanged. Volumetric CT analyses reveal no trade-offs between the feeding apparatus and the remaining cranial structures. The few trade-offs were isolated to the nasal capsule volume's inverse correlation with braincase, chondrocranial, and total cephalofoil volume. Eye volume also decreased as cephalofoil width increased. These data indicate that despite considerable changes in head shape, much of the head is morphologically conserved through sphyrnid phylogeny, particularly the jaw cartilages and their associated feeding muscles, with shape change and morphological trade-offs being primarily confined to the lateral wings of the cephalofoil and their associated sensory structures.

Jaw function in Smilodon fatalis: a reevaluation of the canine shear-bite and a proposal for a new forelimb-powered Class 1 Lever Model

The jaw function of Smilodon fatalis has long been a source of debate. Although modern-day lions subdue large prey through the use of a suffocating throat bite, the dramatically elongated maxillary canines of S. fatalis suggest an alternative bite mechanism. The current literature favors a "canine shear-bite," in which the depression of the cranium by the ventral neck flexors assists the mandibular adductors in closing the jaws. Although the model makes intuitive sense and appears to be supported by scientific data, the mechanical feasibility of "neck-powered" biting has not been experimentally demonstrated. In the present study, the computer-assisted manipulation of digitized images of a high-quality replica of an S. fatalis neck and skull shows that a rotation of the cranium by the ventral neck flexors will not result in jaw closure. Instead, the cranium and mandible rotate ventrally together (at the atlantooccipital joint), and the jaws remain in an open configuration. The only manner by which rotation of the cranium can simultaneously result in jaw closure is by an anterior rotation at the temporomandibular joint. Based on this finding, the author proposes a new Class 1 lever mechanism for S. fatalis jaw function. In this model, the mandible is immobilized against the neck of the prey and a dorsally directed force from the extension of the forelimbs rotates the cranium anteriorly at the temporomandibular joint. The maxillary canines pierce the prey's neck and assist in clamping the ventral neck structures. The model is based on a maximum gape angle of approximately 90° and incorporates a secondary virtual point of rotation located slightly anteroventral to the temporomandibular joint. The Class 1 Lever Model is mechanically feasible, consistent with current data on S. fatalis anatomy and ecology, and may provide a basis for similar studies on other fossil taxa.

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

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

Cranial myology and bite force performance of Erlikosaurus andrewsi: a novel approach for digital muscle reconstructions

The estimation of bite force and bite performance in fossil and extinct animals is a challenging subject in palaeontology and is highly dependent on the reconstruction of the cranial myology. Furthermore, the morphology and arrangement of the adductor muscles considerably affect feeding processes and mastication and thus also have important dietary and ecological ramifications. However, in the past, the reconstruction of the (cranial) muscles was restricted to the identification of muscle attachment sites or simplified computer models. This study presents a detailed reconstruction of the adductor musculature of the Cretaceous therizinosaur Erlikosaurus andrewsi based on a stepwise and iterative approach. The detailed, three-dimensional models of the individual muscles allow for more accurate measurements of the muscle properties (length, cross-section, attachment angle and volume), from which muscle and bite force estimates are calculated. Bite force estimations are found to be the lowest at the tip of the snout (43-65 N) and respectively higher at the first (59-88 N) and last tooth (90-134 N) position. Nevertheless, bite forces are comparatively low for E. andrewsi, both in actual numbers as well as in comparison with other theropod dinosaurs. The results further indicate that the low bite performance was mainly used for leaf-stripping and plant cropping, rather than active mastication or chewing processes. Muscle and thus bite force in E. andrewsi (and most likely all therizinosaurs) is considerably constrained by the cranial anatomy and declines in derived taxa of this clade. This trend is reflected in the changes of dietary preferences from carnivory to herbivory in therizinosaurs.

Convergence and divergence in the evolution of cat skulls: temporal and spatial patterns of morphological diversity

BACKGROUND

Studies of biological shape evolution are greatly enhanced when framed in a phylogenetic perspective. Inclusion of fossils amplifies the scope of macroevolutionary research, offers a deep-time perspective on tempo and mode of radiations, and elucidates life-trait changes. We explore the evolution of skull shape in felids (cats) through morphometric analyses of linear variables, phylogenetic comparative methods, and a new cladistic study of saber-toothed cats.

METHODOLOGY/PRINCIPAL FINDINGS

A new phylogenetic analysis supports the monophyly of saber-toothed cats (Machairodontinae) exclusive of Felinae and some basal felids, but does not support the monophyly of various saber-toothed tribes and genera. We quantified skull shape variation in 34 extant and 18 extinct species using size-adjusted linear variables. These distinguish taxonomic group membership with high accuracy. Patterns of morphospace occupation are consistent with previous analyses, for example, in showing a size gradient along the primary axis of shape variation and a separation between large and small-medium cats. By combining the new phylogeny with a molecular tree of extant Felinae, we built a chronophylomorphospace (a phylogeny superimposed onto a two-dimensional morphospace through time). The evolutionary history of cats was characterized by two major episodes of morphological divergence, one marking the separation between saber-toothed and modern cats, the other marking the split between large and small-medium cats.

CONCLUSIONS/SIGNIFICANCE

Ancestors of large cats in the 'Panthera' lineage tend to occupy, at a much later stage, morphospace regions previously occupied by saber-toothed cats. The latter radiated out into new morphospace regions peripheral to those of extant large cats. The separation between large and small-medium cats was marked by considerable morphologically divergent trajectories early in feline evolution. A chronophylomorphospace has wider applications in reconstructing temporal transitions across two-dimensional trait spaces, can be used in ecophenotypical and functional diversity studies, and may reveal novel patterns of morphospace occupation.

Exploring patterns of interspecific variation in quantitative traits using sequential phylogenetic eigenvector regressions

A number of metrics have been developed for estimating phylogenetic signal in data and to evaluate correlated evolution, inferring broad-scale evolutionary and ecological processes. Here, we proposed an approach called phylogenetic signal-representation (PSR) curve, built upon phylogenetic eigenvector regression (PVR). In PVR, selected eigenvectors extracted from a phylogenetic distance matrix are used to model interspecific variation. In the PSR curve, sequential PVR models are fitted after successively increasing the number of eigenvectors and plotting their R(2) against the accumulated eigenvalues. We used simulations to show that a linear PSR curve is expected under Brownian motion and that its shape changes under alternative evolutionary models. The PSR area, expressing deviations from Brownian motion, is strongly correlated (r= 0.873; P < 0.01) with Blomberg's K-statistics, so nonlinear PSR curves reveal if traits are evolving at a slower or higher rate than expected by Brownian motion. The PSR area is also correlated with phylogenetic half-life under an Ornstein-Uhlenbeck process, suggesting how both methods describe the shape of the relationship between interspecific variation and time since divergence among species. The PSR curve provides an elegant exploratory method to understand deviations from Brownian motion, in terms of acceleration or deceleration of evolutionary rates occurring at large or small phylogenetic distances.

Jaw biomechanics and the evolution of biting performance in theropod dinosaurs

Despite the great diversity in theropod craniomandibular morphology, the presence and distribution of biting function types across Theropoda has rarely been assessed. A novel method of biomechanical profiling using mechanical advantage computed for each biting position along the entirety of the tooth row was applied to 41 extinct theropod taxa. Multivariate ordination on the polynomial coefficients of the profiles reveals the distribution of theropod biting performance in function space. In particular, coelophysoids are found to occupy a unique region of function space, while tetanurans have a wide but continuous function space distribution. Further, the underlying phylogenetic structure and evolution of biting performance were investigated using phylogenetic comparative methods. There is a strong phylogenetic signal in theropod biomechanical profiles, indicating that evolution of biting performance does not depart from Brownian motion evolution. Reconstructions of ancestral function space occupation conform to this pattern, but phylogenetically unexpected major shifts in function space occupation can be observed at the origins of some clades. However, uncertainties surround ancestor estimates in some of these internal nodes, so inferences on the nature of these evolutionary changes must be viewed with caution.