EARLY BURSTS OF BODY SIZE AND SHAPE EVOLUTION ARE RARE IN COMPARATIVE DATA

phyddle: software for phylogenetic model exploration with deep learning

Many realistic phylogenetic models lack tractable likelihood functions, prohibiting their use with standard inference methods. We present phyddle, a pipeline-based toolkit for performing phylogenetic modeling tasks using likelihood-free deep learning approaches. phyddle coordinates modeling tasks through five analysis steps (Simulate, Format, Train, Estimate, and Plot) that transform raw phylogenetic datasets as input into numerical and visualized model-based output. Benchmarks show that phyddle accurately performs a range of inference tasks, such as estimating macroevolutionary parameters, selecting among continuous trait evolution models, and passing coverage tests for epidemiological models, even for models that lack tractable likelihoods. phyddle has a flexible command-line interface, making it easy to integrate deep learning approaches for phylogenetics into research workflows. Learn more about phyddle at https://phyddle.org.

Navigating uncertainty in maximum body size in marine metazoans

Body size is a fundamental biological trait shaping ecological interactions, evolutionary processes, and our understanding of the structure and dynamics of marine communities on a global scale. Accurately defining a species' body size, despite the ease of measurement, poses significant challenges due to varied methodologies, tool usage, and subjectivity among researchers, resulting in multiple, often discrepant size estimates. These discrepancies, stemming from diverse measurement approaches and inherent variability, could substantially impact the reliability and precision of ecological and evolutionary studies reliant on body size data across extensive species datasets. This study examines the variation in reported maximum body sizes across 69,570 individual measurements of maximum size, ranging from <0.2 μm to >45 m, for 27,271 species of marine metazoans. The research aims to investigate how reported maximum size variations within species relate to organism size, taxonomy, habitat, and the presence of skeletal structures. The investigation particularly focuses on understanding why discrepancies in maximum size estimates arise and their potential implications for broader ecological and evolutionary studies relying on body size data. Variation in reported maximum sizes is zero for 38% of species, and low for most species, although it exceeds two orders of magnitude for some species. The likelihood of zero variation in maximum size decreased with more measurements and increased in larger species, though this varied across phyla and habitats. Pelagic organisms consistently had low maximum size range values, while small species with unspecified habitats had the highest variation. Variations in maximum size within a species were notably smaller than interspecific variation at higher taxonomic levels. Significant variation in maximum size estimates exists within marine species, and partially explained by organism size, taxonomic group, and habitat. Variation in maximum size could be reduced by standardized measurement protocols and improved meta-data. Despite the variation, egregious errors in published maximum size measurements are rare, and their impact on comparative macroecological and macroevolutionary research is likely minimal.

A new abelisaurid dinosaur from the end Cretaceous of Patagonia and evolutionary rates among the Ceratosauria

Gondwanan dinosaur faunae during the 20 Myr preceding the Cretaceous-Palaeogene (K/Pg) extinction included several lineages that were absent or poorly represented in Laurasian landmasses. Among these, the South American fossil record contains diverse abelisaurids, arguably the most successful groups of carnivorous dinosaurs from Gondwana in the Cretaceous, reaching their highest diversity towards the end of this period. Here we describe Koleken inakayali gen. et sp. n., a new abelisaurid from the La Colonia Formation (Maastrichtian, Upper Cretaceous) of Patagonia. Koleken inakayali is known from several skull bones, an almost complete dorsal series, complete sacrum, several caudal vertebrae, pelvic girdle and almost complete hind limbs. The new abelisaurid shows a unique set of features in the skull and several anatomical differences from Carnotaurus sastrei (the only other abelisaurid known from the La Colonia Formation). Koleken inakayali is retrieved as a brachyrostran abelisaurid, clustered with other South American abelisaurids from the latest Cretaceous (Campanian-Maastrichtian), such as Aucasaurus, Niebla and Carnotaurus. Leveraging our phylogeny estimates, we explore rates of morphological evolution across ceratosaurian lineages, finding them to be particularly high for elaphrosaurine noasaurids and around the base of Abelisauridae, before the Early Cretaceous radiation of the latter clade. The Noasauridae and their sister clade show contrasting patterns of morphological evolution, with noasaurids undergoing an early phase of accelerated evolution of the axial and hind limb skeleton in the Jurassic, and the abelisaurids exhibiting sustained high rates of cranial evolution during the Early Cretaceous. These results provide much needed context for the evolutionary dynamics of ceratosaurian theropods, contributing to broader understanding of macroevolutionary patterns across dinosaurs.

Climatic and biogeographic processes underlying the diversification of the pantropical flowering plant family Annonaceae

Tropical forests harbor the richest biodiversity among terrestrial ecosystems, but few studies have addressed the underlying processes of species diversification in these ecosystems. We use the pantropical flowering plant family Annonaceae as a study system to investigate how climate and biogeographic events contribute to diversification. A super-matrix phylogeny comprising 835 taxa (34% of Annonaceae species) based on eight chloroplast regions was used in this study. We show that global temperature may better explain the recent rapid diversification in Annonaceae than time and constant models. Accelerated accumulation of niche divergence (around 15 Ma) lags behind the increase of diversification rate (around 25 Ma), reflecting a heterogeneous transition to recent diversity increases. Biogeographic events are related to only two of the five diversification rate shifts detected. Shifts in niche evolution nevertheless appear to be associated with increasingly seasonal environments. Our results do not support the direct correlation of any particular climatic niche shifts or historical biogeographical event with shifts in diversification rate. Instead, we suggest that Annonaceae diversification can lead to later niche divergence as a result of increasing interspecific competition arising from species accumulation. Shifts in niche evolution appear to be associated with increasingly seasonal environments. Our results highlight the complexity of diversification in taxa with long evolutionary histories.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Using phylogeographic link-prediction in primates to prioritize human parasite screening

