Increasing the impact of vertebrate scientific collections through 3D imaging: The openVertebrate (oVert) Thematic Collections Network

The impact of preserved museum specimens is transforming and increasing by three-dimensional (3D) imaging that creates high-fidelity online digital specimens. Through examples from the openVertebrate (oVert) Thematic Collections Network, we describe how we created a digitization community dedicated to the shared vision of making 3D data of specimens available and the impact of these data on a broad audience of scientists, students, teachers, artists, and more. High-fidelity digital 3D models allow people from multiple communities to simultaneously access and use scientific specimens. Based on our multiyear, multi-institution project, we identify significant technological and social hurdles that remain for fully realizing the potential impact of digital 3D specimens.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Specimen collection is essential for modern science

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

Ecological correlates of cranial evolution in the megaradiation of dipsadine snakes

BACKGROUND

Dipsadine snakes represent one of the most spectacular vertebrate radiations that have occurred in any continental setting, with over 800 species in South and Central America. Their species richness is paralleled by stunning ecological diversity, ranging from arboreal snail-eating and aquatic eel-eating specialists to terrestrial generalists. Despite the ecological importance of this clade, little is known about the extent to which ecological specialization shapes broader patterns of phenotypic diversity within the group. Here, we test how habitat use and diet have influenced morphological diversification in skull shape across 160 dipsadine species using micro-CT and 3-D geometric morphometrics, and we use a phylogenetic comparative approach to test the contributions of habitat use and diet composition to variation in skull shape among species.

RESULTS

We demonstrate that while both habitat use and diet are significant predictors of shape in many regions of the skull, habitat use significantly predicts shape in a greater number of skull regions when compared to diet. We also find that across ecological groupings, fossorial and aquatic behaviors result in the strongest deviations in morphospace for several skull regions. We use simulations to address the robustness of our results and describe statistical anomalies that can arise from the application of phylogenetic generalized least squares to complex shape data.

CONCLUSIONS

Both habitat and dietary ecology are significantly correlated with skull shape in dipsadines; the strongest relationships involved skull shape in snakes with aquatic and fossorial lifestyles. This association between skull morphology and multiple ecological axes is consistent with a classic model of adaptive radiation and suggests that ecological factors were an important component in driving morphological diversification in the dipsadine megaradiation.

Diversification dynamics in Caribbean rain frogs (Eleutherodactylus) are uncoupled from the anuran community and consistent with adaptive radiation

Adaptive radiation is proposed to play a key role in generating differences in species richness among lineages and geographical regions. Due to the importance of ecological divergence in adaptive radiation, species richness is predicted to be influenced by equilibrium diversity dynamics, although the concept continues to generate much debate. An additional important question is whether radiating clades have intrinsic biological characteristics that make them particularly prone to diversify. We tackle these questions by analysing (i) the temporal patterns of diversification of Caribbean Eleutherodactylus frogs, and (ii) assembly of the complete native anuran community of the Caribbean archipelago (197 species), testing for the presence of equilibrium dynamics and whether diversification patterns of Eleutherodactylus differ from those of the rest of the Caribbean anurans. Diversification rates follow the predicted pattern of rapid diversification early in the radiation which gradually decreases towards the present. Eleutherodactylus diversification is significantly faster than that of the Caribbean anuran community, and although equilibrium dynamics influence richness of all Caribbean anurans, Eleutherodactylus shows higher carrying capacity. Our results indicate that ecological opportunity per se is not sufficient for adaptive radiation and that diverse lineages present intrinsic characteristics that enable them to make the most of available opportunity.

Evolutionary time and species diversity in aquatic ecosystems worldwide

The latitudinal diversity gradient (LDG) is frequently described as the most dramatic biodiversity pattern on Earth, yet ecologists and biogeographers have failed to reach consensus on its primary cause. A key problem in explaining the LDG involves collinearity between multiple factors that are predicted to affect species richness in the same direction. In terrestrial systems, energy input, geographic area, and evolutionary time for species accumulation tend to covary positively with species richness at the largest spatial scales, such that their individual contributions to the LDG are confounded in global analyses. I review three diversity patterns from marine and freshwater systems that break this collinearity and which may thus provide stronger tests of the influence of time on global richness gradients. Specifically, I contrast biodiversity patterns along oceanic depth gradients, in geologically young versus ancient lakes, and in the north versus south polar marine biomes. I focus primarily on fishes due to greater data availability but synthesize patterns for invertebrates where possible. I find that regional-to-global species richness generally declines with depth in the oceans, despite the great age and stability of the deep-sea biome. Geologically ancient lakes generally do not contain more species than young lakes, and the Antarctic marine biome is not appreciably more species rich than the much younger Arctic marine biome. However, endemism is consistently higher in older systems. Patterns for invertebrate groups are less clear than for fishes and reflect a critical need for primary biodiversity data. In summary, the available data suggest that species richness is either decoupled from or only weakly related to the amount of time for diversification. These results suggest that energy, productivity, or geographic area are the primary drivers of large-scale diversity gradients. To the extent that marine and terrestrial diversity gradients result from similar processes, these examples provide evidence against a primary role for evolutionary time as the cause of the LDG.

