The role of tropical dry forest as a long-term barrier to dispersal: a comparative phylogeographical analysis of dry forest tolerant and intolerant frogs

Unraveling the complex phylogeographic history of freshwater fishes in Lower Central America: A study of the electric fish Brachyhypopomus occidentalis

Lower Central America (LCA) has a complex biogeographic history shaped by the rise of the Isthmus of Panama and the global climatic oscillations of the Pleistocene. These events have been crucial in structuring biodiversity in LCA, but their consequences for the distribution and partitions of genetic diversity across the region remain to be elucidated. We combined complete mitochondrial genomes and nuclear ultraconserved elements (UCEs) to study the phylogeographic history and population genetic structure of the electric fish Brachyhypopomus occidentalis in LCA. Our results are consistent with the known phylogeographic history of B. occidentalis in LCA, but we update this history in several important ways that help illuminate the phylogeographic history of freshwater fishes in the region. We provide: i) support for three waves of colonization, two of which occurred prior to the final closure of the Panama Isthmus; ii) a more precise understanding of each colonization event, with evidence for a larger footprint of the first event, as well as genetic exchange across the continental divide in subsequent events; and iii) evidence for high levels of previously unrecognized population genetic structure across LCA. This updated model of colonization and diversification of B. occidentalis consists of three waves of dispersal and colonization, which triggered the evolution of geographic breaks in both nuclear and mitochondrial genomes across LCA. These processes are tightly linked to the dynamic uplift of the Isthmus, recent volcanic activity in the region, and the sea-level oscillations of the Pleistocene. These results improve previous phylogeographic inferences regarding the distribution and diversification of freshwater fishes in LCA, and generate testable hypotheses to guide future research exploring the factors shaping biodiversity in the region.

Amphibian diversity across three adjacent ecosystems in Área de Conservación Guanacaste, Costa Rica

Amphibians are the most threatened species-rich vertebrate group, with species extinctions and population declines occurring globally, even in protected and seemingly pristine habitats. These 'enigmatic declines' are generated by climate change and infectious diseases. However, the consequences of these declines are undocumented as no baseline ecological data exists for most affected areas. Like other neotropical countries, Costa Rica, including Área de Conservación Guanacaste (ACG) in north-western Costa Rica, experienced rapid amphibian population declines and apparent extinctions during the past three decades. To delineate amphibian diversity patterns within ACG, a large-scale comparison of multiple sites and habitats was conducted. Distance and time constrained visual encounter surveys characterised species richness at five sites-Murciélago (dry forest), Santa Rosa (dry forest), Maritza (mid-elevation dry-rain forest intersect), San Gerardo (rainforest) and Cacao (cloud forest). Furthermore, species-richness patterns for Cacao were compared with historic data from 1987-8, before amphibians declined in the area. Rainforests had the highest species richness, with triple the species of their dry forest counterparts. A decline of 45% (20 to 11 species) in amphibian species richness was encountered when comparing historic and contemporary data for Cacao. Conservation efforts sometimes focus on increasing the resilience of protected areas, by increasing their range of ecosystems. In this sense ACG is unique containing many tropical ecosystems compressed in a small geographic space, all protected and recognised as a UNESCO world heritage site. It thus provides an extraordinary platform to understand changes, past and present, and the resilience of tropical ecosystems and assemblages, or lack thereof, to climate change.

Cryptic species diversity in a widespread neotropical tree genus: The case of Cedrela odorata

PREMISE

Reconciling the use of taxonomy to partition morphological variation and describe genetic divergence within and among closely related species is a persistent challenge in phylogenetics. We reconstructed phylogenetic relationships among Cedrela odorata (Meliaceae) and five closely allied species to test the genetic basis for the current model of species delimitation in this economically valuable and threatened genus.

METHODS

We prepared a nuclear species tree with the program SNPhylo and 16,000 single-nucleotide polymorphisms from 168 Cedrela specimens. Based on clades present and ancestral patterns ADMIXTURE, we designed nine species delimitation models and compared each model to current taxonomy with Bayes factor delimitation. Timing of major lineage divergences was estimated with the program SNAPP.

RESULTS

The resulting analysis revealed that modern C. odorata evolved from two genetically distinct ancestral sources. All species delimitation models tested better fit the data than the model representing current taxonomic delimitation. Models with the greatest marginal likelihoods separated Mesoamerican C. odorata and South American C. odorata into two species and lumped C. angustifolia and C. montana as a single species. We estimated that Cedrela diversified in South America within the last 19 million years following one or more dispersal events from Mesoamerican lineages.

CONCLUSIONS

Our analyses show that the present taxonomic understanding within the genus obscures divergent lineages in C. odorata due in part to morphological differentiation and taxonomic distinctions that are not predictably associated with genetic divergence. A more accurate application of taxonomy to C. odorata and related species may aid in its conservation, management, and restoration efforts.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

Evolutionary history of the oriental fire-bellied toad (Bombina orientalis) in Northeast China

The evolutionary history of a species is generally affected by the combination of geological events and climate fluctuations. By analyzing the population features, genetic structure and the effective population historical dynamics of existing species, the population evolutionary history can be reestablished. In recent years, geological evidence shows that the Yilan-Yitong fault zone located in Northeast Asia experienced strong and frequent geological changes in the late Quaternary period. Species population history has been shaped by the combination of the complex climatic conditions of the Quaternary and Pleistocene glacial interglacial cycles and palaeogeological events in Northeast Asia and it has become a research focus for evolutionary biology researchers. In this study, mitochondrial and microsatellite molecular markers were used to reveal the population features, genetic structure, and the effective population historical dynamics of the Oriental fire-bellied toad (Bombina orientalis). The results showed that the strong seismic activity of the Yilan-Yitong fault zone in the late Quaternary period was the main reason for the population differentiation of Oriental fire-bellied toad in northeast China. The Quaternary Pleistocene glacial interglacial cycles led to the significant bottleneck effect of the western population located in the Maoer mountain area. As a result, the western population has low genetic diversity. Recent gene flow between eastern and western populations and historical evidence of population expansion proved that the dispersal behavior of the western populations was the main cause of the low genetic diversity and mitochondrial and nuclear discordance. Human economic activity may be the mainly driving factor. These evidences showed that the comprehensive influence of geology, climate, human activities and other factors should be considered in the process of exploring the evolutionary history of species.

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.

