The geography of climate and the global patterns of species diversity.

Climate's effect on global biodiversity is typically viewed through the lens of temperature, humidity and resulting ecosystem productivity1-6. However, it is not known whether biodiversity depends solely on these climate conditions, or whether the size and fragmentation of these climates are also crucial. Here we shift the common perspective in global biodiversity studies, transitioning from geographic space to a climate-defined multidimensional space. Our findings suggest that larger and more isolated climate conditions tend to harbour higher diversity and species turnover among terrestrial tetrapods, encompassing more than 30,000 species. By considering both the characteristics of climate itself and its geographic attributes, we can explain almost 90% of the variation in global species richness. Half of the explanatory power (45%) may be attributed either to climate itself or to the geography of climate, suggesting a nuanced interplay between them. Our work evolves the conventional idea that larger climate regions, such as the tropics, host more species primarily because of their size7,8. Instead, we underscore the integral roles of both the geographic extent and degree of isolation of climates. This refined understanding presents a more intricate picture of biodiversity distribution, which can guide our approach to biodiversity conservation in an ever-changing world.

Landscape dynamics and diversification of the megadiverse South American freshwater fish fauna.

Landscape dynamics are widely thought to govern the tempo and mode of continental radiations, yet the effects of river network rearrangements on dispersal and lineage diversification remain poorly understood. We integrated an unprecedented occurrence dataset of 4,967 species with a newly compiled, time-calibrated phylogeny of South American freshwater fishes-the most species-rich continental vertebrate fauna on Earth-to track the evolutionary processes associated with hydrogeographic events over 100 Ma. Net lineage diversification was heterogeneous through time, across space, and among clades. Five abrupt shifts in net diversification rates occurred during the Paleogene and Miocene (between 30 and 7 Ma) in association with major landscape evolution events. Net diversification accelerated from the Miocene to the Recent (c. 20 to 0 Ma), with Western Amazonia having the highest rates of in situ diversification, which led to it being an important source of species dispersing to other regions. All regional biotic interchanges were associated with documented hydrogeographic events and the formation of biogeographic corridors, including the Early Miocene (c. 23 to 16 Ma) uplift of the Serra do Mar and Serra da Mantiqueira and the Late Miocene (c. 10 Ma) uplift of the Northern Andes and associated formation of the modern transcontinental Amazon River. The combination of high diversification rates and extensive biotic interchange associated with Western Amazonia yielded its extraordinary contemporary richness and phylogenetic endemism. Our results support the hypothesis that landscape dynamics, which shaped the history of drainage basin connections, strongly affected the assembly and diversification of basin-wide fish faunas.

gen3sis: A general engine for eco-evolutionary simulations of the processes that shape Earth's biodiversity.

Understanding the origins of biodiversity has been an aspiration since the days of early naturalists. The immense complexity of ecological, evolutionary, and spatial processes, however, has made this goal elusive to this day. Computer models serve progress in many scientific fields, but in the fields of macroecology and macroevolution, eco-evolutionary models are comparatively less developed. We present a general, spatially explicit, eco-evolutionary engine with a modular implementation that enables the modeling of multiple macroecological and macroevolutionary processes and feedbacks across representative spatiotemporally dynamic landscapes. Modeled processes can include species' abiotic tolerances, biotic interactions, dispersal, speciation, and evolution of ecological traits. Commonly observed biodiversity patterns, such as α, ß, and γ diversity, species ranges, ecological traits, and phylogenies, emerge as simulations proceed. As an illustration, we examine alternative hypotheses expected to have shaped the latitudinal diversity gradient (LDG) during the Earth's Cenozoic era. Our exploratory simulations simultaneously produce multiple realistic biodiversity patterns, such as the LDG, current species richness, and range size frequencies, as well as phylogenetic metrics. The model engine is open source and available as an R package, enabling future exploration of various landscapes and biological processes, while outputs can be linked with a variety of empirical biodiversity patterns. This work represents a key toward a numeric, interdisciplinary, and mechanistic understanding of the physical and biological processes that shape Earth's biodiversity.

