mv morph : an r package for fitting multivariate evolutionary models to morphometric data

Ecological diversification preceded geographical expansion during the evolutionary radiation of Cataglyphis desert ants

Biological diversity often arises as organisms adapt to new ecological conditions (i.e., ecological opportunities) or colonize suitable areas (i.e., spatial opportunities). Cases of geographical expansion followed by local ecological divergence are well described; they result in clades comprising ecologically heterogeneous subclades. Here, we show that the desert ant genus Cataglyphis likely originated in open grassland habitats in the Middle East ∼18 million years ago and became a taxon of diverse species specializing in prey of different masses. The genus then colonized the Mediterranean Basin around 9 million years ago. The result was the rapid accumulation of species, and the appearance of local assemblages containing species from different lineages that still displayed ancestral foraging specialties. These findings highlight that, in Cataglyphis, ecological diversification preceded geographical expansion, resulting in a clade composed of ecologically homogeneous subclades.

The rise of pelagic sharks and adaptive evolution of pectoral fin morphology during the Cretaceous

The emergence and subsequent evolution of pectoral fins is a key point in vertebrate evolution, as pectoral fins are dominant control surfaces for locomotion in extant fishes.1,2,3 However, major gaps remain in our understanding of the diversity and evolution of pectoral fins among cartilaginous fishes (Chondrichthyes), a group with an evolutionary history spanning over 400 million years with current selachians (modern sharks) appearing about 200 million years ago.4,5,6 Modern sharks are a charismatic group of vertebrates often thought to be predators roaming the open ocean and coastal areas, but most extant species occupy the seafloor.4 Here we use an integrative approach to understand what facilitated the expansion to the pelagic realm and what morphological changes accompanied this shift. On the basis of comparative analyses in the framework of a time-calibrated molecular phylogeny,7 we show that modern sharks expanded to the pelagic realm no later than the Early Cretaceous (Barremian). The pattern of pectoral fin aspect ratios across selachians is congruent with adaptive evolution, and we identify an increase of the subclade disparity of aspect ratio at a time when sea surface temperatures were at their highest.8 The expansion to open ocean habitats likely involved extended bouts of sustained fast swimming, which led to the selection for efficient movement via higher aspect ratio pectoral fins. Swimming performance was likely enhanced in pelagic sharks during this time due to the elevated temperatures in the sea, highlighting that shark evolution has been greatly impacted by climate change.

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instan-taneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny [30]. Here, we introduce a new family of models - the so-called "Phylogenetic Integrated Velocity" (PIV) models - that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the U.S.A. indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Robust Phylogenetic Regression

Modern comparative biology owes much to phylogenetic regression. At its conception, this technique sparked a revolution that armed biologists with phylogenetic comparative methods (PCMs) for disentangling evolutionary correlations from those arising from hierarchical phylogenetic relationships. Over the past few decades, the phylogenetic regression framework has become a paradigm of modern comparative biology that has been widely embraced as a remedy for shared ancestry. However, recent evidence has shown doubt over the efficacy of phylogenetic regression, and PCMs more generally, with the suggestion that many of these methods fail to provide an adequate defense against unreplicated evolution-the primary justification for using them in the first place. Importantly, some of the most compelling examples of biological innovation in nature result from abrupt lineage-specific evolutionary shifts, which current regression models are largely ill equipped to deal with. Here we explore a solution to this problem by applying robust linear regression to comparative trait data. We formally introduce robust phylogenetic regression to the PCM toolkit with linear estimators that are less sensitive to model violations than the standard least-squares estimator, while still retaining high power to detect true trait associations. Our analyses also highlight an ingenuity of the original algorithm for phylogenetic regression based on independent contrasts, whereby robust estimators are particularly effective. Collectively, we find that robust estimators hold promise for improving tests of trait associations and offer a path forward in scenarios where classical approaches may fail. Our study joins recent arguments for increased vigilance against unreplicated evolution and a better understanding of evolutionary model performance in challenging-yet biologically important-settings.

Evolutionary trends in the elasmobranch neurocranium

The neurocranium (braincase) is one of the defining vertebrate characters. Housing the brain and other key sensory organs, articulating with the jaws and contributing to the shape of the anteriormost portion of the body, the braincase is undoubtedly of great functional importance. Through studying relationships between braincase shape and ecology we can gain an improved understanding of form-function relationships in extant and fossil taxa. Elasmobranchii (sharks and rays) represent an important case study of vertebrate braincase diversity as their neurocranium is simplified and somewhat decoupled from other components of the cranium relative to other vertebrates. Little is known about the associations between ecology and braincase shape in this clade. In this study we report patterns of mosaic cranial evolution in Elasmobranchii that differ significantly from those present in other clades. The degree of evolutionary modularity also differs between Selachii and Batoidea. In both cases innovation in the jaw suspension appears to have driven shifts in patterns of integration and modularity, subsequently facilitating ecological diversification. Our results confirm the importance of water depth and biogeography as drivers of elasmobranch cranial diversity and indicate that skeletal articulation between the neurocranium and jaws represents a major constraint upon the evolution of braincase shape in vertebrates.

Body size shapes song in honeyeaters

Birdsongs are among the most distinctive animal signals. Their evolution is thought to be shaped simultaneously by habitat structure and by the constraints of morphology. Habitat structure affects song transmission and detectability, thus influencing song (the acoustic adaptation hypothesis), while body size and beak size and shape necessarily constrain song characteristics (the morphological constraint hypothesis). Yet, support for the acoustic adaptation and morphological constraint hypotheses remains equivocal, and their simultaneous examination is infrequent. Using a phenotypically diverse Australasian bird clade, the honeyeaters (Aves: Meliphagidae), we compile a dataset consisting of song, environmental, and morphological variables for 163 species and jointly examine predictions of these two hypotheses. Overall, we find that body size constrains song frequency and pace in honeyeaters. Although habitat type and environmental temperature influence aspects of song, that influence is indirect, likely via effects of environmental variation on body size, with some evidence that elevation constrains the evolution of song peak frequency. Our results demonstrate that morphology has an overwhelming influence on birdsong, in support of the morphological constraint hypothesis, with the environment playing a secondary role generally via body size rather than habitat structure. These results suggest that changing body size (a consequence of both global effects such as climate change and local effects such as habitat transformation) will substantially influence the nature of birdsong.

