Iterative adaptive radiations of fossil canids show no evidence for diversity-dependent trait evolution

RRmorph-a new R package to map phenotypic evolutionary rates and patterns on 3D meshes

The study of evolutionary rates and patterns is the key to understand how natural selection shaped the current and past diversity of phenotypes. Phylogenetic comparative methods offer an array of solutions to undertake this challenging task, and help understanding phenotypic variation in full in most circumstances. However, complex, three-dimensional structures such as the skull and the brain serve disparate goals, and different portions of these phenotypes often fulfil different functions, making it hard to understand which parts truly were recruited by natural selection. In the recent past, we developed tools apt to chart evolutionary rate and patterns directly on three-dimensional shapes, according to their magnitude and direction. Here, we present further developments of these tools, which now allow to restitute the mapping of rates and patterns with full biological realism. The tools are condensed in a new R software package.

Mosaic evolution underlies feliform morphological disparity

Constraint is a fundamental concept in evolutionary theory. Morphology and ecology both are limited by functional, historical and developmental factors to a subset of the theoretical range species could occupy. Cat-like carnivorans (Feliformia) offer a unique opportunity to investigate phenotypic constraint, as several feliform clades are purported to be limited to generalized ecomorphological roles, while others possess extremely specialized durophagous (bone-crushing) and sabretooth morphology. We investigated the evolutionary history of feliforms by considering their phylogeny, morphological disparity and rates of evolution. We recover results that show a mosaic pattern exists in the degree of morphological disparity per anatomical region per clade and ecology. Non-hypercarnivores, such as viverrids (civets and genets), Malagasy euplerids and lophocyonids (extinct hypocarnivores), have the greatest dental disparity, while hypercarnivores (felids, nimravids, many hyaenids) have the lowest dental disparity but highest cranial and mandibular disparity (excluding dentition). However, high disparity is not necessarily associated with high rates of evolution, but instead with ecological radiations. We reveal that relationships between specialization and disparity are not as simple as past research has concluded. Instead, morphological disparity results from an anatomical mosaic of evolution, where different ecologies correlate with and likely channel unique patterns/combinations of disparity per anatomical partition.

Limited ecological opportunity influences the tempo of morphological evolution in birds

According to classic models of lineage diversification and adaptive radiation, phenotypic evolution should accelerate in the context of ecological opportunity and slow down when niches become saturated.1,2 However, only weak support for these ideas has been found in nature, perhaps because most analyses make the biologically unrealistic assumption that clade members contribute equally to reducing ecological opportunity, even when they occur in different continents or specialize on different habitats and diets. To view this problem through a different lens, we adapted a new phylogenetic modeling approach that accounts for the fact that competition for ecological opportunity only occurs between species that coexist and share similar habitats and diets. Applying this method to trait data for nearly all extant species of landbirds,3 we find a widespread signature of decelerating trait evolution in lineages adapted to similar habitats or diets. The strength of this pattern was consistent across latitudes when comparing tropical and temperate assemblages. Our results provide little support for the idea that increased diversity and tighter packing of niches accentuates evolutionary slowdowns in the tropics and instead suggest that limited ecological opportunity can be an important factor determining the rate of morphological diversification at a global scale.

The Major Features of Macroevolution

Evolutionary dynamics operating across deep time leave footprints in the shapes of phylogenetic trees. For the last several decades, researchers have used increasingly large and robust phylogenies to study the evolutionary history of individual clades and to investigate the causes of the glaring disparities in diversity among groups. Whereas typically not the focal point of individual clade-level studies, many researchers have remarked on recurrent patterns that have been observed across many different groups and at many different time scales. Whereas previous studies have documented various such regularities in topology and branch length distributions, they have typically focused on a single pattern and used a disparate collection (oftentimes, of quite variable reliability) of trees to assess it. Here we take advantage of modern megaphylogenies and unify previous disparate observations about the shapes embedded in the Tree of Life to create a catalog of the "major features of macroevolution." By characterizing such a large swath of subtrees in a consistent way, we hope to provide a set of phenomena that process-based macroevolutionary models of diversification ought to seek to explain.

Bayesian Prediction of Multivariate Ecology from Phenotypic Data Yields New Insights into the Diets of Extant and Extinct Taxa

AbstractMorphology often reflects ecology, enabling the prediction of ecological roles for taxa that lack direct observations, such as fossils. In comparative analyses, ecological traits, like diet, are often treated as categorical, which may aid prediction and simplify analyses but ignores the multivariate nature of ecological niches. Furthermore, methods for quantifying and predicting multivariate ecology remain rare. Here, we ranked the relative importance of 13 food items for a sample of 88 extant carnivoran mammals and then used Bayesian multilevel modeling to assess whether those rankings could be predicted from dental morphology and body size. Traditional diet categories fail to capture the true multivariate nature of carnivoran diets, but Bayesian regression models derived from living taxa have good predictive accuracy for importance ranks. Using our models to predict the importance of individual food items, the multivariate dietary niche, and the nearest extant analogs for a set of data-deficient extant and extinct carnivoran species confirms long-standing ideas for some taxa but yields new insights into the fundamental dietary niches of others. Our approach provides a promising alternative to traditional dietary classifications. Importantly, this approach need not be limited to diet but serves as a general framework for predicting multivariate ecology from phenotypic traits.

A macroevolutionary pathway to megaherbivory

Several scenarios have been proposed to explain rapid net size increases in some early Cenozoic mammalian lineages: sustained and gradual directional change, successive occupation of adaptive zones associated with progressively larger body sizes, and nondirectional evolution associated with branching events in combination with higher diversification potential of the larger lineages. We test these hypotheses in brontotheres, which are among the first radiations of mammals that consistently evolved multitonne sizes. Body-mass evolution in brontotheres mainly occurred during speciation and had no preferential direction. Long-term directional change stemmed from the higher survival of larger lineages in less-saturated herbivore guilds. Our study emphasizes the role of differential species proliferation in explaining the long-term phenotypic trends observed in the fossil record, which are more than an accumulation of steady microevolutionary changes.

