RATES OF MORPHOLOGICAL EVOLUTION ARE CORRELATED WITH SPECIES RICHNESS IN SALAMANDERS

Size is not everything: rates of genome size evolution, not C-value, correlate with speciation in angiosperms

Angiosperms represent one of the key examples of evolutionary success, and their diversity dwarfs other land plants; this success has been linked, in part, to genome size and phenomena such as whole genome duplication events. However, while angiosperms exhibit a remarkable breadth of genome size, evidence linking overall genome size to diversity is equivocal, at best. Here, we show that the rates of speciation and genome size evolution are tightly correlated across land plants, and angiosperms show the highest rates for both, whereas very slow rates are seen in their comparatively species-poor sister group, the gymnosperms. No evidence is found linking overall genome size and rates of speciation. Within angiosperms, both the monocots and eudicots show the highest rates of speciation and genome size evolution, and these data suggest a potential explanation for the megadiversity of angiosperms. It is difficult to associate high rates of diversification with different types of polyploidy, but it is likely that high rates of evolution correlate with a smaller genome size after genome duplications. The diversity of angiosperms may, in part, be due to an ability to increase evolvability by benefiting from whole genome duplications, transposable elements and general genome plasticity.

How Metamorphosis Is Different in Plethodontids: Larval Life History Perspectives on Life-Cycle Evolution

Plethodontid salamanders exhibit biphasic, larval form paedomorphic, and direct developing life cycles. This diversity of developmental strategies exceeds that of any other family of terrestrial vertebrate. Here we compare patterns of larval development among the three divergent lineages of biphasic plethodontids and other salamanders. We discuss how patterns of life-cycle evolution and larval ecology might have produced a wide array of larval life histories. Compared with many other salamanders, most larval plethodontids have relatively slow growth rates and sometimes exceptionally long larval periods (up to 60 mo). Recent phylogenetic analyses of life-cycle evolution indicate that ancestral plethodontids were likely direct developers. If true, then biphasic and paedomorphic lineages might have been independently derived through different developmental mechanisms. Furthermore, biphasic plethodontids largely colonized stream habitats, which tend to have lower productivity than seasonally ephemeral ponds. Consistent with this, plethodontid larvae grow very slowly, and metamorphic timing does not appear to be strongly affected by growth history. On the basis of this, we speculate that feeding schedules and stress hormones might play a comparatively reduced role in governing the timing of metamorphosis of stream-dwelling salamanders, particularly plethodontids.

Rapid postglacial diversification and long-term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.

The impact of rate heterogeneity on inference of phylogenetic models of trait evolution

Rates of trait evolution are known to vary across phylogenies; however, standard evolutionary models assume a homogeneous process of trait change. These simple methods are widely applied in small‐scale phylogenetic studies, whereas models of rate heterogeneity are not, so the prevalence and patterns of potential rate variation in groups up to hundreds of species remain unclear. The extent to which trait evolution is modelled accurately on a given phylogeny is also largely unknown because studies typically lack absolute model fit tests. We investigated these issues by applying both rate‐static and variable‐rates methods on (i) body mass data for 88 avian clades of 10–318 species, and (ii) data simulated under a range of rate‐heterogeneity scenarios. Our results show that rate heterogeneity is present across small‐scaled avian clades, and consequently applying only standard single‐process models prompts inaccurate inferences about the generating evolutionary process. Specifically, these approaches underestimate rate variation, and systematically mislabel temporal trends in trait evolution. Conversely, variable‐rates approaches have superior relative fit (they are the best model) and absolute fit (they describe the data well). We show that rate changes such as single internal branch variations, rate decreases and early bursts are hard to detect, even by variable‐rates models. We also use recently developed absolute adequacy tests to highlight misleading conclusions based on relative fit alone (e.g. a consistent preference for constrained evolution when isolated terminal branch rate increases are present). This work highlights the potential for robust inferences about trait evolution when fitting flexible models in conjunction with tests for absolute model fit.

Analyzing endocrine system conservation and evolution

