A high‐density linkage map enables a second‐generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high‐density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best‐order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM . Facilitated by the array being designed to include markers from most scaffolds, we obtained a second‐generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super‐scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM /Mb and correlated closely with chromosome size, from 2 cM /Mb for chromosomes >100 Mb to >10 cM /Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.     

Habitat preference and flowering‐time variation contribute to reproductive isolation between diploid and autotetraploid Anacamptis pyramidalis

Tetraploid lineages are typically reproductively isolated from their diploid ancestors by post‐zygotic isolation via triploid sterility. Nevertheless, polyploids often also exhibit ecological divergence that could contribute to reproductive isolation from diploid ancestors. In this study, we disentangled the contribution of different forms of reproductive isolation between sympatric diploid and autotetraploid individuals of the food‐deceptive orchid Anacamptis pyramidalis by quantifying the strength of seven reproductive barriers: three prepollination, one post‐pollination prezygotic and three post‐zygotic. The overall reproductive isolation between the two cytotypes was found very high, with a preponderant contribution of two prepollination barriers, that is phenological and microhabitat differences. Although the contribution of post‐zygotic isolation (triploid sterility) is confirmed in our study, these results highlight that prepollination isolation, not necessarily involving pollinator preference, can represent a strong component of reproductive isolation between different cytotypes. Thus, in the context of polyploidy as quantum speciation, that generates reproductive isolation via triploid sterility, ecological divergence can strengthen the reproductive isolation between cytotypes, reducing the waste of gametes in low fitness interploidy crosses and thus favouring the initial establishment of the polyploid lineage. Under this light, speciation by polyploidy involves ecological processes and should not be strictly considered as a nonecological form of speciation.     

The effectiveness of pre‐ and post‐zygotic barriers in avoiding hybridization between two snapdragons (Antirrhinum L.: Plantaginaceae)

Reproductive barriers play an important role in the maintenance of species boundaries. However, to date, few studies have provided a detailed analysis of reproductive isolation barriers between species or examined their importance in maintaining species identity. This is the first detailed study into pre‐ and post‐zygotic reproductive isolation barriers in Antirrhinum, based on a mixed population with two species that rarely co‐occur. The study revealed that pollinator constancy and preference and poor hybrid seed viability were the most important reproductive isolating mechanisms. Reproductive isolation was practically complete by both pre‐ and post‐zygotic barriers. Average pre‐zygotic isolation was greater than post‐zygotic isolation, in accordance with the trend observed in flowering plants in which reproductive isolation is principally caused by pre‐zygotic mechanisms. However, average post‐zygotic isolation was also high, in contrast to what was expected among Antirrhinum spp. This case highlights the importance of quantifying the reproductive isolation barriers thoroughly to understand how and why species boundaries are maintained. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 176 , 159–172.     

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.     

Multiple sources of reproductive isolation in a bimodal butterfly hybrid zone

An important evolutionary question concerns whether one or many barriers are involved in the early stages of speciation. We examine pre‐ and post‐zygotic reproductive barriers between two species of butterflies (Heliconius erato chestertonii and H. e. venus) separated by a bimodal hybrid zone in the Cauca Valley, Colombia. We show that there is both strong pre‐ and post‐mating reproductive isolation, together leading to a 98% reduction in gene flow between the species. Pre‐mating isolation plays a primary role, contributing strongly to this isolation (87%), similar to previous examples in Heliconius. Post‐mating isolation was also strong, with absence of Haldane’s rule, but an asymmetric reduction in fertility (< 11%) in inter‐specific crosses depending on maternal genotype. In summary, this is one of the first examples of post‐zygotic reproductive isolation playing a significant role in early stages of parapatric speciation in Heliconius and demonstrates the importance of multiple barriers to gene flow in the speciation process.     

Evolution of pre‐zygotic and post‐zygotic barriers to gene flow among three cryptic species within the Anastrepha fraterculus complex

