Geoclimatic drivers of diversification in the largest arid and semi-arid environment of the Neotropics: Perspectives from phylogeography

The South American Dry Diagonal, also called the Diagonal of Open Formations, is a large region of seasonally dry vegetation extending from northeastern Brazil to northern Argentina, comprising the Caatinga, Cerrado, and Chaco subregions. A growing body of phylogeography literature has determined that a complex history of climatic changes coupled with more ancient geological events has produced a diverse and endemic-rich Dry Diagonal biota. However, the exact drivers are still under investigation, and their relative strengths and effects are controversial. Pleistocene climatic fluctuations structured lineages via vegetation shifts, refugium formation, and corridors between the Amazon and Atlantic forests. In some taxa, older geological events, such as the reconfiguration of the São Francisco River, uplift of the Central Brazilian Plateau, or the Miocene inundation of the Chaco by marine incursions, were more important. Here, we review the Dry Diagonal phylogeography literature, discussing each hypothesized driver of diversification and assessing degree of support. Few studies statistically test these hypotheses, with most support drawn from associating encountered phylogeographic patterns such as population structure with the timing of ancient geoclimatic events. Across statistical studies, most hypotheses are well supported, with the exception of the Pleistocene Arc Hypothesis. However, taxonomic and regional biases persist, such as a proportional overabundance of herpetofauna studies, and the under-representation of Chaco studies. Overall, both Pleistocene climate change and Neogene geological events shaped the evolution of the Dry Diagonal biota, though the precise effects are regionally and taxonomically varied. We encourage further use of model-based analyses to test evolutionary scenarios, as well as interdisciplinary collaborations to progress the field beyond its current focus on the traditional set of geoclimatic hypotheses.

The desert wolf-spider genus Xenoctenus: two endemic species from the Brazilian Caatinga, and a redescription of the type-species, X. unguiculatus (Araneae: Xenoctenidae)

The Caatinga is a nucleus of seasonally dry tropical forest (SDTF) known as a hotspot of diversification and endemism. Despite its importance, this biome is still insufficiently sampled, resulting in extensive knowledge gaps regarding its species richness and composition. In this study we report two species of Xenoctenidae that are endemic to, and widely distributed in the Caatinga. We redescribe and illustrate Odo vittatus (Mello-Leito, 1936), the only xenoctenid species previously known from the Caatinga. We transfer this species to Xenoctenus Mello-Leito,1938, a genus currently known from six species restricted to Argentina, Bolivia, and Colombia. We also newly describe the male of Xenoctenus vittatus comb. nov. and provide new records of this species, which was hitherto known only from the type-locality, throughout the Caatinga and nearby semiarid vegetation formations. We also describe and illustrate a new species, Xenoctenus kaatinga sp. nov., based on males and female specimens collected throughout the Caatinga. Additionally, we propose diagnostic characters for Xenoctenus and redescribe the type-species, X. unguiculatus.

Assessing model adequacy for Bayesian Skyline plots using posterior predictive simulation

Bayesian skyline plots (BSPs) are a useful tool for making inferences about demographic history. For example, researchers typically apply BSPs to test hypotheses regarding how climate changes have influenced intraspecific genetic diversity over time. Like any method, BSP has assumptions that may be violated in some empirical systems (e.g., the absence of population genetic structure), and the naïve analysis of data collected from these systems may lead to spurious results. To address these issues, we introduce P2C2M.Skyline, an R package designed to assess model adequacy for BSPs using posterior predictive simulation. P2C2M.Skyline uses a phylogenetic tree and the log file output from Bayesian Skyline analyses to simulate posterior predictive datasets and then compares this null distribution to statistics calculated from the empirical data to check for model violations. P2C2M.Skyline was able to correctly identify model violations when simulated datasets were generated assuming genetic structure, which is a clear violation of BSP model assumptions. Conversely, P2C2M.Skyline showed low rates of false positives when models were simulated under the BSP model. We also evaluate the P2C2M.Skyline performance in empirical systems, where we detected model violations when DNA sequences from multiple populations were lumped together. P2C2M.Skyline represents a user-friendly and computationally efficient resource for researchers aiming to make inferences from BSP.

