Variation in contact zone dynamics between two species of topminnows, Fundulus notatus and F. olivaceus, across isolated drainage systems

Spatially variable selection pressure within heterogeneous environments can result in the evolution of specialist phenotypes that facilitate co-occurrence of closely related species and limit genetic exchange. If divergent selection pressures maintain reproductive isolation, hybridization is expected to correlate with the strength of underlying ecological gradients and the traits shaped by adaptive processes. We sampled ten replicate topminnow (Fundulus olivaceus and Fundulus notatus) hybrid zones in isolated drainages throughout central and southern North America. In all drainages, species were distributed in an upstream–downstream manner with contact zones localized at confluences featuring abrupt shifts from tributary to river habitat. In two drainages, the typical up and downstream positions of species were reversed. Phenotype differences between the species reflect predicted selection differences along stream gradients. Downstream populations (lower food availability and greater predator pressure) generally showed larger investment in reproduction (higher gonadal somatic index), smaller body size and lower somatic condition compared to upstream populations. Phenotypic differences between the species in the two reversed drainages were consistent with convergence of life history traits in the respective habitats. Phenotypes of individuals of hybrid origin (F₁ hybrids or backcrosses) were not significantly different from the average of the two parental forms, though there were trends towards reduced fitness. The prevalence of hybridization among drainages ranged from no hybrids in two drainages to near random mating. The strongest correlates of hybridization rate among replicate hybrid zones were similarity in body shape and the homogeneity of habitat through tributary-river confluences. The two reversed orientation hybrid zones also exhibited high prevalence of hybrids suggesting that phenotypic convergence could lead to increased hybridization.     

Ecological and genetic assessment of spatial structure among replicate contact zones between two topminnow species

The spatial structure of contact zones is often described as disjunct, diffuse or mosaic and presumed to be related to underlying ecological gradients. However, the ecology of contact zones, how they are structured, and if that structure is predictable based on the strength and nature of ecological gradients is unknown. Large spatial scales and the unreplicated nature of many of the best studied contact zones has made it difficult to codify broader ecological patterns. Freshwater stream fish contact zones have the advantage of being potentially replicated with well defined boundaries and predictable linear gradients (river continuum concept). We sampled four replicate topminnow (Fundulus olivaceus and F. notatus) contact zones in Gulf of Mexico drainages. In each, we quantified contact zone spatial structure and the strength of ecological gradients (habitat, physicochemical variables and fish community functional traits). All three types of contact zone structure were represented. Systems with weaker gradients had diffuse contact zones, low species richness and were numerically dominated by generalist species. Rates of hybridization were also variable among systems. There was no hybridization detected in the mosaic zone while hybrids were found at most of the co-occurrence sites in the diffuse and disjunct zones. Overall, local ecology clearly influences contact zone structure and the two species interact in fundamentally different ways in these four systems.     

High gene flow despite opposite chirality in hybrid zones between enantiomorphic door snails

We studied differentiation and geneflow patterns between enantiomorphic door‐snail species in two hybrid zones in the Bucegi Mountains (Romania) to investigate the effects of intrinsic barriers (complications in copulation) and extrinsic selection by environmental factors. A mitochondrial gene tree confirmed the historical separation of the examined populations into the dextral Alopia livida and the sinistral Alopia straminicollis in accordance with the morphological classification, but also indicated gene flow between the species. By contrast, a network based on amplified fragment length polymorphisms (AFLP) markers revealed local groups of populations as units independent of their species affiliation. Admixture analyses based on AFLP data showed that the genomes of most individuals in the hybrid zones are composed of parts of the genomes of both parental taxa. The introgression patterns of a notable fraction of the examined markers deviated from neutral introgression. However, the patterns of most non‐neutral markers were not concordant between the two hybrid zones. There was also no concordance between non‐neutral markers in the two genomic clines and markers that were correlated with environmental variables or markers that were correlated with the proportion of dextral individuals in the populations. Neither extrinsic selection by environmental factors nor intrinsic barriers resulting from positive frequency‐dependent selection of the prevailing coiling direction were sufficient to maintain the distinctness of A. straminicollis and A. livida. Despite being historically separated units, we conclude that these taxa now merge where they come into contact.     

