Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs

Widespread sympatry in a species-rich clade of marine fishes (Carangoidei)

A universal paradigm describing patterns of speciation across the tree of life has been debated for decades. In marine organisms, inferring patterns of speciation using contemporary and historical patterns of biogeography is challenging due to the deficiency of species-level phylogenies and information on species' distributions, as well as conflicting relationships between species' dispersal, range size and co-occurrence. Most research on global patterns of marine fish speciation and biogeography has focused on coral reef or pelagic species. Carangoidei is an ecologically important clade of marine fishes that use coral reef and pelagic environments. We used sequence capture of 1314 ultraconserved elements (UCEs) from 154 taxa to generate a time-calibrated phylogeny of Carangoidei and its parent clade, Carangiformes. Age-range correlation analyses of the geographical distributions and divergence times of sister species pairs reveal widespread sympatry, with 73% of sister species pairs exhibiting sympatric geographical distributions, regardless of node age. Most species pairs coexist across large portions of their ranges. We also observe greater disparity in body length and maximum depth between sympatric relative to allopatric sister species. These and other ecological or behavioural attributes probably facilitate sympatry among the most closely related carangoids.

Integrative revision of the Blommersia wittei complex, with description of a new species of frog from western and north-western Madagascar

Frogs of the Blommersia wittei complex are widespread in western and northern Madagascar, and are one of two clades of the family Mantellidae that have colonized the Comoran island of Mayotte. Based on a comprehensive set of DNA sequences of the mitochondrial 16S rRNA gene and the nuclear-encoded RAG1 and SACS genes, integrated with morphological and bioacoustic data, we here analyze the genetic differentiation of populations of this complex across Madagascar. We confirm that a candidate species named B. sp. Ca5 in previous studies represents a genetically well-defined evolutionary lineage distributed over much of western Madagascar, which we describe herein as Blommeria bara sp. nov. based on its molecular and bioacoustic differentiation. Blommersia wittei occurs across northern Madagascar but its type locality Ambanja, at the lower Sambirano river, is very close to the range of another, newly discovered microendemic lineage that was only found at two sites along the upper Sambirano river (here named as candidate species B. sp. Ca12). The B. wittei complex thus provides an example of a clade of closely related Malagasy frogs that contains species widespread over hundreds of kilometers, as well as extreme microendemics. For a full resolution of this species complex, more data need to be collected on the geographical contact among these two lineages, on the morphology and bioacoustics of B. sp. Ca12, and on the north-eastern populations of B. wittei at Sambava, which are weakly differentiated in mitochondrial genes but differ in bioacoustics and possibly in the extent of foot webbing.

Body size predicts the rate of contemporary morphological change in birds

Variation in evolutionary rates among species is a defining characteristic of the tree of life and may be an important predictor of species' capacities to adapt to rapid environmental change. It is broadly assumed that generation length is an important determinant of microevolutionary rates, and body size is often used as a proxy for generation length. However, body size has myriad biological correlates that could affect evolutionary rates independently from generation length. We leverage two large, independently collected datasets on recent morphological change in birds (52 migratory species breeding in North America and 77 South American resident species) to test how body size and generation length are related to the rates of contemporary morphological change. Both datasets show that birds have declined in body size and increased in wing length over the past 40 y. We found, in both systems, a consistent pattern wherein smaller species declined proportionally faster in body size and increased proportionally faster in wing length. By contrast, generation length explained less variation in evolutionary rates than did body size. Although the mechanisms warrant further investigation, our study demonstrates that body size is an important predictor of contemporary variation in morphological rates of change. Given the correlations between body size and a breadth of morphological, physiological, and ecological traits predicted to mediate phenotypic responses to environmental change, the relationship between body size and rates of phenotypic change should be considered when testing hypotheses about variation in adaptive responses to climate change.

Analysis of intrinsic evolutionary factors leading to microendemic distributions in New Caledonian leaf beetles

Microendemicity, or the condition of some species having local ranges, is a relatively common pattern in nature. However, the factors that lead to this pattern are still largely unknown. Most studies addressing this issue tend to focus on extrinsic factors associated with microendemic distributions, such as environmental conditions, hypothesising a posteriori about underlying potential speciation mechanisms, linked or not to these conditions. Here, we use a multi-faceted approach mostly focusing on intrinsic factors instead, namely diversification dynamics and speciation modes in two endemic sibling genera of leaf beetles with microendemic distributions, Taophila and Tricholapita, in a microendemicity hotspot, New Caledonia. Results suggest that the diversification rate in this lineage slowed down through most of the Neogene and consistently with a protracted speciation model possibly combined with several ecological and environmental factors potentially adding rate-slowing effects through time. In turn, species accumulated following successive allopatric speciation cycles, possibly powered by marked geological and climatic changes in the region in the last 25 million years, with daughter species ranges uncorrelated with the time of speciation. In this case, microendemicity seems to reflect a mature state for the system, rather than a temporary condition for recent species, as suggested for many microendemic organisms.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Slow and steady wins the race: Diversification rate is independent from body size and lifestyle in Malagasy skinks (Squamata: Scincidae: Scincinae)

