The signature of competition in ecomorphological traits across the avian radiation

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

Ecomorphological variation in male and female wall lizards and the macroevolution of sexual dimorphism in relation to habitat use

Understanding how phenotypic diversity evolves is a major interest of evolutionary biology. Habitat use is an important factor in the evolution of phenotypic diversity of many animal species. Interestingly, male and female phenotypes have been frequently shown to respond differently to environmental variation. At the macroevolutionary level, this difference between the sexes is frequently analysed using phylogenetic comparative tools to assess variation in sexual dimorphism (SD) across taxa in relation to habitat. A shortcoming of such analyses is that they evaluate the degree of dimorphism itself and therefore they do not provide access to the evolutionary trajectories of each sex. As such, the relative contribution of male and female phenotypes on macroevolutionary patterns of sexual dimorphism cannot be directly assessed. Here, we investigate how habitat use shapes phenotypic diversity in wall lizards using phylogenetic comparative tools to simultaneously assess the tempo and mode of evolution in males, females and the degree of sexual dimorphism. We find that both sexes have globally diversified under similar, but not identical, processes, where habitat use seems to drive macroevolutionary variation in head shape, but not in body size or relative limb length. However, we also observe small differences in the evolutionary dynamics of male and female phenotypes that have a marked impact on macroevolutionary patterns of SD, with important implications for our interpretation of what drives phenotypic diversification within and between the sexes.     

How traits shape trees: new approaches for detecting character state‐dependent lineage diversification

Biologists have long sought to understand the processes underlying disparities in clade size across the tree of life and the extent to which such clade size differences can be attributed to the evolution of particular traits. The association of certain character states with species‐rich clades suggests that trait evolution can lead to increased diversification, but such a pattern could also arise due other processes, such as directional trait evolution. Recent advances in phylogenetic comparative methods have provided new statistical approaches for distinguishing between these intertwined and potentially confounded macroevolutionary processes. Here, we review the historical development of methods for detecting state‐dependent diversification and explore what new methods have revealed about classic examples of traits that affect diversification, including evolutionary dead ends, key innovations and geographic traits. Applications of these methods thus far collectively suggest that trait diversity commonly arises through the complex interplay between transition, speciation and extinction rates and that long hypothesized evolutionary dead ends and key innovations are instead often cases of directional trends in trait evolution.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

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

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

Macroevolution of arboreality in salamanders

Evolutionary theory predicts that selection in distinct microhabitats generates correlations between morphological and ecological traits, and may increase both phenotypic and taxonomic diversity. However, some microhabitats exert unique selective pressures that act as a restraining force on macroevolutionary patterns of diversification. In this study, we use phylogenetic comparative methods to investigate the evolutionary outcomes of inhabiting the arboreal microhabitat in salamanders. We find that arboreality has independently evolved at least five times in Caudata and has arisen primarily from terrestrial ancestors. However, the rate of transition from arboreality back to terrestriality is 24 times higher than the converse. This suggests that macroevolutionary trends in microhabitat use tend toward terrestriality over arboreality, which influences the extent to which use of the arboreal microhabitat proliferates. Morphologically, we find no evidence for an arboreal phenotype in overall body proportions or in foot shape, as variation in both traits overlaps broadly with species that utilize different microhabitats. However, both body shape and foot shape display reduced rates of phenotypic evolution in arboreal taxa, and evidence of morphological convergence among arboreal lineages is observed. Taken together, these patterns suggest that arboreality has played a unique role in the evolution of this family, providing neither an evolutionary opportunity, nor an evolutionary dead end.     

The evolution of floral signals in relation to range overlap in a clade of California Jewelflowers (Streptanthus s.l.)