Unexpectedly uneven distribution of functional trade-offs explains cranial morphological diversity in carnivores

Functional trade-offs can affect patterns of morphological and ecological evolution as well as the magnitude of morphological changes through evolutionary time. Using morpho-functional landscape modelling on the cranium of 132 carnivore species, we focused on the macroevolutionary effects of the trade-off between bite force and bite velocity. Here, we show that rates of evolution in form (morphology) are decoupled from rates of evolution in function. Further, we found theoretical morphologies optimising for velocity to be more diverse, while a much smaller phenotypic space was occupied by shapes optimising force. This pattern of differential representation of different functions in theoretical morphological space was highly correlated with patterns of actual morphological disparity. We hypothesise that many-to-one mapping of cranium shape on function may prevent the detection of direct relationships between form and function. As comparatively only few morphologies optimise bite force, species optimising this function may be less abundant because they are less likely to evolve. This, in turn, may explain why certain clades are less variable than others. Given the ubiquity of functional trade-offs in biological systems, these patterns may be general and may help to explain the unevenness of morphological and functional diversity across the tree of life.

The hierarchical radiation of phyllostomid bats as revealed by adaptive molar morphology

Adaptive radiations are bursts in biodiversity that generate new evolutionary lineages and phenotypes. However, because they typically occur over millions of years, it is unclear how their macroevolutionary dynamics vary through time and among groups of organisms. Phyllostomid bats radiated extensively for diverse diets-from insects to vertebrates, fruit, nectar, and blood-and we use their molars as a model system to examine the dynamics of adaptive radiations. Three-dimensional shape analyses of lower molars of Noctilionoidea (Phyllostomidae and close relatives) indicate that different diet groups exhibit distinct morphotypes. Comparative analyses further reveal that phyllostomids are a striking example of a hierarchical radiation; phyllostomids' initial, higher-level diversification involved an "early burst" in molar morphological disparity as lineages invaded new diet-affiliated adaptive zones, followed by subsequent lower-level diversifications within adaptive zones involving less dramatic morphological changes. We posit that strong selective pressures related to initial shifts to derived diets may have freed molars from morpho-functional constraints associated with the ancestral molar morphotype. Then, lineages with derived diets (frugivores and nectarivores) diversified within broad adaptive zones, likely reflecting finer-scale niche partitioning. Importantly, the observed early burst pattern is only evident when examining molar traits that are strongly linked to diet, highlighting the value of ecomorphological traits in comparative studies. Our results support the hypothesis that adaptive radiations are commonly hierarchical and involve different tempos and modes at different phylogenetic levels, with early bursts being more common at higher levels.

Opening the Museum's Vault: Historical Field Records Preserve Reliable Ecological Data

AbstractMuseum specimens have long served as foundational data sources for ecological, evolutionary, and environmental research. Continued reimagining of museum collections is now also generating new types of data associated with but beyond physical specimens, a concept known as "extended specimens." Field notes penned by generations of naturalists contain firsthand ecological observations associated with museum collections and comprise a form of extended specimens with the potential to provide novel ecological data spanning broad geographic and temporal scales. Despite their data-yielding potential, however, field notes remain underutilized in research because of their heterogeneous, unstandardized, and qualitative nature. We introduce an approach for transforming descriptive ecological notes into quantitative data suitable for statistical analysis. Tests with simulated and real-world published data show that field notes and our transformation approach retain reliable quantitative ecological information under a range of sample sizes and evolutionary scenarios. Unlocking the wealth of data contained within field records could facilitate investigations into the ecology of clades whose diversity, distribution, or other demographic features present challenges to traditional ecological studies, improve our understanding of long-term environmental and evolutionary change, and enhance predictions of future change.

The origins and drivers of sexual size dimorphism in sharks

While sexual size dimorphism (SSD) is abundant in nature, there is huge variation in both the intensity and direction of SSD. SSD results from a combination of sexual selection for large male size, fecundity selection for large female size and ecological selection for either. In most vertebrates, it is variation in the intensity of male-male competition that primarily underlies variation in SSD. In this study, we test four hypotheses regarding the adaptive value of SSD in sharks-considering the potential for each of fecundity, sexual, ecological selection and reproductive mode as the primary driver of variation in SSD between species. We also estimate past macroevolutionary shifts in SSD direction/intensity through shark phylogeny. We were unable to find evidence of significant SSD in early sharks and hypothesise that SSD is a derived state in this clade, that has evolved independently of SSD observed in other vertebrates. Moreover, there is no significant relationship between SSD and fecundity, testes mass or oceanic depth in sharks. However, there is evidence to support previous speculation that reproductive mode is an important determinant of interspecific variation in SSD in sharks. This is significant as in most vertebrates sexual selection is thought to be the primary driver of SSD trends, with evidence for the role of fecundity selection in other clades being inconsistent at best. While the phylogenetic distribution of SSD among sharks is superficially similar to that observed in other vertebrate clades, the relative importance of selective pressures underlying its evolution appears to differ.

Analytical advances alleviate model misspecification in non-Brownian multivariate comparative methods

Adams and Collyer argue that contemporary multivariate (Gaussian) phylogenetic comparative methods are prone to favouring more complex models of evolution and sometimes rotation invariance can be an issue. Here we dissect the concept of rotation invariance and point out that, depending on the understanding, this can be an issue with any method that relies on numerical instead of analytical estimation approaches. We relate this to the ongoing discussion concerning phylogenetic principal component analysis. Contrary to what Adams and Collyer found, we do not observe a bias against the simpler Brownian motion process in simulations when we use the new, improved, likelihood evaluation algorithm employed by mvSLOUCH, which allows for studying much larger phylogenies and more complex model setups.

Predatory synapsid ecomorphology signals growing dynamism of late Palaeozoic terrestrial ecosystems

