- Open Access
Avian introgression in the genomic era
© The Author(s) 2017
Received: 20 July 2017
Accepted: 16 November 2017
Published: 30 November 2017
Introgression, the incorporation of genetic material from one (sub)species into the gene pool of another by means of hybridization and backcrossing, is a common phenomenon in birds and can provide important insights into the speciation process. In the last decade, the toolkit for studying introgression has expanded together with the development of molecular markers. In this review, we explore how genomic data, the most recent step in this methodological progress, impacts different aspects in the study of avian introgression. First, the detection of hybrids and backcrosses has improved dramatically. The most widely used software package is STRUCTURE. Phylogenetic discordance (i.e. different loci resulting in discordant gene trees) is another means for the detection of introgression, although it should be regarded as a starting point for further analyses, not as a definitive proof of introgression. Specifically, disentangling introgression from other biological processes, such as incomplete lineage sorting, remains a challenging endeavour, although new techniques, such as the D-statistic, are being developed. In addition, phylogenetics might require a shift from trees to networks. Second, the study of hybrid zones by means of geographical or genomic cline analysis has led to important insights into the complex interplay between hybridization and speciation. However, because each hybrid zone study is just a single snapshot of a complex and continuously changing interaction, hybrid zones should be studied across different temporal and/or spatial scales. A third powerful tool is the genome scan. The debate on which evolutionary processes underlie the genomic landscape is still ongoing, as is the question whether loci involved in reproductive isolation cluster together in ‘islands of speciation’ or whether they are scattered throughout the genome. Exploring genomic landscapes across the avian tree of life will be an exciting field for further research. Finally, the findings from these different methods should be incorporated into specific speciation scenarios, which can consequently be tested using a modelling approach. All in all, this genomic perspective on avian hybridization and speciation will further our understanding in evolution in general.
In this review, we focus on the final step in this methodological progress: genomic data, which has become standard practise in ornithology (Kraus and Wink 2015; Jarvis 2016; Oyler-McCance et al. 2016; Toews et al. 2016). Following Toews et al. (2016), we consider the following next generation sequencing techniques as genomic tools: genome sequencing and resequencing, reduced representation techniques (genotype-by-sequencing [GBS] and restriction-site-associated DNA sequencing [RADseq]), sequence capture and RNA sequencing. In this review, we explore how genomic data impacts the different aspects in the study of avian introgression, such as detecting introgression, hybrid zone studies, and the role of introgression in genome architecture.
One of the first steps in avian hybridization research is the identification of admixed individuals. To tackle this issue several software packages have been developed, such as NewHybrids (Anderson and Thompson 2002), AFLPOP (Duchesne and Bernatchez 2002), BAPS (Corander et al. 2004) and HYBRIDLAB (Nielsen et al. 2006). But the most widely used software package is STRUCTURE (Pritchard et al. 2000). Indeed, 62 out of 87 studies (about 71%) on avian introgression applied this software package (Additional file 1: Table S1). STRUCTURE uses multilocus genotype data and is based on a clustering algorithm that assigns individuals to populations. Early studies used microsatellites, often in combination with mitochondrial markers (e.g., Barilani et al. 2005). The rapid progress in sequencing techniques introduced the application of SNPs and other genome-wide markers to this method (Saetre et al. 2003; Kraus et al. 2012). A new and updated version, fastSTRUCTURE, allows for the analysis of large SNP datasets (Raj et al. 2014; Elgvin et al. 2017).
STRUCTURE, however popular it remains since over a decade, does have limitations with respect to the underlying population genetic model, such as adherence to Hardy–Weinberg and linkage equilibria (Jombart et al. 2010). Alternatives to STRUCTURE have been developed to address these issues, such as ADMIXTURE (Alexander et al. 2009) and Discriminant Analysis of Principle Components or DAPC (Jombart et al. 2010). Eventually, STRUCTURE, ADMIXTURE and DAPC (and other similar software) are best utilized alongside each other (Frosch et al. 2014). However, all genomic studies in our literature search relied solely on STRUCTURE.
Hybridization and consequent exchange of genetic material can become apparent in phylogenetic analyses with different loci resulting in discordant gene trees (Degnan and Rosenberg 2009). Hence, this gene tree discordance, also referred to as phylogenetic incongruence, can be used to detect introgression. For example, in a phylogenetic analysis of woodpeckers, Fuchs et al. (2013) attributed conflicting topologies between several loci to an ancient hybridization event between members of the Campephilus and the melanerpine (genera Melanerpes and Sphyrapicus) lineages.
However, introgression is not the only process culminating in phylogenetic incongruence. Other biological processes, such as incomplete lineage sorting and gene duplication, can amount to similar patterns (Maddison 1997; Degnan and Rosenberg 2009). Particularly, incomplete lineage sorting, when lineages fail to coalesce in the ancestral population of two species, seems to be pervasive (Pamilo and Nei 1988). Hence, phylogenetic discordance should not be seen as definitive proof of introgression, but rather as a starting point for further analyses. It is advised to provide other lines of evidence to show that gene tree discordance is the outcome of introgressive hybridization, not incomplete lineage sorting.
However, an excess of ABBA or BABA can also arise from other processes, such as non-random mating in the ancestral population due to population structure (Eriksson and Manica 2012). Also, the D-statistic was originally developed to infer introgression on a genome-wide or chromosome-wide scale (Green et al. 2010). Calculating this statistic for small genomic regions or specific loci in order to characterize patterns of introgression across the genome can lead to unreliable results, because significant D-statistics tend to cluster in regions of reduced genetic diversity (Martin et al. 2014). It is, therefore, advised to apply the D-statistic for the detection of genome-wide or chromosome-wide introgression, not to single out possible introgressed regions.
We found five studies on avian introgression that adopted the D-statistic to study introgression patterns in Zimmerius flycatchers (Rheindt et al. 2014), Corvus crows (Poelstra et al. 2014), Darwin’s Finches (Lamichhaney et al. 2015), Aphelocoma jays (Zarza et al. 2016), and Passer sparrows (Elgvin et al. 2017). Most studies applied the D-statistic in the appropriate way (i.e. to infer genome-wide introgression), only the study on Zimmerius flycatchers seems to fall into the trap of using the ABBA-BABA-test to pinpoint introgressed loci. In this study, Rheindt et al. (2014) attempt to reconstruct the evolutionary history of a phenotypically mosaic population of the Peruvian Tyrannulet (Z. viridiflavus), which shares plumage characteristics with the Golden-faced Tyrannulet (Z. chrysops). Mapping the results from the ABBA-BABA-test to the Zebra Finch (Taeniopygia guttata) genome, they find that some ABBA-favoured SNPs are close to genes involved in cell projection and plasma membranes. Given the connection between cell membranes and plumage coloration, the authors speculate that there might have been introgression of plumage alleles between these populations. However, as discussed above, these loci have not necessarily introgressed, they might represent genomic regions of low genetic diversity due to purifying selection. Examining the level of divergence in putative introgressed regions is recommended here.
Another way to deal with introgression and incomplete lineage sorting is to switch from phylogenetic trees to phylogenetic networks (Huson and Bryant 2006; Ottenburghs et al. 2016, 2017). Some methods for constructing phylogenetic networks allow for hybridization while ignoring incomplete lineage sorting (Beiko and Hamilton 2006), whereas other methods take into account incomplete lineage sorting while ignoring hybridization (Maddison and Knowles 2006; Than and Nakhleh 2009). New methods are being developed to reconstruct phylogenetic networks taking both incomplete lineage sorting and hybridization into account (Joly et al. 2009; Kubatko 2009; Meng and Kubatko 2009; Yu et al. 2013; Wen et al. 2016). This network approach has not been applied to avian genomic data yet, but promises to be a fruitful strategy for the detection of introgression.
Genomic clines in hybrid zones
Apart from the descriptive nature of detecting introgression, patterns of interspecific genetic exchange can also be used to study the genetics of reproductive isolation and speciation (Harrison and Larson 2014; Payseur and Rieseberg 2016). The extent of introgression of a specific allele depends on several factors, such as hybrid fitness, reproductive isolation and genetic linkage (Barton 1979; Wu 2001; Payseur 2010). Alleles can be roughly divided into three categories: (1) neutral alleles that are free to flow between species, (2) alleles that confer an adaptive advantage and introgress quickly, and (3) alleles that lead to reduced fitness and inhibit gene flow. Hybrid genomes are a mosaic of these three categories, mingled by migration and recombination (Payseur 2010; Wang et al. 2011). Hence, most species boundaries are semipermeable: some genomic regions (e.g., those leading to reduced hybrid fitness) show restricted gene flow while other regions (e.g., comprising neutral or advantageous alleles) are allowed to flow freely.
Hybrid zones, regions where two genetically distinct populations interbreed, are excellent natural laboratories to explore these locus-specific patterns of introgression (Hewitt 1988; Harrison 1990; Harrison and Larson 2016), especially in combination with the application of geographical cline theory (Barton and Hewitt 1985). Cline theory provides a framework to analyse changes in traits or allele frequencies as a function of geographic distance across a hybrid zone transect. Several characteristics of the observed clines can be used to make inferences about hybrid zone dynamics. For instance, cline width in combination with dispersal rates allows for estimation of selection pressures (Barton and Gale 1993). Alleles and traits under similar selective pressures will show concordant cline widths and centres, whereas those subject to different selection pressures will show displaced cline centres compared to the majority of the other clines (Barton 1983).
Genomic studies on avian hybrid zones
Parchman et al. (2013)
Alcaide et al. (2014)
Baldassarre et al. (2014)
Taylor et al. (2014)
Seneviratne et al. (2016)
Geographic and Genomic
Grossen et al. (2016)
Kingston et al. (2017)
Characterizing an avian hybrid zone by means of geographical or genomic cline analyses can provide important insights into the genetic underpinnings of reproductive isolation between the hybridizing species. But each hybrid zone study is just a single snapshot of a complex and continuously changing interaction. To capture the dynamic nature of hybrid zones, one can study a particular hybrid zone across different temporal and/or spatial scales. For example, by comparing historical (2000‒2002) and recent (2010‒2012) genetic data, Taylor et al. (2014) showed that the hybrid zone between Black-capped Chickadee (Poecile atricapillus) and Carolina Chickadee (P. carolensis) in Pennsylvania has moved north due to changing winter temperatures. However, a temporal comparison is only possible when historical data is available or can be obtained from museum specimens (Spurgin et al. 2014; Linck et al. 2017). An alternative strategy is to compare different hybrid zones between the same species (Schaefer et al. 2016; Kingston et al. 2017; Lackey and Boughman 2017). On the one hand, concordant patterns of introgression and differentiation across multiple hybrid zones might pinpoint genes that are important in reproductive isolation or adaptation, regardless of any environmental differences between the hybrid zones. On the other hand, discordant patterns of introgression and adaptation might be used to identify loci under environment-specific selection or drift within the respective hybrid zones. Finally, the spatial comparison of different hybrid zones can be extended to multiple pairs of closely related species. By analysing contact zones between these hybridizing species, patterns of introgression and differentiation can be related to the age of the interacting species, providing insights into the build-up of reproductive isolation and differentiation over time (Grossen et al. 2016; Vijay et al. 2016).
Exploring the genomic landscape
Genomic data allows researchers to zoom out from locus-specific introgression and differentiation uncovered in hybrid zones and study patterns across the whole genome. Here, genome scans provide a powerful approach (Haasl and Payseur 2016): align the genomes of two species, slide a window across them and calculate a divergence statistic—mostly F ST. The resulting picture is a genomic landscape with islands of differentiated regions in a sea of neutral variation. These islands are commonly referred to as ‘differentiation islands’ (Harr 2006) or ‘speciation islands’ (Turner et al. 2005). As the latter term suggests, these islands have often been related to speciation, in particular to the genic view of speciation (Wu 2001) in combination with divergence-with-gene-flow (Pinho and Hey 2010). This view holds that divergent selection against gene flow is initially restricted to a few loci. These loci contribute to reproductive isolation and are less likely to introgress compared to selectively neutral loci. Hence, these loci and closely linked genomic regions are expected to diverge while gene flow homogenizes the remainder of the genome (Feder et al. 2012; Via 2012). Furthermore, if genomic islands of divergence are the outcome of reduced gene flow, one expects to find genes contributing to reproductive isolation within them. An alternative explanation for the formation of genomic islands is that they arose through positive and purifying selection (including background selection at linked sites, referred to as ‘linked selection’) in allopatry, independent of reproductive isolation and gene flow (Cruickshank and Hahn 2014).
Avian studies on genomic islands
Ellegren et al. (2012)
Parchman et al. (2013)
Migration, plumage, and song (−)
Ruegg et al. (2014)
Delmore et al. (2015)
Lavretsky et al. (2015)
4 species of Ficedula
Burri et al. (2015)
Plumage coloration (+)
Poelstra et al. (2014)
Plumage coloration (+)
Vijay et al. (2016)
16 species of Darwin’s finches
Bill morphology (+)
Lamichhaney et al. (2015)
Grossen et al. (2016)
Elgvin et al. (2017)
If ‘islands of differentiation’ play a pivotal role in the origin of new species, one expects to find genes involved in reproductive isolation within these islands. This expectation has been confirmed for Carrion Crows (Corvus corone corone) and Hooded Crows (C. c. cornix). A genome scan uncovered a highly differentiated genomic region which contains genes involved in pigmentation and visual perception. This result was further corroborated by differential patterns of gene expression (Poelstra et al. 2014). Genomic analyses of contact zones between other Corvus subspecies (cornix, corone, orientalis and pectoralis) revealed clustering of certain pigmentation genes, albeit in different genomic islands (Vijay et al. 2016). Other studies, however, have shown that genes potentially involved in reproductive isolation do not always cluster together, instead they are scattered across the genome (Parchman et al. 2013; Ruegg et al. 2014; but see Delmore et al. 2015). These contrasting findings—candidate ‘speciation genes’ clustered in genomic islands versus scattered throughout the genome—suggests that the genetic basis of speciation is highly species-specific and context-dependent.
The approaches discussed above (e.g., phylogenetic discordance, genome scans and genomic cline analysis) are best utilized alongside each other to achieve a complete picture of the speciation and hybridization process. For instance, Parchman et al. (2013) combined genomic cline analyses and genome scans (based on F ST) to investigate the Panamanian hybrid zone between White-collared Manakin (Manacus candei) and Golden-collared Manakin (M. vitellinus). They showed that differentiated loci with high F ST and loci with discordant clines did not cluster together in the genome. This is in line with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. Most studies on avian introgression, however, relied on one approach (only 58 of the 165 studies [35%] applied multiple methods).
To discover and describe hybrid zones in the past, many studies relied on ‘randomly’ sampled specimens, mostly collected within and outside the putative hybrid zone. Sampling strategies that take into account the resolution of genomic data will need to be developed. Moreover, each hybrid zone study is just a single snapshot of a complex and continuously changing interaction. To capture the dynamic nature of hybrid zones in a genomic context, sampling should be carried out across different temporal and spatial scales.
Due to the increasing amount of genetic data, biology is moving from a hypothesis-driven to a data-driven science. The same pattern can be observed in ornithological research. This shift brings many computational challenges. Most software packages that are routinely used in avian introgression studies (e.g., STRUCTURE, IMa) cannot cope with large genomic datasets (Darriba et al. 2015). Therefore, new computational tools will need to be developed to handle the increasing amount of genomic data. For instance, ExaML, a computationally more efficient version of the maximum likelihood program RAxML, was developed to analyse 48 complete avian genomes (Kozlov et al. 2015). Jarvis et al. (2014) indicated that ‘these computationally intensive analyses were conducted on more than 9 supercomputer centers and required the equivalent of > 400 years of computing using a single processor.’
Apart from developing new computational methods, analyses should be conducted more efficiently. For example, the newest computer hardware technology, such as graphics processing units (GPUs) and multicore central processing units (CPUs), can be implemented to parallelize calculations (Ayres et al. 2011). In addition, new techniques from Artificial Intelligence and deep learning could be applied to ‘text mine’ genomes which may yield cases of hybridization where we had not expected them (Fogel 2008; Angermueller et al. 2016; Leung et al. 2016).
In this review, we showed how genomic data can be applied to the study of avian introgression. First, the detection of hybrids and backcrosses has improved dramatically, although the monopoly of STRUCTURE should be broken up by applying other software in concert with it (e.g., DAPC or ADMIXTURE). Another way of detecting introgressive hybridization, phylogenetic discordance (i.e. different loci resulting in discordant gene trees), should be regarded as a starting point for further analyses, not as a definitive proof of introgression. Specifically, disentangling introgression from incomplete lineage sorting remains a challenging endeavour, although new techniques, such as the D-statistic, are being developed. Furthermore, with the advent of genomic data, phylogenetics might require a shift from trees to networks (Edwards et al. 2016; Ottenburghs et al. 2016).
The study of hybrid zones has led to important insights into the complex interplay between hybridization and speciation (Harrison and Larson 2016). Genomic data provide the opportunity to augment the resolution of geographical cline analysis. In addition, genomic cline analysis provides a fresh perspective on hybrid zone dynamics, circumventing several limitations of geographical cline analysis (Gompert and Buerkle 2009). It is, however, important to keep in mind that each hybrid zone study is just a single snapshot of a complex and continuously changing interaction. To capture the dynamic nature of hybrid zones, they should be studied across different temporal and spatial scales.
Genome scans, which uncovered a highly heterogeneous genomic landscape, have become a powerful tool in the genomic toolbox (Nosil and Feder 2012). When performing these genome-wide comparisons, one should realize that there is more to life than F ST. Other summary statistics provide different perspectives on the processes that sculpted the genomic landscape (Cruickshank and Hahn 2014; Wolf and Ellegren 2017). The debate which evolutionary processes underlie the genomic landscape, linked selection or reduced gene flow, is still ongoing. Also, the question whether loci involved in reproductive isolation cluster together in ‘islands of speciation’ or whether they are scattered throughout the genome remains to be answered. Exploring the genomic landscapes across the avian tree of life represents an exciting field for further research.
Finally, the findings from different methods should be incorporated into specific speciation scenarios, which can consequently be tested using a modelling approach. Especially, the application of ABC modelling, which allows for the comparison of different scenarios, will increase our understanding in avian speciation, which turns out to be more complex than the classical Mayrian triumvirate of allopatric, sympatric and parapatric speciation. Although there are many computation challenges ahead, this genomic perspective on avian hybridization and speciation will further our understanding in evolution in general.
JO conceived the idea and wrote the manuscripts. RHSK, PvH, SvW, RCY and HHTP commented on previous versions of the manuscript. All authors read and approved the final manuscript.
We are grateful to Leo Joseph and Peter Grant for their useful comments on early versions of this manuscript. We also thank the anonymous reviewer whose suggestions greatly improved the manuscript.
The authors declare that they have no competing interests.
This research was funded by Stichting de Eik.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Alcaide M, Scordato ES, Price TD, Irwin DE. Genomic divergence in a ring species complex. Nature. 2014;511:83–5.PubMedView ArticleGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Allendorf FW, Leary RF, Spruell P, Wenburg JK. The problems with hybrids: setting conservation guidelines. Trends Ecol Evol. 2001;16:613–22.View ArticleGoogle Scholar
- Anderson E, Thompson EA. A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002;160:1217–29.PubMedPubMed CentralGoogle Scholar
- Angermueller C, Pärnamaa T, Parts L, Stegle O. Deep learning for computational biology. Mol Syst Biol. 2016;12:878.PubMedPubMed CentralView ArticleGoogle Scholar
- Arnold ML. Evolution through genetic exchange. Oxford: Oxford University Press; 2006.Google Scholar
- Ayres DL, Darling A, Zwickl DJ, Beerli P, Holder MT, Lewis PO, Huelsenbeck JP, Ronquist F, Swofford DL, Cummings MP. BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics. Syst Biol. 2011;61:170–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldassarre DT, White TA, Karubian J, Webster MS. Genomic and morphological analysis of a semipermeable avian hybrid zone suggests asymmetrical introgression of a sexual signal. Evolution. 2014;68:2644–57.PubMedView ArticleGoogle Scholar
- Barilani M, Deregnaucourt S, Gallego S, Galli L, Mucci N, Piombo R, Puigcerver M, Rimondi S, Rodriguez-Teijeiro JD, Spano S, Randi E. Detecting hybridization in wild (Coturnix c. coturnix) and domesticated (Coturnix c. japonica) quail populations. Biol Conserv. 2005;126:445–55.View ArticleGoogle Scholar
- Barton NH. Gene flow past a cline. Heredity. 1979;43:333–9.View ArticleGoogle Scholar
- Barton NH. Multilocus clines. Evolution. 1983;37:454–71.PubMedView ArticleGoogle Scholar
- Barton NH, Gale KS. Genetic analysis of hybrid zones. In: Harrison RG, editor. Hybrid zones and the evolutionary process. New York: Oxford University Press; 1993. p. 13–45.Google Scholar
- Barton NH, Hewitt GM. Analysis of hybrid zones. Annu Rev Ecol Syst. 1985;16:113–48.View ArticleGoogle Scholar
- Beaumont MA. Approximate Bayesian computation in evolution and ecology. Annu Rev Ecol Evol Syst. 2010;41:379–406.View ArticleGoogle Scholar
- Beiko RG, Hamilton N. Phylogenetic identification of lateral genetic transfer events. BMC Evol Biol. 2006;6:15.PubMedPubMed CentralView ArticleGoogle Scholar
- Borge T, Lindroos K, Nadvornik P, Syvanen AC, Saetre GP. Amount of introgression in flycatcher hybrid zones reflects regional differences in pre and post-zygotic barriers to gene exchange. J Evol Biol. 2005;18:1416–24.PubMedView ArticleGoogle Scholar
- Burri R, Nater A, Kawakami T, Mugal CF, Olason PI, Smeds L, Suh A, Dutoit L, Bureš S, Garamszegi LZ. Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers. Genome Res. 2015;25:1656–65.PubMedPubMed CentralView ArticleGoogle Scholar
- Corander J, Waldmann P, Marttinen P, Sillanpää MJ. BAPS 2: enhanced possibilities for the analysis of genetic population structure. Bioinformatics. 2004;20:2363–9.PubMedView ArticleGoogle Scholar
- Cruickshank TE, Hahn MW. Reanalysis suggests that genomic islands of speciation are due to reduced diversity, not reduced gene flow. Mol Ecol. 2014;23:3133–57.PubMedView ArticleGoogle Scholar
- Darriba D, Flouri T, Stamatakis A. The state of software in evolutionary biology. bioRxiv. 2015;031930.Google Scholar
- Degnan JH, Rosenberg NA. Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol Evol. 2009;24:332–40.PubMedView ArticleGoogle Scholar
- Delmore KE, Hubner S, Kane NC, Schuster R, Andrew RL, Camara F, Guigo R, Irwin DE. Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation. Mol Ecol. 2015;24:1873–88.PubMedView ArticleGoogle Scholar
- Duchesne P, Bernatchez L. AFLPOP: a computer program for simulated and real population allocation, based on AFLP data. Mol Ecol Notes. 2002;2:380–3.View ArticleGoogle Scholar
- Durand EY, Patterson N, Reich D, Slatkin M. Testing for ancient admixture between closely related populations. Mol Biol Evol. 2011;28:2239–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Eaton DAR, Ree RH. Inferring phylogeny and introgression using RADseq data: an example from flowering plants (Pedicularis: Orobanchaceae). Syst Biol. 2013;62:689–706.PubMedPubMed CentralView ArticleGoogle Scholar
- Edwards SV, Potter S, Schmitt CJ, Bragg JG, Moritz C. Reticulation, divergence, and the phylogeography–phylogenetics continuum. PNAS. 2016;113:8025–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Elgvin TO, Trier CN, Tørresen OK, Hagen IJ, Lien S, Nederbragt AJ, Ravinet M, Jensen H, Sætre G-P. The genomic mosaicism of hybrid speciation. Sci Adv. 2017;3:e1602996.PubMedPubMed CentralView ArticleGoogle Scholar
- Ellegren H, Smeds L, Burri R, Olason PI, Backström N, Kawakami T, Kunstner A, Makinen H, Nadachowska-Brzyska K, Qvarnström A, Uebbing S, Wolf JBW. The genomic landscape of species divergence in Ficedula flycatchers. Nature. 2012;491:756–60.PubMedGoogle Scholar
- Eriksson A, Manica A. Effect of ancient population structure on the degree of polymorphism shared between modern human populations and ancient hominins. PNAS. 2012;109:13956–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Feder JL, Egan SP, Nosil P. The genomics of speciation-with-gene-flow. Trends Genet. 2012;28:342–50.PubMedView ArticleGoogle Scholar
- Fitzpatrick BM. Rates of evolution of hybrid inviability in birds and mammals. Evolution. 2004;58:1865–70.PubMedView ArticleGoogle Scholar
- Fitzpatrick BM. Alternative forms for genomic clines. Ecol Evol. 2013;3:1951–66.PubMedPubMed CentralView ArticleGoogle Scholar
- Fogel GB. Computational intelligence approaches for pattern discovery in biological systems. Brief Bioinform. 2008;9:307–16.PubMedView ArticleGoogle Scholar
- Frosch C, Kraus RHS, Angst C, Allgower R, Michaux J, Teubner J, Nowak C. The genetic legacy of multiple beaver reintroductions in central Europe. PLoS ONE. 2014;9:e97619.PubMedPubMed CentralView ArticleGoogle Scholar
- Fuchs J, Pons JM, Liu L, Ericson PGP, Couloux A, Pasquet E. A multi-locus phylogeny suggests an ancient hybridization event between Campephilus and melanerpine woodpeckers (Aves: Picidae). Mol Phylogenet Evol. 2013;67:578–88.PubMedView ArticleGoogle Scholar
- Gay L, Neubauer G, Zagalska-Neubauer M, Debain C, Pons JM, David P, Crochet PA. Molecular and morphological patterns of introgression between two large white-headed gull species in a zone of recent secondary contact. Mol Ecol. 2007;16:3215–27.PubMedView ArticleGoogle Scholar
- Gompert Z, Buerkle CA. A powerful regression-based method for admixture mapping of isolation across the genome of hybrids. Mol Ecol. 2009;18:1207–24.PubMedView ArticleGoogle Scholar
- Gompert Z, Buerkle CA. Bayesian estimation of genomic clines. Mol Ecol. 2011;20:2111–27.PubMedView ArticleGoogle Scholar
- Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai WW, Fritz MHY, Hansen NF, Durand EY, Malaspinas AS, Jensen JD, Marques-Bonet T, Alkan C, Prufer K, Meyer M, Burbano HA, Good JM, Schultz R, Aximu-Petri A, Butthof A, Hober B, Hoffner B, Siegemund M, Weihmann A, Nusbaum C, Lander ES, Russ C, Novod N, Affourtit J, Egholm M, Verna C, Rudan P, Brajkovic D, Kucan Z, Gusic I, Doronichev VB, Golovanova LV, Lalueza-Fox C, de la Rasilla M, Fortea J, Rosas A, Schmitz RW, Johnson PLF, Eichler EE, Falush D, Birney E, Mullikin JC, Slatkin M, Nielsen R, Kelso J, Lachmann M, Reich D, Paabo SA. Draft sequence of the Neandertal genome. Science. 2010;328:710–22.PubMedPubMed CentralView ArticleGoogle Scholar
- Grossen C, Seneviratne SS, Croll D, Irwin DE. Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact. Mol Ecol. 2016;25:4247–66.PubMedView ArticleGoogle Scholar
- Haasl RJ, Payseur BA. Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication. Mol Ecol. 2016;25:5–23.PubMedView ArticleGoogle Scholar
- Harr B. Genomic islands of differentiation between house mouse subspecies. Genome Res. 2006;16:730–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Harrison RG. Hybrid zones: windows on evolutionary process. Oxf Surv Evol Biol. 1990;7:69–129.Google Scholar
- Harrison RG, Larson EL. Hybridization, introgression, and the nature of species boundaries. J Hered. 2014;105:795–809.PubMedView ArticleGoogle Scholar
- Harrison RG, Larson EL. Heterogeneous genome divergence, differential introgression, and the origin and structure of hybrid zones. Mol Ecol. 2016;25:2454–66.PubMedPubMed CentralView ArticleGoogle Scholar
- Hewitt GM. Hybrid zones—natural laboratories for evolutionary studies. Trends Ecol Evol. 1988;3:158–67.PubMedView ArticleGoogle Scholar
- Hey J. Isolation with migration models for more than two populations. Mol Biol Evol. 2010;27:905–20.PubMedView ArticleGoogle Scholar
- Hey J, Nielsen R. Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics. 2004;167:747–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Hey J, Chung YJ, Sethuraman A. On the occurrence of false positives in tests of migration under an isolation-with-migration model. Mol Ecol. 2015;24:5078–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–67.PubMedView ArticleGoogle Scholar
- Jarvis ED. Perspectives from the avian phylogenomics project: questions that can be answered with sequencing all genomes of a vertebrate class. Annu Rev Anim Biosci. 2016;4:45–59.PubMedView ArticleGoogle Scholar
- Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, Ho SYW, Faircloth BC, Nabholz B, Howard JT, Suh A, Weber CC, da Fonseca RR, Li JW, Zhang F, Li H, Zhou L, Narula N, Liu L, Ganapathy G, Boussau B, Bayzid MS, Zavidovych V, Subramanian S, Gabaldon T, Capella-Gutierrez S, Huerta-Cepas J, Rekepalli B, Munch K, Schierup M, Lindow B, Warren WC, Ray D, Green RE, Bruford MW, Zhan XJ, Dixon A, Li SB, Li N, Huang YH, Derryberry EP, Bertelsen MF, Sheldon FH, Brumfield RT, Mello CV, Lovell PV, Wirthlin M, Schneider MPC, Prosdocimi F, Samaniego JA, Velazquez AMV, Alfaro-Nunez A, Campos PF, Petersen B, Sicheritz-Ponten T, Pas A, Bailey T, Scofield P, Bunce M, Lambert DM, Zhou Q, Perelman P, Driskell AC, Shapiro B, Xiong ZJ, Zeng YL, Liu SP, Li ZY, Liu BH, Wu K, Xiao J, Yinqi X, Zheng QM, Zhang Y, Yang HM, Wang J, Smeds L, Rheindt FE, Braun M, Fjeldsa J, Orlando L, Barker FK, Jonsson KA, Johnson W, Koepfli KP, O’Brien S, Haussler D, Ryder OA, Rahbek C, Willersle E, Graves GR, Glenn TC, McCormack J, Burt D, Ellegren H, Alström P, Edwards SV, Stamatakis A, Mindell DP, Cracraft J, Braun EL, Warnow T, Jun W, Gilbert MTP, Zhang GJ. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346:1320–31.PubMedPubMed CentralView ArticleGoogle Scholar
- Joly S, McLenachan PA, Lockhart PJ. A statistical approach for distinguishing hybridization and incomplete lineage sorting. Am Nat. 2009;174:E54–70.PubMedView ArticleGoogle Scholar
- Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.PubMedPubMed CentralView ArticleGoogle Scholar
- Kingston SE, Parchman TL, Gompert Z, Buerkle CA, Braun MJ. Heterogeneity and concordance in locus-specific differentiation and introgression between species of towhees. J Evol Biol. 2017;30:474–85.PubMedView ArticleGoogle Scholar
- Kozlov AM, Aberer AJ, Stamatakis A. ExaML version 3: a tool for phylogenomic analyses on supercomputers. Bioinformatics. 2015;31:2577–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Kraus RH, Wink M. Avian genomics: fledging into the wild! J Ornithol. 2015;156:1–15.View ArticleGoogle Scholar
- Kraus RHS, Kerstens HHD, van Hooft P, Megens HJ, Elmberg J, Tsvey A, Sartakov D, Soloviev SA, Crooijmans RPMA, Groenen MAM, Ydenberg RC, Prins HHT. Widespread horizontal genomic exchange does not erode species barriers among sympatric ducks. BMC Evol Biol. 2012;12:45.PubMedPubMed CentralView ArticleGoogle Scholar
- Kubatko LS. Identifying hybridization events in the presence of coalescence via model selection. Syst Biol. 2009;58:478–88.PubMedView ArticleGoogle Scholar
- Lackey AC, Boughman JW. Evolution of reproductive isolation in stickleback fish. Evolution. 2017;71:357–72.PubMedView ArticleGoogle Scholar
- Lamichhaney S, Berglund J, Almen MS, Maqbool K, Grabherr M, Martinez-Barrio A, Promerova M, Rubin CJ, Wang C, Zamani N, Grant BR, Grant PR, Webster MT, Andsersson L. Evolution of Darwin’s finches and their beaks revealed by genome sequencing. Nature. 2015;518:371–5.PubMedView ArticleGoogle Scholar
- Lavretsky P, Engilis A, Eadie JM, Peters JL. Genetic admixture supports an ancient hybrid origin of the endangered Hawaiian duck. J Evol Biol. 2015;28:1005–15.PubMedView ArticleGoogle Scholar
- Leung MK, Delong A, Alipanahi B, Frey BJ. Machine learning in genomic medicine: a review of computational problems and data sets. Proc IEEE. 2016;104:176–97.View ArticleGoogle Scholar
- Linck EB, Hanna ZR, Sellas A, Dumbacher JP. Evaluating hybridization capture with RAD probes as a tool for museum genomics with historical bird specimens. Ecol Evol. 2017;7:4755–67.PubMedPubMed CentralView ArticleGoogle Scholar
- Maddison WP. Gene trees in species trees. Syst Biol. 1997;46:523–36.View ArticleGoogle Scholar
- Maddison WP, Knowles LL. Inferring phylogeny despite incomplete lineage sorting. Syst Biol. 2006;55:21–30.PubMedView ArticleGoogle Scholar
- Martin SH, Davey JW, Jiggins CD. Evaluating the use of ABBA–BABA statistics to locate introgressed loci. Mol Biol Evol. 2014;32:244–57.PubMedPubMed CentralView ArticleGoogle Scholar
- Meng C, Kubatko LS. Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: a model. Theor Popul Biol. 2009;75:35–45.PubMedView ArticleGoogle Scholar
- Nadachowska-Brzyska K, Li C, Smeds L, Zhang GJ, Ellegren H. Temporal dynamics of avian populations during Pleistocene revealed by whole-genome sequences. Curr Biol. 2015;25:1375–80.PubMedPubMed CentralView ArticleGoogle Scholar
- Nater A, Burri R, Kawakami T, Smeds L, Ellegren H. Resolving evolutionary relationships in closely related species with whole-genome sequencing data. Syst Biol. 2015;64:1000–17.PubMedPubMed CentralView ArticleGoogle Scholar
- Nielsen EE, Bach LA, Kotlicki P. HYBRIDLAB (version 1.0): a program for generating simulated hybrids from population samples. Mol Ecol Notes. 2006;6:971–3.View ArticleGoogle Scholar
- Nosil P, Feder JL. Widespread yet heterogeneous genomic divergence. Mol Ecol. 2012;21:2829–32.PubMedView ArticleGoogle Scholar
- Ottenburghs J, Ydenberg RC, van Hooft P, van Wieren SE, Prins HHT. The avian hybrids project: gathering the scientific literature on avian hybridization. Ibis. 2015;157:892–4.View ArticleGoogle Scholar
- Ottenburghs J, van Hooft P, Van Wieren SE, Ydenberg RC, Prins HHT. Birds in a bush: toward an avian phylogenetic network. Auk. 2016;133:577–82.View ArticleGoogle Scholar
- Ottenburghs J, Megens H-J, Kraus RHS, van Hooft P, van Wieren SE, Crooijmans RPMA, Ydenberg RC, Groenen MAM, Prins HHT. A history of hybrids? Genomic patterns of introgression in the True Geese. BMC Evol Biol. 2017;17:201.PubMedPubMed CentralView ArticleGoogle Scholar
- Oyler-McCance SJ, Oh KP, Langin KM, Aldridge CL. A field ornithologist’s guide to genomics: practical considerations for ecology and conservation. Auk. 2016;133:626–48.View ArticleGoogle Scholar
- Pamilo P, Nei M. Relationships between gene trees and species trees. Mol Biol Evol. 1988;1988(5):568–83.Google Scholar
- Parchman TL, Gompert Z, Braun MJ, Brumfield RT, Mcdonald DB, Uy JAC, Zhang G, Jarvis ED, Schlinger BA, Buerkle CA. The genomic consequences of adaptive divergence and reproductive isolation between species of manakins. Mol Ecol. 2013;22:3304–17.PubMedView ArticleGoogle Scholar
- Parsons TJ, Olson SL, Braun MJ. Unidirectional spread of secondary sexual plumage traits across an avian hybrid zone. Science. 1993;260:1643–6.PubMedView ArticleGoogle Scholar
- Payseur BA. Using differential introgression in hybrid zones to identify genomic regions involved in speciation. Mol Ecol Resour. 2010;10:806–20.PubMedView ArticleGoogle Scholar
- Payseur BA, Rieseberg LH. A genomic perspective on hybridization and speciation. Mol Ecol. 2016;25:2337–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Pease JB, Hahn MW. Detection and polarization of introgression in a five-taxon phylogeny. Syst Biol. 2015;64:651–62.PubMedView ArticleGoogle Scholar
- Pinho C, Hey J. Divergence with gene flow: models and data. Annu Rev Ecol Evol Syst. 2010;41:215–30.View ArticleGoogle Scholar
- Poelstra JW, Vijay N, Bossu CM, Lantz H, Ryll B, Muller I, Baglione V, Unneberg P, Wikelski M, Grabherr MG, Wolf JBW. The genomic landscape underlying phenotypic integrity in the face of gene flow in crows. Science. 2014;344:1410–4.PubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics. 2014;197:573–89.PubMedPubMed CentralView ArticleGoogle Scholar
- Raposo do Amaral F, Albers PK, Edwards SV, Miyaki CY. Multilocus tests of Pleistocene refugia and ancient divergence in a pair of Atlantic forest antbirds (Myrmeciza). Mol Ecol. 2013;22:3996–4013.PubMedView ArticleGoogle Scholar
- Rheindt FE, Edwards SV. Genetic introgression: an integral but neglected component of speciation in birds. Auk. 2011;128:620–32.View ArticleGoogle Scholar
- Rheindt FE, Fujita MK, Wilton PR, Edwards SV. Introgression and phenotypic assimilation in Zimmerius flycatchers (Tyrannidae): population genetic and phylogenetic inferences from genome-wide SNPs. Syst Biol. 2014;63:134–52.PubMedView ArticleGoogle Scholar
- Ruegg K, Anderson EC, Boone J, Pouls J, Smith TB. A role for migration-linked genes and genomic islands in divergence of a songbird. Mol Ecol. 2014;23:4757–69.PubMedView ArticleGoogle Scholar
- Saetre GP, Borge T, Lindell J, Moum T, Primmer CR, Sheldon BC, Haavie J, Johnsen A, Ellegren H. Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers. Mol Ecol. 2001;10:737–49.PubMedView ArticleGoogle Scholar
- Saetre GP, Borge T, Lindroos K, Haavie J, Sheldon BC, Primmer C, Syvanen AC. Sex chromosome evolution and speciation in Ficedula flycatchers. Proc R Soc B. 2003;270:53–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Schaefer J, Duvernell D, Campbell DC. Hybridization and introgression in two ecologically dissimilar Fundulus hybrid zones. Evolution. 2016;70:1051–63.PubMedView ArticleGoogle Scholar
- Seneviratne SS, Davidson P, Martin K, Irwin DE. Low levels of hybridization across two contact zones among three species of woodpeckers (Sphyrapicus sapsuckers). J Avian Biol. 2016;47:887–98.View ArticleGoogle Scholar
- Smyth JF, Patten MA, Pruett CL. The evolutionary ecology of a species ring: a test of alternative models. Folia Zool. 2015;64:233–44.Google Scholar
- Spurgin LG, Wright DJ, van der Velde M, Collar NJ, Komdeur J, Burke T, Richardson DS. Museum DNA reveals the demographic history of the endangered Seychelles warbler. Evol Appl. 2014;7:1134–43.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor SA, White TA, Hochachka WM, Ferretti V, Curry RL, Lovette I. Climate-mediated movement of an avian hybrid zone. Curr Biol. 2014;24:671–6.PubMedView ArticleGoogle Scholar
- Tegelström H, Gelter HP. Haldane rule and sex biased gene flow between two hybridizing flycatcher species (Ficedula albicollis and F. hypoleuca, Aves, Muscicapidae). Evolution. 1990;44:2012–21.PubMedGoogle Scholar
- Than C, Nakhleh L. Species tree inference by minimizing deep coalescences. PLoS Comput Biol. 2009;5:e1000501.PubMedPubMed CentralView ArticleGoogle Scholar
- Toews DPL, Campagna L, Taylor SA, Balakrishnan CN, Baldassarre DT, Deane-Coe PE, Harvey MG, Hooper DM, Irwin DE, Judy CD, Mason NA, McCormack JE, McCracken KG, Oliveros CH, Safran RJ, Scordato ESC, Faust Stryjweksi K, Tigano A, Uy JAC, Winger BM. Genomic approaches to understanding population divergence and speciation in birds. Auk. 2016;133:13–30.View ArticleGoogle Scholar
- Turner TL, Hahn MW, Nuzhdin SV. Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005;3:1572–8.View ArticleGoogle Scholar
- Via S. Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow. Philos Trans R Soc Lond B Biol. 2012;367:451–60.View ArticleGoogle Scholar
- Vijay N, Bossu CM, Poelstra JW, Weissensteiner MH, Suh A, Kryukov AP, Wolf JBW. Evolution of heterogeneous genome differentiation across multiple contact zones in a crow species complex. Nat Commun. 2016;7:13195.PubMedPubMed CentralView ArticleGoogle Scholar
- Walsh J, Rowe RJ, Olsen BJ, Shriver WG, Kovach AI. Genotype-environment associations support a mosaic hybrid zone between two tidal marsh birds. Ecol Evol. 2016;6:279–94.PubMedView ArticleGoogle Scholar
- Wang LY, Luzynski K, Pool JE, Janousek V, Dufkova P, Vyskocilova MM, Teeter KC, Nachman MW, Munclinger P, Macholan M, Pialek J, Tucker PK. Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone. Mol Ecol. 2011;20:2985–3000.PubMedView ArticleGoogle Scholar
- Wen D, Yu Y, Nakhleh L. Bayesian inference of reticulate phylogenies under the multispecies network coalescent. PLoS Genet. 2016;12:e1006006.PubMedPubMed CentralView ArticleGoogle Scholar
- Wolf JB, Ellegren H. Making sense of genomic islands of differentiation in light of speciation. Nat Rev Genet. 2017;18:87–100.PubMedView ArticleGoogle Scholar
- Wu CI. The genic view of the process of speciation. J Evol Biol. 2001;14:851–65.View ArticleGoogle Scholar
- Yeung CKL, Tsai PW, Chesser RT, Lin RC, Yao CT, Tian XH, Li SH. Testing founder effect speciation: divergence population genetics of the spoonbills Platalea regia and P. minor (Threskiornithidae, Aves). Mol Biol Evol. 2011;28:473–82.PubMedView ArticleGoogle Scholar
- Yu Y, Barnett RM, Nakhleh L. Parsimonious inference of hybridization in the presence of incomplete lineage sorting. Syst Biol. 2013;62:738–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Zarza E, Faircloth BC, Tsai WL, Bryson RW Jr, Klicka J, McCormack JE. Hidden histories of gene flow in highland birds revealed with genomic markers. Mol Ecol. 2016;25:5144–57.PubMedView ArticleGoogle Scholar