Genetic structure and landscape effects on gene flow in the Neotropical lizard Norops brasiliensis (Squamata: Dactyloidae)

One key research goal of evolutionary biology is to understand the origin and maintenance of genetic variation. In the Cerrado, the South American savanna located primarily in the Central Brazilian Plateau, many hypotheses have been proposed to explain how landscape features (e.g., geographic distance, river barriers, topographic compartmentalization, and historical climatic fluctuations) have promoted genetic structure by mediating gene flow. Here, we asked whether these landscape features have influenced the genetic structure and differentiation in the lizard species Norops brasiliensis (Squamata: Dactyloidae). To achieve our goal, we used a genetic clustering analysis and estimate an effective migration surface to assess genetic structure in the focal species. Optimized isolation-by-resistance models and a simulation-based approach combined with machine learning (convolutional neural network; CNN) were then used to infer current and historical effects on population genetic structure through 12 unique landscape models. We recovered five geographically distributed populations that are separated by regions of lower-than-expected gene flow. The results of the CNN showed that geographic distance is the sole predictor of genetic variation in N. brasiliensis, and that slope, rivers, and historical climate had no discernible influence on gene flow. Our novel CNN approach was accurate (89.5%) in differentiating each landscape model. CNN and other machine learning approaches are still largely unexplored in landscape genetics studies, representing promising avenues for future research with increasingly accessible genomic datasets.

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

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

Species limits and diversification of the Dendropsophus rubicundulus subgroup (Anura, Hylidae) in Neotropical savannas

Understanding the processes that generate and maintain biodiversity at and below the species level is a central goal of evolutionary biology. Here we explore the spatial and temporal drivers of diversification of the treefrog subgroup Dendropsophus rubicundulus, a subgroup of the D. microcephalus species group, over periods of pronounced geological and climatic changes in the Neotropical savannas that they inhabit. This subgroup currently comprises 11 recognized species distributed across the Brazilian and Bolivian savannas, but the taxonomy has been in a state of flux, necessitating reexamination. Using newly generated single nucleotide polymorphism (SNP) data from restriction-site associated DNA sequencing (RADseq) and mitochondrial 16S sequence data for ∼150 specimens, we inferred phylogenetic relationships, tested species limits using a model-based approach, and estimated divergence times to gain insights into the geographic and climatic events that affected the diversification of this subgroup. Our results recognized at least nine species: D. anataliasiasi, D. araguaya, D. cerradensis, D. elianeae, D. jimi, D. rubicundulus, D. tritaeniatus, D. rozenmani, and D. sanborni. Although we did not collect SNP data for the latter two species, they are likely distinct based on mitochondrial data. In addition, we found genetic structure within the widespread species D. rubicundulus, which comprises three allopatric lineages connected by gene flow upon secondary contact. We also found evidence of population structure and perhaps undescribed diversity in D. elianeae, which warrants further study. The D. rubicundulus subgroup is estimated to have originated in the Late Miocene (∼5.45 million years ago), with diversification continuing through the Pliocene and Early Pleistocene, followed by the most recent divergence of D. rubicundulus lineages in the Middle Pleistocene. The epeirogenic uplift followed by erosion and denudation of the central Brazilian plateau throughout the Pliocene and Pleistocene, in combination with the increasing frequency and amplitude of climatic fluctuations during the Pleistocene, was important for generating and structuring diversity at or below the species level in the D. rubicundulus subgroup.

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

Peripheral isolation and demographic stability are reflected in the genetic diversity of the populations of the Helmeted Manakin: a bird endemic to the gallery forests

