Genetic divergence in tropical anurans: deeper phylogeographic structure in forest specialists and in topographically complex regions

Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas

Despite their limited vagility and pronounced habitat heterogeneity in the tropics, many anuran species have unexpectedly extensive geographic ranges. One prominent example of this phenomenon is Pithecopus hypochondrialis, which is found in the Cerrado, Guianan savanna, and Llanos domains, as well as isolated tracts of savanna and open habitat within the Amazon Forest. The present study employs an integrative species delimitation approach to test the hypothesis that P. hypochondrialis is in fact a species complex. We also reconstruct the relationships among the lineages delimited here and other Pithecopus species. In this study, we employ Ecological Niche Modelling (ENM) and spatiotemporal phylogeographic reconstruction approaches to evaluate a multitude of scenarios of connectivity across the Neotropical savannas. We identified three divergent lineages, two of which have been described previously. The lineages were allocated to a lowland Pithecopus clade, although the relationships among these lineages are weakly supported. Both the ENM and the phylogeographic reconstruction highlight the occurrence of periods of connectivity among the Neotropical savannas over the course of the Pliocene and Pleistocene epochs. These processes extended from eastern Amazonia to the northern coast of Brazil. The findings of the present study highlight the presence of hidden diversity within P. hypochondrialis, and reinforce the need for a comprehensive taxonomic review. These findings also indicate intricate and highly dynamic patterns of connectivity across the Neotropical savannas that date back to the Pliocene.

Landscape genetics reveals contrasting patterns of connectivity in two newt species (Lissotriton montandoni and L. vulgaris)

Ecologically distinct species may respond to landscape changes in different ways. In addition to basic ecological data, the extent of the geographic range has been successfully used as an indicator of species sensitivity to anthropogenic landscapes, with widespread species usually found to be less sensitive compared to range-restricted species. In this study, we investigate connectivity patterns of two closely related but ecologically distinct newt species - the range-restricted, Lissotriton montandoni and the widespread, L. vulgaris - using genomic data, a highly replicated setting (six geographic regions per species), and tools from landscape genetics. Our results show the importance of forest for connectivity in both species, but at the same time suggest differential use of forested habitat, with L. montandoni and L. vulgaris showing the highest connectivity at forest-core and forest-edges, respectively. Anthropogenic landscapes (i.e., higher crop- or urban-cover) increased resistance in both species, but the effect was one to three orders of magnitude stronger in L. montandoni than in L. vulgaris. This result is consistent with a view of L. vulgaris as an ecological generalist. Even so, currently, the negative impact of anthropogenic landscapes is mainly seen in connectivity among L. vulgaris populations, which show significantly stronger isolation and lower effective sizes relative to L. montandoni. Overall, this study emphasizes how habitat destruction is compromising genetic connectivity not only in endemic, range-restricted species of conservation concern but also in widespread generalist species, despite their comparatively lower sensitivity to anthropogenic landscape changes.

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.

Unravelling the cryptic diversity and evolution of the dwarf swamp frog Pseudopaludicola mystacalis (Anura, Leptodactylidae) in open habitats of South America

Many studies on the diversity of the South American biota support the role of ecological and geological events as main drivers of species diversification. For many groups, geomorphological events are the key drivers of diversification, while the influence of Pleistocene climate oscillations is prominent for others. To precisely indicate which events were key for the development of the astonishing biodiversity in South America, studies on widely distributed species are paramount. One such species, the dwarf swamp frog Pseudopaludicola mystacalis (Leptodactylidae, Leiuperinae), is widely distributed in open habitats of South America and we herein investigate population differentiation and diversification in this species across its geographic range. We sequenced a 1374 bp mtDNA fragment from 64 specimens across 25 localities. We used population assignment and species delimitation methods to assess genetic structure and lineage limits across the species distribution. We estimated, for each lineage, intraspecific diversity, divergence times, and demographic histories. Our results recovered ten lineages with up to 5% of genetic divergence among them. Diversification occurred mainly during the Tertiary, suggesting that Miocene-Pliocene topographic events had a major influence on the diversification of P. mystacalis. Pleistocene climatic oscillations also played a role on evolutionary history of P. mystacalis, causing demographic changes in one lineage.

