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

Systematics, distribution patterns and historical biogeography of the Central America wandering spider genus Kiekie Polotow & Brescovit, 2018 (Araneae: Ctenidae)

Kiekie Polotow & Brescovit, 2018 is a Neotropical genus of Ctenidae, with most of its species occuring in Central America. In this study, we review the systematics of Kiekie and describe five new species and the unknown females of K. barrocolorado Polotow & Brescovit, 2018 and K. garifuna Polotow & Brescovit, 2018, and the unknown male of K. verbena Polotow & Brescovit, 2018. In addition, we described the female of K. montanense which was wrongly assigned as K. griswoldi Polotow & Brescovit, 2018 (both species are sympatric). We provided a modified diagnosis for previously described species based on the morphology of the newly discovered species and in situ photographs of living specimens. We inferred a molecular phylogeny using four nuclear (histone H3, 28S rRNA, 18S rRNA and ITS-2) and three mitochondrial genes (cytochrome c oxidase subunit I or COI, 12S rRNA and 16S rRNA) to test the monophyly of the genus and the evolutionary relationships of its species. Lastly, we reconstruct the historical biogeography and map diversity and endemism distributional patterns of the different species. This study increased the number of known species of Kiekie from 13 to 18, and we describe a new genus, Eldivo which is sister lineage of Kiekie. Most of the diversity and endemism of the genus Kiekie is located in the montane ecosystems of Costa Rica followed by the lowland rainforest of the Pacific side (Limon Basin). Kiekie originated in the North America Tropical region, this genus started diversifying in the Late Miocene and spread to Lower Central America and South America. In that region, Kiekie colonized independently several times the montane ecosystems corresponding to periods of uplifting of Talamanca and Central Cordilleras.

Impact of persistent barrier to gene flow and catastrophic events on red algae evolutionary history along the Chilean coast

Historical vicariance events, linked to the existence of stable physical barriers to gene flow, generate concordant genetic breaks in co-distributed species while stochastic processes (e.g., costal uplift) could cause species-specific genetic breaks as a result of local strong demographic bottlenecks or extinction. In Chile, previous studies show that the area of the 30°S-33°S could correspond to a stable barrier to gene flow that have affected the genetic structure of various algae and marine invertebrates. Here we sequenced two organellar genes (COI and rbcL) in four taxonomically accepted co-distributed red seaweeds species characterized by a low dispersal potential: Mazzaella laminarioides, M. membranacea, Asterfilopsis disciplinalis, and Ahnfeltiopsis vermicularis. Our results revealed the existence of ten strongly differentiated linages in the taxa studied. Strong genetic breaks, concordant in both space and time (divergence estimated to have occurred some 2.9-12.4 million years ago), were observed between taxa distributed across the 33°S. Conversely, in the Central/South part of the Chilean coast, the localization of the genetic breaks/sub-structure observed varied widely (36°S, 38°S, 39°S, and 40°S). These results suggest that a major historical vicariance event has modeled the genetic structure of several Chilean marine organisms in the north of the Chilean coast during the mid-Miocene, while more recent stochastic events and genetic drift could be the driving forces of genetic divergence/structuration in the central-southern part of the coast.

Genetic analyses reveal panmixia in Indian waters and population subdivision across Indian Ocean and Indo-Malay Archipelago for Decapterus russelli

The Indian Scad, Decapterus russelli is an important pelagic carangid widely distributed throughout the Indian Ocean and the Indo-West Pacific. Despite being widely distributed in the Indian Ocean, the information regarding genetic structuring and diversity has been lacking compared to its Indo Malay Archipelago counterparts. The present study was conducted to investigate the genetic stock structure of D. russelli based on mitochondrial (Cyt b) and nuclear (DrAldoB1) markers along Indian waters. The results indicated the presence of a single panmictic stock across the Indian Ocean region. High haplotype diversity associated with low nucleotide diversity suggested a population bottleneck followed by rapid population growth. Phylogenetic analysis revealed the absence of geographical clustering of lineages with the most common haplotype distributed globally. The pelagic life style, migratory capabilities, and larval dispersal may be the contributing factors to the observed spatial homogeneity of D. russelli. However, significant genetic differentiation was observed between the populations from Indian Ocean and Indo-Malay Archipelago. Hierarchical molecular variance analysis (AMOVA), pairwise FST comparisons and SAMOVA showed existence of two distinct genetic stocks of D. russelli in the Indian Ocean and IMA. The observed interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the Indo-Malay Archipelago and the Indian Ocean regions suggest the influence of historic isolation, ocean surface currents and biotic and abiotic features of the ocean. Also, there was a significant relationship between genetic distance and geographical distance between population pairs in a manner consistent with isolation-by-distance. These resulted in the evolution of a phylogeographic break for this species between these regions. The findings of these results suggest that D. russelli from the Indian Ocean shall be managed in its entire area of distribution as a single stock. Further, the Indian Ocean and Indo-Malayan stocks can be managed separately.

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.

