Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise F_ST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.     

Patterns of reproductive isolation in a haplodiploid – strong post‐mating,prezygotic barriers among three forms of a social spider mite

In speciation research, much attention is paid to the evolution of reproductive barriers, preventing diverging groups from hybridizing back into one gene pool. The prevalent view is that reproductive barriers evolve gradually as a by‐product of genetic changes accumulated by natural selection and genetic drift in groups that are segregated spatially and/or temporally. Reproductive barriers, however, can also be reinforced by natural selection against maladaptive hybridization. These mutually compatible theories are both empirically supported by studies, analysing relationships between intensity of reproductive isolation and genetic distance in sympatric taxa and allopatric taxa. Here, we present the – to our knowledge – first comparative study in a haplodiploid organism, the social spider mite Stigmaeopsis miscanthi, by measuring premating and post‐mating, pre‐ and post‐zygotic components of reproductive isolation, using three recently diverged forms of the mite that partly overlap in home range. We carried out cross‐experiments and measured genetic distances (mitochondrial DNA and nuclear DNA) among parapatric and allopatric populations of the three forms. Our results show that the three forms are reproductively isolated, despite the absence of premating barriers, and that the post‐mating, prezygotic component contributes most to reproductive isolation. As expected, the strength of post‐mating reproductive barriers positively correlated with genetic distance. We did not find a clear pattern of prezygotic barriers evolving faster in parapatry than in allopatry, although one form did show a trend in line with the ecological and behavioural relationships between the forms. Our study advocates the versatility of haplodiploid animals for investigating the evolution of reproductive barriers.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

Multiple strong postmating and intrinsic postzygotic reproductive barriers isolate florally diverse species of Jaltomata (Solanaceae)

Divergence in phenotypic traits often contributes to premating isolation between lineages, but could also promote isolation at postmating stages. Phenotypic differences could directly result in mechanical isolation or hybrids with maladapted traits; alternatively, when alleles controlling these trait differences pleiotropically affect other components of development, differentiation could indirectly produce genetic incompatibilities in hybrids. Here, we determined the strength of nine postmating and intrinsic postzygotic reproductive barriers among 10 species of Jaltomata (Solanaceae), including species with highly divergent floral traits. To evaluate the relative importance of floral trait diversification for the strength of these postmating barriers, we assessed their relationship to floral divergence, genetic distance, geographical context, and ecological differences, using conventional tests and a new linear‐mixed modeling approach. Despite close evolutionary relationships, all species pairs showed moderate to strong isolation. Nonetheless, floral trait divergence was not a consistent predictor of the strength of isolation; instead this was best explained by genetic distance, although we found evidence for mechanical isolation in one species, and a positive relationship between floral trait divergence and fruit set isolation across species pairs. Overall, our data indicate that intrinsic postzygotic isolation is more strongly associated with genome‐wide genetic differentiation, rather than floral divergence.     

Intraspecific variation in reproductive barriers between two closely related Arabidopsis sister species

Reproductive isolation (RI) is a critical component of speciation and varies strongly in timing and strength among different sister taxa, depending on, for example the geography of speciation and divergence time. However, these factors may also produce variation in timing and strength among populations within species. Here we tested for variation in the expression of RI among replicate population pairs between the sister taxa Arabidopsis lyrata subsp. lyrata and A. arenicola. While the former is predominantly outcrossing, the latter is predominantly selfing. We focused on intrinsic prezygotic and postzygotic RI as both species occur largely in allopatry. We assessed RI by performing within-population crosses and interspecific between-population crosses, and by raising offspring. RI was generally high between all interspecific population pairs, but it varied in timing and strength depending on population history. Prezygotic isolation was strongest between the closest-related population pair, while early postzygotic isolation was high for all other population pairs. Furthermore, the timing and strength of RI depended strongly on cross direction. Our study provides empirical support that reproductive barriers between species are highly variable among population pairs and asymmetric within population pairs, and this variation seems to follow patterns typically described across species pairs.     

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.     

Allopatric origin of cryptic butterfly species that were discovered feeding on distinct host plants in sympatry