Subspecies at crossroads: the evolutionary significance of genomic and phenotypic variation in a wide-ranging Australian lizard (Ctenotus pantherinus)

Many subspecies were described to capture phenotypic variation in wide-ranging taxa, with some later being found to correspond to divergent genetic lineages. We investigate whether currently recognized subspecies correspond to distinctive and coherent evolutionary lineages in the widespread Australian lizard Ctenotus pantherinus based on morphological, mitochondrial and genome-wide nuclear variation. We find weak and inconsistent correspondence between morphological patterns and the presumed subspecies ranges, with character polymorphism within regions and broad morphological overlap across regions. Phylogenetic analyses suggest paraphyly of populations assignable to each subspecies, mitonuclear discordance and little congruence between subspecies ranges and the distribution of inferred clades. Genotypic clustering supports admixture across regions. These results undermine the presumed phenotypic and genotypic coherence and distinctiveness of C. pantherinus subspecies. Based on our findings, we comment on the operational and conceptual shortcomings of morphologically defined subspecies and discuss practical challenges in applying the general notion of subspecies as incompletely separated population lineages. We conclude by highlighting a historical asymmetry that has implications for ecology, evolution and conservation: subspecies proposed in the past are difficult to falsify even in the face of new data that challenge their coherence and distinctiveness, whereas modern researchers appear hesitant to propose new subspecies.

Assembling microbial communities: a genomic analysis of a natural experiment in neotropical bamboo internodes

Microbes participate in ecological communities, much like multicellular organisms. However, microbial communities lack the centuries of observation and theory describing and predicting ecological processes available for multicellular organisms. Here, we examine early bacterial community assembly in the water-filled internodes of Amazonian bamboos from the genus Guadua. Bamboo stands form distinct habitat patches within the lowland Amazonian rainforest and provide habitat for a suite of vertebrate and invertebrate species. Guadua bamboos develop sealed, water-filled internodes as they grow. Internodes are presumed sterile or near sterile while closed, but most are eventually opened to the environment by animals, after which they are colonized by microbes. We find that microbial community diversity increases sharply over the first few days of environmental exposure, and taxonomic identity of the microbes changes through this time period as is predicted for early community assembly in macroscopic communities. Microbial community taxonomic turnover is consistent at the bacteria phylum level, but at the level of Operational Taxonomic Units (OTUs), internode communities become increasingly differentiated through time. We argue that these tropical bamboos form an ideal study system for microbial community ecology due to their near-sterile condition prior to opening, relatively consistent environment after opening, and functionally limitless possibilities for replicates. Given the possible importance of opened internode habitats as locations of transmission for both pathogenic and beneficial microbes among animals, understanding the microbial dynamics of the internode habitat is a key conservation concern for the insect and amphibian species that use this microhabitat.

Genetic variability and the ecology of geographic range: A test of the central-marginal hypothesis in Australian scincid lizards

For many species, both local abundance and regional occupancy are highest near the centre of their geographic distributions. One hypothesis for this pattern is that niche suitability declines with increasing distance from a species geographic centre, such that populations near range margins are characterized by reduced density and increased patchiness. In these smaller edge populations, genetic drift is more powerful, leading to the loss of genetic diversity. This simple verbal model has been formalized as the central-marginal hypothesis, which predicts that core populations should have greater genetic diversity than edge populations. Here, we tested the central-marginal hypothesis using a genomic data set of 25 species-level taxa of Australian scincid lizards in the genera Ctenotus and Lerista. A majority of taxa in our data set showed range-wide patterns of genetic variation consistent with central-marginal hypothesis, and eight of 25 taxa showed significantly greater genetic diversity in the centre of their range. We then explored biological, historical, and methodological factors that might predict which taxa support the central-marginal hypothesis. We found that taxa with the strongest evidence for range expansion were the least likely to follow predictions of the central-marginal hypothesis. The majority of these taxa had range expansions that originated at the range edge, which led to a gradient of decreasing genetic diversity from the range edge to the core, contrary to the central-marginal hypothesis.

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

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

No link between population isolation and speciation rate in squamate reptiles

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Biodiversity across space and time in the fossil record