Late Pleistocene Altitudinal Segregation and Demography Define Future Climate Change Distribution of the Peromyscus mexicanus Species Group: Conservation Implications

Mountains harbor a significant number of the World's biodiversity, both on tropical and temperate regions. Notably, one crucial gap in conservation is the consideration of historical and contemporary patterns influencing differential distribution in small mammal mountain species and how climate change will affect their distribution and survival. The mice Peromyscus mexicanus species group is distributed across mountains in Guatemala-Chiapas and Central America, which experienced significant effects of glacial and interglacial cycles. We determined phylogeographic and demographic patterns of lowlands and highlands mountain lineages, revealing that the radiation of modern P. mexicanus lineages occurred during the Pleistocene (ca. 2.6 mya) along Nuclear Central America. In concert with climatic cycles and the distribution of habitats, lowland and highland lineages showed recent population size increase and decrease, respectively. We also estimated the current and future distribution ranges for six lineages, finding marked area size increase for two lineages for which vegetation type and distribution would facilitate migrating towards higher elevations. Contrastingly, three lineages showed range size decrease; their ecological requirements make them highly susceptible to future habitat loss. Our findings are clear evidence of the negative impacts of future climate change, while our ability to manage and conserve these vulnerable ecosystems and mountain species is contingent on our understanding of the implications of climate change on the distribution, ecology, and genetics of wildlife populations.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

Phylogenomics of trans-Andean tetras of the genus Hyphessobrycon Durbin 1908 (Stethaprioninae: Characidae) and colonization patterns of Middle America

Hyphessobrycon is one of the most species rich and widely distributed genera in the family Characidae, with more than 160 species ranging from Veracruz, Mexico to Mar Chiquita Lagoon in Buenos Aires, Argentina. The majority of Hyphessobrycon diversity shows a cis-Andean distribution; only nine species are trans-Andean including H. compressus (Meek 1908). It is well established that Hyphessobrycon is not monophyletic but it has been suggested that natural groups can be identified within the larger Hyphessobrycon species group. In this study, we tested the monophyly of trans-Andean species of Hyphessobrycon and investigated the placement of H. compressus. We inferred the first phylogenomic hypothesis of trans-Andean Hyphessobrycon that includes nearly complete taxonomic sampling (eight of nine valid species) using ultraconserved elements (UCEs). We analyzed 75% (1682 UCEs), 90% (1258 UCEs), and 95% (838 UCEs) complete data matrices, and inferred phylogenomic hypotheses under concatenation and coalescent approaches. In all cases, we recovered the monophyly of trans-Andean Hyphessobrycon inclusive of H. compressus, strong support for three species groups, and evidence of cryptic diversity within the widespread H. compressus and H. condotensis. We used our phylogenomic hypothesis to investigate the biogeographic history of Hyphessobrycon in Middle America. Our ancestral range estimation analysis suggests a single event of cis- to trans-Andean colonization followed by stepwise colonization from the Pacific slope of northwestern South America (Chocó block) to northern Middle America (Maya block). Our work supports the recognition of the trans-Andean species as Hyphessobrycon sensu stricto and provides an evolutionary template to examine morphological characters that will allow us to better understand the diversity of Hyphessobrycon in Middle America.

Filling the knowledge gap of Middle American freshwater fish parasite biodiversity: metazoan parasite fauna of Nicaragua

The heterogeneous landscape of Nicaragua harbours a large diversity of freshwater fishes. The great Nicaraguan lakes, Managua and Nicaragua, and several adjacent crater lakes harbour numerous endemic fish species. However, information about their parasite fauna is still fragmentary. Here, we surveyed the great Nicaraguan lakes and four crater lakes and provide data for 17 metazoan parasite taxa infecting seven fish host species. We also gathered all the published records from the literature on the parasites reported from Nicaraguan freshwater fishes, as well as those for Costa Rica and Panama to discuss the region of Lower Central America as a whole. With this information we built a parasite-host and a host-parasite checklist. With data from near 50% of the native and endemic freshwater fishes in Nicaragua, the parasite fauna comprises 101 taxa in 51 fish species allocated in 11 families. Cichlids are the most diverse group of fishes in this region and have been the most extensively surveyed for their metazoan parasites. Helminths are the best-represented groups of metazoan parasites, with 42 trematodes, five cestodes, 24 monogeneans, two acanthocephalans, 20 nematodes and one hirudinean. Additionally, freshwater fishes are parasitized by copepods, branchiurans and oribatid mites. Even though the inventory is not yet complete, the patterns of diversity uncovered revealed promising information about the origin, biogeography and evolutionary history of the Nicaraguan freshwater fish parasite fauna. More studies are necessary to complete our knowledge about the diversity, host association and distribution of metazoan parasites in Nicaragua and other Central American countries.

Testing spatial and environmental factors to explain body shape variation in the widespread Central American Blackbelt cichlid Vieja maculicauda (Teleostei: Cichlidae)

