Females drive postmating reproductive trait evolution across Drosophila species, but not via remating rate

While traits that contribute to premating sexual interactions are known to be wildly diverse, much less is known about the diversity of postmating (especially female) reproductive traits and the mechanisms shaping this diversity. To assess the rate, pattern, and potential drivers of postmating reproductive trait evolution, we analyzed male and female traits across up to 30 Drosophila species within a phylogenetic comparative framework. In addition to postmating reproductive morphology (e.g., sperm length, reproductive tract length and mass), we also quantified mating behaviors including female remating rate-a common proxy for the strength of postmating sexual selection. We found evidence for strong coevolution between male and female postmating traits (specifically sperm length and sperm storage organ size). However, remating rate was not associated with the rate of evolution or exaggeration of either male or female postmating reproductive morphology, once phylogenetic relatedness was accounted for. We infer that female-mediated and intersexual selection predominantly drive the evolution of our postmating morphological traits, including via divergent male and female interests in controlling paternity. In comparison, remating rate has a complex and likely secondary role in shaping this evolution, in part because this trait can be both a driver and a product of postmating selection.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

The Sequential Direct and Indirect Effects of Mountain Uplift, Climatic Niche, and Floral Trait Evolution on Diversification Dynamics in an Andean Plant Clade

Why and how organismal lineages radiate is commonly studied through either assessing abiotic factors (biogeography, geomorphological processes, and climate) or biotic factors (traits and interactions). Despite increasing awareness that both abiotic and biotic processes may have important joint effects on diversification dynamics, few attempts have been made to quantify the relative importance and timing of these factors, and their potentially interlinked direct and indirect effects, on lineage diversification. We here combine assessments of historical biogeography, geomorphology, climatic niche, vegetative, and floral trait evolution to test whether these factors jointly, or in isolation, explain diversification dynamics of a Neotropical plant clade (Merianieae, Melastomataceae). After estimating ancestral areas and the changes in niche and trait disparity over time, we employ Phylogenetic Path Analyses as a synthesis tool to test eleven hypotheses on the individual direct and indirect effects of these factors on diversification rates. We find strongest support for interlinked effects of colonization of the uplifting Andes during the mid-Miocene and rapid abiotic climatic niche evolution in explaining a burst in diversification rate in Merianieae. Within Andean habitats, later increases in floral disparity allowed for the exploitation of wider pollination niches (i.e., shifts from bee to vertebrate pollinators), but did not affect diversification rates. Our approach of including both vegetative and floral trait evolution, rare in assessments of plant diversification in general, highlights that the evolution of woody habit and larger flowers preceded the colonization of the Andes, but was likely critical in enabling the rapid radiation in montane environments. Overall, and in concert with the idea that ecological opportunity is a key element of evolutionary radiations, our results suggest that a combination of rapid niche evolution and trait shifts was critical for the exploitation of newly available niche space in the Andes in the mid-Miocene. Further, our results emphasize the importance of incorporating both abiotic and biotic factors into the same analytical framework if we aim to quantify the relative and interlinked effects of these processes on diversification.

Opportunity begets opportunity to drive macroevolutionary dynamics of a diverse lizard radiation

Evolution proceeds unevenly across the tree of life, with some lineages accumulating diversity more rapidly than others. Explaining this disparity is challenging as similar evolutionary triggers often do not result in analogous shifts across the tree, and similar shifts may reflect different evolutionary triggers. We used a combination of approaches to directly consider such context-dependency and untangle the complex network of processes that shape macroevolutionary dynamics, focusing on Pleurodonta, a diverse radiation of lizards. Our approach shows that some lineage-wide signatures are lost when conditioned on sublineages: while viviparity appears to accelerate diversification, its effect size is overestimated by its association with the Andean mountains. Conversely, some signals that erode at broader phylogenetic scales emerge at shallower ones. Mountains, in general, do not affect speciation rates; rather, the occurrence in the Andean mountains specifically promotes diversification. Likewise, the evolution of larger sizes catalyzes diversification rates, but only within certain ecological and geographical settings. We caution that conventional methods of fitting models to entire trees may mistakenly assign diversification heterogeneity to specific factors despite evidence against their plausibility. Our study takes a significant stride toward disentangling confounding factors and identifying plausible sources of ecological opportunities in the diversification of large evolutionary radiations.

Language follows a distinct mode of extra-genomic evolution

As one of the most specific, yet most diverse of human behaviors, language is shaped by both genomic and extra-genomic evolution. Sharing methods and models between these modes of evolution has significantly advanced our understanding of language and inspired generalized theories of its evolution. Progress is hampered, however, by the fact that the extra-genomic evolution of languages, i.e. linguistic evolution, maps only partially to other forms of evolution. Contrasting it with the biological evolution of eukaryotes and the cultural evolution of technology as the best understood models, we show that linguistic evolution is special by yielding a stationary dynamic rather than stable solutions, and that this dynamic allows the use of language change for social differentiation while maintaining its global adaptiveness. Linguistic evolution furthermore differs from technological evolution by requiring vertical transmission, allowing the reconstruction of phylogenies; and it differs from eukaryotic biological evolution by foregoing a genotype vs phenotype distinction, allowing deliberate and biased change. Recognising these differences will improve our empirical tools and open new avenues for analyzing how linguistic, cultural, and biological evolution interacted with each other when language emerged in the hominin lineage. Importantly, our framework will help to cope with unprecedented scientific and ethical challenges that presently arise from how rapid cultural evolution impacts language, most urgently from interventional clinical tools for language disorders, potential epigenetic effects of technology on language, artificial intelligence and linguistic communicators, and global losses of linguistic diversity and identity. Beyond language, the distinctions made here allow identifying variation in other forms of biological and cultural evolution, developing new perspectives for empirical research.

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.

