Climate change and carnivores: shifts in the distribution and effectiveness of protected areas in the Amazon.

Background: Carnivore mammals are animals vulnerable to human interference, such as climate change and deforestation. Their distribution and persistence are affected by such impacts, mainly in tropical regions such as the Amazon. Due to the importance of carnivores in the maintenance and functioning of the ecosystem, they are extremely important animals for conservation. We evaluated the impact of climate change on the geographic distribution of carnivores in the Amazon using Species Distribution Models (SDMs). Do we seek to answer the following questions: (1) What is the effect of climate change on the distribution of carnivores in the Amazon? (2) Will carnivore species lose or gain representation within the Protected Areas (PAs) of the Amazon in the future? Methods: We evaluated the distribution area of 16 species of carnivores mammals in the Amazon, based on two future climate scenarios (RCP 4.5 and RCP 8.5) for the year 2070. For the construction of the SDMs we used bioclimatic and vegetation cover variables (land type). Based on these models, we calculated the area loss and climate suitability of the species, as well as the effectiveness of the protected areas inserted in the Amazon. We estimated the effectiveness of PAs on the individual persistence of carnivores in the future, for this, we used the SDMs to perform the gap analysis. Finally, we analyze the effectiveness of PAs in protecting taxonomic richness in future scenarios. Results: The SDMs showed satisfactory predictive performance, with Jaccard values above 0.85 and AUC above 0.91 for all species. In the present and for the future climate scenarios, we observe a reduction of potencial distribution in both future scenarios (RCP4.5 and RCP8.5), where five species will be negatively affected by climate change in the RCP 4.5 future scenario and eight in the RCP 8.5 scenario. The remaining species stay stable in terms of total area. All species in the study showed a loss of climatic suitability. Some species lost almost all climatic suitability in the RCP 8.5 scenario. According to the GAP analysis, all species are protected within the PAs both in the current scenario and in both future climate scenarios. From the null models, we found that in all climate scenarios, the PAs are not efficient in protecting species richness.

Kill site database: A unified dataset on human-megafauna interactions across time and space.

The database presented in this data article is related to the paper "Megafauna kill sites in South America: a critical review" [1]. It includes a list of 134 publications on human-megafauna interaction, with 69 archaeological sites showing human-megafauna interaction. From these sites, 44 present a minimum human-megafauna association, from which up to 17 megafauna kill sites were classified, with up to 15 exploited extinct megafauna taxa. It also provides a list of current taxonomic classifications of extinct megafauna that humans have exploited according to empirical evidence presented in the related paper. The megafauna kill sites were classified based on five restrictive criteria according to Grayson and Meltzer's (2015, 2002), Borrero's (2009) and Mothé et al.'s (2020) protocol. The kill sites database reflects the empirical evidence on megafauna exploitation by humans available in scientific literature and is useful to understand the human-megafauna interactions in the late Quaternary. Finally, we also provide our online repository (www.killsitedatabase.com), an initiative to unify the evidence on megafauna kill sites (and their related data) worldwide, starting in South America.

Rugged relief and climate promote isolation and divergence between two neotropical cold-associated birds.

The role of historical factors in establishing patterns of diversity in tropical mountains is of interest to understand the buildup of megadiverse biotas. In these regions, the historical processes of range fragmentation and contraction followed by dispersal are thought to be mediated by the interplay between rugged relief (complex topography) and climate fluctuations and likely explain most of the dynamics of diversification in plants and animals. Although empirical studies addressing the interaction between climate and topography have provided invaluable insights into population divergence and speciation patterns in tropical montane organisms, a more detailed and robust test of such processes in an explicit spatio-temporal framework is still lacking. Consequently, our ability to gain insights into historical range shifts over time and the genomic footprint left by them is limited. Here, we used niche modeling and subgenomic population-level datasets to explore the evolution of two species of warbling finches (genus Microspingus) disjunctly distributed across the Montane Atlantic Forest, a Neotropical region with complex geological and environmental histories. Population structure inferences suggest a scenario of three genetically differentiated populations, which are congruent with both geography and phenotypic variation. Demographic simulations support asynchronous isolation of these populations as recently as ∼40,000 years ago, relatively stable population sizes over recent time, and past gene flow subsequent to divergence. Throughout the last 800,000 years, niche models predicted extensive expansion into lowland areas with increasing overlap of species distributions during glacial periods, with prominent retractions and isolation into higher altitudes during interglacials, which are in line with signs of introgression of currently isolated populations. These results support a dual role of cyclical climatic changes: population divergence and persistence in mountain tops during warm periods followed by periods of expansion and admixture in lower elevations during cold periods. Our results underscore the role of the interplay between landscape and climate as an important mechanism in the evolution of the Neotropical montane biota.