Tropical tephritids are ideally suited for studies on population divergence and speciation because they include species groups undergoing rapid radiation, in which morphologically cryptic species and sister species are abundant. The fraterculus species group in the Neotropical genus Anastrepha is a case in point, as it is composed of a complex of up to seven A. fraterculus morphotypes proposed to be cryptic species. Here, we document pre‐ and post‐zygotic barriers to gene flow among adults of the Mexican A. fraterculus morphotype and three populations (Argentina, Brazil, and Peru) belonging to two separate morphotypes (Brazilian 1 and Peruvian). We unveiled three forms of pre‐zygotic reproductive isolation resulting in strong assortative mating. In field cages, free‐ranging male and female A. fraterculus displayed a strong tendency to form couples with members of the opposite sex belonging to their own morphotype, suggesting that male pheromone emission, courtship displays, or both intervene in shaping female choice before actual contact and coupling. In addition, males and females of the Peruvian morphotype became receptive and mated significantly later than adults of the Mexican and Brazilian 1 morphotypes. After contact, Mexican females exhibited greater mating discrimination than males when facing adults of the opposite sex belonging to either the Peruvian or the Brazilian 1 morphotype as evidenced by vigorous resistance to penetration once they had been forcefully mounted by heterotypic males. Forced copulations resulted in production of F1 hybrids that were either less viable (and partially fertile) than parental crosses or even sterile. Our results suggest that the Mexican morphotype is a distinct biological entity and that pre‐zygotic reproductive isolation through divergence in courtship or male‐produced pheromone and other mechanisms appear to evolve faster than post‐zygotic isolation in the fraterculus species group.     

Abundant recent activity of retrovirus‐like retrotransposons within and among flycatcher species implies a rich source of structural variation in songbird genomes

Transposable elements (TEs) are genomic parasites capable of inserting virtually anywhere in the host genome, with manifold consequences for gene expression, DNA methylation and genomic stability. Notably, they can contribute to phenotypic variation and hence be associated with, for example, local adaptation and speciation. However, some organisms such as birds have been widely noted for the low densities of TEs in their genomes and this has been attributed to a potential dearth in transposition during their evolution. Here, we show that avian evolution witnessed diverse and abundant transposition on very recent timescales. First, we made an in‐depth repeat annotation of the collared flycatcher genome, including identification of 23 new, retrovirus‐like LTR retrotransposon families. Then, using whole‐genome resequencing data from 200 Ficedula flycatchers, we detected 11,888 polymorphic TE insertions (TE presence/absence variations, TEVs) that segregated within and among species. The density of TEVs was one every 1.5–2.5 Mb per individual, with heterozygosities of 0.12–0.16. The majority of TEVs belonged to some 10 different LTR families, most of which are specific to the flycatcher lineage. TEVs were validated by tracing the segregation of hundreds of TEVs across a three‐generation pedigree of collared flycatchers and also by their utility as markers recapitulating the phylogenetic relationships among flycatcher species. Our results suggest frequent germline invasions of songbird genomes by novel retroviruses as a rich source of structural variation, which may have had underappreciated phenotypic consequences for the diversification of this species‐rich group of birds.     

Ongoing speciation within the Anastrepha fraterculus cryptic species complex: the case of the Andean morphotype

The Anastrepha fraterculus (Wiedemann) (Diptera: Tephritidae) cryptic species complex is currently composed of seven taxonomically recognized morphotypes. Both, pre‐ and post‐zygotic isolation has been documented among four of these morphotypes, revealing that in fact they appear to be distinct biological entities. In order to progress in the full delimitation of species within the complex, we examined reproductive isolation between a Colombian population of the Andean morphotype and populations belonging to four other morphotypes spanning from Mexico to Argentina. Flies from the Andean morphotype exhibited strong pre‐zygotic mating isolation through temporal partitioning of mating activity. Post‐zygotic isolation was observed for crosses of males of all morphotypes and Andean morphotype females, yet most of the F1 hybrid ♂ × F1 hybrid ♀ self‐crosses showed normal levels of fertility, a finding suggesting a nuclear–cytoplasmic interaction according to previous studies. Overall, the Andean morphotype within the complex also appears to be a distinct biological entity. We discuss the implications of these findings for the understanding of speciation mechanisms in the Neotropical genus Anastrepha.     

