Evolutionary dynamics of hybridization and introgression following the recent colonization of Glossy Ibis (Aves:Plegadis falcinellus) into the New World

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Estimation of species divergence times in presence of cross-species gene flow

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent].

Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.

Rapid introgression of non-native alleles following hybridization between a native Anolis lizard species and a cryptic invader across an urban landscape

Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced-representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non-native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non-native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non-native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non-native genetic material even in the absence of ongoing immigration, indicating that selection favouring non-native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non-native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long-term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.

Distribución geográfica y abundancia poblacional de Plegadis ridgwayi, el ibis de la Puna (Threskiornithidae) con énfasis en las poblaciones del Perú

El ibis de la puna Plegadis ridgwayi, es una especie de Threskiornithidae que habita humedales andinos y realiza migraciones altitudinales hacia la costa. Datos propios, de GBIF, información bibliográfica y del Censo Neotropical de Aves Acuáticas (1992 a 2015) muestran que el ibis de la puna Plegadis ridgwayi se distribuye en Ecuador, Perú, Bolivia, Argentina y Chile, con las mayores densidades poblacionales en Perú y Bolivia en siete y tres localidades respectivamente, que acumulan más del 1% de la población biogeográfica. Se encuentran de 0 a 5000 m de altitud, con las mayores densidades entre 3000 a 4500 m y 0 a 500 m. La mayor incidencia de registros ocurre al sur y centro del Perú, así como costa del centro y norte del Perú. La ampliación de la distribución hacia el norte y costa peruana puede deberse a la disponibilidad ambiental y al deterioro de su hábitat andino. En cuatro humedales costeros del centro del Perú se registraron hasta 818 ibis en 2006, la gran mayoría en Pantanos de Villa y Paraíso. El número de migrantes costeros parece relacionado a la intensidad de sequías en la sierra del Perú central. La abundancia de ibis en el lago altoandino de Junín muestra una disminución histórica, con énfasis después de la sequía de 2004-2005. La expansión distribucional requiere investigar la posible hibridación con las otras especies del género antes alopátridas.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

Genomic Evidence for Rare Hybridization and Large Demographic Changes in the Evolutionary Histories of Four North American Dove Species

Introductions and invasions provide opportunities for interaction and hybridization between colonists and closely related native species. We investigate this phenomenon using the mitochondrial DNA COI and 81,416 base-pairs of overlapping nuclear variation to examine the evolutionary histories and signatures of hybridization among introduced feral Rock Pigeon and Eurasian Collared-Dove and native White-winged and Mourning doves in southwestern North America. First, we report all four species to be highly divergent across loci (overall pair-wise species Φ_ST range = 0.17-0.70) and provide little evidence for gene flow at evolutionary timescales. Despite this, evidence from multiple population genetics analyses supports the presence of six putative contemporary late-stage hybrids among the 182 sampled individuals. These putative hybrids contain various ancestry combinations, but all involve the most populous species, the Mourning Dove. Next, we use a novel method to reconstruct demographic changes through time using partial genome sequence data. We identify recent, species-specific fluctuations in population size that are likely associated with changing environments since the Miocene and suggest that these fluctuations have influenced the genetic diversity of each dove species in ways that may impact their future persistence. Finally, we discuss the importance of using multiple marker types when attempting to infer complex evolutionary histories and propose important considerations when analyzing populations that were recently established or of domestic origins.

The rates of introgression and barriers to genetic exchange between hybridizing species: sex chromosomes vs autosomes

Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Locus-specific introgression in young hybrid swarms: Drift may dominate selection

Closely related species that have previously inhabited geographically separated ranges are hybridizing at an increasing rate due to human disruptions. These human-mediated hybrid zones can be used to study reproductive isolation between species at secondary contact, including examining locus-specific rates of introgression. Introgression is expected to be heterogenous across the genome, reflecting variation in selection. Those loci that introgress especially slowly are good candidates for being involved in reproductive isolation, while those loci that introgress quickly may be involved in adaptive introgression. In the context of conservation, policy makers are especially concerned about introduced alleles moving quickly into the background of a native or endemic species, as these alleles could replace the native alleles in the population, leading to extinction via hybridization. We applied genomic cline analyses to 44,997 SNPs to identify loci introgressing more or less when compared to the genome wide expectation in a human-mediated hybridizing population of red deer and sika in Kintyre Scotland. We found 11.4% of SNPs had cline centres that were significantly different from the genome wide expectation, and 17.6% of all SNPs had excess rates of introgression. Based on simulations, we believe that many of these markers have diverged from the genome-wide average due to drift, rather than because of selection, and we suggest that these simulations can be useful as a null distribution for future studies of genomic clines. Future work on red deer and sika could determine the policy implications of allelic-replacement due to drift rather than selection, and could use replicate, geographically distinct hybrid zones to narrow down those loci that are responding to selection.

