A mid‐ P leistocene rainforest corridor enabled synchronous invasions of the A tlantic F orest by A mazonian anole lizards

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.

Complete phylogeny of Micrathena spiders suggests multiple dispersal events among Neotropical rainforests, islands and landmasses, and indicates that Andean orogeny promotes speciation

The Neotropical region is the most diverse on the planet, largely owing to its mosaic of tropical rainforests. Multiple tectonic and climatic processes have been hypothesized to contribute to generating this diversity, including Andean orogeny, the closure of the Isthmus of Panama, the GAARlandia land bridge and historical connections among currently isolated forests. Micrathena spiders are diverse and widespread in the region, and thus a complete phylogeny of this genus allows the testing of hypotheses at multiple scales. We estimated a complete, dated phylogeny using morphological data for 117 Micrathena species and molecular data of up to five genes for a subset of 79 species. Employing eventc-based approaches and biogeographic stochastic mapping while considering phylogenetic uncertainty, we estimated ancestral distributions, the timing and direction of dispersal events and diversification rates among areas. The phylogeny is generally robust, with uncertainty in the position of some of the species lacking sequences. Micrathena started diversifying around 25 Ma. Andean cloud forests show the highest in-situ speciation, while the Amazon is the major dispersal source for adjacent areas. The Dry Diagonal generated few species and is a sink of diversity. Species exchange between Central and South America involved approximately 23 dispersal events and started ~20 Ma, which is consistent with a Miocene age for the Isthmus of Panama closure. We inferred four dispersal events from Central America to the Antilles in the last 20 Myr, indicating the spiders did not reach the islands through the GAARlandia land bridge. We identified important species exchange routes among the Amazon, Andean cloud forests and Atlantic forests during the Plio-Pleistocene. Sampling all species of the genus was fundamental to the conclusions above, especially in identifying the Andean forests as the area that generated the majority of species. This highlights the importance of complete taxonomic sampling in biogeographic studies.

Biogeographic factors contributing to the diversification of Euphoniinae (Aves, Passeriformes, Fringillidae): a phylogenetic and ancestral areas analysis

Factors such as the Andean uplift, Isthmus of Panama, and climate changes have influenced bird diversity in the Neotropical region. Studying bird species that are widespread in Neotropical highlands and lowlands can help us understand the impact of these factors on taxa diversification. Our main objectives were to determine the biogeographic factors that contributed to the diversification of Euphoniinae and re-evaluate their phylogenetic relationships. The nextRAD and mitochondrial data were utilized to construct phylogenies. The ancestral distribution range was then estimated using a time-calibrated phylogeny, current species ranges, and neotropical regionalization. The phylogenies revealed two main Euphoniinae clades, Chlorophonia and Euphonia, similar to previous findings. Furthermore, each genus has distinctive subclades corresponding to morphology and geography. The biogeographic results suggest that the Andean uplift and the establishment of the western Amazon drove the vicariance of Chlorophonia and Euphonia during the Miocene. The Chlorophonia lineage originated in the Andes mountains and spread to Central America and the Mesoamerican highlands after the formation of the Isthmus of Panama. Meanwhile, the ancestral area of Euphonia was the Amazonas, from which it spread to trans-Andean areas during the Pliocene and Pleistocene due to the separation of the west lowlands from Amazonas due to the Northern Andean uplift. Chlorophonia and Euphonia species migrated to the Atlantic Forest during the Pleistocene through corridors from the East Andean Humid Forest and Amazonas. These two genera had Caribbean invasions with distinct geographic origins and ages. Finally, we suggested taxonomic changes in the genus Euphonia based on the study's phylogenetic, morphological, and biogeographic findings.

Amazonian colonization from the Atlantic Forest: New perspectives on the connections of South American tropical forests