Analyzing variation in rates of evolution can provide important insights into the factors that constrain trait evolution, as well as those that promote diversification. Metazoan endocrine systems exhibit apparent variation in evolutionary rates of their constituent components at multiple levels, yet relatively few studies have quantified these patterns and analyzed them in a phylogenetic context. This may be in part due to historical and current data limitations for many endocrine components and taxonomic groups. However, recent technological advancements such as high-throughput sequencing provide the opportunity to collect large-scale comparative data sets for even non-model species. Such ventures will produce a fertile data landscape for evolutionary analyses of nucleic acid and amino acid based endocrine components. Here I summarize evolutionary rate analyses that can be applied to categorical and continuous endocrine traits, and also those for nucleic acid and protein-based components. I emphasize analyses that could be used to test whether other variables (e.g., ecology, ontogenetic timing of expression, etc.) are related to patterns of rate variation and endocrine component diversification. The application of phylogenetic-based rate analyses to comparative endocrine data will greatly enhance our understanding of the factors that have shaped endocrine system evolution.

What explains patterns of species richness? The relative importance of climatic-niche evolution, morphological evolution, and ecological limits in salamanders

A major goal of evolutionary biology and ecology is to understand why species richness varies among clades. Previous studies have suggested that variation in richness among clades might be related to variation in rates of morphological evolution among clades (e.g., body size and shape). Other studies have suggested that richness patterns might be related to variation in rates of climatic-niche evolution. However, few studies, if any, have tested the relative importance of these variables in explaining patterns of richness among clades. Here, we test their relative importance among major clades of Plethodontidae, the most species-rich family of salamanders. Earlier studies have suggested that climatic-niche evolution explains patterns of diversification among plethodontid clades, whereas rates of morphological evolution do not. A subsequent study stated that rates of morphological evolution instead explained patterns of species richness among plethodontid clades (along with "ecological limits" on richness of clades, leading to saturation of clades with species, given limited resources). However, they did not consider climatic-niche evolution. Using phylogenetic multiple regression, we show that rates of climatic-niche evolution explain most variation in richness among plethodontid clades, whereas rates of morphological evolution do not. We find little evidence that ecological limits explain patterns of richness among plethodontid clades. We also test whether rates of morphological and climatic-niche evolution are correlated, and find that they are not. Overall, our results help explain richness patterns in a major amphibian group and provide possibly the first test of the relative importance of climatic niches and morphological evolution in explaining diversity patterns.

Evolution of body morphology and beak shape revealed by a morphometric analysis of 14 Paridae species

BACKGROUND

Morphological characters of birds reflect their adaptive evolution and ecological requirements and are also relevant to phylogenetic relationships within a group of related species. The tits (Paridae) are known to be outwardly homogeneous in shape, with one aberrant member, the Ground Tit (Pseudopodoces humilis), which is quite different from its relatives in both body morphology and beak shape. We combined traditional measurements and geometric morphometrics to quantify the variation in body morphology and beak shape of 14 Paridae species distributed in China. Based on these results, we sought to assess the contribution of phylogeny, altitude and species interactions to the evolution of morphological traits.

RESULTS

The basic features for discriminating among the 14 species studied here were overall body size, the ratio of body and tail length to culmen and tarsus length, and beak shape (long/slender/pointy vs. short/robust/blunt). These dimensions clearly separate Ps. humilis and Melanochlora sultanea from the other species in shape space. Body length and PC3 of beak shape (round outline vs. straight outline) show significant phylogenetic signals. Across 14 species, altitude is related to tarsus, culmen length and PC1 of beak shape. Within Parus major, altitude is related to body weight, body length, culmen length and PC1 of body morphology. Morphological distances and geographic distances among species are positively correlated.

CONCLUSIONS