Surveys of tropical insects are increasingly uncovering cryptic species – morphologically similar yet reproductively isolated taxa once thought to comprise a single interbreeding entity. The vast majority of such species are described from a single location. This leaves us with little information on geographic range and intraspecific variation and limits our ability to infer the forces responsible for generating such diversity. For example, in herbivorous and parasitic insects, multiple specialists are often discovered within what were thought to be single more generalized species. Host shifts are likely to have contributed to speciation in these cases. But when and where did those shifts occur, and were they facilitated by geographic isolation? We attempted to answer these questions for two cryptic species within the butterfly Cymothoe egesta that were recently discovered on different host plants in central Cameroon. We first used mtDNA markers to separate individuals collected on the two hosts within Cameroon and then extended our analysis to incorporate individuals collected across the entire pan‐Afrotropical range of the original taxon. To our surprise, we found that the species are almost entirely allopatric, dividing the original range and overlapping only in the narrow zone of West‐Central Africa where they were first discovered in sympatry. This finding, combined with analyses of genetic variation within each butterfly species, strongly suggests that speciation occurred in allopatry, probably during the Pleistocene. We discuss the implications of our results for understanding speciation among other cryptic species recently discovered in the tropics and argue that more work is needed on geographic patterns and host usage in such taxa.     

Asymmetric reproductive barriers and mosaic reproductive isolation: insights from Misty lake–stream stickleback

Ecological speciation seems to occur readily but is clearly not ubiquitous – and the relative contributions of different reproductive barriers remain unclear in most systems. We here investigate the potential importance of selection against migrants in lake/stream stickleback (Gasterosteus aculeatus) from the Misty Lake system, Canada. This system is of particular interest because one population contrast (Lake vs. Outlet stream) shows very low genetic and morphological divergence, whereas another population contrast (Lake vs. Inlet stream) shows dramatic genetic and morphological divergence apparently without strong and symmetric reproductive barriers. To test whether selection against migrants might solve this “conundrum of missing reproductive isolation”, we performed a fully factorial reciprocal transplant experiment using 225 individually marked stickleback collected from the wild. Relative fitness of the different ecotypes (Lake, Inlet, and Outlet) was assessed based on survival and mass change in experimental enclosures. We found that Inlet fish performed poorly in the lake (selection against migrants in that direction), whereas Lake fish outperformed Inlet fish in all environments (no selection against migrants in the opposite direction). As predicted from their phenotypic and genetic similarity, Outlet and Lake fish performed similarly in all environments. These results suggest that selection against migrants is asymmetric and, together with previous work, indicates that multiple reproductive barriers contribute to reproductive isolation. Similar mosaic patterns of reproductive isolation are likely in other natural systems.     

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species‐level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near‐East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.     

The effects of heterospecific mating frequency on the strength of cryptic reproductive barriers

Heterospecific mating frequency is critical to hybrid zone dynamics and can directly impact the strength of reproductive barriers and patterns of introgression. The effectiveness of post‐mating prezygotic (PMPZ) reproductive barriers, which include reduced fecundity via heterospecific matings and conspecific sperm precedence, may depend on the number, identity and order of mates. Studies of PMPZ barriers suggest that they may be important in many systems, but whether these barriers are effective at realistic heterospecific mating frequencies has not been tested. Here, we evaluate the strength of cryptic reproductive isolation in two leaf beetles (Chrysochus auratus and C. cobaltinus) in the context of a range of heterospecific mating frequencies observed in natural populations. We found both species benefited from multiple matings, but the benefits were greater in C. cobaltinus and extended to heterospecific matings. We found that PMPZ barriers greatly limited hybrid production by C. auratus females with moderate heterospecific mating frequencies, but that their effectiveness diminished at higher heterospecific mating frequencies. In contrast, there was no evidence for PMPZ barriers in C. cobaltinus females at any heterospecific mating frequency. We show that integrating realistic estimates of cryptic isolation with information on relative abundance and heterospecific mating frequency in the field substantially improves our understanding of the strong directional bias in F1 production previously documented in the Chrysochus hybrid zone. Our results demonstrate that heterospecific mating frequency is critical to understanding the impact of cryptic post‐copulatory barriers on hybrid zone structure and dynamics, and that future studies of such barriers should incorporate field‐relevant heterospecific mating frequencies.     

Speciation in ninespine stickleback: reproductive isolation and phenotypic divergence among cryptic species of Japanese ninespine stickleback