The fossil record is the primary source of information on how biodiversity has varied in deep time, providing unique insight on the long-term dynamics of diversification and their drivers. However, interpretations of fossil record diversity patterns have been much debated, with a traditional focus on global diversity through time. Problems arise because the fossil record is spatially and temporally patchy, so 'global' diversity estimates actually represent the summed diversity across a set of geographically and environmentally distinct regions that vary substantially in number and identity through time. Furthermore, a focus on global diversity lumps the signal of ecological drivers at local and regional scales with the signal of global-scale processes, including variation in the distribution of environments and in provincialism (the extent of subdivision into distinct biogeographic regions). These signals cannot be untangled by studying global diversity measures alone. These conceptual and empirical concerns necessitate a shift away from the study of 'biodiversity through time' and towards the study of 'biodiversity across time and space'. Spatially explicit investigations, including analyses of local- and regional-scale datasets, are central to achieving this and allow analysis of geographic scale, location and the environmental parameters directly experienced by organisms. So far, research in this area has revealed the stability of species richness variation among environments through time, and the potential climatic and Earth-system drivers of changing biodiversity. Ultimately, this research program promises to address key questions regarding the assembly of biodiversity, and the contributions of local-, regional- and global-scale processes to the diversification of life on Earth.

Macroevolutionary thermodynamics: Temperature and the tempo of evolution in the tropics

An influential hypothesis proposes that the tempo of evolution is faster in the tropics. Emerging evidence, including a study in this issue of PLOS Biology, challenges this view, raising new questions about the causes of Earth’s iconic latitudinal diversity gradient (LDG).

A return-on-investment approach for prioritization of rigorous taxonomic research needed to inform responses to the biodiversity crisis

Global biodiversity loss is a profound consequence of human activity. Disturbingly, biodiversity loss is greater than realized because of the unknown number of undocumented species. Conservation fundamentally relies on taxonomic recognition of species, but only a fraction of biodiversity is described. Here, we provide a new quantitative approach for prioritizing rigorous taxonomic research for conservation. We implement this approach in a highly diverse vertebrate group-Australian lizards and snakes. Of 870 species assessed, we identified 282 (32.4%) with taxonomic uncertainty, of which 17.6% likely comprise undescribed species of conservation concern. We identify 24 species in need of immediate taxonomic attention to facilitate conservation. Using a broadly applicable return-on-investment framework, we demonstrate the importance of prioritizing the fundamental work of identifying species before they are lost.

A test for rate-coupling of trophic and cranial evolutionary dynamics in New World bats

Morphological evolution is often assumed to be causally related to underlying patterns of ecological trait evolution. However, few studies have directly tested whether evolutionary dynamics of-and major shifts in-ecological resource use are coupled with morphological shifts that may facilitate trophic innovation. Using diet and multivariate cranial (microCT) data, we tested whether rates of trophic and cranial evolution are coupled in the radiation of New World bats. We developed a generalizable information-theoretic method for describing evolutionary rate heterogeneity across large candidate sets of multirate evolutionary models, without relying on a single best-fitting model. We found considerable variation in trophic evolutionary dynamics, in sharp contrast to a largely homogeneous cranial evolutionary process. This dichotomy is surprising given established functional associations between overall skull morphology and trophic ecology. We suggest that assigning discrete trophic states may underestimate trophic generalism and opportunism, and that this radiation could be characterized by labile crania and a homogeneous dynamic of generally high morphological rates. Overall, we discuss how trophic classifications could substantively impact our interpretation of how these dynamics covary in adaptive radiations.

Detecting Lineage-Specific Shifts in Diversification: A Proper Likelihood Approach

The branching patterns of molecular phylogenies are generally assumed to contain information on rates of the underlying speciation and extinction processes. Simple birth-death models with constant, time-varying, or diversity-dependent rates have been invoked to explain these patterns. They have one assumption in common: all lineages have the same set of diversification rates at a given point in time. It seems likely, however, that there is variability in diversification rates across subclades in a phylogenetic tree. This has inspired the construction of models that allow multiple rate regimes across the phylogeny, with instantaneous shifts between these regimes. Several methods exist for calculating the likelihood of a phylogeny under a specified mapping of diversification regimes and for performing inference on the most likely diversification history that gave rise to a particular phylogenetic tree. Here, we show that the likelihood computation of these methods is not correct. We provide a new framework to compute the likelihood correctly and show, with simulations of a single shift, that the correct likelihood indeed leads to parameter estimates that are on average in much better agreement with the generating parameters than the incorrect likelihood. Moreover, we show that our corrected likelihood can be extended to multiple rate shifts in time-dependent and diversity-dependent models. We argue that identifying shifts in diversification rates is a nontrivial model selection exercise where one has to choose whether shifts in now-extinct lineages are taken into account or not. Hence, our framework also resolves the recent debate on such unobserved shifts. [Diversification; macroevolution; phylogeny; speciation].

Estimating Diversification Rates on Incompletely Sampled Phylogenies: Theoretical Concerns and Practical Solutions