Temporal and spatial diversification along the Amazonia-Cerrado transition in Neotropical treefrogs of the Boana albopunctata species group

Despite extensive research on biodiversity in Neotropical forests, biodiversity in seasonally dry, open biomes in South America has been underestimated until recently. We leverage a widespread group, Boana albopunctata, to uncover cryptic lineages and investigate the timing of diversification in Neotropical anurans with a focus on dry diagonal biomes (Cerrado, Caatinga and Chaco) and the ecotone between Amazonia and the Cerrado. We inferred a multilocus phylogeny of the B. albopunctata species group that includes 15 of 18 described species, recovered two cryptic species, and reconstructed the timing of diversification among species distributed across multiple South American biomes. One new potential species (B. aff. steinbachi), sampled in the Amazonian state of Acre, clustered within the B. calcara-fasciata species complex and is close to B. steinbachi. A second putative new species (B. aff. multifasciata), sampled in the Amazonia-Cerrado ecotone, is closely related to B. multifasciata. Lastly, we place a recently identified Cerrado lineage (B. aff. albopuncata) into the B. albopunctata species group phylogeny for the first time. Our ancestral range reconstruction showed that species in the B. albopuctata group likely dispersed from Amazonia-Cerrado into the dry-diagonal and Atlantic Forest. Intraspecies demography showed, for both B. raniceps and B. albopunctata, signs of rapid expansion across the dry diagonal. Similarly, for one clade of B. multifasciata, our analyses support an invasion of the Cerrado from Amazonia, followed by a rapid expansion across the open diagonal biomes. Thus, our study recovers several recent divergences along the Amazonia-Cerrado ecotone in northern Brazil. Tectonic uplift and erosion in the late Miocene and climate oscillations in the Pleistocene corresponded with estimated divergence times in the dry diagonal and Amazonia-Cerrado ecotone. Our study highlights the importance of these threatened open formations in the generation of biodiversity in the Neotropics.

Deep Genomic Divergence and Phenotypic Admixture of the Treefrog Dendropsophus elegans (Hylidae: Amphibia) Coincide With Riverine Boundaries at the Brazilian Atlantic Forest

The Atlantic Forest (AF) domain is one of the Earth’s biodiversity hotspots, known for its high levels of species diversity and endemism. Factors related to palaeoenvironmental dynamics, such as the establishment of vegetational refugia and river basins, have different impacts on biological communities and biodiversity patterns in this domain. Here, we sample genome-wide RADseq data from a widespread treefrog (Dendropsophus elegans), inhabiting natural and human-impacted ecosystems at the Brazilian AF to test the impact of riverine boundaries and climatic refugia on population structure and diversification. We estimate divergence times and migration rate across identified genetic breaks related to the rivers Doce, Paraíba do Sul, Ribeira de Iguape, and Paraguaçu, known to represent barriers to gene flow for other AF endemic species, and test the role of climatic refugia. Finally, we investigate the impact of spatio-temporal population history on morphological variation in this species. We recovered a phylogeographic history supporting three distinct clades separated into two geographically structured populations, corresponding to the north and south of AF. In addition, we identified an admixture zone between north and south populations in the latitude close to the Doce River. Our findings support a pattern of isolation-by-distance and the existence of a secondary contact zone between populations, which might have been promoted by gene flow during population expansion. Further, we found support for models considering migration parameters for all the tested rivers with different population divergence times. Based on the species history and the AF palaeoenvironmental dynamics, we corroborate the role of forest refugia impacting population structure for this species through recent range expansion after the Last Glacial Maximum (LGM). The Doce and Paraíba do Sul Rivers coincide with the main genetic breaks, suggesting they might also have played a role in the diversification processes. Finally, despite finding subtle correlations for phenotypic data among different populations, variation is not strongly detectable and does not seem associated with speciation-level processes that could warrant taxonomic changes. Such results can be explained by phenotypic plasticity of the evaluated traits and by recent divergence times, where there has been insufficient time and weak selective pressures to accumulate enough phenotypic differences.