The Spatial Signature of Introgression After a Biological Invasion With Hybridization

The accumulation of genome-wide molecular data has emphasized the important role of hybridization in the evolution of many organisms, which may carry introgressed genomic segments resulting from past admixture events with other taxa. Despite a number of examples of hybridization occurring during biological invasions, the resulting spatial patterns of genomic introgression remain poorly understood. Preliminary simulation studies have suggested a heterogeneous spatial level of introgression for invasive taxa after range expansion. We investigated in detail the robustness of this pattern and its persistence over time for both invasive and local organisms. Using spatially explicit simulations, we explored the spatial distribution of introgression across the area of colonization of an invasive taxon hybridizing with a local taxon. The general pattern for neutral loci supported by our results is an increasing introgression of local genes into the invasive taxon with the increase in the distance from the source of the invasion and a decreasing introgression of invasive genes into the local taxon. However, we also show there is some variation in this general trend depending on the scenario investigated. Spatial heterogeneity of introgression within a given taxon is thus an expected neutral pattern in structured populations after a biological invasion with a low to moderate amount of hybridization. We further show that this pattern is consistent with published empirical observations. Using additional simulations, we argue that the spatial pattern of Neanderthal introgression in modern humans, which has been documented to be higher in Asia than in Europe, can be explained by a model of hybridization with Neanderthals in Eurasia during the range expansion of modern humans from Africa. Our results support the view that weak hybridization during range expansion may explain spatially heterogeneous introgression patterns without the need to invoke selection.

New approaches to species delimitation and population structure of anthozoans: Two case studies of octocorals using ultraconserved elements and exons

As coral populations decline worldwide in the face of ongoing environmental change, documenting their distribution, diversity and conservation status is now more imperative than ever. Accurate delimitation and identification of species is a critical first step. This task, however, is not trivial as morphological variation and slowly evolving molecular markers confound species identification. New approaches to species delimitation in corals are needed to overcome these challenges. Here, we test whether target enrichment of ultraconserved elements (UCEs) and exons can be used for delimiting species boundaries and population structure within species of corals by focusing on two octocoral genera, Alcyonium and Sinularia, as exemplary case studies. We designed an updated bait set (29,181 baits) to target-capture 3,023 UCE and exon loci, recovering a mean of 1,910 ± 168 SD per sample with a mean length of 1,055 ± 208 bp. Similar numbers of loci were recovered from Sinularia (1,946 ± 227 SD) and Alcyonium (1,863 ± 177 SD). Species-level phylogenies were highly supported for both genera. Clustering methods based on filtered single nucleotide polymorphisms delimited species and populations that are congruent with previous allozyme, DNA barcoding, reproductive and ecological data for Alcyonium, and offered further evidence of hybridization among species. For Sinularia, results were congruent with those obtained from a previous study using restriction site associated DNA sequencing. Both case studies demonstrate the utility of target-enrichment of UCEs and exons to address a wide range of evolutionary and taxonomic questions across deep to shallow timescales in corals.

Projecting introgression from domestic cats into European wildcats in the Swiss Jura

Hybridization between wild and domesticated organisms is a worldwide conservation issue. In the Jura Mountains, threatened European wildcats (Felis silvestris) have been demographically spreading for approximately the last 50 years, but this recovery is coupled with hybridization with domestic cats (Felis catus). Here, we project the pattern of future introgression using different spatially explicit scenarios to model the interactions between the two species, including competition and different population sizes. We project the fast introgression of domestic cat genes into the wildcat population under all scenarios if hybridization is not severely restricted. If the current hybridization rate and population sizes remain unchanged, we expect the loss of genetic distinctiveness between wild and domestic cats at neutral nuclear, mitochondrial and Y chromosome markers in one hundred years. However, scenarios involving a competitive advantage for wildcats and a future increase in the wildcat population size project a slower increase in introgression. We recommend that future studies assess the fitness of these hybrids and better characterize their ecological niche and their ecological interactions with parental species to elucidate effective conservation measures.

Introgression between non-sister species of honeyeaters (Aves: Meliphagidae) several million years after speciation

We reassessed whether two parapatric non-sister Australian honeyeater species (Aves: Meliphagidae), varied and mangrove honeyeaters (Gavicalis versicolor and G. fasciogularis, respectively), that diverged from a common ancestor c. 2.5 Mya intergrade in the Townsville area of north-eastern Queensland. Consistent with a previous specimen-based study, by using genomics methods we show one-way gene flow for autosomal but not Z-linked markers from varied into mangrove honeyeaters. Introgression barely extends south of the area of parapatry in and around the city of Townsville. While demonstrating the long-term porosity of species boundaries over several million years, our data also suggest a clear role of sex chromosomes in maintaining reproductive isolation.