Abstract Body-shape variability may respond to the interaction between history and environment, particularly in species whose range includes known biogeographical barriers. Central America has traditionally been regarded to have functioned as a complete land-bridge connecting two continents since the Plio-Pleistocene and as an incomplete one since much earlier. This history has helped shape species diversification and distributions. Vieja maculicauda is a widely distributed cichlid, found throughout most of the Atlantic slope of Central America, across three of the four geological blocks that make up this region and whose borders are known dispersal barriers. We asked whether there is an effect of geologic blocks on body shape of V. maculicauda as determined by geometric morphometric assessment of 151 individuals from across its range. We asked how much variance in body shape is explained by environmental variables. We used a CVA and a MANCOVA to determine whether a relationship between body shape and geologic blocks exists. For testing body shape correlation against geographic and environmental variables, we used a Mantel test. Our results suggest that body shape in V. maculicauda can be segregated by geologic blocks, while the correlation results showed no strong correlation between our sets of variables.

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon-specific attributes.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Genetic Introgression and Morphological Variation in Naked-Back Bats (Chiroptera: Mormoopidae: Pteronotus Species) along Their Contact Zone in Central America

Two sibling bare-backed bat species (Pteronotus fulvus and P. gymnonotus) have been traditionally differentiated by their size. However, intermediate specimens between the two species have been found in sympatric populations along southern Mexico and it has been suggested that they may be the outcome of a hybridization process between the two species. We used one mitochondrial (COI), three nuclear markers (PRKCL, STAT5A and RAG2) and 13 microsatellites to explore the evolutionary relationships between these two species and elucidate whether the intermediate morphotypes correspond to hybrid individuals. These markers have been analyzed in sympatric and allopatric populations of the two species plus the closely related species Pteronotus davyi. We confirmed the species-level differentiation of the three lineages (P. fulvus, P. davyi and P. gymnonotus), but the phylogenetic hypotheses suggested by the nuclear and mitochondrial markers were discordant. We confirm that the discordance between markers is due to genetic introgression through the mitochondrial capture of P. fulvus in P. gymnonotus populations. Such introgression was found in all P. gymnonotus specimens across its sympatric distribution range (Mexico to Costa Rica) and is related to expansion/retraction species distribution pulses associated with changes in forest distribution during the Quaternary climate cycles. Microsatellite analyses showed contemporary genetic contact between the two sympatric species and 3.0% of the samples studied were identified as hybrids. In conclusion, we found a historical and asymmetric genetic introgression (through mitochondrial capture) of P. fulvus into P. gymnonotus in Mexico and Central America and a limited contemporary gene exchange between the two species. However, no relationship was found between hybridization and the intermediate-sized specimens from southern Mexico, which might likely result from a clinal variation with latitude. These results confirm the need for caution when using forearm size to identify these species in the field and when differentiating them in the laboratory based on mitochondrial DNA alone.

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.

A new vine snake (Reptilia, Colubridae, Oxybelis) from Peru and redescription of O. acuminatus

The Brown Vine Snake, Oxybelis aeneus, was until recently considered a single species, distributed from southern Arizona through the Neotropics into southeastern Brazil. However, newly conducted research restructured the species with a substantial taxonomic revision, recognizing five additional taxa (i.e. O. koehleri, O. microphthalmus, O. potosiensis, O. rutherfordi, O. vittatus) in this species complex. This revision focused on populations in North America, Central America, and northern South America while neglecting the southern portion of its distribution. Here, we examine the taxonomic history of the complex and use it along with specimen data to resurrect O. acuminatus from southeastern Brazil. Finally, we describe a new species from the Peruvian Amazon based on morphological characters. This work increases the species diversity of the O. aeneus complex to eight, and we expect further increases in biodiversity discoveries with continued exploration of the New World vine snakes.

Towards a synthesis of the Caribbean biogeography of terrestrial arthropods

BACKGROUND

The immense geologic and ecological complexity of the Caribbean has created a natural laboratory for interpreting when and how organisms disperse through time and space. However, competing hypotheses compounded with this complexity have resulted in a lack of unifying principles of biogeography for the region. Though new data concerning the timing of geologic events and dispersal events are emerging, powerful new analytical tools now allow for explicit hypothesis testing. Arthropods, with varying dispersal ability and high levels of endemism in the Caribbean, are an important, albeit understudied, biogeographic model system. Herein, we include a comprehensive analysis of every publicly available genetic dataset (at the time of writing) of terrestrial Caribbean arthropod groups using a statistically robust pipeline to explicitly test the current extent of biogeographic hypotheses for the region.

RESULTS

Our findings indicate several important biogeographic generalizations for the region: the South American continent is the predominant origin of Caribbean arthropod fauna; GAARlandia played a role for some taxa in aiding dispersal from South America to the Greater Antilles; founder event dispersal explains the majority of dispersal events by terrestrial arthropods, and distance between landmasses is important for dispersal; most dispersal events occurred via island hopping; there is evidence of 'reverse' dispersal from islands to the mainland; dispersal across the present-day Isthmus of Panama generally occurred prior to 3 mya; the Greater Antilles harbor more lineage diversity than the Lesser Antilles, and the larger Greater Antilles typically have greater lineage diversity than the smaller islands; basal Caribbean taxa are primarily distributed in the Greater Antilles, the basal-most being from Cuba, and derived taxa are mostly distributed in the Lesser Antilles; Jamaican taxa are usually endemic and monophyletic.

