Consequences of divergence and introgression for speciation in Andean cloud forest birds

An Elevational Phylogeographic Diversity Gradient in Neotropical Birds Is Decoupled from Speciation Rates

AbstractA key question about macroevolutionary speciation rates is whether they are controlled by microevolutionary processes operating at the population level. For example, does spatial variation in population genetic differentiation underlie geographical gradients in speciation rates? Previous work suggests that speciation rates increase with elevation in Neotropical birds, but underlying population-level gradients remain unexplored. Here, we characterize elevational phylogeographic diversity between montane and lowland birds in the megadiverse Andes-Amazonian system and assess its relationship to speciation rates to evaluate the link between population-level differentiation and species-level diversification. We aggregated and georeferenced nearly 7,000 mitochondrial DNA sequences across 103 species or species complexes in the Andes and Amazonia and used these sequences to describe phylogeographic differentiation across both regions. Our results show increased levels of both discrete and continuous metrics of population structure in the Andean mountains compared with the Amazonian lowlands. However, higher levels of population differentiation do not predict higher rates of speciation in our dataset. Multiple potential factors may lead to our observed decoupling of initial population divergence and speciation rates, including the ephemerality of incipient species and the multifaceted nature of the speciation process, as well as methodological challenges associated with estimating rates of population differentiation and speciation.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

The dynamics of introgression across an avian radiation

Hybridization and resulting introgression can play both a destructive and a creative role in the evolution of diversity. Thus, characterizing when and where introgression is most likely to occur can help us understand the causes of diversification dynamics. Here, we examine the prevalence of and variation in introgression using phylogenomic data from a large (1300+ species), geographically widespread avian group, the suboscine birds. We first examine patterns of gene tree discordance across the geographic distribution of the entire clade. We then evaluate the signal of introgression in a subset of 206 species triads using Patterson's D‐statistic and test for associations between introgression signal and evolutionary, geographic, and environmental variables. We find that gene tree discordance varies across lineages and geographic regions. The signal of introgression is highest in cases where species occur in close geographic proximity and in regions with more dynamic climates since the Pleistocene. Our results highlight the potential of phylogenomic datasets for examining broad patterns of hybridization and suggest that the degree of introgression between diverging lineages might be predictable based on the setting in which they occur.

A new species of wood lizard (Hoplocercinae, Enyalioides) from the Río Huallaga Basin in Central Peru

We report the discovery of a new species of Enyalioides from the premontane forest of the Río Huallaga basin in central Peru. The most similar and phylogenetically related species are E. binzayedi and E. rudolfarndti. However, the new species differs from E. binzayedi (state of character in parentheses) by having dorsal scales strongly keeled on paravertebral region and feebly keeled or smooth elsewhere (prominent medial keel on each dorsal scale), more dorsals in transverse row between dorsolateral crests at midbody 26–39, x̄ = 30.44 (22–31, x̄ = 27.57), and a conspicuous posteromedial black patch in the gular region of males (absent). Contrarily, adult males of the new species and E. rudolfarndti are readily distinguished by having a conspicuous posteromedial black patch in the gular region, absent in E. rudolfarndti, and by lacking a conspicuous orange blotch (faint if present) on the antehumeral region, as in E. rudolfarndti. We also present an updated molecular phylogenetic tree of hoplocercines, which strongly supports both referral of the newly discovered species to Enyalioides and its status as distinct from those recognized previously.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

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.

Rugged relief and climate promote isolation and divergence between two neotropical cold-associated birds