Terrestrial ecosystems evolved substantially through the Palaeozoic, especially the Permian, gaining much new complexity, especially among predators. Key among these predators were non-mammalian synapsids. Predator ecomorphology reflect interactions with prey and competitors, which are key controls on carnivore diversity and ecology. Therefore, carnivorous synapsids may offer insight on wider ecological evolution as the first complex, tetrapod-dominated, terrestrial ecosystems formed through the late Palaeozoic. Using morphometric and phylogenetic comparative methods, we chart carnivorous synapsid trophic morphology from the latest Carboniferous to the earliest Triassic (307-251.2 Ma). We find a major morphofunctional shift in synapsid carnivory between the early and middle Permian, via the addition of new feeding modes increasingly specialised for greater biting power or speed that captures the growing antagonism and dynamism of terrestrial tetrapod predator-prey interactions. The further evolution of new hypo- and hypercarnivorous synapsids highlight the nascent intrinsic pressures and complexification of terrestrial ecosystems across the mid-late Permian.

Exploring the Climatic Niche Evolution of the Genus Falco (Aves: Falconidae) in Europe

By integrating species distribution modeling techniques, phylogenetic comparative methods, and climatic data, we analyzed how European falcon climatic niches have changed over evolutionary time in order to understand their tempo and mode of evolution and gain phylogenetic insights related to the ecological context of falcon evolution. For this purpose, we tested the relative contributions of niche conservatism, convergent evolution, and divergent evolution in the evolutionary history of this group of species in Europe. The occupation of climatic niche spaces by falcon species in Europe was not similar, considering that their climatic niche evolution was characterized by heterotachy, especially after ca. 4 Mya. Our results indicate that convergent evolution and niche divergence played an important role in the evolutionary history of these species, with no significant evidence of closely related species retaining their fundamental niche over time (phylogenetic niche conservatism). In most analyses, less closely related falcon species occupied similar climatic environments. We found that speciation in the European genus Falco was influenced by climatic niche differentiation, more prevalent in the last 4 million years, with the main climatic niche shifts occurring between closely related falcon species.

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.

Drivers of morphological evolution in the toothed whale jaw

Toothed whales (odontocetes) emit high-frequency underwater sounds (echolocate)-an extreme and unique innovation allowing them to sense their prey and environment. Their highly specialized mandible (lower jaw) allows high-frequency sounds to be transmitted back to the inner ear. Echolocation is evident in the earliest toothed whales, but little research has focused on the evolution of mandibular form regarding this unique adaptation. Here, we use a high-density, three-dimensional geometric morphometric analysis of 100 living and extinct cetacean species spanning their ∼50-million-year evolutionary history. Our analyses demonstrate that most shape variation is found in the relative length of the jaw and the mandibular symphysis. The greatest morphological diversity was obtained during two periods of rapid evolution: the initial evolution of archaeocetes (stem whales) in the early to mid-Eocene as they adapted to an aquatic lifestyle, representing one of the most extreme adaptive transitions known, and later on in the mid-Oligocene odontocetes as they became increasingly specialized for a range of diets facilitated by increasingly refined echolocation. Low disparity in the posterior mandible suggests the shape of the acoustic window, which receives sound, has remained conservative since the advent of directional hearing in the aquatic archaeocetes, even as the earliest odontocetes began to receive sounds from echolocation. Diet, echolocation, feeding method, and dentition type strongly influence mandible shape. Unlike in the toothed whale cranium, we found no significant asymmetry in the mandible. We suggest that a combination of refined echolocation and associated dietary specializations have driven morphology and disparity in the toothed whale mandible.

Predicting evolutionary targets and parameters of gene deletion from expression data

Motivation

Gene deletion is traditionally thought of as a nonadaptive process that removes functional redundancy from genomes, such that it generally receives less attention than duplication in evolutionary turnover studies. Yet, mounting evidence suggests that deletion may promote adaptation via the "less-is-more" evolutionary hypothesis, as it often targets genes harboring unique sequences, expression profiles, and molecular functions. Hence, predicting the relative prevalence of redundant and unique functions among genes targeted by deletion, as well as the parameters underlying their evolution, can shed light on the role of gene deletion in adaptation.

Results

Here, we present CLOUDe, a suite of machine learning methods for predicting evolutionary targets of gene deletion events from expression data. Specifically, CLOUDe models expression evolution as an Ornstein-Uhlenbeck process, and uses multi-layer neural network, extreme gradient boosting, random forest, and support vector machine architectures to predict whether deleted genes are "redundant" or "unique", as well as several parameters underlying their evolution. We show that CLOUDe boasts high power and accuracy in differentiating between classes, and high accuracy and precision in estimating evolutionary parameters, with optimal performance achieved by its neural network architecture. Application of CLOUDe to empirical data from Drosophila suggests that deletion primarily targets genes with unique functions, with further analysis showing these functions to be enriched for protein deubiquitination. Thus, CLOUDe represents a key advance in learning about the role of gene deletion in functional evolution and adaptation.

Availability and implementation

CLOUDe is freely available on GitHub (https://github.com/anddssan/CLOUDe).

phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things)

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Uncovering the mosaic evolution of the carnivoran skeletal system

The diversity of vertebrate skeletons is often attributed to adaptations to distinct ecological factors such as diet, locomotion, and sensory environment. Although the adaptive evolution of skull, appendicular skeleton, and vertebral column is well studied in vertebrates, comprehensive investigations of all skeletal components simultaneously are rarely performed. Consequently, we know little of how modes of evolution differ among skeletal components. Here, we tested if ecological and phylogenetic effects led to distinct modes of evolution among the cranial, appendicular and vertebral regions in extant carnivoran skeletons. Using multivariate evolutionary models, we found mosaic evolution in which only the mandible, hindlimb and posterior (i.e. last thoracic and lumbar) vertebrae showed evidence of adaptation towards ecological regimes whereas the remaining skeletal components reflect clade-specific evolutionary shifts. We hypothesize that the decoupled evolution of individual skeletal components may have led to the origination of distinct adaptive zones and morphologies among extant carnivoran families that reflect phylogenetic hierarchies. Overall, our work highlights the importance of examining multiple skeletal components simultaneously in ecomorphological analyses. Ongoing work integrating the fossil and palaeoenvironmental record will further clarify deep-time drivers that govern the carnivoran diversity we see today and reveal the complexity of evolutionary processes in multicomponent systems.