Fossil bacula of five species of Borophaginae (Family: Canidae): Implications for their reproductive biology

The baculum of mammals offers the opportunity to study the reproductive biology of extinct species given that it is a fossilizable part of the male genitalia and that its size and shape correlate with several aspects of the reproductive biology of extant mammals. Fossil bacula, however, are rare. Currently, bacula have been described from only two extinct species of canids, one from the subfamily Caninae and the other from the extinct subfamily Hesperocyoninae. Here, I describe the bacula of five extinct species of Borophaginae, each of which was found with other skeletal elements that have enabled identification to the species level. Two specimens (Aelurodon ferox and Aelurodon stirtoni) are largely complete, while the baculum from Carpocyon compressus is complete but still embedded in matrix that obscures some of its features. The bacula of Paratomarctus euthos and Desmocyon thomsoni are incomplete, but they provide useful information nonetheless. These borophagine bacula are similar to extant canines in being robust, having a urethral groove, and a simple distal end. These features suggest that the Borophaginae had long-lasting copulation and possibly spontaneous ovulation, similar to the extant canines. However, unlike the straight baculum of extant canines, borophagine bacula are ventrally curved (arched), which is also observed in the hesperocyonine baculum. The implication of this curvature for the reproductive biology of these animals remains unknown.

Decoupled evolution of the cranium and mandible in carnivoran mammals

The relationship between skull morphology and diet is a prime example of adaptive evolution. In mammals, the skull consists of the cranium and the mandible. Although the mandible is expected to evolve more directly in response to dietary changes, dietary regimes may have less influence on the cranium because additional sensory and brain-protection functions may impose constraints on its morphological evolution. Here, we tested this hypothesis by comparing the evolutionary patterns of cranium and mandible shape and size across 100+ species of carnivoran mammals with distinct feeding ecologies. Our results show decoupled modes of evolution in cranial and mandibular shape; cranial shape follows clade-based evolutionary shifts, whereas mandibular shape evolution is linked to broad dietary regimes. These results are consistent with previous hypotheses regarding hierarchical morphological evolution in carnivorans and greater evolutionary lability of the mandible with respect to diet. Furthermore, in hypercarnivores, the evolution of both cranial and mandibular size is associated with relative prey size. This demonstrates that dietary diversity can be loosely structured by craniomandibular size within some guilds. Our results suggest that mammal skull morphological evolution is shaped by mechanisms beyond dietary adaptation alone.

Modeling the Evolution of Rates of Continuous Trait Evolution

Rates of phenotypic evolution vary markedly across the tree of life, from the accelerated evolution apparent in adaptive radiations to the remarkable evolutionary stasis exhibited by so-called "living fossils". Such rate variation has important consequences for large-scale evolutionary dynamics, generating vast disparities in phenotypic diversity across space, time, and taxa. Despite this, most methods for estimating trait evolution rates assume rates vary deterministically with respect to some variable of interest or change infrequently during a clade's history. These assumptions may cause underfitting of trait evolution models and mislead hypothesis testing. Here, we develop a new trait evolution model that allows rates to vary gradually and stochastically across a clade. Further, we extend this model to accommodate generally decreasing or increasing rates over time, allowing for flexible modeling of "early/late bursts" of trait evolution. We implement a Bayesian method, termed "evolving rates" (evorates for short), to efficiently fit this model to comparative data. Through simulation, we demonstrate that evorates can reliably infer both how and in which lineages trait evolution rates varied during a clade's history. We apply this method to body size evolution in cetaceans, recovering substantial support for an overall slowdown in body size evolution over time with recent bursts among some oceanic dolphins and relative stasis among beaked whales of the genus Mesoplodon. These results unify and expand on previous research, demonstrating the empirical utility of evorates.

Topographically distinct adaptive landscapes for teeth, skeletons, and size explain the adaptive radiation of Carnivora (Mammalia)

Models of adaptive radiation were originally developed to explain the early, rapid appearance of distinct modes of life within diversifying clades. Phylogenetic tests of this hypothesis have yielded limited support for temporally declining rates of phenotypic evolution across diverse clades, but the concept of an adaptive landscape that links form to fitness, while also crucial to these models, has received more limited attention. Using methods that assess the temporal accumulation of morphological variation and estimate the topography of the underlying adaptive landscape, I found evidence of an early partitioning of mandibulo-dental morphological variation in Carnivora (Mammalia) that occurs on an adaptive landscape with multiple peaks, consistent with classic ideas about adaptive radiation. Although strong support for this mode of adaptive radiation is present in traits related to diet, its signal is not present in body mass data or for traits related to locomotor behavior and substrate use. These findings suggest that adaptive radiations may occur along some axes of ecomorphological variation without leaving a signal in others and that their dynamics are more complex than simple univariate tests might suggest.

Dietary niches of creodonts and carnivorans of the late Eocene Cypress Hills Formation

