Genetic and morphological evidence of a geographically widespread hybrid zone between two crocodile species,Crocodylus acutus and Crocodylus moreletii

Hybrid zones represent natural laboratories to study gene flow, divergence and the nature of species boundaries between closely related taxa. We evaluated the level and extent of hybridization between Crocodylus moreletii and Crocodylus acutus using genetic and morphological data on 300 crocodiles from 65 localities. To our knowledge, this is the first genetic study that includes the entire historic range and sympatric zone of the two species. Contrary to expectations, Bayesian admixture proportions and maximum‐likelihood estimates of hybrid indexes revealed that most sampled crocodiles were admixed and that the hybrid zone is geographically extensive, extending well beyond their historical region of sympatry. We identified a few geographically isolated, nonadmixed populations of both parental species. Hybrids do not appear to be F₁s or recent backcrosses, but rather are more likely later‐generation hybrids, suggesting that hybridization has been going on for several to many generations and is mostly the result of natural processes. Crocodylus moreletii is not the sister species of C. acutus, suggesting that the hybrid zone formed from secondary contact rather than primary divergence. Nonadmixed individuals from the two species were distinguishable based on morphological characters, whereas hybrids had a complex mosaic of morphological characters that hinders identification in the wild. Very few nonadmixed C. acutus and C. moreletii populations exist in the wild. Consequently, the last nonadmixed C. moreletii populations have become critically endangered. Indeed, not only the parental species but also the naturally occurring hybrids should be considered for their potential conservation value.     

Introgressive hybridization and morphological transgression in the contact zone between two Mediterranean Solea species

Hybrid zones provide natural experiments where new combinations of genotypes and phenotypes are produced. Studying the reshuffling of genotypes and remodeling of phenotypes in these zones is of particular interest to document the building of reproductive isolation and the possible emergence of transgressive phenotypes that can be a source of evolutionary novelties. Here, we specifically investigate the morphological variation patterns associated with introgressive hybridization between two species of sole, Solea senegalensis and Solea aegyptiaca. The relationship between genetic composition at nuclear loci and individual body shape variation was studied in four populations sampled across the hybrid zone located in northern Tunisia. A strong correlation between genetic and phenotypic variation was observed among all individuals but not within populations, including the two most admixed ones. Morphological convergence between parental species was observed close to the contact zone. Nevertheless, the samples taken closest to the hybrid zone also displayed deviant segregation of genotypes and phenotypes, as well as transgressive phenotypes. In these samples, deviant body shape variation could be partly attributed to a reduced condition index, and the distorted genetic composition was most likely due to missing allelic combinations. These results were interpreted as an indication of hybrid breakdown, which likely contributes to postmating reproductive isolation between the two species.     

Hybridization following recent secondary contact results in asymmetric genotypic and phenotypic introgression between island species of Myzomela honeyeaters

Hybridization and introgression can have important evolutionary consequences for speciation, especially during early stages of secondary contact when reproductive barriers may be weak. Few studies, however, have quantified dynamics of hybridization and introgression in systems in which recent natural dispersal across a geographic barrier resulted in secondary contact. We investigated patterns of hybridization and introgression between two Myzomela honeyeaters (M. tristrami and M. cardinalis) that recently achieved secondary contact on Makira in the Solomon Islands. Hybridization in this system was hypothesized to be a byproduct of conspecific mate scarcity during early stages of colonization. Our research, however, provides evidence of ongoing hybridization more than a century after secondary contact. Mitochondrial sequencing revealed strongly asymmetric reproductive isolation that is most likely driven by postzygotic incompatibilities rather than prezygotic behavioral barriers. Nuclear introgression was observed from the native species (M. tristrami) to the colonizing species (M. cardinalis). Nuclear introgression in the reverse direction is almost exclusively limited to birds that are phenotypically M. tristrami but possess M. cardinalis mitochondrial haplotypes, consistent with introgression of plumage‐related alleles into the genomic background of M. cardinalis. These results provide unique insight into the dynamics and consequences of hybridization and introgression during early stages of secondary contact.     

