Flycatcher song in allopatry and sympatry--convergence, divergence and reinforcement

The theory of reinforcement predicts that natural selection against the production of unfit hybrids favours traits that increase assortative mating. Whether culturally inherited traits, such as bird song, can increase assortative mating by reinforcement is largely unknown. We compared songs of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) from two hybrid zones of different ages with songs from allopatric populations. Previously, a character divergence in male plumage traits has been shown to reinforce premating isolation in sympatric flycatchers. In contrast, we find that the song of the pied flycatcher has converged towards that of the collared flycatcher (mixed singing). However, a corresponding divergence in the collared flycatcher shows that the species differences in song characters are maintained in sympatry. Genetic analyses suggest that mixed song is not caused by introgression from the collared flycatcher, but rather due to heterospecific copying. Circumstantial evidence suggests that mixed song may increase the rate of maladaptive hybridization. In the oldest hybrid zone where reinforcement on plumage traits is most pronounced, the frequency of mixed singing and hybridization is also lowest. Thus, we suggest that reinforcement has reduced the frequency of mixed singing in the pied flycatcher and caused a divergence in the song of the collared flycatcher. Whether a culturally inherited trait promotes or opposes speciation in sympatry may depend on its plasticity. The degree of plasticity may be genetically determined and accordingly under selection by reinforcement.     

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (P_ST) was higher than that in neutral genetic markers (F_ST), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non‐independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.     

PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (N_e) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in N_e over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. N_e curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an N_e of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate N_e trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in N_e differed among species, and for pied flycatcher, the temporal dynamics of N_e differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.     

Malaria infections reinforce competitive asymmetry between two Ficedula flycatchers in a recent contact zone

Parasites may influence the outcome of interspecific competition between closely related host species through lower parasite virulence in the host with which they share the longer evolutionary history. We tested this idea by comparing the prevalence of avian malaria (Haemosporidia) lineages and their association with survival in pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breeding in a recent contact zone on the Swedish island of Öland. A nested PCR protocol amplifying haemosporidian fragments of mtDNA was used to screen the presence of malaria lineages in 1048 blood samples collected during 6 years. Competitively inferior pied flycatchers had a higher prevalence of blood parasites, including the lineages that were shared between the two flycatcher species. Multistate mark–recapture models revealed a lower survival of infected versus uninfected female pied flycatchers, while no such effects were detected in male pied flycatchers or in collared flycatchers of either sex. Our results show that a comparatively new host, the collared flycatcher, appears to be less susceptible to a local northern European malarial lineage where the collared flycatchers have recently expanded their distribution. Pied flycatchers experience strong reproductive interference from collared flycatchers, and the additional impact of species‐specific blood parasite effects adds to this competitive exclusion. These results support the idea that parasites can strongly influence the outcome of interspecific competition between closely related host species, but that the invading species need not necessarily be more susceptible to local parasites.     

Contrasting patterns of polymorphism and divergence on the Z chromosome and autosomes in two Ficedula flycatcher species

In geographic areas where pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breed in sympatry, hybridization occurs, leading to gene flow (introgression) between the two recently diverged species. Notably, while such introgression is observable at autosomal loci it is apparently absent at the Z chromosome, suggesting an important role for genes on the Z chromosome in creating reproductive isolation during speciation. To further understand the role of Z-linked loci in the formation of new species, we studied genetic variation of the two species from regions where they live in allopatry. We analyzed patterns of polymorphism and divergence in introns from 9 Z-linked and 23 autosomal genes in pied and collared flycatcher males. Average variation on the Z chromosome is greatly reduced compared to neutral expectations based on autosomal diversity in both species. We also observe significant heterogeneity between patterns of polymorphism and divergence at Z-linked loci and a relative absence of polymorphisms that are shared by the two species on the Z chromosome compared to the autosomes. We suggest that these observations may indicate the action of recurrent selective sweeps on the Z chromosome during the evolution of the two species, which may be caused by sexual selection acting on Z-linked genes. Alternatively, reduced variation on the Z chromosome could result from substantially higher levels of introgression at autosomal than at Z-linked loci or from a complex demographic history, such as a population bottleneck.     

Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome

Although there is growing interest to take genomics into the complex realms of natural populations, there is a general shortage of genomic resources and tools available for wild species. This applies not at least to birds, for which genomic approaches should be helpful to questions such as adaptation, speciation and population genetics. In this study, we describe a genome-wide reference set of conserved avian gene markers, broadly applicable across birds. By aligning protein-coding sequences from the recently assembled chicken genome with orthologous sequences in zebra finch, we identified particularly conserved exonic regions flanking introns of suitable size for subsequent amplification and sequencing. Primers were designed for 242 gene markers evenly distributed across the chicken genome, with a mean inter-marker interval of 4.2 Mb. Between 78% and 93% of the markers amplified a specific product in five species tested (chicken, peregrine falcon, collared flycatcher, great reed warbler and blue tit). Two hundred markers were sequenced in collared flycatcher, yielding a total of 122.41 kb of genomic DNA sequence (12096 bp coding sequence and 110 314 bp noncoding). Intron size of collared flycatcher and chicken was highly correlated, as was GC content. A polymorphism screening using these markers in a panel of 10 unrelated collared flycatchers identified 871 single nucleotide polymorphisms (pi = 0.0029) and 33 indels (mainly very short). Avian genome characteristics such as uniform genome size and low rate of syntenic rearrangements suggest that this marker set will find broad utility as a genome-wide reference resource for molecular ecological and population genomic analysis of birds. We envision that it will be particularly useful for obtaining large-scale orthologous targets in different species--important in, for instance, phylogenetics--and for large-scale identification of evenly distributed single nucleotide polymorphisms needed in linkage mapping or in studies of gene flow and hybridization.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

Song similarity predicts hybridization in flycatchers

Given that population divergence in sexual signals is an important prerequisite for reproductive isolation, a key prediction is that cases of signal convergence should lead to hybridization. However, empirical studies that quantitatively demonstrate links between phenotypic characters of individuals and their likelihood to hybridize are rare. Here we show that song convergence between sympatric pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) influence social and sexual interactions between the two species. In sympatry, the majority of male pied flycatchers (65%) include various parts of collared flycatcher song in their song repertoire (but not vice versa). Playback experiments on male interactions demonstrate that male collared flycatchers respond similarly to this 'mixed' song as to conspecific song. Long-term data on pairing patterns show that males singing a converged song attract females of the other species: female collared flycatchers only pair with male pied flycatchers if the males sing the mixed song type. From the perspective of a male pied flycatcher, singing a mixed song type is associated with 30% likelihood of hybridization. This result, combined with our estimates of the frequency of mixed singers, accurately predicts the observed occurrence of hybridization among male pied flycatchers in our study populations (20.45% of 484 pairs; predicted 19.5%). Our results support the suggestion that song functions as the most important prezygotic isolation mechanism in many birds.     

Genetic evidence on the demography of speciation in allopatric dolphin species

Under a neutral model, the stochastic lineage sorting that leads to gene monophyly proceeds slowly in large populations. Therefore, in many recent species with large population size, the genome will have mixed support for monophyly unless historical bottlenecks have accelerated coalescence. We use genealogical patterns in mitochondrial DNA and in introns of four nuclear loci to test for historical bottlenecks during the speciation and divergence of two temperate Lagenorhynchus dolphin species isolated by tropical Pacific waters (an antitropical distribution). Despite distinct morphologies, foraging behaviors, and mitochondrial DNAs, these dolphin species are polyphyletic at all four nuclear loci. The abundance of shared polymorphisms between these sister taxa is most consistent with the maintenance of large effective population sizes (5.09 x 10(4) to 10.9 x 10(4)) during 0.74-1.05 million years of divergence. A variety of population size histories are possible, however. We used gene tree coalescent probabilities to explore the rejection region for historical bottlenecks of different intensity given best estimates of effective population size under a strict isolation model of divergence. In L. obliquidens the data are incompatible with a colonization propagule of an effective size of 10 or fewer individuals. Although the ability to reject less extreme historical bottlenecks will require data from additional loci, the intermixed genealogical patterns observed between these dolphin sister species are highly probable only under an extended history of large population size. If similar demographic histories are inferred for other marine antitropical taxa, a parsimonious model for the Pleistocene origin of these distributions would not involve rare breaches of a constant dispersal barrier by small colonization propagules. Instead, a history of large population size in L. obliquidens and L. obscurus contributes to growing biological and environmental evidence that the equatorial barrier became permeable during glacial/interglacial cycles, leading to vicariant isolation of antitropical populations.     