Molecular phylogenies are a key source of information about the tempo and mode of species diversification. However, most empirical phylogenies do not contain representatives of all species, such that diversification rates are typically estimated from incompletely sampled data. Most researchers recognize that incomplete sampling can lead to biased rate estimates, but the statistical properties of methods for accommodating incomplete sampling remain poorly known. In this point of view, we demonstrate theoretical concerns with the widespread use of analytical sampling corrections for sparsely sampled phylogenies of higher taxonomic groups. In particular, corrections based on "sampling fractions" can lead to low statistical power to infer rate variation when it is present, depending on the likelihood function used for inference. In the extreme, the sampling fraction correction can lead to spurious patterns of diversification that are driven solely by unbalanced sampling across the tree in concert with low overall power to infer shifts. Stochastic polytomy resolution provides an alternative to sampling fraction approaches that avoids some of these biases. We show that stochastic polytomy resolvers can greatly improve the power of common analyses to estimate shifts in diversification rates. We introduce a new stochastic polytomy resolution method (Taxonomic Addition for Complete Trees [TACT]) that uses birth-death-sampling estimators across an ultrametric phylogeny to estimate branching times for unsampled taxa, with taxonomic information to compatibly place new taxa onto a backbone phylogeny. We close with practical recommendations for diversification inference under several common scenarios of incomplete sampling. [Birth-death process; diversification; incomplete sampling; phylogenetic uncertainty; rate heterogeneity; rate shifts; stochastic polytomy resolution.].

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

What makes a fang? Phylogenetic and ecological controls on tooth evolution in rear-fanged snakes

BACKGROUND

Fangs are a putative key innovation that revolutionized prey capture and feeding in snakes, and - along with their associated venom phenotypes - have made snakes perhaps the most medically-significant vertebrate animals. Three snake clades are known for their forward-positioned fangs, and these clades (Elapidae, Viperidae, and Atractaspidinae) contain the majority of snakes that are traditionally considered venomous. However, many other snakes are "rear-fanged": they possess potentially venom-delivering teeth situated at the rear end of the upper jaw. Quantification of fang phenotypes - and especially those of rear-fanged species - has proved challenging or impossible owing to the small size and relative rarity of many such snakes. Consequently, it has been difficult to understand the evolutionary history of both venom and prey-capture strategies across extant snakes. We quantified variation in the dentition of 145 colubriform ("advanced") snake species using microCT scanning and compared dental characters with ecological data on species' diet and prey capture method(s) to understand broader patterns in snake fang evolution.

RESULTS

Dental traits such as maxilla length, tooth number, and fang size show strong phylogenetic signal across Colubriformes. We find extreme heterogeneity and evolutionary lability in the rear-fanged phenotype in colubrid (colubrine, dipsadine, and natricine lineages) and lamprophiid snakes, in contrast to relative uniformity in the front fanged phenotypes of other groups (vipers and, to a lesser extent, elapids). Fang size and position are correlated with venom-use in vipers, elapids, and colubrid snakes, with the latter group shifting fangs anteriorly by shortening the entire maxillary bone. We find that maxilla length and tooth number may also be correlated with the evolution of dietary specialization. Finally, an ancestral state reconstruction suggests that fang loss is a recurring phenomenon in colubrid snakes, likely accompanied by shifts in diet and prey capture mode.

CONCLUSIONS

Our study provides a framework for quantifying the complex morphologies associated with venom use in snakes. Our results suggest that fang phenotypes, and particularly the rear-fanged phenotype, in snakes are both diverse and labile, facilitating a wide range of ecological strategies and contributing to spectacular radiations of these organisms in tropical and subtropical biomes worldwide.

Complex Ecological Phenotypes on Phylogenetic Trees: A Markov Process Model for Comparative Analysis of Multivariate Count Data

The evolutionary dynamics of complex ecological traits—including multistate representations of diet, habitat, and behavior—remain poorly understood. Reconstructing the tempo, mode, and historical sequence of transitions involving such traits poses many challenges for comparative biologists, owing to their multidimensional nature. Continuous-time Markov chains are commonly used to model ecological niche evolution on phylogenetic trees but are limited by the assumption that taxa are monomorphic and that states are univariate categorical variables. A necessary first step in the analysis of many complex traits is therefore to categorize species into a predetermined number of univariate ecological states, but this procedure can lead to distortion and loss of information. This approach also confounds interpretation of state assignments with effects of sampling variation because it does not directly incorporate empirical observations for individual species into the statistical inference model. In this study, we develop a Dirichlet-multinomial framework to model resource use evolution on phylogenetic trees. Our approach is expressly designed to model ecological traits that are multidimensional and to account for uncertainty in state assignments of terminal taxa arising from effects of sampling variation. The method uses multivariate count data across a set of discrete resource categories sampled for individual species to simultaneously infer the number of ecological states, the proportional utilization of different resources by different states, and the phylogenetic distribution of ecological states among living species and their ancestors. The method is general and may be applied to any data expressible as a set of observational counts from different categories. [Comparative methods; Dirichlet multinomial; ecological niche evolution; macroevolution; Markov model.]

Beyond Reproductive Isolation: Demographic Controls on the Speciation Process