Climatic dynamics and topography control genetic variation in Atlantic Forest montane birds.

Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity.

A genome-wide scan shows evidence for local adaptation in a widespread keystone Neotropical forest tree.

The role of natural selection in shaping patterns of diversity is still poorly understood in the Neotropics. We carried out the first genome-wide population genomics study in a Neotropical tree, Handroanthus impetiginosus (Bignoniaceae), sampling 75,838 SNPs by sequence capture in 128 individuals across 13 populations. We found evidences for local adaptation using Bayesian correlations of allele frequency and environmental variables (32 loci in 27 genes) complemented by an analysis of selective sweeps and genetic hitchhiking events using SweepFinder2 (81 loci in 47 genes). Fifteen genes were identified by both approaches. By accounting for population genetic structure, we also found 14 loci with selection signal in a STRUCTURE-defined lineage comprising individuals from five populations, using Outflank. All approaches pinpointed highly diverse and structurally conserved genes affecting plant development and primary metabolic processes. Spatial interpolation forecasted differences in the expected allele frequencies at loci under selection over time, suggesting that H. impetiginosus may track its habitat during climate changes. However, local adaptation through natural selection may also take place, allowing species persistence due to niche evolution. A high genetic differentiation was seen among the H. impetiginosus populations, which, together with the limited power of the experiment, constrains the improved detection of other types of soft selective forces, such as background, balanced, and purifying selection. Small differences in allele frequency distribution among widespread populations and the low number of loci with detectable adaptive sweeps advocate for a polygenic model of adaptation involving a potentially large number of small genome-wide effects.

Unravelling the genetic differentiation among varieties of the Neotropical savanna tree Hancornia speciosa Gomes.

Background and Aims: Spatial distribution of species genetic diversity is often driven by geographical distance (isolation by distance) or environmental conditions (isolation by environment), especially under climate change scenarios such as Quaternary glaciations. Here, we used coalescent analyses coupled with ecological niche modelling (ENM), spatially explicit quantile regression analyses and the multiple matrix regression with randomization (MMRR) approach to unravel the patterns of genetic differentiation in the widely distributed Neotropical savanna tree, Hancornia speciosa (Apocynaceae). Due to its high morphological differentiation, the species was originally classified into six botanical varieties by Monachino, and has recently been recognized as only two varieties by Flora do Brasil 2020. Thus, H. speciosa is a good biological model for learning about evolution of phenotypic plasticity under genetic and ecological effects, and predicting their responses to changing environmental conditions. Methods: We sampled 28 populations (777 individuals) of Monachino's four varieties of H. speciosa and used seven microsatellite loci to genotype them. Key Results: Bayesian clustering showed five distinct genetic groups (K = 5) with high admixture among Monachino's varieties, mainly among populations in the central area of the species geographical range. Genetic differentiation among Monachino's varieties was lower than the genetic differentiation among populations within varieties, with higher within-population inbreeding. A high historical connectivity among populations of the central Cerrado shown by coalescent analyses may explain the high admixture among varieties. In addition, areas of higher climatic suitability also presented higher genetic diversity in such a way that the wide historical refugium across central Brazil might have promoted the long-term connectivity among populations. Yet, FST was significantly related to geographic distances, but not to environmental distances, and coalescent analyses and ENM predicted a demographical scenario of quasi-stability through time. Conclusions: Our findings show that demographical history and isolation by distance, but not isolation by environment, drove genetic differentiation of populations. Finally, the genetic clusters do not support the two recently recognized botanical varieties of H. speciosa, but partially support Monachino's classification at least for the four sampled varieties, similar to morphological variation.

Correction: Pollination Mode and Mating System Explain Patterns in Genetic Differentiation in Neotropical Plants.

[This corrects the article DOI: 10.1371/journal.pone.0158660.].

Two sides of a coin: Effects of climate change on the native and non-native distribution of Colossoma macropomum in South America.