Recent and Rapid Radiation of the Highly Endangered Harlequin Frogs (Atelopus) into Central America Inferred from Mitochondrial DNA Sequences

Populations of amphibians are experiencing severe declines worldwide. One group with the most catastrophic declines is the Neotropical genus Atelopus (Anura: Bufonidae). Many species of Atelopus have not been seen for decades and all eight Central American species are considered “Critically Endangered”, three of them very likely extinct. Nonetheless, the taxonomy, phylogeny, and biogeographic history of Central American Atelopus are still poorly known. In this study, the phylogenetic relationships among seven of the eight described species in Central America were inferred based on mitochondrial DNA sequences from 103 individuals, including decades-old museum samples and two likely extinct species, plus ten South American species. Among Central American samples, we discovered two candidate species that should be incorporated into conservation programs. Phylogenetic inference revealed a ladderized topology, placing species geographically furthest from South America more nested in the tree. Model-based ancestral area estimation supported either one or two colonization events from South America. Relaxed-clock analysis of divergence times indicated that Atelopus colonized Central America prior to 4 million years ago (Ma), supporting a slightly older than traditional date for the closure of the Isthmus. This study highlights the invaluable role of museum collections in documenting past biodiversity, and these results could guide future conservation efforts. An abstract in Spanish (Resumen) is available as supplementary material.

Contrasting Patterns of Population History and Seed-mediated Gene Flow in Two Endemic Costa Rican Oak Species

Lower Central America is an important area to study recent population history and diversification of Neotropical species due to its complex and dynamic geology and climate. Phylogeographic studies in this region are few in comparison with other regions and even less for tree species. The aim of the present study was to characterize the phylogeographic structure in 2 partially co-distributed endemic oak species (Quercus costaricensis and Q. bumelioides) of the Costa Rican mountains using chloroplast short sequence repeats (cpSSRs), and to test for the effect of geological and palaeoclimatic processes on their population history. Genetic diversity and structure, haplotype networks, patterns of seed-mediated gene flow and historical demography were estimated for both species. Results suggested contrasting patterns. Quercus costaricensis exhibited high values of genetic diversity, a marked phylogeographic structure, a north-to-south genetic diversity gradient and evidence of a demographic expansion during the Quaternary. Quercus bumelioides did not show significant genetic structure and the haplotype network and historical demography estimates suggested a recent population expansion probably during the Pleistocene-Holocene transition. The phylogeographic structure of Q. costaricensis seems to be related to Pleistocene altitudinal migration due to its higher altitudinal distribution. Meanwhile, historical seed-mediated gene flow through the lower altitudinal distribution of Q. bumelioides may have promoted the homogenization of genetic variation. Population expansion and stable availability of suitable climatic areas in both species probably indicate that palaeoclimatic changes promoted downwards altitudinal migration and formation of continuous forests allowing oak species to expand their distribution into the Panamanian mountains during glacial stages.

Cryptic diversity in Lithobateswarszewitschii (Amphibia, Anura, Ranidae)

Lithobateswarszewitschii is a species of ranid frog distributed from southern Honduras to Panama. This species suffered severe population declines at higher elevations (above 500 m a.s.l.) from the 1980s to early 1990s, but there is more recent evidence of recovery in parts of its range. Here we advocate for the status of Lithobateswarszewitschii as a candidate cryptic species complex based on sequence data from mitochondrial genes CO1 and 16S. Using concatenated phylogenies, nucleotide diversity (K2P-π), net between group mean distance (NBGMD) (πnet) and species delimitation methods, we further elucidate cryptic diversity within this species. All phylogenies display polyphyletic lineages within Costa Rica and Panama. At both loci, observed genetic polymorphism (K2P-π) is also high within and between geographic populations, surpassing proposed species threshold values for amphibians. Additionally, patterns of phylogeographic structure are complicated for this species, and do not appear to be explained by geographic barriers or isolation by distance. These preliminary findings suggest L.warszewitschii is a wide-ranging species complex. Therefore, we propose further research within its wider range, and recommend integrative taxonomic assessment is merited to assess species status.

The role of Central American barriers in shaping the evolutionary history of the northernmost glassfrog, Hyalinobatrachium fleischmanni (Anura: Centrolenidae)

The complex geological history of Central America has been useful for understanding the processes influencing the distribution and diversity of multiple groups of organisms. Anurans are an excellent choice for such studies because they typically exhibit site fidelity and reduced movement. The objective of this work was to identify the impact of recognized geographic barriers on the genetic structure, phylogeographic patterns and divergence times of a wide-ranging amphibian species, Hyalinobatrachium fleischmanni. We amplified three mitochondrial regions, two coding (COI and ND1) and one ribosomal (16S), in samples collected from the coasts of Veracruz and Guerrero in Mexico to the humid forests of Chocó in Ecuador. We examined the biogeographic history of the species through spatial clustering analyses (Geneland and sPCA), Bayesian and maximum likelihood reconstructions, and spatiotemporal diffusion analysis. Our data suggest a Central American origin of H. fleischmanni and two posterior independent dispersals towards North and South American regions. The first clade comprises individuals from Colombia, Ecuador, Panama and the sister species Hyalinobatrachium tatayoi; this clade shows little structure, despite the presence of the Andes mountain range and the long distances between sampling sites. The second clade consists of individuals from Costa Rica, Nicaragua, and eastern Honduras with no apparent structure. The third clade includes individuals from western Honduras, Guatemala, and Mexico and displays deep population structure. Herein, we synthesize the impact of known geographic areas that act as barriers to glassfrog dispersal and demonstrated their effect of differentiating H. fleischmanni into three markedly isolated clades. The observed genetic structure is associated with an initial dispersal event from Central America followed by vicariance that likely occurred during the Pliocene. The southern samples are characterized by a very recent population expansion, likely related to sea-level and climatic oscillations during the Pleistocene, whereas the structure of the northern clade has probably been driven by dispersal through the Isthmus of Tehuantepec and isolation by the Motagua–Polochic–Jocotán fault system and the Mexican highlands.

