Testing hypotheses of convergence with multivariate data: morphological and functional convergence among herbivorous lizards

Abstract Despite its importance to evolutionary theory, convergence remains an understudied phenomenon and is usually investigated using qualitative data. This paper advances a new, multidimensional view of convergence. Three patterns indicative of convergence are discussed, and techniques to discover and test convergent patterns in a quantitative framework are developed. These concepts and methods are applied to a dataset of digitized coordinates on 1554 lizard skulls and 1292 lower jaws to test hypotheses of convergence among herbivorous lizards. Encompassing seven independent acquisitions of herbivory, this lizard sample provides an ideal natural experiment for exploring ideas of convergence among different systems (here, morphological and functional). Three related questions are addressed: (1) Do herbivorous lizards show evidence of convergence in skull and lower jaw morphology? (2) What, if any, is the morphospace pattern associated with this convergence? (3) Is it possible to predict the direction of convergence using functional models? Relative warp analysis and permutation tests reveal that the skulls and lower jaws of herbivorous lizards do show evidence of convergence. Herbivore skulls deviate from their carnivorous or omnivorous sister groups toward the same area of morphospace. Without a phylogenetic perspective, this pattern would not be recognizable. Lower jaws of herbivores are not convergent in morphology but are convergent in function: herbivores deviate away from their carnivorous sister groups toward higher values of mechanical advantage. These results illustrate the desirability of quantitative methods, informed by phylogenetic information, in the study of convergence.     

Molecular phylogenetics of squamata: the position of snakes, amphisbaenians, and dibamids, and the root of the squamate tree

Squamate reptiles (snakes, lizards, and amphisbaenians) serve as model systems for evolutionary studies of a variety of morphological and behavioral traits, and phylogeny is crucial to many generalizations derived from such studies. Specifically, the traditional dichotomy between Iguania (anoles, iguanas, chameleons, etc.) and Scleroglossa (skinks, geckos, snakes, etc.) has been correlated with major evolutionary shifts within Squamata. We present a molecular phylogenetic study of 69 squamate species using approximately 4600 (2876 parsimony-informative) base pairs (bp) of DNA sequence data from the nuclear genes RAG-1(approximately 2750 bp) and c-mos(approximately 360 bp) and the mitochondrial ND2 region (approximately 1500 bp), sampling all major clades and most major subclades. Under our hypothesis, species previously placed in Iguania, Anguimorpha, and almost all recognized squamate families form strongly supported monophyletic groups. However, species previously placed in Scleroglossa, Varanoidea, and several other higher taxa do not form monophyletic groups. Iguania, the traditional sister group of Scleroglossa, is actually highly nested within Scleroglossa. This unconventional rooting does not seem to be due to long-branch attraction, base composition biases among taxa, or convergence caused by similar selective forces acting on nonsister taxa. Studies of functional tongue morphology and feeding mode have contrasted the similar states found in Sphenodon(the nearest outgroup to squamates) and Iguania with those of Scleroglossa, but our findings suggest that similar states in Sphenodonand Iguania result from homoplasy. Snakes, amphisbaenians, and dibamid lizards, limbless forms whose phylogenetic positions historically have been impossible to place with confidence, are not grouped together and appear to have evolved this condition independently. Amphisbaenians are the sister group of lacertids, and dibamid lizards diverged early in squamate evolutionary history. Snakes are grouped with iguanians, lacertiforms, and anguimorphs, but are not nested within anguimorphs.     

Lipophilic compounds in femoral secretions of male collared lizards,Crotaphytus bicinctores (Iguania,Crotaphytidae)

In many lizards, chemical compounds from the femoral gland secretions are used in intraspecific communication, but most studies describing these chemicals are for lizard species included in the Scleroglossa clade, whereas lizards within the Iguanian clade have been much less studied, probably because these lizards were considered to rely more on visual cues. However, many iguanian lizards have abundant femoral secretions and are able of chemosensory conspecific recognition, which might be based on compounds secreted by femoral glands. By using GC–MS analyses, we found 58 lipophilic compounds in femoral gland secretions of male Great Basin collared lizard, Crotaphytus bicinctores (Iguania, Crotaphytidae). Main compounds were steroids (mainly two triunsaturated steroids and cholesterol), carboxylic acids (mainly hexadecanoic acid), waxy esters of long chain fatty acids, alcohols (mainly hexadecanol), aldehydes and other minor compounds. We compared these compounds with those found in other lizard species and discussed the potential signaling function of some compounds and how the xeric habitat of this lizard could have conditioned the composition of secretions.     

Shape at the cross‐roads: homoplasy and history in the evolution of the carnivoran skull towards herbivory