Ecological divergence and synchronous Pleistocene diversification in the widespread South American butter frog complex

Phylogeographic studies primarily focus on the major role of landscape topography in driving lineage diversification. However, populational phylogeographic breaks may also occur as a result of either niche conservatism or divergence, in the absence of geographic barriers to gene flow. Furthermore, these two factors are not mutually exclusive and can act in concert, making it challenging to evaluate their relative importance on explaining genetic variation in nature. Herein, we use sequences of two mitochondrial and four nuclear genes to investigate the timing and diversification patterns of species pertaining to the Leptodactylus latrans complex, which harbors four morphologically cryptic species with broad distributions across environmental gradients in eastern South America. The origin of this species complex dates back to the late Miocene (ca. 5.5 Mya), but most diversification events occurred synchronically during the late Pleistocene likely as the result of ecological divergence driven by Quaternary climatic oscillations. Further, significant patterns of environmental niche divergences among species in the L. latrans complex imply that ecological isolation is the primary mode of genetic diversification, mostly because phylogenetic breaks are associated with environmental transitions rather than topographic barriers at both species and populational scale. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America.

Subtle environmental variation affects phenotypic differentiation of shallow divergent treefrog lineages in Amazonia

Amazonia harbours a vast biotic and ecological diversity, enabling investigation of the effects of microevolutionary processes and environmental variation on species diversification. Integrative approaches combining phenotypic and genetic variation can improve our knowledge on diversification processes in megadiverse regions. Here, we investigate the influence of environmental and geographic variation on the genetic and morphological differentiation in the Amazonian Boana calcarata-fasciata (Anura: Hylidae) species complex. We analysed the variation of one mtDNA gene from individuals of different forest environments, and assessed their phylogenetic relationships and species limits to define the lineages to perform a phenotypic-environmental approach. We collected morphological data (head shape and size) using 3D models and investigated the phylogenetic signal, evolutionary model and influence of environmental variables on morphology. We verified associations between environmental and geographical distances with morphological and genetic variation using distance-based redundancy analyses and Mantel tests. We found an even higher cryptic diversity than already recognized within the species complex. Body size and head shape varied among specimens, but did not present phylogenetic signal, diverging under a selective evolutionary model. Our results show that diverse factors have influenced morphological and genetic variation, but environmental conditions such as vegetation cover, precipitation and climate change velocity influenced morphological diversification. Possible population-level mechanisms such as parallel morphological evolution or plastic responses to similar environments could account for such patterns in these typical Amazonian treefrogs.

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.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

Predicting amphibian intraspecific diversity with machine learning: Challenges and prospects for integrating traits, geography, and genetic data

The growing availability of genetic data sets, in combination with machine learning frameworks, offers great potential to answer long-standing questions in ecology and evolution. One such question has intrigued population geneticists, biogeographers, and conservation biologists: What factors determine intraspecific genetic diversity? This question is challenging to answer because many factors may influence genetic variation, including life history traits, historical influences, and geography, and the relative importance of these factors varies across taxonomic and geographic scales. Furthermore, interpreting the influence of numerous, potentially correlated variables is difficult with traditional statistical approaches. To address these challenges, we analysed repurposed data using machine learning and investigated predictors of genetic diversity, focusing on Nearctic amphibians as a case study. We aggregated species traits, range characteristics, and >42,000 genetic sequences for 299 species using open-access scripts and various databases. After identifying important predictors of nucleotide diversity with random forest regression, we conducted follow-up analyses to examine the roles of phylogenetic history, geography, and demographic processes on intraspecific diversity. Although life history traits were not important predictors for this data set, we found significant phylogenetic signal in genetic diversity within amphibians. We also found that salamander species at northern latitudes contained low genetic diversity. Data repurposing and machine learning provide valuable tools for detecting patterns with relevance for conservation, but concerted efforts are needed to compile meaningful data sets with greater utility for understanding global biodiversity.