Genetic structuring and secondary contact in the white‐chested Amazilia hummingbird species complex

Pleistocene climate cycles have been recognized to be a major driver of postglacial northward range expansion of North American bird populations. During glacial maxima, allopatric lineages that were reproductively isolated might have come into secondary contact with one another during expansion periods and the genetic signatures of past hybridization as a result of secondary contact events should produce detectable hybrid zones. The white‐chested hummingbirds Amazilia violiceps and A. viridifrons, constitute a species complex showing phenotypic similarity across its range. One exception is the subspecies found in the Central Depression of Chiapas (A. viridifrons villadai), which shares some plumage traits with the endemic but allopatric green‐fronted populations in Oaxaca. Phylogenetic relationships, taxonomy and species limits among violiceps, viridifrons and villadai have been controversial for decades. We assessed genetic structure of populations and introgression in this species complex by analysing 95 individuals at ten nuclear microsatellites and morphology. Bayesian analysis yielded four clusters. However, only two clusters generally match previously described mtDNA haplogroups, one parental taxon in the south (villadai) and a cluster with two admixed taxa (viridifrons and violiceps) that cannot be attributed to any pure parental population. High genetic admixture was recorded in the violiceps/viridifrons range, probably as a consequence of a postglacial northern expansion of violiceps. Signs of admixture and gene flow between violiceps/viridifrons and villadai were low. Historical and contemporary migration rates and Approximate Bayesian computations support a scenario of divergence with gene flow: a Pleistocene basal split separating A. violiceps and the other two clades are derived from a second split (villadai and viridifrons) or from a merger of violiceps and villadai into viridifrons due to gene flow.     

Differential introgression of a female competitive trait in a hybrid zone between sex‐role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male–male competition in maintaining or eroding species barriers. Although female–female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex‐role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female–female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex‐role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.     

Divergent character clines across a recent secondary contact zone in a Hispaniolan lizard

Studies of genetic contact zones provide valuable information regarding the processes of population divergence, adaptation and speciation. In this paper, I examine transitions in morphology, mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) haplotypes across a recent secondary contact zone in a Hispaniolan lizard Ameiva chrysolaema . Maximum likelihood cline fitting analyses suggest non-coincidence of cline centers and that the mtDNA cline is significantly displaced to the west of the remaining clines. nDNA and morphological clines are coincident and tend to be associated with the prevailing environmental gradient. The lack of cytonuclear disequilibrium near the center of the contact zone and the non-coincidence of character clines suggest that this zone does not conform to a tension zone model of hybridization; thus, gene flow across the zone does not seem to be impeded. The extremely narrow width of the dorsal scale size cline and the close association of this cline with the steepness of the environmental (precipitation) gradient suggest that this character may be under environmental selection. Taken together, this contact zone appears to be structured by a combination of mtDNA introgression, possibly associated with eastward movement of the zone, and environmental selection on some characters.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Hybridization between subspecies of waterbuck (Kobus ellipsiprymnus) in zones of overlap with limited introgression

Two subspecies of waterbuck (Kobus ellipsiprymnus), common (Kobus ellipsiprymnus ellipsiprymnus) and defassa (Kobus ellipsiprymnus defassa), are recognized based on differences in rump pattern, coat colour and geographical distribution. These forms are parapatrically distributed with an area of range overlap in East Africa, where phenotypically intermediate populations occur. Variation in 478 bp of the mitochondrial DNA control region and 14 polymorphic microsatellite loci were used to describe the genetic structure and phylogeographical pattern of the species, and to assess if the intermediate populations are the results of hybridization. In total, 186 individuals from 11 localities were analysed. A high degree of genetic differentiation was found between subspecies, although this was most evident from the microsatellite data. Hybridization was suggested in the phenotypically and geographically intermediate Nairobi NP population in Kenya. A neighbour-joining (NJ) tree based on microsatellite population genetic distances grouped Nairobi between the common and defassa populations, and a Bayesian analysis clearly showed introgression. Individuals sampled in Samburu NP, Kenya, had a common waterbuck phenotype, but introgression was suggested by both markers. Although a high degree of maternal defassa input was indicated from the sequence data, the Samburu population grouped with the common waterbuck in the microsatellite population genetic distance tree, with high support. Analyses of linkage disequilibrium and maximum-likelihood estimates of genetic drift suggested that admixture between subspecies is a recent event. The fact that introgression is limited between subspecies could be caused by chromosomal differences, hindering gene flow between common and defassa waterbuck.     

