Genomic data reveal a protracted window of introgression during the diversification of a neotropical woodcreeper radiation*

Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

Diversification and introgression in four chromosomal taxa of the Pearson's horseshoe bat (Rhinolophus pearsoni) group

Chromosomal variation among closely related taxa is common in both plants and animals, and can reduce rates of introgression as well as promote reproductive isolation and speciation. In mammals, studies relating introgression to chromosomal variation have tended to focus on a few model systems and typically characterized levels of introgression using small numbers of loci. Here we took a genome-wide approach to examine how introgression rates vary among four closely related horseshoe bats (Rhinolophus pearsoni group) that possess different diploid chromosome numbers (2n = 42, 44, 46, and 60) resulting from Robertsonian (Rb) changes (fissions/fusions). Using a sequence capture we obtained orthologous loci for thousands of nuclear loci, as well as mitogenomes, and performed phylogenetic and population genetic analyses. We found that the taxon with 2n = 60 was the first to diverge in this group, and that the relationships among the three other taxa (2n = 42, 44 and 46) showed discordance across our different analyses. Our results revealed signatures of multiple ancient introgression events between the four taxa, with evidence of mitonuclar discordance in phylogenetic trees and reticulation events in their evolutionary history. Despite this, we found no evidence of recent and/or ongoing introgression between taxa. Overall, our results indicate that the effects of Rb changes on the reduction of introgression are complicated and that these may contribute to reproductive isolation and speciation in concert with other factors (e.g. phenotypic and genic divergence).

Phylogeny and historical biogeography of the water boatmen (Insecta: Hemiptera: Heteroptera: Nepomorpha: Corixoidea)

The water boatmen of Corixoidea, a group of aquatic bugs with more than 600 extant species, is one of the largest superfamilies of Nepomorpha. Contrary to the other nepomorphan lineages, the Corixoidea are most diverse in the Laurasian remnant Holarctic region. To explicitly test whether the present-day Holarctic distribution of diverse corixids is associated with the arising of the Laurasian landmass that was separated from Gondwana, we investigated the phylogeny, divergence times and historical biogeography of Corixoidea based on morphological and molecular characters sampled from 122 taxa representing all families, subfamilies, tribes and approximately 54 % of the genera. Our results were largely congruent with the phylogenetic relationships within the established nepomorphan phylogenetic context. The fossil calibrated chronogram, diversification analysis and ancestral ranges reconstruction indicated that Corixoidea began to diversify in Gondwana in the late Triassic approximately at 224 Ma and the arising of the most diverse subfamily Corixinae in Corixidae in the Holarctic region was largely congruent with the time of separation of Laurasia from Gondwana. The large-scale expansion of the temperate and cold zones on the northward-moving Laurasian landmass after the breakup of the Pangea provided new aquatic niches and ecological opportunities for promoting rapid diversification for the Holarctic corixid lineage.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Ecology and the origin of non-ephemeral species

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.

Gene flow, genomic homogenization and the timeline to speciation in Amazonian manakins

Phylogeographic studies of the most species rich region of the planet - the Amazon basin - repeatedly uncover 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 datasets together with detailed phylogeographic 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 geographic 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 1.4 Ma basal node 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 phylogeographic studies and would help determine species classification of the plethora of newly identified lineages in the Amazon basin and other regions.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Genome-wide sequence data show no evidence of hybridization and introgression among pollinator wasps associated with a community of Panamanian strangler figs

The specificity of pollinator host choice influences opportunities for reproductive isolation in their host plants. Similarly, host plants can influence opportunities for reproductive isolation in their pollinators. For example, in the fig and fig wasp mutualism, offspring of fig pollinator wasps mate inside the inflorescence that the mothers pollinate. Although often host specific, multiple fig pollinator species are sometimes associated with the same fig species, potentially enabling hybridization between wasp species. Here, we study the 19 pollinator species (Pegoscapus spp.) associated with an entire community of 16 Panamanian strangler fig species (Ficus subgenus Urostigma, section Americanae) to determine whether the previously documented history of pollinator host switching and current host sharing predicts genetic admixture among the pollinator species, as has been observed in their host figs. Specifically, we use genome‐wide ultraconserved element (UCE) loci to estimate phylogenetic relationships and test for hybridization and introgression among the pollinator species. In all cases, we recover well‐delimited pollinator species that contain high interspecific divergence. Even among pairs of pollinator species that currently reproduce within syconia of shared host fig species, we found no evidence of hybridization or introgression. This is in contrast to their host figs, where hybridization and introgression have been detected within this community, and more generally, within figs worldwide. Consistent with general patterns recovered among other obligate pollination mutualisms (e.g. yucca moths and yuccas), our results suggest that while hybridization and introgression are processes operating within the host plants, these processes are relatively unimportant within their associated insect pollinators.

OUP accepted manuscript

Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve “ghosts,” that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA–BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.]

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.