Modern North American carnivorous mammal assemblages consist of species from a single clade: the Carnivora. Carnivorans once coexisted with members of other meat-eating clades, including the creodonts (Hyaenodontida and Oxyaenida). Creodonts, however, went extinct in North America during the late Eocene and early Oligocene, potentially due to niche overlap and resource competition with contemporary carnivorans. In this study, we employ a community ecology approach to understand whether the dietary niches of coexisting creodonts and carnivorans overlapped during the late Eocene (Chadronian North American Land Mammal Age), a time when creodonts were dwindling and carnivorans were diversifying. We quantify niche overlap based on inferences of diet from carnassial tooth shape estimated using Orientation Patch Count, Dirichlet's Normal Surface Energy, and linear dental measurements as well as from body mass for all species in the Calf Creek Local Fauna of Cypress Hills, Saskatchewan (Treaty 4 land). Although creodonts and carnivorans shared characteristics of their carnassial tooth shape, suggesting similar chewing mechanics and feeding habits, we find that marked differences in body size likely facilitated niche partitioning, at least between the largest creodonts and carnivorans. Calculations of prey focus masses and prey mass spectra indicate that only the smallest creodont may have experienced significant competition for prey with the coeval carnivorans. We suggest that the ultimate extinction of creodonts from North America during the late Eocene and Oligocene was unlikely to have been driven by factors related to niche overlap with carnivorans.

Distributions of extinction times from fossil ages and tree topologies: the example of mid-Permian synapsid extinctions

Given a phylogenetic tree that includes only extinct, or a mix of extinct and extant taxa, where at least some fossil data are available, we present a method to compute the distribution of the extinction time of a given set of taxa under the Fossilized-Birth-Death model. Our approach differs from the previous ones in that it takes into account (i) the possibility that the taxa or the clade considered may diversify before going extinct and (ii) the whole phylogenetic tree to estimate extinction times, whilst previous methods do not consider the diversification process and deal with each branch independently. Because of this, our method can estimate extinction times of lineages represented by a single fossil, provided that they belong to a clade that includes other fossil occurrences. We assess and compare our new approach with a standard previous one using simulated data. Results show that our method provides more accurate confidence intervals. This new approach is applied to the study of the extinction time of three Permo-Carboniferous synapsid taxa (Ophiacodontidae, Edaphosauridae, and Sphenacodontidae) that are thought to have disappeared toward the end of the Cisuralian (early Permian), or possibly shortly thereafter. The timing of extinctions of these three taxa and of their component lineages supports the idea that the biological crisis in the late Kungurian/early Roadian consisted of a progressive decline in biodiversity throughout the Kungurian.

Rapid increase in snake dietary diversity and complexity following the end-Cretaceous mass extinction

The Cenozoic marked a period of dramatic ecological opportunity in Earth history due to the extinction of non-avian dinosaurs as well as to long-term physiographic changes that created new biogeographic theaters and new habitats. Snakes underwent massive ecological diversification during this period, repeatedly evolving novel dietary adaptations and prey preferences. The evolutionary tempo and mode of these trophic ecological changes remain virtually unknown, especially compared with co-radiating lineages of birds and mammals that are simultaneously predators and prey of snakes. Here, we assemble a dataset on snake diets (34,060 observations on the diets of 882 species) to investigate the history and dynamics of the multidimensional trophic niche during the global radiation of snakes. Our results show that per-lineage dietary niche breadths remained remarkably constant even as snakes diversified to occupy disparate outposts of dietary ecospace. Rapid increases in dietary diversity and complexity occurred in the early Cenozoic, and the overall rate of ecospace expansion has slowed through time, suggesting a potential response to ecological opportunity in the wake of the end-Cretaceous mass extinction. Explosive bursts of trophic innovation followed colonization of the Nearctic and Neotropical realms by a group of snakes that today comprises a majority of living snake diversity. Our results indicate that repeated transformational shifts in dietary ecology are important drivers of adaptive radiation in snakes and provide a framework for analyzing and visualizing the evolution of complex ecological phenotypes on phylogenetic trees.

The Cenozoic marked a period of dramatic ecological opportunity in Earth history due to the extinction of non-avian dinosaurs and long-term physiographic changes. This phylogenetic natural history study offers new insights into the evolution of snake ecological diversity after the end-Cretaceous mass extinction, as they took advantage of these new opportunities.

Testing for adaptive radiation: A new approach applied to Madagascar frogs

Adaptive radiation is a key topic at the intersection of ecology and evolutionary biology. Yet the definition and identification of adaptive radiation both remain contentious. Here, we introduce a new approach for identifying adaptive radiations that combines key aspects of two widely used definitions. Our approach compares evolutionary rates in morphology, performance, and diversification between the candidate radiation and other clades. We then apply this approach to a putative adaptive radiation of frogs from Madagascar (Mantellidae). We present new data on morphology and performance from mantellid frogs, then compare rates of diversification and multivariate evolution of size, shape, and performance between mantellids and other frogs. We find that mantellids potentially pass our test for accelerated rates of evolution for shape, but not for size, performance, or diversification. Our results demonstrate that clades can have accelerated phenotypic evolution without rapid diversification (dubbed "adaptive non-radiation"). We also highlight general issues in testing for adaptive radiation, including taxon sampling and the problem of including another adaptive radiation among the comparison clades. Finally, we suggest that similar tests should be conducted on other putative adaptive radiations on Madagascar, comparing their evolutionary rates to those of related clades outside Madagascar. Based on our results, we speculate that older Madagascar clades may show evolutionary patterns more similar to those on a continent than an island.

Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales

Evolutionary models account for either population- or species-level processes but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth–death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology.[Caninae; fossilized birth–death; molecular clock; multispecies coalescent; phylogenetics; species trees.]

Ecological Drivers of Carnivoran Body Shape Evolution