The Cerrado encompasses a complex network of hydrographic basins, which is responsible for the formation and maintenance of the riparian and gallery forests. Alterations in the vegetation resulting from the paleoclimatic changes that occurred during the Pleistocene influenced the current distribution of these humid forests. To understand of the evolutionary dynamics of this landscape on the population structure of the associated organisms, we studied the population genetics of the Antilophia galeata (Pipridae), a bird endemic to the gallery forests of the Cerrado. For this, we evaluated the variability of the mitochondrial control region of 71 A. galeata specimens from 18 localities distributed throughout the Cerrado. The results indicated the existence of significant substructure among the populations located in the peripheral areas of the Cerrado, in comparison with the populations found in the central portion, reflecting historical changes in the environment. We also found high levels of diversity in all the populations, while the analysis of the demographic history revealed a scenario of stability. Overall, then, these findings indicate that the historical modifications of the gallery forest distribution may have been most accentuated in the peripheral zones, with a greater stability and connectivity persisting in the central portion of the biome.

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

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

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

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

Molecular data reveal multiple lineages of Scinax nebulosus (Spix, 1824) (Anura: Hylidae) with Plio-Pleistocene diversification in different Brazilian regions

To understand the organism's history, we can start assessing the complexity of the biome where they occur. In this study, we used a region of the mitochondrial genome, the rRNA 16S, to evaluate the genetic differentiation in Scinax nebulosus along with its geographical range highlighting important Brazilian biomes as Restinga, Cerrado, Amazon, and Atlantic Forest. Geographically structured genetic divergence was observed within the species S. nebulosus. The values of the fixation index (Фst) and the pairwise Fst index were high and significant regarding this structuring. Besides, the haplotype network corroborates these results with the haplotypes arrangement found by separating the S. nebulosus populations in two major groups: North and Northeast. The lineage delimitation analyses indicate the occurrence of several lineages with divergence mainly between the samples from the Northeast group. Thus, we can suggest that S. nebulosus may present itself as a group of cryptic species due to the genetic characteristics found. The existence of a mosaic of heterogeneous habitats may explain the genetic divergence found, which justifies the existence of cryptic species in this group. However, this hypothesis needs more detail in molecular studies, including large sample sizes and other population and demographic analyses.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

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.

Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes

Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.

Phylogeographic model selection using convolutional neural networks

The discipline of phylogeography has evolved rapidly in terms of the analytical toolkit used to analyse large genomic data sets. Despite substantial advances, analytical tools that could potentially address the challenges posed by increased model complexity have not been fully explored. For example, deep learning techniques are underutilized for phylogeographic model selection. In non-model organisms, the lack of information about their ecology and evolution can lead to uncertainty about which demographic models are appropriate. Here, we assess the utility of convolutional neural networks (CNNs) for assessing demographic models in South American lizards in the genus Norops. Three demographic scenarios (constant, expansion, and bottleneck) were considered for each of four inferred population-level lineages, and we found that the overall model accuracy was higher than 98% for all lineages. We then evaluated a set of 26 models that accounted for evolutionary relationships, gene flow, and changes in effective population size among the four lineages, identifying a single model with an estimated overall accuracy of 87% when using CNNs. The inferred demography of the lizard system suggests that gene flow between non-sister populations and changes in effective population sizes through time, probably in response to Pleistocene climatic oscillations, have shaped genetic diversity in this system. Approximate Bayesian computation (ABC) was applied to provide a comparison to the performance of CNNs. ABC was unable to identify a single model among the larger set of 26 models in the subsequent analysis. Our results demonstrate that CNNs can be easily and usefully incorporated into the phylogeographer's toolkit.

Taxonomic revision of the genus Xenopholis Peters, 1869 (Serpentes: Dipsadidae): Integrating morphology with ecological niche