Evidence for historical introgression along a contact zone between two species of char (Pisces: Salmonidae) in northwestern North America

Phylogeographic analyses can yield valuable insights into the geographic and historical contexts of contact and hybridization between taxa. Two species of char (Salmonidae), Dolly Varden (Salvelinus malma) and bull trout (S. confluentus) have largely parapatric distributions in watersheds of northwestern North America. They are, however, sympatric in several localities and hybridization and some introgression occurs across a broad area of contact. We conducted a comparative phylogenetic analysis of Dolly Varden and bull trout to gain a historical perspective of hybridization between these species and to test for footprints of historical introgression. We resolved two major Dolly Varden mitochondrial DNA (mtDNA) clades (with 1.4-2.2% sequence divergence between haplotypes) that had different geographical distributions. Clade N is distributed across most of the range of Dolly Varden, from southern British Columbia through to the Kuril Islands in Asia. Clade S had a much more limited distribution, from Washington state, at the southern limit of the Dolly Varden range, to the middle of Vancouver Island. The distribution and inferred ages of the mtDNA clades suggested that Dolly Varden survived the Wisconsinan glaciation in a previously unsuspected refuge south of the ice sheet, and that Dolly Varden and bull trout were probably in continuous contact over most of the last 100,000 years. When bull trout were included in the phylogenetic analysis, however, the mtDNA of neither species was monophyletic: Clade S Dolly Varden clustered within the bull trout mtDNA clade. This pattern was discordant with two nuclear phylogenies produced (growth hormone 2 and rRNA internal transcribed sequence 1), in which Dolly Varden and bull trout were reciprocally monophyletic. This discordance between mtDNA- and nDNA-based phylogenies indicates that historical introgression of bull trout mtDNA into Dolly Varden occurred. Percent sequence divergence within these introgressed Dolly Varden (clade S) was 0.2-0.6%, implying that the introgression occurred prior to the most recent glaciation. Our analysis and other evidence of contact between divergent lineages in northwestern North America strongly suggests that the area may be the site of previously unsuspected suture zones of aquatic biotas.     

Appearance of morphological novelty in a hybrid zone between two species of land snail

On the small oceanic island of Chichijima, two endemic species of land snails, Mandarina mandarina and M. chichijimana, have discrete distributions separated by a hybrid zone. This study investigates the potential of hybridization as a source of morphological novelty in these snails. Mandarina mandarina possesses a shell with a higher whorl expansion rate and a smaller protoconch than M. chichijimana, relative to shell size. The number of whorls and shell size of M. mandarina do not differ from those of M. chichijimana, because the effect of higher expansion rate on number of whorls and size of the former is compensated for by its smaller protoconch. The whorl expansion rate and protoconch diameter of the individuals from the hybrid populations are intermediate or typical of either of the two species, and their average values show clinal changes along the hybrid zone. However, the hybrid populations include exceptionally high shells with many whorls and flat shells with few whorls, which are never found in the pure populations of either species. In addition, gradual increases in variance in shell height and number of whorls were found from the edges to the center of the hybrid zone. A combination of low expansion rate (typical of M. chichijimana) and a small protoconch (typical of M. mandarina) produces a shell with an extremely large number of whorls because of the geometry of shell coiling. However, the combination of high expansion rate and a large protoconch produces a shell with an extremely small number of whorls. Because of the correlation between the number of whorls and shell height, shells with an exceptional number of whorls possess an extraordinarily high or flat spire. Hybrids can inherit a mosaic of characters that, as they play out during growth, lead to novel adult morphologies. These findings emphasize the importance of hybridization as a source of morphological variation and evolutionary novelty in land snails.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Morphological and molecular evidence for hybridization and introgression in a willow (Salix) hybrid zone