Testing Hypotheses of Diversification in Panamanian Frogs and Freshwater Fishes Using Hierarchical Approximate Bayesian Computation with Model Averaging

Most Neotropical frog and freshwater fish species sampled to date show phylogeographic breaks along the Pacific coast of the Isthmus of Panama, with lineages in Costa Rica and western Panama isolated from central Panama. We examine temporal patterns of diversification of taxa across this ‘western Panama isthmus’ (WPI) break to test hypotheses about the origin of species geographical distributions and genetic structuring in this region. We tested for synchronous diversification of four codistributed frog taxon-pairs and three fish taxon-pairs sharing the WPI break using hierarchical approximate Bayesian computation with model averaging based on mitochondrial DNA sequences. We also estimated lineage divergence times using full-Bayesian models. Several of our results supported synchronous divergences within the frog and freshwater fish assemblages; however, Bayes factor support was equivocal for or against synchronous or asynchronous diversification. Nevertheless, we infer that frog populations were likely isolated by one or multiple Pliocene–Pleistocene events more recently than predicted by previous models, while fish genetic diversity was structured by Pleistocene events. By integrating our results with external information from geology and elevational sea level modeling, we discuss the implications of our findings for understanding the biogeographical scenario of the diversification of Panamanian frogs and fishes. Consistent with the ‘Bermingham/Martin model’ (Molecular Ecology 1998, 7, 499–517), we conclude that the regional fish assemblage was fractured by processes shaping isthmian landscapes during the Pleistocene glaciations, including drainage basin isolation during lowered sea levels.

Inferring responses to climate dynamics from historical demography in neotropical forest lizards

We apply a comparative framework to test for concerted demographic changes in response to climate shifts in the neotropical lowland forests, learning from the past to inform projections of the future. Using reduced genomic (SNP) data from three lizard species codistributed in Amazonia and the Atlantic Forest (Anolis punctatus, Anolis ortonii, and Polychrus marmoratus), we first reconstruct former population history and test for assemblage-level responses to cycles of moisture transport recently implicated in changes of forest distribution during the Late Quaternary. We find support for population shifts within the time frame of inferred precipitation fluctuations (the last 250,000 y) but detect idiosyncratic responses across species and uniformity of within-species responses across forest regions. These results are incongruent with expectations of concerted population expansion in response to increased rainfall and fail to detect out-of-phase demographic syndromes (expansions vs. contractions) across forest regions. Using reduced genomic data to infer species-specific demographical parameters, we then model the plausible spatial distribution of genetic diversity in the Atlantic Forest into future climates (2080) under a medium carbon emission trajectory. The models forecast very distinct trajectories for the lizard species, reflecting unique estimated population densities and dispersal abilities. Ecological and demographic constraints seemingly lead to distinct and asynchronous responses to climatic regimes in the tropics, even among similarly distributed taxa. Incorporating such constraints is key to improve modeling of the distribution of biodiversity in the past and future.

Phylogeography of the Central American lancehead Bothrops asper (SERPENTES: VIPERIDAE)

The uplift and final connection of the Central American land bridge is considered the major event that allowed biotic exchange between vertebrate lineages of northern and southern origin in the New World. However, given the complex tectonics that shaped Middle America, there is still substantial controversy over details of this geographical reconnection, and its role in determining biogeographic patterns in the region. Here, we examine the phylogeography of Bothrops asper, a widely distributed pitviper in Middle America and northwestern South America, in an attempt to evaluate how the final Isthmian uplift and other biogeographical boundaries in the region influenced genealogical lineage divergence in this species. We examined sequence data from two mitochondrial genes (MT-CYB and MT-ND4) from 111 specimens of B. asper, representing 70 localities throughout the species’ distribution. We reconstructed phylogeographic patterns using maximum likelihood and Bayesian methods and estimated divergence time using the Bayesian relaxed clock method. Within the nominal species, an early split led to two divergent lineages of B. asper: one includes five phylogroups distributed in Caribbean Middle America and southwestern Ecuador, and the other comprises five other groups scattered in the Pacific slope of Isthmian Central America and northwestern South America. Our results provide evidence of a complex transition that involves at least two dispersal events into Middle America during the final closure of the Isthmus.

Dispersal corridors for plant species in the Poyang Lake Basin of southeast China identified by integration of phylogeographic and geospatial data

Measuring the dispersal of wildlife through landscapes is notoriously difficult. Recently, the categorical least cost path algorithm that integrates population genetic data with species distribution models has been applied to reveal population connectivity. In this study, we use this method to identify the possible dispersal corridors of five plant species (Castanopsis tibetana, Schima superba, Cyclocarya paliurus, Sargentodoxa cuneata, Eomecon chionantha) in the Poyang Lake Basin (PLB, largely coinciding with Jiangxi Province), China, in the late Quaternary. The results showed that the strongest population connectivity for the five species occurred in the Wuyi Mountains and the Yu Mountains of the eastern PLB (East Corridor) during the late Quaternary. In the western PLB, populations of the five species were connected by the Luoxiao Mountains and the Jiuling Mountains (West Corridor) but with a lower degree of connectivity. There were some minor connections between the eastern and the western populations across the Gannan Hills. When the corridors of five species were overlaid, the East Corridor and the West Corridor were mostly shared by multiple species. These results indicate that plant species in the PLB could have responded to the Quaternary climate changes by moving along the East Corridor and the West Corridor. Given that dispersal corridors have seldom been considered in the governmental strategies of biodiversity conservation in the PLB, preserving and restoring natural vegetation along these corridors should be prioritized to mitigate the effects of anthropogenic climate change by facilitating migration of plant species and other biota.