Because of their function as reproductive signals in plants, floral traits experience distinct selective pressures related to their role in speciation, reinforcement, and prolonged coexistence with close relatives. However, few studies have investigated whether population‐level processes translate into detectable signatures at the macroevolutionary scale. Here, we ask whether patterns of floral trait evolution and range overlap across a clade of California Jewelflowers reflect processes hypothesized to shape floral signal differentiation at the population level. We found a pattern of divergence in floral scent composition across the clade such that close relatives had highly disparate floral scents given their age. Accounting for range overlap with close relatives explained additional variation in floral scent over time, with sympatric species pairs having diverged more than allopatric species pairs given their age. However, three other floral traits (flower size, scent complexity and flower color) did not fit these patterns, failing to deviate from a null Brownian motion model of evolution. Together, our results suggest that selection for divergence among close relatives in the composition of floral scents may play a key, sustained role in mediating speciation and coexistence dynamics across this group, and that signatures of these dynamics may persist at the macroevolutionary scale.     

Correlates of rate heterogeneity in avian ecomorphological traits

Heterogeneity in rates of trait evolution is widespread, but it remains unclear which processes drive fast and slow character divergence across global radiations. Here, we test multiple hypotheses for explaining rate variation in an ecomorphological trait (beak shape) across a globally distributed group (birds). We find low support that variation in evolutionary rates of species is correlated with life history, environmental mutagenic factors, range size, number of competitors, or living on islands. Indeed, after controlling for the negative effect of species' age, 80% of variation in species‐specific evolutionary rates remains unexplained. At the clade level, high evolutionary rates are associated with unusual phenotypes or high species richness. Taken together, these results imply that macroevolutionary rates of ecomorphological traits are governed by both ecological opportunity in distinct adaptive zones and niche differentiation among closely related species.     

Trait independence primes convergent trait loss

The repeated, independent evolution of traits (convergent evolution) is often attributed to shared environmental selection pressures. However, developmental dependencies among traits can limit the phenotypic variation available to selection and bias evolutionary outcomes. Here, we determine how changes in developmentally correlated traits may impact convergent loss of the tympanic middle ear, a highly labile trait within toads that currently lack adaptive explanation. The middle ear's lability could reflect evolutionary trade‐offs with other skull features under selection, or the middle ear may evolve independently of the rest of the skull, allowing it to be modified by active or passive processes without pleiotropic trade‐offs with other skull features. We compare the skulls of 55 species (39 eared, 16 earless) within the family Bufonidae, spanning six hypothesized independent middle ear transitions. We test whether shared or lineage‐specific changes in skull shape distinguish earless species from eared species and whether earless skulls lack other late‐forming skull bones. We find no evidence for pleiotropic trade‐offs between the middle ear and other skull structures. Instead, middle ear loss in anurans may provide a rare example of developmental independence contributing to evolutionary lability of a sensory system.     

Approaches to Macroevolution: 1. General Concepts and Origin of Variation

Approaches to macroevolution require integration of its two fundamental components, i.e. the origin and the sorting of variation, in a hierarchical framework. Macroevolution occurs in multiple currencies that are only loosely correlated, notably taxonomic diversity, morphological disparity, and functional variety. The origin of variation within this conceptual framework is increasingly understood in developmental terms, with the semi-hierarchical structure of gene regulatory networks (GRNs, used here in a broad sense incorporating not just the genetic circuitry per se but the factors controlling the timing and location of gene expression and repression), the non-linear relation between magnitude of genetic change and the phenotypic results, the evolutionary potential of co-opting existing GRNs, and developmental responsiveness to nongenetic signals (i.e. epigenetics and plasticity), all requiring modification of standard microevolutionary models, and rendering difficult any simple definition of evolutionary novelty. The developmental factors underlying macroevolution create anisotropic probabilities—i.e., an uneven density distribution—of evolutionary change around any given phenotypic starting point, and the potential for coordinated changes among traits that can accommodate change via epigenetic mechanisms. From this standpoint, “punctuated equilibrium” and “phyletic gradualism” simply represent two cells in a matrix of evolutionary models of phenotypic change, and the origin of trends and evolutionary novelty are not simply functions of ecological opportunity. Over long timescales, contingency becomes especially important, and can be viewed in terms of macroevolutionary lags (the temporal separation between the origin of a trait or clade and subsequent diversification); such lags can arise by several mechanisms: as geological or phylogenetic artifacts, or when diversifications require synergistic interactions among traits, or between traits and external events. The temporal and spatial patterns of the origins of evolutionary novelties are a challenge to macroevolutionary theory; individual events can be described retrospectively, but a general model relating development, genetics, and ecology is needed. An accompanying paper (Jablonski in Evol Biol 2017) reviews diversity dynamics and the sorting of variation, with some general conclusions.     