The body morphology of Paridae species shows extensive evolutionary changes, while their beak has mainly evolved along the long/slender/pointy vs. short/robust/blunt dimension. Only body length and beak curvature show a phylogenetic signal. Altitude correlates with multiple traits both across and within species, suggesting that altitude is an important factor in promoting morphological divergence. The deviant appearance of Ps. humilis corresponds to its foraging and feeding adaptations to high-altitude steppe habitats. Our results also show a higher level of morphological divergence with greater difference in distribution ranges among the Paridae species involved in this study.

Diversification rates and phenotypic evolution in venomous snakes (Elapidae)

The relationship between rates of diversification and of body size change (a common proxy for phenotypic evolution) was investigated across Elapidae, the largest radiation of highly venomous snakes. Time-calibrated phylogenetic trees for 175 species of elapids (more than 50% of known taxa) were constructed using seven mitochondrial and nuclear genes. Analyses using these trees revealed no evidence for a link between speciation rates and changes in body size. Two clades (Hydrophis, Micrurus) show anomalously high rates of diversification within Elapidae, yet exhibit rates of body size evolution almost identical to the general elapid 'background' rate. Although correlations between speciation rates and rates of body size change exist in certain groups (e.g. ray-finned fishes, passerine birds), the two processes appear to be uncoupled in elapid snakes. There is also no detectable shift in diversification dynamics associated with the colonization of Australasia, which is surprising given that elapids appear to be the first clade of venomous snakes to reach the continent.

Challenges in the estimation of extinction from molecular phylogenies: A response to Beaulieu and O'Meara

Time-calibrated phylogenies that contain only living species have been widely used to study the dynamics of speciation and extinction. Concerns about the reliability of phylogenetic extinction estimates were raised by Rabosky (2010), where I suggested that unaccommodated heterogeneity in speciation rate could lead to positively biased extinction estimates. In a recent article, Beaulieu and O'Meara (2015a) correctly point out several technical errors in the execution of my 2010 study and concluded that phylogenetic extinction estimates are robust to speciation rate heterogeneity under a range of model parameters. I demonstrate that Beaulieu and O'Meara underestimated the magnitude of speciation rate variation in real phylogenies and consequently did not incorporate biologically meaningful levels of rate heterogeneity into their simulations. Using parameter values drawn from the recent literature, I find that modest levels of heterogeneity in speciation rate result in a consistent, positive bias in extinction estimates that are exacerbated by phylogenetic tree size. This bias, combined with the inherent lack of information about extinction in molecular phylogenies, suggests that extinction rate estimates from phylogenies of extant taxa only should be treated with caution.

The Evolution of Bony Vertebrate Enhancers at Odds with Their Coding Sequence Landscape

Enhancers lie at the heart of transcriptional and developmental gene regulation. Therefore, changes in enhancer sequences usually disrupt the target gene expression and result in disease phenotypes. Despite the well-established role of enhancers in development and disease, evolutionary sequence studies are lacking. The current study attempts to unravel the puzzle of bony vertebrates’ conserved noncoding elements (CNE) enhancer evolution. Bayesian phylogenetics of enhancer sequences spotlights promising interordinal relationships among placental mammals, proposing a closer relationship between humans and laurasiatherians while placing rodents at the basal position. Clock-based estimates of enhancer evolution provided a dynamic picture of interspecific rate changes across the bony vertebrate lineage. Moreover, coelacanth in the study augmented our appreciation of the vertebrate cis-regulatory evolution during water–land transition. Intriguingly, we observed a pronounced upsurge in enhancer evolution in land-dwelling vertebrates. These novel findings triggered us to further investigate the evolutionary trend of coding as well as CNE nonenhancer repertoires, to highlight the relative evolutionary dynamics of diverse genomic landscapes. Surprisingly, the evolutionary rates of enhancer sequences were clearly at odds with those of the coding and the CNE nonenhancer sequences during vertebrate adaptation to land, with land vertebrates exhibiting significantly reduced rates of coding sequence evolution in comparison to their fast evolving regulatory landscape. The observed variation in tetrapod cis-regulatory elements caused the fine-tuning of associated gene regulatory networks. Therefore, the increased evolutionary rate of tetrapods’ enhancer sequences might be responsible for the variation in developmental regulatory circuits during the process of vertebrate adaptation to land.