Climate change and species invasions interact in nature, disrupting biological communities. Based on this knowledge, we simultaneously assessed the effects of climate change on the native distribution of the Amazonian fish Colossoma macropomum as well as on its invasiveness across river basins of South America, using ecological niche modeling. We used six niche models within the ensemble forecast context to predict the geographical distribution of C. macropomum for the present time, 2050 and 2080. Given that this species has been continuously introduced into non-native South American basins by fish farming activities, we added the locations of C. macropomum farms into the modeling process to obtain a more realistic scenario of its invasive potential. Based on modelling outputs we mapped climate refuge areas at different times. Our results showed that a plenty of climatically suitable areas for the occurrence of C. macropomum occurrence are located outside the original basins at the present time and that its invasive potential is greatly amplified by fish farms. Simulations of future geographic ranges revealed drastic range contraction in the native region, implying concerns not only with respect to the species conservation but also from a socio-economic perspective since the species is a cornerstone of artisanal and commercial fisheries in the Amazon. Although the invasive potential is projected to decrease in the face of climate change, climate refugia will concentrate in Paraná River, Southeast Atlantic and East Atlantic basins, putting intense, negative pressures on the native fish fauna these regions. Our findings show that short and long-term management actions are required for: i) the conservation of natural stocks of C. macropomum in the Amazon, and ii) protecting native fish fauna in the climate refuges of the invaded regions.

Climatic changes can drive the loss of genetic diversity in a Neotropical savanna tree species.

The high rates of future climatic changes, compared with the rates reported for past changes, may hamper species adaptation to new climates or the tracking of suitable conditions, resulting in significant loss of genetic diversity. Trees are dominant species in many biomes and because they are long-lived, they may not be able to cope with ongoing climatic changes. Here, we coupled ecological niche modelling (ENM) and genetic simulations to forecast the effects of climatic changes on the genetic diversity and the structure of genetic clusters. Genetic simulations were conditioned to climatic variables and restricted to plant dispersal and establishment. We used a Neotropical savanna tree as species model that shows a preference for hot and drier climates, but with low temperature seasonality. The ENM predicts a decreasing range size along the more severe future climatic scenario. Additionally, genetic diversity and allelic richness also decrease with range retraction and climatic genetic clusters are lost for both future scenarios, which will lead genetic variability to homogenize throughout the landscape. Besides, climatic genetic clusters will spatially reconfigure on the landscape following displacements of climatic conditions. Our findings indicate that climate change effects will challenge population adaptation to new environmental conditions because of the displacement of genetic ancestry clusters from their optimal conditions.

Demographic stability and high historical connectivity explain the diversity of a savanna tree species in the Quaternary.

Background and Aims: Cyclic glaciations were frequent throughout the Quaternary and this affected species distribution and population differentiation worldwide. The present study reconstructed the demographic history and dispersal routes of Eugenia dysenterica lineages and investigated the effects of Quaternary climate change on its spatial pattern of genetic diversity. Methods: A total of 333 individuals were sampled from 23 populations and analysed by sequencing four regions of the chloroplast DNA and the internal transcribed spacer of the nuclear DNA. The analyses were performed using a multi-model inference approach based on ecological niche modelling and statistical phylogeography. Key Results: Coalescent simulation showed that population stability through time is the most likely scenario. The palaeodistribution dynamics predicted by the ecological niche models revealed that the species was potentially distributed across a large area, extending over Central-Western Brazil through the last glaciation. The lineages of E. dysenterica dispersed from Central Brazil towards populations at the northern, western and south-eastern regions. A historical refugium through time may have favoured lineage dispersal and the maintenance of genetic diversity. Conclusions: The results suggest that the central region of the Cerrado biome is probably the centre of distribution of E. dysenterica and that the spatial pattern of its genetic diversity may be the outcome of population stability throughout the Quaternary. The lower genetic diversity in populations in the south-eastern Cerrado biome is probably due to local climatic instability during the Quaternary.

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.

Demographical history and palaeodistribution modelling show range shift towards Amazon Basin for a Neotropical tree species in the LGM.