Incorporating Topological and Age Uncertainty into Event-Based Biogeography of Sand Spiders Supports Paleo-Islands in Galapagos and Ancient Connections among Neotropical Dry Forests

Event-based biogeographic methods, such as dispersal-extinction-cladogenesis, have become increasingly popular for attempting to reconstruct the biogeographic history of organisms. Such methods employ distributional data of sampled species and a dated phylogenetic tree to estimate ancestral distribution ranges. Because the input tree is often a single consensus tree, uncertainty in topology and age estimates are rarely accounted for, even when they may affect the outcome of biogeographic estimates. Even when such uncertainties are taken into account for estimates of ancestral ranges, they are usually ignored when researchers compare competing biogeographic hypotheses. We explore the effect of incorporating this uncertainty in a biogeographic analysis of the 21 species of sand spiders (Sicariidae: Sicarius) from Neotropical xeric biomes, based on a total-evidence phylogeny including a complete sampling of the genus. Using a custom R script, we account for uncertainty in ages and topology by estimating ancestral ranges over a sample of trees from the posterior distribution of a Bayesian analysis, and for uncertainty in biogeographic estimates by using stochastic maps. This approach allows for counting biogeographic events such as dispersal among areas, counting lineages through time per area, and testing biogeographic hypotheses, while not overestimating the confidence in a single topology. Including uncertainty in ages indicates that Sicarius dispersed to the Galapagos Islands when the archipelago was formed by paleo-islands that are now submerged; model comparison strongly favors a scenario where dispersal took place before the current islands emerged. We also investigated past connections among currently disjunct Neotropical dry forests; failing to account for topological uncertainty underestimates possible connections among the Caatinga and Andean dry forests in favor of connections among Caatinga and Caribbean + Mesoamerican dry forests. Additionally, we find that biogeographic models including a founder-event speciation parameter (“+J”) are more prone to suffer from the overconfidence effects of estimating ancestral ranges using a single topology. This effect is alleviated by incorporating topological and age uncertainty while estimating stochastic maps, increasing the similarity in the inference of biogeographic events between models with or without a founder-event speciation parameter. We argue that incorporating phylogenetic uncertainty in biogeographic hypothesis-testing is valuable and should be a commonplace approach in the presence of rogue taxa or wide confidence intervals in age estimates, and especially when using models including founder-event speciation.

Cryptic diversity within three South American whip spider species (Arachnida, Amblypygi)

Cryptic diversity (CD), the presence of highly divergent phylogenetic lineages within closed morphological species, has been documented for many taxa. Great arachnid orders such as Araneae or Scorpiones are well studied and many cases of CD have been described therein; to date, however, related research on smaller arachnid orders, such as whip spiders (Amblypygi), remains lacking. In the current study, we investigated CD based on cytochrome oxidase 1 (COI) in three nominal species of the genus Heterophrynus (H. alces, H. batesii, and H. longicornis), represented by 65 specimens. The sequences were compared using three different methods. All three species showed geographically structured CD. Thus, given its existence in this genus, it is important that CD and its spatial distribution be considered in future studies and possible conservation projects.

Genetics and Distribution Modeling: The Demographic History of the Cactophilic Drosophila buzzatii Species Cluster in Open Areas of South America

Biodiversity is the result of historical and recurrent events acting on populations and species. The Drosophila buzzatii species cluster is distributed along a diagonal of open areas in South America. Combining genetic analyses with species distribution models we evaluated the influence of climatic changes in the demography history of this cluster. We performed a Bayesian Skyline analysis and reconstructed the ancestral areas based on mitochondrial cytochrome oxidase I (mtCOI) gene. We modeled the current distribution and projected it to past (mid-Holocene and Last Glacial Maximum) and future. Our results demonstrate that climate change plays a critical role in historical demography and in defining the current and future geographic ranges of these species. The inter-Andean dry valleys and the Chiquitano Seasonally Dry Tropical Forests (SDTF) in Bolivia are considered the ancestral area for the D. buzzatii cluster. From this area, the migration route was through a west-east corridor from central Andes, throughout Bolivia, Paraguay and Argentina toward eastern and northeastern Brazil, along the edges of rain forest. The responses of the species to the climatic changes differ from the dominant assumptions of expansion during dry/cold weather and contraction during wet/warm weather that characterized the glacial cycles. We suggest that the influence of ecology and the potential responsiveness of each taxon to the environmental dynamics should be considered as well. Predictions for the future suggest a decline of suitable areas for the cluster, threatening biodiversity of these habitats. This work showed the importance of an integrative analysis of genetics and geography information to improve the inferences about demographic history hypotheses for the cluster.