An open and dry vegetation belt separates Amazonia (AM) and the Atlantic Forest (AF). Evidence from palaeoclimatic and phylogenetic studies suggests past connections between these forests during cycles of increased humidity through the formation of forest corridors. The distinctive northern AF avifauna is known to have affinities both with AM and the southern AF. Still, the extent of how these two regions contributed to the assemblage of this avifauna remains poorly understood. Using historical demographic analyses and comparative phylogeography based on sub-genomic genetic sampling, we assessed how past connections between AM and AF led to shared vicariance and colonization events in four avian AF endemic taxa. Our results supported the occurrence of humid forest corridors promoting the contact between AF and AM populations and suggested two vicariant events and two colonization events from AF to AM. Population divergences were mostly non-synchronous and occurred multiple times during the Pleistocene. Historical gene flow was prevalent across study groups, supporting migration flows after the initial separation between AM and AF - a pattern previously unknown in birds between these regions. Idiosyncratic histories and divergent demographic syndromes suggest that organisms' responses to climate-driven habitat shifts broadly depend on their ecological attributes. This study strengthened our knowledge of past connections between AM and AF and provided demographic scenarios amenable for testing in other groups of co-distributed organisms.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

Phylogeography of a Typical Forest Heliothermic Lizard Reveals the Combined Influence of Rivers and Climate Dynamics on Diversification in Eastern Amazonia

The formation of the Amazon drainage basin has been considered an important driver of speciation of several taxa, promoting vicariant events or reinforcement of barriers that restrict gene flow between opposite river margins. Several recent studies reported a set of miscellaneous events involving climatic fluctuations, geomorphological changes, and dispersal mechanisms as propellers of diversification of Amazonian rainforest taxa. Here, we show the results of dated phylogenetic, biogeographic, and populational analyses to investigate which events could better explain the current distribution of a heliothermic, active foraging lizard in the central and eastern portions of the Amazonian rainforest (besides a disjunct distribution in part of the Atlantic Forest). We sampled Kentropyx calcarata from most of its area of occurrence in Amazonia and used mitochondrial and nuclear markers to evaluate if the genetic structure agrees with evolutionary scenarios previously proposed for Amazonia. We performed phylogenetic and populational analyses to better understand the dynamics of this species in the Amazonia rainforest over time. Phylogenetic inference recovered ten K. calcarata structured lineages in eastern Amazonia, some of them limited by the Amazon River and its southern tributaries (Tapajós, Xingu, and Tocantins), although we detected occasional haplotype sharing across some of the river banks. According to molecular dating, K. calcarata diversified since Miocene–Pliocene, and some of the lineages presented signs of demographic expansion during the Pleistocene, supposedly triggered by climatic dynamics. The putative ancestral lineage of K. calcarata was distributed on the Guiana Shield, later spreading south and southeastward by dispersion. Our results indicate that Amazonian rivers acted as barriers to the dispersal of Kentropyx calcarata, but they were not the sole drivers of diversification.

A multidisciplinary framework for biodiversity prediction in the Brazilian Atlantic Forest hotspot

Abstract: We briefly describe selected results from our thematic project focused on the biodiversity of the Atlantic Forest (“AF BIOTA”), which was jointly funded by FAPESP’s BIOTA Program, the U.S. National Science Foundation Dimensions of Biodiversity Program, and the National Aeronautics and Space Administration (NASA). As one of the five most important hotspots of biodiversity in the world, the Atlantic Forest (AF) holds less than 16% of its vegetation cover, yet, amongst the hotspots, it still harbors one of the highest numbers of species, including endemics. By gathering specialists across multiple disciplines (biology, geology, engineering), we aimed to understand how this megabiodiversity was built through time, informing biodiversity science and conservation. Among the results, we trained 18 Master’s and 26 Ph.D. students, published more than 400 peer-reviewed papers that improved our knowledge about the forest’s biologic and climatic diversity and dynamics through time, developed new analytical methods, produced outreach videos and articles, and provided data to help define biodiversity conservation policies.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Interspecific Gene Flow and Mitochondrial Genome Capture During the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae)

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, A. opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting. Here we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene-tree, species-tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene-tree recovers the same relationships previously inferred for this group, which is strikingly different from the species-tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model, so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group, however coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent adaptive radiations.

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon-specific attributes.