Reproductive barriers between two sympatric beetle species specialized on different host plants

Knowledge on interspecific pre‐ and post‐zygotic isolation mechanisms provides insights into speciation patterns. Using crosses (F₁ and backcrosses) of two closely related flea beetles species, Altica fragariae and A. viridicyanea, specialized on different hosts in sympatry, we measured: (a) the type of reproductive isolation and (b) the inheritance mode of preference and host‐specific performance, using a joint‐scaling test. Each species preferred almost exclusively its host plant, creating strong prezygotic isolation between them, and suggesting that speciation may occur at least partly in sympatry. Reproductive isolation was intrinsic between females of A. fragariae and either A. viridicyanea or F₁ males, whereas the other crosses showed ecologically dependent reproductive isolation, suggesting ecological speciation. The genetic basis of preference and performance was at least partially independent, and several loci coded for preference, which limits the possibility of sympatric speciation. Hence, both ecological and intrinsic factors may contribute to speciation between these species.     

Incipient post‐zygotic barrier in a model system of ecological speciation with gene flow

The role of post‐zygotic isolation in nonallopatric ecological speciation is still mostly unknown and information on the nature and strength of these barriers in well‐known speciation models is essential for a deeper understanding of such processes. The Galician ecotypes of the marine snail Littorina saxatilis represent one of the best studied cases of nonallopatric ecological speciation. Here, we test the existence of incipient post‐zygotic isolation by comparing the fertility of male hybrids with that of both pure forms [ridged and banded (RB) and smooth and unbanded (SU) ecotypes]. We analysed the degree of sperm DNA fragmentation (SDF) of individuals morphologically classified as RB, SU and hybrids, sampled from two locations. SDF analyses were chosen to study sperm quality because, in other animal species, SDF rates correlate with important parameters for speciation research, such as fertilization and abortion rates and viability of adult progeny. In the present work, hybrids showed significantly higher SDF rates than RB and SU males in one location and significantly higher variances in both locations. These results suggest the existence of an incipient post‐zygotic barrier, the strength of which may vary across the Galician shore, and highlight the potential of SDF analyses for speciation research.     

Discordant divergence times among Z-chromosome regions between two ecologically distinct swallowtail butterfly species

We investigate multilocus patterns of differentiation between parental populations of two swallowtail butterfly species that differ at a number of ecologically important sex-linked traits. Using a new coalescent-based approach, we show that there is significant heterogeneity in estimated divergence times among five Z-linked markers, rejecting a purely allopatric speciation model. We infer that the Z chromosome is a mosaic of regions that differ in the extent of historical gene flow, potentially due to isolating barriers that prevent the introgression of species-specific traits that result in hybrid incompatibilities. Surprisingly, a candidate region for a strong barrier to introgression, Ldh, does not show a significantly deeper divergence time than other markers on the Z chromosome. Our approach can be used to test alternative models of speciation and can potentially assign chronological order to the appearance of factors contributing to reproductive isolation between species.     

Asymmetric Assortative Mating Behaviour Reflects Incomplete Pre‐zygotic Isolation in the Nasonia Species Complex

