Genetic divergence within frog species is greater in topographically more 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.

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.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Amphibian Speciation Rates Support a General Role of Mountains as Biodiversity Pumps

AbstractContinental mountain areas cover <15% of global land surface, yet these regions concentrate >80% of global terrestrial diversity. One prominent hypothesis to explain this pattern proposes that high mountain diversities could be explained by higher diversification rates in regions of high topographic complexity (HTC). While high speciation in mountains has been detected for particular clades and regions, the global extent to which lineages experience faster speciation in mountains remains unknown. Here we addressed this issue using amphibians as a model system (>7,000 species), and we found that families showing high speciation rates contain a high proportion of species distributed in mountains. Moreover, we found that lineages inhabiting areas of HTC speciate faster than lineages occupying areas that are topographically less complex. When comparing across regions, we identified the same pattern in five biogeographical realms where higher speciation rates are associated with higher levels of complex topography. Low-magnitude differences in speciation rates between some low and high complex topographies suggest that high mountain diversity is also affected by low extinction and/or high colonization rates. Nevertheless, our results bolster the importance of mountains as engines of speciation at different geographical scales and highlight their importance for the conservation of global biodiversity.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

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.

Patterns of speciation are similar across mountainous and lowland regions for a Neotropical plant radiation (Costaceae: Costus)

High species richness and endemism in tropical mountains are recognized as major contributors to the latitudinal diversity gradient. The processes underlying mountain speciation, however, are largely untested. The prevalence of steep ecogeographic gradients and the geographic isolation of populations by topographic features are predicted to promote speciation in mountains. We evaluate these processes in a species-rich Neotropical genus of understory herbs that range from the lowlands to montane forests and have higher species richness in topographically complex regions. We ask whether climatic niche divergence, geographic isolation, and pollination shifts differ between mountain-influenced and lowland Amazonian sister pairs inferred from a 756-gene phylogeny. Neotropical Costus ancestors diverged in Central America during a period of mountain formation in the last 3 million years with later colonization of Amazonia. Although climatic divergence, geographic isolation, and pollination shifts are prevalent in general, these factors do not differ between mountain-influenced and Amazonian sister pairs. Despite higher climatic niche and species diversity in the mountains, speciation modes in Costus appear similar across regions. Thus, greater species richness in tropical mountains may reflect differences in colonization history, diversification rates, or the prevalence of rapidly evolving plant life forms, rather than differences in speciation mode.

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.

Anurans of the Caparaó National Park and surroundings, southeast Brazil

Abstract The Atlantic Forest (AF) is one of the biodiversity hotspots of the world, and the most fragmented biome of Brazil. This biome includes different phytophysiognomies, as riparian, slope, cloudy forests, and grasslands. Such complexity, allied to huge latitudinal and high elevational range, provides diverse habitats and conditions for amphibian speciation. As a result, there are over 600 amphibian species known to occur in the AF. Within this biome the Caparaó National Park (CNP) is relevant, as it includes the highest peak of the biome, the Pico da Bandeira at almost 3,000 m above sea level, as well as different phytophysiognomies as rocky fields and humid forests. In spite of that, its amphibian fauna is still poorly described. We inventoried amphibians at the CNP and surrounding areas from 2016 to 2018 and recorded 47 anuran species, of which two are locally endemic and at least six have not been described yet. Additionally, we compiled data from previous surveys (2004 to 2008) and secondary data from scientific collections. All together, we registered a total of 61 anuran species from 12 families for the CNP and surroundings, placing this area among the 10 amphibian richest sites in the AF. Some of these species are represented by only one or two collected specimens and have not been registered in the CNP since the 1980’s, such as Thoropa lutzi and Hylodes vanzolinii. These species could be examples of population declines or even past local extinctions, highlighting the need of further sampling efforts in that highly biodiverse site.

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.

Everything you always wanted to know about gene flow in tropical landscapes (but were afraid to ask)

The bulk of the world’s biodiversity is found in tropical regions, which are increasingly threatened by the human-led degradation of natural habitats. Yet, little is known about tropical biodiversity responses to habitat loss and fragmentation. Here we review all available literature assessing landscape effects on gene flow in tropical species, aiming to help unravel the factors underpinning functional connectivity in the tropics. We map and classify studies by focus species, the molecular markers employed, statistical approaches to assess landscape effects on gene flow, and the evaluated landscape and environmental variables. We then compare qualitatively and quantitatively landscape effects on gene flow across species and units of analysis. We found 69 articles assessing landscape effects on gene flow in tropical organisms, most of which were published in the last five years, were concentrated in the Americas, and focused on amphibians or mammals. Most studies employed population-level approaches, microsatellites were the preferred type of markers, and Mantel and partial Mantel tests the most common statistical approaches used. While elevation, land cover and forest cover were the most common gene flow predictors assessed, habitat suitability was found to be a common predictor of gene flow. A third of all surveyed studies explicitly assessed the effect of habitat degradation, but only 14 of these detected a reduced gene flow with increasing habitat loss. Elevation was responsible for most significant microsatellite-based isolation by resistance effects and a single study reported significant isolation by non-forested areas in an ant. Our study reveals important knowledge gaps on the study of landscape effects on gene flow in tropical organisms, and provides useful guidelines on how to fill them.

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.