CONCLUSIONS

Given the diversity and deep history of terrestrial arthropods, incongruence of biogeographic patterns is expected, but focusing on both similarities and differences among divergent taxa with disparate life histories emphasizes the importance of particular qualities responsible for resulting diversification patterns. Furthermore, this study provides an analytical toolkit that can be used to guide researchers interested in answering questions pertaining to Caribbean biogeography using explicit hypothesis testing.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Phylogeography of the Neotropical epiphytic orchid, Brassavola nodosa: evidence for a secondary contact zone in northwestern Costa Rica

Spatial patterns of genetic variation can reveal otherwise cryptic evolutionary and landscape processes. In northwestern Costa Rica, an approximately concordant genetic discontinuity occurs among populations of several plant species. We conducted phylogeographic analyses of an epiphytic orchid, Brassavola nodosa, to test for genetic discontinuity and to explore its underlying causes. We genotyped 18 populations with 19 nuclear loci and two non-coding chloroplast sequence regions. We estimated genetic diversity and structure, relative importance of pollen and seed dispersal, and divergence time to understand how genetic diversity was spatially partitioned. Nuclear genetic diversity was high with little differentiation among populations (GSTn = 0.065). In contrast, chloroplast haplotypes were highly structured (GSTc = 0.570) and reveal a discontinuity between northwestern and southeastern populations within Costa Rica. Haplotype differences suggest two formerly isolated lineages that diverged ~10,000-100,000 YBP. Haplotype mixing and greater genetic diversity occur in an intermediate transition zone. Patterns of nuclear and chloroplast data were consistent. Different levels of genetic differentiation for the two genomes reflect the relative effectiveness of biotic versus abiotic dispersers of pollen and seeds, respectively. Isolation of the two lineages likely resulted from the complex environmental and geophysical history of the region. Our results suggest a recent cryptic seed dispersal barrier and/or zone of secondary contact. We hypothesize that powerful northeasterly trade winds hinder movement of wind-borne seeds between the two regions, while the multi-directional dispersal of pollen by strong-flying sphinx moths resulted in lower differentiation of nuclear loci.

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.

DNA barcode and minibarcode identification of freshwater fishes from Cerrado headwater streams in Central Brazil

The extraordinary species diversity of the Neotropical freshwater fish fauna is world renown. Yet, despite rich species diversity, taxonomic and genetic resources for its Cerrado ichthyofauna remain poorly developed. We provide a reference library of 149 DNA barcodes for 39 species/lineages of Cerrado headwater stream fishes from the Brazilian Distrito Federal and nearby areas and test the utility of distance-based criteria, tree-based criteria and minibarcodes for specimen identification. Mean Kimura 2-parameter genetic distances within species to orders ranged 1·8-12·1%. However, mean intraspecific v. congeneric-interspecific distances (0·9-1·3%) overlapped extensively and distance-based barcoding failed to achieve correct identifications due to c. 4-12·1% error rates and 19·5% ambiguous identifications related to the presence of singletons. Overlap was reduced and best-match success rates improved drastically to 83·5% when Characidium barcodes representing potential misidentifications or undescribed species were removed. Tree-based monophyly criteria generally performed similarly to distance methods, correctly differentiating up to c. 85% of species/lineages despite neighbour-joining and Bayesian tree errors (random lineage-branching events, long-branch attraction). Five clusters (Ancistrus aguaboensis, Characidium spp., Eigenmannia trilineata, Hasemania hanseni and Hypostomus sp. 2) exhibited deep intraspecific divergences or para-/polyphyly and multiple Barcode Index Number assignments indicative of putative candidate species needing taxonomic re-examination. Sliding-window analyses also indicated that a 200 bp minibarcode region performed just as well at specimen identification as the entire barcode gene. Future DNA barcoding studies of Distrito Federal-Cerrado freshwater fishes will benefit from increased sampling coverage, as well as consideration of minibarcode targets for degraded samples and next-generation sequencing.

Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level

Batrachochytrium dendrobatidis (Bd) has been associated with the severe declines and extinctions of amphibians in Costa Rica that primarily occurred during the 1980s and 1990s. However, the current impact of Bd infection on amphibian species in Costa Rica is unknown. We aimed to update the list of amphibian species in Costa Rica and evaluate the prevalence and infection intensity of Bd infection across the country to aid in the development of effective conservation strategies for amphibians. We reviewed taxonomic lists and included new species descriptions and records for a total of 215 amphibian species in Costa Rica. We also sampled for Bd at nine localities from 2015–2018 and combined these data with additional Bd occurrence data from multiple studies conducted in amphibian communities across Costa Rica from 2005–2018. With this combined dataset, we found that Bd was common (overall infection rate of 23%) across regions and elevations, but infection intensity was below theoretical thresholds associated with mortality. Bd was also more prevalent in Caribbean lowlands and in terrestrial amphibians with an aquatic larval stage; meanwhile, infection load was the highest in direct-developing species (forest and stream-dwellers). Our findings can be used to prioritize regions and taxonomic groups for conservation strategies.

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.