BACKGROUND: We studied the phylogeography and demographical history of Tabebuia serratifolia (Bignoniaceae) to understand the disjunct geographical distribution of South American seasonally dry tropical forests (SDTFs). We specifically tested if the multiple and isolated patches of SDTFs are current climatic relicts of a widespread and continuously distributed dry forest during the last glacial maximum (LGM), the so called South American dry forest refugia hypothesis, using ecological niche modelling (ENM) and statistical phylogeography. We sampled 235 individuals of T. serratifolia in 17 populations in Brazil and analysed the polymorphisms at three intergenic chloroplast regions and ITS nuclear ribosomal DNA. RESULTS: Coalescent analyses showed a demographical expansion at the last c. 130 ka (thousand years before present). Simulations and ENM also showed that the current spatial pattern of genetic diversity is most likely due to a scenario of range expansion and range shift towards the Amazon Basin during the colder and arid climatic conditions associated with the LGM, matching the expected for the South American dry forest refugia hypothesis, although contrasting to the Pleistocene Arc hypothesis. Populations in more stable areas or with higher suitability through time showed higher genetic diversity. Postglacial range shift towards the Southeast and Atlantic coast may have led to spatial genome assortment due to leading edge colonization as the species tracks suitable environments, leading to lower genetic diversity in populations at higher distance from the distribution centroid at 21 ka. CONCLUSION: Haplotype sharing or common ancestry among populations from Caatinga in Northeast Brazil, Atlantic Forest in Southeast and Cerrado biome and ENM evince the past connection among these biomes.

Pollination Mode and Mating System Explain Patterns in Genetic Differentiation in Neotropical Plants.

We studied genetic diversity and differentiation patterns in Neotropical plants to address effects of life history traits (LHT) and ecological attributes based on an exhaustive literature survey. We used generalized linear mixed models (GLMMs) to test the effects as fixed and random factors of growth form, pollination and dispersal modes, mating and breeding systems, geographical range and habitat on patterns of genetic diversity (HS, HeS, π and h), inbreeding coefficient (FIS), allelic richness (AR) and differentiation among populations (FST) for both nuclear and chloroplast genomes. In addition, we used phylogenetic generalized least squares (pGLS) to account for phylogenetic independence on predictor variables and verify the robustness of the results from significant GLMMs. In general, GLMM revealed more significant relationships among LHTs and genetic patterns than pGLS. After accounting for phylogenetic independence (i.e., using pGLS), FST for nuclear microsatellites was significantly related to pollination mode, mating system and habitat. Plants specifically with outcrossing mating system had lower FST. Moreover, AR was significantly related to pollination mode and geographical range and HeS for nuclear dominant markers was significantly related to habitat. Our findings showed that different results might be retrieved when phylogenetic non-independence is taken into account and that LHTs and ecological attributes affect substantially the genetic pattern in Neotropical plants, hence may drive key evolutionary processes in plants.

Coalescent Simulation and Paleodistribution Modeling for Tabebuia rosealba Do Not Support South American Dry Forest Refugia Hypothesis.

Studies based on contemporary plant occurrences and pollen fossil records have proposed that the current disjunct distribution of seasonally dry tropical forests (SDTFs) across South America is the result of fragmentation of a formerly widespread and continuously distributed dry forest during the arid climatic conditions associated with the Last Glacial Maximum (LGM), which is known as the modern-day dry forest refugia hypothesis. We studied the demographic history of Tabebuia rosealba (Bignoniaceae) to understand the disjunct geographic distribution of South American SDTFs based on statistical phylogeography and ecological niche modeling (ENM). We specifically tested the dry forest refugia hypothesis; i.e., if the multiple and isolated patches of SDTFs are current climatic relicts of a widespread and continuously distributed dry forest during the LGM. We sampled 235 individuals across 18 populations in Central Brazil and analyzed the polymorphisms at chloroplast (trnS-trnG, psbA-trnH and ycf6-trnC intergenic spacers) and nuclear (ITS nrDNA) genomes. We performed coalescence simulations of alternative hypotheses under demographic expectations from two a priori biogeographic hypotheses (1. the Pleistocene Arc hypothesis and, 2. a range shift to Amazon Basin) and other two demographic expectances predicted by ENMs (3. expansion throughout the Neotropical South America, including Amazon Basin, and 4. retraction during the LGM). Phylogenetic analyses based on median-joining network showed haplotype sharing among populations with evidence of incomplete lineage sorting. Coalescent analyses showed smaller effective population sizes for T. roseoalba during the LGM compared to the present-day. Simulations and ENM also showed that its current spatial pattern of genetic diversity is most likely due to a scenario of range retraction during the LGM instead of the fragmentation from a once extensive and largely contiguous SDTF across South America, not supporting the South American dry forest refugia hypothesis.

Relaxed random walk model coupled with ecological niche modeling unravel the dispersal dynamics of a Neotropical savanna tree species in the deeper Quaternary.