Objectives

The ongoing risk of emerging infectious disease has renewed calls for understanding the origins of zoonoses and identifying future zoonotic disease threats. Given their close phylogenetic relatedness and geographic overlap with humans, non-human primates (NHPs) have been the source of many infectious diseases throughout human evolution. NHPs harbor diverse parasites, with some infecting only a single host species while others infect species from multiple families.

Materials and Methods

We applied a novel link-prediction method to predict undocumented instances of parasite sharing between humans and NHPs. Our model makes predictions based on phylogenetic distances and geographic overlap among NHPs and humans in six countries with high NHP diversity: Columbia, Brazil, Democratic Republic of Congo, Madagascar, China and Indonesia.

Results

Of the 899 human parasites documented in the Global Infectious Diseases and Epidemiology Network (GIDEON) database for these countries, 12% were shared with at least one other NHP species. The link prediction model identified an additional 54 parasites that are likely to infect humans but were not reported in GIDEON. These parasites were mostly host generalists, yet their phylogenetic host breadth varied substantially.

Discussion

As human activities and populations encroach on NHP habitats, opportunities for parasite sharing between human and non-human primates will continue to increase. Our study identifies specific infectious organisms to monitor in countries with high NHP diversity, while the comparative analysis of host generalism, parasite taxonomy, and transmission mode provides insights to types of parasites that represent high zoonotic risk.

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.

Disparification and extinction trade-offs shaped the evolution of Permian to Jurassic Odonata

Owing to their prevalence in nowadays terrestrial ecosystems, insects are a relevant group to assess the impact of mass extinctions on emerged land. However, limitations of the insect fossil record make it difficult to assess the impact of such events based on taxonomic diversity alone. Therefore, we documented trends in morphological diversity, i.e., disparity, using wings of Permian to Jurassic Odonata as model. Our results show a decreasing trend in disparity while species richness increased. Both the Permian-Triassic and Triassic-Jurassic transitions are revealed as important events, associated with strong morphospace restructuring due to selective extinction. In each case, a recovery was assured by the diversification of new forms compensating the loss of others. Early representatives of Odonata continuously evolved new shapes, a pattern contrasting with the classical assertion of a morphospace fulfilled early and followed by selective extinctions and specialization within it.

Evaluating the Performance of Widely Used Phylogenetic Models for Gene Expression Evolution

Phylogenetic comparative methods are increasingly used to test hypotheses about the evolutionary processes that drive divergence in gene expression among species. However, it is unknown whether the distributional assumptions of phylogenetic models designed for quantitative phenotypic traits are realistic for expression data and importantly, the reliability of conclusions of phylogenetic comparative studies of gene expression may depend on whether the data is well-described by the chosen model. To evaluate this, we first fit several phylogenetic models of trait evolution to 8 previously published comparative expression datasets, comprising a total of 54,774 genes with 145,927 unique gene-tissue combinations. Using a previously developed approach, we then assessed how well the best model of the set described the data in an absolute (not just relative) sense. First, we find that Ornstein-Uhlenbeck models, in which expression values are constrained around an optimum, were the preferred model for 66% of gene-tissue combinations. Second, we find that for 61% of gene-tissue combinations, the best fit model of the set was found to perform well; the rest were found to be performing poorly by at least one of the test statistics we examined. Third, we find that when simple models do not perform well, this appears to be typically a consequence of failing to fully account for heterogeneity in the rate of the evolution. We advocate that assessment of model performance should become a routine component of phylogenetic comparative expression studies; doing so can improve the reliability of inferences and inspire the development of novel models.

Body Size Poorly Predicts Host-Associated Microbial Diversity in Wild Birds

The composition and diversity of avian microbiota are shaped by many factors, including host ecologies and environmental variables. In this study, we examine microbial diversity across 214 bird species sampled in Malawi at five major body sites: blood, buccal cavity, gizzard, intestinal tract, and cloaca. Microbial community dissimilarity differed significantly across body sites. Ecological theory predicts that as area increases, so does diversity. We tested the hypothesis that avian microbiota diversity is correlated with body size, used as a proxy for area, using comparative phylogenetic methods. Using Pagel's lambda, we found that few microbial diversity metrics had significant phylogenetic signals. Phylogenetic generalized least squares identified a significant but weak negative correlation between host size and microbial diversity of the blood and a similarly significant but weakly positive correlation between the cloacal microbiota and host size among birds within the order Passeriformes. Phylosymbiosis, or a congruent branching pattern between host phylogeny and their associated microbiota similarity, was tested and found to be weak or not significant in four of the body sites with sufficient sample size (blood, buccal, cloaca, and intestines). Taken together, these results suggest that the avian microbiome is highly variable, with microbiota diversity demonstrating few clear associations with bird size. Finally, the blood microbiota have a unique relationship with host size. IMPORTANCE All animals coexist and interact with microorganisms, including bacteria, archaea, microscopic eukaryotes, and viruses. These microorganisms can have an enormous influence on the biology and health of macro-organisms. However, the general rules that govern these host-associated microbial communities are poorly described, especially in wild animals. In this paper, we investigate the microbial communities of over 200 species of birds from Malawi and characterize five body site bacterial microbiota in depth. Because the evolutionary relationships of the host underlie the relationship between any host-associated microbiota relationships, we use phylogenetic comparative methods to account for this relationship. We find that the size of a host (the bird) and the diversity and composition of the microbiota are largely uncorrelated. We also find that the general pattern of similarity between host phylogeny and microbiota similarity is weak. Together, we see that bird microbiota are not strongly tied to host size or evolutionary history.

