THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

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.     

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.     

GENOMIC AND MORPHOLOGICAL ANALYSIS OF A SEMIPERMEABLE AVIAN HYBRID ZONE SUGGESTS ASYMMETRICAL INTROGRESSION OF A SEXUAL SIGNAL

Hybrid zones are geographic regions where differentiated taxa meet and potentially exchange genes. Increasingly, genomic analyses have demonstrated that many hybrid zones are semipermeable boundaries across which introgression is highly variable. In some cases, certain alleles penetrate across the hybrid zone in only one direction, recombining into the alternate genome. We investigated this phenomenon using genomic (genotyping‐by‐sequencing) and morphological (plumage reflectance spectrophotometry) analyses of the hybrid zone between two subspecies of the red‐backed fairy‐wren (Malurus melanocephalus) that differ conspicuously in a sexual signal, male back plumage color. Geographic cline analyses revealed a highly variable pattern of differential introgression, with many narrow coincident clines combined with several significantly wider clines, suggesting that the hybrid zone is a semipermeable tension zone. The plumage cline was shifted significantly into the genomic background of the orange subspecies, consistent with sexual selection driving asymmetrical introgression of red plumage alleles across the hybrid zone. This interpretation is supported by previous experimental work demonstrating an extra‐pair mating advantage for red males, but the role of genetic dominance in driving this pattern remains unclear. This study highlights the potential for sexual selection to erode taxonomic boundaries and promote gene flow, particularly at an intermediate stage of divergence.     

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.     

Homoploid hybridization of plants in the Hengduan mountains region

The Hengduan Mountains Region (HMR) is a major global biodiversity hotspot. Complex tectonic and historical climatic conditions created opportunities for natural interspecific hybridization. Likewise, anthropogenic disturbance potentially raises the frequency of hybridization. Among species studies to date, the frequency of homoploid hybridization appears in the HMR. Of nine taxa in which natural hybridization has been detected, three groups are involved in homoploid hybrid speciation, and species pairs from the remaining six genera suggest that continuous gene flow occurs in hybrid zones. Reproductive isolation may greatly affect the dynamic and architecture of hybrid zones in the HMR. Asymmetrical hybridization and introgression can primarily be attributed to both prezygotic and postzygotic barriers. The frequent observation of such asymmetry may imply that reproductive barrier contributes to maintaining species boundaries in the alpine region. Ecological isolations with environmental disturbance may promote breeding barriers between parental species and hybrids. Hybrid zones may be an important phase for homoploid hybrid speciation. Hybrid zones potentially provided abundant genetic resources for the diversification of the HMR flora. The ecological and molecular mechanisms of control and mediation for natural hybridization will help biologists to understand the formation of biodiversity in the HMR. More researches from ecological and molecular aspects were required in future studies.     

Genetic differentiation and habitat connectivity across towhee hybrid zones in Mexico

Hybridization can either reinforce or erode species boundaries; therefore, hybrid zones offer a natural experimental setting in which to assess the dynamics of reproductive isolation. Secondary contact zones, in particular, present a partial separation of the original divergence mechanisms and the subsequent genomic architecture of reproductive isolation (or lack thereof). The spatial context of secondary contact and its consequent effect on dispersal play vital roles on the contact’s outcome. In a hybrid complex between two towhee species in Mexico, Pipilo maculatus and Pipilo ocai, two major hybrid gradients provide natural replicates for comparison. However, genetic analyses demonstrate significant divergence between geographically separate parental populations of each species and divergence of populations within each hybrid zone. The two hybrid transects (Teziutlán and Transvolcanic) are distinct and evidence suggests allelic introgression both across the species boundary and between the two transects. Habitat corridors for dispersal represent functional connectivity hotspots where the two transects meet. Both habitat connectivity and genetic differentiation between geographically disparate parental populations appear to influence the dynamics of gene flow across the hybrid gradients. In southern sympatric populations (Mt. Orizaba and Oaxaca) where morphological evidence for hybridization is scarce, opposing species’ alleles appear to traverse through the hybrid zones rather than arising from cryptic local hybridization. These results illustrate the importance of environmentally mediated gene flow in the context of secondary contact as an important force influencing evolutionary trajectory.     

Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

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.     

Variable patterns of introgression in two sculpin hybrid zones suggest that genomic isolation differs among populations

Theory predicts that reproductive isolation may be due to intrinsic genetic incompatibilities or extrinsic ecological factors. Therefore, an understanding of the genetic basis of isolation may require analyses of evolutionary processes in situ to include environmental factors. Here we study genetic isolation between populations of sculpins ( Cottus ) at 168 microsatellites. Genomic clines were fit using 480 individuals sampled across independent natural hybrid zones that have formed between one invading species and two separate populations of a resident species. Our analysis tests for deviations from neutral patterns of introgression at individual loci based on expectations given genome-wide admixture. Roughly 51% of the loci analysed displayed significant deviations. An overall deficit of interspecific heterozygotes in 26% and 21% of the loci suggests that widespread underdominance drives genomic isolation. At the same time, selection promotes introgression of almost 30% of the markers, which implies that hybridization may increase the fitness of admixed individuals. Cases of overdominance or epistatic interactions were relatively rare. Despite the similarity of the two hybrid zones in their overall genomic composition, patterns observed at individual loci show little correlation between zones and many fit different genotypic models of fitness. At this point, it remains difficult to determine whether these results are due to differences in external selection pressures or cryptic genetic differentiation of distinct parental populations. In the future, data from mapped genetic markers and on variation of ecological factors will provide additional insights into the contribution of these factors to variation in the evolutionary consequences of hybridization.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

Low intraspecific variation for genomic isolation between hybridizing sunflower species