Most of the unique and diverse vertebrate fauna that inhabits Madagascar derives from in situ diversification from colonisers that reached this continental island through overseas dispersal. The endemic Malagasy Scincinae lizards are amongst the most species-rich squamate groups on the island. They colonised all bioclimatic zones and display many ecomorphological adaptations to a fossorial (burrowing) lifestyle. Here we propose a new phylogenetic hypothesis for their diversification based on the largest taxon sampling so far compiled for this group. We estimated divergence times and investigated several aspects of their diversification (diversification rate, body size and fossorial lifestyle evolution, and biogeography). We found that diversification rate was constant throughout most of the evolutionary history of the group, but decreased over the last 6-4 million years and independently from body size and fossorial lifestyle evolution. Fossoriality has evolved from fully quadrupedal ancestors at least five times independently, which demonstrates that even complex morphological syndromes - in this case involving traits such as limb regression, body elongation, modification of cephalic scalation, depigmentation, and eyes and ear-opening regression - can evolve repeatedly and independently given enough time and eco-evolutionary advantages. Initial diversification of the group likely occurred in forests, and the divergence of sand-swimmer genera around 20 Ma appears linked to a period of aridification. Our results show that the large phenotypic variability of Malagasy Scincinae has not influenced diversification rate and that their rich species diversity results from a constant accumulation of lineages through time. By compiling large geographic and trait-related datasets together with the computation of a new time tree for the group, our study contributes important insights on the diversification of Malagasy vertebrates.

Low microbiome diversity in threatened amphibians from two biodiversity hotspots