Adaptations to a cold climate promoted social evolution in Asian colobine primates

The biological mechanisms that underpin primate social evolution remain poorly understood. Asian colobines display a range of social organizations, which makes them good models for investigating social evolution. By integrating ecological, geological, fossil, behavioral, and genomic analyses, we found that colobine primates that inhabit colder environments tend to live in larger, more complex groups. Specifically, glacial periods during the past 6 million years promoted the selection of genes involved in cold-related energy metabolism and neurohormonal regulation. More-efficient dopamine and oxytocin pathways developed in odd-nosed monkeys, which may have favored the prolongation of maternal care and lactation, increasing infant survival in cold environments. These adaptive changes appear to have strengthened interindividual affiliation, increased male-male tolerance, and facilitated the stepwise aggregation from independent one-male groups to large multilevel societies.

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.

Testing drivers of acoustic divergence in cicadas (Cicadidae: Tettigettalna)

Divergence in acoustic signals may have a crucial role in the speciation process of animals that rely on sound for intra-specific recognition and mate attraction. The acoustic adaptation hypothesis (AAH) postulates that signals should diverge according to the physical properties of the signalling environment. To be efficient, signals should maximize transmission and decrease degradation. To test which drivers of divergence exert the most influence in a speciose group of insects, we used a phylogenetic approach to the evolution of acoustic signals in the cicada genus Tettigettalna, investigating the relationship between acoustic traits (and their mode of evolution) and body size, climate and micro-/macro-habitat usage. Different traits showed different evolutionary paths. While acoustic divergence was generally independent of phylogenetic history, some temporal variables' divergence was associated with genetic drift. We found support for ecological adaptation at the temporal but not the spectral level. Temporal patterns are correlated with micro- and macro-habitat usage and temperature stochasticity in ways that run against the AAH predictions, degrading signals more easily. These traits are likely to have evolved as an anti-predator strategy in conspicuous environments and low-density populations. Our results support a role of ecological selection, not excluding a likely role of sexual selection in the evolution of Tettigettalna calling songs, which should be further investigated in an integrative approach.

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.

Evidence for an evo-devo-derived hypothesis on three-dimensional flower shape modularity in a tropical orchid clade

Covarying suites of phenotypic traits, or modules, are increasingly recognized to promote morphological evolution. However, information on how modularity influences flower diversity is rare and lacking for Orchidaceae. Here, we combine high-resolution X-ray computed tomography scanning with three-dimensional geometric morphometrics and phylogenetic comparative methods to test various hypotheses about three-dimensional patterns of flower evolutionary modularity in Malagasy Bulbophyllum orchids and examine rates and modes of module evolution. Based on the four evolutionary modules identified (i.e., sepals, lateral petals, labellum + column-foot, and column-part), our data support the hypothesis that both genetic-developmental and functional adaptive factors shaped evolutionary flower trait covariation in these tropical orchids. In line with "evo-devo" studies, we also find that the labellum evolved independently from the rest of the petal whorl. Finally, we show that modules evolved with different rates, and either in a neutral fashion (only column-part) or under selective constraints, as likely imposed by pollinators. Overall, this study supports current views that modular units can enhance the range and rate of morphological evolution.

Attenuated evolution of mammals through the Cenozoic

The Cenozoic diversification of placental mammals is the archetypal adaptive radiation. Yet, discrepancies between molecular divergence estimates and the fossil record fuel ongoing debate around the timing, tempo, and drivers of this radiation. Analysis of a three-dimensional skull dataset for living and extinct placental mammals demonstrates that evolutionary rates peak early and attenuate quickly. This long-term decline in tempo is punctuated by bursts of innovation that decreased in amplitude over the past 66 million years. Social, precocial, aquatic, and herbivorous species evolve fastest, especially whales, elephants, sirenians, and extinct ungulates. Slow rates in rodents and bats indicate dissociation of taxonomic and morphological diversification. Frustratingly, highly similar ancestral shape estimates for placental mammal superorders suggest that their earliest representatives may continue to elude unequivocal identification.

Floral evolution and pollinator diversification in Hedychium: Revisiting Darwin's predictions using an integrative taxonomic approach

PREMISE

Hedychium J. Koenig (Zingiberaceae) is endemic to the Indo-Malayan Realm and is known for its colorful and fragrant flowers. Historically, two different pollination syndromes characterize Hedychium: diurnal or bird pollination, and nocturnal or moth pollination. In this study, we aim to understand the evolution of nocturnal and diurnal flowers, and to test its putative association with lineage diversification in Hedychium.

METHODS

A molecular tree of Hedychium was used as a scaffold upon which we estimated ancestral character states, phylogenetic signals, and correlations for certain categorical and continuous floral traits. Furthermore, we used phylomorphospace and trait-dependent diversification rate estimation analyses to understand phenotypic evolution and associated lineage diversification in Hedychium.

RESULTS

Although floral color and size lacked any association with specific pollinators, white or pale flowers were most common in the early branching clades when compared to bright-colored flowers, which were more widely represented in the most-derived clade IV. Five categorical and two continuous characters were identified to have informative evolutionary patterns, which also emphasized that ecology may have played a critical role in the diversification of Hedychium.

