Brian Tilston Smith

Biological evidence supports an early and complex emergence of the Isthmus of Panama

Significance The formation of the Isthmus of Panama, which linked North and South America, is key to understanding the biodiversity, oceanography, atmosphere, and climate in the region. Despite its importance across multiple disciplines, the timing of formation and emergence of the Isthmus and the biological patterns it created have been controversial. Here, we analyze molecular and fossil data, including terrestrial and marine organisms, to show that biotic migrations across the Isthmus of Panama began several million years earlier than commonly assumed. An earlier evolution of the Isthmus has broad implications for the mechanisms driving global climate (e.g., Pleistocene glaciations, thermohaline circulation) as well as the rich biodiversity of the Americas. The linking of North and South America by the Isthmus of Panama had major impacts on global climate, oceanic and atmospheric currents, and biodiversity, yet the timing of this critical event remains contentious. The Isthmus is traditionally understood to have fully closed by ca. 3.5 million years ago (Ma), and this date has been used as a benchmark for oceanographic, climatic, and evolutionary research, but recent evidence suggests a more complex geological formation. Here, we analyze both molecular and fossil data to evaluate the tempo of biotic exchange across the Americas in light of geological evidence. We demonstrate significant waves of dispersal of terrestrial organisms at approximately ca. 20 and 6 Ma and corresponding events separating marine organisms in the Atlantic and Pacific oceans at ca. 23 and 7 Ma. The direction of dispersal and their rates were symmetrical until the last ca. 6 Ma, when northern migration of South American lineages increased significantly. Variability among taxa in their timing of dispersal or vicariance across the Isthmus is not explained by the ecological factors tested in these analyses, including biome type, dispersal ability, and elevation preference. Migration was therefore not generally regulated by intrinsic traits but more likely reflects the presence of emergent terrain several millions of years earlier than commonly assumed. These results indicate that the dramatic biotic turnover associated with the Great American Biotic Interchange was a long and complex process that began as early as the Oligocene–Miocene transition.

The drivers of tropical speciation

Since the recognition that allopatric speciation can be induced by large-scale reconfigurations of the landscape that isolate formerly continuous populations, such as the separation of continents by plate tectonics, the uplift of mountains or the formation of large rivers, landscape change has been viewed as a primary driver of biological diversification. This process is referred to in biogeography as vicariance. In the most species-rich region of the world, the Neotropics, the sundering of populations associated with the Andean uplift is ascribed this principal role in speciation. An alternative model posits that rather than being directly linked to landscape change, allopatric speciation is initiated to a greater extent by dispersal events, with the principal drivers of speciation being organism-specific abilities to persist and disperse in the landscape. Landscape change is not a necessity for speciation in this model. Here we show that spatial and temporal patterns of genetic differentiation in Neotropical birds are highly discordant across lineages and are not reconcilable with a model linking speciation solely to landscape change. Instead, the strongest predictors of speciation are the amount of time a lineage has persisted in the landscape and the ability of birds to move through the landscape matrix. These results, augmented by the observation that most species-level diversity originated after episodes of major Andean uplift in the Neogene period, suggest that dispersal and differentiation on a matrix previously shaped by large-scale landscape events was a major driver of avian speciation in lowland Neotropical rainforests.

Target Capture and Massively Parallel Sequencing of Ultraconserved Elements for Comparative Studies at Shallow Evolutionary Time Scales

Comparative genetic studies of non-model organisms are transforming rapidly due to major advances in sequencing technology. A limiting factor in these studies has been the identification and screening of orthologous loci across an evolutionarily distant set of taxa. Here, we evaluate the efficacy of genomic markers targeting ultraconserved DNA elements (UCEs) for analyses at shallow evolutionary timescales. Using sequence capture and massively parallel sequencing to generate UCE data for five co-distributed Neotropical rainforest bird species, we recovered 776-1516 UCE loci across the five species. Across species, 53-77% of the loci were polymorphic, containing between 2.0 and 3.2 variable sites per polymorphic locus, on average. We performed species tree construction, coalescent modeling, and species delimitation, and we found that the five co-distributed species exhibited discordant phylogeographic histories. We also found that species trees and divergence times estimated from UCEs were similar to the parameters obtained from mtDNA. The species that inhabit the understory had older divergence times across barriers, contained a higher number of cryptic species, and exhibited larger effective population sizes relative to the species inhabiting the canopy. Because orthologous UCEs can be obtained from a wide array of taxa, are polymorphic at shallow evolutionary timescales, and can be generated rapidly at low cost, they are an effective genetic marker for studies investigating evolutionary patterns and processes at shallow timescales.

Sequence Capture Versus Restriction Site Associated DNA Sequencing for Shallow Systematics