Studies of speciation typically investigate the evolution of reproductive isolation between populations, but several other processes can serve as key steps limiting the formation of species. In particular, the probability of successful speciation can be influenced by factors that affect the frequency with which population isolates form as well as their persistence through time. We suggest that population isolation and persistence have an inherently spatial dimension that can be profitably studied using a conceptual framework drawn from metapopulation ecology. We discuss models of speciation that incorporate demographic processes and highlight the need for a broader application of phylogenetic comparative approaches to evaluate the general importance of population isolation, persistence, and reproductive isolation in speciation. We review diverse and nontraditional data sources that can be leveraged to study isolation and persistence in a comparative framework. This incorporation of spatial demographic information facilitates the integration of perspectives on speciation across disciplines and timescales.

Metabolically similar cohorts of bacteria exhibit strong cooccurrence patterns with diet items and eukaryotic microbes in lizard guts

Gut microbiomes perform essential services for their hosts, including helping them to digest food and manage pathogens and parasites. Performing these services requires a diverse and constantly changing set of metabolic functions from the bacteria in the microbiome. The metabolic repertoire of the microbiome is ultimately dependent on the outcomes of the ecological interactions of its member microbes, as these interactions in part determine the taxonomic composition of the microbiome. The ecological processes that underpin the microbiome's ability to handle a variety of metabolic challenges might involve rapid turnover of the gut microbiome in response to new metabolic challenges, or it might entail maintaining sufficient diversity in the microbiome that any new metabolic demands can be met from an existing set of bacteria. To differentiate between these scenarios, we examine the gut bacteria and resident eukaryotes of two generalist-insectivore lizards, while simultaneously identifying the arthropod prey each lizard was digesting at the time of sampling. We find that the cohorts of bacteria that occur significantly more or less often than expected with arthropod diet items or eukaryotes include bacterial species that are highly similar to each other metabolically. This pattern in the bacterial microbiome could represent an early step in the taxonomic shifts in bacterial microbiome that occur when host lineages change their diet niche over evolutionary timescales.

Thermal physiological traits in tropical lowland amphibians: Vulnerability to climate warming and cooling

Climate change is affecting biodiversity and ecosystem function worldwide, and the lowland tropics are of special concern because organisms living in this region experience temperatures that are close to their upper thermal limits. However, it remains unclear how and whether tropical lowland species will be able to cope with the predicted pace of climate warming. Additionally, there is growing interest in examining how quickly thermal physiological traits have evolved across taxa, and whether thermal physiological traits are evolutionarily conserved or labile. We measured critical thermal maximum (CTmax) and minimum (CTmin) in 56 species of lowland Amazonian frogs to determine the extent of phylogenetic conservatism in tolerance to heat and cold, and to predict species' vulnerability to climate change. The species we studied live in sympatry and represent ~65% of the known alpha diversity at our study site. Given that critical thermal limits may have evolved differently in response to different temperature constraints, we tested whether CTmax and CTmin exhibit different rates of evolutionary change. Measuring both critical thermal traits allowed us to estimate species' thermal breadth and infer their potential to respond to abrupt changes in temperature (warming and cooling). Additionally, we assessed the contribution of life history traits and found that both critical thermal traits were correlated with species' body size and microhabitat use. Specifically, small direct-developing frogs in the Strabomantidae family appear to be at highest risk of thermal stress while tree frogs (Hylidae) and narrow mouthed frogs (Microhylidae) tolerate higher temperatures. While CTmax and CTmin had considerable variation within and among families, both critical thermal traits exhibited similar rates of evolutionary change. Our results suggest that 4% of lowland rainforest frogs assessed will experience temperatures exceeding their CTmax, 25% might be moderately affected and 70% are unlikely to experience pronounced heat stress under a hypothetical 3°C temperature increase.

Is genomic diversity a useful proxy for census population size? Evidence from a species-rich community of desert lizards

Species abundance data are critical for testing ecological theory, but obtaining accurate empirical estimates for many taxa is challenging. Proxies for species abundance can help researchers circumvent time and cost constraints that are prohibitive for long-term sampling. Under simple demographic models, genetic diversity is expected to correlate with census size, such that genome-wide heterozygosity may provide a surrogate measure of species abundance. We tested whether nucleotide diversity is correlated with long-term estimates of abundance, occupancy and degree of ecological specialization in a diverse lizard community from arid Australia. Using targeted sequence capture, we obtained estimates of genomic diversity from 30 species of lizards, recovering an average of 5,066 loci covering 3.6 Mb of DNA sequence per individual. We compared measures of individual heterozygosity to a metric of habitat specialization to investigate whether ecological preference exerts a measurable effect on genetic diversity. We find that heterozygosity is significantly correlated with species abundance and occupancy, but not habitat specialization. Demonstrating the power of genomic sampling, the correlation between heterozygosity and abundance/occupancy emerged from considering just one or two individuals per species. However, genetic diversity does no better at predicting abundance than a single day of traditional sampling in this community. We conclude that genetic diversity is a useful proxy for regional-scale species abundance and occupancy, but a large amount of unexplained variation in heterozygosity suggests additional constraints or a failure of ecological sampling to adequately capture variation in true population size.

Inferring Diversification Rate Variation From Phylogenies With Fossils