Shape variability in topminnows (Fundulus notatus species complex) along the river continuum

We examined intra‐ and interspecific variability in shape of three topminnow species (Funduluidae: Fundulus notatus, F. olivaceus, and F. euryzonus) across ten drainages. Within each drainage, five or more adult male topminnows were digitized at multiple sites (83 total sites) along the river continuum representing a range of stream sizes (cumulative drainage area) and hydrological conditions. Nine of the ten drainages contained two Fundulus species that were longitudinally separated along the river continuum with narrow areas of coexistence. Upstream–downstream distribution patterns were variable by drainage, allowing us to examine patterns repeated across ecologically similar species. More variability in shape was explained by drainage (19.7%) than by species (7.4%) differences. Populations of F. notatus from headwaters (three drainages) converged on a deep‐bodied form similar to F. olivaceus which was typically sampled in headwaters. Fundulus notatus shape was more closely related to stream size than in the other two species. Headwater populations of F. notatus and F. olivaceus had fineness ratios near the hydrodynamic optima of 4.5 whereas downstream populations of F. notatus had shallower bodies. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 612–621.     

Altitudinal gradients,plant hybrid zones and evolutionary novelty

Altitudinal gradients are characterized by steep changes of the physical and biotic environment that present challenges to plant adaptation throughout large parts of the world. Hybrid zones may form where related species inhabit different neighbouring altitudes and can facilitate interspecific gene flow and potentially the breakdown of species barriers. Studies of such hybrid zones can reveal much about the genetic basis of adaptation to environmental differences stemming from changes in altitude and the maintenance of species divergence in the face of gene flow. Furthermore, owing to recombination and transgressive effects, such hybrid zones can be sources of evolutionary novelty. We document plant hybrid zones associated with altitudinal gradients and emphasize similarities and differences in their structure. We then focus on recent studies of a hybrid zone between two Senecio species that occur at high and low altitude on Mount Etna, Sicily, showing how adaptation to local environments and intrinsic selection against hybrids act to maintain it. Finally, we consider the potential of altitudinal hybrid zones for generating evolutionary novelty through adaptive introgression and hybrid speciation. Examples of homoploid hybrid species of Senecio and Pinus that originated from altitudinal hybrid zones are discussed.     

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.     

Formation of a recent hybrid zone offers insight into the geographic puzzle and maintenance of species boundaries in musk turtles

Speciation is the result of an accumulation of reproductive barriers between populations, but pinpointing the factors involved is often difficult. However, hybrid zones can form when these barriers are not complete, especially when lineages come into contact in intermediate or modified habitats. We examine a hybrid zone between two closely related riverine turtle species, Sternotherus depressus and S. peltifer, and use dual‐digest RAD sequencing to understand how this hybrid zone formed and elucidate genomic patterns of reproductive isolation. First, the geographical extent and timing of formation of the hybrid zone is established to provide context for understanding the role of extrinsic and intrinsic reproductive isolating mechanisms in this system. The strength of selection on taxon‐specific contributions to maintenance of the hybrid zone is then inferred using a Bayesian genomic cline model. These analyses identify a role for selection inhibiting introgression in some genomic regions at one end of the hybrid zone and promoting introgression in many loci at the other. When selective pressures necessary to generate outliers to the genomic cline are considered with the geographical and temporal context of this hybrid zone, we conclude that habitat‐specific selection probably limits introgression from S. depressus to S. peltifer in the direction of river flow. However, selection is mediating rapid, unidirectional introgression from S. peltifer to S. depressus, which is probably facilitated by anthropogenic habitat alteration. These findings indicate a potentially imminent threat of population‐level genomic extinction for an already imperiled species due to ongoing human‐caused habitat alteration.     

Hybridization and reproductive isolation among syntopic populations of the topminnows Fundulus notatus and F. olivaceus