Diversification of the rainfrog Pristimantis ornatissimus in the lowlands and Andean foothills of Ecuador

Geographic barriers and elevational gradients have long been recognized as important in species diversification. Here, we illustrate an example where both mechanisms have shaped the genetic structure of the Neotropical rainfrog, Pristimantis ornatissimus, which has also resulted in speciation. This species was thought to be a single evolutionary lineage distributed throughout the Ecuadorian Chocó and the adjacent foothills of the Andes. Based on recent sampling of P. ornatissimus sensu lato, we provide molecular and morphological evidence that support the validity of a new species, which we name Pristimantis ecuadorensis sp. nov. The sister species are elevational replacements of each other; the distribution of Pristimantis ornatissimus sensu stricto is limited to the Ecuadorian Chocó ecoregion (< 1100 m), whereas the new species has only been found at Andean localities between 1450-1480 m. Given the results of the Multiple Matrix Regression with Randomization analysis, the genetic difference between P. ecuadorensis and P. ornatissimus is not explained by geographic distance nor environment, although environmental variables at a finer scale need to be tested. Therefore this speciation event might be the byproduct of stochastic historic extinction of connected populations or biogeographic events caused by barriers to dispersal such as rivers. Within P. ornatissimus sensu stricto, morphological patterns and genetic structure seem to be related to geographic isolation (e.g., rivers). Finally, we provide an updated phylogeny for the genus, including the new species, as well as other Ecuadorian Pristimantis.

Phylogenetic relationships within the Callicebus cupreus species group (Pitheciidae: Primates): Biogeographic and taxonomic implications

The genus Callicebus (Thomas, 1903) is one of the most diverse of Neotropical primate genera and the only extant member of the Callicebinae subfamily. It has a widespread distribution from Colombia to Brazil, Bolivia, Peru and northern Paraguay. Coat colouring and colour pattern vary substantially within the genus, and this has led to the description of numerous species and subspecies, as well as numerous species groups. However, a lack of molecular phylogenetic analyses on the genus means that phylogenetic relationships and biogeographic history of species are poorly understood. Here, we examined phylogenetic relationships and patterns of diversification within the Callicebus cupreus species Group (sensu Kobayashi, 1995) using complete mitochondrial DNA cytochrome b gene sequence. Analyses indicate that the Callicebus cupreus Group underwent recent and extensive diversification. The common ancestor appears to have emerged some 2.3 million years ago (Ma) from a centre of origin in the western Amazon region, followed by diversification of the group between about 1.5 and 1.2Ma. Phylogenetic analyses were able to recover most previously described species (including the recently described Colombian endemic Callicebus caquetensis). However, there are some notable inconsistences between the obtained phylogeny and current taxonomy. Some previously recognized taxa were not separated by our data (e.g., Callicebus caligatus and Callicebus dubius), while currently unrecognized species diversity was uncovered within C. cupreus in the form of two divergent lineages: one of which exhibited greater phylogenetic similarity to species from the C. moloch Group. Based on the present study, we challenge current taxonomic arrangements for the C. cupreus species Group and call for a thorough taxonomic revision within the genus Callicebus.

Integrative Taxonomic Approach for Describing a New Cryptic Species of Bush Frog (Raorchestes: Anura: Rhacophoridae) from the Western Ghats, India

A new cryptic species of bush frog Raorchestes honnametti sp. nov. is described from the south-eastern part of the Western Ghats, India. This newly described species belongs to the Charius clade and is morphologically similar to other clade members--R. charius and R. griet. Therefore, an integrative taxonomic approach based on molecular and bioacoustic analysis along with morphology was used to delimit the new species. Raorchestes honnametti sp. nov., is currently known only from Biligiri Rangaswamy Temple Tiger Reserve, a part of Biligiri Rangaswamy horst mountain range (a mountain formed due movement of two faults) formed during the Late Quaternary period (1.8-2.58 Ma). Discovery of cryptic species from a highly speciose and well-studied genus Raorchestes hints at the possible existence of several more cryptic species in this genus. We discuss the possible reasons for crypsis and emphasize the need for continued systematic surveys of amphibians across the Western Ghats.

Climate as a driver of tropical insular diversity: comparative phylogeography of two ecologically distinctive frogs in Puerto Rico

The effects of late Quaternary climate on distributions and evolutionary dynamics of insular species are poorly understood in most tropical archipelagoes. We used ecological niche models under past and current climate to derive hypotheses regarding how stable climatic conditions shaped genetic diversity in two ecologically distinctive frogs in Puerto Rico. Whereas the Mountain Coquí, Eleutherodactylus portoricensis, is restricted to montane forest in the Cayey and Luquillo Mountains, the Red-eyed Coquí, E. antillensis, is a habitat generalist distributed across the entire Puerto Rican Bank (Puerto Rico and the Virgin Islands, excluding St. Croix). To test our hypotheses, we conducted phylogeographic and population genetic analyses based on mitochondrial and nuclear loci of each species across their range in Puerto Rico. Patterns of population differentiation in E. portoricensis, but not in E. antillensis, supported our hypotheses. For E. portoricensis, these patterns include: individuals isolated by long-term unsuitable climate in the Río Grande de Loíza Basin in eastern Puerto Rico belong to different genetic clusters; past and current climate strongly predicted genetic differentiation; and Cayey and Luquillo Mountains populations split prior to the last interglacial. For E. antillensis, these patterns include: genetic clusters did not fully correspond to predicted long-term unsuitable climate; and past and current climate weakly predicted patterns of genetic differentiation. Genetic signatures in E. antillensis are consistent with a recent range expansion into western Puerto Rico, possibly resulting from climate change and anthropogenic influences. As predicted, regions with a large area of long-term suitable climate were associated with higher genetic diversity in both species, suggesting larger and more stable populations. Finally, we discussed the implications of our findings for developing evidence-based management decisions for E. portoricensis, a taxon of special concern. Our findings illustrate the role of persistent suitable climatic conditions in promoting the persistence and diversification of tropical island organisms.