Understanding the dispersal routes of Neotropical savanna tree species is an essential step to unravel the effects of past climate change on genetic patterns, species distribution and population demography. Here we reconstruct the demographic history and dispersal dynamics of the Neotropical savanna tree species Tabebuia aurea to understand the effects of Quaternary climate change on its current spatial patterns of genetic diversity. We sampled 285 individuals from 21 populations throughout Brazilian savannas and sequenced all individuals for three chloroplast intergenic spacers and ITS nrDNA. We analyzed data using a multi-model inference framework by coupling the relaxed random walk model (RRW), ecological niche modeling (ENM) and statistical phylogeography. The most recent common ancestor of T. aurea lineages dated from ~4.0 ± 2.5 Ma. T. aurea lineages cyclically dispersed from the West toward the Central-West Brazil, and from the Southeast toward the East and Northeast Brazil, following the paleodistribution dynamics shown by the ENMs through the last glacial cycle. A historical refugium through time may have allowed dispersal of lineages among populations of Central Brazil, overlapping with population expansion during interglacial periods and the diversification of new lineages. Range and population expansion through the Quaternary were, respectively, the most frequent prediction from ENMs and the most likely demographic scenario from coalescent simulations. Consistent phylogeographic patterns among multiple modeling inferences indicate a promising approach, allowing us to understand how cyclical climate changes through the Quaternary drove complex population dynamics and the current patterns of species distribution and genetic diversity.

A Short Guide to the Climatic Variables of the Last Glacial Maximum for Biogeographers.

Ecological niche models are widely used for mapping the distribution of species during the last glacial maximum (LGM). Although the selection of the variables and General Circulation Models (GCMs) used for constructing those maps determine the model predictions, we still lack a discussion about which variables and which GCM should be included in the analysis and why. Here, we analyzed the climatic predictions for the LGM of 9 different GCMs in order to help biogeographers to select their GCMs and climatic layers for mapping the species ranges in the LGM. We 1) map the discrepancies between the climatic predictions of the nine GCMs available for the LGM, 2) analyze the similarities and differences between the GCMs and group them to help researchers choose the appropriate GCMs for calibrating and projecting their ecological niche models (ENM) during the LGM, and 3) quantify the agreement of the predictions for each bioclimatic variable to help researchers avoid the environmental variables with a poor consensus between models. Our results indicate that, in absolute values, GCMs have a strong disagreement in their temperature predictions for temperate areas, while the uncertainties for the precipitation variables are in the tropics. In spite of the discrepancies between model predictions, temperature variables (BIO1-BIO11) are highly correlated between models. Precipitation variables (BIO12-BIO19) show no correlation between models, and specifically, BIO14 (precipitation of the driest month) and BIO15 (Precipitation Seasonality (Coefficient of Variation)) show the highest level of discrepancy between GCMs. Following our results, we strongly recommend the use of different GCMs for constructing or projecting ENMs, particularly when predicting the distribution of species that inhabit the tropics and the temperate areas of the Northern and Southern Hemispheres, because climatic predictions for those areas vary greatly among GCMs. We also recommend the exclusion of BIO14 and BIO15 from ENMs because those variables show a high level of discrepancy between GCMs. Thus, by excluding them, we decrease the level of uncertainty of our predictions. All the climatic layers produced for this paper are freely available in http://ecoclimate.org/.

Multi-model inference in comparative phylogeography: an integrative approach based on multiple lines of evidence.

Comparative phylogeography has its roots in classical biogeography and, historically, relies on a pattern-based approach. Here, we present a model-based framework for comparative phylogeography. Our framework was initially developed for statistical phylogeography based on a multi-model inference approach, by coupling ecological niche modeling, coalescent simulation and direct spatio-temporal reconstruction of lineage diffusion using a relaxed random walk model. This multi-model inference framework is particularly useful to investigate the complex dynamics and current patterns in genetic diversity in response to processes operating on multiple taxonomic levels in comparative phylogeography. In addition, because of the lack, or incompleteness of fossil record, the understanding of the role of biogeographical events (vicariance and dispersal routes) in most regions worldwide is barely known. Thus, we believe that the expansion of that framework for multiple species under a comparative approach may give clues on genetic legacies in response to Quaternary climate changes and other biogeographical processes.

Recovering species demographic history from multi-model inference: the case of a Neotropical savanna tree species.