Nested radiations and the pulse of angiosperm diversification: increased diversification rates often follow whole genome duplications

Our growing understanding of the plant tree of life provides a novel opportunity to uncover the major drivers of angiosperm diversity. Using a time-calibrated phylogeny, we characterized hot and cold spots of lineage diversification across the angiosperm tree of life by modeling evolutionary diversification using stepwise AIC (MEDUSA). We also tested the whole-genome duplication (WGD) radiation lag-time model, which postulates that increases in diversification tend to lag behind established WGD events. Diversification rates have been incredibly heterogeneous throughout the evolutionary history of angiosperms and reveal a pattern of 'nested radiations' - increases in net diversification nested within other radiations. This pattern in turn generates a negative relationship between clade age and diversity across both families and orders. We suggest that stochastically changing diversification rates across the phylogeny explain these patterns. Finally, we demonstrate significant statistical support for the WGD radiation lag-time model. Across angiosperms, nested shifts in diversification led to an overall increasing rate of net diversification and declining relative extinction rates through time. These diversification shifts are only rarely perfectly associated with WGD events, but commonly follow them after a lag period.

Potential rapid evolution of foot morphology in Italian plethodontid salamanders (Hydromantes strinatii) following the colonization of an artificial cave

How organisms respond to environmental change is a long-standing question in evolutionary biology. Species invading novel habitats provide an opportunity to examine contemporary evolution in action and decipher the pace of evolutionary change over short timescales. Here, we characterized phenotypic evolution in the Italian plethodontid salamander, Hydromantes strinatii, following the recent colonization of an artificial cave by a forest floor population. When compared with a nearby and genetically related population in the natural forest floor and a nearby cave population, the artificial cave population displayed significant differences in overall foot shape, with more interdigital webbing relative to the other populations. Further, this population evolved significantly larger feet, which corresponded more closely to those found in other cave populations than to forest floor populations to which the cave population is closely related. Finally, we quantified the rate of evolution for both foot shape and foot area, and found that both traits displayed large and significant evolutionary rates, at levels corresponding to other classic cases of rapid evolution in vertebrates. Together, these findings reveal that the response to novel environmental pressures can be large and rapid and that the anatomical shifts observed in the artificial cave population of H. strinatii may represent a case of rapid evolution in response to novel environmental pressures.

Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation

Several evolutionary theories predict that rates of morphological change should be positively associated with the rate at which new species arise. For example, the theory of punctuated equilibrium proposes that phenotypic change typically occurs in rapid bursts associated with speciation events. However, recent phylogenetic studies have found little evidence linking these processes in nature. Here we demonstrate that rates of species diversification are highly correlated with the rate of body size evolution across the 30,000+ living species of ray-finned fishes that comprise the majority of vertebrate biological diversity. This coupling is a general feature of fish evolution and transcends vast differences in ecology and body-plan organization. Our results may reflect a widespread speciational mode of character change in living fishes. Alternatively, these findings are consistent with the hypothesis that phenotypic 'evolvability'-the capacity of organisms to evolve-shapes the dynamics of speciation through time at the largest phylogenetic scales.

Quantifying and comparing phylogenetic evolutionary rates for shape and other high-dimensional phenotypic data

