Jeremy M. Beaulieu

Three keys to the radiation of angiosperms into freezing environments

Early flowering plants are thought to have been woody species restricted to warm habitats. This lineage has since radiated into almost every climate, with manifold growth forms. As angiosperms spread and climate changed, they evolved mechanisms to cope with episodic freezing. To explore the evolution of traits underpinning the ability to persist in freezing conditions, we assembled a large species-level database of growth habit (woody or herbaceous; 49,064 species), as well as leaf phenology (evergreen or deciduous), diameter of hydraulic conduits (that is, xylem vessels and tracheids) and climate occupancies (exposure to freezing). To model the evolution of species’ traits and climate occupancies, we combined these data with an unparalleled dated molecular phylogeny (32,223 species) for land plants. Here we show that woody clades successfully moved into freezing-prone environments by either possessing transport networks of small safe conduits and/or shutting down hydraulic function by dropping leaves during freezing. Herbaceous species largely avoided freezing periods by senescing cheaply constructed aboveground tissue. Growth habit has long been considered labile, but we find that growth habit was less labile than climate occupancy. Additionally, freezing environments were largely filled by lineages that had already become herbs or, when remaining woody, already had small conduits (that is, the trait evolved before the climate occupancy). By contrast, most deciduous woody lineages had an evolutionary shift to seasonally shedding their leaves only after exposure to freezing (that is, the climate occupancy evolved before the trait). For angiosperms to inhabit novel cold environments they had to gain new structural and functional trait solutions; our results suggest that many of these solutions were probably acquired before their foray into the cold.

Detecting Hidden Diversification Shifts in Models of Trait-Dependent Speciation and Extinction

The distribution of diversity can vary considerably from clade to clade. Attempts to understand these patterns often employ state-dependent speciation and extinction models to determine whether the evolution of a particular novel trait has increased speciation rates and/or decreased extinction rates. It is still unclear, however, whether these models are uncovering important drivers of diversification, or whether they are simply pointing to more complex patterns involving many unmeasured and co-distributed factors. Here we describe an extension to the popular state-dependent speciation and extinction models that specifically accounts for the presence of unmeasured factors that could impact diversification rates estimated for the states of any observed trait, addressing at least one major criticism of BiSSE (Binary State Speciation and Extinction) methods. Specifically, our model, which we refer to as HiSSE (Hidden State Speciation and Extinction), assumes that related to each observed state in the model are hidden states that exhibit potentially distinct diversification dynamics and transition rates than the observed states in isolation. We also demonstrate how our model can be used as character-independent diversification models that allow for a complex diversification process that is independent of the evolution of a character. Under rigorous simulation tests and when applied to empirical data, we find that HiSSE performs reasonably well, and can at least detect net diversification rate differences between observed and hidden states and detect when diversification rate differences do not correlate with the observed states. We discuss the remaining issues with state-dependent speciation and extinction models in general, and the important ways in which HiSSE provides a more nuanced understanding of trait-dependent diversification.

Extinction can be estimated from moderately sized molecular phylogenies

Hundreds of studies have been dedicated to estimating speciation and extinction from phylogenies of extant species. Although it has long been known that estimates of extinction rates using trees of extant organisms are often uncertain, an influential paper by Rabosky (2010) suggested that when birth rates vary continuously across the tree, estimates of the extinction fraction (i.e., extinction rate/speciation rate) will appear strongly bimodal, with a peak suggesting no extinction and a peak implying speciation and extinction rates are approaching equality. On the basis of these results, and the realistic nature of this form of rate variation, it is now generally assumed by many practitioners that extinction cannot be understood from molecular phylogenies alone. Here, we reevaluated and extended the analyses of Rabosky (2010) and come to the opposite conclusion—namely, that it is possible to estimate extinction from molecular phylogenies, even with model violations due to heritable variation in diversification rate. Note that while it may be tempting to interpret our study as advocating the application of simple birth–death models, our goal here is to show how a particular model violation does not necessitate the abandonment of an entire field: use prudent caution, but do not abandon all hope.

Behavioural changes and the adaptive diversification of pigeons and doves.

What factors determine the extent of evolutionary diversification remains a major question in evolutionary biology. Behavioural changes have long been suggested to be a major driver of phenotypic diversification by exposing animals to new selective pressures. Nevertheless, the role of behaviour in evolution remains controversial because behavioural changes can also retard evolutionary change by hiding genetic variation from selection. In the present study, we apply recently implemented Ornstein–Uhlenbeck evolutionary models to show that behavioural changes led to associated evolutionary responses in functionally relevant morphological traits of pigeons and doves (Columbiformes). Specifically, changes from terrestrial to arboreal foraging behaviour reconstructed in a set of phylogenies brought associated shorter tarsi and longer tails, consistent with functional predictions. Interestingly, the transition to arboreality accelerated the rates of evolutionary divergence, leading to an increased morphological specialization that seems to have subsequently constrained reversals to terrestrial foraging. Altogether, our results support the view that behaviour may drive evolutionary diversification, but they also highlight that its evolutionary consequences largely depend on the limits imposed by the functional demands of the adaptive zone.

