Viviparity imparts a macroevolutionary signature of ecological opportunity in the body size of female Liolaemus lizards.

Viviparity evolved ~115 times across squamate reptiles, facilitating the colonization of cold habitats, where oviparous species are scarce or absent. Whether the ecological opportunity furnished by such colonization reconfigures phenotypic diversity and accelerates evolution is unclear. We investigated the association between viviparity and patterns and rates of body size evolution in female Liolaemus lizards, the most species-rich tetrapod genus from temperate regions. Here, we discover that viviparous species evolve ~20% larger optimal body sizes than their oviparous relatives, but exhibit similar rates of body size evolution. Through a causal modeling approach, we find that viviparity indirectly influences body size evolution through shifts in thermal environment. Accordingly, the colonization of cold habitats favors larger body sizes in viviparous species, reconfiguring body size diversity in Liolaemus. The catalyzing influence of viviparity on phenotypic evolution arises because it unlocks access to otherwise inaccessible sources of ecological opportunity, an outcome potentially repeated across the tree of life.

There and back again: when and how the world's richest snake family (Dipsadidae) dispersed and speciated across the Neotropical region

Aim The widespread megadiverse Neotropical snake family Dipsadidae occurs in a large range of diverse habitats. Therefore, it represents an excellent model to study the diversification of Neotropical biota. Herein, by generating a time-calibrated species-level phylogeny, we investigate the origin and historical biogeography of Dipsadidae and test if its two main Neotropical subfamilies, Xenodontinae and Dipsadinae, have different geographical origins. Location Neotropical region. Taxon Dipsadidae (Serpentes). Methods We generated a new Bayesian time-calibrated phylogeny based on published sequences from six genes for 344 species, including 287 species of Dipsadidae. We subsequently estimated ancestral areas of distribution by comparing models in BioGeoBEARS: DEC (subset sympatry, narrow vicariance), DIVALIKE (narrow and wide vicariance), BAYAREALIKE (no vicariance and widespread sympatry), also testing jump dispersal. We also estimated shifts in the diversification of this group using BAMM, exploring possible relationships with its historical biogeography. Results The best models show that Dipsadidae likely originated approximately 50 million years ago (mya) in Asia. Dispersal was a fundamental process in its historical biogeography. The DEC model with jump dispersal indicated that this family underwent a range extension from Asia and posterior vicariance of North and Central America ancestors. Both Xenodontinae and Dipsadinae originated in Central America and dispersed to South America during Middle Eocene, but did so to different regions (cis and trans-Andean South America, respectively). Xenodontinae entered cis-Andean South America around 39 mya and jump dispersed to the West Indies around 33 mya, while Dipsadinae entered trans-Andean South America multiple times 20–38 mya. The diversification rate decreased through time, except for a clade within Dipsadinae composed of the Dipsadini tribe and the Atractus and Geophis genera. Main Conclusions Our results show that Dipsadidae has an Asian origin and that the two main Neotropical subfamilies originated in Central America, later dispersing to South America in different time periods. This difference is also reflected in the higher diversification rate for the ‘goo-eaters’ in the Dipsadinae subfamily. The current biogeographical patterns of the family Dipsadidae, the most species-rich snake family in the world, have likely been shaped by complex evolutionary and geological processes such as Eocene land bridges, Andean uplift and the formation of the Panama isthmus.

VenomMaps: Updated species distribution maps and models for New World pitvipers (Viperidae: Crotalinae).

Beyond providing critical information to biologists, species distributions are useful for naturalists, curious citizens, and applied disciplines including conservation planning and medical intervention. Venomous snakes are one group that highlight the importance of having accurate information given their cosmopolitan distribution and medical significance. Envenomation by snakebite is considered a neglected tropical disease by the World Health Organization and venomous snake distributions are used to assess vulnerability to snakebite based on species occurrence and antivenom/healthcare accessibility. However, recent studies highlighted the need for updated fine-scale distributions of venomous snakes. Pitvipers (Viperidae: Crotalinae) are responsible for >98% of snakebites in the New World. Therefore, to begin to address the need for updated fine-scale distributions, we created VenomMaps, a database and web application containing updated distribution maps and species distribution models for all species of New World pitvipers. With these distributions, biologists can better understand the biogeography and conservation status of this group, researchers can better assess vulnerability to snakebite, and medical professionals can easily discern species found in their area.

The role of vicariance and dispersal on the temporal range dynamics of forest vipers in the Neotropical region.

