Search

Ever since Alfred Russel Wallace's nineteenth-century observation that related terrestrial species are often separated on opposing riverbanks, major Amazonian rivers have been recognized as key drivers of speciation. However, rivers are dynamic entities whose widths and courses may vary through time. It thus remains unknown how effective rivers are at reducing gene flow and promoting speciation over long timescales. We fit demographic models to genomic sequences to reconstruct the history of gene flow in three pairs of avian taxa fully separated by different Amazonian rivers, and whose geographic ranges do not make contact in headwater regions where rivers may cease to be barriers. Models with gene flow were best fit but still supported an initial period without any gene flow, which ranged from 187 000 to over 959 000 years, suggesting that rivers are capable of initiating speciation through long stretches of allopatric divergence. Allopatry was followed by either bursts or prolonged episodes of gene flow that retarded genomic differentiation but did not fully homogenize populations. Our results support Amazonian rivers as key barriers that promoted speciation and the build-up of species richness, but they also suggest that river barriers are often leaky, with genomic divergence accumulating slowly owing to episodes of substantial gene flow.

The "Amazon tipping point" is a global change scenario resulting in replacement of upland terra-firme forests by large-scale "savannization" of mostly southern and eastern Amazon. Reduced rainfall accompanying the Last Glacial Maximum (LGM) has been proposed to have acted as such a tipping point in the past, with the prediction that terra-firme inhabiting species should have experienced reductions in population size as drier habitats expanded. Here, we use whole-genomes of an Amazonian endemic organism (Scale-backed antbirds - Willisornis spp.) sampled from nine populations across the region to test this historical demography scenario. Populations from southeastern Amazonia and close to the Amazon- Cerrado ecotone exhibited a wide range of demographic patterns, while most of those from northern and western Amazonia experienced uniform expansions between 400 kya and 80-60 kya, with gradual declines toward 20 kya. Southeastern populations of Willisornis were the last to diversify and showed smaller heterozygosity and higher runs of homozygosity values than western and northern populations. These patterns support historical population declines throughout the Amazon that affected more strongly lineages in the southern and eastern areas, where historical "tipping point" conditions existed due to the widespread replacement of humid forest by drier and open vegetation during the LGM., Competing Interests: The authors declare that the research was conducted without any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2024 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