CONCLUSIONS

From our phylogenetic analyses and ecological observations, we conclude that specializations in pollinator interactions are rare in the hyperdiverse clade IV, thus challenging the role of both moth-specialization and bird-specialization as central factors in the diversification of Hedychium. However, our results also suggest that clade III (predominantly island clade) may show specializations, and future studies should investigate ecological and pollinator interactions, along with inclusion of new traits such as floral fragrance and anthesis time.

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.

Glaucoma through Animal’s Eyes: Insights from the Evolution of Intraocular Pressure in Mammals and Birds

Simple Summary

Understanding how a disease evolved across the animal kingdom could help us better understand the disease and might lead to novel methods for treatment. Here, we studied the evolution of glaucoma, an irreversible eye disease, in mammals and birds, by studying the evolution of intraocular pressure (IOP), a central driver of glaucoma, and searching for associations between life history traits and IOP. Our results revealed that IOP is a taxa-specific trait that is higher in some species than in others. Higher IOPs appear to have evolved multiple times in mammals and birds. Higher IOPs were found in mammals with higher body mass and in aquatic birds. We also found that higher IOPs evolved through stabilizing selection, with the optimum IOP in mammals and birds being 17.67 and 14.31 mmHg, respectively. This supports the hypothesis that higher IOPs may be an adaptive trait for certain animals. Focusing on species with higher IOPs but no evidence of glaucoma may help identify glaucoma-resistant adaptations, which could be developed into human therapies.

Abstract

Glaucoma, an eye disorder caused by elevated intraocular pressure (IOP), is the leading cause of irreversible blindness in humans. Understanding how IOP levels have evolved across animal species could shed light on the nature of human vulnerability to glaucoma. Here, we studied the evolution of IOP in mammals and birds and explored its life history correlates. We conducted a systematic review, to create a dataset of species-specific IOP levels and reconstructed the ancestral states of IOP using three models of evolution (Brownian, Early burst, and Ornstein–Uhlenbeck (OU)) to understand the evolution of glaucoma. Furthermore, we tested the association between life history traits (e.g., body mass, blood pressure, diet, longevity, and habitat) and IOP using phylogenetic generalized least squares (PGLS). IOP in mammals and birds evolved under the OU model, suggesting stabilizing selection toward an optimal value. Larger mammals had higher IOPs and aquatic birds had higher IOPs; no other measured life history traits, the type of tonometer used, or whether the animal was sedated when measuring IOP explained the significant variation in IOP in this dataset. Elevated IOP, which could result from physiological and anatomical processes, evolved multiple times in mammals and birds. However, we do not understand how species with high IOP avoid glaucoma. While we found very few associations between life history traits and IOP, we suggest that more detailed studies may help identify mechanisms by which IOP is decoupled from glaucoma. Importantly, species with higher IOPs (cetaceans, pinnipeds, and rhinoceros) could be good model systems for studying glaucoma-resistant adaptations.

Bridging macroecology and macroevolution in the radiation of sigmodontine rodents

Investigations of phenotypic disparity across geography often ignore macroevolutionary processes. As a corollary, the random null expectations to which disparity is compared and interpreted may be unrealistic. We tackle this issue by representing, in geographical space, distinct processes of phenotypic evolution underlying ecological disparity. Under divergent natural selection, assemblages in a given region should have empirical disparity higher than expected under an evolutionarily oriented null model, whereas the opposite may indicate constraints on phenotypic evolution. We gathered phylogenies, biogeographic distributions, and data on the skull morphology of sigmodontine rodents to discover which regions of the Neotropics were more influenced by divergent, neutral, or constrained phenotypic evolution. We found that regions with higher disparity than expected by the evolutionary-oriented null model, in terms of both size and shape, were concentrated in the Atlantic Forest, suggesting a larger role for divergent natural selection there. Phenotypic disparity in the rest of South America, mainly the Amazon basin, northeastern Brazil, and Southern Andes, was constrained-lower than predicted by the evolutionary model. We also demonstrated equivalence between the disparity produced by randomization-based null models and constrained-evolution null models. Therefore, including evolutionary simulations into the null modeling framework used in ecophylogenetics can strengthen inferences on the processes underlying phenotypic evolution.

Breaking constraints: The development and evolution of extreme fin morphology in the Bramidae

The developmental process establishes the foundation upon which natural selection may act. In that same sense, it is inundated with numerous constraints that work to limit the directions in which a phenotype may respond to selective pressures. Extreme phenotypes have been used in the past to identify tradeoffs and constraints and may aid in recognizing how alterations to the Baupläne can influence the trajectories of lineages. The Bramidae, a family of Scombriformes consisting of 20 extant species, are unique in that five species greatly deviate from the stout, ovaloid bodies that typify the bramids. The Ptericlinae, or fanfishes, are instead characterized by relatively elongated body plans and extreme modifications to their medial fins. Here, we explore the development of Bramidae morphologies and examine them through a phylogenetic lens to investigate the concepts of developmental and evolutionary constraints. Contrary to our predictions that the fanfishes had been constrained by inherited properties of an ancestral state, we find that the fanfishes exhibit both increased rates of trait evolution and differ substantially from the other bramids in their developmental trajectories. Conversely, the remaining bramid genera differ little, both among one another and in comparison, to the sister family Caristiidae. In all, our data suggest that the fanfishes have broken constraints, thereby allowing them to mitigate trade-offs on distinctive aspects of morphology.

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26-taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