Comparison of the mitochondrial phylogeographical structure of a generalist and two specialist frog species reveals contrasting patterns in the Eastern and Western Cape provinces of South Africa

In this study, we examined the phylogeographical structure of three frog species (Anhydrophryne rattrayi, Arthroleptis wageri and Cacosternum nanum) in the Eastern and Western Cape provinces of South Africa. The first two species are forest dwelling and exhibit direct development, whereas the last species is a habitat generalist, breeding in open freshwater bodies and exhibiting rapid metamorphosis. Evolutionary relationships were inferred using combined mitochondrial DNA (mtDNA) sequence data [16S rRNA and cytochrome b (Cytb)]. Divergence times were estimated for each species using the combined mtDNA dataset, and population genetic structuring was inferred using haplotype networks and analysis of molecular variance using the rapidly evolving Cytb locus. The two forest-dwelling species (Anhydrophryne rattrayi and Arthroleptis wageri) each revealed two statistically well-supported clades and were characterized by marked genetic differentiation and the general absence of shared maternal haplotypes, indicating no maternal gene flow between conspecific populations; a result corroborated by moderate values of pairwise genetic distance. In contrast, for the generalist species C. nanum, two shallow clades were observed and several haplotypes were shared between localities, indicating moderate maternal gene flow. Diversification in the two forest-dwelling species occurred during the Plio-Pleistocene climatic oscillations and was associated with increased xeric conditions, whereas in C. nanum a recent, Holocene divergence was inferred. Ancient climatic factors thought to be causal in the divergences within each of the three species are discussed. The marked mtDNA differentiation in the two forest specialist species suggests the presence of distinct management units that should be considered in future conservation management decisions.

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.

Allometric escape from acoustic constraints is rare for frog calls

Allometric constraint is a product of natural selection and physical laws, particularly with respect to body size and traits constrained by properties thereof, such as metabolism, longevity, and vocal frequency. Allometric relationships are often conserved across lineages, indicating that physical constraints dictate scaling patterns in deep time, despite substantial genetic and ecological divergence among organisms. In particular, acoustic allometry (sound frequency ~ body size) is conserved across frogs, in defiance of massive variation in both body size and frequency. Here, we ask how many instances of allometric escape have occurred across the frog tree of life using a Bayesian framework that estimates the location, number, and magnitude of shifts in the adaptive landscape of acoustic allometry. Moreover, we test whether ecology in terms of calling site could affect these relationships. We find that calling site has a major influence on acoustic allometry. Despite this, we identify only four major instances of allometric escape, potentially deriving from ecomorphological adaptations to new signal modalities. In these instances of allometric escape, the optima and strength of the scaling relationship are different than expected for most other frog species, representing new adaptive regimes of body size ~ call frequency. Allometric constraints on frog calls are highly conserved and have rarely allowed escape, despite frequent invasions of new adaptive regimes and dramatic ecomorphological divergence. Our results highlight the rare instances in which natural and sexual selection combined can overcome physical constraints on sound production.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

Functional traits explain amphibian distribution in the Brazilian Atlantic Forest

AIM

Species distributions are one of the most important ways to understand how communities interact through macroecological relationships. The functional abilities of a species, such as its plasticity in various environments, can determine its distribution, species richness and beta diversity patterns. In this study, we evaluate how functional traits influence the distribution of amphibians, and hypothesize which functional traits explain the current pattern of amphibian species composition.

LOCATION

Atlantic Forest, Brazil.

TAXON

Amphibia (Anura and Gymnophiona).

METHODS

Using potential distributions of Brazilian amphibians from Atlantic Forest based on their functional traits, we analysed the influence of biotic and abiotic factors on species richness, endemism (with permutation multivariate analysis) and beta diversity components (i.e. total, turnover and nestedness dissimilarities).