Prevalence and Adaptive Impact of Introgression

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus)

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Most Genomic Loci Misrepresent the Phylogeny of an Avian Radiation Because of Ancient Gene Flow

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis, and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the "winner-takes-all democratic majority tree" is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger data sets, especially whole-genome data sets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.[Bifurcation; introgression; recombination; reticulation; Phylloscopus.].

Genome-wide SNPs clarify a complex radiation and support recognition of an additional cat species

Phylogenetic reconstruction and species delimitation are often challenging in the case of recent evolutionary radiations, especially when postspeciation gene flow is present. Leopardus is a Neotropical cat genus that has a long history of recalcitrant taxonomic problems, along with both ancient and current episodes of interspecies admixture. Here, we employ genome-wide SNP data from all presently recognized Leopardus species, including several individuals from the tigrina complex (representing Leopardus guttulus and two distinct populations of Leopardus tigrinus), to investigate the evolutionary history of this genus. Our results reveal that the tigrina complex is paraphyletic, containing at least three distinct species. While one can be assigned to L. guttulus, the other two remain uncertain regarding their taxonomic assignment. Our findings indicate that the “tigrina” morphology may be plesiomorphic within this group, which has led to a longstanding taxonomic trend of lumping these poorly known felids into a single species.

Holocene population expansion of a tropical bee coincides with early human colonization of Fiji rather than climate change

There is substantial debate about the relative roles of climate change and human activities on biodiversity and species demographies over the Holocene. In some cases, these two factors can be resolved using fossil data, but for many taxa such data are not available. Inferring historical demographies of taxa has become common, but the methodologies are mostly recent and their shortcomings often unexplored. The bee genus Homalictus is developing into a tractable model system for understanding how native bee populations in tropical islands have responded to past climate change. We greatly expand on previous studies using sequences of the mitochondrial gene COI from 474 specimens and between 171 and 3928 autosomal (DArTSeq) single nucleotide polymorphism loci from 19 specimens of the native Fijian bee, Homalictus fijiensis, to explore its historical demography using coalescent and mismatch analyses. We ask whether past changes in demography were human- or climate-driven, while considering analytical assumptions. We show that inferred changes in population sizes are too recent to be explained by past climate change. Instead we find that a dramatic increase in population size for the main island of Viti Levu coincides with increasing occupation by humans and their modification of the environment. We found no corresponding change in bee population size for another major island, Kadavu, where human populations and agricultural activities have been historically very low. Our analyses indicate that molecular approaches can be used to disentangle the impacts of humans and climate change on a major tropical pollinator and that stringent analytical approaches are required for reliable interpretation of results.

Molecular systematics and phylogeography of a widespread Neotropical avian lineage: evidence for cryptic speciation with protracted gene flow throughout the Late Quaternary

Here we use an integrative approach, including coalescent-based methods, isolation–migration and species distribution models, to infer population structure, divergence times and diversification in the two species of the genus Cymbilaimus (Aves, Thamnophilidae). Our results support a recent and rapid diversification with both incomplete lineage sorting and gene flow shaping the evolutionary history of Cymbilaimus. The spatio-temporal pattern of cladogenesis suggests that Cymbilaimus originated in the north/western portion of cis-Andean South America and then diversified into the Brazilian Shield and Central America after consolidation of the modern Amazonian drainage and the Andean range. This evolutionary scenario is explained by cycles of range expansion and dispersal, followed by isolation, and recurrent gene flow, during the last 1.2 Myr. Our results agree with those recently reported for other closely related suboscine lineages, whereby the window of introgression between closely related taxa remains open for up to a few million years after their original split. In Cymbilaimus, introgression was recurrent between C. lineatus and C. sanctaemariae, even after they acquired vocal and ecological differentiation, supporting the claim that at least in Neotropical suboscines, full reproductive compatibility may take millions of years to evolve and cannot be interpreted as synonymous with a lack of speciation.

Avian Diversity: Speciation, Macroevolution, and Ecological Function

The origin, distribution, and function of biological diversity are fundamental themes of ecology and evolutionary biology. Research on birds has played a major role in the history and development of these ideas, yet progress was for many decades limited by a focus on patterns of current diversity, often restricted to particular clades or regions. Deeper insight is now emerging from a recent wave of integrative studies combining comprehensive phylogenetic, environmental, and functional trait data at unprecedented scales. We review these empirical advances and describe how they are reshaping our understanding of global patterns of bird diversity and the processes by which it arises, with implications for avian biogeography and functional ecology. Further expansion and integration of data sets may help to resolve longstanding debates about the evolutionary origins of biodiversity and offer a framework for understanding and predicting the response of ecosystems to environmental change.

Digest: Avian genomes are permeable to introgression for a few million years

How wide is the time window for introgression after divergence? Pulido-Santacruz et al. addressed this question by studying the evolutionary history of the bird genus Dendrocincla. They found five introgression events that occurred between a few hundred thousand and around 2.5 million years after divergence. The introgressed genomic proportion declined exponentially with the age of the hybridizing taxa. Hence, this study suggests that species boundaries are permeable up to 2.5 million years after divergence.