A reliable identification and delimitation of species is an essential pre-requisite for many fields of science and conservation. The Neotropical herpetofauna is the world's most diverse, including many taxa of uncertain or debated taxonomy. Here we tackle one such species complex, by evaluating the taxonomic status of species currently allocated in the snake genus Xenopholis (X. scalaris, X. undulatus, and X. werdingorum). We base our conclusions on concordance between quantitative (meristic and morphometric) and qualitative (color pattern, hemipenes and skull features) analyses of morphological characters, in combination with ecological niche modeling. We recognize all three taxa as valid species and improve their respective diagnosis, including new data on color in life, pholidosis, bony morphology, and male genitalia. We find low overlap among the niches of each species, corroborating the independent source of phenotypic evidence. Even though all three species occur in the leaf litter of distinct forested habitats, Xenopholis undulatus is found in the elevated areas of the Brazilian Shield (Caatinga, Cerrado and Chaco), whereas X. scalaris occurs in the Amazon and Atlantic rainforests, and X. werdingorum in the Chiquitanos forest and Pantanal wetlands. We discuss the disjunct distribution between Amazonian and Atlantic Forest snake species in the light of available natural history and ecological aspects. This study shows the advantages of combining multiple data sources for reliable identification and circumscription of ecologically similar species.

Dispersal and local persistence shape the genetic structure of a widespread Neotropical plant species with a patchy distribution

BACKGROUND AND AIMS

Isolated populations constitute an ideal laboratory to study the consequences of intraspecific divergence, because intrinsic incompatibilities are more likely to accumulate under reduced gene flow. Here, we use a widespread bromeliad with a patchy distribution, Pitcairnia lanuginosa, as a model to infer processes driving Neotropical diversification and, thus, to improve our understanding of the origin and evolutionary dynamics of biodiversity in this highly speciose region.

METHODS

We assessed the timing of lineage divergence, genetic structural patterns and historical demography of P. lanuginosa, based on microsatellites, and plastid and nuclear sequence data sets using coalescent analyses and an Approximate Bayesian Computation framework. Additionally, we used species distribution models (SDMs) to independently estimate potential changes in habitat suitability.

KEY RESULTS

Despite morphological uniformity, plastid and nuclear DNA data revealed two distinct P. lanuginosa lineages that probably diverged through dispersal from the Cerrado to the Central Andean Yungas, following the final uplift of the Andes, and passed through long-term isolation with no evidence of migration. Microsatellite data indicate low genetic diversity and high levels of inbreeding within populations, and restricted gene flow among populations, which are likely to be a consequence of bottlenecks (or founder events), and high selfing rates promoting population persistence in isolation. SDMs showed a slight expansion of the suitable range for P. lanuginosa lineages during the Last Glacial Maximum, although molecular data revealed a signature of older divergence. Pleistocene climatic oscillations thus seem to have played only a minor role in the diversification of P. lanuginosa, which probably persisted through adverse conditions in riparian forests.

CONCLUSIONS

Our results imply drift as a major force shaping the evolution of P. lanuginosa, and suggest that dispersal events have a prominent role in connecting Neotropical open and forest biomes.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

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

Historical demography of four gecko species specializing in boulder cave habitat: Implications in the evolutionary dead end hypothesis and conservation

Specialization in narrow ecological niches may not only help species to survive in competitive or unique environments but also contribute to their extermination over evolutionary time. Although the "evolutionary dead end" hypothesis has long been debated, empirical evidence from species with detailed information on niche specialization and evolutionary history remains rare. Here we use a group of four closely related Cnemaspis gecko species that depend highly on granite boulder caves in the Mekong Delta to investigate the potential impact of ecological specialization on their evolution and population dynamics. Isolated by unsuitable floodplain habitats, these boulder-dwelling geckos are among the most narrowly distributed Squamata in the world. We applied several coalescence-based approaches combined with the RAD-seq technique to estimate their divergence times, gene flow and demographic fluctuations during the speciation and population differentiation processes. Our results reveal long-term population shrinkage in the four geckos and limited gene flow during their divergence. The results suggest that the erosion and fragmentation of the granite boulder hills have greatly impacted population divergence and declines. The habitat specialization of these geckos has led to fine-scaled speciation in these granite rocky hills; in contrast, specialization might also have pushed these species toward the edge of extinction. Our study also emphasizes the conservation urgency of these vulnerable, cave-dependent geckos.

Diversification of the widespread neotropical frog Physalaemus cuvieri in response to Neogene-Quaternary geological events and climate dynamics