The emergence of the diagonal of open/dry vegetations, including Chaco, Cerrado and Caatinga, is suggested to have acted as a dispersal barrier for terrestrial organisms by fragmenting a single large forest that existed in South America into the present Atlantic and Amazon forests. Here we tested the hypothesis that the expansion of the South American diagonal of open/dry landscapes acted as a vicariant process for forest lanceheads of the genus Bothrops, by analyzing the temporal range dynamics of those snakes. We estimated ancestral geographic ranges of the focal lancehead clade and its sister clade using a Bayesian dated phylogeny and the BioGeoBEARS package. We compared nine Maximum Likelihood models to infer ancestral range probabilities and their related biogeographic processes. The best fitting models (DECTS and DIVALIKETS) recovered the ancestor of our focal clade in the Amazon biogeographic region of northwestern South America. Vicariant processes in two different subclades resulted in disjunct geographic distributions in the Amazon and the Atlantic Forest. Dispersal processes must have occurred mostly within the Amazon and the Atlantic Forest and not between them. Our results suggest the fragmentation of a single ancient large forest into the Atlantic and Amazon forests acting as a driver of vicariant processes for the snake lineage studied, highlighting the importance of the diagonal of open/dry landscapes in shaping distribution patterns of terrestrial biota in South America.

Do key innovations unlock diversification? A case-study on the morphological and ecological impact of pharyngognathy in acanthomorph fishes.

Key innovations may allow lineages access to new resources and facilitate the invasion of new adaptive zones, potentially influencing diversification patterns. Many studies have focused on the impact of key innovations on speciation rates, but far less is known about how they influence phenotypic rates and patterns of ecomorphological diversification. We use the repeated evolution of pharyngognathy within acanthomorph fishes, a commonly cited key innovation, as a case study to explore the predictions of key innovation theory. Specifically, we investigate whether transitions to pharyngognathy led to shifts in the rate of phenotypic evolution, as well as shifts and/or expansion in the occupation of morphological and dietary space, using a dataset of 8 morphological traits measured across 3,853 species of Acanthomorpha. Analyzing the 6 evolutionarily independent pharyngognathous clades together, we found no evidence to support pharyngognathy as a key innovation; however, comparisons between individual pharyngognathous lineages and their sister clades did reveal some consistent patterns. In morphospace, most pharyngognathous clades cluster in areas that correspond to deeper-bodied morphologies relative to their sister clades, whereas occupying greater areas in dietary space that reflects a more diversified diet. Additionally, both Cichlidae and Labridae exhibited higher univariate rates of phenotypic evolution compared with their closest relatives. However, few of these results were exceptional relative to our null models. Our results suggest that transitions to pharyngognathy may only be advantageous when combined with additional ecological or intrinsic factors, illustrating the importance of accounting for lineage-specific effects when testing key innovation hypotheses. Moreover, the challenges we experienced formulating informative comparisons, despite the ideal evolutionary scenario of multiple independent evolutionary origins of pharyngognathous clades, illustrates the complexities involved in quantifying the impact of key innovations. Given the issues of lineage specific effects and rate heterogeneity at macroevolutionary scales we observed, we suggest a reassessment of the expected impacts of key innovations may be warranted.

Exploring the drivers of population structure across desert snakes can help to link micro and macroevolution.

To understand the underlying mechanisms generating population genetic divergence and structure is a critical step towards understanding how biodiversity evolves at both micro- and macroevolutionary scales. At the population-level, geographic isolation as well as adaptation to local environmental conditions can generate different patterns of spatial genetic variation among populations. Specific organismal traits as well as the characteristics of the environment might influence the process under which populations become spatially structured. In a From the Cover article in this issue of Molecular Ecology, Myers et al. (2019) present an integrative approach to investigate if the Cochise filter barrier (CFB), lying between the Sonoran and Chihuahuan Deserts, and the surrounding river networks were relevant in driving the population structure of 13 snake species. While local environmental conditions seem to predominantly contribute to lineage divergence, traditionally studied vicariant barriers seem to have played a minor role in shaping population structure across the studied species. This study brings insights into how population-level processes could contribute to the formation of incipient species, which ultimately might affect the speciation rates measured at macroevolutionary scales. Hence, Myers et al. (2019) not only represents an integrative study aiming to understand the drivers of population genetic divergence, but also a potentially important contribution to our ongoing challenge in linking micro- and macroevolution.

How Well Can We Estimate Diversity Dynamics for Clades in Diversity Decline?

