Concordance between phylogeographical and biogeographical patterns in the Brazilian Cerrado: diversification of the endemic tree Dalbergia miscolobium (Fabaceae)

Phylogeography of Eriotheca species complex: insights into the origin and range expansion of apomictic and polyploid trees in Neotropical Savannas

Polyploidy and whole genome duplication are major evolutionary drivers in plants. Climate variations during the Pleistocene have influenced distribution and range expansion worldwide. Similar trends have been reported for Cerrado plants, but no attempt has been made to link phylogeography with ploidy and breeding changes. Thus, we aimed to (i) assess ploidy and genome size of Eriotheca estevesiae Carv.-Sobr., and compare it with E. pubescens (Mart.) Schott & Endl. (Both included into the Eriotheca Stellate Trichome Species Complex - ESTSC). (ii) Subsequently, we investigated their phylogeography to see whether genetic structure and range expansion trends were similar to those previously described for the Cerrado biome. Finally (iii), we discuss whether ESTSC phylogeographic patterns could be associated with geographic parthenogenesis processes. Common cytogenetic techniques and flow cytometry were used to confirm chromosome number and genome size of E. estevesiae. We used three cpDNA regions to analyse 14 ESTSC Cerrado populations, for which we also obtained ploidy level and breeding information. We investigated haplotype diversity, population structure and tested neutrality, aiming to reconstruct phylogeographic scenarios. We found three ploidy levels and eight cpDNA haplotypes in ESTSC, one shared by most populations. Haplotype and ploidy distribution corroborated that E. pubescens, the widely distributed polyploid and apomictic species, may have originated from northern diploid and probably sexual E. estevesiae. Matrilinear cpDNA links support the idea that apomixis and polyploidy in ESTSC may have allowed range expansion during the Pleistocene, in a process analogous to the geographic parthenogenesis described elsewhere.

Evolutionary history of Nasutitermes kemneri (Termitidae, Nasutitermitinae), a termite from the South American diagonal of open formations

Little is known about the phylogeography of termites in the Neotropical region. Here, we explored the genetic patterns and phylogeographical processes in the evolutionary history of Nasutitermes kemneri, an endemic termite of the South American diagonal of open formations (DOF) formed by the Chaco, Cerrado, and Caatinga phytogeographic domains. We sampled 60 individuals across the three domains of the DOF, and using the mitochondrial genes 16S, COI, and COII, as well as the nuclear gene ITS, evaluated the genetic diversity and divergence time of the populations, along with their genetic structure. The results show a strong genetic and spatial structure within the samples, evidencing the existence of two well-differentiated genetic groups: the Northeastern and the Southwestern populations, which diverged about 2.5 Mya, during the Pliocene-Pleistocene boundary. The Northeastern population, which encompasses Caatinga and northern portions of Cerrado, has an intricate structure and seems to have suffered repetitive retraction-expansion events due to climactic fluctuations during the Quaternary. The Southwestern population, which ranges from central-south Cerrado to the northeast peripherical portions of the Chaco, displays a star-shaped haplotype structure, indicating that this region may have acted as a refugia during interglacial periods.

Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree Plathymenia reticulata (n = 233 individuals). Genetic structure analysis of P. reticulata was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F CT = 0.222 and F ST = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F2) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in P. reticulata. The low genetic diversity of the P. reticulata forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

Stomatal size, ploidy and polyembryony in Eriotheca Stellate Trichome Species Complex (Bombacoideae - Malvaceae) in the Cerrados of Brazil

Geographic parthenogenesis, range expansion of apomictic plants after climate changes, has been described for Northern Hemisphere gametophytic apomicts. But similar trends have been observed for sporophytic apomicts of Cerrado, the savannas in Brazil. Eriotheca pubescens is a common Cerrado tree, an agamic complex of either hexaploid/polyembryonic apomicts or tetraploid/monoembryonic sexual individuals. Some populations have been described as a new species, Eriotheca estevesiae, all included in the Eriotheca Stellate Trichome Species Complex (ESTSC). Since breeding systems and ploidy are clearly associated with polyembryony and stomatal size, we used these ancillary features to map the reproductive and ploidy level traits of E. pubescens and E. estevesiae. Leaves and seeds were collected from individuals of 19 populations. Seeds were evaluated for the presence of polyembryony and leaves for stomatal measurements. Eight populations were monoembryonic while another eight were polyembryonic and for other three, the embryonic pattern was not readily verified. E. pubescens polyembryonic and hexaploid populations formed a homogeneous group, but monoembryonic plants were more variable. E. estevesiae populations were monoembryonic with smaller stomata. In contrast, some E. pubescens monoembryonic populations further south presented larger stomata. Despite these outliers, possibly mixed populations, stomatal size and embryonic pattern differed from northern to southern populations. Embryonic pattern and stomatal size indicated that northernmost populations of Eriotheca STSC (E. estevesiae) are diploid and sexual. Southernmost populations, mostly polyembryonic and with large stomata, are hexaploid and apomictic. This is in agreement with geographic parthenogenesis and range expansion of apomictic lineages to southern habitats available after the last glacial maximum.