Explaining the distribution of breeding and dispersal syndromes in conifers

The evolution of plants exhibiting different sexes, or dioecy, is correlated with a number of ecological and life-history traits such as woody growth form and animal-dispersed seeds, but the underlying causes of these associations are unclear. Previous work in seed plants has suggested that the evolution of fleshy cones or seeds may favour dioecy. In this study, we use a well-sampled molecular phylogeny of conifers to show that although dioecy and fleshiness strongly co-occur at the species level, this relationship has not resulted from numerous separate origins of this trait combination or from differential rates of diversification. Instead, we suggest that two character combinations—the ancestral dry-monoecious condition and the derived fleshy-dioecious condition—have persisted in conifers longer than other combinations over evolutionary time. The persistence of these trait combinations appears to reflect differences in the rate of successful transition into and out of these character states over time, as well as the geographical restriction of species with rare combinations and their consequent vulnerability to extinction. In general, we argue that such persistence explanations should be considered alongside ‘key innovation’ hypotheses in explaining the phylogenetic distribution of traits.

The impact of shifts in marine biodiversity hotspots on patterns of range evolution: Evidence from the Holocentridae (squirrelfishes and soldierfishes)

One of the most striking biodiversity patterns is the uneven distribution of marine species richness, with species diversity in the Indo‐Australian Archipelago (IAA) exceeding all other areas. However, the IAA formed fairly recently, and marine biodiversity hotspots have shifted across nearly half the globe since the Paleogene. Understanding how lineages have responded to shifting biodiversity hotspots represents a necessary historic perspective on the formation and maintenance of global marine biodiversity. Such evolutionary inferences are often challenged by a lack of fossil evidence that provide insights into historic patterns of abundance and diversity. The greatest diversity of squirrelfishes and soldierfishes (Holocentridae) is in the IAA, yet these fishes also represent some of the most numerous fossil taxa in deposits of the former West Tethyan biodiversity hotspot. We reconstruct the pattern of holocentrid range evolution using time‐calibrated phylogenies that include most living species and several fossil lineages, demonstrating the importance of including fossil species as terminal taxa in ancestral area reconstructions. Holocentrids exhibit increased range fragmentation following the West Tethyan hotspot collapse. However, rather than originating within the emerging IAA hotspot, the IAA has acted as a reservoir for holocentrid diversity that originated in adjacent regions over deep evolutionary time scales.

Non-null Effects of the Null Range in Biogeographic Models: Exploring Parameter Estimation in the DEC Model

Historical biogeography seeks to understand the distribution of biodiversity in space and time. The dispersal-extinction-cladogenesis (DEC) model, a likelihood-based model of geographic range evolution, is widely used in assessing the biogeography of clades. Robust inference of dispersal and local extinction parameters is crucial for biogeographic inference, and yet a major caveat to its use is that the DEC model severely underestimates local extinction. We suggest that this is mainly due to the way in which the model is constructed to allow observed species to transition into being present in no areas (i.e., null range). By prohibiting transitions into the null range in the transition rate matrix, we were able to better infer local extinction and support this with simulations. This modified model, DEC*, has higher model fit and model adequacy than DEC, suggesting this modification should be considered for DEC and other models of geographic range evolution.

Modelling Stabilizing Selection: The Attraction of Ornstein–Uhlenbeck Models

Ornstein–Uhlenbeck models are a generalization of Brownian motion models that allow trait values to evolve to follow optima. They have become broadly popular in evolutionary studies due to their ability to better fit empirical data as well as for the biological conclusions which can be drawn based on their parameter estimates, especially optimum trait values. We include a survey of available software implementing these models in phylogenetics as well as cautions regarding the use of this software.

Hidden Markov Models for Studying the Evolution of Binary Morphological Characters

Biologists now have the capability of building large phylogenetic trees consisting of tens of thousands of species, from which important comparative questions can be addressed. However, to the extent that biologists have applied these large trees to comparative data, it is clear that current methods, such as those that deal with the evolution of binary morphological characters, make unrealistic assumptions about how these characters are modeled. As phylogenies increase both in size and scope, it is likely that the lability of a binary character will differ significantly among lineages. In this chapter, we describe how a new generalized model, which we refer to as the “hidden rates model” (HRM), can be used to identify different rates of evolution in a discrete binary character along different branches of a phylogeny. The HRM is part of a class of models that are more broadly known as Hidden Markov models because it presupposes that unobserved “hidden” rate classes underlie each observed state and that each rate class represents potentially different transition rates to and from these observed states. As we discuss, the recognition and accommodation of this heterogeneity can provide a robust picture of binary character evolution.