Here we reconstructed the demographical history and the dispersal dynamics of Physalemus cuvieri through the Neogene-Quaternary periods by coupling DNA regions with different mutation rates, ecological niche modelling, reconstruction of spatio-temporal lineage dispersal and coalescent simulations. Still, to test alternative diversification scenarios we used approximate Bayesian computation. Molecular phylogenetic analysis recovered four deep and strongly supported clades, which we interpret as population lineages. The ancestral location reconstruction placed the root in southcentral Amazonia, and the dispersal events indicate that spatial displacement was widespread early in the diversification of this species. The demographical scenario of "Multiple Refugia" with recent lineage admixture was the most likely hypothesis to predict the observed genetic parameters of P. cuvieri. Our results revealed that Neogene orogenic events might have played a prominent role in the early diversification of P. cuvieri. The species shows deep divergences with strong regional population structure, despite its widespread distribution. Final uplift of the central Brazilian Plateau and formation of the river basins in Central South America played an important role in the origin, diversification and the maintenance of P. cuvieri lineages.

Phylogeographic analyses of the pampas cat (Leopardus colocola; Carnivora, Felidae) reveal a complex demographic history

The pampas cat is a small felid that occurs in open habitats throughout much of South America. Previous studies have revealed intriguing patterns of morphological differentiation and genetic structure among its populations, as well as molecular evidence for hybridization with the closely related L. tigrinus. Here we report phylogeographic analyses encompassing most of its distribution (focusing particularly on Brazilian specimens, which had been poorly sampled in previous studies), using a novel dataset comprising 2,143 bp of the mitogenome, along with previously reported mtDNA sequences. Our data revealed strong population strutucture and supported a west-to-east colonization process in this species' history. We detected two population expansion events, one older (ca. 200 thousand years ago [kya]) in western South America and another more recent (ca. 60-50 kya) in eastern areas, coinciding with the expansion of savanna environments in Brazil. Analyses including L. tigrinus individuals bearing introgressed mtDNA from L. colocola showed a complete lack of shared haplotypes between species, indicating that their hybridization was ancient. Finally, we observed a close relationship between Brazilian/Uruguayan L. colocola haplotypes and those sampled in L. tigrinus, indicating that their hybridization was likely related to the demographic expansion of L. colocola into eastern South America.

Three new species of the killifish genus Melanorivulus from the central Brazilian Cerrado savanna (Cyprinodontiformes, Aplocheilidae)

Three new species are described from the Neotropical region comprising the Cerrado savannas of the central Brazilian plateaus, which is among the most important biodiversity centres in the world. These species are considered closely related to Melanorivulusdapazi from the same region, with which they share the presence of a rudimentary interarcual cartilage and a dark reddish brown distal margin on the male anal fin. The group comprising Melanorivulusdapazi and the three new species is here named as the Melanorivulusdapazi species group. Melanorivulusignescenssp. n., from the upper Rio Araguaia basin, is distinguished from all other species of the Melanorivulusdapazi group by the anal-fin colour pattern in males; Melanorivulusflavipinnissp. n. and Melanorivulusregularissp. n. from the Rio Paraguai basin are distinguished from all other congeners of the Melanorivulusdapazi group by the colour pattern of the caudal fin and number of scales in the longitudinal series, respectively. All the new species are further unambiguously diagnosed by unique combinations of morphological characters, including meristic and morphometric data, and colour patterns. This study reinforces the importance of using live colour patterns to diagnose species and species groups of the genus Melanorivulus, but also indicates that osteological characters may be informative for species diagnosis. This study confirms the high diversity of species of Melanorivulus in the central Brazilian Cerrado plateaus already reported in previous studies, indicating that endemic species are often restricted to short segments of a single river drainage.

A large historical refugium explains spatial patterns of genetic diversity in a Neotropical savanna tree species

BACKGROUND AND AIMS