The fossil record shows that the vast majority of all species that ever existed are extinct and that most lineages go through an expansion and decline in diversity. However, macroevolutionary analyses based upon molecular phylogenies have difficulty inferring extinction dynamics, raising questions about whether the neontological record can contribute to an understanding of the decline phenomenon. Two recently developed diversification methods for molecular phylogenies (RPANDA and BAMM) incorporate models that theoretically have the capacity to capture decline dynamics by allowing extinction to be higher than speciation. However, the performance of these frameworks over a wide range of decline scenarios has not been studied. Here, we investigate the behavior of these methods under decline scenarios caused by decreasing speciation and increasing extinction through time on simulated trees at fixed intervals over diversity trajectories with expansion and decline phases. We also compared method performance over a comprehensive data set of 214 empirical trees. Our results show that both methods perform equally well when varying speciation rates control decline. When decline was only caused by an increase in extinction rates both methods wrongly assign the variation in net diversification to a drop in speciation, even though the positive gamma values of those trees would suggest otherwise. We also found a tendency for RPANDA to favor increasing extinction and BAMM to favor decreasing speciation as the most common cause of decline in empirical trees. Overall our results shed light on the limitations of both methods, encouraging researchers to carefully interpret the results from diversification studies.

Diversification in vipers: Phylogenetic relationships, time of divergence and shifts in speciation rates.

Snakes of the cosmopolitan family Viperidae comprise around 329 venomous species showing a striking heterogeneity in species richness among lineages. While the subfamily Azemiopinae comprises only two species, 70% of all viper species are arranged in the subfamily Crotalinae or the "pit vipers". The radiation of the pit vipers was marked by the evolution of the heat-sensing pits, which has been suggested to be a key innovation for the successful diversification of the group. Additionally, only crotalines were able to successfully colonize the New World. Here, we present the most complete molecular phylogeny for the family to date that comprises sequences from nuclear and mitochondrial genes representing 79% of all living vipers. We also investigated the time of divergence between lineages, using six fossils to calibrate the tree, and explored the hypothesis that crotalines have undergone an explosive radiation. Our phylogenetic analyses retrieved high support values for the monophyly of the family Viperidae, subfamilies Viperinae and Crotalinae, and 22 out of 27 genera, as well as well-supported intergeneric relationships throughout the family. We were able to recover a strongly supported sister clade to the New World pit vipers that comprises Gloydius, Ovophis, Protobothrops and Trimeresurus gracilis. Our results agree in many aspects with other studies focusing on the phylogenetics of vipers, but we recover new relationships as well. Despite the addition of new sequences we were not able to resolve some of the poor supported relationships previously suggested. Time of divergence estimates suggested that vipers started to radiate around the late Paleocene to middle Eocene with subfamilies most likely dating back to the Eocene. The invasion of the New World might have taken place sometime close to the Oligocene/Miocene boundary. Diversification analyses suggested a shift in speciation rates during the radiation of a sub-clade of pit vipers where speciation rates rapidly increased but slowed down toward the present. Thus, the evolution of the loreal pits alone does not seem to explain their explosive speciation rates. We suggest that climatic and geological changes in Asia and the invasion of the New World may have also contributed to the speciation shift found in vipers.

Seasonal Variation in Population Abundance and Chytrid Infection in Stream-Dwelling Frogs of the Brazilian Atlantic Forest.

Enigmatic amphibian declines were first reported in southern and southeastern Brazil in the late 1980s and included several species of stream-dwelling anurans (families Hylodidae and Cycloramphidae). At that time, we were unaware of the amphibian-killing fungus Batrachochytrium dendrobatidis (Bd); therefore, pollution, habitat loss, fragmentation and unusual climatic events were hypothesized as primary causes of these declines. We now know that multiple lineages of Bd have infected amphibians of the Brazilian Atlantic forest for over a century, yet declines have not been associated specifically with Bd outbreaks. Because stream-dwelling anurans occupy an environmental hotspot ideal for disease transmission, we investigated temporal variation in population and infection dynamics of three stream-adapted species (Hylodes asper, H. phyllodes, and Cycloramphus boraceiensis) on the northern coast of São Paulo state, Brazil. We surveyed standardized transects along streams for four years, and show that fluctuations in the number of frogs correlate with specific climatic variables that also increase the likelihood of Bd infections. In addition, we found that Bd infection probability in C. boraceiensis, a nocturnal species, was significantly higher than in Hylodes spp., which are diurnal, suggesting that the nocturnal activity may either facilitate Bd zoospore transmission or increase susceptibility of hosts. Our findings indicate that, despite long-term persistence of Bd in Brazil, some hosts persist with seasonally variable infections, and thus future persistence in the face of climate change will depend on the relative effect of those changes on frog recruitment and pathogen proliferation.