Hybrid zones provide biologists with the opportunity to examine genetic and ecological interactions between differentiated populations. Accurate identification of hybrid genealogies is considered a necessary prerequisite to understanding observed patterns of hybridization-related phenomena. We analysed molecular and morphological data from individuals in a hybrid zone between two species of willows (Salix sericea Marshall and S. eriocephala Michaux) and report the use of randomly amplified polymorphic DNA (RAPD), chloroplast DNA (cpDNA), and ribosomal DNA (rDNA) markers, as well as vegetative morphology and foliar chemistry data to identify individuals in terms of hybrid genealogy and to infer the direction and extent of backcrossing and introgression within the hybrid zone. A novel version of a maximum likelihood estimate approach (developed for this study) was used to calculate hybrid index scores from RAPD marker data; this method produced results similar to those obtained using traditional arithmetic methods. Distribution of rDNA, cpDNA, and chemistry data were examined within the graphical context of RAPD-based hybrid index score histograms and principal component analyses (PCA) on RAPD and morphology data. Seven of the 21 plants classified as S. eriocephala in the field were possible introgressants. Another plant presented an unequivocal example of backcrossed S. sericea chemistry and RAPD markers. Inter- and intraspecific chloroplast diversity found within the hybrid zone suggests both historic introgression (perhaps in a glacial refugium), and contemporary hybridization. Patterns of inheritance and expression within the hybrid zone suggest that morphological characters are often not expressed in a simple additive fashion, and problems associated with both morphological and molecular data are considered.     

Reproductive isolation between oviparous and viviparous lineages of the Eurasian common lizard Zootoca vivipara in a contact zone

Contact zones between two evolutionary lineages are often useful for understanding the process of speciation because the observed genetic pattern reflects the history of differentiation. The Eurasian lacertid lizard Zootoca vivipara is a potentially interesting model for studying the role of reproductive mode in the speciation of squamate reptiles because it has both oviparous (Zootoca vivipara carniolica) and viviparous (Zootoca vivipara vivipara) populations that have recently been shown to be genetically distinct. We studied a newly‐discovered syntopic area of these two Zootoca subspecies in the central Italian Alps using genetic markers to investigate the level of introgression between them. Patterns of genetic differentiation in a fragment of the mitochondrial DNA cytb gene and a set of nuclear microsatellites show that the speciation process is complete in this area, with no evidence of recent introgression. Phylogenetic and genotypic divergence suggests that the two subspecies have experienced long independent evolutionary histories, during which genetic and phenotypic differences evolved. The possible roles of biogeography, reproductive mode, and cytogenetic differentiation in this speciation process are discussed. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 566–573.     

Interspecific competition promotes habitat and morphological divergence in a secondary contact zone between two hybridizing songbirds

Interspecific competition is assumed to play an important role in the ecological differentiation of species and speciation. However, empirical evidence for competition's role in speciation remains surprisingly scarce. Here, we studied the role of interspecific competition in the ecological differentiation and speciation of two closely related songbird species, the Common Nightingale (Luscinia megarhynchos) and the Thrush Nightingale (Luscinia luscinia). Both species are insectivorous and ecologically very similar. They hybridize in a secondary contact zone, which is a mosaic of sites where both species co‐occur (syntopy) and sites where only one species is present (allotopy). We analysed fine‐scale habitat data for both species in both syntopic and allotopic sites and looked for associations between habitat use and bill morphology, which have been previously shown to be more divergent in sympatry than in allopatry. We found that the two nightingale species differ in habitat use in allotopic sites, where L. megarhynchos occurred in drier habitats and at slightly higher elevations, but not in syntopic sites. Birds from allotopic sites also showed higher interspecific divergence in relative bill size compared to birds from syntopic sites. Finally, we found an association between bill morphology and elevation. Our results are consistent with the view that interspecific competition in nightingales has resulted in partial habitat segregation in sympatry and that the habitat‐specific food supply has in turn very likely led to bill size divergence. Such ecological divergence may enhance prezygotic as well as extrinsic postzygotic isolation and thus accelerate the completion of the speciation process.     