The role of vicariance and dispersal on the temporal range dynamics of forest vipers in the Neotropical region

The emergence of the diagonal of open/dry vegetations, including Chaco, Cerrado and Caatinga, is suggested to have acted as a dispersal barrier for terrestrial organisms by fragmenting a single large forest that existed in South America into the present Atlantic and Amazon forests. Here we tested the hypothesis that the expansion of the South American diagonal of open/dry landscapes acted as a vicariant process for forest lanceheads of the genus Bothrops, by analyzing the temporal range dynamics of those snakes. We estimated ancestral geographic ranges of the focal lancehead clade and its sister clade using a Bayesian dated phylogeny and the BioGeoBEARS package. We compared nine Maximum Likelihood models to infer ancestral range probabilities and their related biogeographic processes. The best fitting models (DECTS and DIVALIKETS) recovered the ancestor of our focal clade in the Amazon biogeographic region of northwestern South America. Vicariant processes in two different subclades resulted in disjunct geographic distributions in the Amazon and the Atlantic Forest. Dispersal processes must have occurred mostly within the Amazon and the Atlantic Forest and not between them. Our results suggest the fragmentation of a single ancient large forest into the Atlantic and Amazon forests acting as a driver of vicariant processes for the snake lineage studied, highlighting the importance of the diagonal of open/dry landscapes in shaping distribution patterns of terrestrial biota in South America.

Genetic variation in neotropical butterflies is associated with sampling scale, species distributions, and historical forest dynamics

Previous studies of butterfly diversification in the Neotropics have focused on Amazonia and the tropical Andes, while southern regions of the continent have received little attention. To address the gap in knowledge about the Lepidoptera of temperate South America, we analysed over 3000 specimens representing nearly 500 species from Argentina for a segment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Representing 42% of the country's butterfly fauna, collections targeted species from the Atlantic and Andean forests, and biodiversity hotspots that were previously connected but are now isolated. We assessed COI effectiveness for species discrimination and identification and how its performance was affected by geographic distances and taxon coverage. COI data also allowed to study patterns of genetic variation across Argentina, particularly between populations in the Atlantic and Andean forests. Our results show that COI discriminates species well, but that identification success is reduced on average by ~20% as spatial and taxonomic coverage rises. We also found that levels of genetic variation are associated with species' spatial distribution type, a pattern which might reflect differences in their dispersal and colonization abilities. In particular, intraspecific distance between populations in the Atlantic and Andean forests was significantly higher in species with disjunct distributions than in those with a continuous range. All splits between lineages in these forests dated to the Pleistocene, but divergence dates varied considerably, suggesting that historical connections between the Atlantic and Andean forests have differentially affected their shared butterfly fauna. Our study supports the fact that large-scale assessments of mitochondrial DNA variation are a powerful tool for evolutionary studies.

Mesoamerica is a cradle and the Atlantic Forest is a museum of Neotropical butterfly diversity: insights from the evolution and biogeography of Brassolini (Lepidoptera: Nymphalidae)