Disentangling the Environmental Factors That Shape Genetic and Phenotypic Leaf Trait Variation in the Tree Qualea grandiflora Across the Brazilian Savanna

Identifying the environmental factors that shape intraspecific genetic and phenotypic diversity of species can provide insights into the processes that generate and maintain divergence in highly diverse biomes such as the savannas of the Neotropics. Here, we sampled Qualea grandiflora, the most widely distributed tree species in the Cerrado, a large Neotropical savanna. We analyzed genetic variation with microsatellite markers in 23 populations (418 individuals) and phenotypic variation of 10 metamer traits (internode, petiole and corresponding leaf lamina) in 36 populations (744 individuals). To evaluate the role of geography, soil, climate, and wind speed in shaping the divergence of genetic and phenotypic traits among populations, we used Generalized Dissimilarity Modelling. We also used multiple regressions to further investigate the contributions of those environmental factors on leaf trait diversity. We found high genetic diversity, which was geographically structured. Geographic distance was the main factor shaping genetic divergence in Qualea grandiflora, reflecting isolation by distance. Genetic structure was more related to past climatic changes than to the current climate. We also found high metamer trait variation, which seemed largely influenced by precipitation, soil bulk density and wind speed during the period of metamer development. The high degree of metamer trait variation seems to be due to both, phenotypic plasticity and local adaptation to different environmental conditions, and may explain the success of the species in occupying all the Cerrado biome.

Megafauna Seed Dispersal in the Neotropics: A Meta-Analysis Shows No Genetic Signal of Loss of Long-Distance Seed Dispersal

Restricted gene flow may lead to the loss of genetic diversity and higher genetic differentiation among populations, but the genetic consequences of megafauna extinction for plant populations still remain to be assessed. We performed a phylogenetic-independent meta-analysis across 102 Neotropical plants to test the hypothesis that plant species with megafaunal seed dispersal syndrome have a lower genetic diversity and a higher genetic differentiation than those without it. We classified as megafauna-dependent plant species those that potentially relied only on megafauna to seed dispersal, and as megafauna-independent those that relied on megafauna and other seed dispersers. Our data comprised 98 studies using microsatellite markers. We found no statistical difference in genetic diversity and differentiation between plants with megafauna and non-megafauna seed dispersal syndrome, although the statistical power to detect differences in genetic differentiation was low. Moreover, we found no statistical difference between megafauna-dependent and megafauna-independent plant species. We then used generalized linear mixed models and phylogenetic generalized least square models to investigate the effects of megafaunal seed dispersal syndromes and reproductive traits on variation in genetic diversity and genetic differentiation. We found no effect of megafaunal syndrome, rather, reproductive traits, such as pollination mode, mating, and breeding systems, showed significant effects. Our findings show that the genetic studies of Neotropical plants performed so far show no difference in genetic diversity and differentiation in plants with megafaunal compared to those with non-megafaunal seed dispersal syndromes. Our results also provide evidence pointing out that plant species with megafaunal seed dispersal syndromes may have used different strategies to counterbalance the extinction of their mutualistic megafauna dispersers, such as the dispersal by extant mammals that may promote long-distance seed dispersal. Our results also reinforce the importance of pollination to long-distance gene flow in Neotropical plants.

Habitat preference differentiates the Holocene range dynamics but not barrier effects on two sympatric, congeneric trees (Tristaniopsis, Myrtaceae)

Niche partitioning can lead to differences in the range dynamics of plant species through its impacts on habitat availability, dispersal, or selection for traits that affect colonization and persistence. We investigated whether niche partitioning into upland and riparian habitats differentiates the range dynamics of two closely related and sympatric eastern Australian trees: the mountain water gum (Tristaniopsis collina) and the water gum (T. laurina). Using genomic data from SNP genotyping of 480 samples, we assessed the impact of biogeographic barriers and tested for signals of range expansion. Circuit theory was used to model isolation-by-resistance across three palaeo-environment scenarios: the Last Glacial Maximum, the Holocene Climate Optimum and present-day (1950-2014). Both trees showed similar genetic structure across historically dry barriers, despite evidence of significant environmental niche differentiation and different post-glacial habitat shifts. Tristaniopsis collina exhibits the signature of serial founder effects consistent with recent or rapid range expansion, whilst T. laurina has genetic patterns consistent with long-term persistence in geographically isolated populations despite occupying a broader bioclimatic niche. We found the minor influence of isolation-by-resistance on both species, though other unknown factors appear to shape genetic variation. We postulate that specialized recruitment traits (adapted to flood-disturbance regimes) rather than habitat availability limited post-glacial range expansion in T. laurina. Our findings indicate that niche breadth does not always facilitate range expansion through colonization and migration across barriers, though it can promote long-term persistence in situ.