Spreading introgression in the wake of a moving contact zone

An increasing number of studies describe moving hybrid zones. This raises the issue of their actual frequency and emphasizes the need for methods that enable the detection of zone movements without historical records. Asymmetric introgression, usually considered as a signature of geographical shift, might be misleading when applied to mitochondrial or potentially non-neutral markers. We investigated mitochondrial and genomic introgression, using 30 AFLP derived markers, in a well-documented moving avian contact zone between two warblers. We found no instances of cross-species transmission of mitochondrial DNA but we detected nuclear introgression. Introgression levels were higher in the expanding species. Highest introgression was observed in populations that recently became allopatric than in current sympatric populations, which suggests that alien genetic material mainly spread at the time just before the receding species became extinct. We propose that either local recruitment or positive selection on some loci contribute to this pattern. Furthermore, we propose that, when the rarefaction of sexual partners drive the hybridization process, movement could be revealed by introgression peaking on the rear edge of the moving zone, or in its close vicinity.     

Distinguishing between two sympatric Acanthopagrus species from Dapeng Bay, Taiwan, using morphometric and genetic characters

Morphometric and genetic data were used to compare two sympatric and morphologically similar species, Acanthopagrus berda and Acanthopagrus taiwanensis , in Dapeng Bay, South-western Taiwan. A principle component analysis of morphological data indicated a distinction between the two species, with pectoral fin length and eye diameter accounting for 32·27% of the variation. Interspecific sequence divergence, based on mtDNA cytochrome b (0·118 ± 0·01), was larger than intraspecific divergences between haplotypes (0·007 for A. taiwanensis and 0·003 for A. berda ). Individuals of the two species clustered into different groups in three phylogenetic trees with 100% bootstrap support. The mean observed heterozygosity for eight microsatellite loci was 0·471 ± 0·202 for A. taiwanensis and 0·637 ± 0·145 for A. berda . Nei's unbiased measure of interspecific genetic distance ( D _S) was 1·334. F _ST (0·134) and R _ST (0·404) values indicated significant differentiation between species. An unrooted neighbour-joining tree was constructed by allele-sharing distances and the factorial correspondence analysis split all specimens into two distinct clusters. The results of morphometric, mtDNA and microsatellite analyses indicated the presence of two species, A. taiwanensis and A. berda .     

Bayesian analysis of hybridization and introgression between the endangered european mink (Mustela lutreola) and the polecat (Mustela putorius)

Human-mediated global change will probably increase the rates of natural hybridization and genetic introgression between closely related species, and this will have major implications for conservation of the taxa involved. In this study, we analyse both mitochondrial and nuclear data to characterize ongoing hybridization and genetic introgression between two sympatric sister species of mustelids, the endangered European mink (Mustela lutreola) and the more abundant polecat (M. putorius). A total of 317 European mink, 114 polecats and 15 putative hybrid individuals were collected from different localities in Europe and genotyped with 13 microsatellite nuclear markers. Recently developed Bayesian methods for assigning individuals to populations and identifying admixture proportions were applied to the genetic data. To identify the direction of hybridization, we additionally sequenced mtDNA and Y chromosomes from 78 individuals and 29 males respectively. We found that both hybridization and genetic introgression occurred at low levels (3% and 0.9% respectively) and indicated that hybridization is asymmetric, as only pure polecat males mate with pure European mink females. Furthermore, backcrossing and genetic introgression was detected only from female first-generation (F1) hybrids of European mink to polecats. This latter result implies that Haldane's rule may apply. Our results suggest that hybridization and genetic introgression between the two species should be considered a rather uncommon event. However, the current low densities of European mink might be changing this trend.     

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.     

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.