Somewhere I belong: phylogeny and morphological evolution in a species-rich lineage of ectoparasitic flatworms infecting cichlid fishes

A substantial portion of biodiversity has evolved through adaptive radiation. However, the effects of explosive speciation on species interactions remain poorly understood. Metazoan parasites infecting radiating host lineages could improve our knowledge because of their intimate host relationships. Yet limited molecular, phenotypic and ecological data discourage multivariate analyses of evolutionary patterns and encourage the use of discrete characters. Here, we assemble new molecular, morphological and host range data widely inferred from a species-rich lineage of parasites (Cichlidogyrus, Platyhelminthes: Monogenea) infecting cichlid fishes to address data scarcity. We infer a multimarker (28S/18S rDNA, ITS1, COI mtDNA) phylogeny of 58 of 137 species and characterize major lineages through synapomorphies inferred from mapping morphological characters. We predict the phylogenetic position of species without DNA data through shared character states, a morphological phylogenetic analysis, and a classification analysis with support vector machines. Based on these predictions and a cluster analysis, we assess the systematic informativeness of continuous characters, search for continuous equivalents for discrete characters, and suggest new characters for morphological traits not analysed to date. We also model the attachment/reproductive organ and host range evolution using the data for 136 of 137 described species and multivariate phylogenetic comparative methods (PCMs). We show that discrete characters not only can mask phylogenetic signals, but also are key for characterizing species groups. Regarding the attachment organ morphology, a divergent evolutionary regime for at least one lineage was detected and a limited morphological variation indicates host and environmental parameters affecting its evolution. However, moderate success in predicting phylogenetic positions, and a low systematic informativeness and high multicollinearity of morphological characters call for a revaluation of characters included in species characterizations.

Rapid diversification of vascular architecture underlies the Carboniferous fern radiation

Vascular plants account for 93% of Earth's terrestrial flora. Xylem and phloem, vital for transporting water and nutrients through the plant, unite this diverse clade. Three-dimensional arrangements of these tissues (vascular architecture) are manifold across living and extinct species. However, the evolutionary processes underlying this variation remain elusive. Using ferns, a diverse clade with multiple radiations over their ca 400-million-year history, we synthesized data across 3339 species to explore the tempo and mode of vascular evolution and to contextualize dynamics of phenotypic innovation during major fern diversification events. Our results reveal three paradigm shifts in our understanding of fern vascular evolution. (i) The canonical theory on the stepwise and unidirectional evolution of vascular architecture does not capture the complexities of character evolution among ferns. Rather, a new model permitting additional transitions, rate heterogeneity and multiple reversions is more likely. (ii) Major shifts in vascular architecture correspond to developmental changes in body size, not regional water availability. (iii) The early Carboniferous radiation of crown-group ferns was characterized by an explosion of phenotypic innovation. By contrast, during the Cretaceous and Cenozoic rise of eupolypods, rates of vascular evolution were dramatically low and seemingly decoupled from lineage diversification.

Brawn before brains in placental mammals after the end-Cretaceous extinction

Mammals are the most encephalized vertebrates, with the largest brains relative to body size. Placental mammals have particularly enlarged brains, with expanded neocortices for sensory integration, the origins of which are unclear. We used computed tomography scans of newly discovered Paleocene fossils to show that contrary to the convention that mammal brains have steadily enlarged over time, early placentals initially decreased their relative brain sizes because body mass increased at a faster rate. Later in the Eocene, multiple crown lineages independently acquired highly encephalized brains through marked growth in sensory regions. We argue that the placental radiation initially emphasized increases in body size as extinction survivors filled vacant niches. Brains eventually became larger as ecosystems saturated and competition intensified.

Contrasting Patterns of Climatic Niche Divergence in Trebouxia-A Clade of Lichen-Forming Algae

Lichen associations are overwhelmingly supported by carbon produced by photosynthetic algal symbionts. These algae have diversified to occupy nearly all climates and continents; however, we have a limited understanding of how their climatic niches have evolved through time. Here we extend previous work and ask whether phylogenetic signal in, and the evolution of, climatic niche, varies across climatic variables, phylogenetic scales, and among algal lineages in Trebouxia—the most common genus of lichen-forming algae. Our analyses reveal heterogeneous levels of phylogenetic signal across variables, and that contrasting models of evolution underlie the evolution of climatic niche divergence. Together these analyses demonstrate the variable processes responsible for shaping climatic tolerance in Trebouxia, and provide a framework within which to better understand potential responses to climate change-associated perturbations. Such predictions reveal a disturbing trend in which the pace at which modern climate change is proceeding will vastly exceed the rate at which Trebouxia climatic niches have previously evolved.

General statistical model shows that macroevolutionary patterns and processes are consistent with Darwinian gradualism

Macroevolution posed difficulties for Darwin and later theorists because species’ phenotypes frequently change abruptly, or experience long periods of stasis, both counter to the theory of incremental change or gradualism. We introduce a statistical model that accommodates this uneven evolutionary landscape by estimating two kinds of historical change: directional changes that shift the mean phenotype along the branches of a phylogenetic tree, and evolvability changes that alter a clade’s ability to explore its trait-space. In mammals, we find that both processes make substantial independent contributions to explaining macroevolution, and are rarely linked. ‘Watershed’ moments of increased evolvability greatly outnumber reductions in evolutionary potentials, and large or abrupt phenotypic shifts are explicable statistically as biased random walks, allowing macroevolutionary theory to engage with the language and concepts of gradualist microevolution. Our findings recast macroevolutionary phenomena, illustrating the necessity of accounting for a variety of evolutionary processes simultaneously.