RESULTS

Environmental variables explained 59.5% of species richness, whereas functional traits explained 15.8% of species distribution (geographical species range) for Anuran and 88.8% for Gymnophiona. Body size had the strongest correlation with species distribution. Results showed that species with medium to large body size, and species that are adapted to living in open areas tended to disperse from west to east direction. Current forest changes directly affected beta diversity patterns (i.e. most species adapted to novel environments increase their ranges). Beta diversity partitioning between humid and dry forests showed decreased nestedness and increased turnover by increasing altitude in the south-eastern region of the Atlantic Forest.

MAIN CONCLUSIONS

Our study shows that functional traits directly influence the ability of the species to disperse. With the alterations of the natural environment, species more apt to these alterations have dispersed or increased their distribution, which consequently changes community structure. As a result, there are nested species distribution patterns and homogenization of amphibian species composition throughout the Brazilian Atlantic Forest.

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.

Morphological and ecological convergence at the lower size limit for vertebrates highlighted by five new miniaturised microhylid frog species from three different Madagascan genera

Miniaturised frogs form a fascinating but poorly understood amphibian ecomorph and have been exceptionally prone to taxonomic underestimation. The subfamily Cophylinae (family Microhylidae), endemic to Madagascar, has a particularly large diversity of miniaturised species which have historically been attributed to the single genus Stumpffia largely based on their small size. Recent phylogenetic work has revealed that several independent lineages of cophyline microhylids evolved towards highly miniaturised body sizes, achieving adult snout-vent lengths under 16 mm. Here, we describe five new species belonging to three clades that independently miniaturised and that are all genetically highly divergent from their relatives: (i) a new genus (Mini gen. nov.) with three new species from southern Madagascar, (ii) one species of Rhombophryne, and (iii) one species of Anodonthyla. Mini mum sp. nov. from Manombo in eastern Madagascar is one of the smallest frogs in the world, reaching an adult body size of 9.7 mm in males and 11.3 mm in females. Mini scule sp. nov. from Sainte Luce in southeastern Madagascar is slightly larger and has maxillary teeth. Mini ature sp. nov. from Andohahela in southeast Madagascar is larger than its congeners but is similar in build. Rhombophryne proportionalis sp. nov. from Tsaratanana in northern Madagascar is unique among Madagascar's miniaturised frogs in being a proportional dwarf, exhibiting far less advanced signs of paedomorphism than other species of similar size. Anodonthyla eximia sp. nov. from Ranomafana in eastern Madagascar is distinctly smaller than any of its congeners and is secondarily terrestrial, providing evidence that miniaturisation and terrestriality may be evolutionarily linked. The evolution of body size in Madagascar's microhylids has been more dynamic than previously understood, and future studies will hopefully shed light on the interplay between ecology and evolution of these remarkably diverse frogs.

Diversity of miniaturized frogs of the genus Adelophryne (Anura: Eleutherodactylidae): A new species from the Atlantic Forest of northeast Brazil

The number of species of frogs in the South American genus Adelophryne has increased in recent years, and it has become apparent that this group contains a substantial amount of undescribed diversity. Currently the genus contains nine described species and five candidate species. Here we describe the tenth species of the genus Adelophryne from the municipality of Igrapiúna, southern Bahia state, Brazil. The new species is characterized by its small body size, indistinct tympanum, and two phalanges in the finger IV. The species of the genus are distributed in three groups, Northern Amazonia Clade, Northern Atlantic Forest Clade and Southern Atlantic Forest Clade. The new species is phylogenetically related to species of the Northern Atlantic Forest Clade of Adelophryne and restricted to forested habitat, as typical for other Adelophryne. The species is restricted to the pristine forests in the type locality, and we consider its conservation status as Near Threatened. New morphological and molecular data of other Adelophryne species are presented, extending the distribution of Adelophryne sp. 2, Adelophryne sp. 4, Adelophryne mucronata and Adelophryne glandulata. However, a more comprehensive revision of the diversity and phylogenetic position of most Adelophryne species is needed, and the evolutionary relationships of A. meridionalis and A. pachydactyla remain unknown.

