Extreme and rapid bursts of functional adaptations shape bite force in amniotes
Abstract
Adaptation is the fundamental driver of functional and biomechanical evolution. Accordingly, the states of biomechanical traits (absolute or relative trait values) have long been used as proxies for adaptations in response to direct selection. However, ignoring evolutionary history, in particular ancestry, passage of time and the rate of evolution, can be misleading. Here, we apply a recently developed phylogenetic statistical approach using significant rate shifts to detect instances of exceptional rates of adaptive changes in bite force in a large group of terrestrial vertebrates, the amniotes. Our results show that bite force in amniotes evolved through multiple bursts of exceptional rates of adaptive changes, whereby whole groups—including Darwin's finches, maniraptoran dinosaurs (group of non-avian dinosaurs including birds), anthropoids and hominins (fossil and modern humans)—experienced significant rate increases compared to the background rate. However, in most parts of the amniote tree of life, we find no exceptional rate increases, indicating that coevolution with body size was primarily responsible for the patterns observed in bite force. Our approach represents a template for future studies in functional morphology and biomechanics, where exceptional rates of adaptive changes can be quantified and potentially linked to specific ecological factors underpinning major evolutionary radiations.
1. Background
Adaptation is the fundamental driver of functional and biomechanical evolution. Measures of biomechanical performance—e.g. bite force—characterize and quantify specific functional performances to fulfil ecological demand—e.g. diet [1,2]. Functions are therefore typically assumed to be under direct selection—i.e. a change in biomechanical performance indicates selection for changes in function. For instance, taxa with higher bite force are often interpreted as having diets requiring powerful bites, driving selection on associated morphological features [3–5]. However, observable trait states can be the result of multiple confounding factors independent of a direct selection for a specific function, including scaling effects [6] and evolutionary history (ancestry and the passage of time). That is, species can often achieve relatively high bite forces owing to the accumulation of evolutionary changes that are within expectation given the passage of time, without invoking exceptionally strong selection pressures towards increased bite forces. Assuming that the states of biomechanical traits can be used as proxies for strong selection/functional demand while ignoring such confounding factors (in particular evolutionary history) can be, and often are, misleading [7,8]. Thus, it is important to identify cases in which trait values exceed the amount of evolutionary change expected given confounding factors—i.e. objectively identifying instances of exceptional rates of adaptive changes [9].
It has long been interpreted that the intensity of natural selection acting on a phenotype is linked to rates of evolution [9,10], and significant shifts in rates can be interpreted as instances of adaptive responses to strong selective pressures—e.g. positive phenotypic selection [9] akin to positive genetic selection [11,12]. Positive selection is invoked as an explanation for trait evolution along a branch on a phylogenetic tree if the amount of evolutionary change exceeds the amount of change expected from the passage of time given the background rate of evolution. In the context of biomechanical evolution, evidence for exceptional rates of adaptive changes in biomechanical traits can be detected as significant rate shifts using this phylogenetic comparative framework. However, in spite of the ever-increasing numbers of comparative biomechanical studies employing a phylogenetic framework [13–27], those explicitly interpreting significant rate shifts in biomechanical trait evolution as responses to strong levels of natural selection are rare (but see, e.g. [23]).
Here, we test the hypothesis that instances of exceptional rates of adaptive changes have shaped the observed diversity in biomechanical traits focusing on bite force evolution. Our hypothesis specifically concerns exceptional rates of adaptive changes and not rate heterogeneity—exceptional rates of adaptation is defined as instances of exceptionally large rate increases (at least twice the background rate, see Materials and methods and [9] for details). Bite force relates to species' niche and feeding ecology, is correlated with several ecological and behavioural traits [28–31] and is widely available from the literature across several fields of study (e.g. biomechanics, ecology, palaeobiology) for a broad sample of the tree of life, making it an ideal biomechanical trait to test our hypothesis using a phylogenetic framework. To this end, we assembled, to the best of our knowledge, the largest dataset of bite forces collected to date for amniotes, both extinct (including non-avian dinosaurs, sabre-toothed cats and fossil hominins) and extant.
2. Material and methods
(a) Data
Bite force and body mass data were collected primarily through the literature, augmented with novel estimates (electronic supplementary material), spanning 434 extant and extinct amniote species (electronic supplementary material, table S1). We used the Time Tree of Life (TTOL) [32] as the backbone phylogeny, with fossil tips/clades inserted at the appropriate phylogenetic and temporal positions (electronic supplementary material) using fossil dates from the Palaeobiology Database (accessed 9 February 2017).
(b) Variable rates phylogenetic regression models
We fitted a variable rates (VR) regression model [9] using BayesTraits [33] on log10 bite force against log10 body mass (single-slope VR model). The VR regression model as implemented in BayesTraits works to modify the branch lengths to ‘detect heterogeneity in rates of phylogenetically structured residual errors’ [9]. That is, once the appropriate level of variance in the response variable—e.g. bite force—is explained by some predictor variable(s)—e.g. body mass—outlying deviations from the regression line will be explained as rate shifts (electronic supplementary material, figure S1). Under Brownian motion, bite force—after accounting for body mass and other confounding variables—evolves at a rate proportional to time (and an estimated background variance) across the phylogeny, and for any evolutionary change along a given branch that is greater/less than the expected amount of change for the duration of time to occur (given body mass), that branch must be stretched/compressed in length in proportion to the observed amount of phenotypic change—corresponding to a rate increase/decrease. The magnitude of branch stretching/compressing is the rate scalar (r).
As previous research indicated that the scaling of bite force is group-specific [34–36], we tested an additional model in which separate slopes were estimated for five different groups (5-group VR model). Each taxon was assigned to one of five groups: mammals excluding bats (hereafter ‘mammals’), bats, finches (Fringillidae, Estrildidae and Darwin's finches); dinosaurs including birds but excluding finches (hereafter ‘dinosaurs’); and non-dinosaurian diapsid reptiles (hereafter ‘reptiles’) (figure 1; electronic supplementary material, table S1). We chose these five groups because they act as good descriptors of the distribution of data (electronic supplementary material, figure S3) as well as conforming to widely recognized taxonomic groups (see the electronic supplementary material, for details). Birds and dinosaurs were grouped together, as fossil evidence points to a blurred distinction in their physiology and biology [37–40], while dinosaurs are very different from other reptiles [41]. Intercept differences were not modelled—groupwise offsets will trade off with rate scalars on branches leading to the last common ancestors of the respective groups (electronic supplementary material, figure S1). In order to account for potential differences owing to bite force acquisition type and biting positions (electronic supplementary material), two additional models were fitted with the confounding variables ‘bite type’ and ‘bite point’ individually added to the regression model as covariates (electronic supplementary material).
We tested for instances of exceptional bursts of evolutionary change in bite force based on rate shifts along branches on the phylogeny [9], based on the premise that phenotypic changes owing to adaptations (potentially as a response to strong selective pressure) would be proportional to r. Thus, we define exceptional change following the criteria of [9]: (i) certainty of rate shifts, the branch in question must be scaled in greater than 95% of the posterior sample of scaled trees; and (ii) magnitude of rate shifts, the r in question must be greater than two. Rate heterogeneity that does not fulfil these two criteria was not considered as instances of exceptional rates of adaptive changes. There are two types of rate shifts: branchwise rate shifts (branch shifts), rate shifts that are detected along individual branches; and cladewise rate shifts (clade shifts), rate shifts that occur across all branches within a clade and represent cases in which rapid divergences in trait values have occurred (electronic supplementary material, figure S1). We determined whether rate shifts constituted exceptional rates of adaptive changes if they satisfied the criteria set out above in all of three independent replicate Markov chain Monte Carlo (MCMC) chains.
In order to determine if rate heterogeneity was statistically significant, we fitted an equal-rate (ER) model (or Brownian motion) as a simpler alternative to each of our VR models. Model selection was performed using the Bayes factor (BF; see the electronic supplementary material, for details): BF is defined as twice the difference in log marginal likelihood (m) between the complex model (model1) and the simple model (model0)—i.e. BF = 2 × (m1 – m0). For instance, we computed BF using m from our 5-group VR model and the simple alternative 5-group ER model and selected the VR model over the ER model when BF value was greater than 2 [42].
We ran our MCMC chains for 109 iterations, with a burn-in period of 108 iterations, sampling every 105 iteration, resulting in a posterior sample of 900 modified VR trees and model estimates, for each regression model. We used stepping stone sampling (over 1000 stones at 105 iterations each) to compute marginal likelihoods from which BF were calculated. Post-processing of the BayesTraits outputs were conducted using an online post-processor (www.evolution.reading.ac.uk/VarRatesWebPP), as well as in R [43].
3. Results
(a) Variable rates regression model
We found strong support for the VR model compared to the ER model (BFVR−ER = 474) for the single-slope regression model. Bite force scales nearly isometrically with body mass, with a slope of 0.674 (pMCMC0 < 0.001, R2mean = 0.79; electronic supplementary material, table S2), which is not significantly different from a theoretical isometric slope of 0.67 (pMCMC0.67 = 0.4) [6].
There is statistical support for favouring the 5-group model over the single-slope model (figure 1); significant differences exist among the slopes of different groups (electronic supplementary material, tables S2 and S3). Finches and bats are not different from each other, but are distinct from mammals, reptiles and dinosaurs (figure 1; electronic supplementary material, table S3); in turn, these three groups are not different from each other (figure 1; electronic supplementary material, table S3). Finches and bats have slopes that deviate from 0.67 (electronic supplementary material, table S4), while the other three groups have slopes that are not significantly different from 0.67 (electronic supplementary material, table S4). Critically, despite allowing for the variation in slopes among taxonomic groups, our 5-group VR model (R2mean = 0.809; electronic supplementary material, table S2) still outperforms a 5-group ER model (BFVR−ER = 429; electronic supplementary material, table S2).
Model selection showed that bite type is not significant (pMCMC greater than 0.05 in both single-slope + bite type and 5-group + bite type models; electronic supplementary material, table S2) while bite point is (pMCMC less than 0.05 in both single-slope + bite point and 5-group + bite point models; electronic supplementary material, table S2). Given a similar body size, bite force is comparable in magnitude between in vivo measurements and indirect estimates, but it differs in magnitude between posterior and anterior bites (posterior positions have higher forces, as expected). There is no slope difference between bite type categories or between bite point categories (electronic supplementary material, figures S6 and S7). We used the 5-group + bite point model as our final model for detecting exceptional rates of adaptive changes in bite force—this enables us to compare evolutionary rates after accounting for effects owing to body size and bite point.
(b) Rate shifts and exceptional rates of adaptive changes
We found a substantial amount of rate heterogeneity (elevated rates in greater than 50% of the posterior sample) in the amniote tree of life, along 439 branches out of 866 branches in the phylogeny (51% of branches; electronic supplementary material, figure S4). Instances of exceptional rates of adaptive changes are found in a far fewer number of branches: in 182 branches (21% of branches) (figures 1 and 2a; electronic supplementary material, table S5 and movie S1). Our results show that bite force evolved through multiple bursts of exceptional rates of adaptive changes, whereby whole groups experienced rate increases of bite force evolution compared to the background rate across the entire amniote tree. We find such clade shifts in Darwin's finches (median r > 55), the hominin lineage excluding Australopithecus anamensis (r > 35), Anthropoidea (r > 6) and maniraptoran theropod dinosaurs (the clade including birds and their closest relatives, here Erlikosaurus and Dromaeosaurus; median r > 3). Thus, most of the scaled branches are because of a smaller number of node scalars that modify all descendant branches.
Two aspects of the distribution of rate shifts are noteworthy. First, cladewise rate shifts show a nested pattern (figure 1). For example, Darwin's finches exhibit an additional level of rate increase above that of the maniraptoran rate increase (figures 1 and 2a,b; electronic supplementary material, figure S5 and movies S1 and S2). The same pattern characterizes hominins (excluding A. anamensis) within anthropoid primates (figures 1 and 2a,c; electronic supplementary material, figure S5 and movies S1 and S2). In particular, Homo species exhibit reductions in bite forces (figure 2c), which is in marked contrast to the apparent increase in hominin body size through time (electronic supplementary material, figures S8) [44]. Thus, humans drastically reduced bite force through time at a rate faster than their anthropoid ancestors and relatives.
The second key aspect of the rate shift distribution is that the branchwise rate increases occur in conjunction with cladewide shifts. We identified an exceptional increase in the rate of bite force evolution along the branch leading to Passeroidea (the clade defined by the last common ancestor of our finches) (median r > 30), followed by a reversal to the ancestral maniraptoran rate. We also recovered a scattering of branch-specific shifts on terminal branches (Proteles cristatus, Panthera onca, Sus scrofa, Stegosaurus and Plateosaurus), marking sudden changes in the biomechanical performance of some species from their close relatives (figures 1 and 2a; electronic supplementary material, figure S5 and movie S1).
4. Discussion
(a) Exceptional rates of adaptive changes in bite force
Working under the premise that the rate of phenotypic trait evolution is proportional to the strength of selection [9,10], we provide, to our knowledge, the first evidence for exceptional rates of adaptive changes shaping the diversity of bite force in both extinct and extant amniotes, using a statistically robust evolutionary framework [9]. We find such instances of exceptional rates of adaptive changes in four clades of amniotes and a handful of independent branches. Conversely, in most lineages of amniotes studied here (79% of branches in the phylogeny), bite force does not undergo exceptional rates of adaptive changes (even in lineages with elevated rates; electronic supplementary material, figure S4), indicating that coevolution with body size is, for the most part, the main factor responsible for bite force variation. A large predator can generate enough bite force to kill its prey just by being large. As an example, Tyrannosaurus rex was most likely capable of ‘pulverizing’ bones [45] simply owing to its colossal size (approx. 5–10 tonnes [46–48]). We did not detect instances of exceptional rates of adaptive changes in bite force in this taxon, and therefore there is no evidence of strong selection on disproportionately high bite force. Similarly, we do not detect signatures of exceptional rates of adaptive changes in classically recognized power-biters such as osteophagous hyenids [49], short-faced hyper-carnivorous felids [50,51] and small-brained carnivorous marsupials [3,5,52], indicating that bite force in these clades was not subjected to strong selection as is often presumed.
Interestingly, contrary to our prior expectations, we do not detect exceptional rates of adaptive changes in sabre-toothed cats (Machairodontinae). Although rates are on average higher in Felidae as a whole (including both sabre- and conical-toothed cats along with the basal cats Proailurus and Hyperailurictis) compared to the background rate in the majority of the posterior sample (i.e. greater than 50%; electronic supplementary material, figure S4), they do not fulfil the criteria for exceptional rates of adaptive changes. Further, there is no difference between rates in machairodontine lineages and other felid lineages. The time elapsed in the lineage leading to sabre-toothed cats since their divergence with conical-toothed cats sufficiently explains the reduction in bite force in sabre-toothed cats.
It is possible that for the majority of taxa in our bite force sample, individuals preferentially seek out and consume food items that can be processed within the naturally generated range of bite forces, and rarely actively seek food items that require maximum biting capacity. This equates to a behavioural adaptation, in which species evolve bite force through correlated evolution with body mass, and preferentially consume food items that fit within their natural range of bite force. If this is true, then selection for improved biting performance with respect to higher bite force may not frequently occur. Additionally, functional–morphological adaptations—e.g. tooth morphology—may facilitate higher biting performance—e.g. piercing, crushing or shearing—without necessitating a more powerful bite.
Similarly, a drastic reduction in bite force—as a trade-off between force and velocity if jaw closing velocity was under strong selection—would constitute an exceptional rate of adaptation significantly below the expected range of bite force given the universal scaling relationship. It is potentially more likely for reductions in bite force to occur over exceptional gains in bite force because muscles are expensive organs to maintain and more so to enlarge. Indeed, we observe more instances of drastic reductions in bite force associated with exceptional rates of adaptive changes than we do exceptional increases in bite force (figure 2; electronic supplementary material, figure S5).
Despite the overall uniform evolution of bite force relative to body mass, we find evidence for exceptional rates of adaptive changes playing a major role in the evolution of bite force (just not as a tree-wide pattern across all major clades). Since the majority of these exceptional rates of adaptive changes occur as cladewide rate shifts, it is possible that they are linked to some biological, ecological or environmental features unique to those clades and shared among constituent members. For instance: the acquisition of a ‘key’ innovation, which allows such clades to rapidly expand and exploit functional niches; a shift into a new environment, habitat or lifestyle that is associated with new opportunities and resources; or an extrinsic environmental event (such as mass extinction events) that results in an abundance of ecological niches available for exploitation. Determining such factors ultimately responsible for exceptional rates of adaptive changes is theoretically possible. Namely, the VR regression framework allows for the inclusion of extrinsic factors such as dietary preference, feeding strategy, sexual display/conflict, etc. as additional covariates. A covariate can be identified as the extrinsic driver of bite force evolution if its inclusion can explain much of the variation in bite force, thereby reducing or eliminating rate shifts. At present, ecological data associated with biting performance are only available for a handful of species, but we hope that future work will considerably augment information on ecological covariates.
(b) Macroevolutionary patterns of exceptional rates of adaptive changes
Overall, our results highlight a combination of cladewise and branchwise rate shifts occurring across the amniote phylogeny. Cladewise rate shifts are characterized by an elevated rate that is homogeneous across all branches within a given clade and are associated with an increase in trait variation in the constituent taxa given the variance in the other taxa in the data [9]. Lineages in such clades continually evolve traits at a faster rate through time compared to other parts of the tree (figure 1). This contrasts with a classic description of an adaptive radiation [53], which is characterized by a rapid initial burst of trait evolution followed by a rate slowdown associated with niche saturation [54]; clade-wise rate shifts see continual changes in functional niche occupation. One implication of such patterns in bite force evolution is that evolutionary lineages (sequence of branches leading to terminal taxa) do not stay in the same regions of function-space (pertaining to biting functional variation), but rather, continue to expand out to unoccupied/unexplored regions of function-space. Evolutionary lineages will be moving through various functional niches as their bite force values change through time. An alternative interpretation is that function-space itself changes through time, i.e. functional/ecological niches are dynamic rather than fixed entities, a constantly moving target [55]. Yet another interpretation is that function-space saturates but convergences occur frequently and rapidly—that is, lineages move in and out of occupied/explored regions of function-space. In this context, our results would support the notion that functional adaptations are relatively labile over evolutionary history and remain responsive to changing environmental conditions and ecological demands.
Branchwise rate shifts (rate shifts associated with single branches only) occurring on branches subtending whole clades (figure 1), such as that observed at the base of Passeroidea in our dataset, can be interpreted as a mean-shift in bite force after accounting for body mass [9]. In our case, this means that there was a rapid shift in the mean bite force value of Passeroidea from the ancestral maniraptoran mean. The total sum of evolutionary changes accumulating along the branch leading to Passeroidea exceeds that expected from the temporal duration of that branch. This is irrespective of any unsampled taxa along the lineage—e.g. other perching birds (Passeriformes, e.g. corvids, shrikes) for which bite force data are not available in the literature as far as we are aware.
Similarly, the two large-bodied herbivorous taxa, Stegosaurus and Plateosaurus, have evolved bite force at excessively high rates (approx. 11 and 35 times background rate, respectively—Stegosaurus and Plateosaurus have extremely small heads, and thus low bite forces, for their body sizes), but these could potentially represent evolutionary patterns within thyreophoran and sauropodomorph dinosaurs respectively, and not specifically associated with these two species. Using a different measure of size such as head length or width may possibly change these results—though body mass has major benefits over head size for its ecological implications. Nonetheless, major changes in bite force relative to body mass did occur along these lineages so the interpretation remains the same: the amount of change in trait values given the duration of time elapsed is exceptionally high compared to the background rate.
(c) Evolution of bite force in maniraptoran dinosaurs
Maniraptoran theropods are perhaps the most diverse among dinosaurs in terms of functional and morphological specializations associated with feeding. Forms like the parrot-like oviraptorosaurs, large herbivorous therizinosaurs, hyper-carnivorous dromaeosaurs with recurved teeth (e.g. Velociraptor), and toothed and toothless avialans are just some typical examples of maniraptoran morpho-functional diversity. High evolutionary rates in maniraptoran bite force indicate that their morphological and presumed ecological diversity are linked with selection on biting performance. Maniraptoran fossils are predominantly known from Cretaceous rocks but are inferred to have originated by the Middle Jurassic (approx. 168 Myr ago; figure 1), with derived members including the avialan Archaeopteryx appearing relatively quickly, by the Late Jurassic (approx. 150 Myr ago; figure 1). This implies that Maniraptora underwent a rapid diversification (both in species diversity but also in bite force variance) early in their evolutionary history, but that they retained high evolutionary rates in bite force throughout the clade's history.
The observation that bite force underwent exceptional rates of adaptive changes in maniraptoran theropods but not uniquely in birds—rates in birds are not distinguishable from those in other non-avian maniraptorans—is consistent with recent findings that the evolution of birds and their immediate close relatives—i.e. paravians—are similar to one another [37], and that many of the features traditionally associated with birds were present in paravians and more broadly in maniraptorans. Here we have demonstrated that this is also the case with bite force evolution (given the available data); the rate of bite force evolution did not change from non-avian maniraptorans to birds. On the other hand, this means that heritable rates of bite force evolution in maniraptoran ancestors possibly contributed to some extent on the subsequent ecological success of birds—the ability to rapidly change bite force in response to changing environmental and ecological pressures would surely have been beneficial for early Cenozoic birds in the post-extinction world.
Our identification of an extreme cladewide rate shift in the Darwin's finches, which is among the highest in the tree (greater than 55 times the background rate; figures 1 and 2; electronic supplementary material, figure S5, movies S1 and S2, table S5), is noteworthy for both historical and biological reasons. Darwin's finches are the classic textbook case of ‘adaptive radiation’, with eco-morphological diversification occurring in a short time interval after the initial colonization of the Galapagos Islands by finches [56–58]. Their diversification in feeding ecology is particularly relevant to the rapid evolution of bite force, as Darwin's finches are well documented to have strong dietary preferences on food types of varying toughness [59] or differences in food manipulations [60]. Within the context of the evolutionary history of amniotes as sampled here (approx. 350 Myr), the radiation of Darwin's finches is comparatively recent with some divergences occurring in a geologically instantaneous manner (figure 1). Compared to their recent divergence times, bite force variance in Darwin's finches is exceptionally high, spanning almost two orders of magnitude (figure 2) accounting for their extraordinarily high evolutionary rates.
(d) Evolution of bite force in humans
The exceptionally high rates of bite force evolution in the hominins excluding A. anamensis (figures 1 and 2)—more than 35 times the background rate and approximately six times those for the branches within other anthropoid primates—highlights an important, recent and rapid evolution in our own lineage. In particular, the decrease in bite force in Homo species (figure 2) is contrary to the increase in hominin body size through time (electronic supplementary material, figure S8)—such a discrepancy between bite force and body mass is indicative of strong directional selection, and coincides with previously documented evolutionary shifts in relative molar sizes, attributed to the reduction in feeding time associated with the introduction of food processing such as cooking [9,25]. Strikingly, Pampush et al. [61] found a rate shift in symphysial angle in the hominin lineage excluding A. anamensis, coinciding with the phylogenetic location where we find a rate shift in bite force (figure 1).
Additionally, the reduction in bite force in the hominin lineage may have occurred as a consequence of an evolutionary trade-off with increasing brain size in this group [62–64] (see the electronic supplementary material, figures S9 and S10). As brain size increases relative to skull size, the temporal fossa (defined as the opening between the braincase and the zygomatic arch) is reduced in dimension, thereby decreasing the amount of space available to house the temporal muscles [63,64], which are critical for achieving hard biting in most animals. This reduction in temporal muscle can be seen in the changing predominance of the sagittal crest (ridge of bone running along the midline of the skull) where the temporal muscles attach through hominin evolution (electronic supplementary material, figure S10). Strikingly, molecular evidence supports the hypothesis that a drastic reduction in the temporal muscle occurred along the lineage to H. sapiens after its divergence from chimpanzees owing to a frameshifting mutation, causing inactivation of the predominant myosin heavy chain expression in masticatory muscles [63]. The mutation has been inferred to have coincided approximately with the enlargement of the brain, presumably concurrent with the origin of the genus Homo [63]—though the timing has been contested [65,66]. Furthermore, H. sapiens relies less on the temporal muscles and more on the masseter muscles to generate bite force [67]. Thus, such a trade-off between temporal muscle size and brain size is a reasonable explanation for the evolutionary reduction in bite force in the hominin lineage through time, with the advent of cooking further accelerating the loss of reliance on high bite force for food processing.
Indeed, an auxiliary phylogenetic regression modelling [68] of bite force on brain size (endocranial volume) [69–71] accounting for body mass on hominids (Pongo, Gorilla, Pan, and hominins), shows that bite force scales negatively with brain size (electronic supplementary material, figure S9). The reduction in bite force is statistically associated with an increase in brain size.
Alternatively, the reduction in bite force may be owing to neoteny, whereby the jaw muscles arrest in development along the hominin lineage. However, neoteny in humans has also been interpreted as being associated with brain enlargement so the precise cause of the reduction in jaw adductor musculature and thus bite force is up for debate.
5. Conclusion
Taken together, our results reveal that the evolution of bite force in amniotes occurred as bursts of accelerated changes across multiple clades, and as the product of repeated and nested pulses of progressively higher rates of change, representing instances of exceptional functional adaptations. Using a phylogenetic evolutionary framework, on a dataset representing the largest taxonomic sample to date, enables us to statistically detect instances of adaptations in biomechanical metrics, not only along specific branches, but also through time, paving the way to better understand how specific ecological niches (feeding ecologies) are occupied. In order to determine whether species' bite force underwent instances of exceptional rates of adaptive changes, it is necessary to demonstrate an exceptionally high rate of evolution associated with that species in a phylogenetic context after accounting for size and the expected evolutionary change associated with divergence time.
Data accessibility
Data are available on the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.q12c06f [72].
Authors' contributions
M.S., M.R. and C.V. conceived of the study, designed the study, coordinated the study and drafted the manuscript. M.S. collected and analysed the data. All authors gave final approval for publication.
Competing interests
We have no competing interests.
Funding
M.S. and C.V. were supported by the Leverhulme Trust (RPG-2013-185 and RPG-2017-071) and M.S. and M.R. were supported by the BBSRC (BB/H007954/1).
Acknowledgements
We thank Joanna Baker, Ciara O'Donovan, Andrew Meade, Jorge Avaria Llautureo and Henry Ferguson-Gow for invaluable discussion, which helped improve the research and manuscript. M.S. thanks the curators and collection managers at the various institutions visited as part of the data collection process, in particular, Rhian Rowson (Bristol City Museum and Arts Gallery), Paolo Viscardi (National Museum of Ireland—Natural History; formerly Horniman Museum), Mark Carnall (Oxford University Museum; formerly Grant Museum, UCL), Lars Werdelin (Swedish Museum of Natural History), Andrew Kitchener (National Museums Scotland—Natural Sciences), Judy Galkin (American Museum of Natural History), Roberto Portela Miguez (Natural History Museum, London) and Milly Farrell (Oxford Brookes University; formerly Hunterian Museum, Royal College of Surgeons). We also thank Larry Witmer (Ohio University) for permission to use his gallery of various theropod skull images. Finally, we thank the reviewers for their constructive comments on this manuscript. Silhouettes in figures 1 and 2 are attributed to Andrew Farke (Stegosaurus, CC-BY 3.0), Michael Keesey (Homo sapiens, CC-BY 3.0; Gorilla, CC0 1.0), Steven Traver (Sphenodon punctatus, CC0 1.0), Sarah Werning (Panthera tigris, CC-BY 3.0) and Yan Wong (Chiropteran, CC0 1.0). Anthropoid is uncredited (CC0 1.0). Silhouettes of finches were drawn using photographs by Vince Smith (Geospiza fuliginosa, CC-BY 2.0) and Brian Gratwicke (Platyspiza crassirostris, CC-BY 2.0). Silhouette of Deinonychus antirrhopus is the author's own work (MS, CC-BY 3.0).
References
- 1.
Grant PR . 1981 The feeding of Darwin's finches on Tribulus cistoides (L.) seeds. Anim. Behav. 29, 785-793. (doi:10.1016/S0003-3472(81)80012-7) Crossref, Web of Science, Google Scholar - 2.
Schluter D, Price TD, Grant PR . 1985 Ecological character displacement in Darwin's finches. Science 227, 1056-1059. (doi:10.1126/science.227.4690.1056) Crossref, PubMed, Web of Science, Google Scholar - 3.
Wroe S, McHenry C, Thomason J . 2005 Bite club: comparative bite force in big biting mammals and the prediction of predatory behaviour in fossil taxa. Proc. R. Soc. B 272, 619-625. (doi:10.1098/rspb.2004.2986) Link, Web of Science, Google Scholar - 4.
Christiansen P, Wroe S . 2007 Bite forces and evolutionary adaptations to feeding ecology in carnivores. Ecology 88, 347-358. (doi:10.1890/0012-9658(2007)88[347:BFAEAT]2.0.CO;2) Crossref, PubMed, Web of Science, Google Scholar - 5.
Wilson GP, Ekdale EG, Hoganson JW, Calede JJ, Vander Linden A . 2016 A large carnivorous mammal from the Late Cretaceous and the North American origin of marsupials. Nat. Commun. 7, 13734. (doi:10.1038/ncomms13734) Crossref, PubMed, Web of Science, Google Scholar - 6.
Schmidt-Nielsen K . 1984 Scaling: why is animal size so important? Cambridge, UK: Cambridge University Press. Crossref, Google Scholar - 7.
Felsenstein J . 1985 Phylogenies and the comparative method. Am. Nat. 125, 1-15. (doi:10.1086/284325) Crossref, Web of Science, Google Scholar - 8.
Harvey PH, Pagel MD . 1991 The comparative method in evolutionary biology. Oxford, UK: Oxford University Press. Google Scholar - 9.
Baker J, Meade A, Pagel M, Venditti C . 2016 Positive phenotypic selection inferred from phylogenies. Biol. J. Linnean Soc. 118, 95-115. (doi:10.1111/bij.12649) Crossref, Web of Science, Google Scholar - 10.
Simpson GG . 1944 Tempo and mode in evolution, p. 237. New York, NY: Columbia University Press. Google Scholar - 11.
Yang Z . 2002 Inference of selection from multiple species alignments. Curr. Opin Genet. Dev. 12, 688-694. (doi:10.1016/S0959-437X(02)00348-9) Crossref, PubMed, Web of Science, Google Scholar - 12.
Yang Z . 2006 Computational molecular evolution, p. 376. Oxford, UK: Oxford University Press. Google Scholar - 13.
Ghalambor CK, Walker JA, Reznick DN . 2003 Multi-trait selection, adaptation, and constraints on the evolution of burst swimming performance. Integr. Comp. Biol. 43, 431–438. (doi:10.1093/icb/43.3.431) Crossref, PubMed, Web of Science, Google Scholar - 14.
Walker J 2007 A general model of functional constraints on phenotypic evolution. Am. Nat. 170, 681-689. (doi:10.1086/521957) Crossref, PubMed, Web of Science, Google Scholar - 15.
Ackerly D . 2009 Conservatism and diversification of plant functional traits: evolutionary rates versus phylogenetic signal. Proc. Natl Acad. Sci. USA 106, 19699-19706. (doi:10.1073/pnas.0901635106) Crossref, PubMed, Web of Science, Google Scholar - 16.
Collar DC, O'Meara BC, Wainwright PC, Near TJ . 2009 Piscivory limits diversification of feeding morphology in centrarchid fishes. Evolution 63, 1557-1573. (doi:10.1111/j.1558-5646.2009.00626.x) Crossref, PubMed, Web of Science, Google Scholar - 17.
Martin CH, Wainwright PC . 2011 Trophic novelty is linked to exceptional rates of morphological diversification in two adaptive radiations of cyprinodon pupfish. Evolution 65, 2197-2212. (doi:10.1111/j.1558-5646.2011.01294.x) Crossref, PubMed, Web of Science, Google Scholar - 18.
Holzman R, Collar DC, Price SA, Hulsey CD, Thomson RC, Wainwright PC . 2012 Biomechanical trade-offs bias rates of evolution in the feeding apparatus of fishes. Proc. R. Soc. B 279, 1287-1292. (doi:10.1098/rspb.2011.1838) Link, Web of Science, Google Scholar - 19.
Wainwright PC, Smith WL, Price SA, Tang KL, Sparks JS, Ferry LA, Kuhn KL, Eytan RI, Near TJ . 2012 The evolution of pharyngognathy: a phylogenetic and functional appraisal of the pharyngeal jaw key innovation in labroid fishes and beyond. Syst. Biol. 61, 1001-1027. (doi:10.1093/sysbio/sys060) Crossref, PubMed, Web of Science, Google Scholar - 20.
Anderson PSL, Claverie T, Patek SN . 2014 Levers and linkages: mechanical trade-offs in a power-amplified system. Evolution 68, 1919-1933. (doi:10.1111/evo.12407) Crossref, PubMed, Web of Science, Google Scholar - 21.
Collar DC, Wainwright PC, Alfaro ME, Revell LJ, Mehta RS . 2014 Biting disrupts integration to spur skull evolution in eels. 5, 5505. (doi:10.1038/ncomms6505) Google Scholar - 22.
Anderson PSL, Patek SN . 2015 Mechanical sensitivity reveals evolutionary dynamics of mechanical systems. Proc. R. Soc. B 282, 20143088. (doi:10.1098/rspb.2014.3088) Link, Web of Science, Google Scholar - 23.
Muñoz MM, Anderson PS, Patek SN . 2017 Mechanical sensitivity and the dynamics of evolutionary rate shifts in biomechanical systems. Proc. R. Soc. B 284, 20162325. (doi:10.1098/rspb.2016.2325) Link, Web of Science, Google Scholar - 24.
Stayton CT . 2011 Biomechanics on the half shell: functional performance influences patterns of morphological variation in the emydid turtle carapace. Zoology 114, 213-223. (doi:10.1016/j.zool.2011.03.002) Crossref, PubMed, Web of Science, Google Scholar - 25.
Organ C, Nunn CL, Machanda Z, Wrangham RW . 2011 Phylogenetic rate shifts in feeding time during the evolution of Homo. Proc. Natl Acad. Sci. USA 108, 14555-14559. (doi:10.1073/pnas.1107806108) Crossref, PubMed, Web of Science, Google Scholar - 26.
Piras P, Sansalone G, Teresi L, Kotsakis T, Colangelo P, Loy A . 2012 Testing convergent and parallel adaptations in talpids humeral mechanical performance by means of geometric morphometrics and finite element analysis. J. Morphol. 273, 696-711. (doi:10.1002/jmor.20015) Crossref, PubMed, Web of Science, Google Scholar - 27.
Hagey TJ, Uyeda JC, Crandell KE, Cheney JA, Autumn K, Harmon LJ . 2017 Tempo and mode of performance evolution across multiple independent origins of adhesive toe pads in lizards. Evolution 71, 2344-2358. (doi:10.1111/evo.13318) Crossref, PubMed, Web of Science, Google Scholar - 28.
van der Meij MA, Bout RG. 2006 Seed husking time and maximal bite force in finches. J. Exp. Biol. 209, 3329-3335. (doi:10.1242/jeb.02379) Crossref, PubMed, Web of Science, Google Scholar - 29.
Verwaijen D, Van Damme R, Herrel A. 2002 Relationships between head size, bite force, prey handling efficiency and diet in two sympatric lacertid lizards. Funct. Ecol. 16, 842-850. (doi:10.1046/j.1365-2435.2002.00696.x) Crossref, Web of Science, Google Scholar - 30.
Vanhooydonck B, Herrel AY, Van Damme R, Irschick DJ. 2005 Does dewlap size predict male bite performance in Jamaican Anolis lizards? Funct. Ecol. 19, 38-42. (doi:10.1111/j.0269-8463.2005.00940.x) Crossref, Web of Science, Google Scholar - 31.
Lappin A, Kristopher BY, Husak JF, Macedonia JM . 2006 Gaping displays reveal and amplify a mechanically based index of weapon performance. Am. Nat. 168, 100-113. (doi:10.1086/505161) Crossref, PubMed, Web of Science, Google Scholar - 32.
Hedges SB, Marin J, Suleski M, Paymer M, Kumar S . 2015 Tree of life reveals clock-like speciation and diversification. Mol. Biol. Evol. 32, 835-845. (doi:10.1093/molbev/msv037) Crossref, PubMed, Web of Science, Google Scholar - 33.
Pagel M, Meade A, Barker D . 2004 Bayesian estimation of ancestral character states on phylogenies. Syst. Biol. 53, 673-684. (doi:10.1080/10635150490522232) Crossref, PubMed, Web of Science, Google Scholar - 34.
van der Meij MA, Bout RG. 2004 Scaling of jaw muscle size and maximal bite force in finches. J. Exp. Biol. 207, 2745-2753. (doi:10.1242/jeb.01091) Crossref, PubMed, Web of Science, Google Scholar - 35.
Herrel A, Podos J, Huber SK, Hendry AP . 2005 Evolution of bite force in Darwin's finches: a key role for head width. J. Evol. Biol. 18, 669-675. (doi:10.1111/j.1420-9101.2004.00857.x) Crossref, PubMed, Web of Science, Google Scholar - 36.
Senawi J, Schmieder D, Siemers B, Kingston T . 2015 Beyond size: morphological predictors of bite force in a diverse insectivorous bat assemblage from Malaysia. Funct. Ecol. 29, 1411-1420. (doi:10.1111/1365-2435.12447) Crossref, Web of Science, Google Scholar - 37.
Brusatte SL, Lloyd GT, Wang SC, Norell MA . 2014 Gradual assembly of avian body plan culminated in rapid rates of evolution across the dinosaur-bird transition. Curr. Biol. 24, 2386-2392. (doi:10.1016/j.cub.2014.08.034) Crossref, PubMed, Web of Science, Google Scholar - 38.
Lee MSY, Cau A, Naish D, Dyke GJ . 2014 Sustained miniaturization and anatomical innovation in the dinosaurian ancestors of birds. Science 345, 562-566. (doi:10.1126/science.1252243) Crossref, PubMed, Web of Science, Google Scholar - 39.
Puttick MN, Thomas GH, Benton MJ . 2014 High rates of evolution preceded the origin of birds. Evolution 68, 1497-1510. (doi:10.1111/evo.12363) Crossref, PubMed, Web of Science, Google Scholar - 40.
Brusatte SL, O'Connor JK, Jarvis ED . 2015 The origin and diversification of birds. Curr. Biol. 25, R888-R898. (doi:10.1016/j.cub.2015.08.003) Crossref, PubMed, Web of Science, Google Scholar - 41.
Brusatte SL, Nesbitt SJ, Irmis RB, Butler RJ, Benton MJ, Norell MA . 2010 The origin and early radiation of dinosaurs. Earth Sci. Rev. 101, 68-100. (doi:10.1016/j.earscirev.2010.04.001) Crossref, Web of Science, Google Scholar - 42.
Raftery AE . 1996 Hypothesis testing and model selection. In Markov chain Monte Carlo in practice (edsGilks WR, Richardson S, Spiegelhalter DJ ), pp. 163-187. London, UK: Chapman & Hall. Google Scholar - 43.
R Core Development Team. 2011 R: a language and environment for statistical computing, (2.10.1 edn). Vienna, Austria: R Foundation for Statistical Computing. Google Scholar
- 44.
Will M, Pablos A, Stock JT . 2017 Long-term patterns of body mass and stature evolution within the hominin lineage. R. Soc. Open Sci. 4, 171339. (doi:10.1098/rsos.171339) Link, Web of Science, Google Scholar - 45.
Gignac PM, Erickson GM . 2017 The biomechanics behind extreme osteophagy in Tyrannosaurus rex. Sci. Rep. 7, 2012. (doi:10.1038/s41598-017-02161-w) Crossref, PubMed, Web of Science, Google Scholar - 46.
Henderson DM, Snively E . 2004 Tyrannosaurus en pointe: allometry minimized rotational inertia of large carnivorous dinosaurs. Proc. R. Soc. Lond. B 271, S57-S60. (doi:10.1098/rsbl.2003.0097) Link, Web of Science, Google Scholar - 47.
Hutchinson JR, Bates KT, Molnar J, Allen V, Makovicky PJ . 2011 A computational analysis of limb and body dimensions in Tyrannosaurus rex with implications for locomotion, ontogeny, and growth. PLoS One 6, e26037. (doi:10.1371/journal.pone.0026037) Crossref, PubMed, Web of Science, Google Scholar - 48.
Meers MB . 2002 Maximum bite force and prey size of Tyrannosaurus rex and their relationships to the inference of feeding behaviour. Hist. Biol. 16, 1-12. (doi:10.1080/0891296021000050755) Crossref, Google Scholar - 49.
Binder WJ, Van Valkenburgh B. 2000 Development of bite strength and feeding behaviour in juvenile spotted hyenas (Crocuta crocuta). J. Zool. 252, 273-283. (doi:10.1111/j.1469-7998.2000.tb00622.x) Crossref, Web of Science, Google Scholar - 50.
Radinsky LB . 1981 Evolution of skull shape in carnivores.1. Representative modern carnivores. Biol. J. Linnean Soc. 15, 369-388. (doi:10.1111/j.1095-8312.1981.tb00770.x) Crossref, Web of Science, Google Scholar - 51.
Radinsky LB . 1981 Evolution of skull shape in carnivores. 2. Additional modern carnivores. Biol. J. Linnean Soc. 16, 337-355. (doi:10.1111/j.1095-8312.1981.tb01657.x) Crossref, Web of Science, Google Scholar - 52.
Wroe S, Troy M, Seebacher F, Ben K, Gillespie A, Mathew C, Steve S . 2003 An alternative method for predicting body mass: the case of the pleistocene marsupial lion. Paleobiology 29, 403-411. (doi:10.1666/0094-8373(2003)029<0403:AAMFPB>2.0.CO;2) Crossref, Web of Science, Google Scholar - 53.
Osborn HF . 1900 The geological and faunal relations of Europe and America during the tertiary period and the theory of the successive invasions of an African fauna. Science 11, 561-574. (doi:10.1126/science.11.276.561) Crossref, PubMed, Google Scholar - 54.
Schluter D . 2000 The ecology of adaptive radiation. Oxford, UK: Oxford University Press. Crossref, Google Scholar - 55.
Venditti C, Meade A, Pagel M . 2011 Multiple routes to mammalian diversity. Nature 479, 393-396. (doi:10.1038/nature10516) Crossref, PubMed, Web of Science, Google Scholar - 56.
Grant BR, Grant PR . 1993 Evolution of Darwin finches caused by a rare climatic event. Proc. R. Soc. Lond. B 251, 111-117. (doi:10.1098/rspb.1993.0016) Link, Web of Science, Google Scholar - 57.
Grant PR, Grant BR . 2006 Evolution of character displacement in Darwin's finches. Science 313, 224-226. (doi:10.1126/science.1128374) Crossref, PubMed, Web of Science, Google Scholar - 58.
Lamichhaney S 2015 Evolution of Darwin's finches and their beaks revealed by genome sequencing. Nature 518, 371-375. (doi:10.1038/nature14181) Crossref, PubMed, Web of Science, Google Scholar - 59.
Schluter D, Grant PR . 1984 Determinants of morphological patterns in communities of Darwin's finches. Am. Nat. 123, 175-196. (doi:10.1086/284196) Crossref, Web of Science, Google Scholar - 60.
Price TD, Grant PR, Gibbs HL, Boag PT . 1984 Recurrent patterns of natural selection in a population of Darwin's finches. Nature 309, 787-789. (doi:10.1038/309787a0) Crossref, PubMed, Web of Science, Google Scholar - 61.
Pampush JD, Scott JE, Robinson CA, Delezene LK . 2018 Oblique human symphyseal angle is associated with an evolutionary rate-shift early in the hominin clade. J. Hum. Evol. 123, 84-95. (doi:10.1016/j.jhevol.2018.06.006) Crossref, PubMed, Web of Science, Google Scholar - 62.
Hublin J-J, Neubauer S, Gunz P . 2015 Brain ontogeny and life history in Pleistocene hominins. Phil. Trans. R. Soc. B 370, 20140062. (doi:10.1098/rstb.2014.0062) Link, Web of Science, Google Scholar - 63.
Stedman HH 2004 Myosin gene mutation correlates with anatomical changes in the human lineage. Nature 428, 415-418. (doi:10.1038/nature02358) Crossref, PubMed, Web of Science, Google Scholar - 64.
Bruner E . 2004 Geometric morphometrics and paleoneurology: brain shape evolution in the genus Homo. J. Hum. Evol. 47, 279-303. (doi:10.1016/j.jhevol.2004.03.009) Crossref, PubMed, Web of Science, Google Scholar - 65.
McCollum MA, Sherwood CC, Vinyard CJ, Lovejoy CO, Schachat F . 2006 Of muscle-bound crania and human brain evolution: the story behind the MYH16 headlines. J. Hum. Evol. 50, 232-236. (doi:10.1016/j.jhevol.2005.10.003) Crossref, PubMed, Web of Science, Google Scholar - 66.
Perry GH, Verrelli BC, Stone AC . 2005 Comparative analyses reveal a complex history of molecular evolution for human MYH16. Mol. Biol. Evol. 22, 379-382. (doi:10.1093/molbev/msi004) Crossref, PubMed, Web of Science, Google Scholar - 67.
Raadsheer MC, van Eijden TM, van Ginkel FC, Prahl-Andersen B. 1999 Contribution of jaw muscle size and craniofacial morphology to human bite force magnitude. J. Dent. Res. 78, 31-42. (doi:10.1177/00220345990780010301) Crossref, PubMed, Web of Science, Google Scholar - 68.
Orme D, Freckleton R, Thomas G, Petzoldt T, Fritz S, Isaac N, Pearse W . 2013 caper: comparative analyses of phylogenetics and evolution in R. R package version 0.5.2. (R package version 0.5.2 edn). Google Scholar - 69.
Rightmire GP . 2004 Brain size and encephalization in early to mid-Pleistocene Homo. Am. J. Phys. Anthropol. 124, 109-123. (doi:10.1002/ajpa.10346) Crossref, PubMed, Web of Science, Google Scholar - 70.
DeSilva JM, Lesnik JJ . 2008 Brain size at birth throughout human evolution: a new method for estimating neonatal brain size in hominins. J. Hum. Evol. 55, 1064-1074. (doi:10.1016/j.jhevol.2008.07.008) Crossref, PubMed, Web of Science, Google Scholar - 71.
Isler K, Christopher Kirk E, Miller JMA, Albrecht GA, Gelvin BR, Martin RD . 2008 Endocranial volumes of primate species: scaling analyses using a comprehensive and reliable data set. J. Hum. Evol. 55, 967-978. (doi:10.1016/j.jhevol.2008.08.004) Crossref, PubMed, Web of Science, Google Scholar - 72.
Sakamoto M, Ruta M, Venditti C . 2019Data from: Extreme and rapid bursts of functional adaptations shape bite force in amniotes . Dryad Digital Repository. (http://dx.doi.org/10.5061/dryad.q12c06f) Google Scholar