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

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

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

Unraveling the genetic diversity of Ceiba pubiflora (Malvaceae) in isolated limestone outcrops: Conservation strategies

Seasonally Dry Tropical Forests (SDTFs) located on limestone outcrops are vulnerable to degradation caused by timber logging and limestone extraction for cement production. Some of these forests represent the last remnants of native vegetation cover, functioning as isolated islands. Ceiba pubiflora (Malvaceae) is a tree frequently found on limestone outcrops in the central region of Brazil. This study aimed to evaluate the genetic diversity and identify suitable populations for the establishment of Management Units (MUs) for conservation. Inter-simple sequence repeat markers were employed to assess the genetic diversity in ten populations sampled from the Caatinga, Cerrado, and Atlantic Forest biomes. The species exhibited substantial genetic diversity (HT = 0.345; PLP = 97.89%). Populations SAH, JAN, and MON demonstrated elevated rates of polymorphic loci (> 84.2%) along with notable genetic diversity (He > 0.325). Additionally, these populations were the primary contributors to gene flow. The analysis of molecular variance (AMOVA) indicated that most genetic variation occurs within populations (91.5%) than between them. In the Bayesian analysis, the ten populations were clustered into five groups, revealing the presence of at least three barriers to gene flow in the landscape: 1) the Central Plateau or Paranã River valley; 2) near the Espinhaço mountains or the São Francisco River valley; and 3) around the Mantiqueira mountain range, Chapada dos Veadeiros plateau, and disturbed areas. A positive and statistically significant correlation was observed between genetic (θB) and geographic distances (r = 0.425, p = 0.008). Based on these findings, we propose the establishment of Management Units in Minas Gerais state, encompassing the (1) southern region (MIN population), (2) central region (SAH population), and (3) north region (MON population), as well as in Goiás state, covering the (4) Central Plateau region. These units can significantly contribute to preserving the genetic diversity of these trees and protecting their habitat against ongoing threats.

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae

Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Evidence for the Pleistocene Arc Hypothesis from genome-wide SNPs in a Neotropical dry forest specialist, the Rufous-fronted Thornbird (Furnariidae: Phacellodomus rufifrons)

South American dry forests have a complex and poorly understood biogeographic history. Based on the fragmented distribution of many Neotropical dry forest species, it has been suggested that this biome was more widely distributed and contiguous under drier climate conditions in the Pleistocene. To test this scenario, known as the Pleistocene Arc Hypothesis, we studied the phylogeography of the Rufous-fronted Thornbird (Phacellodomus rufifrons), a widespread dry forest bird with a disjunct distribution closely matching that of the biome itself. We sequenced mtDNA and used ddRADseq to sample 7,167 genome-wide single-nucleotide polymorphisms from 74 P. rufifrons individuals across its range. We found low genetic differentiation over two prominent geographic breaks - particularly across a 1,000 km gap between populations in Bolivia and Northern Peru. Using demographic analyses of the joint site frequency spectrum, we found evidence of recent divergence without subsequent gene flow across those breaks. By contrast, parapatric morphologically distinct populations in northeastern Brazil show high genetic divergence with evidence of recent gene flow. These results, in combination with our paleoclimate species distribution modelling, support the idea that currently disjunct patches of dry forest were more connected in the recent past, probably during the Middle and Late Pleistocene. This notion fits the major predictions of the Pleistocene Arc Hypothesis and illustrates the importance of comprehensive genomic and geographic sampling for examining biogeographic and evolutionary questions in complex ecosystems like Neotropical dry forests.

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.

Distributional patterns of Neotropical seasonally dry forest birds: a biogeographical regionalization

Neotropical seasonally dry forests (NSDFs) are widely distributed and possess high levels of species richness and endemism; however, their biogeography remains only partially understood. Using species distribution modelling and parsimony analysis of endemicity, we analysed the distributional patterns of the NSDF avifauna in order to identify their areas of endemism and provide a better understanding of the historical relationships among those areas. The strict consensus trees revealed 17 areas of endemism for NSDFs, which involve four large regions: Baja California, Caribbean-Antilles islands, Mesoamerica and South America. These well-resolved clades are circumscribed by geographical and ecological barriers associated with the Gulf of California, the leading edge of the Caribbean plate, the Tehuantepec Isthmus, the Polochic-Motagua fault, the Nicaragua Depression, the Chocó forest, the Amazon basin and the Andean Cordillera. Relationships among groups of NSDFs found here suggest that evolution of their avifauna involved a mixture of vicariance and dispersal events. Our results support the idea of independent diversification patterns and biogeographical processes in each region, including those previously associated with the Pleistocene Arc Hypothesis for NSDFs of south-eastern South America. This study provides a biogeographical framework to open new lines of research related to the biotic diversification of NSDFs.

