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

Evolutionary bursts drive morphological novelty in the world's largest skinks

Animal phenotypes evolve and diverge as a result of differing selective pressures and drift. These processes leave unique signatures in patterns of trait evolution, impacting the tempo and mode of morphological macroevolution. While there is a broad understanding of the history of some organismal traits (e.g., body size), there is little consensus about the evolutionary mode of most others. This includes the relative contribution of prolonged (Darwinian gradualist) and episodic (Simpsonian jump) changes toward the evolution of novel morphologies. Here, we use new exon-capture and linear morphological datasets to investigate the tempo and mode of morphological evolution in Australo-Melanesian Tiliquini skinks. We generate a well-supported time-calibrated phylogenomic tree from ∼400 nuclear markers for more than 100 specimens, including undescribed diversity, and provide unprecedented resolution of the rapid Miocene diversification of these lizards. By collecting a morphological dataset that encompasses the lizard body plan (19 traits across the head, body, limb, and tail), we are able to identify that most traits evolve conservatively, but infrequent evolutionary bursts result in morphological novelty. These phenotypic discontinuities occur via rapid rate increases along individual branches, inconsistent with both gradualistic and punctuated equilibrial evolutionary modes. Instead, this "punctuated gradualism" has resulted in the rapid evolution of blue-tongued giants and armored dwarves in the ∼20 million years since colonizing Australia. These results outline the evolutionary pathway toward new morphologies and highlight the heterogeneity of evolutionary tempo and mode, even within individual traits.

Co-evolutionary dynamics of mammalian brain and body size

Despite decades of comparative studies, puzzling aspects of the relationship between mammalian brain and body mass continue to defy satisfactory explanation. Here we show that several such aspects arise from routinely fitting log-linear models to the data: the correlated evolution of brain and body mass is in fact log-curvilinear. This simultaneously accounts for several phenomena for which diverse biological explanations have been proposed, notably variability in scaling coefficients across clades, low encephalization in larger species and the so-called taxon-level problem. Our model implies a need to revisit previous findings about relative brain mass. Accounting for the true scaling relationship, we document dramatically varying rates of relative brain mass evolution across the mammalian phylogeny, and we resolve the question of whether there is an overall trend for brain mass to increase through time. We find a trend in only three mammalian orders, which is by far the strongest in primates, setting the stage for the uniquely rapid directional increase ultimately producing the computational powers of the human brain.

Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

Primate social organization evolved from a flexible pair-living ancestor

Explaining the evolution of primate social organization has been fundamental to understand human sociality and social evolution more broadly. It has often been suggested that the ancestor of all primates was solitary and that other forms of social organization evolved later, with transitions being driven by various life history traits and ecological factors. However, recent research showed that many understudied primate species previously assumed to be solitary actually live in pairs, and intraspecific variation in social organization is common. We built a detailed database from primary field studies quantifying the number of social units expressing different social organizations in each population. We used Bayesian phylogenetic models to infer the probability of each social organization, conditional on several socioecological and life history predictors. Here, we show that when intraspecific variation is accounted for, the ancestral social organization of primates was inferred to be variable, with the most common social organization being pair-living but with approximately 10 to 20% of social units of the ancestral population deviating from this pattern by being solitary living. Body size and activity patterns had large effects on transitions between types of social organizations. As in other mammalian clades, pair-living is closely linked to small body size and likely more common in ancestral species. Our results challenge the assumption that ancestral primates were solitary and that pair-living evolved afterward emphasizing the importance of focusing on field data and accounting for intraspecific variation, providing a flexible statistical framework for doing so.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

How Important Is Budding Speciation for Comparative Studies?