Microbial diversity positively influences community resilience of the host microbiome. However, extinction risk factors such as habitat specialization, narrow environmental tolerances, and exposure to anthropogenic disturbance may homogenize host-associated microbial communities critical for stress responses including disease defense. In a dataset containing 43 threatened and 90 non-threatened amphibian species across two biodiversity hotspots (Brazil's Atlantic Forest and Madagascar), we found that threatened host species carried lower skin bacterial diversity, after accounting for key environmental and host factors. The consistency of our findings across continents suggests the broad scale at which low bacteriome diversity may compromise pathogen defenses in species already burdened with the threat of extinction.

The evolution of reproductive modes and life cycles in amphibians

Amphibians have undergone important evolutionary transitions in reproductive modes and life-cycles. We compare large-scale macroevolutionary patterns in these transitions across the three major amphibian clades: frogs, salamanders, and caecilians. We analyse matching reproductive and phylogenetic data for 4025 species. We find that having aquatic larvae is ancestral for all three groups and is retained by many extant species (33-44%). The most frequent transitions in each group are to relatively uncommon states: live-bearing in caecilians, paedomorphosis in salamanders, and semi-terrestriality in frogs. All three groups show transitions to more terrestrial reproductive modes, but only in caecilians have these evolved sequentially from most-to-least aquatic. Diversification rates are largely independent of reproductive modes. However, in salamanders direct development accelerates diversification whereas paedomorphosis decreases it. Overall, we find a widespread retention of ancestral modes, decoupling of trait transition rates from patterns of species richness, and the general independence of reproductive modes and diversification.

Melanesia holds the world's most diverse and intact insular amphibian fauna

Identifying hotspots of biological diversity is a key step in conservation prioritisation. Melanesia-centred on the vast island of New Guinea-is increasingly recognised for its exceptionally species-rich and endemic biota. Here we show that Melanesia has the world's most diverse insular amphibian fauna, with over 7% of recognised global frog species in less than 0.7% of the world's land area, and over 97% of species endemic. We further estimate that nearly 200 additional candidate species have been discovered but remain unnamed, pointing to a total fauna in excess of 700 species. Nearly 60% of the Melanesian frog fauna is in a lineage of direct-developing microhylids characterised by smaller distributions than co-occurring frog families, suggesting lineage-specific high beta diversity is a key driver of Melanesian anuran megadiversity. A comprehensive conservation status assessment further highlights geographic concentrations of recently described range-restricted threatened taxa that warrant urgent conservation actions. Nonetheless, by world standards, the Melanesian frog fauna is relatively intact, with 6% of assessed species listed as threatened and no documented extinctions; and thus it provides an unparalleled opportunity to understand and conserve a megadiverse and relatively intact insular biota.

Coalescent simulations indicate that the São Francisco River is a biogeographic barrier for six vertebrates in a seasonally dry South American forest

The riverine barrier hypothesis has been extensively explored in Neotropical rainforests, while its importance in drier regions such as the Caatinga, a seasonally dry tropical forest in northeastern Brazil, has only recently received more attention. The Caatinga is bisected by the São Francisco River (SFR), which has long been suggested to be an important biogeographic feature in the region. However, recent studies have found mixed support for the role of the SFR as a hard barrier, most of them relying on the presence or absence of genetic breaks congruent with its course. Here, we used published multilocus and next-generation data from six vertebrate species to test the SFR’s strength as a barrier. Using model-based approaches (approximate Bayesian computation and supervised machine learning), we tested demographic models incorporating full, intermediate, and zero migration across the SFR, estimating divergence times and migration rates for each species. We found support for the SFR’s role as a barrier, allowing gene flow for some species. Estimated divergence times varied among species but are limited to the late Pleistocene, coherent with one of several proposed paleocourse changes in the river’s geological history. Contrary to the mixed results of previous studies, our study supports the SFR as an important phylogeographic barrier across different taxonomic groups, driving diversification in the Caatinga.

Understanding the relationship between dispersal and range size

The drivers of variability in species range sizes remain an outstanding enigma in ecology. The theoretical expectation of a positive dispersal-range size relationship has received mixed empirical support, despite dispersal being one of the most prominent hypothesised predictors of range size. Here, we synthesised results from 86 studies examining the dispersal-range size relationship for plants and animals in marine, terrestrial and freshwater realms. Overall, our meta-analysis showed that dispersal positively affects range size, but its effect is dependent on the clade and dispersal proxy studied. Moreover, despite potential differences in habitat connectivity, we did not find an effect of realm on the dispersal-range size relationship. Finally, the strength of the dispersal-range size relationship was dependent on latitude, range size metric and the taxonomic breadth of the study clade. Our synthesis emphasizes the importance of developing a mechanistic understanding of the trait to dispersal to range size relationship, considering the complexity of dispersal departure, transfer and settlement, as well as evolutionary components such as time for range expansion, speciation and past geological-environmental dynamics. We, therefore, call for a more integrative view of the dispersal process and its causal relationship with range size.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

An unexpected new red-bellied Stumpffia (Microhylidae) from forest fragments in central Madagascar highlights remaining cryptic diversity

The Madagascan endemic subfamily Cophylinae in the family Microhylidae, is an example of a taxonomic group for which much is still to be discovered. Indeed, the cophyline frogs present a large portion of Madagascar's cryptic and microendemic amphibian diversity, yet they remain understudied. A new red-bellied species of the microhylid frog genus Stumpffia is described from the central plateau of Madagascar. Visual encounter surveys in Ambohitantely and Anjozorobe in 2019 and 2020 identified this previously unknown Stumpffia species, which closely resembles Stumpffiakibomena known from Andasibe in the east. Stumpffialynnae sp. nov. adds another species to the red-bellied species complex, differing from S.kibomena by genetic differentiation in the mitochondrial 16S rRNA gene (3.6-3.9%) and distinct nuclear RAG1 haplotypes, as well as strongly by its advertisement call. The new species is known from across Ambohitantely Special Reserve and Anjozorobe Angavo protected area, but is known only from one complete specimen and eight individual tissue samples. Based on the rarity of the species, the small number of locations in which it has been found, and its disappearing forest habitat, its IUCN Red List classification is suggested as "Endangered". This species is the first Stumpffia described from Madagascar's central plateau, highlighting the importance of conserving the remnant forest fragments in this area and the ongoing need to survey and protect this threatened habitat type.

Population diversification in the frog Mantidactylus bellyi on an isolated massif in northern Madagascar based on genetic, morphological, bioacoustic and ecological evidence

In the processes that give rise to new species, changes first occur at the population level. But with the continuous nature of the divergence process, change in biological properties delimiting the shift from “individuals of divergent populations” towards “individuals of distinct species”, as well as abiotic factors driving the change, remain largely ambivalent. Here we study diversification processes at the population level in a semi-aquatic frog, Mantidactylus (Brygoomantis) bellyi, across the diverse vegetation types of Montagne d’Ambre National Park (MANP), Madagascar. Genetic diversity was assessed with seven newly developed microsatellite markers as well as mitochondrial DNA sequences and concordance with patterns of ecological, morphological, and bioacoustic divergence evaluated. We found M. bellyi lacking mitochondrial differentiation within MANP, while microsatellite datasets partitioned them into three highly differentiated, geographically separated subpopulations (with indications for up to five subpopulations). The molecular grouping–primarily clustering individuals by geographic proximity–was coincident with differences in mean depth and width of waters, suggesting a possible role of fluvial characteristics in genetic exchange in this stream-breeding species. Genetic clustering not consistent with differences in call properties, except for dominant call frequencies under the two-subpopulations model. Morphological divergence was mostly consistent with the genetic clustering; subpopulations strongly differed by their snout-vent length, with individuals from high-elevation subpopulations smaller than those from populations below 1000 m above sea level. These results exemplify how mountains and environmental conditions might primarily shape genetic and morphological divergence in frog populations, without strongly affecting their calls.

Comparative cytogenetics on eight Malagasy Mantellinae (Anura, Mantellidae) and a synthesis of the karyological data on the subfamily

We performed a molecular and cytogenetic analysis on different Mantellinae species and revised the available chromosomal data on this group to provide an updated assessment of its karyological diversity and evolution. Using a fragment of the mitochondrial 16S rRNA, we performed a molecular taxonomic identification of the samples that were used for cytogenetic analyses. A comparative cytogenetic analysis, with Giemsa's staining, Ag-NOR staining and sequential C-banding + Giemsa + CMA + DAPI was performed on eight species: Gephyromantis sp. Ca19, G.striatus (Vences, Glaw, Andreone, Jesu et Schimmenti, 2002), Mantidactylus (Chonomantis) sp. Ca11, M. (Brygoomantis) alutus (Peracca, 1893), M. (Hylobatrachus) cowanii (Boulenger, 1882), Spinomantispropeaglavei "North" (Methuen et Hewitt, 1913), S.phantasticus (Glaw et Vences, 1997) and S. sp. Ca3. Gephyromantisstriatus, M. (Brygoomantis) alutus and Spinomantispropeaglavei "North" have a karyotype of 2n = 24 chromosomes while the other species show 2n = 26 chromosomes. Among the analysed species we detected differences in the number and position of telocentric elements, location of NOR loci (alternatively on the 6th, 7th or 10th pair) and in the distribution of heterochromatin, which shows species-specific patterns. Merging our data with those previously available, we propose a karyotype of 2n = 26 with all biarmed elements and loci of NORs on the 6th chromosome pair as the ancestral state in the whole family Mantellidae. From this putative ancestral condition, a reduction of chromosome number through similar tandem fusions (from 2n = 26 to 2n = 24) occurred independently in Mantidactylus Boulenger, 1895 (subgenus Brygoomantis Dubois, 1992), Spinomantis Dubois, 1992 and Gephyromantis Methuen, 1920. Similarly, a relocation of NORs, from the putative primitive configuration on the 6th chromosome, occurred independently in Gephyromantis, Blommersia Dubois, 1992, Guibemantis Dubois, 1992, Mantella Boulenger, 1882 and Spinomantis. Chromosome inversions of primitive biarmed elements likely generated a variable number of telocentric elements in Mantellanigricans Guibé, 1978 and a different number of taxa of Gephyromantis (subgenera Duboimantis Glaw et Vences, 2006 and Laurentomantis Dubois, 1980) and Mantidactylus (subgenera Brygoomantis, Chonomantis Glaw et Vences, 1994, Hylobatrachus Laurent, 1943 and Ochthomantis Glaw et Vences, 1994).

A fantastic new species of secretive forest frog discovered from forest fragments near Andasibe, Madagascar

We describe a fantastic new species of forest frog (Mantellidae: Gephyromantis: subgenus: Laurentomantis) from moderately high elevations in the vicinity of Andasibe, Madagascar. This region has been surveyed extensively and has a remarkably high anuran diversity with many undocumented species still being discovered. Surprisingly, by exploring areas around Andasibe that lacked biodiversity surveys, we discovered a spectacular and clearly morphologically distinct species, previously unknown to science, Gephyromantis marokorokosp. nov., documented for the first time in 2015. The new species is well characterised by a very rugose and granular dorsum, dark brown skin with bright red mottling, sparse light orange to white spots on the ventre, vibrant red eyes and femoral glands present only in males that consist of eight medium-sized granules. Bioacoustically, the new species has a quiet advertisement call that differs from related species by having a moderate call duration, 2–4 strongly pulsed notes and a slow note repetition rate. Furthermore, it has substantial differentiation in mitochondrial DNA, with pairwise distances of 7–9% to all other related species in sequences of the mitochondrial 16S rRNA marker. Additional evidence is given through a combined four mitochondrial markers and four nuclear exons concatenated species tree, strongly supporting G. striatus as the sister species of the new species in both analyses. The discovery of this new species highlights the need for continued inventory work in high elevation rainforests of Madagascar, even in relatively well-studied regions.

The Pristimantis trachyblepharis species group, a clade of miniaturized frogs: description of four new species and insights into the evolution of body size in the genus

Species richness in the genus Pristimantis is underestimated due to the existence of morphologically cryptic species. This is worsened by the low sampling effort and the lack of studies using genetic markers. Here, we use molecular and morphological data to determine the phylogenetic relationships of a clade of Pristimantis distributed throughout montane tropical forests in the eastern Andes, from central Ecuador to northern Perú. We name this clade the Pristimantis trachyblepharis species group. Our results show that it comprises nine species, of which four are formally described and five are new. Four of these undescribed species are formally described here. The group is composed of miniaturized species, such as Pristimantis nanus sp. nov., currently the smallest known species of the genus and the smallest vertebrate in Ecuador. As a first approach to understanding the evolutionary origin and implications of body-size reduction in Pristimantis, we here study the phylogenetic signal and evolutionary trends of body size within the genus. We also provide the first record of P. aquilonaris in Ecuador and we show, for the first time, the phylogenetic position of P. albujai, P. aquilonaris, P. minimus and P. trachyblepharis, which are also members of the P. trachyblepharis species group.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA's diverse biotas.

Machine learning approaches identify male body size as the most accurate predictor of species richness

BACKGROUND

A major challenge in biodiversity science is to understand the factors contributing to the variability of species richness -the number of different species in a community or region - among comparable taxonomic lineages. Multiple biotic and abiotic factors have been hypothesized to have an effect on species richness and have been used as its predictors, but identifying accurate predictors is not straightforward. Spiders are a highly diverse group, with some 48,000 species in 120 families; yet nearly 75% of all species are found within just the ten most speciose families. Here we use a Random Forest machine learning algorithm to test the predictive power of different variables hypothesized to affect species richness of spider genera.

RESULTS

We test the predictive power of 22 variables from spiders' morphological, genetic, geographic, ecological and behavioral landscapes on species richness of 45 genera selected to represent the phylogenetic and biological breath of Araneae. Among the variables, Random Forest analyses find body size (specifically, minimum male body size) to best predict species richness. Multiple Correspondence analysis confirms this outcome through a negative relationship between male body size and species richness. Multiple Correspondence analyses furthermore establish that geographic distribution of congeneric species is positively associated with genus diversity, and that genera from phylogenetically older lineages are species poorer. Of the spider-specific traits, neither the presence of ballooning behavior, nor sexual size dimorphism, can predict species richness.

CONCLUSIONS

We show that machine learning analyses can be used in deciphering the factors associated with diversity patterns. Since no spider-specific biology could predict species richness, but the biologically universal body size did, we believe these conclusions are worthy of broader biological testing. Future work on other groups of organisms will establish whether the detected associations of species richness with small body size and wide geographic ranges hold more broadly.

Two-hundred million years of anuran body-size evolution in relation to geography, ecology and life history

Surprisingly, little is known about body-size evolution within the most diverse amphibian order, anurans (frogs and toads), despite known effects of body size on the physiological, ecological and life-history traits of animals more generally. Here, we examined anuran body-size evolution among 2,434 species with over 200 million years of shared evolutionary history. We found clade-specific evolutionary shifts to new body-size optima along with numerous independent transitions to gigantic and miniature body sizes, despite the upper limits of anuran body size remaining quite consistent throughout the fossil record. We found a weak, positive correlation between a species' body size and maximum latitude and elevation, including a dearth of small species at higher elevations and broader latitudinal and elevational ranges in larger anurans. Although we found modest differences in mean anuran body size among microhabitats, there was extensive overlap in the range of body sizes across microhabitats. Finally, we found that larger anurans are more likely to consume vertebrate prey than smaller anurans are and that species with a free-swimming larval phase during development are larger on average than those in which development into a froglet occurs within the egg. Overall, anuran body size does not conform to geographic and ecological patterns observed in other tetrapods but is perhaps more notable for variation in body size within geographic regions, ecologies and life histories. Here, we document this variation and propose target clades for detailed studies aimed at disentangling how and why variation in body size was generated and is maintained in anurans.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

The geographical diversification in varanid lizards: the role of mainland versus island in driving species evolution

Monitor lizards (Varanidae) inhabit both the mainland and islands of all geological types and have diversified into an exceptionally wide range of body sizes, thus providing an ideal model for examining the role of mainland versus island in driving species evolution. Here we use phylogenetic comparative methods to examine whether a link exists between body size-driven diversification and body size-frequency distributions in varanid lizards and to test the hypothesis that island lizards differ from mainland species in evolutionary processes, body size, and life-history traits (offspring number and size). We predict that: 1) since body size drives rapid diversification in groups, a link exists between body size-driven diversification and body size-frequency distributions; 2) because of various environments on island, island species will have higher speciation, extinction, and dispersal rates, compared with mainland species; 3) as a response to stronger intraspecific competition, island species will maximize individual ability associated with body size to outcompete closely-related species, and island species will produce smaller clutches of larger eggs to increase offspring quality. Our results confirm that the joint effect of differential macroevolutionary rates shapes the species richness pattern of varanid lizards. There is a link between body size-driven diversification and body size-frequency distributions, and the speciation rate is maximized at medium body sizes. Island species will have higher speciation, equal extinction, and higher dispersal rates compared with mainland species. Smaller clutch size and larger hatchling in the island than in mainland species indicate that offspring quality is more valuable than offspring quantity for island varanids.