A likelihood inference of historical biogeography in the world's most diverse terrestrial vertebrate genus: diversification of direct-developing frogs (Craugastoridae: Pristimantis) across the Neotropics

The geology of the northern Andean region has driven the evolutionary history of Neotropical fauna through the creation of barriers and connections that have resulted in speciation and dispersal events, respectively. One of the most conspicuous groups of anuran fauna in the Andes and surrounding areas is the direct-developing species of the genus Pristimantis. We investigated the molecular phylogenetic placement of 12 species from the montane Andes of Colombia in a broader geographical context with a new genus-level phylogeny in order to identify the role of Andean orogeny over the last 40million years and the effect of elevational differences in diversification of Pristimantis. We examined the biogeographic history of the genus using ancestral range reconstruction by biogeographic regions and elevational ranges. We recognized the middle elevational band (between 1000 and 3000m) in the Northwestern Andes region of Colombia and Ecuador as a focal point for the origin and radiation of Pristimantis species. Additionally, we found several Andean migrations toward new habitats in Central Andes and Merida Andes for some species groups. We suggest that the paleogeological changes in the Northwestern Andes were the main promoter of speciation in Pristimantis, and may have served as a corridor for the dispersion of lowland species.

Extremely low prevalence of Batrachochytrium dendrobatidis in frog populations from neotropical dry forest of Costa Rica supports the existence of a climatic refuge from disease

Population declines and extinctions of numerous species of amphibians, especially stream-breeding frogs, have been linked to the emerging infectious disease chytridiomycosis, caused by the chytrid fungus, Batrachochytrium dendrobatidis. In Central America, most of the 34 species of the Craugastor punctariolus species group have disappeared in recent years in high- and low-elevation rainforests. Distribution models for B. dendrobatidis and the continuous presence of the extirpated stream-dwelling species, Craugastor ranoides, in the driest site of Costa Rica (Santa Elena Peninsula), suggest that environmental conditions might restrict the growth and development of B. dendrobatidis, existing as a refuge from chytridiomycosis-driven extinction. We conducted field surveys to detect and quantify the pathogen using Real-time PCR in samples from 15 species of frogs in two locations of tropical dry forest. In Santa Elena Peninsula, we swabbed 310 frogs, and only one sample of the species, C. ranoides, tested positive for B. dendrobatidis (prevalence <0.1%). In Santa Rosa Station, we swabbed 100 frogs, and nine samples from three species tested positive (prevalence = 9.0%). We failed to detect signs of chytridiomycosis in any of the 410 sampled frogs, and low quantities of genetic equivalents (between 0 and 1073) were obtained from the ten positive samples. The difference in the prevalence between locations might be due not only to the hotter and drier conditions of Santa Elena Peninsula but also to the different compositions of species in both locations. Our results suggest that B. dendrobatidis is at the edge of its distribution in these dry and hot environments of tropical dry forest. This study supports the existence of climatic refuges from chytridiomycosis and highlights the importance of tropical dry forest conservation for amphibians in the face of epidemic disease.

Phylogeography and biogeography of the lower Central American Neotropics: diversification between two continents and between two seas

Lower Central America (LCA) provides a geologically complex and dynamic, richly biodiverse model for studying the recent assembly and diversification of a Neotropical biota. Here, we review the growing literature of LCA phylogeography studies and their contribution to understanding the origins, assembly, and diversification of the LCA biota against the backdrop of regional geologic and climatic history, and previous biogeographical inquiry. Studies to date reveal that phylogeographical signal within taxa of differing distributions reflects a diversity of patterns and processes rivalling the complexities of LCA landscapes themselves. Even so, phylogeography is providing novel insights into regional diversification (e.g. cryptic lineage divergences), and general evolutionary patterns are emerging. Congruent multi-taxon phylogeographic breaks are found across the Nicaraguan depression, Chorotega volcanic front, western and central Panama, and the Darién isthmus, indicating that a potentially shared history of responses to regional-scale (e.g. geological) processes has shaped the genetic diversity of LCA communities. By contrast, other species show unique demographic histories in response to overriding historical events, including no phylogeographic structure at all. These low-structure or incongruent patterns provide some evidence for a role of local, ecological factors (e.g. long-distance dispersal and gene flow in plants and bats) in shaping LCA communities. Temporally, comparative phylogeographical structuring reflects Pliocene-Pleistocene dispersal and vicariance events consistent with the timeline of emergence of the LCA isthmus and its major physiographic features, e.g. cordilleras. We emphasise the need to improve biogeographic inferences in LCA through in-depth comparative phylogeography projects capitalising on the latest statistical phylogeographical methods. While meeting the challenges of reconstructing the biogeographical history of this complex region, phylogeographers should also take up the critical service to society of applying their work to the conservation of its fascinating biodiversity.

Phylogeny and biogeography of the Poecilia sphenops species complex (Actinopterygii, Poeciliidae) in Central America