Barriers to gene flow between species result from selection against foreign linkage blocks in hybrids. When the geographic ranges of taxa meet at multiple locations, the opportunity exists for variation in the genetic architecture of isolating barriers. Hybrid zones between two sunflower species (Helianthus annuus and H. petiolaris) in Nebraska and California exhibited remarkably similar patterns of introgression of mapped molecular markers. Congruence among hybrid zones may result from limited intraspecific variation at loci contributing to isolation and from similar selective effects of alleles in the heterospecific genetic background. The observed consistency of introgression patterns across distantly separated hybrid zones suggests that intrinsic forces predominate in determining hybrid zone dynamics and boundaries between these sunflower species.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION IN LOUISIANA IRISES: HYBRID FITNESS IN NATURE

Negative epistasis in hybrid genomes commonly results in postzygotic isolation between divergent lineages. However, some genomic regions may be selectively neutral or adaptive in hybrids and thus may potentially cross species barriers. We examined postzygotic isolation between ecologically similar species of Louisiana Iris: Iris brevicaulis and I. fulva to determine the potential for adaptive introgression in nature. Line-cross analyses allowed us a general overview of the gene action responsible for fitness-related traits. We then used a QTL mapping approach to detect genomic regions responsible for variation in these traits. Although hybrid classes suffered reduced fitness for many traits, hybrid means were equivalent to at least one of the parental species in overall estimates of maternal and paternal fitness during the two years of the field study. The genetic architecture underlying the fitness-related traits varied across field site and year of the study, thus emphasizing the importance of the environment in determining the degree of postzygotic isolation and potential for introgression across natural hybrid zones.     

The strength of reproductive isolation in two hybridizing food-deceptive orchid species

Reproductive isolation is of fundamental importance for maintaining species boundaries in sympatry. In orchids, the wide variety of pollination systems and highly diverse floral traits have traditionally suggested a prominent role for pollinator isolation, and thus for prezygotic isolation, as an effective barrier to gene flow among species. Here, we examined the nature of reproductive isolation between Anacamptis morio and Anacamptis papilionacea, two sister species of Mediterranean food-deceptive orchids, in two natural hybrid zones. Comparative analyses of the two hybrid zones that are located on soils with volcanic origin and have different and well-dated ages consistently revealed that all hybrid individuals were morphologically and genetically intermediate between the parental species, but had strongly reduced fitness. Molecular analyses based on nuclear ITS1 and (amplified fragment length polymorphism) AFLP markers clearly showed that all examined hybrids were F1 hybrids, and that no introgression occurred between parental species. The maternally inherited plastid DNA markers indicated that hybridization between A. morio and A. papilionacea was bidirectional, as confirmed by the molecular analysis of seed families. The genetic architecture of the two hybrid zones suggests that the two parental species easily and frequently hybridize in sympatry as a consequence of partial pollinator overlap but that strong postzygotic barriers reduce hybrid fitness and prevent gene introgression. These results corroborate that chromosomal divergence is instrumental for reproductive isolation between these food-deceptive orchids and suggest that hybridization is of limited importance for their diversification.     

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.     

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.     

Parallel introgression and selection on introduced alleles in a native species

As humans cause the redistribution of species ranges, hybridization between previously allopatric species is on the rise. Such hybridization can have complex effects on overall fitness of native species as new allelic combinations are tested. Widespread species introductions provide a unique opportunity to study selection on introgressed alleles in independent, replicated populations. We examined selection on alleles that repeatedly introgressed from introduced rainbow trout (Oncorhynchus mykiss) into native westslope cutthroat trout (Oncorhynchus clarkii lewisi) populations in western Canada. We found that the degree of introgression of individual single nucleotide polymorphisms from the invasive species into the native is correlated between independent watersheds. A number of rainbow trout alleles have repeatedly swept to high frequency in native populations, suggesting parallel adaptive advantages. Using simulations, we estimated large selection coefficients up to 0.05 favoring several rainbow trout alleles in the native background. Although previous studies have found reduced hybrid fitness and genome‐wide resistance to introgression in westslope cutthroat trout, our results suggest that some introduced genomic regions are strongly favored by selection. Our study demonstrates the utility of replicated introductions as case studies for understanding parallel adaptation and the interactions between selection and introgression across the genome. We suggest that understanding this variation, including consideration of beneficial alleles, can inform management strategies for hybridizing species.     

Tracing species replacement in Iberian marbled newts

Secondary contact between closely related species can lead to the formation of hybrid zones, allowing for interspecific gene flow. Hybrid zone movement can take place if one of the species possesses a competitive advantage over the other, ultimately resulting in species replacement. Such hybrid zone displacement is predicted to leave a genomic footprint across the landscape in the form of asymmetric gene flow (or introgression) of selectively neutral alleles from the displaced to the advancing species. Hybrid zone movement has been suggested for marbled newts in the Iberian Peninsula, supported by asymmetric gene flow and a distribution relict (i.e., an enclave) of Triturus marmoratus in the range of T. pygmaeus. We developed a panel of nuclear and mitochondrial SNP markers to test for the presence of a T. marmoratus genomic footprint in the Lisbon peninsula, south of the enclave. We found no additional populations of T. marmoratus. Analysis with the software Structure showed no genetic traces of T. marmoratus in T. pygmaeus. A principal component analysis showed some variation within the local T. pygmaeus, but it is unclear if this represents introgression from T. marmoratus. The results may be explained by (a) species replacement without introgressive hybridization and (b) displacement with hybridization followed by the near‐complete erosion of the footprint by purifying selection. We predict that testing for a genomic footprint north of the reported enclave would confirm that species replacement in these marbled newts occurred with hybridization.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

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.