Patterns of skull shape in Carnivora provide examples of parallel and convergent evolution for similar ecomorphological adaptations. However, although most researchers report on skull homoplasies among hypercarnivorous taxa, evolutionary trends towards herbivory remain largely unexplored. In this study, we analyse the skull of the living herbivorous carnivorans to evaluate the importance of natural selection and phylogenetic legacy in shaping the skulls of these peculiar species. We quantitatively estimated shape variability using geometric morphometrics. A principal components analysis of skull shape incorporating all families of arctoid carnivorans recognized several common adaptations towards herbivory. Ancestral state reconstructions of skull shape and the reconstructed phylogenetic history of morphospace occupation more explicitly reveal the true patterns of homoplasy among the herbivorous carnivorans. Our results indicate that both historical constraints and adaptation have interplayed in the evolution towards herbivory of the carnivoran skull, which has resulted in repeated patterns of biomechanical homoplasy.     

Absence of Prey Chemical Discrimination by Tongue-flicking in an Ambush-foraging Lizard Having Actively Foraging Ancestors

Lingually mediated detection of prey chemicals is widespread in one major clade of lizards, Scleroglossa, but rare in the other, Iguania. It is absent in all ambush-foraging families tested and present in all actively foraging families. In Iguania, prey chemical discrimination is known only in the herbivorous Iguanidae; in Scleroglossa, it was heretofore known to be absent only in ambush-foraging gekkonids. Because ambush foraging precludes lingual sampling of a wide area and tongue-flicking would disrupt the crypticity ambushers maintain by immobility, we predicted that prey chemical discrimination would be absent in scleroglossans that have secondarily adopted ambush foraging. The Cape girdled lizard, Cordylus cordylus, is member of Cordylidae, a family of ambush foragers considered derived from active foragers in the Autarchoglossa, a group of scleroglossan families having highly developed lingual chemosensory behaviours. As predicted, this species did not discriminate surface chemicals of three prey species from control substances in a series of standardized experiments in which prey chemicals were presented on cotton-tipped applicators. Thus, even in taxa having highly developed prey chemical discrimination, adoption of ambush foraging may induce loss of prey chemical discrimination, providing further and stronger evidence that prey chemical discrimination is adaptively adjusted to foraging mode.     

Determinants of assemblage structure in Neotropical dry forest lizards

We investigated the structure of a lizard assemblage from Seasonally Dry Tropical Forest enclaves in the Brazilian Cerrado biome, by testing the roles of ecological and historical components. We analysed data from 469 individuals, belonging to 18 lizard species, sampled by a combination of pitfall, funnel and glue traps, as well as by haphazard sampling. Null model analyses and Canonical Phylogenetic Ordination analysis, coupled with Monte Carlo simulations, revealed lack of both ecological and phylogenetic structure in microhabitat use. Conversely, these analyses revealed a mean overlap in diet composition significantly smaller than expected by chance and significant historical structure. Structure in diet composition was due to phylogenetic effects corresponding to the most basal divergence of the squamate phylogeny (Iguania/Scleroglossa) and the clades Teiidae and Gymnophthalmidae. Among lizards, evolutionary constraints on microhabitat use appear less than on prey use, suggesting that the availability of historically preferred prey types moderates microhabitat selection. The lack of structure in microhabitat use suggests absence of competitive interactions on the spatial component. On the other hand, food preferences have a deep historical basis and do not reflect current competitive interactions.     

How phylogeny and foraging ecology drive the level of chemosensory exploration in lizards and snakes

The chemical senses are crucial for squamates (lizards and snakes). The extent to which squamates utilize their chemosensory system, however, varies greatly among taxa and species’ foraging strategies, and played an influential role in squamate evolution. In lizards, ‘Scleroglossa’ evolved a state where species use chemical cues to search for food (active foragers), whereas ‘Iguania’ retained the use of vision to hunt prey (ambush foragers). However, such strict dichotomy is flawed as shifts in foraging modes have occurred in all clades. Here, we attempted to disentangle effects of foraging ecology from phylogenetic trait conservatism as leading cause of the disparity in chemosensory investment among squamates. To do so, we used species’ tongue‐flick rate (TFR) in the absence of ecological relevant chemical stimuli as a proxy for its fundamental level of chemosensory investigation, that is baseline TFR. Based on literature data of nearly 100 species and using phylogenetic comparative methods, we tested whether and how foraging mode and diet affect baseline TFR. Our results show that baseline TFR is higher in active than ambush foragers. Although baseline TFRs appear phylogenetically stable in some lizard taxa, that is a consequence of concordant stability of foraging mode: when foraging mode shifts within taxa, so does baseline TFR. Also, baseline TFR is a good predictor of prey chemical discriminatory ability, as we established a strong positive relationship between baseline TFR and TFR in response to prey. Baseline TFR is unrelated to diet. Essentially, foraging mode, not phylogenetic relatedness, drives convergent evolution of similar levels of squamate chemosensory investigation.     

Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.     

Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach

Understanding how and why certain clades diversify greatly in morphology whereas others do not remains a major theme in evolutionary biology. Projecting families of phylogenies into multivariate morphospaces can distinguish two scenarios potentially leading to unequal morphological diversification: unequal magnitude of change per phylogenetic branch, and unequal efficiency in morphological innovation. This approach is demonstrated using a case study of skulls in sister clades within the South American fish superfamily Anostomoidea. Unequal morphological diversification in this system resulted not from the morphologically diverse clade changing more on each phylogenetic branch, but from that clade distributing an equal amount of change more widely through morphospace and innovating continually. Although substantial morphological evolution occurred throughout the less diverse clade's history, most of that clade's expansion in morphospace occurred in the most basal branches, and more derived portions of that radiation oscillated within previously explored limits. Because simulations revealed that there is a maximum 2.7% probability of generating two clades that differ so greatly in the density of lineages within morphospace under a null Brownian model, the observed difference in pattern likely reflects a difference in the underlying evolutionary process. Clade-specific factors that may have promoted or arrested morphological diversification are discussed.     

The skull evolution of oviraptorosaurian dinosaurs: the role of niche partitioning in diversification

Oviraptorosaurs are bird‐like theropod dinosaurs that thrived in the final pre‐extinction ecosystems during the latest Cretaceous, and the beaked, toothless skulls of derived species are regarded as some of the most peculiar among dinosaurs. Their aberrant morphologies are hypothesized to have been caused by rapid evolution triggered by an ecological/biological driver, but little is known about how their skull shapes and functional abilities diversified. Here, we use quantitative techniques to study oviraptorosaur skull form and mandibular function. We demonstrate that the snout is particularly variable, that mandibular form and upper/lower beak form are significantly correlated with phylogeny, and that there is a strong and significant correlation between mandibular function and mandible/lower beak shape, suggesting a form–function association. The form–function relationship and phylogenetic signals, along with a moderate allometric signal in lower beak form, indicate that similar mechanisms governed beak shape in oviraptorosaurs and extant birds. The two derived oviraptorosaur clades, oviraptorids and caenagnathids, are significantly separated in morphospace and functional space, indicating that they partitioned niches. Oviraptorids coexisting in the same ecosystem are also widely spread in morphological and functional space, suggesting that they finely partitioned feeding niches, whereas caenagnathids exhibit extreme disparity in beak size. The diversity of skull form and function was likely key to the diversification and evolutionary success of oviraptorosaurs in the latest Cretaceous.     

GEOMETRIC MORPHOMETRICS OF THE SKULL ROOF OF STEREOSPONDYLS (AMPHIBIA: TEMNOSPONDYLI)

Abstract:  Geometric morphometric analysis using relative warps is applied to the skull roof of 62 species of stereospondyls and their closest outgroups (i.e. basal archegosauriforms) from among temnospondyl amphibians. Twenty-one landmarks and five taxonomic groups are used for comparisons. Their skull evolution is quantified in a morphospace defined by two relative warps axes. The majority of groups show poor concordance between morphological and phylogenetic distances. The only exception is represented by Yates and Warren's study of stereospondyl relationships, in which concordance is high. Only basal archegosauriforms and rhinesuchids show significant overlap in morphospace, although this might be due to low sample sizes. Regression of estimated mean disparity against taxon sample size shows that species within both the trematosauroid and the rhytidostean groups are more widely dispersed in morphospace than species belonging to any of the remaining stereospondyl groups. Stereospondyl skull evolution was characterized by divergence between major clades and convergence within those clades. Changes in patterns of morphospace occupation through time agree with the hypothesis of an 'explosive' radiation in the early Early Triassic, after the extinction of basal archegosauriforms at the end of the Permian.     

Soft anatomy, diffuse homoplasy, and the relationships of lizards and snakes