Evolutionary paths toward multi-level convergence of lactic acid bacteria in fructose-rich environments

Convergence provides clues to unveil the non-random nature of evolution. Intermediate paths toward convergence inform us of the stochasticity and the constraint of evolutionary processes. Although previous studies have suggested that substantial constraints exist in microevolutionary paths, it remains unclear whether macroevolutionary convergence follows stochastic or constrained paths. Here, we performed comparative genomics for hundreds of lactic acid bacteria (LAB) species, including clades showing a convergent gene repertoire and sharing fructose-rich habitats. By adopting phylogenetic comparative methods we showed that the genomic convergence of distinct fructophilic LAB (FLAB) lineages was caused by parallel losses of more than a hundred orthologs and the gene losses followed significantly similar orders. Our results further suggested that the loss of adhE, a key gene for phenotypic convergence to FLAB, follows a specific evolutionary path of domain architecture decay and amino acid substitutions in multiple LAB lineages sharing fructose-rich habitats. These findings unveiled the constrained evolutionary paths toward the convergence of free-living bacterial clades at the genomic and molecular levels.

Urban birds' tolerance towards humans was largely unaffected by COVID-19 shutdown-induced variation in human presence

The coronavirus disease 2019 (COVID-19) pandemic and respective shutdowns dramatically altered human activities, potentially changing human pressures on urban-dwelling animals. Here, we use such COVID-19-induced variation in human presence to evaluate, across multiple temporal scales, how urban birds from five countries changed their tolerance towards humans, measured as escape distance. We collected 6369 escape responses for 147 species and found that human numbers in parks at a given hour, day, week or year (before and during shutdowns) had a little effect on birds' escape distances. All effects centered around zero, except for the actual human numbers during escape trial (hourly scale) that correlated negatively, albeit weakly, with escape distance. The results were similar across countries and most species. Our results highlight the resilience of birds to changes in human numbers on multiple temporal scales, the complexities of linking animal fear responses to human behavior, and the challenge of quantifying both simultaneously in situ.

Maternal investment evolves with larger body size and higher diversification rate in sharks and rays

Across vertebrates, live bearing evolved at least 150 times from ancestral egg laying into diverse forms and degrees of prepartum maternal investment.1,2 A key question is how reproductive diversity arose and whether reproductive diversification underlies species diversification.3,4,5,6,7,8,9,10,11 To test this, we evaluate the most basal jawed vertebrates: the sharks, rays, and chimaeras, which have one of the greatest ranges of reproductive and ecological diversity among vertebrates.2,12 We reconstruct the sequence of reproductive mode evolution across a phylogeny of 610 chondrichthyans.13 We reveal egg laying as ancestral, with live bearing evolving at least seven times. Matrotrophy evolved at least 15 times, with evidence of one reversal. In sharks, transitions to live bearing and matrotrophy are more prevalent in larger-bodied tropical species. Further, the evolution of live bearing is associated with a near doubling of the diversification rate, but there is only a small increase associated with the appearance of matrotrophy. Although pre-copulatory sexual selection is associated with increased rates of speciation in teleosts,3 sexual size dimorphism in chondrichthyans does not appear to be related to sexual selection,14,15 and instead we find increased rates of speciation associated with the colonization of novel habitats. This highlights a potential key difference between chondrichthyans and other fishes, specifically a slower rate of evolution of reproductive isolation following speciation, suggesting different rate-limiting mechanisms for diversification between these clades.16 The chondrichthyan diversification and radiation, particularly throughout shallow tropical shelf seas and oceanic pelagic habitats, appear to be associated with the evolution of live bearing and proliferation of a wide range of maternal investment in developing offspring.

The respiratory system influences flight mechanics in soaring birds

The subpectoral diverticulum (SPD) is an extension of the respiratory system in birds that is located between the primary muscles responsible for flapping the wing1,2. Here we survey the pulmonary apparatus in 68 avian species, and show that the SPD was present in virtually all of the soaring taxa investigated but absent in non-soarers. We find that this structure evolved independently with soaring flight at least seven times, which indicates that the diverticulum might have a functional and adaptive relationship with this flight style. Using the soaring hawks Buteo jamaicensis and Buteo swainsoni as models, we show that the SPD is not integral for ventilation, that an inflated SPD can increase the moment arm of cranial parts of the pectoralis, and that pectoralis muscle fascicles are significantly shorter in soaring hawks than in non-soaring birds. This coupling of an SPD-mediated increase in pectoralis leverage with force-specialized muscle architecture produces a pneumatic system that is adapted for the isometric contractile conditions expected in soaring flight. The discovery of a mechanical role for the respiratory system in avian locomotion underscores the functional complexity and heterogeneity of this organ system, and suggests that pulmonary diverticula are likely to have other undiscovered secondary functions. These data provide a mechanistic explanation for the repeated appearance of the SPD in soaring lineages and show that the respiratory system can be co-opted to provide biomechanical solutions to the challenges of flight and thereby influence the evolution of avian volancy.