The Cauchy Process on Phylogenies: A Tractable Model for Pulsed Evolution

Phylogenetic comparative methods use random processes, such as the Brownian Motion, to model the evolution of continuous traits on phylogenetic trees. Growing evidence for non-gradual evolution motivated the development of complex models, often based on Lévy processes. However, their statistical inference is computationally intensive and currently relies on approximations, high-dimensional sampling, or numerical integration. We consider here the Cauchy Process (CP), a particular pure-jump Lévy process in which the trait increment along each branch follows a centered Cauchy distribution with a dispersion proportional to its length. In this work, we derive an exact algorithm to compute both the joint probability density of the tip trait values of a phylogeny under a CP and the ancestral trait values and branch increments posterior densities in quadratic time. A simulation study shows that the CP generates patterns in comparative data that are distinct from any Gaussian process, and that restricted maximum likelihood parameter estimates and root trait reconstruction are unbiased and accurate for trees with 200 tips or less. The CP has only two parameters but is rich enough to capture complex-pulsed evolution. It can reconstruct posterior ancestral trait distributions that are multimodal, reflecting the uncertainty associated with the inference of the evolutionary history of a trait from extant taxa only. Applied on empirical datasets taken from the Evolutionary Ecology and Virology literature, the CP suggests nuanced scenarios for the body size evolution of Greater Antilles Lizards and for the geographical spread of the West Nile Virus epidemics in North America, both consistent with previous studies using more complex models. The method is efficiently implemented in C with an R interface in package cauphy, which is open source and freely available online.

PhyloJunction: a computational framework for simulating, developing, and teaching evolutionary models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, through its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This paper describes the features of PhyloJunction - which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models - and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

How hibernation in frogs drives brain and reproductive evolution in opposite directions

Environmental seasonality can promote the evolution of larger brains through cognitive and behavioral flexibility but can also hamper it when temporary food shortage is buffered by stored energy. Multiple hypotheses linking brain evolution with resource acquisition and allocation have been proposed for warm-blooded organisms, but it remains unclear how these extend to cold-blooded taxa whose metabolism is tightly linked to ambient temperature. Here, we integrated these hypotheses across frogs and toads in the context of varying brumation (hibernation) durations and their environmental correlates. We showed that protracted brumation covaried negatively with brain size but positively with reproductive investment, likely in response to brumation-dependent changes in the socio-ecological context and associated selection on different tissues. Our results provide novel insights into resource allocation strategies and possible constraints in trait diversification, which may have important implications for the adaptability of species under sustained environmental change.

Inferring the ecological and evolutionary determinants of community genetic diversity

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Cochlea development shapes bat sensory system evolution

Sensory organs must develop alongside the skull within which they are largely encased, and this relationship can manifest as the skull constraining the organs, organs constraining the skull, or organs constraining one another in relative size. How this interplay between sensory organs and the developing skull plays out during the evolution of sensory diversity; however, remains unknown. Here, we examine the developmental sequence of the cochlea, the organ responsible for hearing and echolocation, in species with distinct diet and echolocation types within the ecologically diverse bat super-family Noctilionoidea. We found the size and shape of the cochlea largely correlates with skull size, with exceptions of Pteronotus parnellii, whose high duty cycle echolocation (nearly constant emission of sound pulses during their echolocation process allowing for detailed information gathering, also called constant frequency echolocation) corresponds to a larger cochlear and basal turn, and Monophyllus redmani, a small-bodied nectarivorous bat, for which interactions with other sensory organs restrict cochlea size. Our findings support the existence of developmental constraints, suggesting that both developmental and anatomical factors may act synergistically during the development of sensory systems in noctilionoid bats.

Phenotypic rate and state are decoupled in response to river-to-lake transitions in cichlid fishes

Geographic access to isolated ecosystems is an important catalyst of adaptive radiation. Cichlid fishes repeatedly colonized rift, crater, and volcanic lakes from surrounding rivers. We test the "lake effect" on the phenotypic rate and state across 253 cichlid species. The rate of evolution was consistently higher (~10-fold) in lakes, and consistent across different dimensions of the phenotype. Rate shifts tended to occur coincident with or immediately following river-to-lake transitions, generally resulting in 2- to 5-fold faster rates than in the founding riverine lineage. By contrast, river- and lake-dwelling cichlids exhibit considerable overlap in phenotypes, generally with less disparity in lakes, but often different evolutionary optima. Taken together, these results suggest that lake radiations rapidly expand into niches largely already represented by ancestral riverine lineages, albeit in different frequencies. Lakes may provide ecological opportunity via ecological release (e.g., from predators/competitors) but need not be coupled with access to novel ecological niches.

Derived faunivores are the forerunners of major synapsid radiations

Evolutionary radiations generate most of Earth's biodiversity, but are there common ecomorphological traits among the progenitors of radiations? In Synapsida (the mammalian total group), 'small-bodied faunivore' has been hypothesized as the ancestral state of most major radiating clades, but this has not been quantitatively assessed across multiple radiations. To examine macroevolutionary patterns in a phylogenetic context, we generated a time-calibrated metaphylogeny ('metatree') comprising 1,888 synapsid species from the Carboniferous through the Eocene (305-34 Ma) based on 269 published character matrices. We used comparative methods to investigate body size and dietary evolution during successive synapsid radiations. Faunivory is the ancestral dietary regime of each major synapsid radiation, but relatively small body size is only established as the common ancestral state of radiations near the origin of Mammaliaformes in the Late Triassic. The faunivorous ancestors of synapsid radiations typically have numerous novel characters compared with their contemporaries, and these derived traits may have helped them to survive faunal turnover events and subsequently radiate.

Fitting and evaluating univariate and multivariate models of within-lineage evolution