BACKGROUND: Glaciations were recurrent throughout the Quaternary and potentially shaped species genetic structure worldwide by affecting population dynamics. Here, we implemented a multi-model inference approach to recover the distribution dynamics and demographic history of a Neotropical savanna tree, Tabebuia aurea (Bignoniaceae). Exploring different algorithms and paleoclimatic simulations, we used ecological niche modelling to generate alternative hypotheses of potential demographic changes through the last glacial cycle and estimated genetic parameters using coalescent modelling. RESULTS: Comparing predictions from demographic hypotheses with genetic parameters of modern populations, our findings revealed a likely scenario of population decline, with spatial displacement towards Northeast Brazil from the last glacial maximum to the mid-Holocene. Subsequently, populations expanded in response to the return of the climatically suitable conditions in Central-West Brazil. Nevertheless, a wide historical refugium across Central Brazil likely maintained large populations connected throughout time. The expected genetic signatures from such predicted distribution dynamics are also corroborated by spatial genetic structure observed in modern populations. CONCLUSION: By exploring uncertainties inherent in multiple working hypotheses, we have shown that multi-model inference is a fruitful and efficient approach to recover the nature, timing and geographical context of the Tabebuia aurea population dynamic in response to the Quaternary climate changes.

Ecological niche and phylogeography elucidate complex biogeographic patterns in Loxosceles rufescens (Araneae, Sicariidae) in the Mediterranean Basin.

BACKGROUND: Understanding the evolutionary history of morphologically cryptic species complexes is difficult, and made even more challenging when geographic distributions have been modified by human-mediated dispersal. This situation is common in the Mediterranean Basin where, aside from the environmental heterogeneity of the region, protracted human presence has obscured the biogeographic processes that shaped current diversity. Loxosceles rufescens (Araneae, Sicariidae) is an ideal example: native to the Mediterranean, the species has dispersed worldwide via cohabitation with humans. A previous study revealed considerable molecular diversity, suggesting cryptic species, but relationships among lineages did not correspond to geographic location. RESULTS: Delimitation analyses on cytochrome c oxidase subunit I identified 11 different evolutionary lineages, presenting two contrasting phylogeographic patterns: (1) lineages with well-structured populations in Morocco and Iberia, and (2) lineages lacking geographic structure across the Mediterranean Basin. Dating analyses placed main diversification events in the Pleistocene, and multiple Pleistocene refugia, identified using ecological niche modeling (ENM), are compatible with allopatric differentiation of lineages. Human-mediated transportation appears to have complicated the current biogeography of this medically important and synanthropic spider. CONCLUSIONS: We integrated ecological niche models with phylogeographic analyses to elucidate the evolutionary history of L. rufescens in the Mediterranean Basin, with emphasis on the origins of mtDNA diversity. We found support for the hypothesis that northern Africa was the center of origin for L. rufescens, and that current genetic diversity originated in allopatry, likely promoted by successive glaciations during the Pleistocene. We corroborated the scenario of multiple refugia within the Mediterranean, principally in northern Africa, combining results from eight atmosphere-ocean general circulation models (AOGCMs) with two different refugium-delimitation methodologies. ENM results were useful for providing general views of putative refugia, with fine-scale details depending on the level of stringency applied for agreement among models.

A coupled phylogeographical and species distribution modelling approach recovers the demographical history of a Neotropical seasonally dry forest tree species.

We investigated here the demographical history of Tabebuia impetiginosa (Bignoniaceae) to understand the dynamics of the disjunct geographical distribution of South American seasonally dry forests (SDFs), based on coupling an ensemble approach encompassing hindcasting species distribution modelling and statistical phylogeographical analysis. We sampled 17 populations (280 individuals) in central Brazil and analysed the polymorphisms at chloroplast (trnS-trnG, psbA-trnH, and ycf6-trnC intergenic spacers) and nuclear (ITS nrDNA) genomes. Phylogenetic analyses based on median-joining network showed no haplotype sharing among population but strong evidence of incomplete lineage sorting. Coalescent analyses showed historical constant populations size, negligible gene flow among populations, and an ancient time to most recent common ancestor dated from ~4.7 ± 1.1 Myr BP. Most divergences dated from the Lower Pleistocene, and no signal of important population size reduction was found in coalescent tree and tests of demographical expansion. Demographical scenarios were built based on past geographical range dynamic models, using two a priori biogeographical hypotheses ('Pleistocene Arc' and 'Amazonian SDF expansion') and on two additional hypotheses suggested by the palaeodistribution modelling built with several algorithms for distribution modelling and palaeoclimatic data. The simulation of these demographical scenarios showed that the pattern of diversity found so far for T. impetiginosa is in consonance with a palaeodistribution expansion during the last glacial maximum (LGM, 21 kyr BP), strongly suggesting that the current disjunct distribution of T. impetiginosa in SDFs may represent a climatic relict of a once more wide distribution.