BBMV: an R package for the estimation of macroevolutionary landscapes

The distribution of traits along phylogenies bears signatures of how ecological and evolutionary processes have interacted to influence phenotypic evolution, which can be deciphered using macroevolutionary models. BBMV implements a model for the evolution of continuous characters on phylogenies that generalizes existing ones, like Brownian motion and the Ornstein‐Uhlenbeck model. In this model quantitative characters evolve under both random diffusion and a deterministic force that can be of any possible shape and strength. The model can be used to infer evolutionary scenarios that remained inaccessible so far, like directional trends, disruptive selection, and even bounded evolution. With this new tool at hand, researchers will be able to test complex hypothesis‐driven scenarios regarding trait evolution, but they will also have the possibility to estimate the shape of the adaptive landscapes in which traits evolved. Ultimately, this will provide a way to infer how ecological processes have influenced phenotypic evolution over long timescales. The BBMV package is implemented in the R statistical language and is freely available on the CRAN repository < https://CRAN.R‐project.org/package=BBMV >. All source code can also be found on < https://github.com/fcboucher/BBMV >, along with a detailed tutorial.     

Ecological constraints from incumbent clades drive trait evolution across the tree‐of‐life of freshwater macroinvertebrates

The rates of species and trait diversification vary across the Tree‐of‐Life and over time. Whereas species richness and clade age generally are decoupled, the correlation of accumulated trait diversity of clades (trait disparity) with clade age remains poorly explored. Total trait disparity may be coupled with clade age if the growth of disparity (disparification) within and across clades is continuous with time in an additive niche expansion process (linear‐cumulative model), or alternatively if the rate of trait disparification varies over time and decreases as ecological space becomes gradually saturated (disparity‐dependent model). Using a clock‐calibrated phylogenetic tree for 143 freshwater macroinvertebrate families and richness and trait databases covering > 6400 species, we measured trait disparity in 18 independent clades that successively transitioned to freshwater ecosystems and analyzed its relation with clade age. We found a positive correlation between clade age and total disparity within clades, but no relationship for most individual traits. Traits unique to freshwater lifestyle were highly variable within older clades, while disparity in younger clades shifted towards partially terrestrial lifestyles and saline tolerance to occupy habitats previously inaccessible or underutilized. These results argue that constraints from incumbent lineages limit trait disparity in younger clades that evolved for filling unoccupied regions of the trait space, which suggests that trait disparification may follow a disparity‐dependent model. Overall, we provide an empirical pattern that reveals the potential of the disparity‐dependent model for understanding fundamental processes shaping trait dynamics across the Tree‐of‐Life.     

When should we expect early bursts of trait evolution in comparative data? Predictions from an evolutionary food web model

Conceptual models of adaptive radiation predict that competitive interactions among species will result in an early burst of speciation and trait evolution followed by a slowdown in diversification rates. Empirical studies often show early accumulation of lineages in phylogenetic trees, but usually fail to detect early bursts of phenotypic evolution. We use an evolutionary simulation model to assemble food webs through adaptive radiation, and examine patterns in the resulting phylogenetic trees and species' traits (body size and trophic position). We find that when foraging trade-offs result in food webs where all species occupy integer trophic levels, lineage diversity and trait disparity are concentrated early in the tree, consistent with the early burst model. In contrast, in food webs in which many omnivorous species feed at multiple trophic levels, high levels of turnover of species' identities and traits tend to eliminate the early burst signal. These results suggest testable predictions about how the niche structure of ecological communities may be reflected by macroevolutionary patterns.     

Rates of morphological evolution are correlated with species richness in salamanders