Positive feedback between ecological and reproductive character displacement in a young avian hybrid zone

Character displacement can reduce costly interspecific interactions between young species. We investigated the mechanisms behind divergence in three key traits-breeding habitat choice, timing of breeding, and plumage coloration-in Ficedula flycatchers. We found that male pied flycatchers became expelled from the preferred deciduous habitat into mixed forest as the superior competitor, collared flycatchers, increased in numbers. The peak in food abundance differs between habitats, and the spatial segregation was paralleled by an increased divergence in timing of breeding between the two species. Male pied flycatchers vary from brown to black with brown coloration being more frequent in sympatry with collared flycatchers, a pattern often proposed to result from selection against hybridization, that is, reinforcement. In contrast to this view, we show that brown male pied flycatchers more often hybridize than black males. Male pied flycatcher plumage coloration influenced the territory obtained in areas of co-occurrence with collared flycatchers, and brown male pied flycatchers experienced higher relative fitness than black males when faced with heterospecific competition. We suggest that allopatric divergence in resource defense ability causes a feedback loop at secondary contact where male pied flycatchers with the most divergent strategy compared to collared flycatchers are favored by selection.     

Single-nucleotide polymorphism characterization in species with limited available sequence information: high nucleotide diversity revealed in the avian genome

As a case study for single-nucleotide polymorphism (SNP) identification in species for which little or no sequence information is available, we investigated several approaches to identifying SNPs in two passerine bird species: pied and collared flycatchers (Ficedula hypoleuca and F. albicollis). All approaches were successful in identifying sequence polymorphism and over 50 candidate SNPs per species were identified from approximately 9.1 kb of sequence. In addition, 17 sites were identified in which the frequency of alternative bases differed by > 50% between species (termed interspecific SNPs). Interestingly, polymorphism of microsatellite/intron loci in the source species appeared to be a positive predictor of nucleotide diversity in homologous flycatcher sequences. The overall nucleotide diversity of flycatchers was 2.3-2.7 x 10(-3), which is approximately 3-6 times higher than observed in recent studies of human SNPs. Higher nucleotide diversity in the avian genome could be due to the relatively older age of flycatcher populations, compared with humans, and/or a higher long-term effective population size.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

No evidence for Z‐chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene‐based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post‐zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z‐chromosome. Could rearrangements of the Z‐chromosome between the species explain their reproductive isolation? We developed high coverage Z‐chromosome linkage maps for both species, using gene‐based markers and large‐scale SNP genotyping. Best order maps contained 57–62 gene markers with an estimated average density of one every 1–1.5 Mb. We estimated the recombination rates in flycatcher Z‐chromosomes to 1.1–1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co‐linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post‐zygotic reproductive isolation between these two Ficedula species.     

Difference in plasticity of resting metabolic rate – the proximate explanation to different niche breadth in sympatric Ficedula flycatchers