Regional species diversity is explained ultimately by speciation, extinction and dispersal. Here, we estimate dispersal and speciation rates of Neotropical butterflies to propose an explanation for the distribution and diversity of extant species. We focused on the tribe Brassolini (owl butterflies and allies), a Neotropical group that comprises 17 genera and 108 species, most of them endemic to rainforest biomes. We inferred a robust species tree using the multispecies coalescent framework and a dataset including molecular and morphological characters. This formed the basis for three changes in Brassolini classification: (1) Naropina syn. nov. is subsumed within Brassolina; (2) Aponarope syn. nov. is subsumed within Narope; and (3) Selenophanes orgetorix comb. nov. is reassigned from Catoblepia to Selenophanes. By applying biogeographical stochastic mapping, we found contrasting species diversification and dispersal dynamics across rainforest biomes, which might be explained, in part, by the geological and environmental history of each bioregion. Our results revealed a mosaic of biome-specific evolutionary histories within the Neotropics, where butterfly species have diversified rapidly (cradles: Mesoamerica), have accumulated gradually (museums: Atlantic Forest) or have diversified and accumulated alternately (Amazonia). Our study contributes evidence from a major butterfly lineage that the Neotropics are a museum and a cradle of species diversity.

Potential mammalian species for investigating the past connections between Amazonia and the Atlantic Forest

Much evidence suggests that Amazonia and the Atlantic Forest were connected through at least three dispersion routes in the past: the Eastern route, the central route, and the Western route. However, few studies have assessed the use of these routes based on multiple species. Here we present a compilation of mammal species that potentially have dispersed between the two forest regions and which may serve to investigate these connections. We evaluate the present-day geographic distributions of mammals occurring in both Amazonia and the Atlantic Forest and the likely connective routes between these forests. We classified the species per habitat occupancy (strict forest specialists, species that prefer forest habitat, or generalists) and compiled the genetic data available for each species. We found 127 mammalian species presently occurring in both Amazonia and the Atlantic Forest for which, substantial genetic data was available. Hence, highlighting their potential for phylogeographic studies investigating the past connections between the two forests. Differently from what was previously proposed, the present-day geographic distribution of mammal species found in both Amazonia and the Atlantic Forest points to more species in the eastern portion of the dry diagonal (and adjoining forested habitats). The Central route was associated with the second most species. Although it remains to be seen how this present-day geography reflects the paleo dispersal routes, our results show the potential of using mammal species to investigate and bring new insights about the past connections between Amazonia and the Atlantic Forest.

Comparative phylogeographic and demographic analyses reveal a congruent pattern of sister relationships between bird populations of the northern and south-central Atlantic Forest

The Pernambuco Center of Endemism (PCE) is the northernmost strip of the Atlantic Forest (AF). Biogeographic affinities among avifaunas in the PCE, the southern-central Atlantic Forest (SCAF), and Amazonia (AM) have not been studied comprehensively, and current patterns of genetic diversity in the PCE remain unclear. The interplay between species' ecological attributes and historical processes, such as Pleistocene climate fluctuations or the appearance of rivers, may have affected population genetic structures in the PCE. Moreover, the role of past connections between the PCE and AM and the elevational distribution of species in assembling the PCE avifauna remain untested. Here, we investigated the biogeographic history of seven taxa endemic to the PCE within a comparative phylogeographic framework based on a mean of 3,618 independent single nucleotide polymorphisms (SNPs) extracted from flanking regions of ultraconserved elements (UCEs) and one mitochondrial gene. We found that PCE populations were more closely related to SCAF populations than they were to those in AM, regardless of their elevational range, with divergence times placed during the Mid-Pleistocene. These splits were consistent with a pattern of allopatric divergence with gene flow until the upper Pleistocene and no signal of rapid changes in population sizes. Our results support the existence of a Pleistocene refugium driving current genetic diversity in the PCE, thereby rejecting the role of the São Francisco River as a primary barrier for population divergence. Additionally, we found that connections with Amazonia also played a significant role in assembling the PCE avifauna through subsequent migration events.

Molecular phylogenetic inference of the howler monkey radiation (Primates: Alouatta)