Repeated Evolution of Unorthodox Feeding Styles Drives a Negative Correlation between Foot Size and Bill Length in Hummingbirds.

AbstractDifferences among hummingbird species in bill length and shape have rightly been viewed as adaptive in relation to the morphology of the flowers they visit for nectar. In this study we examine functional variation in a behaviorally related but neglected feature: hummingbird feet. We gathered records of hummingbirds clinging by their feet to feed legitimately as pollinators or illegitimately as nectar robbers-"unorthodox" feeding behaviors. We measured key features of bills and feet for 220 species of hummingbirds and compared the 66 known "clinger" species (covering virtually the entire scope of hummingbird body size) with the 144 presumed "non-clinger" species. Once the effects of phylogenetic signal, body size, and elevation above sea level are accounted for statistically, hummingbirds display a surprising but functionally interpretable negative correlation. Clingers with short bills and long hallux (hind-toe) claws have evolved-independently-more than 20 times and in every major clade. Their biomechanically enhanced feet allow them to save energy by clinging to feed legitimately on short-corolla flowers and by stealing nectar from long-corolla flowers. In contrast, long-billed species have shorter hallux claws, as plant species with long-corolla flowers enforce hovering to feed, simply by the way they present their flowers.

Salivary annexin A1: A candidate biomarker for periodontitis.

AIM: To compare the salivary proteomic profile of periodontitis-affected (PA) parents and their offspring to periodontally healthy (PH) dyads in the pursuit of possible biomarkers for early diagnosis of this disease. MATERIALS AND METHODS: Unstimulated saliva samples collected from 17 pairs of PA or PH individuals and their children were submitted to mass spectrometric analyses followed by proteomic analyses. Primary PA fibroblasts were triggered towards having an inflammatory response, and an immunoenzymatic assay of its supernatant was performed to validate the obtained data. RESULTS: ANXA1, KRT4, GSTP1, HPX, A2M and KRT13 were lower in PA parents and their children, and IGHG1, CSTB, KRT9, SMR3B, IGHG4 and SERPINA1 were higher. ANXA1 presented the highest fold change, 7.1 times less produced in children of PA parents, and was selected as a potential biomarker for periodontitis. The in vitro assay also showed lower ANXA1 production by cells of PA patients. CONCLUSION: Before any clinical sign of periodontal loss, descendants of PA patients have an altered proteomic profile compared to PH individuals, presenting a lower abundance of ANXA1. This protein is suggested as a potential biomarker for periodontitis.

Smoking negatively impacts the clinical, microbiological, and immunological treatment response of young adults with Grade C periodontitis.

OBJECTIVE: This study aimed to investigate the influence of smoking on clinical, microbiological and immunological parameters in young adult with stage III-IV Grade C periodontitis after full-mouth ultrasonic debridement (FMUD) associated with Amoxicillin and Metronidazole (AMX + MTZ), comparing smokers (PerioC-Y-Smk) with non-smokers (PerioC-Y-NSmk). MATERIALS AND METHODS: Fifteen PerioC-Y-NSmk and 14 PerioC-Y-Smk patients underwent FMUD associated with AMX + MTZ for 10 days. All parameters were collected at baseline and 3 and 6 months after treatment. Plaque index (PI), bleeding on probing (BoP), probing depth (PD), clinical attachment level (CAL)- the primary variable-, and gingival recession (GR) were clinically assessed. The impact of PI on CAL change at 6-month was verified by a regression analysis. Samples of the subgingival biofilm was collected for detection of levels of Aggregatibacter actinomycetemcomitans (A. actinomycetemcomitans), Porphyromonas gingivalis (P.gingivalis), Tannerella forsythia (T. forsythia), and Fusobacterium nucleatum ssp (F. nucleatum), and were analyzed by real-time qPCR; gingival crevicular fluid was collected for detection of levels of interleukin (IL)-1ß, IL-4, IL-6, IL-10, tumor necrosis factor (TNF)-α, and interferon (IFN)-γ, which were analyzed using an enzyme immunoassay. RESULTS: PerioC-Y-Smk had significantly higher PI, BOP, and GR at baseline compared to non-smokers (p < .05). PerioC-Y-Smk presented higher PD, CAL, and GR at 3 and 6 months (p < .05) compared with PerioC-Y-NSmk in the same periods; PI negatively affected CAL gain in PerioC-Y-NSmk at 6-month follow-up (p = .052) and did not impact on clinical response in PerioC-Y-Smk (p = .882). Lower levels of IFN-γ, IL1-ß, and IL-4 were observed at 3 months in the PerioC-Y-NSmk (p < .05) compared with PerioC-Y-Smk. Lower proportions of P. gingivalis were observed in PerioC-Y-NSmk at baseline and at 3 months (p < .05) and lower proportions of F. nucleatum were observed at 6 months, in the PerioC-Y-NSmk (p < .05). CONCLUSIONS: PerioC-Y-Smk presents an unfavorable clinical, microbiological, and immunological response after 3 and 6 months after FMUD associated with AMX + MTZ. CLINICAL RELEVANCE: Smoking worsens periodontal condition of young treated adults presenting stage III/IV Grade C periodontitis.