Small- to medium-sized mammals show greater morphological disparity in cervical than lumbar vertebrae across different terrestrial modes of locomotion

During mammalian terrestrial locomotion, body flexibility facilitated by the vertebral column is expected to be correlated with observed modes of locomotion, known as gait (e.g., sprawl, trot, hop, bound, gallop). In small- to medium-sized mammals (average weight up to 5 kg), the relationship between locomotive mode and vertebral morphology is largely unexplored. Here we studied the vertebral column from 46 small- to medium-sized mammals. Nine vertebrae across cervical, thoracic, and lumbar regions were chosen to represent the whole vertebral column. Vertebra shape was analysed using three-dimensional geometric morphometrics with the phylogenetic comparative method. We also applied the multi-block method, which can consider all vertebrae as a single structure for analysis. We calculated morphological disparity, phylogenetic signal, and evaluated the effects of allometry and gait on vertebral shape. We also investigated the pattern of integration in the column. We found the cervical vertebrae show the highest degree of morphological disparity, and the first thoracic vertebra shows the highest phylogenetic signal. A significant effect of gait type on vertebrae shape was found, with the lumbar vertebrae having the strongest correlation; but this effect was not significant after taking phylogeny into account. On the other hand, allometry has a significant effect on all vertebrae regardless of the contribution from phylogeny. The regions showed differing degrees of integration, with cervical vertebrae most strongly correlated. With these results, we have revealed novel information that cannot be captured from study of a single vertebra alone: although the lumbar vertebrae are the most correlated with gait, the cervical vertebrae are more morphologically diverse and drive the diversity among species when considering whole column shape.

Connecting Species-Specific Extents of Genome Reduction in Mitochondria and Plastids

Mitochondria and plastids have both dramatically reduced their genomes since the endosymbiotic events that created them. The similarities and differences in the evolution of the two organelle genome types have been the target of discussion and investigation for decades. Ongoing work has suggested that similar mechanisms may modulate the reductive evolution of the two organelles in a given species, but quantitative data and statistical analyses exploring this picture remain limited outside of some specific cases like parasitism. Here, we use cross-eukaryote organelle genome data to explore evidence for coevolution of mitochondrial and plastid genome reduction. Controlling for differences between clades and pseudoreplication due to relatedness, we find that extents of mtDNA and ptDNA gene retention are related to each other across taxa, in a generally positive correlation that appears to differ quantitatively across eukaryotes, for example, between algal and nonalgal species. We find limited evidence for coevolution of specific mtDNA and ptDNA gene pairs, suggesting that the similarities between the two organelle types may be due mainly to independent responses to consistent evolutionary drivers.

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.

Facing the facts: adaptive trade-offs along body size ranges determine mammalian craniofacial scaling

The mammalian cranium (skull without lower jaw) is representative of mammalian diversity and is thus of particular interest to mammalian biologists across disciplines. One widely retrieved pattern accompanying mammalian cranial diversification is referred to as 'craniofacial evolutionary allometry' (CREA). This posits that adults of larger species, in a group of closely related mammals, tend to have relatively longer faces and smaller braincases. However, no process has been officially suggested to explain this pattern, there are many apparent exceptions, and its predictions potentially conflict with well-established biomechanical principles. Understanding the mechanisms behind CREA and causes for deviations from the pattern therefore has tremendous potential to explain allometry and diversification of the mammalian cranium. Here, we propose an amended framework to characterise the CREA pattern more clearly, in that 'longer faces' can arise through several kinds of evolutionary change, including elongation of the rostrum, retraction of the jaw muscles, or a more narrow or shallow skull, which all result in a generalised gracilisation of the facial skeleton with increased size. We define a standardised workflow to test for the presence of the pattern, using allometric shape predictions derived from geometric morphometrics analysis, and apply this to 22 mammalian families including marsupials, rabbits, rodents, bats, carnivores, antelopes, and whales. Our results show that increasing facial gracility with size is common, but not necessarily as ubiquitous as previously suggested. To address the mechanistic basis for this variation, we then review cranial adaptations for harder biting. These dictate that a more gracile cranium in larger species must represent a structural sacrifice in the ability to produce or withstand harder bites, relative to size. This leads us to propose that facial gracilisation in larger species is often a product of bite force allometry and phylogenetic niche conservatism, where more closely related species tend to exhibit more similar feeding ecology and biting behaviours and, therefore, absolute (size-independent) bite force requirements. Since larger species can produce the same absolute bite forces as smaller species with less effort, we propose that relaxed bite force demands can permit facial gracility in response to bone optimisation and alternative selection pressures. Thus, mammalian facial scaling represents an adaptive by-product of the shifting importance of selective pressures occurring with increased size. A reverse pattern of facial 'shortening' can accordingly also be found, and is retrieved in several cases here, where larger species incorporate novel feeding behaviours involving greater bite forces. We discuss multiple exceptions to a bite force-mediated influence on facial proportions across mammals which lead us to argue that ecomorphological specialisation of the cranium is likely to be the primary driver of facial scaling patterns, with some developmental constraints as possible secondary factors. A potential for larger species to have a wider range of cranial functions when less constrained by bite force demands might also explain why selection for larger sizes seems to be prevalent in some mammalian clades. The interplay between adaptation and constraint across size ranges thus presents an interesting consideration for a mechanistically grounded investigation of mammalian cranial allometry.