Sequence capture and restriction site associated DNA sequencing (RAD-Seq) are two genomic enrichment strategies for applying next-generation sequencing technologies to systematics studies. At shallow timescales, such as within species, RAD-Seq has been widely adopted among researchers, although there has been little discussion of the potential limitations and benefits of RAD-Seq and sequence capture. We discuss a series of issues that may impact the utility of sequence capture and RAD-Seq data for shallow systematics in non-model species. We review prior studies that used both methods, and investigate differences between the methods by re-analyzing existing RAD-Seq and sequence capture data sets from a Neotropical bird (Xenops minutus). We suggest that the strengths of RAD-Seq data sets for shallow systematics are the wide dispersion of markers across the genome, the relative ease and cost of laboratory work, the deep coverage and read overlap at recovered loci, and the high overall information that results. Sequence captures benefits include flexibility and repeatability in the genomic regions targeted, success using low-quality samples, more straightforward read orthology assessment, and higher per-locus information content. The utility of a method in systematics, however, rests not only on its performance within a study, but on the comparability of data sets and inferences with those of prior work. In RAD-Seq data sets, comparability is compromised by low overlap of orthologous markers across species and the sensitivity of genetic diversity in a data set to an interaction between the level of natural heterozygosity in the samples examined and the parameters used for orthology assessment. In contrast, sequence capture of conserved genomic regions permits interrogation of the same loci across divergent species, which is preferable for maintaining comparability among data sets and studies for the purpose of drawing general conclusions about the impact of historical processes across biotas. We argue that sequence capture should be given greater attention as a method of obtaining data for studies in shallow systematics and comparative phylogeography.

As Old as the Hills: Montane Scorpions in Southwestern North America Reveal Ancient Associations between Biotic Diversification and Landscape History

Background The age of lineages has become a fundamental datum in studies exploring the interaction between geological transformation and biotic diversification. However, phylogeographical studies are often biased towards lineages that are younger than the geological features of the landscapes they inhabit. A temporally deeper historical biogeography framework may be required to address episodes of biotic diversification associated with geologically older landscape changes. Signatures of such associations may be retained in the genomes of ecologically specialized (stenotopic) taxa with limited vagility. In the study presented here, genetic data from montane scorpions in the Vaejovis vorhiesi group, restricted to humid rocky habitats in mountains across southwestern North America, were used to explore the relationship between scorpion diversification and regional geological history. Results Strong phylogeographical signal was evident within the vorhiesi group, with 27 geographically cohesive lineages inferred from a mitochondrial phylogeny. A time-calibrated multilocus species tree revealed a pattern of Miocene and Pliocene (the Neogene period) lineage diversification. An estimated 21 out of 26 cladogenetic events probably occurred prior to the onset of the Pleistocene, 2.6 million years ago. The best-fit density-dependent model suggested diversification rate in the vorhiesi group gradually decreased through time. Conclusions Scorpions of the vorhiesi group have had a long history in the highlands of southwestern North America. Diversification among these stenotopic scorpions appears to have occurred almost entirely within the Neogene period, and is temporally consistent with the dynamic geological history of the Basin and Range, and Colorado Plateau physiographical provinces. The persistence of separate lineages at small spatial scales suggests that a combination of ecological stenotopy and limited vagility may make these scorpions particularly valuable indicators of geomorphological evolution.

Identifying biases at different spatial and temporal scales of diversification: a case study in the Neotropical parrotlet genus Forpus

The temporal origins of the extraordinary biodiversity of the Neotropical region are highly debated. Recent empirical work has found support for alternative models on the tempo of speciation in Neotropical species further fuelling the debate. However, relationships within many Neotropical lineages are poorly understood, and it is unclear how this uncertainty impacts inferences on the evolution of taxa in the region. We examined the robustness of diversification patterns in the avian genus Forpus by testing whether the use of different units of biodiversity (i.e. biological species and statistically inferred species) impacted diversification rates and inferences regarding important biogeographic breaks in the genus. We found that the best‐fit model of diversification for the biological species data set was a declining rate of diversification; whereas a model of constant diversification was the best‐fit model for statistically inferred species or subspecies. Moreover, the relative importance of different landscape features in delimiting genetic structure across the landscape varied across data sets with differing units of biodiversity. Patterns based on divergence times among biological species indicated old speciation events across major geographic and river barriers. In contrast, data sets more inclusive of the diversity in Forpus illustrate the role of both old divergence across major landscape features and more recent divergences that are possibly attributed to Pleistocene climatic changes. Overall, these results indicate that conflicting models on the temporal origins of Neotropical birds may be attributable to sampling biases.

Evaluating the role of contracting and expanding rainforest in initiating cycles of speciation across the Isthmus of Panama