Historical biogeography and multi-trait evolution in miniature toadlets of the genus Brachycephalus (Anura: Brachycephalidae)

Evolutionary changes towards a miniaturized body plan may directly affect other important phenotypic traits related to the physiology, behaviour and ecology of organisms. The frog genus Brachycephalus is an outstanding example of a radiation of miniaturized species endemic to the Brazilian Atlantic Forest. We inferred ancestral states and historical changes in body size, body colour and hyperossification to test hypotheses about diversification and selective environmental mechanisms leading to the evolution of these specialized traits. The ancestral distribution was associated with high-elevation regions in the northern Serra do Mar mountain range, and diversification in the genus was coincident with important geological and climatic events during the history of the Atlantic Forest. The dynamic historical changes provided an opportunity for multiple lowland lineages and for speciation via dispersal and vicariance in multiple invasions of the highlands. The ancestral Brachycephalus was reconstructed as miniaturized and dull coloured, without hyperossification in the skin, skull or postcranial skeleton. A parallel evolution of phenotypic traits has occurred in northern and southern Atlantic Forest lineages, beginning in the Miocene. Shifts in body size are not related to elevation range or latitude. However, we found a significant correlation between the evolution of hyperossification and aposematism with increasing body size.

Mode of miniaturisation influences body shape evolution in New World anchovies (Engraulidae)