‘Macroevolution posed difficulties for Darwin and later theorists because species frequently change abruptly, or experience long periods of stasis, both counter to the theory of incremental change or gradualism. Here, the authors propose a macroevolutionary statistical model that accommodates this uneven evolutionary landscape, and shows how even abrupt macroevolutionary changes are compatible with gradualist microevolutionary processes.’

Competition and geography underlie speciation and morphological evolution in Indo-Australasian monitor lizards

How biotic and abiotic factors act together to shape biological diversity is a major question in evolutionary biology. The recent availability of large datasets and development of new methodological approaches provide new tools to evaluate the predicted effects of ecological interactions and geography on lineage diversification and phenotypic evolution. Here, we use a near complete phylogenomic-scale phylogeny and a comprehensive morphological dataset comprising more than a thousand specimens to assess the role of biotic and abiotic processes in the diversification of monitor lizards (Varanidae). This charismatic group of lizards shows striking variation in species richness among its clades and multiple instances of endemic radiation in Indo-Australasia (i.e., the Indo-Australian Archipelago and Australia), one of Earth's most biogeographically complex regions. We found heterogeneity in diversification dynamics across the family. Idiosyncratic biotic and geographic conditions appear to have driven diversification and morphological evolution in three endemic Indo-Australasian radiations. Furthermore, incumbency effects partially explain patterns in the biotic exchange between Australia and New Guinea. Our results offer insight into the dynamic history of Indo-Australasia, the evolutionary significance of competition, and the long-term consequences of incumbency effects.

Phylogenetic analysis of adaptation in comparative physiology and biomechanics: overview and a case study of thermal physiology in treefrogs

Comparative phylogenetic studies of adaptation are uncommon in biomechanics and physiology. Such studies require data collection from many species, a challenge when this is experimentally intensive. Moreover, researchers struggle to employ the most biologically appropriate phylogenetic tools for identifying adaptive evolution. Here, we detail an established but greatly underutilized phylogenetic comparative framework - the Ornstein-Uhlenbeck process - that explicitly models long-term adaptation. We discuss challenges in implementing and interpreting the model, and we outline potential solutions. We demonstrate use of the model through studying the evolution of thermal physiology in treefrogs. Frogs of the family Hylidae have twice colonized the temperate zone from the tropics, and such colonization likely involved a fundamental change in physiology due to colder and more seasonal temperatures. However, which traits changed to allow colonization is unclear. We measured cold tolerance and characterized thermal performance curves in jumping for 12 species of treefrogs distributed from the Neotropics to temperate North America. We then conducted phylogenetic comparative analyses to examine how tolerances and performance curves evolved and to test whether that evolution was adaptive. We found that tolerance to low temperatures increased with the transition to the temperate zone. In contrast, jumping well at colder temperatures was unrelated to biogeography and thus did not adapt during dispersal. Overall, our study shows how comparative phylogenetic methods can be leveraged in biomechanics and physiology to test the evolutionary drivers of variation among species.

The evolution of the traplining pollinator role in hummingbirds: specialization is not an evolutionary dead end

Trapliners are pollinators that visit widely dispersed flowers along circuitous foraging routes. The evolution of traplining in hummingbirds is thought to entail morphological specialization through the reciprocal coevolution of longer bills with the long-tubed flowers of widely dispersed plant species. Specialization, such as that exhibited by traplining hummingbirds, is often viewed as both irreversible and an evolutionary dead end. We tested these predictions in a macroevolutionary framework. Specifically, we assessed the relationship between beak morphology and foraging and tested whether transitions to traplining are irreversible and lead to lower rates of diversification as predicted by the hypothesis that specialization is an evolutionary dead end. We find that there have been multiple independent transitions to traplining across the hummingbird phylogeny, but reversals have been rare or incomplete at best. Multiple independent lineages of trapliners have become morphologically specialized, convergently evolving relatively large bills for their body size. Traplining is not an evolutionary dead end however, since trapliners continue to give rise to new traplining species at a rate comparable to non-trapliners.

High Blood Flow Into the Femur Indicates Elevated Aerobic Capacity in Synapsids Since the Synapsida-Sauropsida Split

Varanids are the only non-avian sauropsids that are known to approach the warm-blooded mammals in stamina. Furthermore, a much higher maximum metabolic rate (MMR) gives endotherms (including birds) higher stamina than crocodiles, turtles, and non-varanid lepidosaurs. This has led researchers to hypothesize that mammalian endothermy evolved as a second step after the acquisition of elevated MMR in non-mammalian therapsids from a plesiomorphic state of low metabolic rates. In recent amniotes, MMR correlates with the index of blood flow into the femur (Qi), which is calculated from femoral length and the cross-sectional area of the nutrient foramen. Thus, Qi may serve as an indicator of MMR range in extinct animals. Using the Qi proxy and phylogenetic eigenvector maps, here we show that elevated MMRs evolved near the base of Synapsida. Non-mammalian synapsids, including caseids, edaphosaurids, sphenacodontids, dicynodonts, gorgonopsids, and non-mammalian cynodonts, show Qi values in the range of recent endotherms and varanids, suggesting that raised MMRs either evolved in synapsids shortly after the Synapsida-Sauropsida split in the Mississippian or that the low MMR of lepidosaurs and turtles is apomorphic, as has been postulated for crocodiles.

Early giant reveals faster evolution of large body size in ichthyosaurs than in cetaceans

[Figure: see text].

Trophic niche shifts and phenotypic trait evolution are largely decoupled in Australasian parrots