Foraging predicts the evolution of warning coloration and mimicry in snakes

Warning coloration and Batesian mimicry are classic examples of Darwinian evolution, but empirical evolutionary patterns are often paradoxical. We test whether foraging costs predict the evolution of striking coloration by integrating genetic and ecological data for aposematic and mimetic snakes (Elapidae and Dipsadidae). Our phylogenetic comparison on a total of 432 species demonstrated that dramatic changes in coloration were well predicted by foraging strategy. Multiple tests consistently indicated that warning coloration and conspicuous mimicry were more likely to evolve in species where foraging costs of conspicuous appearance were relaxed by poor vision of their prey, concealed habitat, or nocturnal activity. Reversion to crypsis was also well predicted by ecology for elapids but not for dipsadids. In contrast to a theoretical prediction and general trends, snakes' conspicuous coloration was correlated with secretive ecology, suggesting that a selection regime underlies evolutionary patterns. We also found evidence that mimicry of inconspicuous models (pitvipers) may have evolved in association with foraging demand for crypsis. These findings demonstrate that foraging is an important factor necessary to understand the evolution, persistence, and diversity of warning coloration and mimicry of snakes, highlighting the significance of additional selective factors in solving the warning coloration paradox.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Functionally mediated cranial allometry evidenced in a genus of rock-wallabies

In assessments of skeletal variation, allometry (disproportionate change of shape with size) is often corrected to examine size-independent variation for hypotheses relating to function. However, size-related trade-offs in functional demands may themselves be an underestimated driver of mammalian cranial diversity. Here, we use geometric morphometrics alongside dental measurements to assess craniodental allometry in the rock-wallaby genus Petrogale (all 17 species, 370 individuals). We identified functional aspects of evolutionary allometry that can be both extensions of, and correlated negatively with, static or ontogenetic allometric patterns. Regarding constraints, larger species tended to have relatively smaller braincases and more posterior orbits, the former of which might represent a constraint on jaw muscle anatomy. However, they also tended to have more anterior dentition and smaller posterior zygomatic arches, both of which support the hypothesis of relaxed bite force demands and accommodation of different selective pressures that favour facial elongation. By contrast, two dwarf species had stouter crania with divergent dental adaptations that together suggest increased relative bite force capacity. This likely allows them to feed on forage that is mechanically similar to that consumed by larger relatives. Our results highlight a need for nuanced considerations of allometric patterns in future research of mammalian cranial diversity.

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.

Transmission mode and dispersal traits correlate with host specificity in mammalian gut microbes

Different host species associate with distinct gut microbes in mammals, a pattern sometimes referred to as phylosymbiosis. However, the processes shaping this host specificity are not well understood. One model proposes that barriers to microbial transmission promote specificity by limiting microbial dispersal between hosts. This model predicts that specificity levels measured across microbes is correlated to transmission mode (vertical vs. horizontal) and individual dispersal traits. Here, we leverage two large publicly available gut microbiota data sets (1490 samples from 195 host species) to test this prediction. We found that host specificity varies widely across bacteria (i.e., there are generalist and specialist bacteria) and depends on transmission mode and dispersal ability. Horizontally-like transmitted bacteria equipped with traits that facilitate switches between host (e.g., tolerance to oxygen) were found to be less specific (more generalist) than microbes without those traits, for example, vertically-like inherited bacteria that are intolerant to oxygen. Altogether, our findings are compatible with a model in which limited microbial dispersal abilities foster host specificity.

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 models 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.

A Comparative Perspective on the Cerebello-Cerebral System and Its Link to Cognition

The longstanding idea that the cerebral cortex is the main neural correlate of human cognition can be elaborated by comparative analyses along the vertebrate phylogenetic tree that support the view that the cerebello-cerebral system is suited to support non-motor functions more generally. In humans, diverse accounts have illustrated cerebellar involvement in cognitive functions. Although the neocortex, and its transmodal association cortices such as the prefrontal cortex, have become disproportionately large over primate evolution specifically, human neocortical volume does not appear to be exceptional relative to the variability within primates. Rather, several lines of evidence indicate that the exceptional volumetric increase of the lateral cerebellum in conjunction with its connectivity with the cerebral cortical system may be linked to non-motor functions and mental operation in primates. This idea is supported by diverging cerebello-cerebral adaptations that potentially coevolve with cognitive abilities across other vertebrates such as dolphins, parrots, and elephants. Modular adaptations upon the vertebrate cerebello-cerebral system may thus help better understand the neuroevolutionary trajectory of the primate brain and its relation to cognition in humans. Lateral cerebellar lobules crura I-II and their reciprocal connections to the cerebral cortical association areas appear to have substantially expanded in great apes, and humans. This, along with the notable increase in the ventral portions of the dentate nucleus and a shift to increased relative prefrontal-cerebellar connectivity, suggests that modular cerebellar adaptations support cognitive functions in humans. In sum, we show how comparative neuroscience provides new avenues to broaden our understanding of cerebellar and cerebello-cerebral functions in the context of cognition.

The Major Features of Macroevolution

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

The evolution of same-sex sexual behaviour in mammals