Most previous phylogenetic analyses of squamates (‘lizards’ and snakes) employing large character sets have focused on osteology. Soft anatomical traits bearing on this problem have usually been considered in small subsets. Here, a comprehensive phylogenetic analysis of squamate soft anatomy is attempted. 126 informative characters are assessed for 23 squamate lineages, representing snakes, amphisbaenians, dibamids, and all the traditionally recognized ‘families’ of lizards. The traditionally recognized groupings Iguania, Scleroglossa, Gekkota, Scincomorpha, Anguimorpha and Varanoidea are corroborated in this analysis. More controversial taxa are resolved as follows. Xantusiids, amphisbaenians and dibamids cluster with gekkotans, and snakes are strongly allied with anguimorphs in general, and varanids in particular. Nearly all these clades are congruent with those found in a recent comprehensive osteological analysis; the strong support for snake‐varanid relationships found in both studies is particularly notable. This congruence is surprising given that previous studies of soft anatomy tended to give differing and often heterodox results. These previous results can be attributed to overrepresentation of misleading characters in small isolated data sets. Such misleading signals are minimized when data sets are combined. For instance, the snake‐varanid clade is contradicted by many characters, and analyses of particular organ systems therefore give differing results. However, characters that are incongruent with the snake‐varanid clade also disagree with each other (diffuse homoplasy), rather than forming coherent support for some particular alternative clade (concerted homoplasy). In a combined analysis these incongruent but diffuse characters cancel each other out to leave a very strong (and orthodox) phylogenetic signal. These results underscore the view that the raw amount of homoplasy — as revealed by consistency and retention indices — is not the only determinant of phylogenetic signal; the distribution of that homoplasy is also important. Thus, questioning a phylogenetic hypothesis (e.g. the snake‐varanid clade) by identifying numerous conflicting characters is insufficient — the structure of the conflicting characters should be assessed in a rigorous phylogenetic analysis.     

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.     

The evolution of cranial form and function in theropod dinosaurs: insights from geometric morphometrics

Theropod dinosaurs, an iconic clade of fossil species including Tyrannosaurus and Velociraptor, developed a great diversity of body size, skull form and feeding habits over their 160+ million year evolutionary history. Here, we utilize geometric morphometrics to study broad patterns in theropod skull shape variation and compare the distribution of taxa in cranial morphospace (form) to both phylogeny and quantitative metrics of biting behaviour (function). We find that theropod skulls primarily differ in relative anteroposterior length and snout depth and to a lesser extent in orbit size and depth of the cheek region, and oviraptorosaurs deviate most strongly from the "typical" and ancestral theropod morphologies. Noncarnivorous taxa generally fall out in distinct regions of morphospace and exhibit greater overall disparity than carnivorous taxa, whereas large-bodied carnivores independently converge on the same region of morphospace. The distribution of taxa in morphospace is strongly correlated with phylogeny but only weakly correlated with functional biting behaviour. These results imply that phylogeny, not biting function, was the major determinant of theropod skull shape.     

Patterns of caudal-autotomy evolution in lizards

Using comparative techniques to account for phylogenetic effects, I examined patterns of evolution of caudal autotomy and foraging in 39 lizard species to test the hypothesis that caudal autotomy has co-evolved with morphology, locomotor performance, and foraging behaviour. There were significant positive associations between evolution of the point on the tail (distance from cloaca) at which tail loss occurs (an indirect measure of caudal autotomy) and evolution of each of the following: tail length, caudifemoralis longus (CFL) muscle length, and jump distance. The correlation with the evolution of sprint speed approached significance. These relationships primarily were due to the influence of tail-length evolution on autotomy-point evolution. With the effect of tail-length evolution removed, autotomy-point evolution was negatively correlated with the evolution of tail-loss frequency. The CFL restricts tail loss to portions of the tail posterior to the most distal point of its insertion in the tail. In addition, with the effect of tail-length evolution removed, CFL length co-evolved with sprint speed. These results indicate that tail morphology has co-evolved with caudal autotomy such that the evolution of the CFL has reduced caudal autotomy in certain groups of lizards.
Ambush foraging, the ability to lose the tail, intermediate CFL length, and low locomotor performance (i.e. slow sprint speed and short jump distance) are hypothesized to be the ancestral conditions in lizards using outgroup rooting. The diversification of lizard taxa has resulted in some lineages moving away from ancestral character states (i.e. family Teiidae, superfamily Varanoidea), while others are very similar or identical to their ancestors (i.e. superfamily Iguania).     

Form and function of damselfish skulls: rapid and repeated evolution into a limited number of trophic niches

Background  