Variation in relative fitness of competing recently formed species across heterogeneous environments promotes coexistence. However, the physiological traits mediating such variation in relative fitness have rarely been identified. Resting metabolic rate (RMR) is tightly associated with life history strategies, thermoregulation, diet use, and inhabited latitude and could therefore moderate differences in fitness responses to fluctuations in local environments, particularly when species have adapted to different climates in allopatry. We work in a long‐term study of collared (Ficedula albicollis) and pied flycatchers (Ficedula hypoleuca) in a recent hybrid zone located on the Swedish island of Öland in the Baltic Sea. Here, we explore whether differences in RMR match changes in relative performance of growing flycatcher nestlings across environmental conditions using an experimental approach. The fitness of pied flycatchers has previously been shown to be less sensitive to the mismatch between the peak in food abundance and nestling growth among late breeders. Here, we find that pied flycatcher nestlings have lower RMR in response to higher ambient temperatures (associated with low food availability). We also find that experimentally relaxed nestling competition is associated with an increased RMR in this species. In contrast, collared flycatcher nestlings did not vary their RMR in response to these environmental factors. Our results suggest that a more flexible nestling RMR in pied flycatchers is responsible for the better adaptation of pied flycatchers to the typical seasonal changes in food availability experienced in this hybrid zone. Generally, subtle physiological differences that have evolved when species were in allopatry may play an important role to patterns of competition, coexistence, or displacements between closely related species in secondary contact.     

Past and potential future population dynamics of three grouse species using ecological and whole genome coalescent modeling

Studying demographic history of species provides insight into how the past has shaped the current levels of overall biodiversity and genetic composition of species, but also how these species may react to future perturbations. Here we investigated the demographic history of the willow grouse (Lagopus lagopus), rock ptarmigan (Lagopus muta), and black grouse (Tetrao tetrix) through the Late Pleistocene using two complementary methods and whole genome data. Species distribution modeling (SDM) allowed us to estimate the total range size during the Last Interglacial (LIG) and Last Glacial Maximum (LGM) as well as to indicate potential population subdivisions. Pairwise Sequentially Markovian Coalescent (PSMC) allowed us to assess fluctuations in effective population size across the same period. Additionally, we used SDM to forecast the effect of future climate change on the three species over the next 50 years. We found that SDM predicts the largest range size for the cold‐adapted willow grouse and rock ptarmigan during the LGM. PSMC captured intraspecific population dynamics within the last glacial period, such that the willow grouse and rock ptarmigan showed multiple bottlenecks signifying recolonization events following the termination of the LGM. We also see signals of population subdivision during the last glacial period in the black grouse, but more data are needed to strengthen this hypothesis. All three species are likely to experience range contractions under future warming, with the strongest effect on willow grouse and rock ptarmigan due to their limited potential for northward expansion. Overall, by combining these two modeling approaches, we have provided a multifaceted examination of the biogeography of these species and how they have responded to climate change in the past. These results help us understand how cold‐adapted species may respond to future climate changes.     

Heterospecific female mimicry in Ficedula flycatchers

Mimicry is a widespread phenomenon. Vertebrate visual mimicry often operates in an intraspecific sexual context, with some males resembling conspecific females. Pied flycatcher (Ficedula hypoleuca) dorsal plumage varies from the ancestral black to female‐like brown. Experimental studies have shown that conspecific and heterospecific (collared flycatcher, F. albicollis) individuals of both sexes respond, at least initially, to brown individuals as if they were female. We quantified the perceptual and biochemical differences between brown feathers and found that brown pied flycatcher males are indistinguishable from heterospecific, but not from conspecific, females in both aspects. To our knowledge, this is the first evidence of a visual mimetic signalling system in a sexual context where the model is heterospecific to the mimic. By only mimicking heterospecific females, brown pied flycatcher males can establish territories next to the more dominant collared flycatcher in sympatry, suffer less aggression by darker conspecifics in allopatry and preserve within‐species sexual recognition throughout the breeding range. A closer look at the evolutionary history and ecology of these two species illustrates how such a mimetic system can evolve. Although likely rare, this phenomenon might not be unique to Ficedula flycatchers.     

HALDANE'S RULE AND SEX BIASSED GENE FLOW BETWEEN TWO HYBRIDIZING FLYCATCHER SPECIES (FICEDULA ALBICOLLIS AND F. HYPOLEUCA,AVES: MUSCICAPIDAE)