Same-sex sexual behaviour has attracted the attention of many scientists working in disparate areas, from sociology and psychology to behavioural and evolutionary biology. Since it does not contribute directly to reproduction, same-sex sexual behaviour is considered an evolutionary conundrum. Here, using phylogenetic analyses, we explore the evolution of same-sex sexual behaviour in mammals. According to currently available data, this behaviour is not randomly distributed across mammal lineages, but tends to be particularly prevalent in some clades, especially primates. Ancestral reconstruction suggests that same-sex sexual behaviour may have evolved multiple times, with its appearance being a recent phenomenon in most mammalian lineages. Our phylogenetically informed analyses testing for associations between same-sex sexual behaviour and other species characteristics suggest that it may play an adaptive role in maintaining social relationships and mitigating conflict.

Detecting natural selection in trait-trait coevolution

No phenotypic trait evolves independently of all other traits, but the cause of trait-trait coevolution is poorly understood. While the coevolution could arise simply from pleiotropic mutations that simultaneously affect the traits concerned, it could also result from multivariate natural selection favoring certain trait relationships. To gain a general mechanistic understanding of trait-trait coevolution, we examine the evolution of 220 cell morphology traits across 16 natural strains of the yeast Saccharomyces cerevisiae and the evolution of 24 wing morphology traits across 110 fly species of the family Drosophilidae, along with the variations of these traits among gene deletion or mutation accumulation lines (a.k.a. mutants). For numerous trait pairs, the phenotypic correlation among evolutionary lineages differs significantly from that among mutants. Specifically, we find hundreds of cases where the evolutionary correlation between traits is strengthened or reversed relative to the mutational correlation, which, according to our population genetic simulation, is likely caused by multivariate selection. Furthermore, we detect selection for enhanced modularity of the yeast traits analyzed. Together, these results demonstrate that trait-trait coevolution is shaped by natural selection and suggest that the pleiotropic structure of mutation is not optimal. Because the morphological traits analyzed here are chosen largely because of their measurability and thereby are not expected to be biased with regard to natural selection, our conclusion is likely general.

Reed bamboos drive skull shape evolution in bush frogs of the Western Ghats, Peninsular India

Reed bamboo is a major ecological and economic resource for many animals, including humans. Nonetheless, the influence of this plant's evolutionary role on the morphology of animal species remains unexplored. Here, we investigated the significance of bamboo habitats as ecological opportunities in shaping the skull morphology of bush frogs (Raorchestes) from the Western Ghats, Peninsular India. We applied a three-dimensional (3D) geometric morphometric approach to capture the skull shape of 55 species of bush frogs. We visualized the skull shape variables in phylomorphospace with principal component analysis and performed phylogenetic generalized least-squares analysis to assess the impact of cranial size (evolutionary allometry) and habitat (bamboo or non-bamboo) on cranial shape. We quantified the morphological disparity between bamboo and non-bamboo bush frogs' skull shape, and employed RRphylo, a phylogenetic ridge regression method, to access the evolutionary rate and rate shifts of skull shape change. The phylomorphospace delineated bamboo and non-bamboo bush frogs. While cranial shape exhibited a significant but smaller association with size, its association with habitat type was non-significant. We detected, however, significant differences in skull shape between the two frog groups, with bamboo frogs showing higher morphological disparity and a remarkable shift in the evolutionary rate of skull shape diversification. These findings underscore the role of reed bamboo in the evolution of skull shape in the radiation of frogs, endemic to the Western Ghats. We demonstrate that the association between the members of two distinct endemic clades (bamboo reeds and bamboo frogs) is the outcome of a deep-time ecological opportunity that dates back to the Miocene.

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.

SOCIAL CONTEXT AND THE EVOLUTION OF DELAYED REPRODUCTION IN BIRDS

Classic life history theory makes generalized predictions about phenotypic correlations across large clades. Modern comparative tests of these correlations account for the underlying structure of phylogenetic trees. Yet neither life history theory nor phylogenetic comparative methods automatically specify how biological mechanisms generate correlations. This problem is evident in comparative analyses of birds. Birds show a correlation between body size and age at first reproduction, but do not actually grow larger if they delay reproduction. Instead, field studies raise the hypothesis that social contexts-especially cooperative breeding, coloniality, and lekking-generate unique demands for behavioral development, which in turn result in delayed reproduction. Here, we support that hypothesis with a comparative dataset spanning 961 species in 155 avian families. Continuous (Ornstein-Uhlenbeck), discrete (hidden state Markov), and phylogenetic regression models revealed delayed reproduction in colonial birds, a weaker signal in cooperative birds, and the consistent evolution of sexual bimaturism in polygynous, lekking birds. These results show an association between diverse social contexts, sex-specific developmental demands, and life history evolution in birds. Considering this diversity, we discuss how even statistically powerful phylogenetic correlations-whether focused on mass, lifespan, or broad social categories-can ultimately fail to model the history of life history evolution.