The relative role of Pleistocene climate changes in driving the geographic distribution and genetic diversity of South American species is not well known, especially from open biomes such as the Cerrado, the most diverse tropical savanna, encompassing high levels of endemism. Here the effects of Quaternary climatic changes on demographic history, distribution dynamics and genetic diversity of Dimorphandra mollis, an endemic tree species widely distributed in the Cerrado, were investigated.

METHODS

A total of 38 populations covering most of the distribution of D. mollis were analysed using internal transcribed spacer (ITS) sequences and nuclear microsatellite variation [simple sequence repeats (SSRs)]. The framework incorporated statistical phylogeography, coalescent analyses and ecological niche modelling (ENM).

KEY RESULTS

Different signatures of Quaternary climatic changes were found for ITS sequences and SSRs corresponding to different time slices. Coalescent analyses revealed large and constant effective population sizes, with high historical connectivity among the populations for ITS sequences and low effective population sizes and gene flow with recent population retraction for SSRs. ENMs indicated a slight geographical range retraction during the Last Glacial Maximum. A large historical refugium across central Brazil was predicted. Spatially explicit analyses showed a spatial cline pattern in genetic diversity related to the paleodistribution of D. mollis and to the centre of its historical refugium.

CONCLUSIONS

The complex genetic patterns found in D. mollis are the result of a slight geographical range retraction during the Last Glacial Maximum followed by population expansion to the east and south from a large refugium in the central part of the Cerrado. This historical refugium is coincident with an area predicted to be climatically stable under future climate scenarios. The identified refugium should be given high priority in conservation polices to safeguard the evolutionary potential of the species under predicted future climatic changes.

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

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

Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America

The Pilosocereus aurisetus complex consists of eight cactus species with a fragmented distribution associated to xeric enclaves within the Cerrado biome in eastern South America. The phylogeny of these species is incompletely resolved, and this instability complicates evolutionary analyses. Previous analyses based on both plastid and microsatellite markers suggested that this complex contained species with inherent phylogeographic structure, which was attributed to recent diversification and recurring range shifts. However, limitations of the molecular markers used in these analyses prevented some questions from being properly addressed. In order to better understand the relationship among these species and make a preliminary assessment of the genetic structure within them, we developed anonymous nuclear loci from pyrosequencing data of 40 individuals from four species in the P. aurisetus complex. The data obtained from these loci were used to identify genetic clusters within species, and to investigate the phylogenetic relationship among these inferred clusters using a species tree methodology. Coupled with a palaeodistributional modelling, our results reveal a deep phylogenetic and climatic disjunction between two geographic lineages. Our results highlight the importance of sampling more regions from the genome to gain better insights on the evolution of species with an intricate evolutionary history. The methodology used here provides a feasible approach to develop numerous genealogical molecular markers throughout the genome for non-model species. These data provide a more robust hypothesis for the relationship among the lineages of the P. aurisetus complex.

Molecular data reveal spatial and temporal patterns of diversification and a cryptic new species of lowland Stenocercus Duméril & Bibron, 1837 (Squamata: Tropiduridae)

Phylogenetic studies have uncovered biogeographic patterns and the associated diversification processes of Neotropical wet forest taxa, yet the extensive open and drier biomes have received much less attention. In the Stenocercus lizard radiation, restricted sampling and phylogenetic information have limited inferences about the timing, spatial context, and environmental drivers of diversification in the open and dry lowland settings of eastern and southern South America. Based on new DNA sequence data of previously unsampled species, we provide an updated historical biogeographic hypothesis of Stenocercus. We infer phylogenetic relationships, estimate divergence times, and track ancestral distributions, asking whether cladogenetic events within the genus correlate to reported shifts in South American landscapes during the past 30millionyears, focusing in the open and drier areas. To examine correlations between genetic and ecological divergence, we extracted environmental data from occurrence records and estimated climatic envelopes occupied by lowland taxa. Our results suggest that Stenocercus began to diversify around the South American Midwest by the late Oligocene. We recovered two main lowland and two main Andean clades within the genus; within both Andean clades, most cladogenetic events date back to the Miocene, synchronously with the most intense phase of Andean uplift. In the western clade of lowland Stenocercus, species ranges and divergence times are consistent with major landscape shifts at the upper Guaporé and Paraguay River basins as a result of Andean orogeny, suggesting vicariant speciation. By contrast, in the 'horned' lowland clade, we find evidence that dispersal and ecological differentiation have shaped species divergences and current ranges in the Brazilian Cerrado, Caatinga, Pampas and Atlantic Forest, possibly under a vanishing refuge scenario. Lastly, our phylogenetic results indicate two divergent clades within the formerly recognized taxon S. sinesaccus, and further evaluation of morphological data corroborates the existence of a distinct, new species of Stenocercus, here described. The new taxon occurs in the Chapada dos Parecis massif in the Brazilian states of Mato Grosso and Rondônia.