The nature of phenotypic evolution within lineages is central to many unresolved questions in paleontology and evolutionary biology. Analyses of evolutionary time-series of ancestor-descendant populations in the fossil record are likely to make important contributions to many of these debates. However, the limited number of models that have been applied to these types of data may restrict our ability to interpret phenotypic evolution in the fossil record. Using uni- and multivariate models of trait evolution that make different assumptions regarding the dynamics of the adaptive landscape, I evaluate contrasting hypotheses to explain evolution of size in the radiolarian Eucyrtidium calvertense and armor in the stickleback Gaserosteus doryssus. Body size evolution in E. calvertense is best explained by a model where the lineage evolves as a consequence of a shift in the adaptive landscape that coincides with the initiation of neosympatry with its sister lineage. Multivariate evolution of armor traits in a stickleback lineage (Gasterosteus doryssus) show evidence of adaptation towards independent optima on the adaptive landscape at the same time as traits change in a correlated fashion. The fitted models are available in a the R package evoTS, which builds on the commonly used paleoTS framework.

Remarkable variation in androecial morphology is closely associated with corolla traits in Western Hemisphere Justiciinae (Acanthaceae: Justicieae)

BACKGROUND AND AIMS

Few studies of angiosperms have focused on androecial evolution in conjunction with evolutionary shifts in corolla morphology and pollinator relationships. The Western Hemisphere clade of Justiciinae (Acanthaceae) presents the rare opportunity to examine remarkable diversity in staminal morphology. We took a phylogenetically informed approach to examine staminal diversity in this hypervariable group and asked whether differences in anther thecae separation is associated with phylogenetically informed patterns of variation in corolla morphology. We further discuss evidence for associations between anther diversity and pollinators in this lineage.

METHODS

For the Dianthera/Sarotheca/Plagiacanthus (DSP) clade of Western Hemisphere Justiciinae, we characterized floral diversity based on a series of corolla measurements and using a model-based clustering approach. We then tested for correlations between anther thecae separation and corolla traits, and for shifts in trait evolution, including evidence for convergence.

KEY RESULTS

There is evolutionary vagility in corolla and anther traits across the DSP clade with little signal of phylogenetic constraint. Floral morphology clusters into four distinct groups that are, in turn, strongly associated with anther thecae separation, a novel result in Acanthaceae and, to our knowledge, across flowering plants. These cluster groups are marked by floral traits that strongly point to associations with pollinating animals. Specifically, species that are known or likely to be hummingbird pollinated have stamens with parallel thecae, whereas those that are likely bee or fly pollinated have stamens with offset, divergent thecae.

CONCLUSIONS

Our results suggest that anther thecae separation is likely under selection in concert with other corolla characters. Significant morphological shifts detected by our analyses corresponded to putative shifts from insect to hummingbird pollination. Results from this study support the hypothesis that floral structures function in an integrated manner and are likely subject to selection as a suite. Further, these changes can be hypothesized to represent adaptive evolution.

Diversity and evolution of cerebellar folding in mammals

The process of brain folding is thought to play an important role in the development and organisation of the cerebrum and the cerebellum. The study of cerebellar folding is challenging due to the small size and abundance of its folia. In consequence, little is known about its anatomical diversity and evolution. We constituted an open collection of histological data from 56 mammalian species and manually segmented the cerebrum and the cerebellum. We developed methods to measure the geometry of cerebellar folia and to estimate the thickness of the molecular layer. We used phylogenetic comparative methods to study the diversity and evolution of cerebellar folding and its relationship with the anatomy of the cerebrum. Our results show that the evolution of cerebellar and cerebral anatomy follows a stabilising selection process. We observed two groups of phenotypes changing concertedly through evolution: a group of 'diverse' phenotypes - varying over several orders of magnitude together with body size, and a group of 'stable' phenotypes varying over less than 1 order of magnitude across species. Our analyses confirmed the strong correlation between cerebral and cerebellar volumes across species, and showed in addition that large cerebella are disproportionately more folded than smaller ones. Compared with the extreme variations in cerebellar surface area, folial anatomy and molecular layer thickness varied only slightly, showing a much smaller increase in the larger cerebella. We discuss how these findings could provide new insights into the diversity and evolution of cerebellar folding, the mechanisms of cerebellar and cerebral folding, and their potential influence on the organisation of the brain across species.

A Cautionary Note on "A Cautionary Note on the Use of Ornstein Uhlenbeck Models in Macroevolutionary Studies"

Models based on the Ornstein-Uhlenbeck process have become standard for the comparative study of adaptation. Cooper et al. (2016) have cast doubt on this practice by claiming statistical problems with fitting Ornstein-Uhlenbeck models to comparative data. Specifically, they claim that statistical tests of Brownian motion may have too high Type I error rates and that such error rates are exacerbated by measurement error. In this note, we argue that these results have little relevance to the estimation of adaptation with Ornstein-Uhlenbeck models for three reasons. First, we point out that Cooper et al. (2016) did not consider the detection of distinct optima (e.g. for different environments), and therefore did not evaluate the standard test for adaptation. Second, we show that consideration of parameter estimates, and not just statistical significance, will usually lead to correct inferences about evolutionary dynamics. Third, we show that bias due to measurement error can be corrected for by standard methods. We conclude that Cooper et al. (2016) have not identified any statistical problems specific to Ornstein-Uhlenbeck models, and that their cautions against their use in comparative analyses are unfounded and misleading. [adaptation, Ornstein-Uhlenbeck model, phylogenetic comparative method.].

Developmental origin underlies evolutionary rate variation across the placental skull

The placental skull has evolved into myriad forms, from longirostrine whales to globular primates, and with a diverse array of appendages from antlers to tusks. This disparity has recently been studied from the perspective of the whole skull, but the skull is composed of numerous elements that have distinct developmental origins and varied functions. Here, we assess the evolution of the skull's major skeletal elements, decomposed into 17 individual regions. Using a high-dimensional morphometric approach for a dataset of 322 living and extinct eutherians (placental mammals and their stem relatives), we quantify patterns of variation and estimate phylogenetic, allometric and ecological signal across the skull. We further compare rates of evolution across ecological categories and ordinal-level clades and reconstruct rates of evolution along lineages and through time to assess whether developmental origin or function discriminate the evolutionary trajectories of individual cranial elements. Our results demonstrate distinct macroevolutionary patterns across cranial elements that reflect the ecological adaptations of major clades. Elements derived from neural crest show the fastest rates of evolution, but ecological signal is equally pronounced in bones derived from neural crest and paraxial mesoderm, suggesting that developmental origin may influence evolutionary tempo, but not capacity for specialisation. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Multiple modes of inference reveal less phylogenetic signal in marsupial basicranial shape compared with the rest of the cranium