Time-calibrated phylogenies of living species have been widely used to study the tempo and mode of species diversification. However, it is increasingly clear that inferences about species diversification-extinction rates in particular-can be unreliable in the absence of paleontological data. We introduce a general framework based on the fossilized birth-death process for studying speciation-extinction dynamics on phylogenies of extant and extinct species. The model assumes that phylogenies can be modeled as a mixture of distinct evolutionary rate regimes and that a hierarchical Poisson process governs the number of such rate regimes across a tree. We implemented the model in BAMM, a computational framework that uses reversible jump Markov chain Monte Carlo to simulate a posterior distribution of macroevolutionary rate regimes conditional on the branching times and topology of a phylogeny. The implementation, we describe can be applied to paleontological phylogenies, neontological phylogenies, and to phylogenies that include both extant and extinct taxa. We evaluate performance of the model on data sets simulated under a range of diversification scenarios. We find that speciation rates are reliably inferred in the absence of paleontological data. However, the inclusion of fossil observations substantially increases the accuracy of extinction rate estimates. We demonstrate that inferences are relatively robust to at least some violations of model assumptions, including heterogeneity in preservation rates and misspecification of the number of occurrences in paleontological data sets.

Phylogenetic tests for evolutionary innovation: the problematic link between key innovations and exceptional diversification

Evolutionary innovation contributes to the spectacular diversity of species and phenotypes across the tree of life. 'Key innovations' are widely operationalized within evolutionary biology as traits that facilitate increased diversification rates, such that lineages bearing the traits ultimately contain more species than closely related lineages lacking the focal trait. In this article, I briefly review the inference, analysis and interpretation of evolutionary innovation on phylogenetic trees. I argue that differential rates of lineage diversification should not be used as the basis for key innovation tests, despite the statistical tractability of such approaches. Under traditional interpretations of the macroevolutionary 'adaptive zone', we should not necessarily expect key innovations to confer faster diversification rates upon lineages that possess them relative to their extant sister clades. I suggest that a key innovation is a trait that allows a lineage to interact with the environment in a fundamentally different way and which, as a result, increases the total diversification-but not necessarily the diversification rate-of the parent clade. Considered alone, branching patterns in phylogenetic trees are poorly suited to the inference of evolutionary innovation due to their inherently low information content with respect to the processes that produce them. However, phylogenies may be important for identifying transformational shifts in ecological and morphological space that are characteristic of innovation at the macroevolutionary scale.This article is part of the themed issue 'Process and pattern in innovations from cells to societies'.

Evolutionary radiation of earless frogs in the Andes: molecular phylogenetics and habitat shifts in high-elevation terrestrial breeding frogs

The loss of hearing structures and loss of advertisement calls in many terrestrial breeding frogs (Strabomantidae) living at high elevations in South America are common and intriguing phenomena. The Andean frog genus Phrynopus Peters, 1873 has undergone an evolutionary radiation in which most species lack the tympanic membrane and tympanic annulus, yet the phylogenetic relationships among species in this group remain largely unknown. Here, we present an expanded molecular phylogeny of Phrynopus that includes 24 nominal species. Our phylogeny includes Phrynopus peruanus, the type species of the genus, and 10 other species for which genetic data were previously unavailable. We found strong support for monophyly of Phrynopus, and that two nominal species-Phrynopus curator and Phrynopus nicoleae-are junior synonyms of Phrynopus tribulosus. Using X-ray computed tomography (CT) imaging, we demonstrate that the absence of external hearing structures is associated with complete loss of the auditory skeletal elements (columella) in at least one member of the genus. We mapped the tympanum condition on to a species tree to infer whether the loss of hearing structures took place once or multiple times. We also assessed whether tympanum condition, body size, and body shape are associated with the elevational distribution and habitat use. We identified a single evolutionary transition that involved the loss of both the tympanic membrane and tympanic annulus, which in turn is correlated with the absence of advertisement calls. We also identified several species pairs where one species inhabits the Andean grassland and the other montane forest. When accounting for phylogenetic relatedness among species, we detected a significant pattern of increasing body size with increasing elevation. Additionally, species at higher elevations tend to develop shorter limbs, shorter head, and shorter snout than species living at lower elevations. Our findings strongly suggest a link between ecological divergence and morphological diversity of terrestrial breeding frogs living in montane gradients.

Genetic diversity is largely unpredictable but scales with museum occurrences in a species-rich clade of Australian lizards

Genetic diversity is a fundamental characteristic of species and is affected by many factors, including mutation rate, population size, life history and demography. To better understand the processes that influence levels of genetic diversity across taxa, we collected genome-wide restriction-associated DNA data from more than 500 individuals spanning 76 nominal species of Australian scincid lizards in the genus Ctenotus To avoid potential biases associated with variation in taxonomic practice across the group, we used coalescent-based species delimitation to delineate 83 species-level lineages within the genus for downstream analyses. We then used these genetic data to infer levels of within-population genetic diversity. Using a phylogenetically informed approach, we tested whether variation in genetic diversity could be explained by population size, environmental heterogeneity or historical demography. We find that the strongest predictor of genetic diversity is a novel proxy for census population size: the number of vouchered occurrences in museum databases. However, museum occurrences only explain a limited proportion of the variance in genetic diversity, suggesting that genetic diversity might be difficult to predict at shallower phylogenetic scales.