The collared flycatcher (Ficedula albicollis) and the pied flycatcher (F. hypoleuca) hybridize where their geographic ranges overlap. Restriction fragment comparison of 5% of the mitochondrial genome showed a sequence divergence of 10% between these flycatcher species. This degree of sequence divergence between a closely related pair of bird species is unusually high and contrasts with the low level of divergence between F. albicollis and F. hypoleuca in nuclear genes (Nei's D = 0.0006) revealed by enzyme electrophoresis. The low nuclear differentiation is explained by sex biassed gene flow and introgression in nuclear genes (via fertile male hybrids), while the high mitochondrial DNA sequence divergence is preserved by sterility of female hybrids, which prevents mitochondrial introgression. This pattern is in accordance with Haldane's rule and is supported by field data on hybrid fertility. The high mtDNA differentiation could be explained by transfer of mitochondrial DNA from a third species during a past period of hybridization.     

Demographic history of a common pioneer tree,Zanthoxylum ailanthoides,reconstructed using isolation-with-migration model

The Japanese prickly ash, Zanthoxylum ailanthoides Siebold & Zucc. (Rutaceae), is one of the most common pioneer tree species that grow in disturbed areas, such as canopy gaps, in Japanese warm-temperate evergreen oak forests. The strong genetic structure of its current population might suggest long-term isolation of subpopulations by demographic events, in addition to ecological features of this species. Analyses of genome-wide nucleotide variation using model-based approaches are necessary to achieve a good understanding of the demographic history of a species. In this study, we analyzed the nucleotide variation in natural populations of Z. ailanthoides using a computer program that applied the isolation-with-migration model to quantitatively infer the demographic history. Nucleotide variations at 10 or 26 nuclear loci in six populations from a wide range of the species distribution were analyzed. The maximum likelihood estimate of the divergence time between the current populations in the two hypothetical refugia of Z. ailanthoides was approximately 24,000 years. Unexpectedly, the estimated size of the ancestral population was larger than the sizes of the two current populations. The results suggested a relatively recent divergence of these two populations and rapid formation of strong genetic structure among subpopulations. The large ancestral population may indicate a more complex demographic history during or after the last glacial period than the simple isolation-with-migration model implies.     

Amount of introgression in flycatcher hybrid zones reflects regional differences in pre and post-zygotic barriers to gene exchange

Introgression is the incorporation of alleles from one species or semispecies into the gene pool of another through hybridization and backcrossing. The rate at which this occurs depends on the frequency of hybridization and the fitness of hybrids and backcrosses compared to 'pure' individuals. The collared flycatcher (Ficedula albicollis) and the pied flycatcher (F. hypoleuca) co-exist and hybridize at low to moderate frequencies in a clinal hybrid zone in Central Europe and on the islands of Gotland and Oland off the Swedish east coast. Data on hatching success suggest that hybrids are less fertile in Central Europe compared to on the islands. Direct fitness estimates using molecular markers to infer paternity are consistent with the demographic data. Applying a tag-array-based minisequencing assay to genotype interspecific substitutions and single nucleotide polymorphisms we demonstrate that the amount of introgression from the pied to the collared flycatcher is higher in the two island populations (Gotland and Oland) than in two geographically distinct areas from the Central European hybrid zone (Czech Republic and Hungary). In all areas the amount of introgression from collared to pied flycatchers is very low or seemingly absent. The different patterns of introgression are consistent with regional differences in rates of hybridization and fitness of hybrids. We suggest that barriers to gene exchange may have been partly broken down on the islands due to asymmetric gene flow from allopatry. Alternatively, or in addition, more pronounced reinforcement of prezygotic isolation in Central Europe might have increased post-zygotic isolation through hitchhiking, since genes affecting pre and post-zygotic isolation are both sex-linked in these birds. One of our genetic markers appears to introgress from pied to collared flycatchers at a much higher rate than the other markers. We discuss the possibility that the introgressed marker may be linked to a gene which is under positive selection in the novel genetic background.