We explored the macroevolutionary dynamics of miniaturisation in New World anchovies by integrating a time-calibrated phylogeny, geometric morphometrics and phylogenetic comparative methods. We found that the paedomorphic species Amazonsprattus scintilla occupies a novel region of shape space, while the dwarf species Anchoviella manamensis has an overall shape consistent with other anchovies. We found that miniaturisation did not increase overall clade disparity in size or shape beyond the expectations of Brownian motion, nor were there differences in rates of size or shape evolution among clades. Overall, our study shows that while the mode of miniaturisation influences shape evolution, the phenotypic novelty produced by the evolution of miniaturisation did not seem to alter macroevolutionary dynamics.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Species complexes and the importance of Data Deficient classification in Red List assessments: The case of Hylobatrachus frogs

Taxonomy is the cornerstone of extinction risk assessments. Currently, the IUCN Red List treats species complexes either under a single overarching species name—resulting in an unhelpfully broad circumscription and underestimated threat assessment that does not apply to any one species lineage—or omits them altogether—resulting in the omission of species that should be assessed. We argue that taxonomic uncertainty alone, as in species complexes, should be grounds for assessment as Data Deficient (DD). Yet, use of the DD category is currently discouraged, resulting in assessments based on poor data quality and dismissal of the importance of taxonomic confidence in conservation. This policy may be leading to volatile and unwarranted assessments of hundreds of species across the world, and needs to be revised. To illustrate this point, we here present a partial taxonomic revision of torrent frogs from eastern Madagascar in the Mantidactylus subgenus Hylobatrachus. Two named species, Mantidactylus (Hylobatrachus) lugubris and M. (H.) cowanii, and several undescribed candidate species are recognised, but the application of the available names has been somewhat ambiguous. In a recent re-assessment of its conservation status, M. (H.) lugubris was assessed including all complex members except M. (H.) cowanii within its distribution, giving it a status of Least Concern and distribution over most of eastern Madagascar. After describing two of the unnamed lineages as Mantidactylus (Hylobatrachus) atsimo sp. nov. (from southeastern Madagascar) and Mantidactylus (Hylobatrachus) petakorona sp. nov. (from the Marojejy Massif in northeastern Madagascar), we show that Mantidactylus (Hylobatrachus) lugubris is restricted to the central east of Madagascar, highlighting the inaccuracy of its current Red List assessment. We propose to re-assess its status under a more restrictive definition that omits well-defined candidate species, thus representing the actual species to which its assessment refers, to the best of current knowledge. We recommend that for species complexes in general, (1) nominal lineages that can be confidently restricted should be assessed under the strict definition, (2) non-nominal species-level lineages and ambiguous names should be prioritised for taxonomic research, and (3) ambiguous names should be assessed as DD to highlight the deficiency in data on their taxonomic status, which is an impediment to their conservation. This would reduce ambiguity and underestimation of threats involved in assessing species complexes, and place the appropriate emphasis on the importance of taxonomy in anchoring conservation.