Phylogeography indicates incomplete genetic divergence among phenotypically differentiated montane forest populations of Atlapetesalbinucha (Aves, Passerellidae)

The White-naped Brushfinch (Atlapetesalbinucha) comprises up to eight allopatric subspecies mainly identified by the color of the underparts (gray vs. yellow belly). Yellow and gray bellied forms were long considered two different species (A.albinucha and A.gutturalis), but they are presently considered as one polytypic species. Previous studies in the genus Atlapetes have shown that the phylogeny, based on molecular data, is not congruent with characters such as coloration, ecology, or distributional patterns. The phylogeography of A.albinucha was analyzed using two mitochondrial DNA regions from samples including 24 different localities throughout montane areas from eastern Mexico to Colombia. Phylogeographic analyses using Bayesian inference, maximum likelihood and haplotype network revealed incomplete geographic structure. The genetic diversity pattern is congruent with a recent process of expansion, which is also supported by Ecological Niche Models (ENM) constructed for the species and projected into three past scenarios. Overall, the results revealed an incomplete genetic divergence among populations of A.albinucha in spite of the species’ ample range, which contrasts with previous results of phylogeographic patterns in other Neotropical montane forest bird species, suggesting idiosyncratic evolutionary histories for different taxa throughout the region.

Phylogeographic and diversification patterns of the white-nosed coati (Nasua narica): Evidence for south-to-north colonization of North America

White-nosed coatis (Nasua narica) are widely distributed throughout North, Central, and South America, but the patterns of temporal and spatial diversification that have contributed to this distribution are unknown. In addition, the biogeographic history of procyonid species in the Americas remains contentious. Using sequences from three mitochondrial loci (Cytochrome b, NAHD5 and 16S rRNA; 2201 bp) and genotypes from 11 microsatellite loci, we analyzed genetic diversity to determine phylogeographic patterns, genetic structure, divergence times, and gene flow among Nasua narica populations throughout the majority of the species' range. We also estimated the ancestral geographic range of N. narica and other procyonid species. We found a high degree of genetic structure and divergence among populations that conform to five evolutionarily significant units. The most southerly distributed population (Panama) branched off much earlier (∼3.8 million years ago) than the northern populations (<1.2 million years ago). Estimated gene flow among populations was low and mostly northwards and westwards. The phylogeographic patterns within N. narica are associated with geographic barriers and habitat shifts likely caused by Pliocene-Pleistocene climate oscillations. Significantly, our findings suggest the dispersal of N. narica was south-to-north beginning in the Pliocene, not in the opposite direction during the Pleistocene as suggested by the fossil record, and that the most recent common ancestor for coati species was most likely distributed in South or Central America six million years ago. Our study implies the possibility that the diversification of Nasua species, and other extant procyonid lineages, may have occurred in South America.

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.

Is local trait variation related to total range size of tropical trees?

The reasons why the range size of closely related species often varies significantly have intrigued scientists for many years. Among other hypotheses, species with high trait variation were suggested to occupy more diverse environments, have more continuity in their distributions, and consequently have larger range sizes. Here, using 34 tree species of lowlands tropical rainforest in southern Costa Rica, we explored whether inherent trait variability expressed at the local scale in functional traits is related to the species’ total geographical range size. We formed 17 congeneric pairs of one narrow endemic and one widespread species, sampled 335 individuals and measured eight functional traits: leaf area, leaf thickness, leaf dry matter content, specific leaf area, leaf nitrogen content, leaf phosphorus content, leaf nitrogen to phosphorus ratio, and wood specific gravity. We tested whether there are significant differences in the locally expressed variation of individual traits or in multidimensional trait variance between the species in congeneric pairs and whether species’ range size could hence be predicted from local trait variability. However, we could not find such differences between widely distributed and narrow range species. We discuss the possible reasons for these findings including the fact that higher trait variability of widespread species may result from successive local adaptations during range expansion and may hence often be an effect rather than the cause of larger ranges.

Accelerated Evolution and Functional Divergence of the Dim Light Visual Pigment Accompanies Cichlid Colonization of Central America