Howler monkeys (Alouatta), comprising between nine and 14 species and ranging from southern Mexico to northern Argentina, are the most widely distributed platyrrhines. Previous phylogenetic studies of howlers have used chromosomal and morphological characters and a limited number of molecular markers; however, branching patterns conflict between studies or remain unresolved. We performed a new phylogenetic analysis of Alouatta using both concatenated and coalescent-based species tree approaches based on 14 unlinked non-coding intergenic nuclear regions. Our taxon sampling included five of the seven South American species (Alouatta caraya, Alouatta belzebul, Alouatta guariba, Alouatta seniculus, Alouatta sara) and the two recognized species from Mesoamerica (Alouatta pigra, Alouatta palliata). Similarly to previous studies, our phylogenies supported a Mesoamerican clade and a South American clade. For the South American howlers, both methods recovered the Atlantic Forest endemic A. guariba as sister to all remaining South American species, albeit with moderate support. Moreover, we found no support for the previously proposed sister relationship between A. guariba and A. belzebul. For the first time, a clade composed of A. sara and A. caraya was identified. The relationships among the other South American howlers, however, were not fully supported. Our estimates for divergence times within Alouatta are generally older compared to estimates in earlier studies. However, they conform to recent studies proposing a Miocene age for the Isthmus of Panama and for the uplift of the northern Andes. Our results also point to an early genetic isolation of A. guariba in the Atlantic Forest, in agreement with the hypothesis of biotic exchange across South American rain forests in the Miocene. Collectively, these findings contribute to a better understanding of the diversification processes among howler monkey species; however, they also suggest that further comprehension of the evolutionary history of the Alouatta radiation will rely on broadened taxonomic, geographic, and genomic sampling.

Diversification history of clown tree frogs in Neotropical rainforests (Anura, Hylidae, Dendropsophus leucophyllatus group)

General consensus emphasizes that no single biological process can explain the patterns of species' distributions and diversification in the Neotropics. Instead, the interplay of several processes across space and time must be taken into account. Here we investigated the phylogenetic relationships and biogeographic history of tree frogs in the Dendropsophus leucophyllatus species group (Amphibia: Hylidae), which is distributed across Amazonia and the Atlantic rainforests. Using Next Generation Sequencing (NGS) and double digest restriction-site associated DNA (ddRADseq), we inferred phylogenetic relationships, species limits, and temporal and geographic patterns of diversification relative to the history of these biomes. Our results indicate that the D. leucophyllatus species group includes at least 14 independent lineages, which are currently arranged into ten described species. Therefore, a significant portion of species in the group are still unnamed. Different processes were associated to the group diversification history. For instance, the Andes uplift likely caused allopatric speciation for Cis-Andean species, whereas it may also be responsible for changes in the Amazonian landscape triggering parapatric speciation by local adaptation to ecological factors. Meanwhile, Atlantic Forest ancestors unable to cross the dry diagonal biomes after rainforest's retraction, evolved in isolation into different species. Diversification in the group began in the early Miocene, when connections between Atlantic Forest and the Andes (Pacific Dominion) by way of a south corridor were possible. The historical scenario in Amazonia, characterized by several speciation events and habitat heterogeneity, helped promoting diversification, resulting in the highest species diversity for the group. This marked species diversification did not happen in Atlantic Forest, where speciation is very recent (late Pliocene and Pleistocene), despite its remarkable climatic heterogeneity.

Why is Amazonia a 'source' of biodiversity? Climate-mediated dispersal and synchronous speciation across the Andes in an avian group (Tityrinae)

Amazonia is a 'source' of biodiversity for other Neotropical ecosystems, but which conditions trigger in situ speciation and emigration is contentious. Three hypotheses for how communities have assembled include (1) a stochastic model wherein chance dispersal events lead to gradual emigration and species accumulation, (2) diversity-dependence wherein successful dispersal events decline through time due to ecological limits, and (3) barrier displacement wherein environmental change facilitates dispersal to other biomes via transient habitat corridors. We sequenced thousands of molecular markers for the Neotropical Tityrinae (Aves) and applied a novel filtering protocol to identify loci with high utility for dated phylogenomics. We used these loci to estimate divergence times and model Tityrinae's evolutionary history. We detected a prominent role for speciation driven by barriers including synchronous speciation across the Andes and found that dispersal increased toward the present. Because diversification was continuous but dispersal was non-random over time, we show that barrier displacement better explains Tityrinae's history than stochasticity or diversity-dependence. We propose that Amazonia is a source of biodiversity because (1) it is a relic of a biome that was once more extensive, (2) environmentally mediated corridors facilitated emigration and (3) constant diversification is attributed to a spatially heterogeneous landscape that is perpetually dynamic through time.

