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, ecological niche modeling (ENM) and statistical phylogeography. The most recent common ancestor of T. aurea lineages dated from ~4.0 ± 2.5 Ma. Tabebuia aurea lineages cyclically dispersed from the West towards the Central-West Brazil, and from the Southeast towards 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.

The use of ecological niche models (ENM) to generate potential geographic distributions of species has rapidly increased in ecology, conservation and evolutionary biology. Many methods are available and the most used are Maximum Entropy Method (MAXENT) and the Genetic Algorithm for Rule Set Production (GARP). Recent studies have shown that MAXENT perform better than GARP. Here we used the statistics methods of ROC - AUC (area under the Receiver Operating Characteristics curve) and bootstrap to evaluate the performance of GARP and MAXENT in generate potential distribution models for 39 species of New World coral snakes. We found that values of AUC for GARP ranged from 0.923 to 0.999, whereas those for MAXENT ranged from 0.877 to 0.999. On the whole, the differences in AUC were very small, but for 10 species GARP outperformed MAXENT. Means and standard deviations for 100 bootstrapped samples with sample sizes ranging from 3 to 30 species did not show any trends towards deviations from a zero difference in AUC values of GARP minus AUC values of MAXENT. Ours results suggest that further studies are still necessary to establish under which circumstances the statistical performance of the methods vary. However, it is also important to consider the possibility that this empirical inductive reasoning may fail in the end, because we almost certainly could not establish all potential scenarios generating variation in the relative performance of models.

• Premise of the study: Wild edible plants (WEPs) have an important cultural and economic role in human population worldwide. Human impacts are quickly converting natural habitats in agricultural, cattle ranch, and urbanized lands, putting native species on peril of risk of extinction, including some WEPs. Moreover, global climate changes also can pose another threat to species persistency. Here, we established conservation priorities for the Cerrado, a neotropical region in South America with high levels of plant endemism and vulnerability, aiming to assure long‐term persistency of 16 most important WEPs. We evaluated these conservation priorities using a conservation biogeography framework using ecological patterns and process at a biogeographical scale to deal with species conservation features.• Methods: We built ecological niche models for 16 WEPs from Cerrado in the neotropics using climate models for preindustrial, past (Last Glacial Maximum) and future (year 2080) time periods to establish climatically stable areas through time, finding refugias for these WEPs. We used a spatial prioritization algorithm based on the spatial pattern of irreplaceability across the neotropics, aiming to ensure the persistence of at least 25% of range size in climatically stable areas for each WEP, using agricultural models as constraints.• Key results: The Southeast Cerrado was the most biotically stable and irreplaceable region for the WEPs compared with other areas across the neotropics.• Conclusions: Our findings strongly suggest that the Southeast Cerrado should be considered a conservation priority, with new protected areas to be sustainably managed and restored, to guarantee the supply of cultural and ecosystem services provided from the Cerrado's WEPs.