The role of historical factors in establishing patterns of diversity in tropical mountains is of interest to understand the buildup of megadiverse biotas. In these regions, the historical processes of range fragmentation and contraction followed by dispersal are thought to be mediated by the interplay between rugged relief (complex topography) and climate fluctuations and likely explain most of the dynamics of diversification in plants and animals. Although empirical studies addressing the interaction between climate and topography have provided invaluable insights into population divergence and speciation patterns in tropical montane organisms, a more detailed and robust test of such processes in an explicit spatio-temporal framework is still lacking. Consequently, our ability to gain insights into historical range shifts over time and the genomic footprint left by them is limited. Here, we used niche modeling and subgenomic population-level datasets to explore the evolution of two species of warbling finches (genus Microspingus) disjunctly distributed across the Montane Atlantic Forest, a Neotropical region with complex geological and environmental histories. Population structure inferences suggest a scenario of three genetically differentiated populations, which are congruent with both geography and phenotypic variation. Demographic simulations support asynchronous isolation of these populations as recently as ∼40,000 years ago, relatively stable population sizes over recent time, and past gene flow subsequent to divergence. Throughout the last 800,000 years, niche models predicted extensive expansion into lowland areas with increasing overlap of species distributions during glacial periods, with prominent retractions and isolation into higher altitudes during interglacials, which are in line with signs of introgression of currently isolated populations. These results support a dual role of cyclical climatic changes: population divergence and persistence in mountain tops during warm periods followed by periods of expansion and admixture in lower elevations during cold periods. Our results underscore the role of the interplay between landscape and climate as an important mechanism in the evolution of the Neotropical montane biota.

DNA barcodes of birds from northern Colombia

DNA barcode datasets are a useful tool for conservation and aid in taxonomic identification, particularly in megadiverse tropical countries seeking to document and describe its biota, which is dropping at an alarming rate during recent decades. Here we report the barcodes for several low elevation bird species from northern Colombia with the goal to provide tools for species identification in this region of South America. We blood-sampled birds in a lowland tropical forest with various degrees of intervention using standard 3 × 12 m mist-nets. We extracted DNA and sequenced the COI barcode gene using standard primers and laboratory methods. We obtained 26 COI sequences from 18 species, 10 families and three orders and found that barcodes largely matched (but not always) phenotypic identification (> 90%) and they also facilitated the identification of several challenging passerine species. Despite our reduced sampling, our study represents the first attempt to document COI barcodes for birds (from blood samples) in this part of Colombia, which fills a considerable gap of sampling in this part of South America.

Pervasive Genomic Signatures of Local Adaptation to Altitude Across Highland Specialist Andean Hummingbird Populations

Populations along steep environmental gradients are subject to differentiating selection that can result in local adaptation, despite countervailing gene flow, and genetic drift. In montane systems, where species are often restricted to narrow ranges of elevation, it is unclear whether the selection is strong enough to influence functional differentiation of subpopulations differing by a few hundred meters in elevation. We used targeted capture of 12 501 exons from across the genome, including 271 genes previously implicated in altitude adaptation, to test for adaptation to local elevations for 2 highland hummingbird species, Coeligena violifer (n = 62) and Colibri coruscans (n = 101). For each species, we described population genetic structure across the complex geography of the Peruvian Andes and, while accounting for this structure, we tested whether elevational allele frequency clines in single nucleotide polymorphisms (SNPs) showed evidence for local adaptation to elevation. Although the 2 species exhibited contrasting population genetic structures, we found signatures of clinal genetic variation with shifts in elevation in both. The genes with SNP-elevation associations included candidate genes previously discovered for high-elevation adaptation as well as others not previously identified, with cellular functions related to hypoxia response, energy metabolism, and immune function, among others. Despite the homogenizing effects of gene flow and genetic drift, natural selection on parts of the genome evidently optimizes elevation-specific cellular function even within elevation range-restricted montane populations. Consequently, our results suggest local adaptation occurring in narrow elevation bands in tropical mountains, such as the Andes, may effectively make them "taller" biogeographic barriers.

A new species of frog (Terrarana, Strabomantidae, Phrynopus) from the Peruvian Andean grasslands