The combined role of dispersal and niche evolution in the diversification of Neotropical lizards

Ecological requirements and environmental conditions can influence diversification across temporal and spatial scales. Understanding the role of ecological niche evolution under phylogenetic contexts provides insights on speciation mechanisms and possible responses to future climatic change. Large-scale phyloclimatic studies on the megadiverse Neotropics, where biomes with contrasting vegetation types occur in narrow contact, are rare. We integrate ecological and biogeographic data with phylogenetic comparative methods, to investigate the relative roles of biogeographic events and niche divergence and conservatism on the diversification of the lizard genus Kentropyx Spix, 1825 (Squamata: Teiidae), distributed in South American rainforests and savannas. Using five molecular markers, we estimated a dated species tree, which recovered three clades coincident with previously proposed species groups diverging during the mid-Miocene. Biogeography reconstruction indicates a role of successive dispersal events from an ancestral range in the Brazilian Shield and western Amazonia. Ancestral reconstruction of climatic tolerances and niche overlap metrics indicates a trend of conservatism during the diversification of groups from the Amazon Basin and Guiana Shield, and a strong signal of niche divergence in the Brazilian Shield savannas. Our results suggest that climatic-driven divergence at dynamic forest-savanna borders might have resulted in adaptation to new environmental niches, promoting habitat shifts and shaping speciation patterns of Neotropical lizards. Dispersal and ecological divergence could have a more important role in Neotropical diversification than previously thought.

Historical signatures in the alpha and beta diversity patterns of Atlantic Forest harvestman communities (Arachnida: Opiliones)

The integration of ecology and historical biogeography is fostering the investigation of diversity patterns. We studied alpha and beta diversity patterns of Brazilian Atlantic Forest harvestman (Arachnida: Opiliones) communities and related them to environmental and historical factors. Our data bank contains 508 species from 63 sites, encompassing almost the entire latitudinal range of Atlantic Forest. Alpha diversity was higher in coastal sites in the south and southeast regions and decreased in sites inland, as well as in sites in the coastal northeast region, especially in northern Bahia state. Alpha diversity was positively influenced by precipitation and altitudinal range, but the low number of species in northeastern coastal sites seems to be more related to the historical distribution of Neotropical harvestman lineages than to recent environmental factors. Geographic distance was the most influential factor for beta diversity. Compositional changes were also remarkably congruent with areas of endemism delimited for Atlantic Forest harvestmen. The percentage of protected areas for each area of endemism was very unbalanced, and Espírito Santo and Pernambuco states were the least protected areas. The turnover process observed in the compositional changes indicates that conservation strategies should include as many reserves as possible because every community presents a unique set of species.

Local adaptation in mainland anole lizards: Integrating population history and genome-environment associations