Although similar patterns of phenotypic diversification are often observed in phylogenetically independent lineages, differences in the magnitude and direction of phenotypic divergence have been also observed among independent lineages, even when exposed to the same ecological gradients. The stickleback family is a good model with which to explore the ecological and genetic basis of parallel and nonparallel patterns of phenotypic evolution, because there are a variety of populations and species that are locally adapted to divergent environments. Although the patterns of phenotypic divergence as well as the genetic and ecological mechanisms have been well characterized in threespine sticklebacks, Gasterosteus aculeatus, we know little about the patterns of phenotypic diversification in other stickleback lineages. In eastern Hokkaido, Japan, there are three species of ninespine sticklebacks, Pungitius tymensis and the freshwater type and the brackish‐water type of the P. pungitius–P. sinensis species complex. They utilize divergent habitats along coast–stream gradients of rivers. Here, we investigated genetic, ecological and phenotypic divergence among three species of Japanese ninespine sticklebacks. Divergence in trophic morphology and salinity tolerance occurred in the direction predicted by the patterns observed in threespine sticklebacks. However, the patterns of divergence in armour plate were different from those previously found in threespine sticklebacks. Furthermore, the genetic basis of plate variation may differ from that in threespine sticklebacks. Because threespine sticklebacks are well‐established model for evolutionary research, the sympatric trio of ninespine sticklebacks will be an invaluable resource for ecological and genetic studies on both common and lineage‐specific patterns of phenotypic diversification.     

An integrative multidisciplinary approach to understanding cryptic divergence in Brazilian species of the Anastrepha fraterculus complex (Diptera: Tephritidae)

An integrative multidisciplinary approach was used to delimit boundaries among cryptic species within the Anastrepha fraterculus complex in Brazil. Sexual compatibility, courtship and sexual acoustic behaviour, female morphometric variability, variation for the mitochondrial gene COI, and the presence of Wolbachia were compared among A. fraterculus populations from the Southern (Vacaria, Pelotas, Bento Gonçalves, São Joaquim) and Southeastern (Piracicaba) regions of Brazil. Our results suggest full mating compatibility among A. fraterculus populations from the Southern region and partial pre‐zygotic reproductive isolation of these populations when compared with the population from the Southeastern region. A. fraterculus populations from both regions differed in the frequency of courtship displays and aspects of the calling phase and mounting acoustic signal. Morphometric analysis showed differences between Southern region and Southeastern region samples. All populations analyzed were infected with Wolbachia. The trees generated from the COI sequencing data are broadly congruent with the behavioural and morphometric data with the exception of one Southern population. The likely mechanisms by which A. fraterculus populations might have diverged are discussed in detail based on behavioural, morphometric, molecular genetics, and biogeographical studies.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Genetic divergence and reproductive barriers among morphologically heterogeneous sympatric clones of Eunotia bilunaris sensu lato (Bacillariophyta)

The study of reproductive isolation between populations, combined with estimates of genetic divergence, provides important insights into mechanisms of speciation. In this study, sixteen morphologically heterogeneous sympatric clones of Eunotia bilunaris sensu lato (Bacillariophyta) were brought into culture to study their phylogenetic relationships and pre- and postzygotic reproductive barriers. An ITS rDNA phylogeny was congruent with morphology and divided the clones into three groups ('slender', 'robust' and 'labile'), pointing to the presence of several species in E. bilunaris. Whereas most strains had a heterothallic mating system, four 'labile' clones displayed apomictic behaviour. A further 'labile' clone had a heterothallic mating behaviour, however, suggesting a very recent origin for apomixis. Despite high sequence divergence, hybridization occurred between clones belonging to different groups, but was 20-400 times less frequent than in intra-group matings. F1 hybrids had an intermediate morphology and were almost completely sterile; gamete formation was generally arrested in the early stages of meiosis I. The ITS divergence of 11.5-12.3% between the 'robust' and 'slender' clones seems to represent an upper limit of divergence in which cell pairing, gamete formation and auxosporulation are still possible but heavily reduced, and where hybrid sterility has already evolved.     

Vocal divergence is concordant with genomic evidence for strong reproductive isolation in grasshopper mice (Onychomys)