Damselfishes (Perciformes, Pomacentridae) are a major component of coral reef communities, and the functional diversity of their trophic anatomy is an important constituent of the ecological morphology of these systems. Using shape analyses, biomechanical modelling, and phylogenetically based comparative methods, we examined the anatomy of damselfish feeding among all genera and trophic groups. Coordinate based shape analyses of anatomical landmarks were used to describe patterns of morphological diversity and determine positions of functional groups in a skull morphospace. These landmarks define the lever and linkage structures of the damselfish feeding system, and biomechanical analyses of this data were performed using the software program JawsModel4 in order to calculate the simple mechanical advantage (MA) employed by different skull elements during feeding, and to compute kinematic transmission coefficients (KT) that describe the efficiency with which angular motion is transferred through the complex linkages of damselfish skulls.     

Evolutionary patterns of shape and functional diversification in the skull and jaw musculature of triggerfishes (Teleostei: Balistidae)

The robust skull and highly subdivided adductor mandibulae muscles of triggerfishes provide an excellent system within which to analyze the evolutionary processes underlying phenotypic diversification. We surveyed the anatomical diversity of balistid jaws using Procrustes‐based geometric morphometric analyses and a phylomorphospace approach to quantifying morphological transformation through evolution. We hypothesized that metrics of interspecific cranial shape would reveal patterns of phylogenetic diversification that are congruent with functional and ecological transformation. Morphological landmarks outlining skull and adductor mandibulae muscle shape were collected from 27 triggerfish species. Procrustes‐transformed skull shape configurations revealed significant phylogenetic and size‐influenced structure. Phylomorphospace plots of cranial shape diversity reveal groupings of shape between different species of triggerfish that are mostly consistent with phylogenetic relatedness. Repeated instances of convergence upon similar cranial shape by genetically disparate taxa are likely due to the functional demands of shared specialized dietary habits. This study shows that the diversification of triggerfish skulls occurs via modifications of cranial silhouette and the positioning of subdivided jaw adductor muscles. Using the morphometric data collected here as input to a biomechanical model of triggerfish jaw function, we find that subdivided jaw adductors, in conjunction with a unique cranial skeleton, have direct biomechanical consequences that are not always congruent with phylomorphospace patterns in the triggerfish lineage. The integration of geometric morphometrics with biomechanical modeling in a phylogenetic context provides novel insight into the evolutionary patterns and ecological role of muscle subdivisions in triggerfishes. J. Morphol. 277:737–752, 2016. © 2016 Wiley Periodicals, Inc.     

Foraging mode of a Madagascan iguanian lizard, Oplurus cuvieri cuvieri

The foraging mode of lizards has been extensively studied in members of most major families. One neglected taxon is Opluridae, a distinct Iguanian group endemic to Madagascar. We studied the foraging mode of Oplurus cuvieri cuvieri in a dry forest in north-western Madagascar. Quantitative data on the number of movements per minute and percentage of time spent in moving indicated that O. c. cuvieri is a typical ambusher, although the lizard occasionally fed on plant materials by active foraging. Possible seasonal changes in foraging activity were suggested. The results support the previous view that insectivorous iguanian lizards are ambush foragers.     

Evolution of Cranial Shape in Caecilians (Amphibia: Gymnophiona)

Insights into morphological diversification can be obtained from the ways the species of a clade occupy morphospace. Projecting a phylogeny into morphospace provides estimates of evolutionary trajectories as lineages diversified information that can be used to infer the dynamics of evolutionary processes that produced patterns of morphospace occupation. We present here a large-scale investigation into evolution of morphological variation in the skull of caecilian amphibians, a major clade of vertebrates. Because caecilians are limbless, predominantly fossorial animals, diversification of their skull has occurred within a framework imposed by the functional demands of head-first burrowing. We examined cranial shape in 141 species, over half of known species, using X-ray computed tomography and geometric morphometrics. Mapping an existing phylogeny into the cranial morphospace to estimate the history of morphological change (phylomorphospace), we find a striking pattern: most species occupy distinct clusters in cranial morphospace that closely correspond to the main caecilian clades, and each cluster is separated by unoccupied morphospace. The empty spaces in shape space are unlikely to be caused entirely by extinction or incomplete sampling. The main caecilian clades have different amounts of morphological disparity, but neither clade age nor number of species account for this variation. Cranial shape variation is clearly linked to phyletic divergence, but there is also homoplasy, which is attributed to extrinsic factors associated with head-first digging: features of caecilian crania that have been previously argued to correlate with differential microhabitat use and burrowing ability, such as subterminal and terminal mouths, degree of temporal fenestration (stegokrotaphy/zygokrotaphy), and eyes covered by bone, have evolved and many combinations occur in modern species. We find evidence of morphological convergence in cranial shape, among species that have eyes covered by bone, resulting in a narrow bullet-shaped head. These results reveal a complex history, including early expansion of morphospace and both divergent and convergent evolution resulting in the diversity we observe today.