Conceptual and empirical bridges between micro- and macroevolution

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

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 (Kmult, 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 Kmult 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'.

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.

Opposing effects of plant traits on diversification

Species diversity can vary dramatically across lineages due to differences in speciation and extinction rates. Here, we explore the effects of several plant traits on diversification, finding that most traits have opposing effects on diversification. For example, outcrossing may increase the efficacy of selection and adaptation but also decrease mate availability, two processes with contrasting effects on lineage persistence. Such opposing trait effects can manifest as differences in diversification rates that depend on ecological context, spatiotemporal scale, and associations with other traits. The complexity of pathways linking traits to diversification suggests that the mechanistic underpinnings behind their correlations may be difficult to interpret with any certainty, and context dependence means that the effects of specific traits on diversification are likely to differ across multiple lineages and timescales. This calls for taxonomically and context-controlled approaches to studies that correlate traits and diversification.

Does haplodiploidy help drive the evolution of insect eusociality?

Understanding the evolution of eusociality in insects has been a long-standing and unsolved challenge in evolutionary biology. For decades, it has been suggested that haplodiploidy plays an important role in the origin of eusociality. However, some researchers have also suggested that eusociality is unrelated to haplodiploidy. Surprisingly, there have been no large-scale phylogenetic tests of this hypothesis (to our knowledge). Here, we test whether haplodiploidy might help explain the origins of eusociality across 874 hexapod families, using three different phylogenetic comparative methods. Two of the methods used support the idea that the evolution of eusociality is significantly associated with haplodiploidy, providing possibly the first phylogenetic support for this decades-old hypothesis across insects. However, some patterns were clearly discordant with this hypothesis, and one phylogenetic test was non-significant. Support for this hypothesis came largely from the repeated origins of eusociality within the haplodiploid hymenopterans (and within thrips). Experimental manipulations of the data show that the non-significant results are primarily explained by the origins of eusociality without haplodiploidy in some groups (i.e., aphids, termites). Overall, our results offer mixed phylogenetic support for the long-standing hypothesis that haplodiploidy helps drive the evolution of eusociality.

Key innovations and the diversification of Hymenoptera

The order Hymenoptera (wasps, ants, sawflies, and bees) represents one of the most diverse animal lineages, but whether specific key innovations have contributed to its diversification is still unknown. We assembled the largest time-calibrated phylogeny of Hymenoptera to date and investigated the origin and possible correlation of particular morphological and behavioral innovations with diversification in the order: the wasp waist of Apocrita; the stinger of Aculeata; parasitoidism, a specialized form of carnivory; and secondary phytophagy, a reversal to plant-feeding. Here, we show that parasitoidism has been the dominant strategy since the Late Triassic in Hymenoptera, but was not an immediate driver of diversification. Instead, transitions to secondary phytophagy (from parasitoidism) had a major influence on diversification rate in Hymenoptera. Support for the stinger and the wasp waist as key innovations remains equivocal, but these traits may have laid the anatomical and behavioral foundations for adaptations more directly associated with diversification.

Evolution of sexually selected traits across animals

Sexual selection is thought to be a major driver of phenotypic diversity and diversification in animals, but large-scale evolutionary patterns in sexually selected (SS) traits remain largely unknown. Here, we survey and analyze the evolution of these traits across animal phylogeny. We find that female mate choice appears to be the most widespread mechanism of sexual selection, but male-male competition appears to be almost as frequent in chordates and male mate choice is also common in arthropods. Among sensory types, tactile traits appear to be most widespread whereas auditory traits are relatively uncommon. Rather than being ubiquitous or randomly distributed across animals, most of these different types of SS traits are confined to clades in arthropods and chordates, which form “hotspots” for the evolution of these diverse trait types. Thus, different sensory types show accelerated rates of evolution in these clades. Moreover, different types of SS traits are strongly correlated with each other in their evolution across animals. Finally, despite the intensive interest in the role of sexual selection in speciation, we find only limited support for the idea that SS traits drive large-scale patterns of diversification and species richness across all animals.

Indigenous and introduced Collembola differ in desiccation resistance but not its plasticity in response to temperature

Biological invasions have significant ecological and economic impacts. Much attention is therefore focussed on predicting establishment and invasion success. Trait-based approaches are showing much promise, but are mostly restricted to investigations of plants. Although the application of these approaches to animals is growing rapidly, it is rare for arthropods and restricted mostly to investigations of thermal tolerance. Here we study the extent to which desiccation tolerance and its phenotypic plasticity differ between introduced (nine species) and indigenous (seven species) Collembola, specifically testing predictions of the 'ideal weed' and 'phenotypic plasticity' hypotheses of invasion biology. We do so on the F2 generation of adults in a full factorial design across two temperatures, to elicit desiccation responses, for the phenotypic plasticity trials. We also determine whether basal desiccation resistance responds to thermal laboratory natural selection. We first show experimentally that acclimation to different temperatures elicits changes to cuticular structure and function that are typically associated with water balance, justifying our experimental approach. Our main findings reveal that basal desiccation resistance differs, on average, between the indigenous and introduced species, but that this difference is weaker at higher temperatures, and is driven by particular taxa, as revealed by phylogenetic generalised least squares approaches. By contrast, the extent or form of phenotypic plasticity does not differ between the two groups, with a 'hotter is better' response being most common. Beneficial acclimation is characteristic of only a single species. Laboratory natural selection had little influence on desiccation resistance over 8-12 generations, suggesting that environmental filtering rather than adaptation to new environments may be an important factor influencing Collembola invasions.

Mating strategy predicts gene presence/absence patterns in a genus of simultaneously hermaphroditic flatworms

Gene repertoire turnover is a characteristic of genome evolution. However, we lack well-replicated analyses of presence/absence patterns associated with different selection contexts. Here, we study ∼100 transcriptome assemblies across Macrostomum, a genus of simultaneously hermaphroditic flatworms exhibiting multiple convergent shifts in mating strategy and associated reproductive morphologies. Many species mate reciprocally, with partners donating and receiving sperm at the same time. Other species convergently evolved to mate by hypodermic injection of sperm into the partner. We find that for orthologous transcripts annotated as expressed in the body region containing the testes, sequences from hypodermically inseminating species diverge more rapidly from the model species, Macrostomum lignano, and have a lower probability of being observed in other species. For other annotation categories, simpler models with a constant rate of similarity decay with increasing genetic distance from M. lignano match the observed patterns well. Thus, faster rates of sequence evolution for hypodermically inseminating species in testis-region genes result in higher rates of homology detection failure, yielding a signal of rapid evolution in sequence presence/absence patterns. Our results highlight the utility of considering appropriate null models for unobserved genes, as well as associating patterns of gene presence/absence with replicated evolutionary events in a phylogenetic context.

Among-species variation in hormone concentrations is associated with urban tolerance in birds

As cities expand across the globe, understanding factors that underlie variation in urban tolerance is vital for predicting changes in patterns of biodiversity. Endocrine traits, like circulating hormone concentrations and regulation of endocrine responses, might contribute to variation in species' ability to cope with urban challenges. For example, variation in glucocorticoid and androgen concentrations has been linked to life-history and behavioural traits that are associated with urban tolerance. However, we lack an understanding of the degree to which evolved differences in endocrine traits predict variation in urban tolerance across species. We analysed 1391 estimates of circulating baseline corticosterone, stress-induced corticosterone, and testosterone concentrations paired with citizen-science-derived urban occurrence scores in a broad comparative analysis of endocrine phenotypes across 71 bird species that differ in their occurrence in urban habitats. Our results reveal context-dependent links between baseline corticosterone and urban tolerance, as well as testosterone and urban tolerance. Stress-induced corticosterone was not related to urban tolerance. These findings suggest that some endocrine phenotypes contribute to a species' tolerance of urban habitats, but also indicate that other aspects of the endocrine phenotype, such as the ability to appropriately attenuate responses to urban challenges, might be important for success in cities.

Modeling the Evolution of Rates of Continuous Trait Evolution

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

Convergence without divergence in North American red-flowering Silene

Combinations of correlated floral traits have arisen repeatedly across angiosperms through convergent evolution in response to pollinator selection to optimize reproduction. While some plant groups exhibit very distinct combinations of traits adapted to specific pollinators (so-called pollination syndromes), others do not. Determining how floral traits diverge across clades and whether floral traits show predictable correlations in diverse groups of flowering plants is key to determining the extent to which pollinator-mediated selection drives diversification. The North American Silene section Physolychnis is an ideal group to investigate patterns of floral evolution because it is characterized by the evolution of novel red floral color, extensive floral morphological variation, polyploidy, and exposure to a novel group of pollinators (hummingbirds). We test for correlated patterns of trait evolution that would be consistent with convergent responses to selection in the key floral traits of color and morphology. We also consider both the role of phylogenic distance and geographic overlap in explaining patterns of floral trait variation. Inconsistent with phenotypically divergent pollination syndromes, we find very little clustering of North American Silene into distinct floral morphospace. We also find little evidence that phylogenetic history or geographic overlap explains patterns of floral diversity in this group. White- and pink-flowering species show extensive phenotypic diversity but are entirely overlapping in morphological variation. However, red-flowering species have much less phenotypic disparity and cluster tightly in floral morphospace. We find that red-flowering species have evolved floral traits that align with a traditional hummingbird syndrome, but that these trait values overlap with several white and pink species as well. Our findings support the hypothesis that convergent evolution does not always proceed through comparative phenotypic divergence, but possibly through sorting of standing ancestral variation.

A simple hierarchical model for heterogeneity in the evolutionary correlation on a phylogenetic tree

Numerous questions in phylogenetic comparative biology revolve around the correlated evolution of two or more phenotypic traits on a phylogeny. In many cases, it may be sufficient to assume a constant value for the evolutionary correlation between characters across all the clades and branches of the tree. Under other circumstances, however, it is desirable or necessary to account for the possibility that the evolutionary correlation differs through time or in different sections of the phylogeny. Here, we present a method designed to fit a hierarchical series of models for heterogeneity in the evolutionary rates and correlation of two quantitative traits on a phylogenetic tree. We apply the method to two datasets: one for different attributes of the buccal morphology in sunfishes (Centrarchidae); and a second for overall body length and relative body depth in rock- and non-rock-dwelling South American iguanian lizards. We also examine the performance of the method for parameter estimation and model selection using a small set of numerical simulations.

The Community Coevolution Model with Application to the Study of Evolutionary Relationships between Genes based on Phylogenetic Profiles

Organismal traits can evolve in a coordinated way, with correlated patterns of gains and losses reflecting important evolutionary associations. Discovering these associations can reveal important information about the functional and ecological linkages among traits. Phylogenetic profiles treat individual genes as traits distributed across sets of genomes and can provide a fine-grained view of the genetic underpinnings of evolutionary processes in a set of genomes. Phylogenetic profiling has been used to identify genes that are functionally linked, and to identify common patterns of lateral gene transfer in microorganisms. However, comparative analysis of phylogenetic profiles and other trait distributions should take into account the phylogenetic relationships among the organisms under consideration. Here we propose the Community Coevolution Model (CCM), a new coevolutionary model to analyze the evolutionary associations among traits, with a focus on phylogenetic profiles. In the CCM, traits are considered to evolve as a community with interactions, and the transition rate for each trait depends on the current states of other traits. Surpassing other comparative methods for pairwise trait analysis, CCM has the additional advantage of being able to examine multiple traits as a community to reveal more dependency relationships. We also develop a simulation procedure to generate phylogenetic profiles with correlated evolutionary patterns that can be used as benchmark data for evaluation purposes. A simulation study demonstrates that CCM is more accurate than other methods including the Jaccard Index and three tree-aware methods. The parameterization of CCM makes the interpretation of the relations between genes more direct, which leads to Darwin's scenario being identified easily based on the estimated parameters. We show that CCM is more efficient and fits real data better than other methods resulting in higher likelihood scores with fewer parameters. An examination of 3786 phylogenetic profiles across a set of 659 bacterial genomes highlights linkages between genes with common functions, including many patterns that would not have been identified under a non-phylogenetic model of common distribution. We also applied the CCM to 44 proteins in the well-studied Mitochondrial Respiratory Complex I and recovered associations that mapped well onto the structural associations that exist in the complex.

A flexible method for estimating tip diversification rates across a range of speciation and extinction scenarios

Estimates of diversification rates at the tips of a phylogeny provide a flexible approach for correlation analyses with multiple traits and to map diversification rates in space while also avoiding the uncertainty of deep time rate reconstructions. Available methods for tip rate estimation make different assumptions, and thus their accuracy usually depends on the characteristics of the underlying model generating the tree. Here, we introduce MiSSE, a trait-free, state-dependent speciation and extinction approach that can be used to estimate varying speciation, extinction, net diversification, turnover, and extinction fractions at the tips of the tree. We compare the accuracy of tip rates inferred by MiSSE against similar methods and demonstrate that, due to certain characteristics of the model, the error is generally low across a broad range of speciation and extinction scenarios. MiSSE can be used alongside regular phylogenetic comparative methods in trait-related diversification hypotheses, and we also describe a simple correction to avoid pseudoreplication from sister tips in analyses of independent contrasts. Finally, we demonstrate the capabilities of MiSSE, with a renewed focus on classic comparative methods, to examine the correlation between plant height and turnover rates in eucalypts, a species-rich lineage of flowering plants.

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.

Fast Likelihood Calculations for Automatic Identification of Macroevolutionary Rate Heterogeneity in Continuous and Discrete Traits

Understanding phenotypic disparity across the tree of life requires identifying where and when evolutionary rates change on phylogeny. A primary methodological challenge in macroevolution is therefore to develop methods for accurate inference of among-lineage variation in rates of phenotypic evolution. Here, we describe a method for inferring among-lineage evolutionary rate heterogeneity in both continuous and discrete traits. The method assumes that the present-day distribution of a trait is shaped by a variable-rate process arising from a mixture of constant-rate processes and uses a single-pass tree traversal algorithm to estimate branch-specific evolutionary rates. By employing dynamic programming optimization techniques and approximate maximum likelihood estimators where appropriate, our method permits rapid exploration of the tempo and mode of phenotypic evolution. Simulations indicate that the method reconstructs rates of trait evolution with high accuracy. Application of the method to datasets on squamate reptile reproduction and turtle body size recovers patterns of rate heterogeneity identified by previous studies but with computational costs reduced by many orders of magnitude. Our results expand the set of tools available for detecting macroevolutionary rate heterogeneity and point to the utility of fast, approximate methods for studying large scale biodiversity dynamics.

How to approach the study of syndromes in macroevolution and ecology

Syndromes, wherein multiple traits evolve convergently in response to a shared selective driver, form a central concept in ecology and evolution. Recent work has questioned the existence of some classic syndromes, such as pollination and seed dispersal syndromes. Here, we discuss some of the major issues that have afflicted research into syndromes in macroevolution and ecology. First, correlated evolution of traits and hypothesized selective drivers is often relied on as the only evidence for adaptation of those traits to those hypothesized drivers, without supporting evidence. Second, the selective driver is often inferred from a combination of traits without explicit testing. Third, researchers often measure traits that are easy for humans to observe rather than measuring traits that are suited to testing the hypothesis of adaptation. Finally, species are often chosen for study because of their striking phenotypes, which leads to the illusion of syndromes and divergence. We argue that these issues can be avoided by combining studies of trait variation across entire clades or communities with explicit tests of adaptive hypotheses and that taking this approach will lead to a better understanding of syndrome-like evolution and its drivers.

Flight hampers the evolution of weapons in birds

Birds are a remarkable example of how sexual selection can produce diverse ornaments and behaviours. Specialised fighting structures like deer's antlers, in contrast, are mostly absent among birds. Here, we investigated if the birds' costly mode of locomotion-powered flight-helps explain the scarcity of weapons among members of this clade. Our simulations of flight energetics predicted that the cost of bony spurs-a specialised avian weapon-should increase with time spent flying. Bayesian phylogenetic comparative analyses using a global spur dataset corroborated this prediction. First, extant species with flight-efficient wings (which presumably fly more frequently) tend to have fewer or no bony spurs. Second, this association likely arose because flying more leads to more frequent evolutionary loss of spurs. Together, these findings suggest that, much like pneumatic bones, absence of weaponry may be another feature of the avian body plan that allows birds to efficiently explore the aerial habitat.

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.