Mitochondrial Introgression, Color Pattern Variation, and Severe Demographic Bottlenecks in Three Species of Malagasy Poison Frogs, Genus Mantella

Madagascar is a biodiversity hotspot particularly rich in amphibian diversity and only a few charismatic Malagasy amphibians have been investigated for their population-level differentiation. The Mantellamadagascariensis group is composed of two rainforest and three swamp forest species of poison frogs. We first confirm the monophyly of this clade using DNA sequences of three nuclear and four mitochondrial genes, and subsequently investigate the population genetic differentiation and demography of the swamp forest species using one mitochondrial, two nuclear and a set of nine microsatellite markers. Our results confirm the occurrence of two main mitochondrial lineages, one dominated by Mantellaaurantiaca (a grouping supported also by our microsatellite-based tree) and the other by Mantellacrocea + Mantellamilotympanum. These two main lineages probably reflect an older divergence in swamp Mantella. Widespread mitochondrial introgression suggests a fairly common occurrence of inter-lineage gene flow. However, nuclear admixture seems to play only a limited role in this group, and the analyses of the RAG-1 marker points to a predominant incomplete lineage sorting scenario between all five species of the group, which probably diverged relatively recently. Our demographic analyses show a common, severe and recent demographic contraction, inferred to be in temporal coincidence with the massive deforestation events that took place in the past 1000 years. Current data do not allow to conclusively delimit independent evolutionary units in these frogs, and we therefore refrain to suggest any taxonomic changes.

Global Rate Variation in Bony Vertebrates

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

Factors affecting interspecific differences in genetic divergence among populations of Anolis lizards in Cuba

BACKGROUND

Geographical patterns and degrees of genetic divergence among populations differ between species, reflecting relative potentials for speciation or cladogenesis and differing capacities for environmental adaptation. Identification of factors that contribute to genetic divergence among populations is important to the understanding of why some species exhibit greater interpopulation genetic divergence. In this study, we calculated the mean pairwise genetic distances among populations as species' average genetic divergence by a phylogeny using nuclear and mitochondrial genes of 303 individuals from 33 Cuban Anolis species and estimated species ages by another phylogeny using nuclear and mitochondrial genes of 51 Cuban and 47 non-Cuban Anolis species. We identified factors that influence species' differences in genetic divergence among 26 species of Anolis lizards from Cuba. Species ages, environmental heterogeneity within species ranges, and ecomorph types were considered as factors affecting average genetic divergences among populations.

RESULTS

The phylogenies presented in this study provide the most comprehensive sampling of Cuban Anolis species to date. The phylogeny showed more conservative evolution of Anolis ecomorphs within Cuba and identified twig anoles as a monophyletic group. Subsequent Phylogenetic Generalized Least Squares (PGLS) analyses showed that species age was positively correlated with species' average genetic divergence among populations.