The tempo and mode of species diversification and phenotypic evolution vary widely across the tree of life, yet the relationship between these processes is poorly known. Previous tests of the relationship between rates of phenotypic evolution and rates of species diversification have assumed that species richness increases continuously through time. If this assumption is violated, simple phylogenetic estimates of net diversification rate may bear no relationship to processes that influence the distribution of species richness among clades. Here, we demonstrate that the variation in species richness among plethodontid salamander clades is unlikely to have resulted from simple time-dependent processes, leading to fundamentally different conclusions about the relationship between rates of phenotypic evolution and species diversification. Morphological evolutionary rates of both size and shape evolution are correlated with clade species richness, but are uncorrelated with simple estimators of net diversification that assume constancy of rates through time. This coupling between species diversification and phenotypic evolution is consistent with the hypothesis that clades with high rates of morphological trait evolution may diversify more than clades with low rates. Our results indicate that assumptions about underlying processes of diversity regulation have important consequences for interpreting macroevolutionary patterns.     

How does mutation affect the distribution of phenotypes?

The potential for mutational processes to influence patterns of neutral or adaptive phenotypic evolution is not well understood. If mutations are directionally biased, shifting trait means in a particular direction, or if mutation generates more variance in some directions of multivariate trait space than others, mutation itself might be a source of bias in phenotypic evolution. Here, we use mutagenesis to investigate the affect of mutation on trait mean and (co)variances in zebrafish, Danio rerio. Mutation altered the relationship between age and both prolonged swimming speed and body shape. These observations suggest that mutational effects on ontogeny or aging have the potential to generate variance across the phenome. Mutations had a far greater effect in males than females, although whether this is a reflection of sex‐specific ontogeny or aging remains to be determined. In males, mutations generated positive covariance between swimming speed, size, and body shape suggesting the potential for mutation to affect the evolutionary covariation of these traits. Overall, our observations suggest that mutation does not generate equal variance in all directions of phenotypic space or in each sex, and that pervasive variation in ontogeny or aging within a cohort could affect the variation available to evolution.     

Calls,colours, shape,and genes: a multi‐trait approach to the study of geographic variation in the Amazonian frog Allobates femoralis

Evolutionary divergence in behavioural traits related to mating may represent the initial stage of speciation. Direct selective forces are usually invoked to explain divergence in mate‐recognition traits, often neglecting a role for neutral processes or concomitant differentiation in ecological traits. We adopted a multi‐trait approach to obtain a deeper understanding of the mechanisms behind allopatric divergence in the Amazonian frog, Allobates femoralis. We tested the null hypothesis that geographic distance between populations correlates with genetic and phenotypic divergence, and compared divergence between mate‐recognition (acoustic) and ecological (coloration, body‐shape) traits. We quantified geographic variation in 39 phenotypic traits and a mitochondrial DNA marker among 125 individuals representing eight populations. Geographic variation in acoustic traits was pronounced and tracked the spatial genetic variation, which appeared to be neutral. Thus, the evolution of acoustic traits tracked the shared history of the populations, which is unexpected for pan‐Amazonian taxa or for mate‐recognition traits. Divergence in coloration appeared uncorrelated with genetic distance, and might be partly attributed to local selective pressures, and perhaps to Batesian mimicry. Divergence in body‐shape traits was low. The results obtained depict a complex evolutionary scenario and emphasize the importance of considering multiple traits when disentangling the forces behind allopatric divergence. ©2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 826–838.     

Body mass evolution and diversification within horses (family Equidae)

Horses (family Equidae) are a classic example of adaptive radiation, exhibiting a nearly 60‐fold increase in maximum body mass and a peak taxonomic diversity of nearly 100 species across four continents. Such patterns are commonly attributed to niche competition, in which increased taxonomic diversity drives increased size disparity. However, neutral processes, such as macroevolutionary ‘diffusion’, can produce similar increases in disparity without increased diversity. Using a comprehensive database of Equidae species size estimates and a common mathematical framework, we measure the contributions of diversity‐driven and diffusion‐driven mechanisms for increased disparity during the Equidae radiation. We find that more than 90% of changes in size disparity are attributable to diffusion alone. These results clarify the role of species competition in body size evolution, indicate that morphological disparity and species diversity may be only weakly coupled in general, and demonstrate that large species may evolve from neutral macroevolutionary diffusion processes alone.     