The acknowledgment of evolutionary dependence among species has fundamentally changed how we ask biological questions. Phylogenetic models became the standard approach for studies with 3 or more lineages, in particular those using extant species. Most phylogenetic comparative methods (PCMs) translate relatedness into covariance, meaning that evolutionary changes before lineages split should be interpreted together whereas after the split lineages are expected to change independently. This clever realization has shaped decades of research. Here, we discuss one element of the comparative method often ignored or assumed as unimportant: if nodes of a phylogeny represent the dissolution of the ancestral lineage into two new ones or if the ancestral lineage can survive speciation events (i.e., budding). Budding speciation is often reported in paleontological studies, due to the nature of the evidence for budding in the fossil record, but it is surprisingly absent in comparative methods. Here, we show that many PCMs assume that divergence happens as a symmetric split, even if these methods do not explicitly mention this assumption. We discuss the properties of trait evolution models for continuous and discrete traits and their adequacy under a scenario of budding speciation. We discuss the effects of budding speciation under a series of plausible evolutionary scenarios and show when and how these can influence our estimates. We also propose that long-lived lineages that have survived through a series of budding speciation events and given birth to multiple new lineages can produce evolutionary patterns that challenge our intuition about the most parsimonious history of trait changes in a clade. We hope our discussion can help bridge comparative approaches in paleontology and neontology as well as foster awareness about the assumptions we make when we use phylogenetic trees.

The macroevolutionary dynamics of mammalian sexual size dimorphism