Do Macrophylogenies Yield Stable Macroevolutionary Inferences? An Example from Squamate Reptiles

Advances in the generation, retrieval, and analysis of phylogenetic data have enabled researchers to create phylogenies that contain many thousands of taxa. These "macrophylogenies"-large trees that typically derive from megaphylogeny, supermatrix, or supertree approaches-provide researchers with an unprecedented ability to conduct evolutionary analyses across broad phylogenetic scales. Many studies have now used these phylogenies to explore the dynamics of speciation, extinction, and phenotypic evolution across large swaths of the tree of life. These trees are characterized by substantial phylogenetic uncertainty on multiple levels, and the stability of macroevolutionary inferences from these data sets has not been rigorously explored. As a case study, we tested whether five recently published phylogenies for squamate reptiles-each consisting of more than 4000 species-yield congruent inferences about the processes that underlie variation in species richness across replicate evolutionary radiations of Australian snakes and lizards. We find discordance across the five focal phylogenies with respect to clade age and several diversification rate metrics, and in the effects of clade age on species richness. We also find that crown clade ages reported in the literature on these Australian groups are in conflict with all of the large phylogenies examined. Macrophylogenies offer an unprecedented opportunity to address evolutionary and ecological questions at broad phylogenetic scales, but accurately representing the uncertainty that is inherent to such analyses remains a critical challenge to our field. [Australia; macroevolution; macrophylogeny; squamates; time calibration.].

Is BAMM Flawed? Theoretical and Practical Concerns in the Analysis of Multi-Rate Diversification Models

Bayesian analysis of macroevolutionary mixtures (BAMM) is a statistical framework that uses reversible jump Markov chain Monte Carlo to infer complex macroevolutionary dynamics of diversification and phenotypic evolution on phylogenetic trees. A recent article by Moore et al. (MEA) reported a number of theoretical and practical concerns with BAMM. Major claims from MEA are that (i) BAMM’s likelihood function is incorrect, because it does not account for unobserved rate shifts; (ii) the posterior distribution on the number of rate shifts is overly sensitive to the prior; and (iii) diversification rate estimates from BAMM are unreliable. Here, we show that these and other conclusions from MEA are generally incorrect or unjustified. We first demonstrate that MEA’s numerical assessment of the BAMM likelihood is compromised by their use of an invalid likelihood function. We then show that “unobserved rate shifts” appear to be irrelevant for biologically plausible parameterizations of the diversification process. We find that the purportedly extreme prior sensitivity reported by MEA cannot be replicated with standard usage of BAMM v2.5, or with any other version when conventional Bayesian model selection is performed. Finally, we demonstrate that BAMM performs very well at estimating diversification rate variation across the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\sim}$\end{document}20% of simulated trees in MEA’s data set for which it is theoretically possible to infer rate shifts with confidence. Due to ascertainment bias, the remaining 80% of their purportedly variable-rate phylogenies are statistically indistinguishable from those produced by a constant-rate birth–death process and were thus poorly suited for the summary statistics used in their performance assessment. We demonstrate that inferences about diversification rates have been accurate and consistent across all major previous releases of the BAMM software. We recognize an acute need to address the theoretical foundations of rate-shift models for phylogenetic trees, and we expect BAMM and other modeling frameworks to improve in response to mathematical and computational innovations. However, we remain optimistic that that the imperfect tools currently available to comparative biologists have provided and will continue to provide important insights into the diversification of life on Earth.

LASER: A Maximum Likelihood Toolkit for Detecting Temporal Shifts in Diversification Rates from Molecular Phylogenies

Rates of species origination and extinction can vary over time during evolutionary radiations, and it is possible to reconstruct the history of diversification using molecular phylogenies of extant taxa only. Maximum likelihood methods provide a useful framework for inferring temporal variation in diversification rates. LASER is a package for the R programming environment that implements maximum likelihood methods based on the birth-death process to test whether diversification rates have changed over time. LASER contrasts the likelihood of phylogenetic data under models where diversification rates have changed over time to alternative models where rates have remained constant over time. Major strengths of the package include the ability to detect temporal increases in diversification rates and the inference of diversification parameters under multiple rate-variable models of diversification. The program and associated documentation are freely available from the R package archive at http://cran.r-project.org.

Positive association between population genetic differentiation and speciation rates in New World birds