Cichlids encompass one of the most diverse groups of fishes in South and Central America, and show extensive variation in life history, morphology, and colouration. While studies of visual system evolution in cichlids have focussed largely on the African rift lake species flocks, Neotropical cichlids offer a unique opportunity to investigate visual system evolution at broader temporal and geographic scales. South American cichlid colonization of Central America has likely promoted accelerated rates of morphological evolution in Central American lineages as they encountered reduced competition, renewed ecological opportunity, and novel aquatic habitats. To investigate whether such transitions have influenced molecular evolution of vision in Central American cichlids, we sequenced the dim-light rhodopsin gene in 101 Neotropical cichlid species, spanning the diversity of the clade. We find strong evidence for increased rates of evolution in Central American cichlid rhodopsin relative to South American lineages, and identify several sites under positive selection in rhodopsin that likely contribute to adaptation to different photic environments. We expressed a Neotropical cichlid rhodopsin protein invitro for the first time, and found that while its spectral tuning properties were characteristic of typical vertebrate rhodopsin pigments, the rate of decay of its active signalling form was much slower, consistent with dim light adaptation in other vertebrate rhodopsins. Using site-directed mutagenesis combined with spectroscopic assays, we found that a key amino acid substitution present in some Central American cichlids accelerates the rate of decay of active rhodopsin, which may mediate adaptation to clear water habitats.

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.

A new species of Trachymyrmex (Hymenoptera, Formicidae) fungus-growing ant from the Sierra Madre Oriental of northeastern Mexico

Here we describe a new species of Trachymyrmex, T. pakawasp. n., from the Gran Sierra Plegada range of the Sierra Madre Oriental, in the states of Coahuila and Nuevo Leon, northeastern Mexico. Trachymyrmex pakawa is a large-sized species compared to other North American Trachymyrmex. Its geographic distribution includes the piedmont of the Gran Sierra Plegada at La Estanzuela, Monterrey, as well as peripheral mountains segregated from the Sierra Madre Oriental (Cerro de las Mitras, Sierra de Zapalinamé, Cañon de San Lorenzo, Cerro de las Letras). The preferred habitats of T. pakawa include oak-pine forest at La Estanzuela, xeric oak forest at Zapalinamé and mesic Chihuahuan desert scrub with sotol (Dasylirion) at other sites. All localities are on slopes, on very rocky, shallow lithosols overlaying large boulders. This species nests under and between large boulders and rocks. It has not been observed on alluvial or better developed, deeper soils, and it is absent from sites with human activity (urban, disturbed, and landscaped areas). It is closely related to and morphologically similar to Trachymyrmex smithi. The known distribution ranges of T. pakawa and T. smithi almost overlap in Saltillo, Coahuila state. The main character that distinguishes the new species from T. smithi is longer antennal scapes in T. pakawa; also, different nesting habits (rocky slopes vs. alluvial sites or deep sand in T. smithi), and geographic distribution. Phylogenetic analysis of DNA sequences from the mitochondrial marker cytochrome c oxidase subunit I (COI) and the first intron of the F1 copy of the nuclear protein-coding gene Elongation Factor 1- α (EF1-α-F1) confirm a sister-species relationship between T. pakawa and T. smithi. Bayesian coalescent analyses indicate a divergence time of about 8.00 million years before present (95% confidence interval: 4.8-11.5 mya) between T. pakawa and T. smithi. The divergence of the lineages of T. pakawa and T. smithi could have been driven by the Pliocene-Holocene desertification of southwestern North America. This process resulted in isolated mesic refugia and forests in the Madrean ranges and piedmonts of northeastern Mexico (the current habitat of T. pakawa) while T. smithi adapted to the deeper, often sandy soils on the drier desert plains of Coahuila and Chihuahua states in Mexico, and New Mexico and Texas in the USA. Within the Nearctic species of the Trachymyrmex septentrionalis species group, T. pakawa is the species that is closest (by geographical distribution) to Neotropical species of Trachymyrmex like T. saussurei.

AFLP diversity and spatial structure of Calycophyllum candidissimum (Rubiaceae), a dominant tree species of Nicaragua's critically endangered seasonally dry forest

The Central American seasonally dry tropical (SDT) forest biome is one of the worlds' most endangered ecosystems, yet little is known about the genetic consequences of its recent fragmentation. A prominent constituent of this biome is Calycophyllum candidissimum, an insect-pollinated and wind-dispersed canopy tree of high socio-economic importance, particularly in Nicaragua. Here, we surveyed amplified fragment length polymorphisms across 13 populations of this species in Nicaragua to elucidate the relative roles of contemporary vs historical factors in shaping its genetic variation. Genetic diversity was low in all investigated populations (mean HE=0.125), and negatively correlated with latitude. Overall population differentiation was moderate (ΦST=0.109, P<0.001), and Bayesian analysis of population structure revealed two major latitudinal clusters (I: 'Pacific North'+'Central Highland'; II: 'Pacific South'), along with a genetic cline between I and II. Population-based cluster analyses indicated a strong pattern of 'isolation by distance' as confirmed by Mantel's test. Our results suggest that (1) the low genetic diversity of these populations reflects biogeographic/population history (colonisation from South America, Pleistocene range contractions) rather than recent human impact; whereas (2) the underlying process of their isolation by distance pattern, which is best explained by 'isolation by dispersal limitation', implies contemporary gene flow between neighbouring populations as likely facilitated by the species' efficient seed dispersal capacity. Overall, these results underscore that even tree species from highly decimated forest regions may be genetically resilient to habitat fragmentation due to species-typical dispersal characteristics, the necessity of broad-scale measures for their conservation notwithstanding.