Climatic and geological changes across time are presumed to have shaped the rich biodiversity of tropical regions. However, the impact climatic drying and subsequent tropical rainforest contraction had on speciation has been controversial because of inconsistent palaeoecological and genetic data. Despite the strong interest in examining the role of climatic change on speciation in the Neotropics there has been few comparative studies, particularly, those that include non-rainforest taxa. We used bird species that inhabit humid or dry habitats that dispersed across the Panamanian Isthmus to characterize temporal and spatial patterns of speciation across this barrier. Here, we show that these two assemblages of birds exhibit temporally different speciation time patterns that supports multiple cycles of speciation. Evidence for these cycles is further corroborated by the finding that both assemblages consist of ‘young’ and ‘old’ species, despite dry habitat species pairs being geographically more distant than pairs of humid habitat species. The matrix of humid and dry habitats in the tropics not only allows for the maintenance of high species richness, but additionally this study suggests that these environments may have promoted speciation. We conclude that differentially expanding and contracting distributions of dry and humid habitats was probably an important contributor to speciation in the tropics.

Examining the role of effective population size on mitochondrial and multilocus divergence time discordance in a songbird.

Estimates of speciation times are subject to a number of potential errors. One source of bias is that effective population size (Ne) has been shown to influence substitution rates. This issue is of particular interest for phylogeographic studies because population sizes can vary dramatically among genetically structured populations across species’ ranges. In this study, we used multilocus data to examine temporal phylogeographic patterns in a widespread North American songbird, the Northern Cardinal (Cardinalis cardinalis). Species tree estimation indicated that the phylogeographic structure of C. cardinalis was comprised of four well-supported mainland lineages with large population sizes (large Ne) and two island lineages comprised of much smaller populations (small Ne). We inferred speciation times from mtDNA and multilocus data and found there was discordance between events that represented island-mainland divergences, whereas both estimates were similar for divergences among mainland lineages. We performed coalescent simulations and found that the difference in speciation times could be attributed to stochasticity for a recently diverged island lineage. However, the magnitude of the change between speciation times estimated from mtDNA and multilocus data of an older island lineage was substantially greater than predicted by coalescent simulations. For this divergence, we found the discordance in time estimates was due to a substantial increase in the mtDNA substitution rate in the small island population. These findings indicate that in phylogeographic studies the relative tempo of evolution between mtDNA and nuclear DNA can become highly discordant in small populations.

Reply to Lessios and Marko et al.: Early and progressive migration across the Isthmus of Panama is robust to missing data and biases

The emergence of the Isthmus of Panama left a major imprint on the biodiversity of the Americas. The connection between South and North America facilitated dispersal of terrestrial and freshwater organisms, while separating marine species between the eastern Pacific and Caribbean seas. Recent geological data have questioned the long-standing view of a Pliocene emergence of the Isthmus (1) and show that the Central American Seaway, defined as the deep oceanic seaway along the tectonic boundary of the South American plate and Panama arc, was already closed by 15–13 Ma (2). Caribbean–Pacific shallow water exchange probably continued, albeit intermittently, until a full closure at 3.5 Ma (1–3). Recently Bacon et al. (3) used molecular and fossil data to evaluate the timing, tempo, and directionality of biotic exchange and vicariance across the Isthmus, and tested whether biological data are congruent with recent geological evidence. Significant increases in terrestrial dispersals were found at ca. 20 and 6 Ma, and increases in marine vicariance at ca. 23 and 7 Ma. Similar patterns prevailed despite intrinsic differences among the taxonomic groups surveyed. This led Bacon et al. (3) to reject the assumption of a single closure of the Isthmus at ca. 3.5 Ma in favor of an older, more complex model of land emergence and biotic interchange.

A latitudinal phylogeographic diversity gradient in birds

High tropical species diversity is often attributed to evolutionary dynamics over long timescales. It is possible, however, that latitudinal variation in diversification begins when divergence occurs within species. Phylogeographic data capture this initial stage of diversification in which populations become geographically isolated and begin to differentiate genetically. There is limited understanding of the broader implications of intraspecific diversification because comparative analyses have focused on species inhabiting and evolving in restricted regions and environments. Here, we scale comparative phylogeography up to the hemisphere level and examine whether the processes driving latitudinal differences in species diversity are also evident within species. We collected genetic data for 210 New World bird species distributed across a broad latitudinal gradient and estimated a suite of metrics characterizing phylogeographic history. We found that lower latitude species had, on average, greater phylogeographic diversity than higher latitude species and that intraspecific diversity showed evidence of greater persistence in the tropics. Factors associated with species ecologies, life histories, and habitats explained little of the variation in phylogeographic structure across the latitudinal gradient. Our results suggest that the latitudinal gradient in species richness originates, at least partly, from population-level processes within species and are consistent with hypotheses implicating age and environmental stability in the formation of diversity gradients. Comparative phylogeographic analyses scaled up to large geographic regions and hundreds of species can show connections between population-level processes and broad-scale species-richness patterns.