AbstractMorphological diversity is often attributed as adaptations to distinct ecologies. Although biologists have long hypothesized that distinct ecologies drive the evolution of body shape, these relationships are rarely tested across macroevolutionary scales in mammals. Here, I tested hypotheses that locomotor, hunting, and dietary ecologies influenced body shape evolution in carnivorans, a morphologically and ecologically diverse clade of mammals. I found that adaptive models with ecological trait regimes were poor predictors of carnivoran body shape and the underlying morphological components that contribute to body shape variation. Instead, the best-supported model exhibited clade-based evolutionary shifts, indicating that the complexity and variation of body shape landscape cannot be effectively captured by a priori ecological regimes. However, ecological adaptations of body shapes cannot be ruled out, as aquatic and terrestrial carnivorans exhibited opposite allometric patterns of body shape that may be driven by different gravitational constraints associated with these different environments. Similar to body size, body shape is a prominent feature of vertebrate morphology that may transcend one-to-one mapping relationships between morphology and ecological traits, enabling species with distinct body shapes to exploit similar resources and exhibit similar ecologies. Together, these results demonstrate that the multidimensionality of both body shape morphology and ecology makes it difficult to disentangle the complex relationship among morphological evolution, ecological diversity, and phylogeny across macroevolutionary scales.

Testing the occurrence of convergence in the craniomandibular shape evolution of living carnivorans

Convergence consists in the independent evolution of similar traits in distantly related species. The mammalian craniomandibular complex constitutes an ideal biological structure to investigate ecomorphological dynamics and the carnivorans, due to their phenotypic variability and ecological flexibility, offer an interesting case study to explore the occurrence of convergent evolution. Here, we applied multiple pattern‐based metrics to test the occurrence of convergence in the craniomandibular shape of extant carnivorans. To this aim, we tested for convergence in many dietary groups and analyzed several cases of carnivoran convergence concerning either ecologically equivalent species or ecologically similar species of different body sizes described in the literature. Our results validate the occurrence of convergence in ecologically equivalent species in a few cases (as well as in the case of giant and red pandas), but almost never support the occurrence of convergent evolution in dietary categories of living carnivorans. Therefore, convergent evolution in this clade appears to be a rare phenomenon. This is probably the consequence of a complex interplay of one‐to‐many, many‐to‐one, and many‐to‐many relationships taking place between ecology, biomechanics, and morphology.

Species diversification and ecomorphological evolution in the radiation of tanagers (Passeriformes: Thraupidae)

Ecological opportunity is hypothesized to cause an early burst of species diversification and trait evolution followed by a slowdown in diversification rates as niches are filled. Nonetheless, few studies have tested these predictions empirically with ecomorphological data at the large spatial scales relevant to most of biodiversity. Tanagers (Passeriformes: Thraupidae), the largest family of songbirds, show an early burst of species diversification and provide an excellent opportunity to test one of the hallmarks of adaptive radiation: rapid ecomorphological evolution. Here, we test for an early-burst pattern of a resource-exploiting trait (bill morphology) across the radiation of tanagers using a time-calibrated molecular phylogeny and high-resolution three-dimensional surface scans of bill structure from museum study skins. Using recently developed methods of multivariate trait evolution, we find evidence for a rapid burst of bill shape evolution early in the radiation of tanagers, followed by a subsequent decrease in rates toward the present. Likewise, we show that morphological disparity is distributed among (rather than within) subclades, indicating that most of the observed bill shape disparity evolved early in the radiation of tanagers and has slowed through time. The diversification dynamics of tanagers match patterns expected from adaptive radiation and the filling of ecomorphospace.

Exploring adaptive landscapes across deep time: A case study using echinoid body size

Adaptive landscapes are a common way of conceptualizing the phenotypic evolution of lineages across deep time. Although multiple approaches exist to implement this concept into operational models of trait evolution, inferring adaptive landscapes from comparative datasets remains challenging. Here, I explore the macroevolutionary dynamics of echinoid body size using data from over 5000 specimens and a phylogenetic framework incorporating a dense fossil sampling and spanning approximately 270 million years. Furthermore, I implement a novel approach of exploring alternative parameterizations of adaptive landscapes that succeeds in finding simpler, yet better-fitting models. Echinoid body size has been constrained to evolve within a single adaptive optimum for much of the clade's history. However, most of the morphological disparity of echinoids was generated by multiple regime shifts that drove the repeated evolution of miniaturized and gigantic forms. Events of body size innovation occurred predominantly in the Late Cretaceous and were followed by a drastic slowdown following the Cretaceous-Paleogene mass extinction. The discovery of these patterns is contingent upon directly sampling fossil taxa. The macroevolution of echinoid body size is therefore characterized by a late increase in disparity (likely linked to an expansion of ecospace), generated by active processes driving lineages toward extreme morphologies.

Why aren't rabbits and hares larger?

Macroevolutionary consequences of competition among large clades have long been sought in patterns of lineage diversification. However, mechanistically clear examples of such effects remain elusive. Here, we postulated that the limited phenotypic diversity and insular gigantism in lagomorphs could be explained at least in part by an evolutionary constraint placed on them by potentially competing ungulate-type herbivores (UTHs). Our analyses yielded three independent lines of evidence supporting this hypothesis: (1) the minimum UTH body mass is the most influential predictor of the maximum lagomorph body mass in modern ecoregions; (2) the scaling patterns of local-population energy use suggest universal competitive disadvantage of lagomorphs weighing over approximately 6.3 kg against artiodactyls, closely matching their observed upper size limit in continental settings; and (3) the trajectory of maximum lagomorph body mass in North America from the late Eocene to the Pleistocene (37.5-1.5 million years ago) was best modeled by the body mass ceiling placed by the smallest contemporary perissodactyl or artiodactyl. Body size evolution in lagomorphs has likely been regulated by the forces of competition within the clade, increased predation in open habitats, and importantly, competition from other ungulate-type herbivores. Our findings suggest conditionally-coupled dynamics of phenotypic boundaries among multiple clades within an adaptive zone, and highlight the synergy of biotic and abiotic drivers of diversity.

A Total-Group Phylogenetic Metatree for Cetacea and the Importance of Fossil Data in Diversification Analyses