Environmental gradients constrain physiological performance and thus species' ranges, suggesting that species occurrence in diverse environments may be associated with local adaptation. Genome-environment association analyses (GEAA) have become central for studies of local adaptation, yet they are sensitive to the spatial orientation of historical range expansions relative to landscape gradients. To test whether potentially adaptive genotypes occur in varied climates in wide-ranged species, we implemented GEAA on the basis of genomewide data from the anole lizards Anolis ortonii and Anolis punctatus, which expanded from Amazonia, presently dominated by warm and wet settings, into the cooler and less rainy Atlantic Forest. To examine whether local adaptation has been constrained by population structure and history, we estimated effective population sizes, divergence times, and gene flow under a coalescent framework. In both species, divergence between Amazonian and Atlantic Forest populations dates back to the mid-Pleistocene, with subsequent gene flow. We recovered eleven candidate genes involved with metabolism, immunity, development, and cell signaling in A. punctatus and found no loci whose frequency is associated with environmental gradients in A. ortonii. Distinct signatures of adaptation between these species are not associated with historical constraints or distinct climatic space occupancies. Similar patterns of spatial structure between selected and neutral SNPs along the climatic gradient, as supported by patterns of genetic clustering in A. punctatus, may have led to conservative GEAA performance. This study illustrates how tests of local adaptation can benefit from knowledge about species histories to support hypothesis formulation, sampling design, and landscape gradient characterization.

Conceptual and empirical advances in Neotropical biodiversity research

The unparalleled biodiversity found in the American tropics (the Neotropics) has attracted the attention of naturalists for centuries. Despite major advances in recent years in our understanding of the origin and diversification of many Neotropical taxa and biotic regions, many questions remain to be answered. Additional biological and geological data are still needed, as well as methodological advances that are capable of bridging these research fields. In this review, aimed primarily at advanced students and early-career scientists, we introduce the concept of "trans-disciplinary biogeography," which refers to the integration of data from multiple areas of research in biology (e.g., community ecology, phylogeography, systematics, historical biogeography) and Earth and the physical sciences (e.g., geology, climatology, palaeontology), as a means to reconstruct the giant puzzle of Neotropical biodiversity and evolution in space and time. We caution against extrapolating results derived from the study of one or a few taxa to convey general scenarios of Neotropical evolution and landscape formation. We urge more coordination and integration of data and ideas among disciplines, transcending their traditional boundaries, as a basis for advancing tomorrow's ground-breaking research. Our review highlights the great opportunities for studying the Neotropical biota to understand the evolution of life.

Deep divergence of Red-crowned Ant Tanager (Habia rubica: Cardinalidae), a multilocus phylogenetic analysis with emphasis in Mesoamerica

Many neotropical species have a complex history of diversification as a result of the influence of geographical, ecological, climatic, and geological factors that determine the distribution of populations within a lineage. Phylogeography identifies such populations, determines their geographic distributions, and quantifies the degree of genetic divergence. In this work we explored the genetic structure of Habia rubica populations, a polytypic taxon with 17 subspecies described, in order to obtain hypotheses about their evolutionary history and processes of diversification. We undertook multilocus analyses using sequences of five molecular markers (ND2, ACOI-I9, MUSK, FGB-I5 and ODC), and sampling from across the species’ distribution range, an area encompassing from Central Mexico throughout much of South America. With these data, we obtained a robust phylogenetic hypothesis, a species delimitation analysis, and estimates of divergence times for these lineages. The phylogenetic hypothesis of concatenated molecular markers shows that H. rubica can be divided in three main clades: the first includes Mexican Pacific coast populations, the second is formed by population from east of Mexico to Panama and the third comprises the South American populations. Within these clades we recognize seven principal phylogroups whose limits have a clear correspondence with important geographical discontinuities including the Isthmus of Tehuantepec in southern Mexico, the Talamanca Cordillera, and the Isthmus of Panama in North America. In South America, we observed a marked separation of two phylogroups that include the populations that inhabit mesic forests in western and central South America (Amazon Forest) and those inhabiting the seasonal forest from the eastern and northern regions of the South America (Atlantic Forest). These areas are separated by an intervening dry vegetation “diagonal” (Chaco, Cerrado and Caatinga). The geographic and genetic structure of these phylogroups describes a history of diversification more active and complex in the northern distribution of this species, producing at least seven well-supported lineages that could be considered species.

Modeling the ecology and evolution of biodiversity: Biogeographical cradles, museums, and graves