Dwarf geckos and giant rivers: the role of the São Francisco River in the evolution of Lygodactylus klugei (Squamata: Gekkonidae) in the semi-arid Caatinga of north-eastern Brazil

Species diversification can be strongly influenced by geomorphological features, such as mountains, valleys and rivers. Rivers can act as hard or soft barriers to gene flow depending on their size, speed of flow, historical dynamics and regional topographical characteristics. The São Francisco River (SFR) is the largest perennial river in the Caatinga biome in north-eastern Brazil and has been considered a barrier to gene flow and dispersal. Herein, we evaluated the role of the SFR on the evolution of Lygodactylus klugei, a small gecko from the Caatinga. Using a single-locus species delimitation method (generalized mixed Yule coalescent), we defined lineages (haploclades). Subsequently, we evaluated the role of the SFR in structuring genetic diversity in this species using a multilocus approach to quantify migration across margins. We also evaluated genetic structure based on nuclear markers, testing the number of populations found through an assignment test (STRUCTURE) across the species distribution. We recovered two mitochondrial lineages structured with respect to the SFR, but only a single population was inferred from nuclear markers. Given that we detected an influence of the SFR only on mitochondrial markers, we suggest that the current river course has acted as a relatively recent geographical barrier for L. klugei, for ~450 000 years.

Signature of the Paleo-Course Changes in the São Francisco River as Source of Genetic Structure in Neotropical Pithecopus nordestinus (Phyllomedusinae, Anura) Treefrog

Historical processes that have interrupted gene flow between distinct evolutionary lineages have played a fundamental role in the evolution of the enormous diversity of species found in the Neotropical region. Numerous studies have discussed the role of geographic barriers and Pleistocene forest refugia in the diversification of the region's biodiversity. In the present study, we investigated the relative contribution of these different factors to the evolutionary history of Pithecopus nordestinus, a Neotropical tree frog, which is amply distributed in the Brazilian Atlantic Forest and adjacent areas of the Caatinga biome. We used an extensive sample and multilocus DNA sequences to provide an overview of the intraspecific genetic diversity of P. nordestinus, characterize historical diversification patterns, and identify possible phylogenetic splits. We tested different scenarios of diversification based on Pleistocene Refugia and river barrier models using approximate Bayesian computation (ABC) and ecological niche modeling (ENM). The phylogenetic approach indicate the occurrence of processes of phylogeographic divergence in both time and space, related to historical shifts in the course of the São Francisco River during Plio-Pleistocene period, resulting in two principal, highly divergent clades. The ABC model provided strong statistical support for this scenario, confirming the hypothesis that the São Francisco River acted as an effective geographical barrier during vicariant events in the evolutionary history of P. nordestinus. We believe that the climatic changes that occurred during the Pleistocene also played a secondary role in the genetic signatures identified, reinforcing the divergence of populations isolated by physical barriers. These findings reinforce the conclusion that the two models of diversification (geographic barriers and refugia) are not mutually exclusive in the Neotropical domain but may interact extensively during the diversification of species on a regional scale.

Phylogeny of Neotropical Sicarius sand spiders suggests frequent transitions from deserts to dry forests despite antique, broad-scale niche conservatism