An implicit assumption of speciation biology is that population differentiation is an important stage of evolutionary diversification, but its significance as a rate-limiting control on phylogenetic speciation dynamics remains largely untested. If population differentiation within a species is related to its speciation rate over evolutionary time, the causes of differentiation could also be driving dynamics of organismal diversity across time and space. Alternatively, geographic variants might be short-lived entities with rates of formation that are unlinked to speciation rates, in which case the causes of differentiation would have only ephemeral impacts. By pairing population genetics datasets from 173 New World bird species (>17,000 individuals) with phylogenetic estimates of speciation rate, we show that the population differentiation rates within species are positively correlated with their speciation rates over long timescales. Although population differentiation rate explains relatively little of the variation in speciation rate among lineages, the positive relationship between differentiation rate and speciation rate is robust to species-delimitation schemes and to alternative measures of both rates. Population differentiation occurs at least three times faster than speciation, which suggests that most populations are ephemeral. Speciation and population differentiation rates are more tightly linked in tropical species than in temperate species, consistent with a history of more stable diversification dynamics through time in the Tropics. Overall, our results suggest that the processes responsible for population differentiation are tied to those that underlie broad-scale patterns of diversity.

Squamate Conserved Loci (SqCL): A unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high-throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches-anchored hybrid enrichment (AHE) and ultraconserved elements (UCE)-and both approaches have proven useful for addressing questions in phylogenomics, phylogeography and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci set, a single integrated probe set that can generate high-quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness and suggest areas for future improvement.

Stable isotope ecology of a hyper-diverse community of scincid lizards from arid Australia

We assessed the utility of stable isotope analysis as a tool for understanding community ecological structure in a species-rich clade of scincid lizards from one of the world's most diverse lizard communities. Using a phylogenetic comparative framework, we tested whether δ¹⁵N and δ¹³C isotopic composition from individual lizards was correlated with species-specific estimates of diet and habitat use. We find that species are highly divergent in isotopic composition with significant correlations to habitat use, but this relationship shows no phylogenetic signal. Isotopic composition corresponds to empirical observations of diet for some species but much variation remains unexplained. We demonstrate the importance of using a multianalytical approach to questions of long-term dietary preference, and suggest that the use of stable isotopes in combination with stomach content analysis and empirical data on habitat use can potentially reveal patterns in ecological traits at finer scales with important implications for community structuring.

Coral snakes predict the evolution of mimicry across New World snakes

Batesian mimicry, in which harmless species (mimics) deter predators by deceitfully imitating the warning signals of noxious species (models), generates striking cases of phenotypic convergence that are classic examples of evolution by natural selection. However, mimicry of venomous coral snakes has remained controversial because of unresolved conflict between the predictions of mimicry theory and empirical patterns in the distribution and abundance of snakes. Here we integrate distributional, phenotypic and phylogenetic data across all New World snake species to demonstrate that shifts to mimetic coloration in nonvenomous snakes are highly correlated with coral snakes in both space and time, providing overwhelming support for Batesian mimicry. We also find that bidirectional transitions between mimetic and cryptic coloration are unexpectedly frequent over both long- and short-time scales, challenging traditional views of mimicry as a stable evolutionary 'end point' and suggesting that insect and snake mimicry may have different evolutionary dynamics.

Unlinked Mendelian inheritance of red and black pigmentation in snakes: Implications for Batesian mimicry

Identifying the genetic basis of mimetic signals is critical to understanding both the origin and dynamics of mimicry over time. For species not amenable to large laboratory breeding studies, widespread color polymorphism across natural populations offers a powerful way to assess the relative likelihood of different genetic systems given observed phenotypic frequencies. We classified color phenotype for 2175 ground snakes (Sonora semiannulata) across the continental United States to analyze morph ratios and test among competing hypotheses about the genetic architecture underlying red and black coloration in coral snake mimics. We found strong support for a two-locus model under simple Mendelian inheritance, with red and black pigmentation being controlled by separate loci. We found no evidence of either linkage disequilibrium between loci or sex linkage. In contrast to Batesian mimicry systems such as butterflies in which all color signal components are linked into a single "supergene," our results suggest that the mimetic signal in colubrid snakes can be disrupted through simple recombination and that color evolution is likely to involve discrete gains and losses of each signal component. Both outcomes are likely to contribute to the exponential increase in rates of color evolution seen in snake mimicry systems over insect systems.

Minimal effects of latitude on present-day speciation rates in New World birds

The tropics contain far greater numbers of species than temperate regions, suggesting that rates of species formation might differ systematically between tropical and non-tropical areas. We tested this hypothesis by reconstructing the history of speciation in New World (NW) land birds using BAMM, a Bayesian framework for modelling complex evolutionary dynamics on phylogenetic trees. We estimated marginal distributions of present-day speciation rates for each of 2571 species of birds. The present-day rate of speciation varies approximately 30-fold across NW birds, but there is no difference in the rate distributions for tropical and temperate taxa. Using macroevolutionary cohort analysis, we demonstrate that clades with high tropical membership do not produce species more rapidly than temperate clades. For nearly any value of present-day speciation rate, there are far more species in the tropics than the temperate zone. Any effects of latitude on speciation rate are marginal in comparison to the dramatic variation in rates among clades.