A Multireference-Based Whole Genome Assembly for the Obligate Ant-Following Antbird, Rhegmatorhina melanosticta (Thamnophilidae)

Current generation high-throughput sequencing technology has facilitated the generation of more genomic-scale data than ever before, thus greatly improving our understanding of avian biology across a range of disciplines. Recent developments in linked-read sequencing (Chromium 10×) and reference-based whole-genome assembly offer an exciting prospect of more accessible chromosome-level genome sequencing in the near future. We sequenced and assembled a genome of the Hairy-crested Antbird (Rhegmatorhina melanosticta), which represents the first publicly available genome for any antbird (Thamnophilidae). Our objectives were to (1) assemble scaffolds to chromosome level based on multiple reference genomes, and report on differences relative to other genomes, (2) assess genome completeness and compare content to other related genomes, and (3) assess the suitability of linked-read sequencing technology for future studies in comparative phylogenomics and population genomics studies. Our R. melanosticta assembly was both highly contiguous (de novo scaffold N50 = 3.3 Mb, reference based N50 = 53.3 Mb) and relatively complete (contained close to 90% of evolutionarily conserved single-copy avian genes and known tetrapod ultraconserved elements). The high contiguity and completeness of this assembly enabled the genome to be successfully mapped to the chromosome level, which uncovered a consistent structural difference between R. melanosticta and other avian genomes. Our results are consistent with the observation that avian genomes are structurally conserved. Additionally, our results demonstrate the utility of linked-read sequencing for non-model genomics. Finally, we demonstrate the value of our R. melanosticta genome for future researchers by mapping reduced representation sequencing data, and by accurately reconstructing the phylogenetic relationships among a sample of thamnophilid species.

Comparison of ultraconserved elements (UCEs) to microsatellite markers for the study of avian hybrid zones: a test in Aphelocoma jays

Objective

Hybrid zones are geographic regions where genetically distinct taxa interbreed, resulting in offspring of mixed ancestry. California Scrub-Jays (Aphelocoma californica) and Woodhouse’s Scrub-Jays (A. woodhouseii) come into secondary contact and hybridize in western Nevada. Although previous work investigated divergence and gene flow between these species using a handful of microsatellite markers, the hybrid zone has not been studied using genome-scale markers, providing an opportunity to assess genome-wide introgression, test for a genetic basis for ecomorphological traits, and compare these estimates to those derived from microsatellites.

Results

Using variant sites flanking ultraconserved elements (UCEs), we performed population assignment and quantified hybrid ancestry for 16 individuals across the zone of secondary contact. Our study included 2468 SNPs distributed throughout the genome, allowing discrimination of genetic affinities of hybrid individuals that were similar to estimates from microsatellites. We show a relationship between bill and wing length and the genetic composition of individuals that was not found in prior work using microsatellites, suggesting a genetic basis for these traits. Our analyses demonstrate the utility of UCEs for the analysis of hybrid zones and provide a basis for future studies to identify the genomic architecture of speciation and phenotypic differences between these incipient species.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4481-z) contains supplementary material, which is available to authorized users.

Climatic shifts drove major contractions in avian latitudinal distributions throughout the Cenozoic

Many higher level avian clades are restricted to Earth's lower latitudes, leading to historical biogeographic reconstructions favoring a Gondwanan origin of crown birds and numerous deep subclades. However, several such "tropical-restricted" clades (TRCs) are represented by stem-lineage fossils well outside the ranges of their closest living relatives, often on northern continents. To assess the drivers of these geographic disjunctions, we combined ecological niche modeling, paleoclimate models, and the early Cenozoic fossil record to examine the influence of climatic change on avian geographic distributions over the last ∼56 million years. By modeling the distribution of suitable habitable area through time, we illustrate that most Paleogene fossil-bearing localities would have been suitable for occupancy by extant TRC representatives when their stem-lineage fossils were deposited. Potentially suitable habitat for these TRCs is inferred to have become progressively restricted toward the tropics throughout the Cenozoic, culminating in relatively narrow circumtropical distributions in the present day. Our results are consistent with coarse-scale niche conservatism at the clade level and support a scenario whereby climate change over geological timescales has largely dictated the geographic distributions of many major avian clades. The distinctive modern bias toward high avian diversity at tropical latitudes for most hierarchical taxonomic levels may therefore represent a relatively recent phenomenon, overprinting a complex biogeographic history of dramatic geographic range shifts driven by Earth's changing climate, variable persistence, and intercontinental dispersal. Earth's current climatic trajectory portends a return to a megathermal state, which may dramatically influence the geographic distributions of many range-restricted extant clades.