- Research article
- Open Access
Patterns of gene flow and selection across multiple species of Acrocephalus warblers: footprints of parallel selection on the Z chromosome
© The Author(s). 2016
Received: 10 March 2016
Accepted: 25 May 2016
Published: 16 June 2016
Understanding the mechanisms and selective forces leading to adaptive radiations and origin of biodiversity is a major goal of evolutionary biology. Acrocephalus warblers are small passerines that underwent an adaptive radiation in the last approximately 10 million years that gave rise to 37 extant species, many of which still hybridize in nature. Acrocephalus warblers have served as model organisms for a wide variety of ecological and behavioral studies, yet our knowledge of mechanisms and selective forces driving their radiation is limited. Here we studied patterns of interspecific gene flow and selection across three European Acrocephalus warblers to get a first insight into mechanisms of radiation of this avian group.
We analyzed nucleotide variation at eight nuclear loci in three hybridizing Acrocephalus species with overlapping breeding ranges in Europe. Using an isolation-with-migration model for multiple populations, we found evidence for unidirectional gene flow from A. scirpaceus to A. palustris and from A. palustris to A. dumetorum. Gene flow was higher between genetically more closely related A. scirpaceus and A. palustris than between ecologically more similar A. palustris and A. dumetorum, suggesting that gradual accumulation of intrinsic barriers rather than divergent ecological selection are more efficient in restricting interspecific gene flow in Acrocephalus warblers. Although levels of genetic differentiation between different species pairs were in general not correlated, we found signatures of apparently independent instances of positive selection at the same two Z-linked loci in multiple species.
Our study brings the first evidence that gene flow occurred during Acrocephalus radiation and not only between sister species. Interspecific gene flow could thus be an important source of genetic variation in individual Acrocephalus species and could have accelerated adaptive evolution and speciation rate in this avian group by creating novel genetic combinations and new phenotypes. Independent instances of positive selection at the same loci in multiple species indicate an interesting possibility that the same loci might have contributed to reproductive isolation in several speciation events.
Interspecific gene flow is an important evolutionary force. It may enrich genetic variation of individual species, facilitate the origin of new phenotypes, lead to the origin of new species, or, conversely, to species fusion . It has been suggested that high speciation rates and rapid phenotypic changes observed during adaptive radiations may be facilitated by increased genetic variation and novel genetic combinations produced by introgressive hybridization among the species [2, 3]. Patterns of interspecific gene flow can also provide an important insight into the genetic basis of reproductive isolation . In species with incomplete reproductive isolation, gene flow occurs at loci with alleles that are neutral or beneficial on the genomic background of the other species, but is limited at loci harboring alleles involved in reproductive isolation .
Most work exploring patterns and rates of interspecific gene flow focused on pairs of closely related sister species  because, until recently , models allowing the analysis of interspecific gene flow were limited to only two species. Several general findings emerged from these studies. First, in the case of speciation driven by intrinsic barriers, interspecific gene flow is often reduced on the sex chromosomes (i.e., XY in heterogametic male organisms and ZW in heterogametic female organisms) [8–11] and in genomic regions with low recombination rates , indicating that genetic incompatibilities tend to preferentially accumulate in these genomic regions. Second, in the case of speciation driven by divergent ecological selection, gene flow is correlated with ecological divergence  and is reduced around genes underlying ecological differences . Interspecific gene flow across larger sections of the phylogenetic tree is, however, much less studied. Such studies would be particularly interesting in instances of rapid radiations, where lineage divergence occurs on a shorter timescale than the completion of reproductive isolation [15–18]. Knowledge of patterns and rates of interspecific gene flow in such systems allows to address questions, which are difficult or even impossible to study using models of only two sister species. One can, for example, examine whether the same genes contribute to reproductive isolation in different speciation events, assess the importance of various reproductive barriers in restricting the gene flow and determine the order in which different reproductive barriers arise during lineage divergence. This information can help us to understand why lineage divergence occurs so quickly in some taxonomic groups, while not in others.
Acrocephalus warblers are small passerines that underwent an adaptive radiation during the last approximately 10 million years , which gave rise to 37 extant species occupying mainly Eurasia, Africa, Australia and Pacific islands . The members of this genus are phenotypically quite uniform. They, however, display a great diversity in ecology and behaviour, and have been used extensively as models in ecological and behavioral research . Interspecific hybridization occasionally occurs in Acrocephalus warblers. It has been reported not only between sister species (e.g., A. scirpaceus × A. palustris, A. arundinaceus × A. stentoreus), but also between more distantly related species (e.g., A. palustris × A. dumetorum, A. palustris × A. schoenobaenus, A. scirpaceus × A. arundinaceus) (see Additional file 1 for all reported cases of interspecific hybridization in Acrocephalus warblers and references). In some of these cases, hybridization occurs between species of very different body size (A. scirpaceus × A. arundinaceus) or even between species belonging to different subgenera (e.g. A. palustris × A. schoenobaenus). So far very little is known about the degree of postzygotic reproductive isolation between the species. In some cases, e.g. between sister species A. arundinaceus and A. stentoreus, postzygotic isolation seems to be quite strong as no backcross hybrids have been observed in a hybridizing population . In other cases (e.g. between A. palustris and A. dumetorum), however, observation of fertile F1 hybrids  suggests that postzygotic isolation is incomplete and interspecific gene flow might occur between the species.
Here we performed the first study exploring patterns and rates of interspecific gene flow in this avian group, focusing on three species of the subgenus Notiocichla: A. scirpaceus, A. palustris and A. dumetorum. All three species are very similar in morphology, but have different habitat requirements, migration strategies and song patterns . A. scirpaceus breeds across Europe, northern Africa and western Asia typically in reed beds and winters in sub-Saharan Africa. Its sister species, A. palustris, has similar geographical distribution, but occupies mostly damp herbaceous vegetation and winters in southeastern Africa. Males of each species sing different songs, nevertheless, in A. palustris, mixed singers that incorporate song phrases from A. scirpaceus to its repertoire are known . Both species occasionally form mixed pairs  and are able to produce viable hybrid offspring (e.g. ; see Additional file 1 for a complete list of references). A. dumetorum, the sister species to the A. scirpaceus and A. palustris clade , breeds in Asia and Northeastern Europe in various herbaceous or bushy vegetation and winters in India. In Eastern Europe and Western Asia, it often co-occurs with A. palustris at the same sites and hybridization between the two species has been observed (e.g. [22, 27]; see Additional file 1 for a complete list of references). The exact frequency of hybridization between A. palustris and A. scirpaceus and between A. palustris and A. dumetorum is hard to estimate owing to morphological similarity of the species and difficulty to distinguish hybrids without genetic analysis. The total number of observed mixed pairs is, however, comparable in both species pairs (Additional file 1). The cases of hybridization between A. dumetorum and A. scirpaceus, which are ecologically more divergent, are not known, although one putative hybrid individual has been reported .
To get the first insight into mechanisms and selective forces driving Acrocephalus radiation, we analyzed nucleotide variation at four autosomal and four loci on the Z chromosome in sympatric populations of A. scirpaceus, A. palustris and A. dumetorum. Using an isolation-with-migration model for multiple populations we tested whether interspecific gene flow occurred among the studied species and if so, whether it was higher between more genetically similar or between more ecologically similar species. In addition, using several tests of neutrality, we explored patterns of positive selection in all three species.
Study area and samples
Molecular sex determination and sequencing
Genomic DNA was purified by DNeasy Tissue Kit (Qiagen) according to manufacturer’s instructions. Sex of individual birds was determined following the method of Griffiths et al.  using P2 and P8 primers. These primers amplify a fragment of two homologous genes located on the Z and W chromosomes. The amplified fragments from the Z and W chromosomes differ in length, which can be visualized on agarose gel. We further amplified and sequenced introns of four autosomal and four Z-linked loci. The four autosomal loci are located on different chromosomes in the Zebra Finch (Taeniopygia guttata) genome. Of the four Z-linked genes, PPWD1 and ADAMTS6 are close physically (68 kb), and therefore may not have independent evolutionary histories. However, the other two loci on the Z chromosome are separated from each other and from PPWD1 and ADAMTS6 by > 19 Mb. Primers for PCR amplification were partially obtained from published studies and partially designed by ourselves. All primers were designed in conserved exonic regions of the Zebra Finch or chicken (Gallus gallus) genome in order to amplify intronic sequences. Primer sequences and lengths of the obtained PCR products are provided in Additional file 3. PCR conditions are same as in Storchová et al. . All PCR products were sequenced in both directions with the PCR primers using Sanger sequencing at the DNA sequencing laboratory at Charles University in Prague.
Polymorphism statistics for eight studied loci
Basic population genetic analyses of polymorphism, divergence, recombination, and tests of neutrality based on the allele frequency spectrum were performed with the program DnaSP . Hudson-Kreitman-Aguadé (HKA) tests of positive selection  were performed using the HKA program (https://bio.cst.temple.edu/~hey/software/software.htm#HKA). This program compares the ratio of polymorphism to divergence at multiple loci. Loci under positive selection are expected to show lower levels of polymorphism relative to divergence. For each locus, we set the information whether the locus is autosomal or Z-linked, based on which the HKA program accounts for lower effective population size and thus the lower nucleotide variation on the Z chromosome. Genealogical relationships among haplotypes for each locus were reconstructed with Network software  using Median joining algorithm.
The data were fit to the isolation-with-migration model (IM) for multiple populations  using the program IMa2. The program estimates several demographic parameters based on Markov chain Monte Carlo simulations of genealogies. These parameters include the effective population size for the current and ancestral populations, migration rates between the populations and the population-split times. Because IMa2 assumes no recombination within loci, we determined the longest region without observed recombination for each locus using the program IMgc . This program removes either sites or haplotypes to produce the most information-rich contiguous DNA sequence segment that passes the four-gamete test. Nonrecombinant regions represented 84 % of the length of each locus, on average, and were used as an input for the IMa2 program. For each locus, we provided the information whether it is autosomal or Z-linked, based on which IMa2 program accounts for different effective population sizes between autosomes and the Z chromosome. We ran the program three times with identical starting conditions, with the exception of the random number seed, to assess convergence. To facilitate mixing of the Markov chains, we used Metropolis coupling with 30 chains and a geometric heating model. Upper bounds for the prior distributions of parameter values were set as suggested in IMa2 documentation (10 for population size parameters, 3 for migration rates parameters and 5 for splitting time). All runs began with a burn-in period of 100,000 steps and were allowed to continue for 7–8 million steps. We were able to achieve adequate mixing of the Markov chains as indicated by trend line plots, effective sample size (ESS) values (for all three runs, most ESS values were higher than 20,000 and no ESS value was lower than 50), and very similar parameter estimates in the first half and the second half of the run. Independent runs converged to the same result (e.g., maximum-likelihood estimates and marginal posterior probability distributions for the demographic parameters were essentially the same for all three runs). To perform likelihood-ratio tests of nested models, we combined results of the three independent IMa2 runs (together containing more than 200,000 genealogies) in a single L mode run.
All estimated parameters of the IM model are scaled to the mutation rate. To convert these parameters to biologically meaningful quantities (i.e., N e , effective population size in number of individuals; m, migration rate per year; 2 Nm, population migration rate; t, divergence time in years) we calculated the neutral mutation rate for each locus (μ) using divergence to the outgroup. The neutral mutation rate per year was calculated for each locus by using the formula D = 2 μt 1 where D is the average pairwise divergence, Dxy, between the three studied species and the outgroup, and t 1 is the divergence time. We assumed that the divergence time between A. scirpaceus/A. palustris/A. dumetorum and A. schoenobaenus is 5.9 million years. This was estimated on the basis of cytochrome b sequence divergence, which is 12.4 % , and assuming approximately 2.1 % sequence divergence per million years . The geometric mean of the locus-specific mutation rates was then used for the parameter conversion. A generation time of one year (Acrocephalus warblers reproduce only once per year) was assumed to estimate mutation rates per generation.
Levels of intraspecific polymorphism and tests of neutrality
HKA test of positive selection
Sum of deviationsb
We also looked for evidence of selection by comparing the distribution of allele frequencies with the expectations under a neutral equilibrium model using Tajima’s D and Fu-Li’s D tests (Table 1). These tests were performed for all studied loci separately in each species. There was a tendency toward negative values of Tajima’s D in all three species (17 of 19 values were negative), which may reflect mild population expansions. Nonetheless, none of the 19 values of Tajima’s D within species were significantly different from the neutral expectation of 0. No significant deviations from neutral expectations were either found in Fu-Li’s D tests. These tests, however, could not be performed for loci with zero within-species nucleotide variation.
Estimation of divergence times, effective population sizes and rates of interspecific gene flow
Maximum-likelihood estimates (MLE) and 95 % highest posterior density (HPD) intervals of demographic parameters of isolation with migration model
Log-likelihood ratio tests of nested models
ms>p = 0, mp>s = 0
ms>d = 0, md>s = 0
mp>d = 0, md>p = 0
mp>sp = 0, msp>d = 0
We were further interested in whether patterns of interspecific gene flow between A. scirpaceus and A. palustris, and between A. palustris and A. dumetorum are similar or not. As IMa2 does not allow to estimate locus-specific migration rates, we compared FST values, which reflect between-species differentiation and are correlated with levels of gene flow in hybridizing taxa (e.g., ). We should, however, note that FST values might be also affected by the action of positive selection that leads to a reduced within-species nucleotide variation . No statistically significant correlation in FST between the two species pairs was found (r = 0.47, p > 0.05), suggesting that the patterns of interspecific gene flow and/or positive selection are in general different for the two species pairs.
Adaptive radiations represent an outstanding system for studying mechanisms of speciation and adaptive divergence. Although divergent ecological selection seems to be an important driver of many adaptive radiations , it has been suggested that special genomic properties, such as gene duplications, activation of transposable elements or presence of gene flow among species, can also contribute to high speciation rates and rapid phenotypic changes during adaptive radiations .
Acrocephalus warblers represent an interesting model system for studying mechanisms of radiation in birds, although their mainly continental radiation was relatively slow compared to the well-known cases of adaptive radiations on islands, such as Darwin finches [13, 17]. Our estimates of divergence times are approximately 1.1 Mya for the sister species A. scirpaceus and A. palustris, and 2.5 Mya for their common ancestor and A. dumetorum. Despite this relatively deep divergence, many currently existing Acrocephalus species still hybridize in nature (Additional file 1). This hybridization might in principle lead to gene flow among the species if F1 hybrids are fertile. However, since very little is known about the degree of postzygotic isolation in this avian group, it is hard to predict whether gene flow could occur among Acrocephalus species and how strong it might be.
Here we provided the first evidence that interspecific gene flow occurred during Acrocephalus radiation and not only between sister species. We found evidence of gene flow between A. scirpaceus and A. palustris, which are the sister species with somewhat different ecological requirements, and between A. palustris and A. dumetorum, which have similar ecological requirements and often co-occur at the same sites where their ranges overlap. Significant gene flow was detected also between A. dumetorum and the ancestor of A. scirpaceus and A. palustris. No gene flow was, however, detected between A. scirpaceus and A. dumetorum. Although the isolation-with-migration model estimates not only contemporary gene flow, but average gene flow among populations since the time of divergence , it is noteworthy that hybridization currently occurs between both pairs of species for which gene flow has been detected, but not between A. scirpaceus and A. palustris for which no clear cases of interspecific hybridization are known.
Interestingly, the estimated levels of gene flow were in all cases asymmetric. Gene flow occurred from A. scirpaceus to A. palustris (2 Nm = 0.238), from A. palustris to A. dumetorum (2 Nm = 0.062), and from the A. scirpaceus/A. palustris ancestor to A. dumetorum (2 Nm = 0.95). Zero gene flow was detected in the opposite directions. Such unidirectional gene flow is expected when a population of one species expands to the area already occupied by a related species and interbreeding is not prevented between the two species. Introgression of neutral alleles then occurs almost exclusively from the local into the invading species . Unfortunately, history of breeding range changes are not well known for the studied species, but historically recent breeding range expansion of A. dumetorum to eastern and northern Europe where it encountered A. palustris  is consistent with this scenario.
An important issue in speciation research is to understand evolutionary forces that are responsible for the establishment of reproductive barriers between the species. Here we tested whether gene flow is higher between genetically more similar, but ecologically somewhat different, A. scirpaceus and A. palustris, or between ecologically more similar, but genetically more distant, A. palustris and A. dumetorum. Our results show that gene flow is more than three times higher between genetically more similar A. scirpaceus and A. palustris than between A. palustris and A. dumetorum, although the latter two species have higher chance to meet in the breeding sites and very likely more often hybridize. Lower levels of gene flow between A. palustris and A. dumetorum might be caused by historically more recent contact between the species , but might also suggest that divergent ecological selection is not the major cause of reproductive isolation and that speciation in Acrocephalus warblers is more likely driven by gradual accumulation of postzygotic barriers, such as hybrid sterility. Hybrid sterility have been thought to evolve quite slowly in birds compared to other taxonomic groups; however, recent findings suggest that at least in some bird species it can evolve very rapidly . It is, however, also possible that very different migratory behaviors of A. dumetorum, wintering in India, and A. palustris, that winters similarly as A. scirpaceus in Africa (although in different regions), might have contributed to their stronger reproductive isolation. Indeed, studies in Swainson’s Thrush (Catharus ustulatus) show that different seasonal migratory behavior of its two subspecies can be a source of selection against hybrids . Analysis of more species with different genetic and ecological similarity would be needed to conclude which selective forces drive radiation of Acrocephalus warblers.
Although the patterns of between species differentiation were in general different for the two species pairs, we detected signatures of apparently independent instances of positive selection at the same two Z-linked loci in multiple species. Selective sweeps reducing the within-species nucleotide variation occurred independently in A. scirpaceus and A. palustris at TG401 locus and even in all three species at TG1505 locus. Although we do not know the physical location of the two loci in the Acrocephalus genome, the two loci lie in different parts of the Z chromosome in the Zebra Finch genome. Given the very conserved karyotypes in birds , it is likely that the two loci are also unlinked in Acrocephalus genome and that selection has led to reduced nucleotide variation independently at the two loci. Similar examples of reduced nucleotide diversity at the same loci in multiple species have been observed for example in a genomic scan of islands of divergence in two hybridizing Ficedula flycatchers  or between two closely related and hybridizing plants, Silene latifolia and S. dioica . Such independent instances of positive selection in multiple species suggest that the same (or closely linked) loci might have been subject to parallel selection in two or even all three Acrocephalus species. This might happen for example in genes involved in intragenomic conflicts, such as meiotic drive, which is expected to be especially common on the sex chromosomes. Whether divergence of these genes could contribute to Acrocephalus speciation, however, still remains to be explored. In addition, it should be tested whether similar patterns of adaptive evolution in the three Acrocephalus species could not result from linked selection in low-recombinant regions rather than parallel adaptive evolution .
Our study provides the first evidence that interspecific gene flow occurred during Acrocephalus radiation and not only between sister species. Interspecific gene flow could thus play an important role in shaping genetic variation in individual species and possibly accelerate adaptive evolution and speciation rates in this group. Interestingly, the gene flow was in all cases asymmetric, suggesting that it might have occurred mainly during range expansions of individual species into areas of closely related and not completely reproductively isolated species. Interspecific competition leading to ecological character displacement and habitat segregation  and/or reinforcement of prezygotic isolation  could then strengthen reproductive barriers between the species. However, in contrast to Darwin finches , our results showed that gene flow was higher between genetically more closely related but ecologically more divergent species suggesting that gradual accumulation of genetic incompatibilities causing intrinsic isolation rather than divergent ecological selection are more efficient in restricting interspecific gene flow in this avian group. Thus in contrast to Darwin finches , the ecological differentiation of Acrocephalus species likely arose in later stages of speciation. Data on more than three Acrocephalus species would be, however, needed to generalize this conclusion to the whole Acrocephalus radiation. In addition, genome-wide data across autosomes and the whole Z chromosome would be needed to understand mechanisms responsible for the parallel adaptive evolution at the same loci in multiple species and its possible role in Acrocephalus speciation. Nevertheless, occurrence of interspecific hybridization among many Acrocephalus species suggests that interspecific gene flow could be common across the whole radiation.
We thank Milica Požgayová for the assistance in the field. The study was supported by junior grant of the Czech Science Foundation (15-10884Y) to RR and PP was supported by the Institutional Research Plan (RVO: 68081766).
Availability of data and material
Sequences have been deposited in GenBank under accession numbers KU094074-KU094575.
RR an PP conceived the study. RR, PP and JR wrote the manuscript. JR, PP, MA, AL collected samples. VM sequenced samples. VM and RR analyzed data. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Manipulation of birds was approved by the National Board on Experimental Animal Research, Helsinki, Finland (permission no. ESAVI-2010-05480/Ym-23).
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.
- Arnold ML. Divergence with genetic exchange. Oxford: Oxford University Press; 2015.View ArticleGoogle Scholar
- Seehausen O. Hybridization and adaptive radiation. Trends Ecol Evol. 2004;19:198–207.View ArticlePubMedGoogle Scholar
- Palmer DH, Kronforst MR. Divergence and gene flow among Darwin’s finches: a genome-wide view of adaptive radiation driven by interspecies allele sharing. Bioessays. 2015;37:968–74.View ArticlePubMedGoogle Scholar
- Seehausen O, Butlin RK, Keller I, Wagner CE, Boughman JW, Hohenlohe PA, et al. Genomics and the origin of species. Nat Rev Genet. 2014;15:176–92.View ArticlePubMedGoogle Scholar
- Wu CI. The genic view of the process of speciation. J Evol Biol. 2001;14:851–65.View ArticleGoogle Scholar
- Pinho C, Hey J. Divergence with gene flow: models and data. Annu Rev Ecol Evol S. 2010;41:215–30.View ArticleGoogle Scholar
- Hey J. Isolation with migration models for more than two populations. Mol Biol Evol. 2010;27:905–20.View ArticlePubMedGoogle Scholar
- Macholán M, Munclinger P, Sugerková M, Dufková P, Bímová B, Bozíková E, et al. Genetic analysis of autosomal and X-linked markers across a mouse hybrid zone. Evolution. 2007;61:746–71.View ArticlePubMedGoogle Scholar
- Storchová R, Reif J, Nachman MW. Female heterogamety and speciation: reduced introgression of the Z chromosome between two species of nightingales. Evolution. 2010;64:456–71.View ArticlePubMedGoogle Scholar
- Ellegren H, Smeds L, Burri R, Olason PI, Backström N, Kawakami T, et al. The genomic landscape of species divergence in Ficedula flycatchers. Nature. 2012;491:756–60.PubMedGoogle Scholar
- Martin SH, Dasmahapatra KK, Nadeau NJ, Salazar C, Walters JR, Simpson F, et al. Genome-wide evidence for speciation with gene flow in Heliconius butterflies. Genome Res. 2013;23:1817–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Nachman MW, Payseur BA. Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice. Phil Trans R Soc B. 2012;367:409–21.View ArticlePubMedPubMed CentralGoogle Scholar
- de León LF, Bermingham E, Podos J, Hendry AP. Divergence with gene flow as facilitated by ecological differences: within-island variation in Darwin’s finches. Philos Trans R Soc Lond B Biol Sci. 2010;365:1041–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Malinsky M, Challis RJ, Tyers AM, Schiffels S, Terai Y, Ngatunga BP, et al. Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake. Science. 2015;350:1493–8.View ArticlePubMedGoogle Scholar
- Keller I, Wagner CE, Greuter L, Mwaiko S, Selz OM, Sivasundar A, et al. Population genomic signatures of divergent adaptation, gene flow and hybrid speciation in the rapid radiation of Lake Victoria cichlid fishes. Mol Ecol. 2013;22:2848–63.View ArticlePubMedGoogle Scholar
- Sullivan J, Demboski JR, Bell KC, Hird S, Sarver B, Reid N, et al. Divergence with gene flow within the recent chipmunk radiation (Tamias). Heredity. 2014;113:185–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Lamichhaney S, Berglund J, Almén MS, Maqbool K, Grabherr M, Martinez-Barrio A, et al. Evolution of Darwin’s finches and their beaks revealed by genome sequencing. Nature. 2015;518:371–5.View ArticlePubMedGoogle Scholar
- Zhang W, Dasmahapatra KK, Mallet J, Moreira GR, Kronforst MR. Genome-wide introgression among distantly related Heliconius butterfly species. Genome Biol. 2016;17:25.View ArticlePubMedPubMed CentralGoogle Scholar
- Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491:444–8.View ArticlePubMedGoogle Scholar
- Leisler B, Schulze-Hagen K. The Reed Warblers: diversity in a uniform bird family. Zeist: KNNV Publishing; 2011.View ArticleGoogle Scholar
- Hansson B, Tarka M, Dawson DA, Horsburgh GJ. Hybridization but no evidence for backcrossing and introgression in a sympatric population of Great Reed Warblers and Clamorous Reed Warblers. PLoS One. 2012;7:e31667.View ArticlePubMedPubMed CentralGoogle Scholar
- Koskimies P. Polygyny in Blyth’s Reed Warbler Acrocephalus dumetorum. Ann Zool Fennici. 1984;21:239–42.Google Scholar
- Cramp S. Handbook of the birds of Europe, the Middle East and North Africa, vol. 6. New York: Oxford University Press; 1992.Google Scholar
- Lemaire F. Mixed song, interspecific competition and hybridisation in the reed and marsh warblers. Behaviour. 1977;63:215–40.View ArticleGoogle Scholar
- Otterbeck A, Dale S, Lindén A, Marthinsen G. A male Reed Warbler and Marsh Warbler hybrid (Acrocephalus scirpaceus × A. palustris) in Norway documented with molecular methods. Ornis Norvegica. 2013;36:6–13.View ArticleGoogle Scholar
- Fregin S, Haase M, Olsson U, Alström P. Multi-locus phylogeny of the family Acrocephalidae (Aves: Passeriformes) – the traditional taxonomy overthrown. Mol Phyl Evol. 2009;52:866–78.View ArticleGoogle Scholar
- Lindholm A, Bensch S, Dowsett-Lemaire F, Forsten A, Karkkainen H. Hybrid Marsh × Blyth’s Reed Warbler with mixed song in Finland in June 2003. Dutch Birding. 2007;29:223–31.Google Scholar
- Hagström B. Om en hybrid mellan rör - och busksångare. Vingspegeln. 1984;3:167.Google Scholar
- Griffiths R, Double MC, Orr K, Dawson RJG. A DNA test to sex most birds. Mol Ecol. 1998;7:1071–5.View ArticlePubMedGoogle Scholar
- Hall TA. BioEdit: a user friendly biological sequence alignment editor and analyses program for Windows 95/98/NT. Nucleic Acids Symp Series. 1999;41:95–8.Google Scholar
- Stephens M, Donnelly P. A comparison of Bayesian methods for haplotype reconstruction from population genetic data. Am J Hum Genet. 2003;73:1162–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–7.View ArticlePubMedGoogle Scholar
- Hudson RR, Kreitman M, Aguadé M. A test of neutral molecular evolution based on nucleotide data. Genetics. 1987;116:153–9.PubMedPubMed CentralGoogle Scholar
- Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.View ArticlePubMedGoogle Scholar
- Woerner AE, Cox MP, Hammer MF. Recombination-filtered genomic datasets by information maximization. Bioinformatics. 2007;23:1851–3.View ArticlePubMedGoogle Scholar
- Fregin S, Haase M, Olsson U, Alström P. Pitfalls in comparisons of genetic distances: a case study of the avian family Acrocephalidae. Mol Phyl Evol. 2012;62:319–28.View ArticleGoogle Scholar
- Weir JT, Schluter D. Calibrating the avian molecular clock. Mol Ecol. 2008;17:2321–8.View ArticlePubMedGoogle Scholar
- Hey J, Nielsen R. Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics. Proc Natl Acad Sci U S A. 2007;104:2785–90.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle Scholar
- Schluter D. Ecological causes of adaptive radiation. Am Nat. 1996;148:S40–64.View ArticleGoogle Scholar
- Berner D, Salzburger W. The genomics of organismal diversification illuminated by adaptive radiations. Trends Genet. 2015;31:491–9.View ArticlePubMedGoogle Scholar
- Currat M, Ruedi M, Petit RJ, Excoffier L. The hidden side of invasions: massive introgression by local genes. Evolution. 2008;62:1908–20.PubMedGoogle Scholar
- Alund M, Immler S, Rice AM, Qvarnström A. Low fertility of wild hybrid male flycatchers despite recent divergence. Biol Lett. 2013;9:20130169.View ArticlePubMedPubMed CentralGoogle Scholar
- Delmore KE, Hübner S, Kane NC, Schuster R, Andrew RL, Câmara F, et al. 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.View ArticlePubMedGoogle Scholar
- Backström N, Karaiskou N, Leder EH, Gustafsson L, Primmer CR, Qvarnström A, et al. A gene-based genetic linkage map of the collared flycatcher (Ficedula albicollis) reveals extensive synteny and gene-order conservation during 100 million years of avian evolution. Genetics. 2008;179:1479–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Muir G, Dixon CJ, Harper AL, Filatov DA. Dynamics of drift, gene flow, and selection during speciation in Silene. Evolution. 2012;66:1447–58.View ArticlePubMedGoogle Scholar
- Burri R, Nater A, Kawakami T, Mugal CF, Olason PI, Smeds L, et al. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Reifová R, Reif J, Antczak M, Nachman MW. Ecological character displacement in the face of gene flow: evidence from two species of nightingales. BMC Evol Biol. 2011;11:138.View ArticlePubMedPubMed CentralGoogle Scholar
- Saetre GP, Moum T, Bures S, Kral M, Adamjan M, Moreno J. A sexually selected character displacement in flycatchers reinforces premating isolation. Nature. 1997;387:589–92.View ArticleGoogle Scholar