Sexual size dimorphism (SSD) is a common phenomenon across the animal kingdom. Mammals are unusual in primarily displaying male-biased SSD, where males of a species are typically larger than females. The driving factors behind the evolution of this SSD have been much debated, with popular hypotheses invoking the influence of mating system and social organization via sexual selection, dietary niche divergence and broad-scale correlations with body size (Rensch's rule). Here, we investigate the macroevolutionary origins and maintenance of SSD among mammals, using phylogenetic general mixed linear models and a comprehensive global dataset to evaluate correlations of diet, body mass, seasonality, social organization and mating system with SSD type. We find that SSD as a whole is lost at a greater rate than it is gained, with female-biased SSD being particularly unstable. Non-monogamous mating systems, vertebrate prey consumption and temperature seasonality correlate with male-biased SSD, while polyandry correlates with female-biased SSD, and both types of SSD are positively correlated with body mass. This is in partial contrast to the predictions of Rensch's rule, which predicts that female-biased SSD would correlate negatively with body size. Taken together, our results highlight the importance of considering multiple ecological and social drivers when evaluating the macroevolutionary trajectory of sex differences in body size.

A novel method for identifying key genes in macroevolution based on deep learning with attention mechanism

Macroevolution can be regarded as the result of evolutionary changes of synergistically acting genes. Unfortunately, the importance of these genes in macroevolution is difficult to assess and hence the identification of macroevolutionary key genes is a major challenge in evolutionary biology. In this study, we designed various word embedding libraries of natural language processing (NLP) considering the multiple mechanisms of evolutionary genomics. A novel method (IKGM) based on three types of attention mechanisms (domain attention, kmer attention and fused attention) were proposed to calculate the weights of different genes in macroevolution. Taking 34 species of diurnal butterflies and nocturnal moths in Lepidoptera as an example, we identified a few of key genes with high weights, which annotated to the functions of circadian rhythms, sensory organs, as well as behavioral habits etc. This study not only provides a novel method to identify the key genes of macroevolution at the genomic level, but also helps us to understand the microevolution mechanisms of diurnal butterflies and nocturnal moths in Lepidoptera.

Innovation and elaboration on the avian tree of life

Widely documented, megaevolutionary jumps in phenotypic diversity continue to perplex researchers because it remains unclear whether these marked changes can emerge from microevolutionary processes. Here, we tackle this question using new approaches for modeling multivariate traits to evaluate the magnitude and distribution of elaboration and innovation in the evolution of bird beaks. We find that elaboration, evolution along the major axis of phenotypic change, is common at both macro- and megaevolutionary scales, whereas innovation, evolution away from the major axis of phenotypic change, is more prominent at megaevolutionary scales. The major axis of phenotypic change among species beak shapes at megaevolutionary scales is an emergent property of innovation across clades. Our analyses suggest that the reorientation of phenotypes via innovation is a ubiquitous route for divergence that can arise through gradual change alone, opening up further avenues for evolution to explore.

The evolution of gestation length in eutherian mammals

Gestation length, or the duration of pregnancy, is a critical component of mammalian reproductive biology1. Eutherian mammals exhibit striking variation in their gestation lengths2-5, which has traditionally been linked to and allometrically scales with variation in other life history traits, including body mass and lifespan5-8. How the phenotypic landscape of gestation length variation, including its associations with body mass and lifespan variation, changed over mammalian evolution remains unknown. Phylogeny-informed analyses of 845 representative extant eutherian mammals showed that gestation length variation substantially differed in both whether and how strongly it was associated with body mass and lifespan across mammalian clades. For example, gestation length variation in Chiroptera and Cetacea was not associated with lifespan or body mass but was strongly associated only with body mass in Carnivora. We also identified 52 adaptive shifts in gestation length variation across the mammal phylogeny and 14 adaptive shifts when considering all three life history traits; the placements of six adaptive shifts are common in the two analyses. Notably, two of these shifts occurred at the roots of Cetacea and Pinnipedia, respectively, coinciding with the transition of these clades to the marine environment. The varying dynamics of the phenotypic landscape of gestation length, coupled with the varying patterns of associations between gestation length and two other major life history traits, raise the hypothesis that evolutionary constraints on gestation length have varied substantially across mammalian phylogeny. This variation in constraints implies that the genetic architecture of gestation length differs between mammal clades.

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.

Extended embryo retention and viviparity in the first amniotes

The amniotic egg with its complex fetal membranes was a key innovation in vertebrate evolution that enabled the great diversification of reptiles, birds and mammals. It is debated whether these fetal membranes evolved in eggs on land as an adaptation to the terrestrial environment or to control antagonistic fetal-maternal interaction in association with extended embryo retention (EER). Here we report an oviparous choristodere from the Lower Cretaceous period of northeast China. The ossification sequence of the embryo confirms that choristoderes are basal archosauromorphs. The discovery of oviparity in this assumed viviparous extinct clade, together with existing evidence, suggests that EER was the primitive reproductive mode in basal archosauromorphs. Phylogenetic comparative analyses on extant and extinct amniotes suggest that the first amniote displayed EER (including viviparity).

The dynamic adaptive landscape of cetacean body size

Adaptive landscapes are central to evolutionary theory, forming a conceptual bridge between micro- and macroevolution.^1,2,3,4 Evolution by natural selection across an adaptive landscape should drive lineages toward fitness peaks, shaping the distribution of phenotypic variation within and among clades over evolutionary timescales.⁵ The location and breadth of these peaks in phenotypic space can also evolve,⁴ but whether phylogenetic comparative methods can detect such patterns has largely remained unexplored.⁶ Here, we characterize the global and local adaptive landscape for total body length in cetaceans (whales, dolphins, and relatives), a trait that spans 5 orders of magnitude, across their ∼53 million year evolutionary history. Using phylogenetic comparative methods, we analyze shifts in long-term mean body length⁷ and directional changes in average trait values⁸ for 345 living and fossil cetacean taxa. Remarkably, we find that the global macroevolutionary adaptive landscape of cetacean body length is relatively flat, with very few peak shifts occurring after cetaceans entered the oceans. Local peaks are more numerous and manifest as trends along branches linked to specific adaptations. These results contrast with previous studies using only extant taxa,⁹ highlighting the vital role of fossil data for understanding macroevolution.^10,11,12 Our results indicate that adaptive peaks are dynamic and are associated with subzones of local adaptations, creating moving targets for species adaptation. In addition, we identify limits in our ability to detect some evolutionary patterns and processes and suggest that multiple approaches are required to characterize complex hierarchical patterns of adaptation in deep time.

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.

Inherency and agency in the origin and evolution of biological functions

Although discussed by 20th century philosophers in terms drawn from the sciences of non-living systems, in recent decades biological function has been considered in relationship to organismal capability and purpose. Bringing two phenomena generally neglected in evolutionary theory (i.e. inherency and agency) to bear on questions of function leads to a rejection of the adaptationist ‘selected effects’ notion of biological function. I review work showing that organisms such as the placozoans can thrive with almost no functional embellishments beyond those of their constituent cells and physical properties of their simple tissues. I also discuss work showing that individual tissue cells and their artificial aggregates exhibit agential behaviours that are unprecedented in the histories of their respective lineages. I review findings on the unique metazoan mechanism of developmental gene expression that has recruited, during evolution, inherent ancestral cellular functionalities into specialized cell types and organs of the different animal groups. I conclude that most essential functions in animal species are inherent to the cells from which they evolved, not selected effects, and that many of the others are optional ‘add-ons’, a status inimical to fitness-based models of evolution positing that traits emerge from stringent cycles of selection to meet external challenges.

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.