We describe a new, medium-sized species of terrestrial frog of the genus Phrynopus from a single locality in the central Andes of Peru (Departamento de Huánuco) at 3,730 meters of elevation. Phylogenetic analyses supported Phrynopus remotum sp. nov. as an independent lineage, sister to most of its congeners. The new species is morphologically distinguishable by the presence of small tubercles on upper eyelids and heels, an areolate venter, and the absence of dorsolateral folds or ridges. This species inhabits the highlands adjacent to the Marañón Dry valley. The only sympatric amphibian species recorded is the marsupial frog Gastrotheca peruana.

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.

Assessing biological factors affecting postspeciation introgression

An increasing number of phylogenomic studies have documented a clear “footprint” of postspeciation introgression among closely related species. Nonetheless, systematic genome‐wide studies of factors that determine the likelihood of introgression remain rare. Here, we propose an a priori hypothesis‐testing framework that uses introgression statistics—including a new metric of estimated introgression, D_p—to evaluate general patterns of introgression prevalence and direction across multiple closely related species. We demonstrate this approach using whole genome sequences from 32 lineages in 11 wild tomato species to assess the effect of three factors on introgression—genetic relatedness, geographical proximity, and mating system differences—based on multiple trios within the “ABBA–BABA” test. Our analyses suggest each factor affects the prevalence of introgression, although our power to detect these is limited by the number of comparisons currently available. We find that of 14 species pairs with geographically “proximate” versus “distant” population comparisons, 13 showed evidence of introgression; in 10 of these cases, this was more prevalent between geographically closer populations. We also find modest evidence that introgression declines with increasing genetic divergence between lineages, is more prevalent between lineages that share the same mating system, and—when it does occur between mating systems—tends to involve gene flow from more inbreeding to more outbreeding lineages. Although our analysis indicates that recent postspeciation introgression is frequent in this group—detected in 15 of 17 tested trios—estimated levels of genetic exchange are modest (0.2–2.5% of the genome), so the relative importance of hybridization in shaping the evolutionary trajectories of these species could be limited. Regardless, similar clade‐wide analyses of genomic introgression would be valuable for disentangling the major ecological, reproductive, and historical determinants of postspeciation gene flow, and for assessing the relative contribution of introgression as a source of genetic variation.

On geographic barriers and Pleistocene glaciations: Tracing the diversification of the Russet-crowned Warbler (Myiothlypis coronata) along the Andes

We studied the phylogeography and plumage variation of the Russet-crowned Warbler (Myiothlypis coronata), from Venezuela to Bolivia, with focus on populations from Ecuador and northern Peru. We analyzed sequences of mitochondrial and nuclear genes, geographic distributions, as well as photographs of specimens deposited at museum collections. Phylogenetic analyses identified three major lineages formed by populations from: Venezuela and Colombia (M. c. regulus), Ecuador and northern Peru (M. elata, M. castaneiceps, M. orientalis, M. c. chapmani), and central Peru and Bolivia (M. c. coronata). We found further population structure within M. c. regulus and M. c. coronata, and population structure and complexity of plumage variation within the Ecuador-northern Peru lineage. Time-calibrated trees estimated that most intraspecific variation originated during the Pleistocene; however, this pattern may not be attributed to an increase in diversification rate during that period. We discuss these results in the context of the importance of geographic-ecological barriers in promoting lineage diversification along the Andes and put forward a preliminary taxonomic proposal for major lineages identified in this study.

Fine scale mapping of genomic introgressions within the Drosophila yakuba clade

The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. However, the amount of interspecific genetic exchange (introgression) mediated by hybridization remains largely unknown, although recent progress in genome sequencing has made measuring introgression more tractable. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of São Tomé, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. In this report, we quantify the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). We sequenced the genomes of 86 individuals from all three species. We also developed and applied a new statistical framework, using a hidden Markov approach, to identify introgression. We found that introgression has occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, we find that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). We additionally show that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, our results indicate that genetic exchange between these species is currently rare.