DNA Barcoding Survey of Anurans across the Eastern Cordillera of Colombia and the Impact of the Andes on Cryptic Diversity

Colombia hosts the second highest amphibian species diversity on Earth, yet its fauna remains poorly studied, especially using molecular genetic techniques. We present the results of the first wide-scale DNA barcoding survey of anurans of Colombia, focusing on a transect across the Eastern Cordillera. We surveyed 10 sites between the Magdalena Valley to the west and the eastern foothills of the Eastern Cordillera, sequencing portions of the mitochondrial 16S ribosomal RNA and cytochrome oxidase subunit 1 (CO1) genes for 235 individuals from 52 nominal species. We applied two barcode algorithms, Automatic Barcode Gap Discovery and Refined Single Linkage Analysis, to estimate the number of clusters or “unconfirmed candidate species” supported by DNA barcode data. Our survey included ~7% of the anuran species known from Colombia. While barcoding algorithms differed slightly in the number of clusters identified, between three and ten nominal species may be obscuring candidate species (in some cases, more than one cryptic species per nominal species). Our data suggest that the high elevations of the Eastern Cordillera and the low elevations of the Chicamocha canyon acted as geographic barriers in at least seven nominal species, promoting strong genetic divergences between populations associated with the Eastern Cordillera.

Evidence for an intrinsic factor promoting landscape genetic divergence in Madagascan leaf-litter frogs

The endemic Malagasy frog radiations are an ideal model system to study patterns and processes of speciation in amphibians. Large-scale diversity patterns of these frogs, together with other endemic animal radiations, led to the postulation of new and the application of known hypotheses of species diversification causing diversity patterns in this biodiversity hotspot. Both extrinsic and intrinsic factors have been studied in a comparative framework, with extrinsic factors usually being related to the physical environment (landscape, climate, river catchments, mountain chains), and intrinsic factors being clade-specific traits or constraints (reproduction, ecology, morphology, physiology). Despite some general patterns emerging from such large-scale comparative analyses, it became clear that the mechanism of diversification in Madagascar may vary among clades, and may be a multifactorial process. In this contribution, I test for intrinsic factors promoting population-level divergence within a clade of terrestrial, diurnal leaf-litter frogs (genus Gephyromantis) that has previously been shown to diversify according to extrinsic factors. Landscape genetic analyses of the microendemic species Gephyromantis enki and its widely distributed, larger sister species Gephyromantis boulengeri over a rugged landscape in the Ranomafana area shows that genetic variance of the smaller species cannot be explained by landscape resistance alone. Both topographic and riverine barriers are found to be important in generating this divergence. This case study yields additional evidence for the probable importance of body size in lineage diversification.

Geographic determinants of gene flow in two sister species of tropical Andean frogs

Complex interactions between topographic heterogeneity, climatic and environmental gradients, and thermal niche conservatism are commonly assumed to indicate the degree of biotic diversification in montane regions. Our aim was to investigate factors that disrupt gene flow between populations and to determine if there is evidence of downslope asymmetric migration in highland frogs with wide elevational ranges and thermal niches. We determined the role of putative impediments to gene flow (as measured by least-cost path (LCP) distances, topographic complexity, and elevational range) in promoting genetic divergence between populations of 2 tropical Andean frog sister species (Dendropsophus luddeckei, N = 114; Dendropsophus labialis, N = 74) using causal modeling and multiple matrix regression. Although the effect of geographic features was species specific, elevational range and LCP distances had the strongest effect on gene flow, with mean effect sizes (Mantel r and regression coefficients β), between 5 and 10 times greater than topographic complexity. Even though causal modeling and multiple matrix regression produced congruent results, the latter provided more information on the contribution of each geographic variable. We found moderate support for downslope migration. We conclude that the climatic heterogeneity of the landscape, the elevational distance between populations, and the inability to colonize suboptimal habitats due to thermal niche conservatism influence the magnitude of gene flow. Asymmetric migration, however, seems to be influenced by life history traits.