Spatial processes and evolutionary models: a critical review

Evolution is a fundamentally population level process in which variation, drift and selection produce both temporal and spatial patterns of change. Statistical model fitting is now commonly used to estimate which kind of evolutionary process best explains patterns of change through time using models like Brownian motion, stabilizing selection (Ornstein–Uhlenbeck) and directional selection on traits measured from stratigraphic sequences or on phylogenetic trees. But these models assume that the traits possessed by a species are homogeneous. Spatial processes such as dispersal, gene flow and geographical range changes can produce patterns of trait evolution that do not fit the expectations of standard models, even when evolution at the local‐population level is governed by drift or a typical OU model of selection. The basic properties of population level processes (variation, drift, selection and population size) are reviewed and the relationship between their spatial and temporal dynamics is discussed. Typical evolutionary models used in palaeontology incorporate the temporal component of these dynamics, but not the spatial. Range expansions and contractions introduce rate variability into drift processes, range expansion under a drift model can drive directional change in trait evolution, and spatial selection gradients can create spatial variation in traits that can produce long‐term directional trends and punctuation events depending on the balance between selection strength, gene flow, extirpation probability and model of speciation. Using computational modelling that spatial processes can create evolutionary outcomes that depart from basic population‐level notions from these standard macroevolutionary models.     

Adaptive processes drive ecomorphological convergent evolution in antwrens (Thamnophilidae)

Phylogenetic niche conservatism (PNC) and convergence are contrasting evolutionary patterns that describe phenotypic similarity across independent lineages. Assessing whether and how adaptive processes give origin to these patterns represent a fundamental step toward understanding phenotypic evolution. Phylogenetic model‐based approaches offer the opportunity not only to distinguish between PNC and convergence, but also to determine the extent that adaptive processes explain phenotypic similarity. The Myrmotherula complex in the Neotropical family Thamnophilidae is a polyphyletic group of sexually dimorphic small insectivorous forest birds that are relatively homogeneous in size and shape. Here, we integrate a comprehensive species‐level molecular phylogeny of the Myrmotherula complex with morphometric and ecological data within a comparative framework to test whether phenotypic similarity is described by a pattern of PNC or convergence, and to identify evolutionary mechanisms underlying body size and shape evolution. We show that antwrens in the Myrmotherula complex represent distantly related clades that exhibit adaptive convergent evolution in body size and divergent evolution in body shape. Phenotypic similarity in the group is primarily driven by their tendency to converge toward smaller body sizes. Differences in body size and shape across lineages are associated to ecological and behavioral factors.     

Burrowing constrains patterns of skull shape evolution in wrasses

The evolution of behavioral and ecological specialization can have marked effects on the tempo and mode of phenotypic evolution. Head-first burrowing has been shown to exert powerful selective pressures on the head and body shapes of many vertebrate and invertebrate taxa. In wrasses, burrowing behaviors have evolved multiple times independently, and are commonly used in foraging and predator avoidance behaviors. While recent studies have examined the kinematics and body shape morphology associated with this behavior, no study to-date has examined the macroevolutionary implications of burrowing on patterns of phenotypic diversification in this clade. Here, we use three-dimensional geometric morphometrics and phylogenetic comparative methods to study the evolution of skull shape in fossorial wrasses and their relatives. We test for skull shape differences between burrowing and non burrowing wrasses and evaluate hypotheses of shape convergence among the burrowing wrasses. We also quantify rates of skull shape evolution between burrowing and non burrowing wrasses to test for whether burrowing constrains or accelerates rates of skull shape evolution in this clade. We find that while burrowing and non burrowing wrasses exhibit similar degrees of morphological disparity, for burrowing wrasses, it took nearly twice as long to amass this disparity. Furthermore, while the disparities between groups are evenly matched, we find that most burrowing species are confined to a particular region of shape space with most species exhibiting narrower heads than many non-burrowing species. These results suggest head-first burrowing constrains patterns of skull shape diversification in wrasses by potentially restricting the range of phenotypes that can perform this behavior.