BACKGROUND

Trophic shifts from one dietary niche to another have played major roles in reshaping the evolutionary trajectories of a wide range of vertebrate groups, yet their consequences for morphological disparity and species diversity differ among groups.

METHODS

Here, we use phylogenetic comparative methods to examine whether the evolution of nectarivory and other trophic shifts have driven predictable evolutionary pathways in Australasian psittaculid parrots in terms of ecological traits such as body size, beak shape, and dispersal capacity.

RESULTS

We found no evidence for an 'early-burst' scenario of lineage or morphological diversification. The best-fitting models indicate that trait evolution in this group is characterized by abrupt phenotypic shifts (evolutionary jumps), with no sign of multiple phenotypic optima correlating with different trophic strategies. Thus, our results point to the existence of weak directional selection and suggest that lineages may be evolving randomly or slowly toward adaptive peaks they have not yet reached.

CONCLUSIONS

This study adds to a growing body of evidence indicating that the relationship between avian morphology and feeding ecology may be more complex than usually assumed and highlights the importance of adding more flexible models to the macroevolutionary toolbox.

The evolution of body size in termites

Termites are social cockroaches. Because non-termite cockroaches are larger than basal termite lineages, which themselves include large termite species, it has been proposed that termites experienced a unidirectional body size reduction since they evolved eusociality. However, the validity of this hypothesis remains untested in a phylogenetic framework. Here, we reconstructed termite body size evolution using head width measurements of 1638 modern and fossil termite species. We found that the unidirectional body size reduction model was only supported by analyses excluding fossil species. Analyses including fossil species suggested that body size diversified along with speciation events and estimated that the size of the common ancestor of modern termites was comparable to that of modern species. Our analyses further revealed that body size variability among species, but not body size reduction, is associated with features attributed to advanced termite societies. Our results suggest that miniaturization took place at the origin of termites, while subsequent complexification of termite societies did not lead to further body size reduction.

The effects of introgression across thousands of quantitative traits revealed by gene expression in wild tomatoes

It is now understood that introgression can serve as powerful evolutionary force, providing genetic variation that can shape the course of trait evolution. Introgression also induces a shared evolutionary history that is not captured by the species phylogeny, potentially complicating evolutionary analyses that use a species tree. Such analyses are often carried out on gene expression data across species, where the measurement of thousands of trait values allows for powerful inferences while controlling for shared phylogeny. Here, we present a Brownian motion model for quantitative trait evolution under the multispecies network coalescent framework, demonstrating that introgression can generate apparently convergent patterns of evolution when averaged across thousands of quantitative traits. We test our theoretical predictions using whole-transcriptome expression data from ovules in the wild tomato genus Solanum. Examining two sub-clades that both have evidence for post-speciation introgression, but that differ substantially in its magnitude, we find patterns of evolution that are consistent with histories of introgression in both the sign and magnitude of ovule gene expression. Additionally, in the sub-clade with a higher rate of introgression, we observe a correlation between local gene tree topology and expression similarity, implicating a role for introgressed cis-regulatory variation in generating these broad-scale patterns. Our results reveal a general role for introgression in shaping patterns of variation across many thousands of quantitative traits, and provide a framework for testing for these effects using simple model-informed predictions.

It is now known from studying large genetic datasets that species often hybridize and cross with each other over many generations – a phenomenon known as introgression. Introgression introduces new genetic variation into a population, and this variation can cause traits to be shared among the introgressing species. When researchers study the evolution of trait variation among species, this source of trait sharing is rarely accounted for. Here, we present a statistical model of the effects of introgression on trait variation. This model predicts that, when averaged across many thousands of traits, introgressing species are consistently more similar than expected from standard approaches. Researchers studying gene expression often consider the expression of many thousands of genes, making this a case where the expected effects of introgression are likely to manifest. We tested our model prediction using ovule gene expression data from the wild tomato genus Solanum, in two groups of species with evidence of historical introgression. We found that patterns of expression similarity in both groups are consistent with their histories of introgression and the predictions from our model. Our results highlight the importance of accounting for introgression as a source of trait variation among species.

Phylogeny of the damselfishes (Pomacentridae) and patterns of asymmetrical diversification in body size and feeding ecology

The damselfishes (family Pomacentridae) inhabit near-shore communities in tropical and temperature oceans as one of the major lineages in coral reef fish assemblages. Our understanding of their evolutionary ecology, morphology and function has often been advanced by increasingly detailed and accurate molecular phylogenies. Here we present the next stage of multi-locus, molecular phylogenetics for the group based on analysis of 12 nuclear and mitochondrial gene sequences from 345 of the 422 damselfishes. The resulting well-resolved phylogeny helps to address several important questions about higher-level damselfish relationships, their evolutionary history and patterns of divergence. A time-calibrated phylogenetic tree yields a root age for the family of 55.5 mya, refines the age of origin for a number of diverse genera, and shows that ecological changes during the Eocene-Oligocene transition provided opportunities for damselfish diversification. We explored the idea that body size extremes have evolved repeatedly among the Pomacentridae, and demonstrate that large and small body sizes have evolved independently at least 40 times and with asymmetric rates of transition among size classes. We tested the hypothesis that transitions among dietary ecotypes (benthic herbivory, pelagic planktivory and intermediate omnivory) are asymmetric, with higher transition rates from intermediate omnivory to either planktivory or herbivory. Using multistate hidden-state speciation and extinction models, we found that both body size and dietary ecotype are significantly associated with patterns of diversification across the damselfishes, and that the highest rates of net diversification are associated with medium body size and pelagic planktivory. We also conclude that the pattern of evolutionary diversification in feeding ecology, with frequent and asymmetrical transitions between feeding ecotypes, is largely restricted to the subfamily Pomacentrinae in the Indo-West Pacific. Trait diversification patterns for damselfishes across a fully resolved phylogeny challenge many recent general conclusions about the evolution of reef fishes.