Area, isolation and climate explain the diversity of mammals on islands worldwide.

Insular biodiversity is expected to be regulated differently than continental biota, but their determinants remain to be quantified at a global scale. We evaluated the importance of physical, environmental and historical factors on mammal richness and endemism across 5592 islands worldwide. We fitted generalized linear and mixed models to accommodate variation among biogeographic realms and performed analyses separately for bats and non-volants. Richness on islands ranged from one to 234 species, with up to 177 single island endemics. Diversity patterns were most consistently influenced by the islands' physical characteristics. Area positively affected mammal diversity, in particular the number of non-volant endemics. Island isolation, both current and past, was associated with lower richness but greater endemism. Flight capacity modified the relative importance of past versus current isolation, with bats responding more strongly to current and non-volant mammals to past isolation. Biodiversity relationships with environmental factors were idiosyncratic, with a tendency for greater effects sizes with endemism than richness. The historical climatic change was positively associated with endemism. In line with theory, we found that area and isolation were among the strongest drivers of mammalian biodiversity. Our results support the importance of past conditions on current patterns, particularly of non-volant species.

Spatial variation in direct and indirect effects of climate and productivity on species richness of terrestrial tetrapods.

AIM: We aimed to dissect the spatial variation of the direct and indirect effects of climate and productivity on global species richness of terrestrial tetrapods. LOCATION: Global. TIME PERIOD: Present. MAJOR TAXA STUDIED: Terrestrial tetrapods. METHODS: We used a geographically weighted path analysis to estimate and map the direct and indirect effects of temperature, precipitation and primary productivity on species richness of terrestrial tetrapods across the globe. RESULTS: We found that all relationships shift in magnitude, and even in direction, among taxonomic groups, geographical regions and connecting paths. Direct effects of temperature and precipitation are generally stronger than both indirect effects mediated by productivity and direct effects of productivity. MAIN CONCLUSIONS: Richness gradients seem to be driven primarily by effects of climate on organismal physiological limits and metabolic rates rather than by the amount of productive energy. Reptiles have the most distinct relationships across tetrapods, with a clear latitudinal pattern in the importance of temperature versus water.

Enamel matrix derivative effects on palatal mucosa wound healing: Randomized clinical trial.