We inferred the phylogenetic relationships among members of the Poecilia sphenops species complex to resolve the colonization process and radiation of this group in Central America. We analyzed 2550 base pairs (bp) of mitochondrial DNA (mtDNA), including ATP synthase 6 and 8, cytochrome oxidase subunit I and NADH dehydrogenase subunit 2 genes, and 906bp of the nuclear S7 ribosomal protein of 86 ingroup individuals from 61 localities spanning most of its distribution from Mexico to Panama. Our mitochondrial data rendered a well-supported phylogeny for the P. sphenops complex that differed with the nuclear data set topology, which did not recover the monophyly of the P. mexicana mitochondrial lineage. Coalescent-based simulations tests indicated that, although hybridization cannot be completely ruled out, this incongruence is most likely due to incomplete lineage sorting in this group, which also showed the widest geographic distribution. A single colonization event of Central America from South America was estimated to have occurred between the early Paleocene and Oligocene (53-22millionyears ago). Subsequently, two largely differentiated evolutionary lineages diverged around the Early Oligocene-Miocene (38-13million years ago), which are considered two separate species complexes: P. sphenops and P. mexicana, which can also be distinguished by their tricuspid and unicuspid inner jaw teeth, respectively. Ultimately, within lineage diversification occurred mainly during the Miocene (22-5million years ago). All major cladogenetic events predated the final closure of the Isthmus of Panama. The allopatric distribution of lineages together with the long basal internodes suggest that vicariance and long term isolations could be the main evolutionary forces promoting radiation in this group, although dispersal through water barriers might also have occurred. Lastly, our results suggest the need to review the current species distribution and taxonomy of the P. sphenops complex sensu lato.

DNA barcoding applied to ex situ tropical amphibian conservation programme reveals cryptic diversity in captive populations

Amphibians constitute a diverse yet still incompletely characterized clade of vertebrates, in which new species are still being discovered and described at a high rate. Amphibians are also increasingly endangered, due in part to disease-driven threats of extinctions. As an emergency response, conservationists have begun ex situ assurance colonies for priority species. The abundance of cryptic amphibian diversity, however, may cause problems for ex situ conservation. In this study we used a DNA barcoding approach to survey mitochondrial DNA (mtDNA) variation in captive populations of 10 species of Neotropical amphibians maintained in an ex situ assurance programme at El Valle Amphibian Conservation Center (EVACC) in the Republic of Panama. We combined these mtDNA sequences with genetic data from presumably conspecific wild populations sampled from across Panama, and applied genetic distance-based and character-based analyses to identify cryptic lineages. We found that three of ten species harboured substantial cryptic genetic diversity within EVACC, and an additional three species harboured cryptic diversity among wild populations, but not in captivity. Ex situ conservation efforts focused on amphibians are therefore vulnerable to an incomplete taxonomy leading to misidentification among cryptic species. DNA barcoding may therefore provide a simple, standardized protocol to identify cryptic diversity readily applicable to any amphibian community.

Sea level, topography and island diversity: phylogeography of the Puerto Rican Red-eyed Coquí, Eleutherodactylus antillensis

Quaternary climatic oscillations caused changes in sea level that altered the size, number and degree of isolation of islands, particularly in land-bridge archipelagoes. Elucidating the demographic effects of these oscillations increases our understanding of the role of climate change in shaping evolutionary processes in archipelagoes. The Puerto Rican Bank (PRB) (Puerto Rico and the Eastern Islands, which comprise Vieques, Culebra, the Virgin Islands and associated islets) in the eastern Caribbean Sea periodically coalesced during glaciations and fragmented during interglacial periods of the quaternary. To explore population-level consequences of sea level changes, we studied the phylogeography of the frog Eleutherodactylus antillensis across the archipelago. We tested hypotheses encompassing vicariance and dispersal narratives by sequencing mtDNA (c. 552 bp) of 285 individuals from 58 localities, and four nuDNA introns (totalling c. 1633 bp) from 173 of these individuals. We found low support for a hypothesis of divergence of the Eastern Islands populations prior to the start of the penultimate interglacial c. 250 kya, and higher support for a hypothesis of colonization of the Eastern Islands from sources in eastern Puerto Rico during the penultimate and last glacial period, when a land bridge united the PRB. The Río Grande de Loíza Basin in eastern Puerto Rico delineates a phylogeographic break. Haplotypes shared between the PRB and St. Croix (an island c. 105 km south-east of this archipelago) likely represent human-mediated introductions. Our findings illustrate how varying degrees of connectivity and isolation influence the evolution of tropical island organisms.

Comparative phylogeography of direct-developing frogs (Anura: Craugastoridae: Pristimantis) in the southern Andes of Colombia

The Andes of South America hosts perhaps the highest amphibian species diversity in the world, and a sizable component of that diversity is comprised of direct-developing frogs of the genus Pristimantis (Anura: Craugastoridae). In order to better understand the initial stages of species formation in these frogs, this study quantified local-scale spatial genetic structuring in three species of Pristimantis. DNA sequences of two mitochondrial gene fragments (16S and COI) were obtained from P. brevifrons, P. palmeri and P. jubatus at different locations in the Cordillera Occidental. We found high levels of genetic diversity in the three species, with highly structured populations (as measured by F(ST)) in P. brevifrons and P. palmeri while P. jubatus showed panmixia. Large effective population sizes, inferred from the high levels of genetic diversity, were found in the three species and two highly divergent lineages were detected within P. jubatus and P. palmeri. Estimated divergence times among populations within P. brevifrons and P. palmeri coincide with the Pleistocene, perhaps due to similar responses to climatic cycling or recent geological history. Such insights have important implications for linking alpha and beta diversity, suggesting regional scale patterns may be associated with local scale processes in promoting differentiation among populations in the Andes.

Evaluating the role of contracting and expanding rainforest in initiating cycles of speciation across the Isthmus of Panama

Climatic and geological changes across time are presumed to have shaped the rich biodiversity of tropical regions. However, the impact climatic drying and subsequent tropical rainforest contraction had on speciation has been controversial because of inconsistent palaeoecological and genetic data. Despite the strong interest in examining the role of climatic change on speciation in the Neotropics there has been few comparative studies, particularly, those that include non-rainforest taxa. We used bird species that inhabit humid or dry habitats that dispersed across the Panamanian Isthmus to characterize temporal and spatial patterns of speciation across this barrier. Here, we show that these two assemblages of birds exhibit temporally different speciation time patterns that supports multiple cycles of speciation. Evidence for these cycles is further corroborated by the finding that both assemblages consist of 'young' and 'old' species, despite dry habitat species pairs being geographically more distant than pairs of humid habitat species. The matrix of humid and dry habitats in the tropics not only allows for the maintenance of high species richness, but additionally this study suggests that these environments may have promoted speciation. We conclude that differentially expanding and contracting distributions of dry and humid habitats was probably an important contributor to speciation in the tropics.