Genetic data improve the assessment of the conservation status based only on herbarium records of a Neotropical tree

Although there is a consensus among conservation biologists about the importance of genetic information, the assessment of extinction risk and conservation decision-making generally do not explicitly consider this type of data. Genetic data can be even more important in species where little other information is available. In this study, we investigated a poorly known legume tree, Dimorphandra exaltata, from the Brazilian Atlantic Forest, a hotspot for conservation. We coupled species distribution models and geospatial assessment based on herbarium records with population genetic analyses to evaluate its genetic status and extinction risk, and to suggest conservation measures. Dimorphandra exaltata shows low genetic diversity, inbreeding, and genetic evidence of decrease in population size, indicating that the species is genetically depleted. Geospatial assessment classified the species as Endangered. Species distribution models projected a decrease in range size in the near future (2050). The genetic status of the species suggests low adaptive potential, which compromises its chances of survival in the face of ongoing climatic change. Altogether, our coupled analyses show that the species is even more threatened than indicated by geospatial analyses alone. Thus, conservation measures that take into account genetic data and the impacts of climate change in the species should be implemented.

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.

Genetic and Historical Colonization Analyses of an Endemic Savanna Tree, Qualea grandiflora, Reveal Ancient Connections Between Amazonian Savannas and Cerrado Core

The evolutionary processes underlying the high diversity and endemism in the Cerrado, the most extensive Neotropical savanna, remain unclear, including the factors promoting the presence and evolution of savanna enclaves in the Amazon forest. In this study, we investigated the effects of past climate changes on genetic diversity, dynamics of species range and the historical connections between the savanna enclaves and Cerrado core for Qualea grandiflora, a tree species widely distributed in the biome. Totally, 40 populations distributed in the Cerrado core and Amazon savannas were analyzed using chloroplast and nuclear DNA sequences. We used phylogeographic, coalescent and ecological niche modeling approaches. Genetic data revealed a phylogeographic structure shaped by Pleistocene climatic oscillations. An eastern-western split in the Cerrado core was observed. The central portion of the Cerrado core harbored most of the sampled diversity for cpDNA. Ecological niche models predicted the presence of a large historical refuge in this region and multiple small refuges in peripheral areas. Relaxed Random Walk (RRW) models indicated the ancestral population in the north-western border of the central portion of the Cerrado core and cyclical dynamics of colonization related to Pleistocene climatic oscillations. Central and western ancient connections between Cerrado core and Amazonian savannas were observed. No evidence of connections among the Amazonian savannas was detected. Our study highlights the importance of Pleistocene climatic oscillations for structuring the genetic diversity of Q. grandiflora and complex evolutionary history of ecotonal areas in the Cerrado. Our results do not support the recent replacement of a large area in the Amazon forest by savanna vegetation. The Amazonian savannas appear to be fragmented and isolated from each other, evolving independently a long ago.

The Calibrated Phylogeny of the Drosophila fasciola Subgroup (D. repleta Group Wasserman) Indicates Neogene Diversification of Its Internal Branches

The species of the Drosophila fasciola subgroup Wasserman represent the dominant section of the Drosophila repleta group Wasserman in the American rainforests and have a broad geographical distribution in the New World. However, despite of its wide range, the D. fasciola subgroup is one of the most overlooked D. repleta subgroups. Here, we report a molecular phylogenetic analysis focused on the D. fasciola subgroup using two mitochondrial [cytochrome oxidase subunit I (COI), cytochrome oxidase subunit II (COII)] and two nuclear [elongation factor-1alpha F1 (EF-alphaF1) and transformer (tra)] genes. Overall, we found that this subgroup is a monophyletic taxon, subdivided into two main internal branches: named Fas1 and Fas2 clades. The diversification of these clades is estimated to have begun in the middle Miocene, around 12 Ma [95% high posterior density (HPD) 9.0-15 Ma], and might be associated with the colonization of South America by Central America populations after the closure of Isthmus of Panama due to the temporal congruence between these events. The terminal branches had their origins estimated to be in the Pliocene or the Plio-Pleistocene transition. For the later estimates, both the geomorphological influences and the climatic oscillations of the Pleistocene may have played a role in shaping the diversification of the D. fasciola group.

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.

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east-west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.

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.