Brazilian Anopheles darlingi Root (Diptera: Culicidae) Clusters by Major Biogeographical Region

The major drivers of the extensive biodiversity of the Neotropics are proposed to be geological and tectonic events together with Pliocene and Pleistocene environmental and climatic change. Geographical barriers represented by the rivers Amazonas/Solimões, the Andes and the coastal mountain ranges in eastern Brazil have been hypothesized to lead to diversification within the primary malaria vector, Anopheles (Nyssorhynchus) darlingi Root, which primarily inhabits rainforest. To test this biogeographical hypothesis, we analyzed 786 single nucleotide polymorphisms (SNPs) in 12 populations of An. darlingi from across the complex Brazilian landscape. Both model-based (STRUCTURE) and non-model-based (Principal Components and Discriminant Analysis) analysis of population structure detected three major genetic clusters that correspond with newly described Neotropical biogeographical regions: 1) Atlantic Forest province (= southeast population); 2) Parana Forest province (= West Atlantic forest population, with one Chacoan population - SP); and 3) Brazilian dominion population (= Amazonian population with one Chacoan population - TO). Significant levels of pairwise genetic divergences were found among the three clusters, allele sharing among clusters was negligible, and geographical distance did not contribute to differentiation. We infer that the Atlantic forest coastal mountain range limited dispersal between the Atlantic Forest province and the Parana Forest province populations, and that the large, diagonal open vegetation region of the Chacoan dominion dramatically reduced dispersal between the Parana and Brazilian dominion populations. We hypothesize that the three genetic clusters may represent three putative species.

Biogeographic history and cryptic diversity of saxicolous Tropiduridae lizards endemic to the semiarid Caatinga

Background

Phylogeographic research has advanced in South America, with increasing efforts on taxa from the dry diagonal biomes. However, the diversification of endemic fauna from the semiarid Caatinga biome in northeastern Brazil is still poorly known. Here we targeted saxicolous lizards of the Tropidurus semitaeniatus species group to better understand the evolutionary history of these endemic taxa and the Caatinga. We estimated a time-calibrated phylogeny for the species group based on two mitochondrial and two nuclear genes and jointly estimated the species limits and species tree within the group. We also devoted a denser phylogeographic sampling of the T. semitaeniatus complex to explore migration patterns, and the spatiotemporal diffusion history to verify a possible role of the São Francisco River as a promoter of differentiation in this saxicolous group of lizards.

Results

Phylogenetic analysis detected high cryptic genetic diversity, occurrence of unique microendemic lineages associated with older highlands, and a speciation history that took place during the Pliocene-Pleistocene transition. Species delimitation detected five evolutionary entities within the T. semitaeniatus species group, albeit with low support. Thus, additional data are needed for a more accurate definition of species limits and interspecific relationships within this group. Spatiotemporal analyses reconstructed the geographic origin of the T. semitaeniatus species complex to be located north of the present-day course of the São Francisco River, followed by dispersal that expanded its distribution towards the northwest and south. Gene flow estimates showed higher migration rates into the lineages located north of the São Francisco River.

Conclusions