BACKGROUND: This study was conducted to evaluate the clinical, immunologic, and patient-centered outcomes of enamel matrix protein derivative (EMD) on excisional wounds in palatal mucosa. MATERIALS: Forty-four patients in need of ridge preservation were randomly allocated into two groups: control group (n = 22): open palatal wound after free gingival graft (FGG) harvest and EMD group (n = 22): open palatal wound after FGG harvest that received 0.3 ml of EMD. Clinical and patient-centered parameters were analyzed for 3 months post-treatment. Wound fluid levels of inflammatory markers were assessed 3 and 7 days postoperatively. RESULTS: No significant inter-group difference was observed in remaining wound area and re-epithelialization. EMD and control groups achieved wound closure and re-epithelialization 30 days postoperatively (p < .001), without inter-group differences. Similarly, number of analgesics and Oral Health Impact Profile scores did not present significant inter-group differences (p > .05). EMD appeared to selectively modulate wound fluid levels of monocyte chemoattractant protein-1, macrophage inflammatory protein-1α, matrix metallopeptidase 9, and tissue inhibitor of metalloproteinases-2. CONCLUSION: Within the limits of the present study, it can be concluded that EMD application to excisional palatal wounds using the investigated protocol does not provide clinical healing benefits, despite an apparent modulation of selected inflammatory markers.

Subgingival endotoxin and lipoteichoic acid modulate cytokine production in diabetic subjects: A Case-control Study.

BACKGROUND: Periodontal disease and diabetes mellitus (DM) are highly prevalent and interrelated diseases, resulting in altered host response microbiota. Thus, this study aimed to evaluate the impact of DM on local levels of lipopolysaccharide (LPS) and lipoteichoic acid (LTA) and their relationship with cytokines and matrix metalloproteinases' (MMPs) profile. METHODS: This case-control study included diabetic (n = 15) and non-diabetic (n = 15) subjects presenting Stage 3-4, Grade C, Periodontitis. Gingival crevicular fluid (GCF) was collected, and LPS and LTA levels were analyzed by enzyme-linked immunosorbent assay (ELISA), while IFN-γ, IL-10, IL-17, IL-1ß, IL-4, MMP-2, and MMP-9 were measured by LUMINEX/MAGpix. Mann-Whitney and Spearman's correlation tests were used to compared and to correlate variables (p < 0.05). RESULTS: Higher levels of LTA, LPS, IL-10, IL-1ß, and MMP-2 (p < 0.05) and lower levels of IL-17 were found in the DM group (p < 0.05). Non-diabetic subjects presented higher LPS, IFN-γ, IL-17, and MMP-2 levels and lower IL-10 concentration (p < 0.05). No significant correlation was seen between LPS and cytokine profile in non-diabetic. Local levels of LTA were positively correlated with IL-17 and MMP-2 and negatively with IL-10. CONCLUSION: LTA and LPS drove the inflammatory profile through the modulation of cytokines and MMPs in a different manner in DM and non-diabetic subjects.

Comparative metal accumulation and toxicogenetic damage induction in three neotropical fish species with distinct foraging habits and feeding preferences.

Environmental pollutants affect fish species differently because their routes of exposure make them more critical to a specific group regarding foraging habits or food preferences. However, the association between local problems and environmental conditions makes a particular species more suitable for use in monitoring programs. Thus, this study compared muscle accumulation of metals (Al, Ba, Cr, Cu, Fe, Mn, Ti and Zn) and toxicogenetic damage to three neotropical species from the Itapemirim River Basin (Brazil) with distinct foraging habits and feeding preferences: Geophagus brasiliensis (benthopelagic - omnivorous); Harttia sp. (benthic - detritivorous); and Leporinus copellandii (pelagic - omnivorous). There was seasonal influence on muscle metal concentrations: Ba and Mn concentrations in G. brasiliensis, Cr in Harttia sp. and Zn, Al and Fe in L. copellandii increased during the rainy season. G. brasiliensis was the most sensitive species, since it had a statistically higher prevalence of erythrocytic nuclear abnormalities, which was influenced by seasonality. This result might be an effect of its benthopelagic habits, which favors contact with a larger group of contaminants due to its interaction with both water column and sediment. The present study showed the differential sensitivity of fish species and that the combination of chemical analysis of pollutants with evaluation of toxicogenetic responses helps to choose the best species for field studies.

Drivers of geographical patterns of North American language diversity.