Incorporating morphological data into modern phylogenies allows integration of fossil evidence, facilitating divergence dating and macroevolutionary inferences. Improvements in the phylogenetic utility of morphological data have been sought via Procrustes-based geometric morphometrics (GMM), but with mixed success and little clarity over what anatomical areas are most suitable. Here, we assess GMM-based phylogenetic reconstructions in a heavily sampled source of discrete characters for mammalian phylogenetics-the basicranium-in 57 species of marsupial mammals, compared with the remainder of the cranium. We show less phylogenetic signal in the basicranium compared with a 'Rest of Cranium' partition, using diverse metrics of phylogenetic signal (K_mult, phylogenetically aligned principal components analysis, comparisons of UPGMA/neighbour-joining/parsimony trees and cophenetic distances to a reference phylogeny) for scaled, Procrustes-aligned landmarks and allometry-corrected residuals. Surprisingly, a similar pattern emerged from parsimony-based analyses of discrete cranial characters. The consistent results across methods suggest that easily computed metrics such as K_mult can provide good guidance on phylogenetic information in a landmarking configuration. In addition, GMM data may be less informative for intricate but conservative anatomical regions such as the basicranium, while better-but not necessarily novel-phylogenetic information can be expected for broadly characterized shapes such as entire bones. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Bridging Performance and Adaptive Landscapes to Understand Long-Term Functional Evolution

AbstractUnderstanding functional adaptation demands an integrative framework that captures the complex interactions between form, function, ecology, and evolutionary processes. In this review, we discuss how to integrate the following two distinct approaches to better understand functional evolution: (1) the adaptive landscape approach (ALA), aimed at finding adaptive peaks for different ecologies, and (2) the performance landscape approach (PLA), aimed at finding performance peaks for different ecologies. We focus on the Ornstein-Uhlenbeck process as the evolutionary model for the ALA and on biomechanical modeling to estimate performance for the PLA. Whereas both the ALA and the PLA have each given insight into functional adaptation, separately they cannot address how much performance contributes to fitness or whether evolutionary constraints have played a role in form-function evolution. We show that merging these approaches leads to a deeper understanding of these issues. By comparing the locations of performance and adaptive peaks, we can infer how much performance contributes to fitness in species' current environments. By testing for the relevance of history on phenotypic variation, we can infer the influence of past selection and constraints on functional adaptation. We apply this merged framework in a case study of turtle shell evolution and explain how to interpret different possible outcomes. Even though such outcomes can be quite complex, they represent the multifaceted relations among function, fitness, and constraints.

Ecological and life-history drivers of avian skull evolution

One of the most famous examples of adaptive radiation is that of the Galápagos finches, where skull morphology, particularly the beak, varies with feeding ecology. Yet increasingly studies are questioning the strength of this correlation between feeding ecology and morphology in relation to the entire neornithine radiation, suggesting that other factors also significantly affect skull evolution. Here, we broaden this debate to assess the influence of a range of ecological and life-history factors, specifically habitat density, migration, and developmental mode, in shaping avian skull evolution. Using 3D geometric morphometric data to robustly quantify skull shape for 354 extant species spanning avian diversity, we fitted flexible phylogenetic regressions and estimated evolutionary rates for each of these factors across the full data set. The results support a highly significant relationship between skull shape and both habitat density and migration, but not developmental mode. We further found heterogenous rates of evolution between different character states within habitat density, migration, and developmental mode, with rapid skull evolution in species that occupy dense habitats, are migratory, or are precocial. These patterns demonstrate that diverse factors affect the tempo and mode of avian phenotypic evolution and that skull evolution in birds is not simply a reflection of feeding ecology.

Early euprimates already had a diverse locomotor repertoire: Evidence from ankle bone morphology

The morphological adaptations of euprimates have been linked to their origin and early evolution in an arboreal environment. However, the ancestral and early locomotor repertoire of this group remains contentious. Although some tarsal bones like the astragalus and the calcaneus have been thoroughly studied, the navicular remains poorly studied despite its potential implications for foot mobility. Here, we evaluate early euprimate locomotion by assessing the shape of the navicular-an important component of the midtarsal region of the foot-using three-dimensional geometric morphometrics in relation to quantified locomotor repertoire in a wide data set of extant primates. We also reconstruct the locomotor repertoire of representatives of the major early primate lineages with a novel phylogenetically informed discriminant analysis and characterize the changes that occurred in the navicular during the archaic primate-euprimate transition. To do so, we included in our study an extensive sample of naviculars (36 specimens) belonging to different species of adapiforms, omomyiforms, and plesiadapiforms. Our results indicate that navicular shape embeds a strong functional signal, allowing us to infer the type of locomotion of extinct primates. We demonstrate that early euprimates displayed a diverse locomotor behavior, although they did not reach the level of specialization of some living forms. Finally, we show that the navicular bone experienced substantial reorganization throughout the archaic primate-euprimate transition, supporting the major functional role of the tarsus during early primate evolution. This study demonstrates that navicular shape can be used as a reliable proxy for primate locomotor behavior. In addition, it sheds light on the diverse locomotor behavior of early primates as well as on the archaic primate-euprimate transition, which involved profound morphological changes within the tarsus, including the navicular bone.

Protein length distribution is remarkably uniform across the tree of life

BACKGROUND

In every living species, the function of a protein depends on its organization of structural domains, and the length of a protein is a direct reflection of this. Because every species evolved under different evolutionary pressures, the protein length distribution, much like other genomic features, is expected to vary across species but has so far been scarcely studied.

RESULTS