The phylogenetic and population structures are intrinsically associated with stable rock surfaces and landscape rearrangements, such as the establishment of drainage basins located to the northern and southern distribution ranges. The T. semitaeniatus complex preserved high genetic diversity during range expansion, possibly as a result of frequent long-distance dispersal events. Our results indicate that both the current course of the São Francisco River and its paleo-courses had an important role in promoting diversification of the Caatinga endemic T. semitaeniatus species group.

Electronic supplementary material

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

Revisiting the vanishing refuge model of diversification

Much of the debate around speciation and historical biogeography has focused on the role of stabilizing selection on the physiological (abiotic) niche, emphasizing how isolation and vicariance, when associated with niche conservatism, may drive tropical speciation. Yet, recent re-emphasis on the ecological dimensions of speciation points to a more prominent role of divergent selection in driving genetic, phenotypic, and niche divergence. The vanishing refuge model (VRM), first described by Vanzolini and Williams (1981), describes a process of diversification through climate-driven habitat fragmentation and exposure to new environments, integrating both vicariance and divergent selection. This model suggests that dynamic climates and peripheral isolates can lead to genetic and functional (i.e., ecological and phenotypic) diversity, resulting in sister taxa that occupy contrasting habitats with abutting distributions. Here, we provide predictions for populations undergoing divergence according to the VRM that encompass habitat dynamics, phylogeography, and phenotypic differentiation across populations. Such integrative analyses can, in principle, differentiate the operation of the VRM from other speciation models. We applied these principles to a lizard species, Coleodactylus meridionalis, which was used to illustrate the model in the original paper. We incorporate data on inferred historic habitat dynamics, phylogeography and thermal physiology to test for divergence between coastal and inland populations in the Atlantic Forest of Brazil. Environmental and genetic analyses are concordant with divergence through the VRM, yet physiological data are not. We emphasize the importance of multidisciplinary approaches to test this and alternative speciation models while seeking to explain the extraordinarily high genetic and phenotypic diversity of tropical biomes.

Karyotypic diversity in seven Amazonian anurans in the genus Hypsiboas (family Hylidae)

Background

Hypsiboas species have been divided into seven groups using morphological and genetic characters, but for most of the species, there is no cytogenetic information available. A cytogenetic analysis using conventional staining, C-banding, silver staining, and fluorescence in situ hybridization (FISH) with telomeric sequence probes were used to investigate the karyotype of seven Amazon species of the genus Hypsiboas belonging to the following intrageneric groups: H. punctatus (H. cinerascens), H. semilineatus (H. boans, H. geographicus, and H. wavrini), and H. albopunctatus (H. lanciformis, H. multifasciatus, and H. raniceps). The aim was to differentiate between the karyotypes and use the chromosomal markers to distinguish between the Hypsiboas groups. The data were compared with a previous phylogenetic proposal for these anurans. In addition, H. lanciformis, H. boans, and H. wavrini are described here for the first time, and we characterize the diploid numbers for H. cinerascens, H. geographicus, H. multifasciatus, and H. raniceps.

Results

The diploid number for all of the species analyzed was 24, with the exception of Hypsiboas lanciformis, which had 2n = 22 chromosomes. The constitutive heterochromatin distribution, nucleolar organizer region locations, and interstitial telomeric sites differed between the species. A hypothesis that the heterochromatic patterns are evolving is proposed, with the divergence of the groups probably involving events such as an increase in the heterochromatin in the species of the H. semilineatus group. The FISH conducted with the telomeric probes detected sites in the terminal regions of all of the chromosomes of all species. Interstitial telomeric sites were detected in three species belonging to the H. semilineatus group: H. boans, H. geographicus, and H. wavrini.

Conclusion

The results of this study reinforce the complexity previously observed within the genus Hypsiboas and in the different groups that compose this taxon. More studies are needed focusing on this group and covering larger sampling areas, especially in the Brazilian Amazon, to improve our understanding of this fascinating and complex group.

Climate, physiological tolerance and sex-biased dispersal shape genetic structure of Neotropical orchid bees