DNA barcoding amphibians and reptiles

Only a few major research programs are currently targeting COI barcoding of amphibians and reptiles (including chelonians and crocodiles), two major groups of tetrapods. Amphibian and reptile species are typically old, strongly divergent, and contain deep conspecific lineages which might lead to problems in species assignment with incomplete reference databases. As far as known, there is no single pair of COI primers that will guarantee a sufficient rate of success across all amphibian and reptile taxa, or within major subclades of amphibians and reptiles, which means that the PCR amplification strategy needs to be adjusted depending on the specific research question. In general, many more amphibian and reptile taxa have been sequenced for 16S rDNA, which for some purposes may be a suitable complementary marker, at least until a more comprehensive COI reference database becomes available. DNA barcoding has successfully been used to identify amphibian larval stages (tadpoles) in species-rich tropical assemblages. Tissue sampling, DNA extraction, and amplification of COI is straightforward in amphibians and reptiles. Single primer pairs are likely to have a failure rate between 5 and 50% if taxa of a wide taxonomic range are targeted; in such cases the use of primer cocktails or subsequent hierarchical usage of different primer pairs is necessary. If the target group is taxonomically limited, many studies have followed a strategy of designing specific primers which then allow an easy and reliable amplification of all samples.

The Great American Biotic Interchange in frogs: multiple and early colonization of Central America by the South American genus Pristimantis (Anura: Craugastoridae)

The completion of the land bridge between North and South America approximately 3.5-3.1 million years ago (Ma) initiated a tremendous biogeographic event called the Great American Biotic Interchange (GABI), described principally from the mammalian fossil record. The history of biotic interchange between continents for taxonomic groups with poor fossil records, however, is not well understood. Molecular and fossil data suggest that a number of plant and animal lineages crossed the Isthmus of Panama well before 3.5 Ma, leading biologists to speculate about trans-oceanic dispersal mechanisms. Here we present a molecular phylogenetic analysis of the frog genus Pristimantis based on 189 individuals of 137 species, including 71 individuals of 31 species from Panama and Colombia. DNA sequence data were obtained from three mitochondrial (COI, 12S, 16S) and two nuclear (RAG-1 and Tyr) genes, for a total of 4074 base pairs. The resulting phylogenetic hypothesis showed statistically significant conflict with most recognized taxonomic groups within Pristimantis, supporting only the rubicundus Species Series, and the Pristimantis myersi and Pristimantis pardalis Species Groups as monophyletic. Inference of ancestral areas based on a likelihood model of geographic range evolution via dispersal, local extinction, and cladogenesis (DEC) suggested that the colonization of Central America by South American Pristimantis involved at least 11 independent events. Relaxed-clock analyses of divergence times suggested that at least eight of these invasions into Central America took place prior to 4 Ma, mainly in the Miocene. These findings contribute to a growing list of molecular-based biogeographic studies presenting apparent temporal conflicts with the traditional GABI model.

Multilocus analysis of the evolutionary dynamics of rainforest bird populations in Southeast Asia

Sundaland has a dynamic geographic history because its landmasses were periodically interconnected when sea levels fell during glacial periods. Superimposed on this geographic dynamism were environmental changes related to climatic oscillations. To investigate how tropical taxa responded to such changes, we studied the divergence and demographic history of two co-distributed rainforest passerine species, Arachnothera longirostra and Malacocincla malaccensis. We sampled birds primarily from Borneo and the Malay Peninsula, which straddle the now-submerged Sunda shelf, and analysed multilocus DNA data with a variety of coalescent and gene genealogy methods. Cross-shelf divergence in both species occurred well before the last glacial maximum, i.e., before the most recent land connection. However, post-divergence gene flow occurred, and it was more pronounced in A. longirostra (a highly vagile nectarivore/insectivore) than in M. malaccensis (an understory insectivore). Despite current habitat continuity on Borneo, the population of M. malaccensis in northeastern Borneo is substantially divergent from that on the rest of the island. The NE population experienced dramatic demographic fluctuations, probably because of competition with the other population, which expanded from western Borneo after the mid-Pleistocene. In contrast, the Borneo population of A. longirostra has little structure and appears to have experienced demographic expansion 16 kya, long after it had diverged from the Malay Peninsula population (630-690 kya). Malay Peninsula populations of both species have remained relatively stable. Overall, the most recent glacial period was not the chief determinant of the evolutionary dynamics of our study species, and in this respect, they are different from temperate species.

Correlates of species richness in the largest Neotropical amphibian radiation

Although tropical environments are often considered biodiversity hotspots, it is precisely in such environments where least is known about the factors that drive species richness. Here, we use phylogenetic comparative analyses to study correlates of species richness for the largest Neotropical amphibian radiation: New World direct-developing frogs. Clade-age and species richness were nonsignificantly, negatively correlated, suggesting that clade age alone does not explain among-clade variation in species richness. A combination of ecological and morphological traits explained 65% of the variance in species richness. A more vascularized ventral skin, the ability to colonize high-altitude ranges, encompassing a large variety of vegetation types, correlated significantly with species richness, whereas larger body size was marginally correlated with species richness. Hence, whereas high-altitude ranges play a role in shaping clade diversity in the Neotropics, intrinsic factors, such as skin structures and possibly body size, might ultimately determine which clades are more speciose than others.

Ecological and genetic divergence between two lineages of middle American túngara frogs Physalaemus (= Engystomops) pustulosus

BACKGROUND