Phylogenetic trees provide a powerful framework for testing macroevolutionary hypotheses, but it is becoming increasingly apparent that inferences derived from extant species alone can be highly misleading. Trees incorporating living and extinct taxa are are needed to address fundamental questions about the origins of diversity and disparity but it has proved challenging to generate robust, species-rich phylogenies that include large numbers of fossil taxa. As a result, most studies of diversification dynamics continue to rely on molecular phylogenies. Here, we extend and apply a recently developed meta-analytic approach for synthesizing previously published phylogenetic studies to infer a well-resolved set of species level, time-scaled phylogenetic hypotheses for extinct and extant cetaceans (whales, dolphins and allies). Our trees extend sampling from the ∼ 90 extant species to over 500 living and extinct species, and therefore allow for more robust inference of macroevolutionary dynamics. While the diversification scenarios we recover are broadly concordant with those inferred from molecular phylogenies they differ in critical ways, notably in the relative contributions of extinction and speciation rate shifts in driving rapid radiations. The metatree approach provides the most immediate route for generating higher level phylogenies of extinct taxa, and opens the door to re-evaluation of macroevolutionary hypotheses derived only from extant taxa.

The signature of competition in ecomorphological traits across the avian radiation

Competition for shared resources represents a fundamental driver of biological diversity. However, the tempo and mode of phenotypic evolution in deep-time has been predominantly investigated using trait evolutionary models which assume that lineages evolve independently from each other. Consequently, the role of species interactions in driving macroevolutionary dynamics remains poorly understood. Here, we quantify the prevalence for signatures of competition between related species in the evolution of ecomorphological traits across the bird radiation. We find that mechanistic trait models accounting for the effect of species interactions on phenotypic divergence provide the best fit for the data on at least one trait axis in 27 out of 59 clades ranging between 21 and 195 species. Where it occurs, the signature of competition generally coincides with positive species diversity-dependence, driven by the accumulation of lineages with similar ecologies, and we find scarce evidence for trait-dependent or negative diversity-dependent phenotypic evolution. Overall, our results suggest that the footprint of interspecific competition is often eroded in long-term patterns of phenotypic diversification, and that other selection pressures may predominantly shape ecomorphological diversity among extant species at macroevolutionary scales.

Ecological Interactions and Macroevolution: A New Field with Old Roots

Linking interspecific interactions (e.g., mutualism, competition, predation, parasitism) to macroevolution (evolutionary change on deep timescales) is a key goal in biology. The role of species interactions in shaping macroevolutionary trajectories has been studied for centuries and remains a cutting-edge topic of current research. However, despite its deep historical roots, classic and current approaches to this topic are highly diverse. Here, we combine historical and contemporary perspectives on the study of ecological interactions in macroevolution, synthesizing ideas across eras to build a zoomed-out picture of the big questions at the nexus of ecology and macroevolution. We discuss the trajectory of this important and challenging field, dividing research into work done before the 1970s, research between 1970 and 2005, and work done since 2005. We argue that in response to long-standing questions in paleobiology, evidence accumulated to date has demonstrated that biotic interactions (including mutualism) can influence lineage diversification and trait evolution over macroevolutionary timescales, and we outline major open questions for future research in the field.

Iterative evolution of large-bodied hypercarnivory in canids benefits species but not clades

Ecological specialization has costs and benefits at various scales: traits benefitting an individual may disadvantage its population, species or clade. In particular, large body size and hypercarnivory (diet over 70% meat) have evolved repeatedly in mammals; yet large hypercarnivores are thought to be trapped in a macroevolutionary "ratchet", marching unilaterally toward decline. Here, we weigh the impact of this specialization on extinction risk using the rich fossil record of North American canids (dogs). In two of three canid subfamilies over the past 40 million years, diversification of large-bodied hypercarnivores appears constrained at the clade level, biasing specialized lineages to extinction. However, despite shorter species durations, extinction rates of large hypercarnivores have been mostly similar to those of all other canids. Extinction was size- and carnivory-selective only at the end of the Pleistocene epoch 11,000 years ago, suggesting that large hypercarnivores were not disadvantaged at the species level before anthropogenic influence.

Selection and Constraints in the Ecomorphological Adaptive Evolution of the Skull of Living Canidae (Carnivora, Mammalia)

The association between phenotype and ecology is essential for understanding the environmental drivers of morphological evolution. This is a particularly challenging task when dealing with complex traits, such as the skull, where multiple selective pressures are at play and evolution might be constrained by ontogenetic and genetic factors. I integrate morphometric tools, comparative methods, and quantitative genetics to investigate how ontogenetic constraints and selection might have interacted during the evolution of the skull in extant Canidae. The results confirm that the evolution of cranial morphology was largely adaptive and molded by changes in diet composition. While the investigation of the adaptive landscape reveals two main selective lines of least resistance (one associated with size and one associated with functional shape features), rates of evolution along size were higher than those found for shape dimensions, suggesting the influence of constraints on morphological evolution. Structural modeling analyses revealed that size, which is the line of most genetic/phenotypic variation, might have acted as a constraint, negatively impacting dietary evolution. Constraints might have been overcome in the case of selection for the consumption of large prey by associating strong selection along both size and shape directions. The results obtained here show that microevolutionary constraints may have played a role in shaping macroevolutionary patterns of morphological evolution.

Exact Distribution of Divergence Times from Fossil Ages and Tree Topologies

Being given a phylogenetic tree of both extant and extinct taxa in which the fossil ages are the only temporal information (namely, in which divergence times are considered unknown), we provide a method to compute the exact probability distribution of any divergence time of the tree with regard to any speciation (cladogenesis), extinction, and fossilization rates under the Fossilized Birth–Death model. We use this new method to obtain a probability distribution for the age of Amniota (the synapsid/sauropsid or bird/mammal divergence), one of the most-frequently used dating constraints. Our results suggest an older age (between about 322 and 340 Ma) than has been assumed by most studies that have used this constraint (which typically assumed a best estimate around 310–315 Ma) and provide, for the first time, a method to compute the shape of the probability density for this divergence time. [Divergence times; fossil ages; fossilized birth–death model; probability distribution.]