Here we evaluate this diversity by comparing protein length distribution across 2326 species (1688 bacteria, 153 archaea, and 485 eukaryotes). We find that proteins tend to be on average slightly longer in eukaryotes than in bacteria or archaea, but that the variation of length distribution across species is low, especially compared to the variation of other genomic features (genome size, number of proteins, gene length, GC content, isoelectric points of proteins). Moreover, most cases of atypical protein length distribution appear to be due to artifactual gene annotation, suggesting the actual variation of protein length distribution across species is even smaller.

CONCLUSIONS

These results open the way for developing a genome annotation quality metric based on protein length distribution to complement conventional quality measures. Overall, our findings show that protein length distribution between living species is more uniform than previously thought. Furthermore, we also provide evidence for a universal selection on protein length, yet its mechanism and fitness effect remain intriguing open questions.

Evolutionary history explains foliar spectral differences between arbuscular and ectomycorrhizal plant species

Leaf spectra are integrated foliar phenotypes that capture a range of traits and can provide insight into ecological processes. Leaf traits, and therefore leaf spectra, may reflect belowground processes such as mycorrhizal associations. However, evidence for the relationship between leaf traits and mycorrhizal association is mixed, and few studies account for shared evolutionary history. We conduct partial least squares discriminant analysis to assess the ability of spectra to predict mycorrhizal type. We model the evolution of leaf spectra for 92 vascular plant species and use phylogenetic comparative methods to assess differences in spectral properties between arbuscular mycorrhizal and ectomycorrhizal plant species. Partial least squares discriminant analysis classified spectra by mycorrhizal type with 90% (arbuscular) and 85% (ectomycorrhizal) accuracy. Univariate models of principal components identified multiple spectral optima corresponding with mycorrhizal type due to the close relationship between mycorrhizal type and phylogeny. Importantly, we found that spectra of arbuscular mycorrhizal and ectomycorrhizal species do not statistically differ from each other after accounting for phylogeny. While mycorrhizal type can be predicted from spectra, enabling the use of spectra to identify belowground traits using remote sensing, this is due to evolutionary history and not because of fundamental differences in leaf spectra due to mycorrhizal type.

The evolution of fast-growing coral reef fishes

Individual growth is a fundamental life history trait^1-4, yet its macroevolutionary trajectories have rarely been investigated for entire animal assemblages. Here we analyse the evolution of growth in a highly diverse vertebrate assemblage-coral reef fishes. We combine state-of-the-art extreme gradient boosted regression trees with phylogenetic comparative methods to detect the timing, number, location and magnitude of shifts in the adaptive regime of somatic growth. We also explored the evolution of the allometric relationship between body size and growth. Our results show that the evolution of fast growth trajectories in reef fishes has been considerably more common than the evolution of slow growth trajectories. Many reef fish lineages shifted towards faster growth and smaller body size evolutionary optima in the Eocene (56-33.9 million years ago), pointing to a major expansion of life history strategies in this Epoch. Of all lineages examined, the small-bodied, high-turnover cryptobenthic fishes shifted most towards extremely high growth optima, even after accounting for body size allometry. These results suggest that the high global temperatures of the Eocene⁵ and subsequent habitat reconfigurations⁶ might have been critical for the rise and retention of the highly productive, high-turnover fish faunas that characterize modern coral reef ecosystems.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

Size And Locomotor Ecology Have Differing Effects on the External and Internal Morphologies of Squirrel (Rodentia: Sciuridae) Limb Bones

Mammals exhibit a diverse range of limb morphologies that are associated with different locomotor ecologies and structural mechanics. Much remains to be investigated, however, about the combined effects of locomotor modes and scaling on the external shape and structural properties of limb bones. Here, we used squirrels (Sciuridae) as a model clade to examine the effects of locomotor mode and scaling on the external shape and structure of the two major limb bones, the humerus and femur. We quantified humeral and femoral morphologies using 3D geometric morphometrics and bone structure analyses on a sample of 76 squirrel species across their four major ecotypes. We then used phylogenetic generalized linear models to test how locomotor ecology, size, and their interaction influenced morphological traits. We found that size and locomotor mode exhibit different relationships with the external shape and structure of the limb bones, and that these relationships differ between the humerus and femur. External shapes of the humerus and, to a lesser extent, the femur are best explained by locomotor ecology rather than by size, whereas structures of both bones are best explained by interactions between locomotor ecology and scaling. Interestingly, the statistical relationships between limb morphologies and ecotype were lost when accounting for phylogenetic relationships among species under Brownian motion. That assuming Brownian motion confounded these relationships is not surprising considering squirrel ecotypes are phylogenetically clustered; our results suggest that humeral and femoral variation partitioned early between clades and their ecomorphologies were maintained to the present. Overall, our results show how mechanical constraints, locomotor ecology, and evolutionary history may enact different pressures on the shape and structure of limb bones in mammals.

Patterns of ontogenetic evolution across extant marsupials reflect different allometric pathways to ecomorphological diversity

The relatively high level of morphological diversity in Australasian marsupials compared to that observed among American marsupials remains poorly understood. We undertake a comprehensive macroevolutionary analysis of ontogenetic allometry of American and Australasian marsupials to examine whether the contrasting levels of morphological diversity in these groups are reflected in their patterns of allometric evolution. We collate ontogenetic series for 62 species and 18 families of marsupials (n = 2091 specimens), spanning across extant marsupial diversity. Our results demonstrate significant lability of ontogenetic allometric trajectories among American and Australasian marsupials, yet a phylogenetically structured pattern of allometric evolution is preserved. Here we show that species diverging more than 65 million years ago converge in their patterns of ontogenetic allometry under animalivorous and herbivorous diets, and that Australasian marsupials do not show significantly greater variation in patterns of ontogenetic allometry than their American counterparts, despite displaying greater magnitudes of extant ecomorphological diversity.

Shark mandible evolution reveals patterns of trophic and habitat-mediated diversification