Fundulus notatus and Fundulus olivaceus are two closely related topminnow species that exhibit similar ecological niches and broad, largely overlapping, North American ranges extending throughout much of the Mississippi River drainage as well as the coastal drainages of the Gulf of Mexico. Previous studies have suggested that these two species are reproductively compatible despite cytogenetic differences and will hybridize when syntopic. We used nuclear and mtDNA loci to assess levels of hybridization and test for introgression in syntopic populations of these two species in four drainages in southern Illinois. Although hybridization was detected in all syntopic populations, an assessment of the proportion of hybrid individuals indicated a deficiency of hybrids relative to expectations under random mating. We determined that, although mtDNA introgression was prevalent and extended beyond the zones of contact, evidence of nuclear introgression was limited to the zone of sympatry.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Introgressive hybridization erodes morphological divergence between lentic and lotic habitats in an endangered minnow

Introgressive hybridization may erode phenotypic divergence along environmental gradients, collapsing locally adapted populations into a hybrid swarm. Alternatively, introgression may promote phenotypic divergence by providing variation on which natural selection can act. In freshwater fishes, water flow often selects for divergent morphological traits in lake versus stream habitats. We tested the effects of introgression on lake–stream morphological divergence in the minnow Owens Tui Chub (Siphateles bicolor snyderi), which has been rendered endangered by introgession from the introduced Lahontan Tui Chub (Siphateles bicolor obesa). Using geometric morphometric analysis of 457 individual Tui Chub from thirteen populations, we found that both native and introgressing parent taxa exhibited divergent body and caudal fin shapes in lake versus stream habitats, but their trajectories of divergence were distinct. In contrast, introgressed populations exhibited intermediate body and caudal fin shapes that were not differentiated by habitat type, indicating that introgression has eroded phenotypic divergence along the lentic–lotic gradient throughout the historic range of the Owens Tui Chub. Individuals within hybrid populations were less morphologically variable than those within parent populations, suggesting hybrid adaptation to selective agents other than water flow or loss of variance by drift.     