Phenotypes in phylogeography: Species' traits, environmental variation, and vertebrate diversification

Almost 30 y ago, the field of intraspecific phylogeography laid the foundation for spatially explicit and genealogically informed studies of population divergence. With new methods and markers, the focus in phylogeography shifted to previously unrecognized geographic genetic variation, thus reducing the attention paid to phenotypic variation in those same diverging lineages. Although phenotypic differences among lineages once provided the main data for studies of evolutionary change, the mechanisms shaping phenotypic differentiation and their integration with intraspecific genetic structure have been underexplored in phylogeographic studies. However, phenotypes are targets of selection and play important roles in species performance, recognition, and diversification. Here, we focus on three questions. First, how can phenotypes elucidate mechanisms underlying concordant or idiosyncratic responses of vertebrate species evolving in shared landscapes? Second, what mechanisms underlie the concordance or discordance of phenotypic and phylogeographic differentiation? Third, how can phylogeography contribute to our understanding of functional phenotypic evolution? We demonstrate that the integration of phenotypic data extends the reach of phylogeography to explain the origin and maintenance of biodiversity. Finally, we stress the importance of natural history collections as sources of high-quality phenotypic data that span temporal and spatial axes.

Evolutionarily significant units of the critically endangered leaf frog Pithecopus ayeaye (Anura, Phyllomedusidae) are not effectively preserved by the Brazilian protected areas network

Protected areas (PAs) are essential for biodiversity conservation, but their coverage is considered inefficient for the preservation of all species. Many species are subdivided into evolutionarily significant units (ESUs) and the effectiveness of PAs in protecting them needs to be investigated. We evaluated the usefulness of the Brazilian PAs network in protecting ESUs of the critically endangered Pithecopus ayeaye through ongoing climate change. This species occurs in a threatened mountaintop ecosystem known as campos rupestres. We used multilocus DNA sequences to delimit geographic clusters, which were further validated as ESUs with a coalescent approach. Ecological niche modeling was used to estimate spatial changes in ESUs' potential distributions, and a gap analysis was carried out to evaluate the effectiveness of the Brazilian PAs network to protect P. ayeaye in the face of climate changes. We tested the niche overlap between ESUs to gain insights for potential management alternatives for the species. Pithecopus ayeaye contains at least three ESUs isolated in distinct mountain regions, and one of them is not protected by any PA. There are no climatic niche differences between the units, and only 4% of the suitable potential area of the species is protected in present and future projections. The current PAs are not effective in preserving the intraspecific diversity of P. ayeaye in its present and future range distributions. The genetic structure of P. ayeaye could represent a typical pattern in campos rupestres endemics, which should be considered for evaluating its conservation status.

Genetic diversity of Morato's Digger Toad, Proceratophrys moratoi: spatial structure, gene flow, effective size and the need for differential management strategies of populations

The Morato's Digger Toad, Proceratophrys moratoi, is a critically endangered toad species with a marked population decline in southern Brazilian Cerrado. Despite this, new populations are being discovered, primarily in the northern part of the distribution range, which raises a number of questions with regard to the conservation status of the species. The present study analyzed the genetic diversity of the species based on microsatellite markers. Our findings permitted the identification of two distinct management units. We found profound genetic structuring between the southern populations, on the left margin of the Tietê River, and all other populations. A marked reduction was observed in the contemporary gene flow among the central populations that are most affected by anthropogenic impacts, such as extensive sugar cane plantations, which presumably decreases habitat connectivity. The results indicated reduced diversity in the southern populations which, combined with a smaller effective population size, may make these populations more susceptible to extinction. We recommend the reclassification of P. moratoi as vulnerable and the establishment of a special protection program for the southern populations. Our results provide important insights about the local extinction of southern populations of this toad.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.