The skuas and jaegers (Stercorariidae) are an enigmatic family of seven seabird species that breed at Arctic and Antarctic latitudes. The phylogenetic relationships amongst the species have been controversial, with one of the biggest enigmas involving the Pomarine Jaeger (Stercorarius pomarinus), which has been proposed to represent a hybrid species originating from the merging of distant lineages within the complex. We inferred a phylogeny for the family using multispecies coalescent methods with whole-genome sequencing for all seven species of Stercorariidae, and document an evolutionary history rich in introgression. We uncover evidence for mitochondrial capture and nuclear introgression between S. pomarinus and Stercorarius skua, providing a potential avenue for adaptive introgression. One candidate for adaptive introgression is the MC1R plumage gene which appears to have introgressed from one of the large skuas into S. pomarinus, where it now forms the basis of the dark-morph color polymorphism of that species. We further highlight a complex biogeographical history of interchange between the Arctic and Antarctic, with unexpected close ancestry between S. skua of the northern hemisphere and Stercorarius antarcticus of the southern hemisphere. These results highlight the dynamic history of introgression during pelagic seabird radiation. [Incomplete lineage sorting; introgression; mitochondrial capture; phylogenomics; skua; species tree; stercorariidae; whole-genome resequencing.]., (© The Author(s) 2022. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

AbstractResearch over the past three decades has shown that ecology-based extrinsic reproductive barriers can rapidly arise to generate incipient species-but such barriers can also rapidly dissolve when environments change, resulting in incipient species collapse. Understanding the evolution of unconditional, "intrinsic" reproductive barriers is therefore important for understanding the longer-term buildup of biodiversity. In this article, we consider ecology's role in the evolution of intrinsic reproductive isolation. We suggest that this topic has fallen into a gap between disciplines: while evolutionary ecologists have traditionally focused on the rapid evolution of extrinsic isolation between co-occurring ecotypes, speciation geneticists studying intrinsic isolation in other taxa have devoted little attention to the ecological context in which it evolves. We argue that for evolutionary ecology to close this gap, the field will have to expand its focus beyond rapid adaptation and its traditional model systems. Synthesizing data from several subfields, we present circumstantial evidence for and against different forms of ecological adaptation as promoters of intrinsic isolation and discuss alternative forces that may be significant. We conclude by outlining complementary approaches that can better address the role of ecology in the evolution of nonephemeral reproductive barriers and, by extension, less ephemeral species.

After decades of debate, biologists today largely agree that most speciation events require an allopatric phase (that is, geographic separation), but the role of adaptive ecological divergence during this critical period is still unknown. Here, we show that relatively few allopatric pairs of birds, mammals, or amphibians exhibit trait differences consistent with models of divergent adaptation in each of many ecologically relevant traits. By fitting new evolutionary models to numerous sets of sister-pair trait differences, we find that speciating and recently speciated allopatric taxa seem to overwhelmingly evolve under similar rather than divergent macro-selective pressures. This contradicts the classical view of divergent adaptation as a prominent driver of the early stages of speciation and helps synthesize two historical controversies regarding the ecology and geography of species formation.

Supervolcanoes are volcanoes capable of mega-colossal eruptions that emit more than 1,000 km 3 of ash and other particles. 1 The earth's most recent mega-colossal eruption was the Oruanui eruption of the Taupo supervolcano 25,580 years before present (YBP) on the central North Island of New Zealand. 2 This eruption blanketed major swaths of the North Island in thick layers of ash and igneous rock, 2 , 3 devastating habitats and likely causing widespread population extinctions. 4-7 An additional devastating super-colossal eruption (>100 km 3 ) of the Taupo supervolcano occurred approximately 1,690 YBP. 8 The impacts of such massive but ephemeral natural disasters on contemporary population genetic structure remain underexplored. Here, we combined data for 4,951 SNPs with spatially explicit demographic and coalescent models within an approximate Bayesian computation framework to test the drivers of genetic structure in brown kiwi (Apteryx mantelli). Our results strongly support the importance of eruptions of the Taupo supervolcano in restructuring pre-existing geographic patterns of population differentiation and genetic diversity. Range shifts due to climatic oscillations-a frequent explanation for genetic structure 9 -are insufficient to fully explain the empirical data. Meanwhile, recent range contraction and fragmentation due to historically documented anthropogenic habitat alteration adds no explanatory power to our models. Our results support a major role for cycles of destruction and post-volcanic recolonization in restructuring the population genomic landscape of brown kiwi and highlight how ancient and ephemeral mega-disasters may leave a lasting legacy on patterns of intraspecific genetic variation., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Phylogeographical studies of the most species-rich region of the planet-the Amazon basin-have repeatedly uncovered genetically distinctive, allopatric lineages within currently named species, but understanding whether such lineages are reproductively isolated species is challenging. Here we harness the power of genome-wide data sets together with detailed phylogeographical sampling to both characterize the number of unique lineages and infer levels of reproductive isolation for three parapatric manakin species that make up the genus Pipra. The mitochondrial and nuclear genomes both support six distinctive lineages. The youngest lineages are now highly admixed with each other across major portions of their geographical ranges with one lineage now extinct in a genomically unadmixed state. In contrast, the oldest sets of lineages-dated to 1.4 million years-exhibit narrow hybrid zones. By fitting demographic models to parapatric lineage pairs we found that levels of gene flow and genomic homogenization decline with increasing evolutionary age. Only lineages descending from the basal node at 1.4 million years ago in the genus experience negligible gene flow, possess genomes resistant to homogenization and are separated by narrow hybrid zones. We conclude that a million years or more were required for Pipra manakins to become reproductively isolated. We suggest the six lineages be reclassified as two or three reproductively isolated species. Our unique approach to quantifying reproductive isolation in parapatric lineages could be applied broadly to other phylogeographical studies and would help determine species classification of the plethora of newly identified lineages in the Amazon basin and other regions., (© 2022 John Wiley & Sons Ltd.)

Small and fragmented populations may become rapidly differentiated due to genetic drift, making it difficult to distinguish whether neutral genetic structure is a signature of recent demographic events, or of long-term evolutionary processes that could have allowed populations to adaptively diverge. We sequenced 52 whole genomes to examine Holocene demographic history and patterns of adaptation in kiwi ( Apteryx ), and recovered 11 strongly differentiated genetic clusters corresponding to previously recognized lineages. Demographic models suggest that all 11 lineages experienced dramatic population crashes relative to early- or mid-Holocene levels. Small population size is associated with low genetic diversity and elevated genetic differentiation ( F ST ), suggesting that population declines have strengthened genetic structure and led to the loss of genetic diversity. However, population size is not correlated with inbreeding rates. Eight lineages show signatures of lineage-specific selective sweeps (284 sweeps total) that are unlikely to have been caused by demographic stochasticity. Overall, these results suggest that despite strong genetic drift associated with recent bottlenecks, most kiwi lineages possess unique adaptations and should be recognized as separate adaptive units in conservation contexts. Our work highlights how whole-genome datasets can address longstanding uncertainty about the evolutionary and conservation significance of small and fragmented populations of threatened species.

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species., (© 2021 John Wiley & Sons Ltd.)

Geographic contact between sister lineages often occurs near the final stages of speciation, but its role in speciation's completion remains debated. Reproductive isolation may be essentially complete prior to secondary contact. Alternatively, costly interactions between partially reproductively isolated species - such as maladaptive hybridization or competition for resources - may select for divergence, increasing reproductive isolation and driving speciation toward completion. Here, we use coalescent demographic modelling and whole-genome data sets to show that a period of contact and elevated hybridization between sympatric eastern North American populations of two cryptic bird species preceded a major increase in reproductive isolation between these populations within the last 10,000 years. In contrast, substantial introgression continues to the present in a western contact zone where geographic overlap is much narrower and probably of more recent origin. In the sympatric eastern region where reproductive isolation has increased, it is not accompanied by character displacement in key morphometric traits, plumage coloration, or ecological traits. While the precise trait and underlying mechanism driving increased reproductive isolation remains unknown, we discuss several possibilities and outline avenues for future research. Overall, our results highlight how demographic models can reveal the geographic context in which reproductive isolation was completed, and demonstrate how contact can accelerate the final stages of speciation., (© 2021 John Wiley & Sons Ltd.)

Coexisting (sympatric) pairs of closely related species are often characterized by exaggerated trait differences. This widespread pattern is consistent with adaptation for reduced similarity due to costly interactions (i.e., "character displacement")-a classic hypothesis in evolutionary theory. But it is equally consistent with a community assembly bias in which lineages with greater trait differences are more likely to establish overlapping ranges in the first place (i.e., "species sorting"), as well as with null expectations of trait divergence through time. Few comparative analyses have explicitly modeled these alternatives, and it remains unclear whether trait divergence is a general prerequisite for sympatry or a consequence of interactions between sympatric species. Here, we develop statistical models that allow us to distinguish the signature of these processes based on patterns of trait divergence in closely related lineage pairs. We compare support for each model using a dataset of bill shape differences in 207 pairs of New World terrestrial birds representing 30 avian families. We find that character displacement models are overwhelmingly supported over species sorting and null expectations, indicating that exaggerated bill shape differences in sympatric pairs result from enhanced divergent selection in sympatry. We additionally detect a latitudinal gradient in character displacement, which appears strongest in the tropics. Our analysis implicates costly species interactions as powerful drivers of trait divergence in a major vertebrate fauna. These results help substantiate a long-standing but equivocally supported linchpin of evolutionary theory., Competing Interests: The authors declare no competing interest.

Collisions with buildings cause up to 1 billion bird fatalities annually in the United States and Canada. However, efforts to reduce collisions would benefit from studies conducted at large spatial scales across multiple study sites with standardized methods and consideration of species- and life-history-related variation and correlates of collisions. We addressed these research needs through coordinated collection of data on bird collisions with buildings at sites in the United States (35), Canada (3), and Mexico (2). We collected all carcasses and identified species. After removing records for unidentified carcasses, species lacking distribution-wide population estimates, and species with distributions overlapping fewer than 10 sites, we retained 269 carcasses of 64 species for analysis. We estimated collision vulnerability for 40 bird species with ≥2 fatalities based on their North American population abundance, distribution overlap in study sites, and sampling effort. Of 10 species we identified as most vulnerable to collisions, some have been identified previously (e.g., Black-throated Blue Warbler [Setophaga caerulescens]), whereas others emerged for the first time (e.g., White-breasted Nuthatch [Sitta carolinensis]), possibly because we used a more standardized sampling approach than past studies. Building size and glass area were positively associated with number of collisions for 5 of 8 species with enough observations to analyze independently. Vegetation around buildings influenced collisions for only 1 of those 8 species (Swainson's Thrush [Catharus ustulatus]). Life history predicted collisions; numbers of collisions were greatest for migratory, insectivorous, and woodland-inhabiting species. Our results provide new insight into the species most vulnerable to building collisions, making them potentially in greatest need of conservation attention to reduce collisions and into species- and life-history-related variation and correlates of building collisions, information that can help refine collision management., (© 2020 Society for Conservation Biology.)

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity 1-4 . Sparse taxon sampling has previously been proposed to confound phylogenetic inference 5 , and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

How species evolve reproductive isolation in the species-rich Amazon basin is poorly understood in vertebrates. Here, we sequenced a reference genome and used a genome-wide sample of SNPs to analyze a hybrid zone between two highly cryptic species of Hypocnemis warbling-antbirds-the Rondonia warbling-antbird (H. ochrogyna) and Spix's warbling-antbird (H. striata)-in a headwater region of southern Amazonia. We found that both species commonly hybridize, producing F 1 s and a variety of backcrosses with each species but we detected only one F 2 -like hybrid. Patterns of heterozygosity, hybrid index, and interchromosomal linkage disequilibrium in hybrid populations closely match expectations under strong postzygotic isolation. Hybrid zone width (15.4 km) was much narrower than expected (211 km) indicating strong selection against hybrids. A remarkably high degree of concordance in cline centers and widths across loci, and a lack of reduced interspecific F st between populations close to versus far from the contact zone, suggest that genetic incompatibilities have rendered most of the genome immune to introgression. These results support intrinsic postzygotic isolation as a driver of speciation in a moderately young cryptic species pair from the Amazon and suggest that species richness of the Amazon may be grossly underestimated., (© 2020 The Authors. Evolution © 2020 The Society for the Study of Evolution.)

Antbirds (Thamnophilidae) are a large neotropical family of passerine bird renowned for the ant-following foraging strategies of several members of this clade. The high diversity of antbirds provides ample opportunity for speciation studies, however these studies can be hindered by the lack of an annotated antbird reference genome. In this study, we produced a high-quality annotated reference genome for the Xingu Scale-backed Antbird (Willisornis vidua nigrigula) using 10X Genomics Chromium linked-reads technology. The assembly is 1.09 Gb, with a scaffold N50 of 12.1 Mb and 17,475 annotated protein coding genes. We compare the proteome of W. v. nigrigula to several other passerines, and produce annotations for two additional antbird genomes in order to identify genes under lineage-specific positive selection and gene families with evidence for significant expansions in antbirds. Several of these genes have functions potentially related to the lineage-specific traits of antbirds, including adaptations for thermoregulation in a humid tropical environment., (Copyright © 2020 Elsevier Inc. All rights reserved.)

AbstractEcological differentiation between lineages is widely considered to be an important driver of speciation, but support for this hypothesis is mainly derived from the detailed study of a select set of model species pairs. Mounting evidence from nonmodel taxa, meanwhile, suggests that speciation often occurs with minimal differentiation in ecology or ecomorphology, calling into question the true contribution of divergent adaptation to species richness in nature. To better understand divergent ecological adaptation and its role in speciation generally, researchers require a comparative approach that can distinguish its signature from alternative processes, such as drift and parallel selection, in data sets containing many species pairs. Here we introduce new statistical models of divergent adaptation in the continuous traits of paired lineages. In these models, ecomorphological characters diverge as two lineages adapt toward alternative phenotypic optima following their departure from a common ancestor. The absolute distance between optima measures the extent of divergent selection and provides a basis for interpretation. We encode the models in the new R package diverge and extend them to allow the distance between optima to vary across continuous and categorical variables. We test model performance using simulation and demonstrate model application using published data sets of trait divergence in birds and mammals. Our framework provides the first explicit test for signatures of divergent selection in trait divergence data sets, and it will enable empiricists from a wide range of fields to better understand the dynamics of divergent adaptation and its prevalence in nature beyond just our best-studied model systems.

The incidence of introgression during the diversification process and the timespan following divergence when introgression is possible are poorly understood in the neotropics where high species richness could provide extensive opportunities for genetic exchange. We used thousands of genome-wide SNPs to infer phylogenetic relationships, calculate ages of splitting, and to estimate the timing of introgression in a widespread avian neotropical genus of woodcreepers. Five distinct introgression events were reconstructed involving taxa classified both as subspecies and species including lineages descending from the basal-most split, dated to 7.3 million years ago. Introgression occurred between just a few hundred thousand to about 2.5 million years following divergence, suggesting substantial portions of the genome are capable of introgressing across taxa boundaries during a protracted time window of a few million years following divergence. Despite this protracted time window, we found that the proportion of the genome introgressing (6-11%) declines with the time of introgression following divergence, suggesting that the genome becomes progressively more immune to introgression as reproductive isolation increases., (© 2019 The Authors. Evolution © 2019 The Society for the Study of Evolution.)

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post-glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water-flow drive the downstream accumulation of genetic diversity, whereas post-glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls - a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls' role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome-wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post-glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species' genetic diversity and highlight the need to better understand post-glacial dispersal pathways., (© 2020 John Wiley & Sons Ltd.)

Genetic data indicate differences in speciation rate across latitudes, but underlying causes have been difficult to assess because a critical phase of the speciation process is initiated in allopatry, in which, by definition, individuals from different taxa do not interact. We conducted song playback experiments between 109 related pairs of mostly allopatric bird species or subspecies in Amazonia and North America to compare the rate of evolution of male discrimination of songs. Relative to local controls, the number of flyovers and approach to the speaker were higher in Amazonia. We estimate that responses to songs of relatives are being lost about 6 times more slowly in Amazonia than in North America. The slow loss of response holds even after accounting for differences in song frequency and song length. Amazonian species with year-round territories are losing aggressive responses especially slowly. We suggest the presence of many species and extensive interspecific territoriality favors recognition of songs sung by sympatric heterospecifics, which results in a broader window of recognition and hence an ongoing response to novel similar songs. These aggressive responses should slow the establishment of sympatry between recently diverged forms. If male responses to novel allopatric taxa reflect female responses, then premating reproductive isolation is also evolving more slowly in Amazonia. The findings are consistent with previously demonstrated slower recent rates of expansion of sister taxa into sympatry, slower rates of evolution of traits important for premating isolation, and slower rates of speciation in general in Amazonia than in temperate North America., Competing Interests: The authors have declared that no competing interests exist

We possess limited understanding of how speciation unfolds in the most species-rich region of the planet-the Amazon basin. Hybrid zones provide valuable information on the evolution of reproductive isolation, but few studies of Amazonian vertebrate hybrid zones have rigorously examined the genome-wide underpinnings of reproductive isolation. We used genome-wide genetic datasets to show that two deeply diverged, but morphologically cryptic sister species of forest understorey birds show little evidence for prezygotic reproductive isolation, but substantial postzygotic isolation. Patterns of heterozygosity and hybrid index revealed that hybrid classes with heavily recombined genomes are rare and closely match simulations with high levels of selection against hybrids. Genomic and geographical clines exhibit a remarkable similarity across loci in cline centres, and have exceptionally narrow cline widths, suggesting that postzygotic isolation is driven by genetic incompatibilities at many loci, rather than a few loci of strong effect. We propose Amazonian understorey forest birds speciate slowly via gradual accumulation of postzygotic genetic incompatibilities, with prezygotic barriers playing a less important role. Our results suggest old, cryptic Amazonian taxa classified as subspecies could have substantial postzygotic isolation deserving species recognition and that species richness is likely to be substantially underestimated in Amazonia., (© 2018 The Author(s).)

Although some taxa are increasing in number due to active management and predator control, the overall number of kiwi ( Apteryx spp.) is declining. Kiwi are cryptic and rare, meaning current monitoring tools, such as call counts, radio telemetry, and surveys using detection dogs are labor-intensive, yield small datasets, and require substantial resources or provide inaccurate estimates of population sizes. A noninvasive genetic approach could help the conservation effort. We optimized a panel of 23 genetic markers (22 autosomal microsatellite loci and an allosomal marker) to discriminate between all species of kiwi and major lineages within species, while simultaneously determining sex. Markers successfully amplified from both fecal and shed feather DNA samples collected in captivity. We found that DNA extraction was more efficient from shed feathers, but DNA quality was greater with feces, although this was sampling dependent. Our microsatellite panel was able to distinguish between contemporary kiwi populations and lineages and provided PI values in the range of 4.3 × 10 -5 to 2.0 × 10 -19 , which in some cases were sufficient for individualization and mark-recapture studies. As such, we have tested a wide-reaching, noninvasive molecular approach that will improve conservation management by providing better parameter estimates associated with population ecology and demographics such as abundance, growth rates, and genetic diversity.

Hybrid speciation is rare in vertebrates, and reproductive isolation arising from hybridization is infrequently demonstrated. Here, we present evidence supporting a hybrid-speciation event involving the genetic admixture of the snow-capped ( Lepidothrix nattereri ) and opal-crowned ( Lepidothrix iris ) manakins of the Amazon basin, leading to the formation of the hybrid species, the golden-crowned manakin ( Lepidothrix vilasboasi ). We used a genome-wide SNP dataset together with analysis of admixture, population structure, and coalescent modeling to demonstrate that the golden-crowned manakin is genetically an admixture of these species and does not represent a hybrid zone but instead formed through ancient genetic admixture. We used spectrophotometry to quantify the coloration of the species-specific male crown patches. Crown patches are highly reflective white (snow-capped manakin) or iridescent whitish-blue to pink (opal-crowned manakin) in parental species but are a much less reflective yellow in the hybrid species. The brilliant coloration of the parental species results from nanostructural organization of the keratin matrix feather barbs of the crown. However, using electron microscopy, we demonstrate that the structural organization of this matrix is different in the two parental species and that the hybrid species is intermediate. The intermediate nature of the crown barbs, resulting from past admixture appears to have rendered a duller structural coloration. To compensate for reduced brightness, selection apparently resulted in extensive thickening of the carotenoid-laden barb cortex, producing the yellow crown coloration. The evolution of this unique crown-color signal likely culminated in premating isolation of the hybrid species from both parental species., Competing Interests: The authors declare no conflict of interest.

Molecular dating largely overturned the paradigm that global cooling during recent Pleistocene glacial cycles resulted in a burst of species diversification although some evidence exists that speciation was commonly promoted in habitats near the expanding and retracting ice sheets. Here, we used a genome-wide dataset of more than half a million base pairs of DNA to test for a glacially induced burst of diversification in kiwi, an avian family distributed within several hundred kilometers of the expanding and retracting glaciers of the Southern Alps of New Zealand. By sampling across the geographic range of the five kiwi species, we discovered many cryptic lineages, bringing the total number of kiwi taxa that currently exist to 11 and the number that existed just before human arrival to 16 or 17. We found that 80% of kiwi diversification events date to the major glacial advances of the Middle and Late Pleistocene. During this period, New Zealand was repeatedly fragmented by glaciers into a series of refugia, with the tiny geographic ranges of many kiwi lineages currently distributed in areas adjacent to these refugia. Estimates of effective population size through time show a dramatic bottleneck during the last glacial cycle in all but one kiwi lineage, as expected if kiwi were isolated in glacially induced refugia. Our results support a fivefold increase in diversification rates during key glacial periods, comparable with levels observed in classic adaptive radiations, and confirm that at least some lineages distributed near glaciated regions underwent rapid ice age diversification., Competing Interests: The authors declare no conflict of interest.

The role of historical factors in driving latitudinal diversity gradients is poorly understood. Here, we used an updated global phylogeny of terrestrial birds to test the role of three key historical factors-speciation, extinction, and dispersal rates-in generating latitudinal diversity gradients for eight major clades. We fit a model that allows speciation, extinction, and dispersal rates to differ, both with latitude and between the New and Old World. Our results consistently support extinction (all clades had lowest extinction where species richness was highest) as a key driver of species richness gradients across each of eight major clades. In contrast, speciation and dispersal rates showed no consistent latitudinal patterns across replicate bird clades, and thus are unlikely to represent general underlying drivers of latitudinal diversity gradients., (© 2016 The Author(s). Evolution © 2016 The Society for the Study of Evolution.)

Evolutionary biologists since Darwin have hypothesized that closely related species compete more intensely and are therefore less likely to coexist. However, recent theory posits that species diverge in two ways: either through the evolution of 'stabilizing differences' that promote coexistence by causing individuals to compete more strongly with conspecifics than individuals of other species, or through the evolution of 'fitness differences' that cause species to differ in competitive ability and lead to exclusion of the weaker competitor. We tested macroevolutionary patterns of divergence by competing pairs of annual plant species that differ in their phylogenetic relationships, and in whether they have historically occurred in the same region or different regions (sympatric versus allopatric occurrence). For sympatrically occurring species pairs, stabilizing differences rapidly increased with phylogenetic distance. However, fitness differences also increased with phylogenetic distance, resulting in coexistence outcomes that were unpredictable based on phylogenetic relationships. For allopatric species, stabilizing differences showed no trend with phylogenetic distance, whereas fitness differences increased, causing coexistence to become less likely among distant relatives. Our results illustrate the role of species' historical interactions in shaping how phylogenetic relationships structure competitive dynamics, and offer an explanation for the evolution of invasion potential of non-native species., (© 2016 The Author(s).)

Many understory birds and other groups form genetically differentiated subspecies or closely related species on opposite sides of major rivers of Amazonia, but are proposed to come into geographic contact in headwater regions where narrower river widths may present less of a dispersal barrier. Whether such forms hybridize in headwater regions is generally unknown, but has important implications to our understanding of the role of rivers as drivers of speciation. We used a dataset of several thousand single nucleotide polymorphisms to show that seven taxon pairs that differentiate across a major Amazonian river come into geographic contact and hybridize in headwater regions. All taxon pairs possessed hybrids with low numbers of loci in which alleles were inherited from both parental species, suggesting they are backcrossed with parentals, and indicating gene flow between parental populations. Ongoing gene flow challenges rivers as the sole cause of in situ speciation, but is compatible with the view that the wide river courses in the heart of Amazonia may have driven interfluvial divergence during episodes of wet forest retraction away from headwater regions. Taxa as old as 4 Ma in our Amazonian dataset continue to hybridize at contact zones, suggesting reproductive isolation evolves at a slow pace., (© 2015 The Author(s). Evolution © 2015 The Society for the Study of Evolution.)

Rapid diversification of sexual traits is frequently attributed to sexual selection, though explicit tests of this hypothesis remain limited. Spermatozoa exhibit remarkable variability in size and shape, and studies report a correlation between sperm morphology (sperm length and shape) and sperm competition risk or female reproductive tract morphology. However, whether postcopulatory processes (e.g., sperm competition and cryptic female choice) influence the speed of evolutionary diversification in sperm form is unknown. Using passerine birds, we quantified evolutionary rates of sperm length divergence among lineages (i.e., species pairs) and determined whether these rates varied with the level of sperm competition (estimated as relative testes mass). We found that relative testes mass was significantly and positively associated with more rapid phenotypic divergence in sperm midpiece and flagellum lengths, as well as total sperm length. In contrast, there was no association between relative testes mass and rates of evolutionary divergence in sperm head size, and models suggested that head length is evolutionarily constrained. Our results are the first to show an association between the strength of sperm competition and the speed of sperm evolution, and suggest that postcopulatory sexual selection promotes rapid evolutionary diversification of sperm morphology., (© 2015 The Author(s).)

Application of high-throughput sequencing platforms in the field of ecology and evolutionary biology is developing quickly with the introduction of efficient methods to reduce genome complexity. Numerous approaches for genome complexity reduction have been developed using different combinations of restriction enzymes, library construction strategies and fragment size selection. As a result, the choice of which techniques to use may become cumbersome, because it is difficult to anticipate the number of loci resulting from each method. We developed SimRAD, an R package that performs in silico restriction enzyme digests and fragment size selection as implemented in most restriction site associated DNA polymorphism and genotyping by sequencing methods. In silico digestion is performed on a reference genome or on a randomly generated DNA sequence when no reference genome sequence is available. SimRAD accurately predicts the number of loci under alternative protocols when a reference genome sequence is available for the targeted species (or a close relative) but may be unreliable when no reference genome is available. SimRAD is also useful for fine-tuning a given protocol to adjust the number of targeted loci. Here, we outline the functionality of SimRAD and provide an illustrative example of the use of the package (available on the CRAN at http://cran.r-project.org/web/packages/SimRAD)., (© 2014 John Wiley & Sons Ltd.)

Despite the importance of divergent selection to the speed of evolution, it remains poorly understood if divergent selection is more prevalent in the tropics (where species richness is highest), or at high latitudes (where paleoclimate change has been most intense). We tested whether the rate of climatic-niche evolution - one proxy for divergent selection - varies with latitude for 111 pairs of bird species. Using Brownian motion and Ornsetin-Ulhenbeck models, we show that evolutionary rates along two important axes of the climatic-niche - temperature and seasonality - have been faster at higher latitudes. We then tested whether divergence of the climatic-niche was associated with evolution in traits important in ecological differentiation (body mass) and reproductive isolation (song), and found that climatic divergence is associated with faster rates in both measures. These results highlight the importance of climate-mediated divergent selection pressures in driving evolutionary divergence and reproductive isolation at high latitudes., (© 2014 John Wiley & Sons Ltd/CNRS.)

Are rates of evolution and speciation fastest where diversity is greatest - the tropics? A commonly accepted theory links the latitudinal diversity gradient to a speciation pump model whereby the tropics produce species at a faster rate than extra-tropical regions. In this issue of Molecular Ecology, Botero et al. () test the speciation pump model using subspecies richness patterns for more than 9000 species of birds and mammals as a proxy for incipient speciation opportunity. Rather than using latitudinal centroids, the authors investigate the role of various environmental correlates of latitude as drivers of subspecies richness. Their key finding points to environmental harshness as a positive predictor of subspecies richness. The authors link high subspecies richness in environmental harsh areas to increased opportunities for geographic range fragmentation and/or faster rates of trait evolution as drivers of incipient speciation. Because environmental harshness generally increases with latitude, these results suggest that opportunity for incipient speciation is lowest where species richness is highest. The authors interpret this finding as incompatible with the view of the tropics as a cradle of diversity. Their results are consistent with a growing body of evidence that reproductive isolation and speciation occur fastest at high latitudes., (© 2013 John Wiley & Sons Ltd.)

Sexual selection is proposed to be an important driver of diversification in animal systems, yet previous tests of this hypothesis have produced mixed results and the mechanisms involved remain unclear. Here, we use a novel phylogenetic approach to assess the influence of sexual selection on patterns of evolutionary change during 84 recent speciation events across 23 passerine bird families. We show that elevated levels of sexual selection are associated with more rapid phenotypic divergence between related lineages, and that this effect is restricted to male plumage traits proposed to function in mate choice and species recognition. Conversely, we found no evidence that sexual selection promoted divergence in female plumage traits, or in male traits related to foraging and locomotion. These results provide strong evidence that female choice and male-male competition are dominant mechanisms driving divergence during speciation in birds, potentially linking sexual selection to the accelerated evolution of pre-mating reproductive isolation.

Freshwater habitats make up only ∼0.01% of available aquatic habitat and yet harbor 40% of all fish species, whereas marine habitats comprise >99% of available aquatic habitat and have only 60% of fish species. One possible explanation for this pattern is that diversification rates are higher in freshwater habitats than in marine habitats. We investigated diversification in marine and freshwater lineages in the New World silverside fish clade Menidiinae (Teleostei, Atherinopsidae). Using a time-calibrated phylogeny and a state-dependent speciation-extinction framework, we determined the frequency and timing of habitat transitions in Menidiinae and tested for differences in diversification parameters between marine and freshwater lineages. We found that Menidiinae is an ancestrally marine lineage that independently colonized freshwater habitats four times followed by three reversals to the marine environment. Our state-dependent diversification analyses showed that freshwater lineages have higher speciation and extinction rates than marine lineages. Net diversification rates were higher (but not significant) in freshwater than marine environments. The marine lineage-through time (LTT) plot shows constant accumulation, suggesting that ecological limits to clade growth have not slowed diversification in marine lineages. Freshwater lineages exhibited an upturn near the recent in their LTT plot, which is consistent with our estimates of high background extinction rates. All sequence data are currently being archived on Genbank and phylogenetic trees archived on Treebase., (© 2013 The Author(s). Evolution © 2013 The Society for the Study of Evolution.)

Through the course of an adaptive radiation, the evolutionary speed of cladogenesis and ecologically relevant trait evolution are expected to slow as species diversity increases, niches become occupied, and ecological opportunity declines. We develop new likelihood-based models to test diversity-dependent evolution in the auks, one of only a few families of seabirds adapted to underwater "flight," and which exhibit a large variety of bill sizes and shapes. Consistent with the expectations of adaptive radiation, we find both a decline in rates of cladogenesis (a sixfold decline) and bill shape (a 64-fold decline) evolution as diversity increased. Bill shape diverged into two clades at the basal cladogenesis event with one clade possessing mostly long, narrow bills used to forage primarily on fish, and the other with short thick bills used to forage primarily on plankton. Following this initial divergence in bill shape, size, a known correlate of both prey size and maximum diving depth, diverged rapidly within each of these clades. These results suggest that adaptive radiation in foraging traits underwent initial divergence in bill shape to occupy different food resources, followed by size differentiation to subdivide each niche along the depth axis of the water column., (© 2012 The Author(s). Evolution© 2012 The Society for the Study of Evolution.)

The GOSLON/5-year-old scoring systems have been used in various national and international studies. This study aims to identify the range of Modified Huddart and Bodenham (MHB) scores, which correlate with each GOSLON and 5-year-old category and thereby create a new scoring system that allows comparison with historical data. Two hundred and eighty-three unilateral cleft lip and palate study models from England and Scotland that had all been previously scored using the 5-year-old and GOSLON indices by calibrated examiners were scored using MHB on two separate occasions a month apart by two examiners. Reliability analysis using intraclass correlation and Bland Altman plots were performed. Ordinal regression was used to define the categories of MHB that correspond to the 5-year-old and GOSLON categories. The results revealed a high level of repeatability for both the 5 and the 10 year old models. The MHB scale was grouped into the five categories of the GOSLON and 5-year-old indices with more precision for the 5 year than the 10 year scores. In conclusion, there was high intra-examiner and inter-examiner agreement using the MHB scoring system. It also correlates well with the 5-year-old and GOSLON indices and has been shown to be a much more sensitive scoring system.

Just as features of the physical and biotic environment constrain evolution of ecological and morphological traits, they may also affect evolution of communication systems. Here we analyze constraints on rates of vocal evolution, using a large dataset of New World avian sister taxa. We show that species breeding in tropical forests sing at generally lower frequencies and across narrower bandwidths than species breeding in open habitats, or at high latitudes. We attribute these restrictions on birdsong frequency to the presence of high-frequency insect noise and greater degradation of high-frequency sounds in tropical forests. We fit Ornstein-Uhlenbeck models to show that recent evolution of song frequency has been more greatly constrained in tropical forests than elsewhere, that is, songs have shown less tendency to diverge over time in tropical forests, consistent with inferred acoustic restrictions. In addition, we find that song frequency has evolved more rapidly overall at high latitudes in both forest and open habitats. Besides a larger available sound window, other factors contributing to more rapid divergence at high latitudes may include an overall increased intensity of sexual selection, occupation of more divergent habitats, and the presence of fewer competing species., (© 2012 The Author(s). Evolution© 2012 The Society for the Study of Evolution.)

Within regions, differences in the number of species among clades must be explained by clade age, net diversification rate, or immigration. We examine these alternatives by assessing historical causes of the low diversity of a bird parvorder in the Himalayas (the core Corvoidea, 57 species present), relative to its more species rich sister clade (the Passerida, ∼400 species present), which together comprise the oscine passerines within this region. The core Corvoidea contain ecologically diverse species spanning a large range of body sizes and elevations. Despite this diversity, on the basis of ecological, morphological, and phylogenetic information, we infer that the best explanation for the low number of species within the Himalayan core Corvoidea is one in which ecology limits diversification and/or immigration. Within the core Corvoidea, body size is correlated with elevation: large species are found at high elevations, and small species at lower elevations. This contrasts with the presence of many small-bodied species spanning all elevations in the Passerida and many large bodied species at low elevations in the other orders of birds (the nonpasserines). Cladogenetic events leading to ecological differences between species in body size and shape mostly occurred millions of years ago, and the rate of evolutionary change has declined toward the present. Elevational distributions have been evolutionarily more labile, but are also associated with ancient cladogenetic events. We suggest the core Corvoidea occupy a restricted volume of ecological space in competition with other bird species, and this has limited in situ diversification and/or immigration., (© 2012 The Author(s). Evolution© 2012 The Society for the Study of Evolution.)

The Andes are known to have influenced speciation patterns in many taxa, yet whether species diversification occurred simultaneously with their uplift or only after uplift was complete remains unknown. We examined both the phylogenetic pattern and dates of branching in Adelomyia hummingbirds in relation to Andean uplift to determine whether diversification coincides with the chronological phases of the uplift or with recent climatic fluctuations after Andean formation. Results suggest that the genus Adelomyia originated in the central Andes in the Miocene and was found to be comprised of six deeply divergent phylogroups dating between 3.5 and 6 Ma. The most basal splits in the tree, corresponding to the most southerly distributed of the six phylogroups, diverged in the late Miocene, whereas the northern phylogroups originated during the early-to-mid-Pliocene, when the northern Andes reached heights sufficient to support Adelomyia populations. Although we provide evidence for a southern origin for the group, the subsequent diversification of the northern phylogroups did not strictly follow the hypothesized south-to-north orogeny of the Andes. Further genetic structure within phylogroups may have resulted from Pleistocene climate fluctuations after the onset of the six lineages during the Mio-Pliocene. We explore the processes that promoted diversification in the Andes and suggest that in at least some groups, divergence was coupled to Andean orogeny., (© 2011 Blackwell Publishing Ltd.)

Andean uplift contributed importantly to the build-up of high Neotropical diversity. Final uplift of the Eastern Cordillera of Colombia separated once-contiguous lowland faunas east and west of the Andes between 5 and 3.5 million years ago (Ma hereafter). We used DNA sequences from several moderate- to fast-evolving mitochondrial and two slow-evolving nuclear genes to generate a well-supported phylogeny of Dendrocincla woodcreepers, a genus with multiple species endemic to lowland regions both east and west of the Andes. A time-calibrated phylogeny and dispersal-vicariance analysis indicated that uplift of the Eastern Cordillera of Colombia resulted in the initial vicariant separation of a widespread lowland form east and west of the Andes at c. 3.6 Ma. This was followed by two separate east-to-west dispersal events over or around the completed Andes, each producing a genetically distinct lineage. Our analysis suggests that Andean uplift promoted the build-up of biodiversity in lowland Neotropical faunas both through vicariance-based speciation during uplift and through dispersal-based speciation following uplift. In contrast to the multiple colonizations of the trans-Andean region by Dendrocincla, the Atlantic Forest was colonized from the Amazon only once, followed by in situ diversification., (© 2011 Blackwell Publishing Ltd.)

We ask whether rates of evolution in traits important for reproductive isolation vary across a latitudinal gradient, by quantifying evolutionary rates of two traits important for pre-mating isolation-avian syllable diversity and song length. We analyse over 2500 songs from 116 pairs of closely related New World passerine bird taxa to show that evolutionary rates for the two main groups of passerines-oscines and suboscines-doubled with latitude in both groups for song length. For syllable diversity, oscines (who transmit song culturally) evolved more than 20 times faster at high latitudes than in low latitudes, whereas suboscines (whose songs are innate in most species and who possess very simple song with few syllable types) show no clear latitudinal gradient in rate. Evolutionary rates in oscines and suboscines were similar at tropical latitudes for syllable complexity as well as for song length. These results suggest that evolutionary rates in traits important to reproductive isolation and speciation are influenced by latitude and have been fastest, not in the tropics where species diversity is highest, but towards the poles.

Range expansions are critical to renewed bouts of allopatric or parapatric speciation. Limits on range expansions-and, by implication, speciation-include dispersal ability and permeability of geographical barriers. In addition, recently diverged taxa may interfere with each other, preventing mutual expansion of each other's range into sympatry, because reproductive isolation is incomplete and/or ecological competition particularly strong. On the basis of geographical distributions and mitochondrial DNA phylogenetic information for 418 recently diverged species of New World birds, we estimate that secondary sympatry takes on the order of millions of years following population splitting and hence could impose an important limit on the rate of range expansion, thereby limiting further rounds of species formation. Average rates of achievement of sympatry have been faster in the temperate region (we estimate 1.7 million years to sympatry at 60°) than in the tropics (3.2 million years to sympatry at the equator). Evidence from the ages of species with hybrid zones implies that one factor associated with the slowed sympatry in the tropics is the rate of accumulation of reproductive isolation.

George Gaylord Simpson famously postulated that much of life's diversity originated as adaptive radiations-more or less simultaneous divergences of numerous lines from a single ancestral adaptive type. However, identifying adaptive radiations has proven difficult due to a lack of broad-scale comparative datasets. Here, we use phylogenetic comparative data on body size and shape in a diversity of animal clades to test a key model of adaptive radiation, in which initially rapid morphological evolution is followed by relative stasis. We compared the fit of this model to both single selective peak and random walk models. We found little support for the early-burst model of adaptive radiation, whereas both other models, particularly that of selective peaks, were commonly supported. In addition, we found that the net rate of morphological evolution varied inversely with clade age. The youngest clades appear to evolve most rapidly because long-term change typically does not attain the amount of divergence predicted from rates measured over short time scales. Across our entire analysis, the dominant pattern was one of constraints shaping evolution continually through time rather than rapid evolution followed by stasis. We suggest that the classical model of adaptive radiation, where morphological evolution is initially rapid and slows through time, may be rare in comparative data.

The sudden exchange of mammals over the land bridge between the previously isolated continents of North and South America is among the most celebrated events in the faunal history of the New World. This exchange resulted in the rapid merging of continental mammalian faunas that had evolved in almost complete isolation from each other for tens of millions of years. Yet, the wider importance of land bridge-mediated interchange to faunal mixing in other groups is poorly known because of the incompleteness of the fossil record. In particular, the ability of birds to fly may have rendered a land bridge unnecessary for faunal merging. Using molecular dating of the unique bird faunas of the two continents, we show that rates of interchange increased dramatically after land bridge completion in tropical forest-specializing groups, which rarely colonize oceanic islands and have poor dispersal abilities across water barriers, but not in groups comprised of habitat generalists. These results support the role of the land bridge in the merging of the tropical forest faunas of North and South America. In contrast to mammals, the direction of traffic across the land bridge in birds was primarily south to north. The event transformed the tropical avifauna of the New World.

Most Neotropical lowland forest taxa occur exclusively on one side of the Andes despite the availability of appropriate habitat on both sides. Almost all molecular phylogenies and phylogenetic analyses of species assemblages (i.e. area cladograms) have supported the hypothesis that Andean uplift during the Late Pliocene created a vicariant barrier affecting lowland lineages in the region. However, a few widespread plant and animal species occurring in lowland forests on both sides of the Andes challenge the generality of this hypothesis. To understand the role of the Andes in the history of such organisms, we reconstructed the phylogeographic history of a widespread Neotropical flycatcher (Mionectes oleagineus) in the context of the other four species in the genus. A molecular phylogeny based on nuclear and mitochondrial sequences unambiguously showed an early basal split between montane and lowland Mionectes. The phylogeographic reconstruction of lowland taxa revealed a complex history, with multiple cases in which geographically proximate populations do not represent sister lineages. Specifically, three populations of M. oleagineus west of the Andes do not comprise a monophyletic clade; instead, each represents an independent lineage with origins east of the Andes. Divergence time estimates suggest that at least two cross-Andean dispersal events post-date Andean uplift.

Molecular clocks are widely used to date phylogenetic events, yet evidence supporting the rate constancy of molecular clocks through time and across taxonomic lineages is weak. Here, we present 90 candidate avian clock calibrations obtained from fossils and biogeographical events. Cross-validation techniques were used to identify and discard 16 inconsistent calibration points. Molecular evolution occurred in an approximately clock-like manner through time for the remaining 74 calibrations of the mitochondrial gene, cytochrome b. A molecular rate of approximately 2.1% (+/- 0.1%, 95% confidence interval) was maintained over a 12-million-year interval and across most of 12 taxonomic orders. Minor but significant variance in rates occurred across lineages but was not explained by differences in generation time, body size or latitudinal distribution as previously suggested.

The Common Bush-Tanager (Chlorospingus ophthalmicus) is distributed in Neotropical cloud-forests from Mexico to Argentina and contains 25 subspecies divided into eight subspecies groups based on biogeography, eye coloration, presence of a postocular spot and chest band. All of Central America is occupied by a single subspecies group; whereas the Andes are believed to be occupied by seven additional subspecies groups. We used five mitochondrial genes to investigate the phylogeography and possible species limits of the ophthalmicus complex. A total of 14 monophyletic lineages were uncovered within the ophthalmicus complex, including three clades currently classified as separate species (C. semifuscus, inornatus and tacarcunae). Divergence estimates for these clades date between 0.8 and 5.2 million years ago (Ma). Contrary to expectations based on morphological diversity, phylogeographic structure was greatest in Mexico and Central America and weakest in the Andes. Morphological and genetic divergences were not significantly correlated and most morphologically defined subspecies groups were not supported. Our evidence suggests the ophthalmicus complex originated in Mexico ca. 6.0 Ma (million years ago) and spread south into the Andes ca. 4.7 Ma before the completion of the Isthmus of Panama. Three genetically divergent lineages of ophthalmicus that formed in the Andes possess a complex checkerboard distribution, with a single lineage represented by disjunct populations from Venezuela and the southern Andes, while intervening populations in Ecuador and Central Peru form two genetically and morphologically divergent lineages.