Although many hypotheses have been proposed to explain why humans speak so many languages and why languages are unevenly distributed across the globe, the factors that shape geographical patterns of cultural and linguistic diversity remain poorly understood. Prior research has tended to focus on identifying universal predictors of language diversity, without accounting for how local factors and multiple predictors interact. Here, we use a unique combination of path analysis, mechanistic simulation modelling, and geographically weighted regression to investigate the broadly described, but poorly understood, spatial pattern of language diversity in North America. We show that the ecological drivers of language diversity are not universal or entirely direct. The strongest associations imply a role for previously developed hypothesized drivers such as population density, resource diversity, and carrying capacity with group size limits. The predictive power of this web of factors varies over space from regions where our model predicts approximately 86% of the variation in diversity, to areas where less than 40% is explained.

Quantitative genetics of body size evolution on islands: an individual-based simulation approach.

According to the island rule, small-bodied vertebrates will tend to evolve larger body size on islands, whereas the opposite happens to large-bodied species. This controversial pattern has been studied at the macroecological and biogeographical scales, but new developments in quantitative evolutionary genetics now allow studying the island rule from a mechanistic perspective. Here, we develop a simulation approach based on an individual-based model to model body size change on islands as a progressive adaptation to a moving optimum, determined by density-dependent population dynamics. We applied the model to evaluate body size differentiation in the pigmy extinct hominin Homo floresiensis, showing that dwarfing may have occurred in only about 360 generations (95% CI ranging from 150 to 675 generations). This result agrees with reports suggesting rapid dwarfing of large mammals on islands, as well as with the recent discovery that small-sized hominins lived in Flores as early as 700 kyr ago. Our simulations illustrate the power of analysing ecological and evolutionary patterns from an explicit quantitative genetics perspective.

Genetic Population Structure and Allele Surfing During Range Expansion in Dynamic Habitats.

Expanding populations may loss genetic diversity because sequential founder events throughout a wave of demographic expansion may cause "allele surfing", as the alleles of founder individuals may propagate rapidly through space. The spatial components of allele surfing have been studied by geneticists, but have never been investigate on dynamic and shifting habitats. Here we used an individual-based-model (IBM) to study how interactions between different habitat restoration scenarios and biological characteristics (dispersal capacity) affect the spatial patterns of the genetic structure of a population during demographic expansion. We found that both habitat dynamics and dispersal capacity, as well as their interaction, were the drivers of emergent pattern of genetic diversity and allele surfing. Specifically, allele surfing is more common when a species with low dispersal capacity colonizes a large geographic area with slow restoration (low carrying capacity). Despite this, we showed that allele surfing can be reduced, or even avoided, by dispersal management through suitable habitat restoration. Thus, investigating how colonization generates a spatial variation in genetic diversity, and which parameters control the emergent genetic pattern, are essential steps to planning assisted gene flow, which is fundamental for an effective planning of habitat restoration.

Neutral Community Dynamics and the Evolution of Species Interactions.

A contemporary goal in ecology is to determine the ecological and evolutionary processes that generate recurring structural patterns in mutualistic networks. One of the great challenges is testing the capacity of neutral processes to replicate observed patterns in ecological networks, since the original formulation of the neutral theory lacks trophic interactions. Here, we develop a stochastic-simulation neutral model adding trophic interactions to the neutral theory of biodiversity. Without invoking ecological differences among individuals of different species, and assuming that ecological interactions emerge randomly, we demonstrate that a spatially explicit multitrophic neutral model is able to capture the recurrent structural patterns of mutualistic networks (i.e., degree distribution, connectance, nestedness, and phylogenetic signal of species interactions). Nonrandom species distribution, caused by probabilistic events of migration and speciation, create nonrandom network patterns. These findings have broad implications for the interpretation of niche-based processes as drivers of ecological networks, as well as for the integration of network structures with demographic stochasticity.

SUNPLIN: simulation with uncertainty for phylogenetic investigations.