Phylogenetic niche conservatism (PNC) shapes the distribution of organisms by constraining lineages to particular climatic conditions. Conversely, if areas with similar climates are geographically isolated, diversification may also be limited by dispersal. Neotropical xeric habitats provide an ideal system to test the relative roles of climate and geography on diversification, as they occur in disjunct areas with similar biotas. Sicariinae sand spiders are intimately associated with these xeric environments, particularly seasonally dry tropical forests (SDTFs) and subtropical deserts/scrublands in Africa (Hexophthalma) and the Neotropics (Sicarius). We explore the role of PNC, geography and biome shifts in their evolution and timing of diversification. We estimated a time-calibrated, total-evidence phylogeny of Sicariinae, and used published distribution records to estimate climatic niche and biome occupancy. Topologies were used for estimating ancestral niches and biome shifts. We used variation partitioning methods to test the relative importance of climate and spatially autocorrelated factors in explaining the spatial variation in phylogenetic structure of Sicarius across the Neotropics. Neotropical Sicarius are ancient and split from their African sister-group around 90 (57-131) million years ago. Most speciation events took place in the Miocene. Sicariinae records can be separated in two groups corresponding to temperate/dry and tropical/seasonally dry climates. The ancestral climatic niche of Sicariinae are temperate/dry areas, with 2-3 shifts to tropical/seasonally dry areas in Sicarius. Similarly, ancestral biomes occupied by the group are temperate and dry (deserts, Mediterranean scrub, temperate grasslands), with 2-3 shifts to tropical, seasonally dry forests and grasslands. Most of the variation in phylogenetic structure is explained by long-distance dispersal limitation that is independent of the measured climatic conditions. Sicariinae have an ancient association to arid lands, suggesting that PNC prevented them from colonizing mesic habitats. However, niches are labile at a smaller scale, with several shifts from deserts to SDTFs. This suggests that PNC and long-distance dispersal limitation played major roles in confining lineages to isolated areas of SDTF/desert over evolutionary history, although shifts between xeric biomes occurred whenever geographical opportunities were presented.

Isolation with asymmetric gene flow during the nonsynchronous divergence of dry forest birds

Dry forest bird communities in South America are often fragmented by intervening mountains and rainforests, generating high local endemism. The historical assembly of dry forest communities likely results from dynamic processes linked to numerous population histories among codistributed species. Nevertheless, species may diversify in the same way through time if landscape and environmental features, or species ecologies, similarly structure populations. Here we tested whether six co-distributed taxon pairs that occur in the dry forests of the Tumbes and Marañón Valley of northwestern South America show concordant patterns and modes of diversification. We employed a genome reduction technique, double-digest restriction site-associated DNA sequencing, and obtained 4407-7186 genomewide SNPs. We estimated demographic history in each taxon pair and inferred that all pairs had the same best-fit demographic model: isolation with asymmetric gene flow from the Tumbes into the Marañón Valley, suggesting a common diversification mode. Overall, we also observed congruence in effective population size (N_e ) patterns where ancestral N_e were 2.9-11.0× larger than present-day Marañón Valley populations and 0.3-2.0× larger than Tumbesian populations. Present-day Marañón Valley N_e was smaller than Tumbes. In contrast, we found simultaneous population isolation due to a single event to be unlikely as taxon pairs diverged over an extended period of time (0.1-2.9 Ma) with multiple nonoverlapping divergence periods. Our results show that even when populations of codistributed species asynchronously diverge, the mode of their differentiation can remain conserved over millions of years. Divergence by allopatric isolation due to barrier formation does not explain the mode of differentiation between these two bird assemblages; rather, migration of individuals occurred before and after geographic isolation.

Three new species of Misionella from northern Brazil (Araneae, Haplogynae, Filistatidae)

Three new species of the genus Misionella are described from Brazil: Misionellacarajassp. n. and Misionellaaikewarasp. n. from caves in the states of Pará and Tocantins and Misionellapallidasp. n. from natural and synanthropic dry areas in the states of Piauí, Maranhão, Rio Grande do Norte and Bahia. These species seem to belong to a distinct group within the genus; the males have an elongate palpal tibia and bulb, a pair of characteristic and hirsute macrosetae in the second metatarsus and the females have internal genitalia with only one pair of spermathecae, with relatively short ducts, lacking the auxiliary receptacles. Their phylogenetic placement and geographic distribution are briefly discussed.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.