Feeding ecology has a stronger evolutionary influence on functional morphology than on body mass in mammals

Ecological specialization is a central driver of adaptive evolution. However, selective pressures may uniquely affect different ecomorphological traits (e.g., size and shape), complicating efforts to investigate the role of ecology in generating phenotypic diversity. Comparative studies can help remedy this issue by identifying specific relationships between ecologies and morphologies, thus elucidating functionally relevant traits. Jaw shape is a dietary correlate that offers considerable insight on mammalian evolution, but few studies have examined the influence of diet on jaw morphology across mammals. To this end, I apply phylogenetic comparative methods to mandibular measurements and dietary data for a diverse sample of mammals. Especially powerful predictors of diet are metrics that capture either the size of the angular process, which increases with greater herbivory, or the length of the posterior portion of the jaw, which decreases with greater herbivory. The size of the angular process likely reflects sizes of attached muscles that produce jaw movements needed to grind plant material. Further, I examine the impact of feeding ecology on body mass, an oft-used ecological surrogate in macroevolutionary studies. Although body mass commonly increases with evolutionary shifts to herbivory, it is outperformed by functional jaw morphology as a predictor of diet. Body mass is influenced by numerous factors beyond diet, and it may be evolutionarily labile relative to functional morphologies. This suggests that ecological diversification events may initially facilitate body mass diversification at smaller taxonomic and temporal scales, but sustained selective pressures will subsequently drive greater trait partitioning in functional morphologies.

Mosaic evolution, preadaptation, and the evolution of evolvability in apes

A major goal in postsynthesis evolutionary biology has been to better understand how complex interactions between traits drive movement along and facilitate the formation of distinct evolutionary pathways. I present analyses of a character matrix sampled across the haplorrhine skeleton that revealed several modules of characters displaying distinct patterns in macroevolutionary disparity. Comparison of these patterns to those in neurological development showed that early ape evolution was characterized by an intense regime of evolutionary and developmental flexibility. Shifting and reduced constraint in apes was met with episodic bursts in phenotypic innovation that built a wide array of functional diversity over a foundation of shared developmental and anatomical structure. Shifts in modularity drove dramatic evolutionary changes across the ape body plan in two distinct ways: (1) an episode of relaxed integration early in hominoid evolution coincided with bursts in evolutionary rate across multiple character suites; (2) the formation of two new trait modules along the branch leading to chimps and humans preceded rapid and dramatic evolutionary shifts in the carpus and pelvis. Changes to the structure of evolutionary mosaicism may correspond to enhanced evolvability that has a "preadaptive" effect by catalyzing later episodes of dramatic morphological remodeling.

A new, fast method to search for morphological convergence with shape data

Morphological convergence is an intensely studied macroevolutionary phenomenon. It refers to the morphological resemblance between phylogenetically distant taxa. Currently available methods to explore evolutionary convergence either: rely on the analysis of the phenotypic resemblance between sister clades as compared to their ancestor, fit different evolutionary regimes to different parts of the tree to see whether the same regime explains phenotypic evolution in phylogenetically distant clades, or assess deviations from the congruence between phylogenetic and phenotypic distances. We introduce a new test for morphological convergence working directly with non-ultrametric (i.e. paleontological) as well as ultrametric phylogenies and multivariate data. The method (developed as the function search.conv within the R package RRphylo) tests whether unrelated clades are morphologically more similar to each other than expected by their phylogenetic distance. It additionally permits using known phenotypes as the most recent common ancestors of clades, taking full advantage of fossil information. We assessed the power of search.conv and the incidence of false positives by means of simulations, and then applied it to three well-known and long-discussed cases of (purported) morphological convergence: the evolution of grazing adaptation in the mandible of ungulates with high-crowned molars, the evolution of mandibular shape in sabertooth cats, and the evolution of discrete ecomorphs among anoles of Caribbean islands. The search.conv method was found to be powerful, correctly identifying simulated cases of convergent morphological evolution in 95% of the cases. Type I error rate is as low as 4–6%. We found search.conv is some three orders of magnitude faster than a competing method for testing convergence.

Fast speciations and slow genes: uncovering the root of living canids

Despite ongoing efforts relying on computationally intensive tree-building methods and large datasets, the deeper phylogenetic relationships between living canid genera remain controversial. We demonstrate that this issue arises fundamentally from the uncertainty of root placement as a consequence of the short length of the branch connecting the major canid clades, which probably resulted from a fast radiation during the early diversification of extant Canidae. Using both nuclear and mitochondrial genes, we investigate the position of the canid root and its consistency by using three rooting methods. We find that mitochondrial genomes consistently retrieve a root node separating the tribe Canini from the remaining canids, whereas nuclear data mostly recover a root that places the Urocyon foxes as the sister lineage of living canids. We demonstrate that, to resolve the canid root, the nuclear segments sequenced so far are significantly less informative than mitochondrial genomes. We also propose that short intervals between speciations obscure the place of the true root, because methods are susceptible to stochastic error in the presence of short internal branches near the root.

The multi-peak adaptive landscape of crocodylomorph body size evolution

BACKGROUND

Little is known about the long-term patterns of body size evolution in Crocodylomorpha, the > 200-million-year-old group that includes living crocodylians and their extinct relatives. Extant crocodylians are mostly large-bodied (3-7 m) predators. However, extinct crocodylomorphs exhibit a wider range of phenotypes, and many of the earliest taxa were much smaller (< 1.2 m). This suggests a pattern of size increase through time that could be caused by multi-lineage evolutionary trends of size increase or by selective extinction of small-bodied species. Here, we characterise patterns of crocodylomorph body size evolution using a model fitting-approach (with cranial measurements serving as proxies). We also estimate body size disparity through time and quantitatively test hypotheses of biotic and abiotic factors as potential drivers of crocodylomorph body size evolution.