Differential introgression across newt hybrid zones: Evidence from replicated transects

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small‐scale (~2 km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.     

Tests of reproductive isolation among species in the <Emphasis Type="Italic">Fundulus notatus</Emphasis> (Cyprinodontiformes: Fundulidae) species complex

We assessed prezygotic (probability of spawning) and postzygotic (hatching success) reproductive isolation among the three ecologically and morphologically similar species in the Fundulus notatus species complex. We employed a multi-generation breeding experiment to test the hypotheses that karyotypic differences, body size differences, or geographic isolation among populations will increase pre or postzygotic reproductive barriers. Overall, prezygotic barriers were strong and postzygotic barriers weak in crosses of non-hybrid heterospecifics (F1 hybrid crosses) while prezygotic barriers were weaker and postzygotic barriers stronger in crosses involving hybrid individuals (F2 hybrid crosses and backcrosses). Prezygotic barriers among the two smaller species (Fundulus notatus and F. euryzonus) broke down rapidly; first generation hybrids spawned (F2 hybrid crosses and backcrosses) as frequently as parental forms in intraspecific crosses. There was no increase in postzygotic barriers among species with cytogenetic differences. There were increased prezygotic, but not postzygotic, barriers among geographically isolated populations of one species. While pure males and females were just as likely to spawn with hybrids, some types of hybrid females suffered from increased sterility, but not inviability, over hybrid males. Female sterility was only seen in hybrids with a Fundulus euryzonus parent, while other female hybrids produced viable eggs.     

Characterizing patterns of introgressive hybridization between two species of Tyrannus following concurrent range expansion

Hybridization between two species typically occurs when allopatric or ecologically dissimilar species expand into areas of secondary contact or habitat transitions. However, as species' ranges shift rapidly in response to environmental change, the potential for novel types of ephemeral hybrid zones exists. Here, we document and describe the occurrence, prevalence and symmetry of a previously undocumented hybrid zone involving two sympatric, ecologically similar sister species that have been expanding their ranges eastward in the central USA: Scissor‐tailed Flycatchers Tyrannus forficatus and Western Kingbirds Tyrannus verticalis. We identified cases of hybridization and introgression using analyses of eight microsatellite DNA loci and a single mitochondrial gene. We also evaluated short‐term reproductive consequences of hybridization for both species by surveying for both species and potential hybrids at the periphery of their ranges in northeastern Arkansas and western Tennessee, USA. Genetic data revealed bi‐directional backcrossing at the periphery of the species’ ranges, including a cryptic hybrid. We also analysed DNA of putative ‘pure’ individuals from other parts of their ranges and detected two cryptic admixed individuals, suggesting backcrossed individuals from the periphery may be dispersing to breed or that hybridization events have occurred in the core. Finally, our results suggest that there are no short‐term reproductive consequences of hybridization for the two species. In total, hybrid zones that occur at the edges of expanding, sympatric ranges may be ephemeral; we suggest they play an important role in introgression and may have long‐standing impacts for sympatric sister species. Exploring the extent of hybrid zones such as this for other range‐expanding taxa will elucidate whether this type of hybrid zone is unique or a common occurrence.     

HYBRIDIZATION IN WESTERN ATLANTIC STONE CRABS (GENUS MENIPPE): EVOLUTIONARY HISTORY AND ECOLOGICAL CONTEXT INFLUENCE SPECIES INTERACTIONS

Two distinct taxa of marine crabs (Menippe mercenaria and M. adina, sensu Williams and Felder [1986]) interbreed extensively in two disjunct areas of the coastal southeastern United States. We used variation in coloration and allele frequencies at three diagnostic loci to examine in detail the structure of the two zones of variation. A narrow hybrid zone in northwest Florida is situated at the junction of the present ranges of the two parental forms, in an area of strong ecological change and hydrological convergence. Despite extensive hybridization in this area, there is a significant deficiency of heterozygotes for the Alp-2 locus, nonrandom associations of alleles for two diagnostic loci, and an absence of certain combinations of phenotypes and genotypes. Along the Atlantic coast (east central Florida into South Carolina), a broad zone of increased variability exists within the range of M. mercenaria. Allele frequencies throughout this zone are similar to those of M. mercenaria but reflect apparent introgression from M. adina. In contrast, color patterns are quite variable, but only in the center of this zone. There is little evidence of a heterozygote deficiency, and the preferred habitat of M. mercenaria is not present. The Atlantic zone of variability is apparently expanding, with alleles at enzyme loci introgressing more rapidly than color characteristics. Despite these differences, certain features are common to both zones. These include 1) asymmetry in terms of the direction of introgression, 2) differential introgression of alleles, and 3) an almost complete absence of M. adina phenotypes that carry high proportions of M. mercenaria alleles. Differences between the two zones illustrate the influence that environmental setting, time of contact in relation to time of divergence, and location of the zone relative to the parental species ranges can have on hybridization events. However, observed similarities between the zones suggest that certain patterns of introgression and recombination may be independent of environmental setting. Thus, we suggest that factors inherent to the organism (intrinsic factors) and factors inherent to the environment (extrinsic factors) both act to structure and maintain the two hybrid zones.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long‐distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present‐day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.     

Extrinsic and intrinsic factors influence fitness in an avian hybrid zone