Strategies for improving approximate Bayesian computation tests for synchronous diversification

Background

Estimating the variability in isolation times across co-distributed taxon pairs that may have experienced the same allopatric isolating mechanism is a core goal of comparative phylogeography. The use of hierarchical Approximate Bayesian Computation (ABC) and coalescent models to infer temporal dynamics of lineage co-diversification has been a contentious topic in recent years. Key issues that remain unresolved include the choice of an appropriate prior on the number of co-divergence events (Ψ), as well as the optimal strategies for data summarization.

Methods

Through simulation-based cross validation we explore the impact of the strategy for sorting summary statistics and the choice of prior on Ψ on the estimation of co-divergence variability. We also introduce a new setting (β) that can potentially improve estimation of Ψ by enforcing a minimal temporal difference between pulses of co-divergence. We apply this new method to three empirical datasets: one dataset each of co-distributed taxon pairs of Panamanian frogs and freshwater fishes, and a large set of Neotropical butterfly sister-taxon pairs.

Results

We demonstrate that the choice of prior on Ψ has little impact on inference, but that sorting summary statistics yields substantially more reliable estimates of co-divergence variability despite violations of assumptions about exchangeability. We find the implementation of β improves estimation of Ψ, with improvement being most dramatic given larger numbers of taxon pairs. We find equivocal support for synchronous co-divergence for both of the Panamanian groups, but we find considerable support for asynchronous divergence among the Neotropical butterflies.

Conclusions

Our simulation experiments demonstrate that using sorted summary statistics results in improved estimates of the variability in divergence times, whereas the choice of hyperprior on Ψ has negligible effect. Additionally, we demonstrate that estimating the number of pulses of co-divergence across co-distributed taxon-pairs is improved by applying a flexible buffering regime over divergence times. This improves the correlation between Ψ and the true variability in isolation times and allows for more meaningful interpretation of this hyperparameter. This will allow for more accurate identification of the number of temporally distinct pulses of co-divergence that generated the diversification pattern of a given regional assemblage of sister-taxon-pairs.

Electronic supplementary material

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

An African Origin of the Eurylaimides (Passeriformes) and the Successful Diversification of the Ground-Foraging Pittas (Pittidae)

The Eurylaimides is one of the few passerine groups with a pantropical distribution. In this study, we generated a multi-calibrated tree with 83% of eurylaimid species diversity based on 30 molecular loci. Particular attention was given to the monotypic Sapayoidae to reconstruct the biogeography of this radiation. We conducted several topological tests including nonoverlapping subsampling of the concatenated alignment and coalescent species tree reconstruction. These tests firmly placed the South American Sapayoidae as the sister group to all other Eurylaimides families (split at ∼28 Ma), with increasing branch support as highly variable sites were removed. This topology is consistent with the breakup of the insular connection between Africa and South America (Atlantogea) that took place between the middle Eocene and the early Oligocene. We recovered Africa as the cradle of the core Eurylaimides, and this result is supported by all African lineages corresponding to the oldest splits within each family in this group. Our timescale suggests that desertification and the uplift of the Tibetan Plateau caused a parallel divergence between African and Asian lineages in all major clades in the core Eurylaimides at 22-9 Ma. We also propose that the ground-foraging behavior in the Pittidae ancestor allowed the pitta lineage to thrive and coexist with the older arboreal lineages of the core Eurylaimides. In contrast, the diversification of pittas in Australia was likely hindered by direct competition with the endemic ground-foraging oscines that had been well established in that continent since the Eocene.

Biodiversity and Topographic Complexity: Modern and Geohistorical Perspectives

Topographically complex regions on land and in the oceans feature hotspots of biodiversity that reflect geological influences on ecological and evolutionary processes. Over geologic time, topographic diversity gradients wax and wane over millions of years, tracking tectonic or climatic history. Topographic diversity gradients from the present day and the past can result from the generation of species by vicariance or from the accumulation of species from dispersal into a region with strong environmental gradients. Biological and geological approaches must be integrated to test alternative models of diversification along topographic gradients. Reciprocal illumination among phylogenetic, phylogeographic, ecological, paleontological, tectonic, and climatic perspectives is an emerging frontier of biogeographic research.

Timing and patterns of diversification in the Neotropical bat genus Pteronotus (Mormoopidae)

We investigate the biogeographic processes related to the origin and current patterns of distribution of the extant species of the genus Pteronotus. This clade of insectivorous bats is widely distributed in the Neotropical Region and has recently gone through a taxonomic update which increased more than twice its diversity. Using six molecular markers of 15 Pteronotus lineages ranging from Mexico to Central Brazil, we reconstruct a time-calibrated tree and infer molecular evolutionary rates for this bat genus. In addition, estimates of range evolution across phylogeny were obtained through statistical model testing among six different biogeographic models. The origin of the genus Pteronotus occurred approximately 16million years ago (Ma), with initial cladogenesis events being evenly distributed across the phylogeny. Divergence between most closely related species is recent, falling in the Pleistocene period less than 2.6Ma. Mainland lineages present congruent patterns of north versus south continent splitting while insular clades differ in their time of arrival in the Caribbean Islands. Temporal and geographic range estimates for early nodes of Pteronotus phylogeny suggest a central role of Neogene tectonic reorganizations of Central America in the group diversification process. Also, South American colonization by Pteronotus occurred early in the genus history. Founder-event speciation was an important mode of lineage splitting in Pteronotus, with two independent dispersal jumps having occurred to the Greater Antilles. Finally, Pleistocenic sea-level variation and climatic oscillations are possibly associated with divergence between sister-species and recent ages of MRCA for Pteronotus species.