CONCLUSION

Although previous studies have focused on factors affecting genetic divergence within species, the present study showed for the first time that species differences in genetic divergence could be largely affected by species age.

Most species are not limited by an Amazonian river postulated to be a border between endemism areas

At broad scales in the Amazon, it is often hypothesized that species distributions are limited by geographical barriers, such as large rivers (river-barrier hypothesis). This hypothesis has been used to explain the spatial-distribution limits of species and to indicate endemism areas for several phylogenetic lineages. We tested the ability of the river-barrier hypothesis to explain patterns of species diversity and spatial-distribution limits for 1952 easily-detected species in 14 taxonomic groups that occur around the Madeira River, and our results indicate that the hypothesis that the Madeira River is the border between endemism areas and explains much of the diversity found in the region is inappropriate for >99% of species. This indicates that alternative hypotheses should be proposed to explain the limits of distributions of species around the Madeira River, as well as a revision of the criteria that are used to determine species-endemism areas.

Syntopic frogs reveal different patterns of interaction with the landscape: A comparative landscape genetic study of Pelophylax nigromaculatus and Fejervarya limnocharis from central China

Amphibians are often considered excellent environmental indicator species. Natural and man‐made landscape features are known to form effective genetic barriers to amphibian populations; however, amphibians with different characteristics may have different species–landscape interaction patterns. We conducted a comparative landscape genetic analysis of two closely related syntopic frog species from central China, Pelophylax nigromaculatus (PN) and Fejervarya limnocharis (FL). These two species differ in several key life history traits; PN has a larger body size and larger clutch size, and reaches sexual maturity later than FL. Microsatellite DNA data were collected and analyzed using conventional (F_ST, isolation by distance (IBD), AMOVA) and recently developed (Bayesian assignment test, isolation by resistance) landscape genetic methods. As predicted, a higher level of population structure in FL (F_ST′ = 0.401) than in PN (F_ST′ = 0.354) was detected, in addition to FL displaying strong IBD patterns (r = .861) unlike PN (r = .073). A general north–south break in FL populations was detected, consistent with the IBD pattern, while PN exhibited clustering of northern‐ and southern‐most populations, suggestive of altered dispersal patterns. Species‐specific resistant landscape features were also identified, with roads and land cover the main cause of resistance to FL, and elevation the main influence on PN. These different species–landscape interactions can be explained mostly by their life history traits, revealing that closely related and ecologically similar species have different responses to the same landscape features. Comparative landscape genetic studies are important in detecting such differences and refining generalizations about amphibians in monitoring environmental changes.

Host Ecology Rather Than Host Phylogeny Drives Amphibian Skin Microbial Community Structure in the Biodiversity Hotspot of Madagascar

Host-associated microbiotas of vertebrates are diverse and complex communities that contribute to host health. In particular, for amphibians, cutaneous microbial communities likely play a significant role in pathogen defense; however, our ecological understanding of these communities is still in its infancy. Here, we take advantage of the fully endemic and locally species-rich amphibian fauna of Madagascar to investigate the factors structuring amphibian skin microbiota on a large scale. Using amplicon-based sequencing, we evaluate how multiple host species traits and site factors affect host bacterial diversity and community structure. Madagascar is home to over 400 native frog species, all of which are endemic to the island; more than 100 different species are known to occur in sympatry within multiple rainforest sites. We intensively sampled frog skin bacterial communities, from over 800 amphibians from 89 species across 30 sites in Madagascar during three field visits, and found that skin bacterial communities differed strongly from those of the surrounding environment. Richness of bacterial operational taxonomic units (OTUs) and phylogenetic diversity differed among host ecomorphs, with arboreal frogs exhibiting lower richness and diversity than terrestrial and aquatic frogs. Host ecomorphology was the strongest factor influencing microbial community structure, with host phylogeny and site parameters (latitude and elevation) explaining less but significant portions of the observed variation. Correlation analysis and topological congruency analyses revealed little to no phylosymbiosis for amphibian skin microbiota. Despite the observed geographic variation and low phylosymbiosis, we found particular OTUs that were differentially abundant between particular ecomorphs. For example, the genus Pigmentiphaga (Alcaligenaceae) was significantly enriched on arboreal frogs, Methylotenera (Methylophilaceae) was enriched on aquatic frogs, and Agrobacterium (Rhizobiaceae) was enriched on terrestrial frogs. The presence of shared bacterial OTUs across geographic regions for selected host genera suggests the presence of core microbial communities which in Madagascar, might be driven more strongly by a species’ preference for specific microhabitats than by the physical, physiological or biochemical properties of their skin. These results corroborate that both host and environmental factors are driving community assembly of amphibian cutaneous microbial communities, and provide an improved foundation for elucidating their role in disease resistance.

Shift from slow- to fast-water habitats accelerates lineage and phenotype evolution in a clade of Neotropical suckermouth catfishes (Loricariidae: Hypoptopomatinae)