Individual processes shaping geographical patterns of biodiversity are increasingly understood, but their complex interactions on broad spatial and temporal scales remain beyond the reach of analytical models and traditional experiments. To meet this challenge, we built a spatially explicit, mechanistic simulation model implementing adaptation, range shifts, fragmentation, speciation, dispersal, competition, and extinction, driven by modeled climates of the past 800,000 years in South America. Experimental topographic smoothing confirmed the impact of climate heterogeneity on diversification. The simulations identified regions and episodes of speciation (cradles), persistence (museums), and extinction (graves). Although the simulations had no target pattern and were not parameterized with empirical data, emerging richness maps closely resembled contemporary maps for major taxa, confirming powerful roles for evolution and diversification driven by topography and climate.

Origins and biogeography of the Anolis crassulus subgroup (Squamata: Dactyloidae) in the highlands of Nuclear Central America

BACKGROUND

Recent studies have begun to reveal the complex evolutionary and biogeographic histories of mainland anoles in Central America, but the origins and relationships of many taxa remain poorly understood. One such group is the Anolis (Norops) crassulus species subgroup, which contains ten morphologically similar highland taxa, the majority of which have restricted distributions. The nominal taxon A. crassulus has a disjunct distribution from Chiapas, Mexico, through Guatemala, in the highlands of El Salvador, and in the Chortís Highlands of Honduras. We test the relationships of these species using multiple mitochondrial and nuclear loci in concatenated and multispecies coalescent frameworks, in an effort to both resolve long-standing taxonomic confusion and present new insights into the evolution and biogeography of these taxa.

RESULTS

Sequences of multiple mitochondrial and nuclear loci were generated for eight of the ten species of the Anolis crassulus species subgroup. We analyzed phylogenetic relationships and estimated divergence times and ancestral ranges of the subgroup, recovering a monophyletic subgroup within Anolis. Within the nominal taxon Anolis crassulus, we recovered multiple genetically distinct lineages corresponding to allopatric populations, and show that the Chortís Highland lineage split from the others over 13 MYA. Additionally, distinct mitochondrial lineages are present within the taxa A. heteropholidotus and A. morazani, and importantly, samples of A. crassulus and A. sminthus previously used in major anole phylogenetic analyses are not recovered as conspecific with those taxa. We infer a Chortís Highland origin for the ancestor of this subgroup, and estimate cladogenesis of this subgroup began approximately 22 MYA.

CONCLUSIONS

Our results provide new insights into the evolution, biogeography, and timing of diversification of the Anolis crassulus species subgroup. The disjunctly distributed Anolis crassulus sensu lato represents several morphologically conserved, molecularly distinct anoles, and several other species in the subgroup contain multiple isolated lineages.

Demographic model selection using random forests and the site frequency spectrum

Phylogeographic data sets have grown from tens to thousands of loci in recent years, but extant statistical methods do not take full advantage of these large data sets. For example, approximate Bayesian computation (ABC) is a commonly used method for the explicit comparison of alternate demographic histories, but it is limited by the "curse of dimensionality" and issues related to the simulation and summarization of data when applied to next-generation sequencing (NGS) data sets. We implement here several improvements to overcome these difficulties. We use a Random Forest (RF) classifier for model selection to circumvent the curse of dimensionality and apply a binned representation of the multidimensional site frequency spectrum (mSFS) to address issues related to the simulation and summarization of large SNP data sets. We evaluate the performance of these improvements using simulation and find low overall error rates (~7%). We then apply the approach to data from Haplotrema vancouverense, a land snail endemic to the Pacific Northwest of North America. Fifteen demographic models were compared, and our results support a model of recent dispersal from coastal to inland rainforests. Our results demonstrate that binning is an effective strategy for the construction of a mSFS and imply that the statistical power of RF when applied to demographic model selection is at least comparable to traditional ABC algorithms. Importantly, by combining these strategies, large sets of models with differing numbers of populations can be evaluated.