BACKGROUND: Phylogenetic comparative analyses usually rely on a single consensus phylogenetic tree in order to study evolutionary processes. However, most phylogenetic trees are incomplete with regard to species sampling, which may critically compromise analyses. Some approaches have been proposed to integrate non-molecular phylogenetic information into incomplete molecular phylogenies. An expanded tree approach consists of adding missing species to random locations within their clade. The information contained in the topology of the resulting expanded trees can be captured by the pairwise phylogenetic distance between species and stored in a matrix for further statistical analysis. Thus, the random expansion and processing of multiple phylogenetic trees can be used to estimate the phylogenetic uncertainty through a simulation procedure. Because of the computational burden required, unless this procedure is efficiently implemented, the analyses are of limited applicability. RESULTS: In this paper, we present efficient algorithms and implementations for randomly expanding and processing phylogenetic trees so that simulations involved in comparative phylogenetic analysis with uncertainty can be conducted in a reasonable time. We propose algorithms for both randomly expanding trees and calculating distance matrices. We made available the source code, which was written in the C++ language. The code may be used as a standalone program or as a shared object in the R system. The software can also be used as a web service through the link: http://purl.oclc.org/NET/sunplin/. CONCLUSION: We compare our implementations to similar solutions and show that significant performance gains can be obtained. Our results open up the possibility of accounting for phylogenetic uncertainty in evolutionary and ecological analyses of large datasets.

Phylogenetic fields of species: cross-species patterns of phylogenetic structure and geographical coexistence.

Differential coexistence among species underlies geographical patterns of biodiversity. Understanding such patterns has relied either on ecological or historical approaches applied separately. Recently, macroecology and community phylogenetics have tried to integrate both ecological and historical approaches. However, macroecology is mostly non-phylogenetic, whereas community phylogenetics is largely focused on local scales. Here, we propose a conceptual framework to link macroecology and community phylogenetics by exploring the evolutionary context of large-scale species coexistence, introducing the phylogenetic field concept. This is defined as the phylogenetic structure of species co-occurrence within a focal species' geographical range. We developed concepts and methods for analysing phylogenetic fields and applied them to study coexistence patterns of the bat family Phyllostomidae. Our analyses showed that phyllostomid bats coexist mostly with closely related species, revealing a north-south gradient from overdispersed to clustered phylogenetic fields. Patterns at different phylogenetic levels (i.e. all species versus close relatives only) presented the same gradient. Results support the tropical niche conservatism hypothesis, potentially mediated by higher speciation rates in the region of origin coupled with shared environmental preferences among species. The phylogenetic field approach enables species-based community phylogenetics, instead of those that are site-based, allowing the description of historical processes at more appropriate macroecological and biogeographic scales.

A new eigenfunction spatial analysis describing population genetic structure.

Several methods of spatial analyses have been proposed to infer the relative importance of evolutionary processes on genetic population structure. Here we show how a new eigenfunction spatial analysis can be used to model spatial patterns in genetic data. Considering a sample of n local populations, the method starts by modeling the response variable (allele frequencies or phenotypic variation) against the eigenvectors sequentially extracted from a geographic distance matrix (n × n). The relationship between the coefficient of determination (R(2)) of the models and the cumulative eigenvalues, which we named the spatial signal-representation (SSR) curve, can be more efficient than Moran's I correlograms in describing different patterns. The SSR curve was also applied to simulated data (under distinct scenarios of population differentiation) and to analyze spatial patterns in alleles from microsatellite data for 25 local populations of Dipteryx alata, a tree species endemic to the Brazilian Cerrado. The SSR curves are consistent with previous phylogeographical patterns of the species, revealing combined effects of isolation-by-distance and range expansion. Our analyses demonstrate that the SSR curve is a useful exploratory tool for describing spatial patterns of genetic variability and for selecting spatial eigenvectors for models aiming to explain spatial responses to environmental variables and landscape features.

Mantel test in population genetics.

The comparison of genetic divergence or genetic distances, estimated by pairwise FST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata ("Baru"), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.

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.