Preference of con‐ over heterospecific mates leading to assortative mating can substantially contribute to pre‐zygotic reproductive isolation and prevent fitness losses if post‐zygotic hybridization barriers already exist. The jewel wasp genus Nasonia displays quite strong and well‐studied post‐zygotic reproductive isolation due to a ubiquitous Wolbachia infection causing cytoplasmic incompatibility between different species. Pre‐zygotic isolation, however, has received far less research attention in this model organism, especially concerning the mechanisms and criteria of mate choice. In the present study, we analysed mate rejection and mate acceptance rates in cross‐comparisons between all four Nasonia species. We put emphasis on observing which sex is more likely to interrupt interspecific matings and how discriminatory behaviour varies across the different species in all possible combinations. We found an asymmetric distribution of assortative mating among the four Nasonia species that appears to be highly influenced by the respective combinations of sex and species. Females appeared to be the main discriminators against heterospecific mating partners, but interestingly, we could also detect mate discrimination and rejection behaviour in males, a widely neglected factor in research on mating behaviour in general and on Nasonia in particular. Moreover, the asymmetry in the assortative mating behaviour was partially reflective of sym‐ or allopatric distributions of natural Nasonia populations.     

Increased divergence but reduced variation on the Z chromosome relative to autosomes in Ficedula flycatchers: differential introgression or the faster-Z effect?

Recent multilocus studies of congeneric birds have shown a pattern of elevated interspecific divergence on the Z chromosome compared to the autosomes. In contrast, intraspecifically, birds exhibit less polymorphism on the Z chromosome relative to the autosomes. We show that the four black-and-white Ficedula flycatcher species show greater genetic divergence on the Z chromosome than on the autosomes, and that the ratios of intraspecific polymorphism at Z-linked versus autosomal markers are below the neutral expectation of 75%. In all species pairs, we found more fixed substitutions and fewer shared polymorphisms on the Z chromosome than on the autosomes. Finally, using isolation with migration (IMa) models we estimated gene flow among the four closely related flycatcher species. The results suggest that different pattern of evolution of Z chromosomes and autosomes is best explained by the faster-Z hypothesis, since the estimated long-term gene flow parameters were close to zero in all comparisons.     

Chromosomal variation segregates within incipient species and correlates with reproductive isolation

Reproductive isolation is a critical step in the process of speciation. Among the most important factors driving reproductive isolation are genetic incompatibilities. Whether these incompatibilities are already present before extrinsic factors prevent gene flow between incipient species remains largely unresolved in natural systems. This question is particularly challenging because it requires that we catch speciating populations in the act before they reach the full‐fledged species status. We measured the extent of intrinsic postzygotic isolation within and between phenotypically and genetically divergent lineages of the wild yeast Saccharomyces paradoxus that have partially overlapping geographical distributions. We find that hybrid viability between lineages progressively decreases with genetic divergence. A large proportion of postzygotic inviability within lineages is associated with chromosomal rearrangements, suggesting that chromosomal differences substantially contribute to the early steps of reproductive isolation within lineages before reaching fixation. Our observations show that polymorphic intrinsic factors may segregate within incipient species before they contribute to their full reproductive isolation and highlight the role of chromosomal rearrangements in speciation. We propose different hypotheses based on adaptation, biogeographical events and life history evolution that could explain these observations.     

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.     

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate ‘genomic island of speciation’

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a ‘genomic island of divergence’ spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X‐island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X‐island despite a lack of autosomal divergence; (ii) individuals with multiple‐recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination‐driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.     

Genomic landscape of reproductive isolation in Lucania killifish: The role of sex loci and salinity