Identifying habitat characteristics that accelerate organismal evolution is essential to understanding both the origins of life on Earth and the ecosystem properties that are most critical to maintaining life into the future. Searching for these characteristics on a large scale has only recently become possible via advances in phylogenetic reconstruction, time-calibration, and comparative analyses. In this study, we combine these tools with habitat and phenotype data for 105 species in a clade of Neotropical suckermouth catfishes commonly known as cascudinhos. Our goal was to determine whether riverine mesohabitats defined by different flow rates (i.e., pools vs. rapids) and substrates (plants vs. rocks) have affected rates of cascudinho cladogenesis and morphological diversification. In contrast to predictions based on general theory related to life in fast-flowing, rocky riverine habitats, Neoplecostomini lineages associated with these habitats exhibited increased body size, head shape diversity, and lineage and phenotype diversification rates. These findings are consistent with a growing understanding of river rapids as incubators of biological diversification and specialization. They also highlight the urgent need to conserve rapids habitats throughout the major rivers of the world.

Transcriptomic and macroevolutionary evidence for phenotypic uncoupling between frog life history phases

Anuran amphibians undergo major morphological transitions during development, but the contribution of their markedly different life-history phases to macroevolution has rarely been analysed. Here we generate testable predictions for coupling versus uncoupling of phenotypic evolution of tadpole and adult life-history phases, and for the underlying expression of genes related to morphological feature formation. We test these predictions by combining evidence from gene expression in two distantly related frogs, Xenopus laevis and Mantidactylus betsileanus, with patterns of morphological evolution in the entire radiation of Madagascan mantellid frogs. Genes linked to morphological structure formation are expressed in a highly phase-specific pattern, suggesting uncoupling of phenotypic evolution across life-history phases. This gene expression pattern agrees with uncoupled rates of trait evolution among life-history phases in the mantellids, which we show to have undergone an adaptive radiation. Our results validate a prevalence of uncoupling in the evolution of tadpole and adult phenotypes of frogs.

Linking speciation to extinction: Diversification raises contemporary extinction risk in amphibians

Many of the traits associated with elevated rates of speciation, including niche specialization and having small and isolated populations, are similarly linked with an elevated risk of extinction. This suggests that rapidly speciating lineages may also be more extinction prone. Empirical tests of a speciation-extinction correlation are rare because assessing paleontological extinction rates is difficult. However, the modern biodiversity crisis allows us to observe patterns of extinction in real time, and if this hypothesis is true then we would expect young clades that have recently diversified to have high contemporary extinction risk. Here, we examine evolutionary patterns of modern extinction risk across over 300 genera within one of the most threatened vertebrate classes, the Amphibia. Consistent with predictions, rapidly diversifying amphibian clades also had a greater share of threatened species. Curiously, this pattern is not reflected in other tetrapod classes and may reflect a greater propensity to speciate through peripheral isolation in amphibians, which is partly supported by a negative correlation between diversification rate and mean geographic range size. This clustered threat in rapidly diversifying amphibian genera means that protecting a small number of species can achieve large gains in preserving amphibian phylogenetic diversity. Nonindependence between speciation and extinction rates has many consequences for patterns of biodiversity and how we may choose to conserve it.

Chemical antipredator defence is linked to higher extinction risk

Many attributes of species may be linked to contemporary extinction risk, though some such traits remain untested despite suggestions that they may be important. Here, I test whether a trait associated with higher background extinction rates, chemical antipredator defence, is also associated with current extinction risk, using amphibians as a model system-a group facing global population declines. I find that chemically defended species are approximately 60% more likely to be threatened than species without chemical defence, although the severity of the contemporary extinction risk may not relate to chemical defence. The results confirm that background and contemporary extinction rates can be predicted from the same traits, at least in certain cases. This suggests that associations between extinction risk and phenotypic traits can be temporally stable over long periods. The results also provide novel insights into the relevance of antipredator defences for species subject to conservation concerns.

High evolutionary constraints limited adaptive responses to past climate changes in toad skulls

Interactions among traits that build a complex structure may be represented as genetic covariation and correlation. Genetic correlations may act as constraints, deflecting the evolutionary response from the direction of natural selection. We investigated the relative importance of drift, selection, and constraints in driving skull divergence in a group of related toad species. The distributional range of these species encompasses very distinct habitats with important climatic differences and the species are primarily distinguished by differences in their skulls. Some parts of the toad skull, such as the snout, may have functional relevance in reproductive ecology, detecting water cues. Thus, we hypothesized that the species skull divergence was driven by natural selection associated with climatic variation. However, given that all species present high correlations among skull traits, our second prediction was of high constraints deflecting the response to selection. We first extracted the main morphological direction that is expected to be subjected to selection by using within- and between-species covariance matrices. We then used evolutionary regressions to investigate whether divergence along this direction is explained by climatic variation between species. We also used quantitative genetics models to test for a role of random drift versus natural selection in skull divergence and to reconstruct selection gradients along species phylogeny. Climatic variables explained high proportions of between-species variation in the most selected axis. However, most evolutionary responses were not in the direction of selection, but aligned with the direction of allometric size, the dimension of highest phenotypic variance in the ancestral population. We conclude that toad species have responded to selection related to climate in their skulls, yet high evolutionary constraints dominated species divergence and may limit species responses to future climate change.