RESULTS

Crocodylomorphs reached an early peak in body size disparity during the Late Jurassic, and underwent an essentially continual decline since then. A multi-peak Ornstein-Uhlenbeck model outperforms all other evolutionary models fitted to our data (including both uniform and non-uniform), indicating that the macroevolutionary dynamics of crocodylomorph body size are better described within the concept of an adaptive landscape, with most body size variation emerging after shifts to new macroevolutionary regimes (analogous to adaptive zones). We did not find support for a consistent evolutionary trend towards larger sizes among lineages (i.e., Cope's rule), or strong correlations of body size with climate. Instead, the intermediate to large body sizes of some crocodylomorphs are better explained by group-specific adaptations. In particular, the evolution of a more aquatic lifestyle (especially marine) correlates with increases in average body size, though not without exceptions.

CONCLUSIONS

Shifts between macroevolutionary regimes provide a better explanation of crocodylomorph body size evolution on large phylogenetic and temporal scales, suggesting a central role for lineage-specific adaptations rather than climatic forcing. Shifts leading to larger body sizes occurred in most aquatic and semi-aquatic groups. This, combined with extinctions of groups occupying smaller body size regimes (particularly during the Late Cretaceous and Cenozoic), gave rise to the upward-shifted body size distribution of extant crocodylomorphs compared to their smaller-bodied terrestrial ancestors.

Conserving evolutionary history does not result in greater diversity over geological time scales

Alternative prioritization strategies have been proposed to safeguard biodiversity over macroevolutionary time scales. The first prioritizes the most distantly related species—maximizing phylogenetic diversity (PD)—in the hopes of capturing at least some lineages that will successfully diversify into the future. The second prioritizes lineages that are currently speciating, in the hopes that successful lineages will continue to generate species into the future. These contrasting schemes also map onto contrasting predictions about the role of slow diversifiers in the production of biodiversity over palaeontological time scales. We consider the performance of the two schemes across 10 dated species-level palaeo-phylogenetic trees ranging from Foraminifera to dinosaurs. We find that prioritizing PD for conservation generally led to fewer subsequent lineages, while prioritizing diversifiers led to modestly more subsequent diversity, compared with random sets of lineages. Importantly for conservation, the tree shape when decisions are made cannot predict which scheme will be most successful. These patterns are inconsistent with the notion that long-lived lineages are the source of new species. While there may be sound reasons for prioritizing PD for conservation, long-term species production might not be one of them.

Evolutionary shifts in extant mustelid (Mustelidae: Carnivora) cranial shape, body size and body shape coincide with the Mid-Miocene Climate Transition

Environmental changes can lead to evolutionary shifts in phenotypic traits, which in turn facilitate the exploitation of novel adaptive landscapes and lineage diversification. The global cooling, increased aridity and expansion of open grasslands during the past 50 Myr are prime examples of new adaptive landscapes that spurred lineage and ecomorphological diversity of several mammalian lineages such as rodents and large herbivorous megafauna. However, whether these environmental changes facilitated evolutionary shifts in small- to mid-sized predator morphology is unknown. Here, I used a complete cranial and body morphological dataset to examine the timing of evolutionary shifts in cranial shape, body size and body shape within extant mustelids (martens, otters, polecats and weasels) during the climatic and environmental changes of the Cenozoic. I found that evolutionary shifts in all three traits occurred within extant mustelid subclades just after the onset of the Mid-Miocene Climate Transition. These mustelid subclades first shifted towards more elongate body plans followed by concurrent shifts towards smaller body sizes and more robust crania. I hypothesize that these cranial and body morphological shifts enabled mustelids to exploit novel adaptive zones associated with the climatic and environmental changes of the Mid to Late Miocene, which facilitated significant increases in clade carrying capacity.

Signatures of echolocation and dietary ecology in the adaptive evolution of skull shape in bats

Morphological diversity may arise rapidly as a result of adaptation to novel ecological opportunities, but early bursts of trait evolution are rarely observed. Rather, models of discrete shifts between adaptive zones may better explain macroevolutionary dynamics across radiations. To investigate which of these processes underlie exceptional levels of morphological diversity during ecological diversification, we use modern phylogenetic tools and 3D geometric morphometric datasets to examine adaptive zone shifts in bat skull shape. Here we report that, while disparity was established early, bat skull evolution is best described by multiple adaptive zone shifts. Shifts are partially decoupled between the cranium and mandible, with cranial evolution more strongly driven by echolocation than diet. Phyllostomidae, a trophic adaptive radiation, exhibits more adaptive zone shifts than all other families combined. This pattern was potentially driven by ecological opportunity and facilitated by a shift to intermediate cranial shapes compared to oral-emitters and other nasal emitters.

What drives changes in morphological diversity? Here, Arbour et al. analyse skull 3D shape evolution across the bat radiation using µCT scan data, finding two phases of skull shape diversification, early adaptive shifts related to echolocation, and more recent shifts related to diet transitions.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Simultaneous detection of macroevolutionary patterns in phenotypic means and rate of change with and within phylogenetic trees including extinct species

Recognizing evolutionary trends in phenotypic means and rates requires the application of phylogenetic comparative methods (PCMs). Most PCMs are unsuited to make full use of fossil information, which is a drawback, given the inclusion of such data improves, and in some cases even corrects, the proper understanding of trait evolution. Here we present a new computer application, written in R, that allows the simultaneous computation of temporal trends in phenotypic mean and evolutionary rate along a phylogeny, and to contrast such patterns among different clades within the tree. By using simulation experiments, we show the new implementation, names search.trend is as powerful as existing PCM tools in discerning macroevolutionary patterns in phenotypic means and rates, but differently from any other PCM allows comparing individual clades to each other, and provides rich information about trait evolution for all lineages in the tree.

Early Arrival and Climatically-Linked Geographic Expansion of New World Monkeys from Tiny African Ancestors

New World Monkeys (NWM) (platyrrhines) are one of the most diverse groups of primates, occupying today a wide range of ecosystems in the American tropics and exhibiting large variations in ecology, morphology, and behavior. Although the relationships among the almost 200 living species are relatively well understood, we lack robust estimates of the timing of origin, ancestral morphology, and geographic range evolution of the clade. Herein, we integrate paleontological and molecular evidence to assess the evolutionary dynamics of extinct and extant platyrrhines. We develop novel analytical frameworks to infer the evolution of body mass, changes in latitudinal ranges through time, and species diversification rates using a phylogenetic tree of living and fossil taxa. Our results show that platyrrhines originated 5–10 million years earlier than previously assumed, dating back to the Middle Eocene. The estimated ancestral platyrrhine was small—weighing 0.4 kg—and matched the size of their presumed African ancestors. As the three platyrrhine families diverged, we recover a rapid change in body mass range. During the Miocene Climatic Optimum, fossil diversity peaked and platyrrhines reached their widest latitudinal range, expanding as far South as Patagonia, favored by warm and humid climate and the lower elevation of the Andes. Finally, global cooling and aridification after the middle Miocene triggered a geographic contraction of NWM and increased their extinction rates. These results unveil the full evolutionary trajectory of an iconic and ecologically important radiation of monkeys and showcase the necessity of integrating fossil and molecular data for reliably estimating evolutionary rates and trends.

Post-Cretaceous bursts of evolution along the benthic-pelagic axis in marine fishes

Ecological opportunity arising in the aftermath of mass extinction events is thought to be a powerful driver of evolutionary radiations. Here, we assessed how the wake of the Cretaceous-Palaeogene (K-Pg) mass extinction shaped diversification dynamics in a clade of mostly marine fishes (Carangaria), which comprises a disparate array of benthic and pelagic dwellers including some of the most astonishing fish forms (e.g. flatfishes, billfishes, remoras, archerfishes). Analyses of lineage diversification show time-heterogeneous rates of lineage diversification in carangarians, with highest rates reached during the Palaeocene. Likewise, a remarkable proportion of Carangaria's morphological variation originated early in the history of the group and in tandem with a marked incidence of habitat shifts. Taken together, these results suggest that all major lineages and body plans in Carangaria originated in an early burst shortly after the K-Pg mass extinction, which ultimately allowed the occupation of newly released niches along the benthic-pelagic habitat axis.

Size does matter: Parallel evolution of adaptive thermal tolerance and body size facilitates adaptation to climate change in domestic cattle

The adaptive potential of livestock under a warming climate is increasingly relevant in relation to the growing pressure of global food security. Studies on heat tolerance demonstrate the interplay of adaptation and acclimatization in functional traits, for example, a reduction in body size and enhanced tolerance in response to a warming climate. However, current lack of understanding of functional traits and phylogenetic history among phenotypically distinct populations constrains predictions of climate change impact. Here, we demonstrate evidence of parallel evolution in adaptive tolerance to heat stress in dwarf cattle breeds (DCB, Bos taurus indicus) and compare their thermoregulatory responses with those in standard size cattle breeds (SCB, crossbred, Bos taurus indicus × Bos taurus taurus). We measured vital physiological, hematological, biochemical, and gene expression changes in DCB and SCB and compared the molecular phylogeny using mitochondrial genome (mitogenome) analysis. Our results show that SCB can acclimatize in the short term to higher temperatures but reach their tolerance limit under prevailing tropical conditions, while DCB is adapted to the warmer climate. Increased hemoglobin concentration, reduced cellular size, and smaller body size enhance thermal tolerance. Mitogenome analysis revealed that different lineages of DCB have evolved reduced size independently, as a parallel adaptation to heat stress. The results illustrate mechanistic ways of dwarfing, body size-dependent tolerance, and differential fitness in a large mammal species under harsh field conditions, providing a background for comparing similar populations during global climate change. These demonstrate the value of studies combining functional, physiological, and evolutionary approaches to delineate adaptive potential and plasticity in domestic species. We thus highlight the value of locally adapted breeds as a reservoir of genetic variation contributing to the global domestic genetic resource pool that will become increasingly important for livestock production systems under a warming climate.

Evolution of morphological integration in the skull of Carnivora (Mammalia): Changes in Canidae lead to increased evolutionary potential of facial traits

Morphological integration refers to the fact that different phenotypic traits of organisms are not fully independent from each other, and tend to covary to different degrees. The covariation among traits is thought to reflect properties of the species' genetic architecture and thus can have an impact on evolutionary responses. Furthermore, if morphological integration changes along the history of a group, inferences of past selection regimes might be problematic. Here, we evaluated the stability and evolution of the morphological integration of skull traits in Carnivora by using evolutionary simulations and phylogenetic comparative methods. Our results show that carnivoran species are able to respond to natural selection in a very similar way. Our comparative analyses show that the phylogenetic signal for pattern of integration is lower than that observed for morphology (trait averages), and that integration was stable throughout the evolution of the group. That notwithstanding, Canidae differed from other families by having higher integration, evolvability, flexibility, and allometric coefficients on the facial region. These changes might have allowed canids to rapidly adapt to different food sources, helping to explain not only the phenotypic diversification of the family, but also why humans were able to generate such a great diversity of dog breeds through artificial selection.