The effects of hybridization on evolutionary processes are primarily determined by the differential between hybrid and parental species fitness. Assessing the impacts of hybridization can be challenging, however, as determining the relationship between individual fitness and the extent of introgression in wild populations is difficult. We evaluated the fitness consequences of hybridization for pure and hybrid females in a hybrid zone between two tidal marsh birds, the saltmarsh sparrow (Ammodramus caudacutus), a salt marsh obligate, and Nelson's sparrow (A. nelsoni), which has a broader ecological niche and a much younger evolutionary association with salt marshes. Biotic stressors associated with nesting in tidal environments suggest an important role for differential adaptation in shaping hybrid zone dynamics, with saltmarsh sparrows predicted to be better adapted to nesting in salt marshes. We collected DNA samples from adults (n = 394) and nestlings (n = 431) to determine the extent of introgression using 12 microsatellite loci and tested for the influence of extrinsic (nest placement) and intrinsic (genotype) factors on female reproductive success. We monitored nests (n = 228), collected data on reproductive output, and estimated daily nest survival rates using female genotype and nest characteristics as covariates. To test for reduced survival of hybrid females, we also used capture data to assess the distribution of admixed male and female individuals across age classes. Reproductive success of females varied by genotypic class, but hybrids did not have intermediate success as predicted. Instead, we found that pure Nelson's sparrows had, on average, 33% lower hatching success than any other genotype, whereas F1/F2 hybrids, backcrossed Nelson's sparrows, and backcrossed and pure saltmarsh sparrows all had similar hatching success. We found no effect of genotype or nest placement on daily nest survival probabilities. However, hybrid individuals with a higher proportion of saltmarsh sparrow alleles exhibit nesting behaviours better suited to nesting successfully in tidal marshes. Further, while the proportion of F1/F2 individuals was similar between nestling and adult males, we found that the proportion of F1/F2 individuals was 2.3 times greater in nestling females compared with adult females, indicating reduced survival of F1 females. We conclude that differences in reproductive success among pure and admixed individuals coupled with intrinsic mechanisms (reduced survival in F1 females) shape hybrid zone dynamics in this system.     

Structure of a mosaic hybrid zone between the field crickets Gryllus firmus and G. pennsylvanicus

Hybrid zones provide insight into the nature of species boundaries and the evolution of barriers to gene exchange. Characterizing multiple regions within hybrid zones is essential for understanding both their history and current dynamics. Here, we describe a previously uncharacterized region of a well‐studied hybrid zone between two species of field crickets, Gryllus pennsylvanicus and G. firmus. We use a combination of mitochondrial DNA sequencing, morphological data, and modeling of environmental variables to identify the ecological factors structuring the hybrid zone and define patterns of hybridization and introgression. We find an association between species distribution and natural habitat; Gryllus pennsylvanicus occupies natural habitat along forest edges and natural clearings, whereas G. firmus occupies more disturbed areas in agricultural and suburban environments. Hybridization and introgression occur across patch boundaries; there is evidence of substantial admixture both in morphological characters and mtDNA, over a broad geographic area. Nonetheless, the distribution of morphological types is bimodal. Given that F₁ hybrids are viable and fertile in the lab, this suggests that strong pre‐zygotic barriers are operating in this portion of the hybrid zone.     

Widespread hybridization and bidirectional introgression in sympatric species of coral reef fish

Coral reefs are highly diverse ecosystems, where numerous closely related species often coexist. How new species arise and are maintained in these high geneflow environments have been long‐standing conundrums. Hybridization and patterns of introgression between sympatric species provide a unique insight into the mechanisms of speciation and the maintenance of species boundaries. In this study, we investigate the extent of hybridization between two closely related species of coral reef fish: the common coral trout (Plectropomus leopardus) and the bar‐cheek coral trout (Plectropomus maculatus). Using a complementary set of 25 microsatellite loci, we distinguish pure genotype classes from first‐ and later‐generation hybrids, identifying 124 interspecific hybrids from a collection of 2,991 coral trout sampled in inshore and mid‐shelf reefs of the southern Great Barrier Reef. Hybrids were ubiquitous among reefs, fertile and spanned multiple generations suggesting both ecological and evolutionary processes are acting to maintain species barriers. We elaborate on these finding to investigate the extent of genomic introgression and admixture from 2,271 SNP loci recovered from a ddRAD library of pure and hybrid individuals. An analysis of genomic clines on recovered loci indicates that 261 SNP loci deviate from a model of neutral introgression, of which 132 indicate a pattern of introgression consistent with selection favouring both hybrid and parental genotypes. Our findings indicate genome‐wide, bidirectional introgression between two sympatric species of coral reef fishes and provide further support to a growing body of evidence for the role of hybridization in the evolution of coral reef fishes.