Environmental controls of species diversity represent a central research focus in evolutionary biology. In the marine realm, sharks are widely distributed, occupying mainly higher trophic levels and varied dietary preferences, mirrored by several morphological traits and behaviours. Recent comparative phylogenetic studies revealed that sharks present a fairly uneven diversification across habitats, from reefs to deep-water. We show preliminary evidence that morphological diversification (disparity) in the feeding system (mandibles) follows these patterns, and we tested hypotheses linking these patterns to morphological specialisation. We conducted a 3D geometric morphometric analysis and phylogenetic comparative methods on 145 specimens representing 90 extant shark species using computed tomography models. We explored how rates of morphological evolution in the jaw correlate with habitat, size, diet, trophic level, and taxonomic order. Our findings show a relationship between disparity and environment, with higher rates of morphological evolution in reef and deep-water habitats. Deep-water species display highly divergent morphologies compared to other sharks. Strikingly, evolutionary rates of jaw disparity are associated with diversification in deep water, but not in reefs. The environmental heterogeneity of the offshore water column exposes the importance of this parameter as a driver of diversification at least in the early part of clade history.

Evolution of cortical geometry and its link to function, behaviour and ecology

Studies in comparative neuroanatomy and of the fossil record demonstrate the influence of socio-ecological niches on the morphology of the cerebral cortex, but have led to oftentimes conflicting theories about its evolution. Here, we study the relationship between the shape of the cerebral cortex and the topography of its function. We establish a joint geometric representation of the cerebral cortices of ninety species of extant Euarchontoglires, including commonly used experimental model organisms. We show that variability in surface geometry relates to species' ecology and behaviour, independent of overall brain size. Notably, ancestral shape reconstruction of the cortical surface and its change during evolution enables us to trace the evolutionary history of localised cortical expansions, modal segregation of brain function, and their association to behaviour and cognition. We find that individual cortical regions follow different sequences of area increase during evolutionary adaptations to dynamic socio-ecological niches. Anatomical correlates of this sequence of events are still observable in extant species, and relate to their current behaviour and ecology. We decompose the deep evolutionary history of the shape of the human cortical surface into spatially and temporally conscribed components with highly interpretable functional associations, highlighting the importance of considering the evolutionary history of cortical regions when studying their anatomy and function.

Diversification of the ruminant skull along an evolutionary line of least resistance

Clarifying how microevolutionary processes scale to macroevolutionary patterns is a fundamental goal in evolutionary biology, but these analyses, requiring comparative datasets of population-level variation, are limited. By analyzing a previously published dataset of 2859 ruminant crania, we find that variation within and between ruminant species is biased by a highly conserved mammalian-wide allometric pattern, CREA (craniofacial evolutionary allometry), where larger species have proportionally longer faces. Species with higher morphological integration and species more biased toward CREA have diverged farther from their ancestors, and Ruminantia as a clade diversified farther than expected in the direction of CREA. Our analyses indicate that CREA acts as an evolutionary "line of least resistance" and facilitates morphological diversification due to its alignment with the browser-grazer continuum. Together, our results demonstrate that constraints at the population level can produce highly directional patterns of phenotypic evolution at the macroevolutionary scale. Further research is needed to explore how CREA has been exploited in other mammalian clades.

Evolutionary convergence of muscle architecture in relation to locomotor ecology in snakes

The epaxial muscles in snakes are responsible for locomotion and as such can be expected to show adaptations in species living in different environments. Here, we tested whether the structural units that comprise the superficial epaxial muscles (semispinalis-spinalis, SSP; longissimus dorsi, LD; iliocostalis, IC) were different in animals occupying similar habitats. To do so, we analyzed and compared the muscle architecture (mass, fiber length, and physiological cross-sectional area) of the superficial epaxial muscle segments in snakes that differ in their habitat use (e.g., arboreal, terrestrial, and aquatic). Our results showed that arboreal species have on average longer muscles and tendons spanning more segments likely important during gap bridging. Moreover, aquatic snakes show relatively heavier semispinalis-spinalis muscles with a greater cross-sectional area. The longissimus dorsi muscles also showed a greater cross-sectional area compared with terrestrial and especially arboreal snakes. Whereas the more strongly developed muscles in aquatic snakes are likely associated with the dense and viscous environment through which they move, the lighter muscles in arboreal snakes may provide an advantage when climbing. Future studies comparing other ecologies (e.g., burrowing snakes) and additional muscle units (e.g., multifidus; hypaxial muscles) are needed to better understand the structural features driving variation in locomotor performance and efficiency in snakes.

Mammalian forelimb evolution is driven by uneven proximal-to-distal morphological diversity

Vertebrate limb morphology often reflects the environment due to variation in locomotor requirements. However, proximal and distal limb segments may evolve differently from one another, reflecting an anatomical gradient of functional specialization that has been suggested to be impacted by the timing of development. Here, we explore whether the temporal sequence of bone condensation predicts variation in the capacity of evolution to generate morphological diversity in proximal and distal forelimb segments across more than 600 species of mammals. Distal elements not only exhibit greater shape diversity, but also show stronger within-element integration and, on average, faster evolutionary responses than intermediate and upper limb segments. Results are consistent with the hypothesis that late developing distal bones display greater morphological variation than more proximal limb elements. However, the higher integration observed within the autopod deviates from such developmental predictions, suggesting that functional specialization plays an important role in driving within-element covariation. Proximal and distal limb segments also show different macroevolutionary patterns, albeit not showing a perfect proximo-distal gradient. The high disparity of the mammalian autopod, reported here, is consistent with the higher potential of development to generate variation in more distal limb structures, as well as functional specialization of the distal elements.

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.

Fertilization mode differentially impacts the evolution of vertebrate sperm components

Environmental change frequently drives morphological diversification, including at the cellular level. Transitions in the environment where fertilization occurs (i.e., fertilization mode) are hypothesized to be a driver of the extreme diversity in sperm morphology observed in animals. Yet how fertilization mode impacts the evolution of sperm components-head, midpiece, and flagellum-each with different functional roles that must act as an integrated unit remains unclear. Here, we test this hypothesis by examining the evolution of sperm component lengths across 1103 species of vertebrates varying in fertilization mode (external vs. internal fertilization). Sperm component length is explained in part by fertilization mode across vertebrates, but how fertilization mode influences sperm evolution varies among sperm components and vertebrate clades. We also identify evolutionary responses not influenced by fertilization mode: midpieces evolve rapidly in both external and internal fertilizers. Fertilization mode thus influences vertebrate sperm evolution through complex component- and clade-specific evolutionary responses.