Behavioral barriers to gene flow often evolve faster than intrinsic incompatibilities and can eliminate the opportunity for hybridization between interfertile species. While acoustic signal divergence is a common driver of premating isolation in birds and insects, its contribution to speciation in mammals is less studied. Here we characterize the incidence of, and potential barriers to, hybridization among three closely related species of grasshopper mice (genus Onychomys). All three species use long‐distance acoustic signals to attract and localize mates; Onychomys arenicola and Onychomys torridus are acoustically similar and morphologically cryptic whereas Onychomys leucogaster is larger and acoustically distinct. We used genotyping‐by‐sequencing (GBS) to test for evidence of introgression in 227 mice from allopatric and sympatric localities in the western United States and northern Mexico. We conducted laboratory mating trials for all species pairs to assess reproductive compatibility, and recorded vocalizations from O. arenicola and O. torridus in sympatry and allopatry to test for evidence of acoustic character displacement. Hybridization was rare in nature and, contrary to prior evidence for O. torridus/O. arenicola hybrids, only involved O. leucogaster and O. arenicola. In contrast, laboratory crosses between O. torridus and O. arenicola produced litters whereas O. leucogaster and O. arenicola crosses did not. Call fundamental frequency in O. torridus and O. arenicola was indistinguishable in allopatry but significantly differentiated in sympatry, a pattern consistent with reproductive character displacement. These results suggest that assortative mating based on a long‐distance signal is an important isolating mechanism between O. torridus and O. arenicola and highlight the importance of behavioral barriers in determining the permeability of species boundaries.     

The accumulation of reproductive barriers during speciation: postmating barriers in two behaviorally isolated species of darters (Percidae: Etheostoma)

Identifying the manner in which reproductive barriers accumulate during lineage divergence is central to establishing general principles of species formation. One outstanding question is which isolating mechanisms form the first complete barrier to gene flow in a given lineage or under a particular set of conditions. To identify these initial reproductive barriers requires examining lineages in very early stages of divergence, before multiple reproductive barriers have evolved to completion. We quantified the strength of three postmating barriers in a pair of darter species and compared these estimates to each other and to the strength of behavioral isolation (BI) reported in a previous study. Results reveal no evidence of gametic incompatibility but intermediate levels of conspecific sperm precedence and hybrid inviability. As BI is nearly complete, our analysis comparing the strength of multiple reproductive barriers implicates the evolution of mate choice as central to both the origin and maintenance of these species. Further examination of ecological isolation and hybrid sterility is necessary to determine the role of these barriers in darter speciation.     

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

Sex-independent transmission ratio distortion system responsible for reproductive barriers between Asian and African rice species

* A sex-independent transmission ratio distortion (siTRD) system detected in the interspecific cross in rice was analyzed in order to understand its significance in reproductive barriers. The S(1) gene, derived from African rice Oryza glaberrima, induced preferential abortion of both male and female gametes possessing its allelic alternative (), from Asian rice O. sativa, only in the heterozygote. * The siTRD was characterized by resolving it into mTRD and fTRD occurring through male and female gametes, respectively, cytological analysis of gametophyte development, and mapping of the S(1) locus using near-isogenic lines. The allelic distribution of the S(1) locus in Asian and African rice species complexes was also analyzed. * The siTRD system involved at least two components affecting male and female gametogeneses, respectively, including a modifier(s) that enhances fTRD. The chromosomal location of the major component causing the mTRD was delimited within an approx. 40 kb region. The S(1) locus induced hybrid sterility in any pairwise combination between Asian and African rice species complexes. * The allelic state of the S(1) locus has diverged between Asian and African rice species complexes, suggesting that the TRD system has a significant role in the reproductive barriers in rice.     

Multiple post‐mating barriers to hybridization in field crickets

Mechanisms that prevent different species from interbreeding are fundamental to the maintenance of biodiversity. Barriers to interspecific matings, such as failure to recognize a potential mate, are often relatively easy to identify. Those occurring after mating, such as differences in the how successful sperm are in competition for fertilisations, are cryptic and have the potential to create selection on females to mate multiply as a defence against maladaptive hybridization. Cryptic advantages to conspecific sperm may be very widespread and have been identified based on the observations of higher paternity of conspecifics in several species. However, a relationship between the fate of sperm from two species within the female and paternity has never been demonstrated. We use competitive microsatellite PCR to show that in two hybridising cricket species, Gryllus bimaculatus and G. campestris, sequential cryptic reproductive barriers are present. In competition with heterospecifics, more sperm from conspecific males is stored by females. Additionally, sperm from conspecific males has a higher fertilisation probability. This reveals that conspecific sperm precedence can occur through processes fundamentally under the control of females, providing avenues for females to evolve multiple mating as a defence against hybridization, with the counterintuitive outcome that promiscuity reinforces isolation and may promote speciation.