Understanding the impact of past climatic events on the demographic history of extant species is critical for predicting species' responses to future climate change. Palaeoclimatic instability is a major mechanism of lineage diversification in taxa with low dispersal and small geographical ranges in tropical ecosystems. However, the impact of these climatic events remains questionable for the diversification of species with high levels of gene flow and large geographical distributions. In this study, we investigate the impact of Pleistocene climate change on three Neotropical orchid bee species (Eulaema bombiformis, E. meriana and E. cingulata) with transcontinental distributions and different physiological tolerances. We first generated ecological niche models to identify species-specific climatically stable areas during Pleistocene climatic oscillations. Using a combination of mitochondrial and nuclear markers, we inferred calibrated phylogenies and estimated historical demographic parameters to reconstruct the phylogeographical history of each species. Our results indicate species with narrower physiological tolerance experienced less suitable habitat during glaciations and currently exhibit strong population structure in the mitochondrial genome. However, nuclear markers with low and high mutation rates show lack of association with geography. These results combined with lower migration rate estimates from the mitochondrial than the nuclear genome suggest male-biased dispersal. We conclude that despite large effective population sizes and capacity for long-distance dispersal, climatic instability is an important mechanism of maternal lineage diversification in orchid bees. Thus, these Neotropical pollinators are susceptible to disruption of genetic connectivity in the event of large-scale climatic changes.

The role of historical barriers in the diversification processes in open vegetation formations during the Miocene/Pliocene using an ancient rodent lineage as a model

The Neotropics harbors a high diversity of species and several hypotheses have been proposed to account for this pattern. However, while species of forested domains are frequently studied, less is known of species from open vegetation formations occupying, altogether, a larger area than the Amazon Forest. Here we evaluate the role of historical barriers and the riverine hypothesis in the speciation patterns of small mammals by analyzing an ancient rodent lineage (Thrichomys, Hystricomorpha). Phylogenetic and biogeographic analyses were carried out with mitochondrial and nuclear DNA markers to analyze the evolutionary relationships between Thrichomys lineages occurring in dry domains along both banks of the Rio São Francisco. This river is one of the longest of South America whose course and water flow have been modified by inland tectonic activities and climate changes. Molecular data showed a higher number of lineages than previously described. The T. inermis species complex with 2n = 26, FN = 48 was observed in both banks of the river showing a paraphyletic arrangement, suggesting that river crossing had occurred, from east to west. A similar pattern was also observed for the T. apereoides complex. Thrichomys speciation occurred in Late Miocene when the river followed a different course. The current geographic distribution of Thrichomys species and their phylogenetic relationships suggested the existence of frequent past connections between both banks in the middle section of the Rio São Francisco. The extensive palaeodune region found in this area has been identified as a centre of endemism of several vertebrate species and is likely to be a center of Thrichomys diversification.

A geographic mosaic of evolutionary lineages within the insular endemic newt Euproctus montanus

Islands are hotspots of biodiversity, with a disproportionately high fraction of endemic lineages, often of ancient origin. Nevertheless, intra-island phylogeographies are surprisingly scarce, leading to a scanty knowledge about the microevolutionary processes induced on island populations by Plio-Pleistocene climatic oscillations, and the manner in which these processes contributed to shape their current genetic diversity. We investigated the phylogeography, historical demography and species distribution models of the Corsican endemic newt Euproctus montanus (north-western Mediterranean). As for many island endemics, the continuous distribution of E. montanus throughout its range has hitherto been considered as evidence for a single large population, a belief that also guided the species' categorization for conservation purposes. Instead, we found a geographic mosaic of ancient evolutionary lineages, with five main clades of likely Pliocene origin (2.6-5.8 My), all but one restricted to northern Corsica. Moreover, the copresence between main lineages in the same population was limited to a single case. As also suggested by growing literature on intra-island phylogeographic variation, it seems that the extensive use of simplifying assumption on the population structure and historical demography of island populations--both in theoretical and applicative studies--should be carefully reconsidered, a claim that is well exemplified by the case presented here.