Evidence for selectively constrained 3D flower shape evolution in a Late Miocene clade of Malagasy Bulbophyllum orchids

Questions concerning the evolution of complex biological structures are central to the field of evolutionary biology. Yet, still little information is known about the modes and temporal dynamics of three-dimensional (3D) flower shape evolution across the history of clades. Here, we combined high-resolution X-ray computed tomography with 3D geometric morphometrics and phylogenetic comparative methods to test models of whole-flower shape evolution in the orchid family, using an early Late Miocene clade (c. 50 spp.) of Malagasy Bulbophyllum as model system. Based on landmark data of 38 species, our high-dimensional model fitting decisively rejects a purely neutral mode of evolution, suggesting instead that flower shapes evolved towards a primary adaptive optimum. Only a small number of recently evolved species/lineages attained alternative shape optima, resulting in an increased rate of phenotypic evolution. Our findings provide evidence of constrained 3D flower shape evolution in a small-sized clade of tropical orchids, resulting in low rates of phenotypic evolution and uncoupled trait-diversification rates. We hypothesise that this deep imprint of evolutionary constraint on highly complex floral structures might reflect long-term (directional and/or stabilizing) selection exerted by the group's main pollinators (flies).

An early burst in brachiopod evolution corresponding with significant climatic shifts during the Great Ordovician Biodiversification Event

We employ modified tip-dating methods to date divergence times within the Strophomenoidea, one of the most abundant and species-rich brachiopod clades to radiate during the Great Ordovician Biodiversification Event (GOBE), to determine if significant environmental changes at this time correlate with the diversification of the clade. Models using origination, extinction and sampling rates to estimate prior probabilities of divergence times strongly support both high rates of anatomical change per million years and rapid divergences shortly before the clade first appears in the fossil record. These divergence times indicate much higher rates of cladogenesis than are typical of brachiopods during this interval. The correspondence of high speciation rates and high anatomical disparity suggests punctuated (speciational) change drove the high frequencies of early anatomical change, which in turn suggests increased ecological opportunities rather than shifting developmental constraints account for high rates of anatomical change. The pulse of rapid evolution began coincident with cooling temperatures, the start of major oscillations in sea level and increased levels of atmospheric oxygen. Our results suggest that these factors permitted major geographical and ecological expansion of strophomenoids with intervals of geographical isolation, resulting in elevated speciation rates and corresponding elevated frequencies of punctuated change.

What does modularity mean?

Modularity is now generally recognized as a fundamental feature of organisms, one that may have profound consequences for evolution. Modularity has recently become a major focus of research in organismal biology across multiple disciplines including genetics, developmental biology, functional morphology, population and evolutionary biology. While the wealth of new data, and also new theory, has provided exciting and novel insights, the concept of modularity has become increasingly ambiguous. That ambiguity is underlain by diverse intuitions about what modularity means, and the ambiguity is not merely about the meaning of the word-the metrics of modularity are measuring different properties and the methods for delimiting modules delimit them by different, sometimes conflicting criteria. The many definitions, metrics and methods can lead to substantial confusion not just about what modularity means as a word but also about what it means for evolution. Here we review various concepts, using graphical depictions of modules. We then review some of the metrics and methods for analyzing modularity at different levels. To place these in theoretical context, we briefly review theories about the origins and evolutionary consequences of modularity. Finally, we show how mismatches between concepts, metrics and methods can produce theoretical confusion, and how potentially illogical interpretations can be made sensible by a better match between definitions, metrics, and methods.

3D skull modelling and description of a new baurusuchid (Crocodyliformes, Mesoeucrocodylia) from the Late Cretaceous (Bauru Basin) of Brazil

Baurusuchidae is one of the most diverse groups of South American notosuchians, unambiguously recorded in Late Cretaceous deposits of Brazil and Argentina. The group is characterized by a reduced tooth formula, a lateromedially compressed rostrum, and a verticalized quadrate, representing one of the top predators of their faunas. Historically, skull morphology is the most employed tool to investigate the relationships of baurusuchids, as most of the species have been primarily based on cranial remains. The present study describes a new baurusuchid species from the Bauru Basin of Brazil, based on the first tridimensional digital reconstruction of individualized skull bones for Notosuchia, and discusses its phylogenetic position within the group. The new species differs from all the other known baurusuchids by a depression on the posterior portion of the nasal bearing a crest, an infraorbital crest of the jugal that extends until the anterior margin of the lacrimal, the dorsal surface of the frontal lacking a longitudinal crest or depression, and the lateral convexity of the squamosal prongs participating in the occipital wall. The new taxon is consistently positioned as sister to the remaining baurusuchines, with Aplestosuchus sordidus and Stratiotosuchus maxhechti, as successive sister-taxa to a monophyletic Baurusuchus (Ba. albertoi, Ba. Salgadoensis, and Ba. pachecoi). Our updated phylogenetic analysis helps to differentiate the two major Baurusuchidae lineages, Baurusuchinae and Pissarrachampsinae. Yet, the new species shares morphological features with both groups, suggesting the occurrence of "Zones of Variability" in the radiation of Baurusuchidae.