What has become of the refugia hypothesis to explain biological diversity in Amazonia?

The spatial distribution of biodiversity and related processes is the core of Biogeography. Amazonia is the world's most diverse rainforest and the primary source of diversity to several Neotropical regions. The origins of such diversity continue to be an unresolved question in evolutionary biology. Among many competing hypotheses to explain the evolution of the Amazonian biodiversity, one stands out as the most influential: the refugia hypothesis by Jürgen Haffer. Here, we provide a chronological overview on how the refugia hypothesis evolved over the decades and how the criticism from different fields affected its acceptance. We conclude that the refugia hypothesis alone cannot explain the diversification of the complex Amazonian diversity, and perhaps it was not the most important diversification mechanism. However, the debate provoked by refugia has produced a great amount of knowledge on Amazonian climatic, geological, and evolutionary processes, as well as on species distributions, movements, and history.

Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

Aridity drives plant biogeographical sub regions in the Caatinga, the largest tropical dry forest and woodland block in South America

Our aims were to quantify and map the plant sub regions of the the Caatinga, that covers 844,453 km² and is the largest block of seasonally dry forest in South America. We performed spatial analyses of the largest dataset of woody plant distributions in this region assembled to date (of 2,666 shrub and tree species; 260 localities), compared these distributions with the current phytogeographic regionalizations, and investigated the potential environmental drivers of the floristic patterns in these sub regions. Phytogeographical regions were identified using quantitative analyses of species turnover calculated as Simpson dissimilarity index. We applied an interpolation method to map NMDS axes of compositional variation over the entire extent of the Caatinga, and then classified the compositional dissimilarity according to the number of biogeographical sub regions identified a priori using k-means analysis. We used multinomial logistic regression models to investigate the influence of contemporary climatic productivity, topographic complexity, soil characteristics, climate stability since the last glacial maximum, and the human footprint in explaining the identified sub regions. We identified nine spatially cohesive biogeographical sub regions. Current productivity, as indicated by an aridity index, was the only explanatory variable retained in the best model, explaining nearly half of the floristic variability between sub regions. The highest rates of endemism within the Caatinga were in the Core and Periphery Chapada Diamantina sub regions. Our findings suggest that the topographic complexity, soil variation, and human footprint in the Caatinga act on woody plant distributions at local scales and not as determinants of broad floristic patterns. The lack of effect of climatic stability since the last glacial maximum probably results from the fact that a single measure of climatic stability does not adequately capture the highly dynamic climatic shifts the region suffered during the Pleistocene. There was limited overlap between our results and previous Caatinga classifications.

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.

Do plant populations on distinct inselbergs talk to each other? A case study of genetic connectivity of a bromeliad species in an Ocbil landscape

Here, we explore the historical and contemporaneous patterns of connectivity among Encholirium horridum populations located on granitic inselbergs in an Ocbil landscape within the Brazilian Atlantic Forest, using both nuclear and chloroplast microsatellite markers. Beyond to assess the E. horridum population genetic structure, we built species distribution models across four periods (current conditions, mid-Holocene, Last Glacial Maximum [LGM], and Last Interglacial) and inferred putative dispersal corridors using a least-cost path analysis to elucidate biogeographic patterns. Overall, high and significant genetic divergence was estimated among populations for both nuclear and plastid DNA (ΦST(n) = 0.463 and ΦST(plastid) = 0.961, respectively, p < .001). For nuclear genome, almost total absence of genetic admixture among populations and very low migration rates were evident, corroborating with the very low estimates of immigration and emigration rates observed among E. horridum populations. Based on the cpDNA results, putative dispersal routes in Sugar Loaf Land across cycles of climatic fluctuations in the Quaternary period revealed that the populations' connectivity changed little during those events. Genetic analyses highlighted the low genetic connectivity and long-term persistence of populations, and the founder effect and genetic drift seemed to have been very important processes that shaped the current diversity and genetic structure observed in both genomes. The genetic singularity of each population clearly shows the need for in situ conservation of all of them.