Many questions in evolutionary biology require the quantification and comparison of rates of phenotypic evolution. Recently, phylogenetic comparative methods have been developed for comparing evolutionary rates on a phylogeny for single, univariate traits (σ(2)), and evolutionary rate matrices (R) for sets of traits treated simultaneously. However, high-dimensional traits like shape remain under-examined with this framework, because methods suited for such data have not been fully developed. In this article, I describe a method to quantify phylogenetic evolutionary rates for high-dimensional multivariate data (σ2 mult), found from the equivalency between statistical methods based on covariance matrices and those based on distance matrices (R-mode and Q-mode methods). I then use simulations to evaluate the statistical performance of hypothesis-testing procedures that compare σ2 mult for two or more groups of species on a phylogeny. Under both isotropic and non-isotropic conditions, and for differing numbers of trait dimensions, the proposed method displays appropriate Type I error and high statistical power for detecting known differences in σ2 mult among groups. In contrast, the Type I error rate of likelihood tests based on the evolutionary rate matrix (R) increases as the number of trait dimensions (p) increases, and becomes unacceptably large when only a few trait dimensions are considered. Further, likelihood tests based on R cannot be computed when the number of trait dimensions equals or exceeds the number of taxa in the phylogeny (i.e., when p ≥ N). These results demonstrate that tests based on σ2 mult provide a useful means of comparing evolutionary rates for high-dimensional data that are otherwise not analytically accessible to methods based on the evolutionary rate matrix. This advance thus expands the phylogenetic comparative toolkit for high-dimensional phenotypic traits like shape. Finally, I illustrate the utility of the new approach by evaluating rates of head shape evolution in a lineage of Plethodon salamanders.

An adaptive radiation of frogs in a southeast Asian island archipelago

Living amphibians exhibit a diversity of ecologies, life histories, and species-rich lineages that offers opportunities for studies of adaptive radiation. We characterize a diverse clade of frogs (Kaloula, Microhylidae) in the Philippine island archipelago as an example of an adaptive radiation into three primary habitat specialists or ecotypes. We use a novel phylogenetic estimate for this clade to evaluate the tempo of lineage accumulation and morphological diversification. Because species-level phylogenetic estimates for Philippine Kaloula are lacking, we employ dense population sampling to determine the appropriate evolutionary lineages for diversification analyses. We explicitly take phylogenetic uncertainty into account when calculating diversification and disparification statistics and fitting models of diversification. Following dispersal to the Philippines from Southeast Asia, Kaloula radiated rapidly into several well-supported clades. Morphological variation within Kaloula is partly explained by ecotype and accumulated at high levels during this radiation, including within ecotypes. We pinpoint an axis of morphospace related directly to climbing and digging behaviors and find patterns of phenotypic evolution suggestive of ecological opportunity with partitioning into distinct habitat specialists. We conclude by discussing the components of phenotypic diversity that are likely important in amphibian adaptive radiations.

Integrated fossil and molecular data reveal the biogeographic diversification of the eastern Asian-eastern North American disjunct hickory genus (Carya Nutt.)

The hickory genus (Carya) contains ca. 17 species distributed in subtropical and tropical regions of eastern Asia and subtropical to temperate regions of eastern North America. Previously, the phylogenetic relationships between eastern Asian and eastern North American species of Carya were not fully confirmed even with an extensive sampling, biogeographic and diversification patterns had thus never been investigated in a phylogenetic context. We sampled 17 species of Carya and 15 species representing all other genera of the Juglandaceae as outgroups, with eight nuclear and plastid loci to reconstruct the phylogeny of Carya. The phylogenetic positions of seven extinct genera of the Juglandaceae were inferred using morphological characters and the molecular phylogeny as a backbone constraint. Divergence times within Carya were estimated with relaxed Bayesian dating. Biogeographic analyses were performed in DIVA and LAGRANGE. Diversification rates were inferred by LASER and APE packages. Our results support two major clades within Carya, corresponding to the lineages of eastern Asia and eastern North America. The split between the two disjunct clades is estimated to be 21.58 (95% HPD 11.07-35.51) Ma. Genus-level DIVA and LAGRANGE analyses incorporating both extant and extinct genera of the Juglandaceae suggested that Carya originated in North America, and migrated to Eurasia during the early Tertiary via the North Atlantic land bridge. Fragmentation of the distribution caused by global cooling in the late Tertiary resulted in the current disjunction. The diversification rate of hickories in eastern North America appeared to be higher than that in eastern Asia, which is ascribed to greater ecological opportunities, key morphological innovations, and polyploidy.

Decoupling of taxonomic diversity and morphological disparity during decline of the Cambrian trilobite family Pterocephaliidae