Mitogenomics phylogenetic relationships of the current sloth's genera and species (Bradypodidae and Megalonychidae)

We sequenced the complete mitogenome of 39 sloths (19 Bradypus variegatus, 4 B. tridactylus, 1 B. pygmaeus, 1 B. torquatus, 4 Choloepus didactylus, and 10 C. hoffmanni). A Bayesian tree (BI) indicated a temporal split between Bradypus and Choloepus around 31 million years ago (MYA, Oligocene) and the other major splits within each genera during the Miocene and Pliocene. A haplotype network (MJN) estimated a lower temporal split between the sloth genera (around 23.5 MYA). Both methods detected the ancestor of B. torquatus as the first to diverge within Bradypus (21 for BI and 19 MJN), followed by that of the ancestor of B. tridactylus. The split of B. pygmaeus from the common ancestor with B. variegatus was around 12 MYA (BI) or 4.3 MYA (MJN). The splits among the previous populations of B. variegatus began around 8 MYA (BI) or 3.6 MYA (MJN). The trans-Andean population was the first to diverge from the remaining cis-Andean populations of B. variegatus. The genetic differentiation of the trans-Andean B. variegatus population relative to the cis-Andean B. variegatus is similar to that found for different species of sloths. The mitogenomic analysis resolved the differentiation of C. hoffmanni from the C. didactylus individuals of the Guiana Shield. However, one C. didactylus from the Colombian Amazon specimen was inside the C. hoffmanni clade. This could be the first example of possible natural hybridization in the Amazon of both Choloepus taxa or the existence of un-differentiable phenotypes of these two species in some Amazonian areas.

Early and dynamic colonization of Central America drives speciation in Neotropical army ants

The emergence of the Isthmus of Panama is one of the most important events in recent geological history, yet its timing and role in fundamental evolutionary processes remain controversial. While the formation of the isthmus was complete around 3 million years ago (Ma), recent studies have suggested prior intercontinental biotic exchange. In particular, the possibility of early intermittent land bridges facilitating colonization constitutes a potential mechanism for speciation and colonization before full closure of the isthmus. To test this hypothesis, we employed genomic methods to study the biogeography of the army ant genus Eciton, a group of keystone arthropod predators in Neotropical rainforests. Army ant colonies are unable to disperse across water and are therefore ideally suited to study the biogeographic impact of land bridge formation. Using a reduced representation genome sequencing approach, we show that all strictly Central American lineages of Eciton diverged from their respective South American sister lineage between 4 and 7 Ma, significantly prior to the complete closure of the isthmus. Furthermore, three of the lineage pairs form extensive and coincident secondary contact zones in Costa Rica and Nicaragua, with no evidence of gene flow. Such a discrete and repeated biogeographic pattern indicates at least two waves of army ant dispersal into Central America that were separated by significant genetic divergence times. Thus, by integrating phylogenomic, population genomic and geographic evidence, we show that early colonization of Central America across the emerging Isthmus of Panamá drove parallel speciation in Eciton army ants.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA₂ proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA₂ molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA₂s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA₂ (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA₂-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers.

SIGNIFICANCE

Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

Unveiling the Hidden Bat Diversity of a Neotropical Montane Forest

Mountain environments, characterized by high levels of endemism, are at risk of experiencing significant biodiversity loss due to current trends in global warming. While many acknowledge their importance and vulnerability, these ecosystems still remain poorly studied, particularly for taxa that are difficult to sample such as bats. Aiming to estimate the amount of cryptic diversity among bats of a Neotropical montane cloud forest in Talamanca Range—south-east Central America—, we performed a 15-night sampling campaign, which resulted in 90 captured bats belonging to 8 species. We sequenced their mitochondrial cytochrome c oxidase subunit I (COI) and screened their inter- and intraspecific genetic variation. Phylogenetic relations with conspecifics and closely related species from other geographic regions were established using Maximum Likelihood and Bayesian inference methods, as well as median-joining haplotype networks. Mitochondrial lineages highly divergent from hitherto characterized populations (> 9% COI dissimilarity) were found in Myotis oxyotus and Hylonycteris underwoodi. Sturnira burtonlimi and M. keaysi also showed distinct mitochondrial structure with sibling species and/or populations. These results suggest that mountains in the region hold a high degree of endemicity potential that has previously been ignored in bats. They also warn of the high extinction risk montane bats may be facing due to climatic change, particularly in isolated mountain systems like Talamanca Range.