Uncovering how populations of a species differ genetically and ecologically is important for understanding evolutionary processes. Here we combine population genetic methods (microsatellites) with phylogenetic information (mtDNA) to define genetic population clusters of the wide-spread Neotropical túngara frog (Physalaemus pustulosus). We measure gene flow and migration within and between population clusters and compare genetic diversity between population clusters. By applying ecological niche modeling we determine whether the two most divergent genetic groups of the túngara frog (1) inhabit different habitats, and (2) are separated geographically by unsuitable habitat across a gap in the distribution.

RESULTS

Most population structure is captured by dividing all sample localities into two allopatric genetic lineages. The Northern genetic lineage (NW Costa Rica) is genetically homogenous while the Southern lineage (SW Costa Rica and Panama) is sub-divided into three population clusters by both microsatellite and mtDNA analyses. Gene flow is higher within the Northern lineage than within the Southern lineage, perhaps due to increased landscape heterogeneity in the South. Niche modeling reveals differences in suitable habitat between the Northern and Southern lineages: the Northern lineage inhabits dry/pine-oak forests, while the Southern lineage is confined to tropical moist forests. Both lineages seem to have had little movement across the distribution gap, which persisted during the last glacial maximum. The lack of movement was more pronounced for the Southern lineage than for the Northern lineage.

CONCLUSIONS

This study confirms the finding of previous studies that túngara frogs diverged into two allopatric genetic lineages north and south of the gap in the distribution in central Costa Rica several million years ago. The allopatric distribution is attributed to unsuitable habitat and probably other unknown ecological factors present across the distribution gap. Niche conservatism possibly contributes to preventing movements across the gap and gene flow between both groups. Genetic and ecological data indicate that there is the potential for ecological divergence in allopatry between lineages. In this context we discuss whether the Northern and Southern lineages should be recognized as separate species, and we conclude that further studies of pre- and post-zygotic isolation are needed for a final assessment. Identified population clusters should motivate future behavioral and ecological research regarding within-species biodiversity and speciation mechanisms.

Genetic diversification, vicariance, and selection in a polytypic frog

Spatial patterns of heritable phenotypic diversity reflect the relative roles of gene flow and selection in determining geographic variation within a species. We quantified color differentiation and genetic divergence among 20 populations of the red-eyed tree frog (Agalychnis callidryas) in lower Central America. Phylogenetic analyses revealed 5 well-supported mitochondrial DNA clades, and we infer from our phylogeny that geographic barriers have played a large role in structuring populations. Two phenotypic characters varied independently among isolated population groups: Flank coloration distinguished Caribbean from Pacific individuals, whereas leg coloration exhibited a more complex geographic pattern. We detected 3 generalized spatial patterns of genetic and phenotypic diversity: 1) phenotypic differentiation in the presence of historical connectivity, 2) phenotypic uniformity across genetically differentiated regions, and 3) codistribution of genetic and phenotypic characters. These patterns indicate that phenotypic diversification is highly regionalized and can result from spatial variation in localized adaptations, geographic isolation, genetic drift, and/or evolutionary stasis. Although the mode of selection underlying color variation was not the focal objective of this study, we discuss the possible roles of natural and sexual selection in mediating population differentiation. Our study underscores the fact that selection gradients vary across relatively small spatial scales, even in species that occupy relatively homogeneous environments.

Rapid color evolution in an aposematic species: a phylogenetic analysis of color variation in the strikingly polymorphic strawberry poison-dart frog

Aposematism is one of the great mysteries of evolutionary biology. The evolution of aposematic coloration is poorly understood, but even less understood is the evolution of polymorphism in aposematic signals. Here, we use a phylogeographic approach to investigate the evolution of color polymorphism in Dendrobates pumilio, a well-known poison-dart frog (family Dendrobatidae), which displays perhaps the most striking color variation of any aposematic species. With over a dozen color morphs, ranging from bright red to dull green, D. pumilio provides an ideal opportunity to examine the evolution of color polymorphism and evolutionary shifts to cryptic coloration in an otherwise aposematic species. We constructed a phylogenetic tree for all D. pumilio color morphs from 3051bp of mtDNA sequence data, reconstructed ancestral states using parsimony and Bayesian methods, and tested the recovered tree against constraint trees using parametric bootstrapping to determine the number of changes to each color type. We find strong evidence for nearly maximal numbers of changes in all color traits, including five independent shifts to dull dorsal coloration. Our results indicate that shifts in coloration in aposematic species may occur more regularly than predicted and that convergence in coloration may indicate that similar forces are repeatedly driving these shifts.

Phylogeography of the Pygmy Rain Frog (Pristimantis ridens) across the lowland wet forests of isthmian Central America

We used a phylogeographic approach to elucidate the evolutionary history of a lineage of frogs, known as Pristimantis (formerly Eleutherodactylus) ridens (Anura: Brachycephalidae), restricted to the wet forests occurring along the Caribbean versant of isthmian Central America as well as the disjunct wet forest on the Pacific slope of Costa Rica. We placed our phylogeographic study of P. ridens within a larger molecular phylogenetic analysis of Central American Pristimantis. All phylogenetic inferences were based on a 1455 base pair fragment of mitochondrial DNA, containing the complete ND2 gene and five flanking tRNA genes. Our reconstruction of the intraspecific phylogeny of P. ridens yielded a basal trichotomy dating to an estimated 12+ million years ago (Ma), consisting of central Panama, western Panama, and Costa Rica plus Honduras. Thus, the presence of P. ridens appears to predate the completion of the Isthmus 3.1Ma. Using a parametric bootstrap (SOWH) test, we evaluated four a priori zoogeographic hypotheses for the origin and spread of P. ridens. This analysis suggested that the P. ridens populations on the Caribbean versant of Costa Rica were established by Pacific versant ancestors only recently, in contrast to the very old lineages found in Panama. Our results support a model of Miocene colonization, long-term geographic stasis, followed by rapid dispersal across the Caribbean lowlands during the Pliocene or Pleistocene.