Although discordance between taxonomic diversity and morphological disparity is common, little is known about the underlying dynamics that drive this decoupling. Early in the history of the Cambrian trilobite family Pterocephaliidae, there was an increase in taxonomic diversity and morphological diversity. As taxonomic diversity declined in the later history of the clade, range of variation stayed high and disparity continued to increase. However, per-branch rates of morphological evolution estimated from a recent phylogeny decreased with time. Neither within-trait nor within-species variation increased or decreased, suggesting that the declining rates of morphological evolution were more likely related to ecological opportunity or niche partitioning, rather than increasing intrinsic constraints. This is further supported by evidence for increased biofacies associations throughout the time period. Thus, the high disparity seen at low taxonomic diversity late in the history of this clade was due to extinction - either random or targeting mean forms - rather than increased rates of morphological evolution. This pattern also provides a scenario that could account for instances of low taxonomic diversity but high morphological disparity in modern groups.

Evidence for determinism in species diversification and contingency in phenotypic evolution during adaptive radiation

Adaptive radiation (AR) theory predicts that groups sharing the same source of ecological opportunity (EO) will experience deterministic species diversification and morphological evolution. Thus, deterministic ecological and morphological evolution should be correlated with deterministic patterns in the tempo and mode of speciation for groups in similar habitats and time periods. We test this hypothesis using well-sampled phylogenies of four squamate groups that colonized the New World (NW) in the Late Oligocene. We use both standard and coalescent models to assess species diversification, as well as likelihood models to examine morphological evolution. All squamate groups show similar early pulses of speciation, as well as diversity-dependent ecological limits on clade size at a continental scale. In contrast, processes of morphological evolution are not easily predictable and do not show similar pulses of early and rapid change. Patterns of morphological and species diversification thus appear uncoupled across these groups. This indicates that the processes that drive diversification and disparification are not mechanistically linked, even among similar groups of taxa experiencing the same sources of EO. It also suggests that processes of phenotypic diversification cannot be predicted solely from the existence of an AR or knowledge of the process of diversification.

Clade age and species richness are decoupled across the eukaryotic tree of life

Explaining the dramatic variation in species richness across the tree of life remains a key challenge in evolutionary biology. At the largest phylogenetic scales, the extreme heterogeneity in species richness observed among different groups of organisms is almost certainly a function of many complex and interdependent factors. However, the most fundamental expectation in macroevolutionary studies is simply that species richness in extant clades should be correlated with clade age: all things being equal, older clades will have had more time for diversity to accumulate than younger clades. Here, we test the relationship between stem clade age and species richness across 1,397 major clades of multicellular eukaryotes that collectively account for more than 1.2 million described species. We find no evidence that clade age predicts species richness at this scale. We demonstrate that this decoupling of age and richness is unlikely to result from variation in net diversification rates among clades. At the largest phylogenetic scales, contemporary patterns of species richness are inconsistent with unbounded diversity increase through time. These results imply that a fundamentally different interpretative paradigm may be needed in the study of phylogenetic diversity patterns in many groups of organisms.

Species richness varies by many orders of magnitude across the evolutionary "tree of life." Some groups, like beetles and flowering plants, contain nearly incomprehensible species diversity, but the overwhelming majority of groups contain far fewer species. Many processes presumably contribute to this variation in diversity, but the most general explanatory variable is the evolutionary age of each group: older groups will simply have had more time for diversity to accumulate than younger groups. We tested whether evolutionary age explains differences in species richness by compiling diversity and age estimates for nearly 1,400 groups of multicellular organisms. Surprisingly, we find no evidence that old groups have more species than young groups. This result appears to hold across the entire tree of life, for taxa as diverse as ferns, fungi, and flies. We demonstrate that this pattern is highly unlikely under simple but widely used evolutionary models that allow diversity to increase through time without bounds. Paleontologists have long contended that diversity-dependent processes have regulated species richness through time, and our results suggest that such processes have left a footprint on the living biota that can even be seen without data from the fossil record.