Adaptation to different environments can directly and indirectly generate reproductive isolation between species. Bluefin killifish (Lucania goodei) and rainwater killifish (L. parva) are sister species that have diverged across a salinity gradient and are reproductively isolated by habitat, behavioural, extrinsic and intrinsic post‐zygotic isolation. We asked if salinity adaptation contributes indirectly to other forms of reproductive isolation via linked selection and hypothesized that low recombination regions, such as sex chromosomes or chromosomal rearrangements, might facilitate this process. We conducted QTL mapping in backcrosses between L. parva and L. goodei to explore the genetic architecture of salinity tolerance, behavioural isolation and intrinsic isolation. We mapped traits relative to a chromosome that has undergone a centric fusion in L. parva (relative to L. goodei). We found that the sex locus appears to be male determining (XX‐XY), was located on the fused chromosome and was implicated in intrinsic isolation. QTL associated with salinity tolerance were spread across the genome and did not overly co‐localize with regions associated with behavioural or intrinsic isolation. This preliminary analysis of the genetic architecture of reproductive isolation between Lucania species does not support the hypothesis that divergent natural selection for salinity tolerance led to behavioural and intrinsic isolation as a by‐product. Combined with previous studies in this system, our work suggests that adaptation as a function of salinity contributes to habitat isolation and that reinforcement may have contributed to the evolution of behavioural isolation instead, possibly facilitated by linkage between behavioural isolation and intrinsic isolation loci on the fused chromosome.     

A molecular study of hybridization and homoploid hybrid speciation in Argyranthemum (Asteraceae) on Tenerife,the Canary Islands

Examples of recurrent homoploid hybrid speciation are few. One often‐cited example is Argyranthemum sundingii. This example includes two described species, A. lemsii and A. sundingii, resulting from reciprocal hybridization between A. broussonetii and A. frutescens on Tenerife. The four species and artificial F₁ and F₂ hybrids have previously been investigated morphologically and cytologically. Here, we examine population differentiation based on amplified fragment length polymorphism to get a better understanding of the genetic relationships among the species and the extent of hybridization. We aim to investigate if there is molecular support for treating the hybrid species as one taxon. Seven parental and four hybrid species populations (149 individuals) were analysed and we scored 85 polymorphic markers. A few (2–5) were private to each species but variably present and mostly rare. Our principal coordinate, STRUCTURE and BAPS analyses and AMOVA resulted in a clear separation of the parental species. The hybrid species were genetically less divergent but not identical. Our data indicate that hybridization and introgression are common in all these species on Tenerife and support the hypothesis that homoploid hybrid speciation has occurred repeatedly. Intrinsic post‐zygotic barriers are notoriously weak in Argyranthemum and reproductive isolation and speciation result primarily from strong ecological selection. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 159 , 19–31.     

GENOME‐WIDE INVESTIGATION OF REPRODUCTIVE ISOLATION IN EXPERIMENTAL LINEAGES AND NATURAL SPECIES OF NEUROSPORA: IDENTIFYING CANDIDATE REGIONS BY MICROARRAY‐BASED GENOTYPING AND MAPPING

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe‐dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky–Muller model of genic incompatibility. We developed microarray‐based, restriction‐site associated DNA (RAD) markers that identified ～1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk‐segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only ～2% of the nuclear genome. This study demonstrates how advances in high‐throughput, genome‐wide genotyping can be applied to mapping reproductive isolation genes and speciation research.     

Possible glimpses into early speciation: the effect of ovarian fluid on sperm velocity accords with post‐copulatory isolation between two guppy populations

Identifying mechanisms of reproductive isolation is key to understanding speciation. Among the putative mechanisms underlying reproductive isolation, sperm–female interactions (post‐mating–prezygotic barriers) are arguably the hardest to identify, not least because these are likely to operate at the cellular or molecular level. Yet sperm–female interactions offer great potential to prevent the transfer of genetic information between different populations at the initial stages of speciation. Here, we provide a preliminary test for the presence of a putative post‐mating–prezygotic barrier operating between three populations of Trinidadian guppies (Poecilia reticulata), an internally fertilizing fish that inhabits streams with different levels of connectivity across Trinidad. We experimentally evaluate the effect of female ovarian fluid on sperm velocity (a predictor of competitive fertilization success) according to whether males and females were from the same (native) or different (foreign) populations. Our results reveal the potential for ovarian fluid to act as a post‐mating–prezygotic barrier between two populations from different drainages, but also that the strength of this barrier is different among populations. This result may explain the previous finding that, in some populations, sperm from native males have precedence over foreign sperm, which could eventually lead to reproductive isolation between these populations.