THE MAINTENANCE OF A CLINE IN THE MARINE SNAIL LITTORINA SAXATILIS: THE ROLE OF HOME SITE ADVANTAGE AND HYBRID FITNESS

Steep clinal transitions in one or several inherited characters between genetically distinct populations are usually referred to as hybrid zones. Essentially two different mechanisms may maintain steep genetic clines. Either selection acts against hybrids that are unfit over the entire zone due to their mixed genetic origin (endogenous selection), or hybrids and parental types attain different fitness values in different parts of the cline (exogenous selection). Survival rate estimates of hybrids and parental forms in different regions of the cline may be used to distinguish between these models to assess how the cline is maintained. We used reciprocal transplants to test the relative survival rates of two parental ecotypes and their hybrids over microscale hybrid zones in the direct-developing marine snail Littorina saxatilis (Olivi) on the rocky shores of Galicia, Spain. One of the parental forms occupies upper and the other lower shores, and the hybrids are found at various proportions (1–38%) along with both parental forms in a midshore zone a few meters wide. The survival rate over one month was 39-52% of the native ecotype on upper shores, but only 2-8% for the lower-shore ecotype. In contrast, only 4-8% of the upper-shore ecotype but 53% of large (> 6 mm) and 8% of small (3-6 mm) native lower-shore ecotype survived in the lower shores. In the midshores, both the two parental ecotypes and the hybrids survived about equally well. Thus there is a considerable advantage for the native ecotypes in the upper and lower shores, while in the hybrid zone none of the morphs, hybrids included, are favored. This indicates that the dimorphism of L. saxatilis is maintained by steep cross-shore selection gradients, thus supporting the selection-gradient model of hybrid zones. We performed field and laboratory experiments that suggest physical factors and predation as important selective agents. Earlier studies indicate assortative mating between the two ecotypes in the midshore. This is unexpected in a hybrid zone maintained by selection gradients, and it seems as if the reproductive barrier compresses the hybrid zone considerably.     

Cytoplasmic and nuclear markers reveal contrasting patterns of spatial genetic structure in a natural Ipomopsis hybrid zone

Spatial variation in natural selection may play an important role in determining the genetic structure of hybridizing populations. Previous studies have found that F1 hybrids between naturally hybridizing Ipomopsis aggregata and Ipomopsis tenuituba in central Colorado differ in fitness depending on both genotype and environment: hybrids had higher survival when I. aggregata was the maternal parent, except in the centre of the hybrid zone where both hybrid types had high survival. Here, we developed both maternally (cpDNA PCR-RFLP) and biparentally inherited (nuclear AFLP) species-diagnostic markers to characterize the spatial genetic structure of the natural Ipomopsis hybrid zone, and tested the prediction that the majority of natural hybrids have I. aggregata cytoplasm, except in areas near the centre of the hybrid zone. Analyses of 352 individuals from across the hybrid zone indicate that cytoplasmic gene flow is bidirectional, but contrary to expectation, most plants in the hybrid zone have I. tenuituba cytoplasm. This cytotype distribution is consistent with a hybrid zone in historical transition, with I. aggregata nuclear genes advancing into the contact zone. Further, nuclear data show a much more gradual cline than cpDNA markers that is consistent with morphological patterns across the hybrid populations. A mixture of environment- and pollinator-mediated selection may contribute to the current genetic structure of this hybrid system.     

A test of the endogenous and exogenous selection hypotheses for the maintenance of a narrow avian hybrid zone

The contributions of genetic and environmental factors to differential reproductive success across hybrid zones have rarely been tested. Here, we report a manipulative experiment that simultaneously tested endogenous (genetic-based) and exogenous (environmental-based) selection within a hybrid zone. We transplanted mated pairs of two chickadee species (Poecile atricapilla and P. carolinensis) and their hybrids into isolated woodlots within their hybrid zone and monitored their reproductive success. Although clutch sizes were similar, based on an estimate of the genetic compatibility of a pair, hybrid pairs produced fewer nestlings and fledglings than did pairs of either parental species. According to a linear model generated from the data, a pure pair of either parental species would be expected to produce 1.91-2.48 times more fledglings per nesting attempt, respectively, than the average or least compatible pair in the experiment. Our result of decreased reproduction for hybrid pairs relative to parental species pairs within same environment (the hybrid zone in this experiment) support the endogenous selection hypothesis for maintenance of this hybrid zone. Because the experiment was conducted entirely within the hybrid zone (i.e., the same environment for parental and hybrid pairings), our data do not support the exogenous selection hypothesis as it predicts either all pairings doing poorly or the hybrid pairs more successful than the parental pairs.     

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.     

STABILITY OF THE NORTHERN FLICKER HYBRID ZONE IN HISTORICAL TIMES: IMPLICATIONS FOR ADAPTIVE SPECIATION THEORY

The hybrid zone between the Red- and Yellow-shafted Flickers has been stable on the United States Great Plains in historical times. This conclusion is based on multivariate comparisons of historical and contemporary collections from 18 locales. Adaptive speciation theory predicts that the hybrid zone should either become broader or narrower as a result of introgressive hybridization or reinforcement of premating isolating mechanisms. Neither of these predictions was borne out. Despite 10,000-13,000 years of hybridization, mating between subspecies remains indiscriminate. The data are also inconsistent with a dynamicequilibrium hypothesis wherein narrow hybrid zones are maintained by hybrid unfitness. According to this hypothesis, the hybrid zone would probably “flow” unless it was trapped by a population density trough. The hybrid zone does not appear to be associated with such a feature. The data are consistent with a bounded hybrid superiority theory of a hybrid zone, but this is more a question of survival in a process of elimination than a resounding corroboration.     

Estuarine habitats protect hybrid mussels from selection

The marine mussels, Mytilus edulis and Mytilus galloprovincialis, form an extensive hybrid zone in Europe where F₂ hybrids and mussels of mixed genetic ancestry are often locally abundant. Hybrid zones are maintained by the interplay of dispersal and selection on hybrid genotypes but there has been vigorous debate on the form of selection that may occur in these systems. Tension zone models argue that selection is against hybrids because of developmental misregulation and is independent of the external environment. Exogenous selection models argue that selection is habitat-dependent and the structure of the hybrid zone is due to the distribution of habitat patches that vary in selection intensity. We test this prediction by comparing the genetic structure of mussel populations in open-coast habitats, where selection on hybrids is strong, to those within two independent estuaries. We show that mussels within these estuaries are protected from selection and thus selection is strongly dependent on habitat, which supports the exogenous selection hypothesis. Hybrid mussel populations on the open-coast experience strong selection against M. edulis-like genotypes, which has been postulated to be the result of differential dislodgment by waves. This hypothesis is supported by our results since mussels within sheltered habitats are protected from selection. There was, despite previous suggestions, no evidence of selection in favor of M. edulis-like mussels within either estuary.     

Movements of adult fish in a hybrid zone revealed by microsatellite genetic analysis and capture-recapture data

• 1 Barbus barbus and Barbus meridionalis (Cyprinidae) form a hybrid zone in the Lergue river (southern France). We attempted to characterise the pattern of adult dispersal with the aim of understanding the evolutionary interactions maintaining this hybrid zone. In a previous capture/recapture study, movements between localities at a distance of 2.5 km or more appeared extremely rare. We therefore decided to investigate adult movement using a different method, based on genetic markers.
• 2 As this hybrid zone has produced allele frequency clines, massive population movements would produce temporal variations in genetic composition at a given locality. In order to determine the relationship between gene frequency and position, we surveyed four diagnostic or semidiagnostic microsatellite loci over the hybrid zone and estimated the cline produced by introgression. Then, we focused on a single locality at the centre of the hybrid zone and established its introgression index over 11 periods during two years.
• 3 The introgression index varied between periods producing significant Fst. A synthetic hybrid index, based on principal component analysis of the logit frequencies, was used for regression analyses and appeared significantly correlated with the river flow. This may be explained by displacement of adult fishes in response to flow increase. Using the information from the cline, we estimate that the most important gene frequency changes among dates, if created exclusively by population movement, correspond to a distance of 1500 m.
• 4 Additionally, we performed recaptures on a finer geographic scale than previously, around the central locality. No recaptured fishes were observed at 1875 m or more downstream, and at 875 m or more upstream, from the central locality. A high proportion of recaptured fishes (20% and 12%) was found in the nearest points (312 m downstream and 437 m upstream, respectively).
• 5 Thus, we established that (i) individual movements appear limited to a few hundred metres and (ii) individuals in a location tend to move in the same direction at the same time, probably in response to the same environmental factor.

     

Alternative forms for genomic clines

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory‐based logit‐logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one‐to‐one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package “HIest” for the R statistical software (CRAN; http://cran.r-project.org/ ).     

Asymmetrical introgression in a freshwater fish hybrid zone as revealed by a morphological index of hybridization

Two cyprinid species, Barbus barbus and Barbus meridionalis , form a hybrid zone in the River Lergue (southern France). Hybrids were distinguished morphologically by discriminant function analysis based on samples from the two parental taxa using only three measurements. This morphological hybridization index showed a strong asymmetry in introgression between the two taxa. B. meridionalis can live in the habitats typical of B. barbus and also in habitats where the latter cannot (small upstream tributaries), while remaining essentially unaffected by the hybridization process. By contrast, B. barbus is strongly introgressed because its habitat is invaded at least sporadically by B. meridionalis migrants. There is no behavioral barrier to hybridization, suggesting that an asymmetrical gene flow is probably partly responsible for the asymmetry in introgression.     

INCIPIENT SPECIATION DESPITE LITTLE ASSORTATIVE MATING: THE YELLOW‐RUMPED WARBLER HYBRID ZONE

Hybrid zones between recently diverged taxa are natural laboratories for speciation research, allowing us to determine whether there is reproductive isolation between divergent forms and the causes of that isolation. We present a study of a classic avian hybrid zone in North America between two subspecies of the yellow‐rumped warbler (Dendroica coronata). Although previous work has shown very little differentiation in mitochondrial DNA across this hybrid zone, we identified two nuclear loci (one sex‐linked and one autosomal) that show fixed differences across the hybrid zone, in a close concordance with patterns of plumage variation. Temporal stability and limited width of the hybrid zone, along with substantial linkage disequilibrium between these two diagnostic markers in the center of the zone, indicate that there is moderate reproductive isolation between these populations, with an estimated strength of selection maintaining the zone of 18%. Pairing data indicate that assortative mating is either very weak or absent, suggesting that this reproductive isolation is largely due to postmating barriers. Thus, despite extensive hybridization the two forms are distinct evolutionary groups carrying genes for divergent adaptive peaks, and this situation appears relatively stable.     

HYBRID DYSFUNCTION IN FIRE-BELLIED TOADS (BOMBINA)

Reproductive isolation between two taxa may be due to endogenous selection, which is generated by incompatibilities between the respective genomes, to exogenous selection, which is generated by differential adaptations to alternative environments, or to both. The continuing debate over the relative importance of either mode of selection has highlighted the need for unambiguous data on the fitness of hybrid genotypes. The hybrid zone between the fire-bellied toad (Bombina bombina) and the yellow-bellied toad (B. variegata) in central Europe involves adaptation to different environments, but evidence of hybrid dysfunction is equivocal. In this study, we followed the development under laboratory conditions of naturally laid eggs collected from a transect across the Bombina hybrid zone in Croatia. Fitness was significantly reduced in hybrid populations: Egg batches from the center of the hybrid zone showed significantly higher embryonic and larval mortality and higher frequencies of morphological abnormalities relative to either parental type. Overall mortality from day of egg collection to three weeks after hatching reached 20% in central hybrid populations, compared to 2% in pure populations. There was no significant difference in fitness between two parental types. Within hybrid populations, there was considerable variation in fitness, with some genotypes showing no evidence of reduced viability. We discuss the implications of these findings for our understanding of barriers to gene flow between species.     

Selective increase of a rare haplotype in a land snail hybrid zone

In hybrid zones, which are regions where genetically differentiated populations of organisms meet and produce hybrids, allozyme studies have often revealed unexpected alleles. The cause of this ''hybrizyme'' or ''rare allele'' phenomenon has been elusive, although it has been variously ascribed to natural selection or increased mutation rates. If the latter hypothesis is correct, selectively neutral markers should demonstrate increased variability in contrast to expressed markers such as allozymes. In this study, we screened selectively neutral variation in an intron of the calmodulin (CaM) gene in a hybrid zone between two subspecies of the Greek land snail Albinaria hippolyti. In previous allozyme studies, this hybrid zone has been shown to exhibit the rare allele phenomenon. We used a variant of the single-strand conformational polymorphism (SSCP) technique to detect seven haplotypes in both parental taxa. In the zone, one of these occurs at unexpectedly high frequencies. Since no additional mutants were found, we concluded that this is the result of selection.     

AN EMPIRICAL TEST OF PREDICTIONS OF TWO COMPETING MODELS FOR THE MAINTENANCE AND FATE OF HYBRID ZONES: BOTH MODELS ARE SUPPORTED IN A HARD-CLAM HYBRID ZONE

Two models developed to discern the mode of selection in hybrid zones differ in some predictions. The tension-zone model predicts that selection acts against hybrids and independently of the environment (endogenous selection) and that selection is invariant throughout the hybrid zone. The ecological selection-gradient, or ecotone, model maintains that fitness of different genotypes varies in response to environmental variation (exogenous selection) and thus, that in a region of the zone, fitness of hybrids is at least equal to that of the parental species. Therefore, to assess the predominant mode of selection operating in a hybrid zone, it is fundamental to evaluate whether selection is acting specifically against hybrid individuals, that is, whether hybridity alone is the basis for deficiencies of hybrids, and to evaluate whether the relative fitness of hybrids versus that of pure species varies across the zone. In a hardclam (genus Mercenaria) hybrid zone located in a polyhaline lagoon in east-central Florida, we used age-specific and location-specific analyses to determine that a hybrid deficit occurrs, that the deficit seems to be due to selection against hybrids, and that selection varies across the zone. Various measures of deviation from Hardy-Weinberg equilibrium, linkage disequilibrium analyses, and shifts in allele frequencies at semidiagnostic loci support the idea that selection is strongest in the northern region of the lagoon, the zone of sympatry and hybridization. Southward, into the range of M. mercenaria (the numerically predominant species), the percentage of hybrids remains relatively high and selection against hybrids decreases. For some genetic linkage groups, selection for M. mercenaria alleles seems to be occurring, but selection seems to be acting principally against alleles characteristic of M. mercenaria and, to a lesser degree, for alleles characteristic of M. campechiensis (the rarer species). These findings and others from previous analyses we have done on this hybrid zone demonstrate that selection in the zone is complex, and that characteristics of both the tension-zone and ecotone models are present. Supporting the tension-zone model, selection against hybrids per se clearly occurs, but specific genotypes seem to be at a selective disadvantage, whereas others have a selective advantage, and selection operates differentially on the two parental species within the zone. Supporting the ecotone model, the strength of overall selection varies throughout the zone, and environmentally mediated selection in which each species and hybrids have an advantage in specific habitats occurs, but some selection against hybrids is invariant throughout the zone. Thus, the structure and genetic architecture of this hybrid zone appear to be products of a complicated interaction between both types of selective forces cited in the two competing models.     

POPULATION STRUCTURE AND HOST USE IN HYBRIDIZING SUBSPECIES OF PAPILIO GLAUCUS (LEPIDOPTERA: PAPILIONIDAE)

Two parapatric subspecies of the eastern tiger swallowtail butterfly, Papilio glaucus glaucus and P. g. canadensis, differ greatly in larval survival and growth on host plants in the Magnoliaceae, Salicaceae, and Betulaceae. In the first part of this study, butterflies were collected from 17 sites along a transect from Georgia to Quebec and used for allozyme electrophoresis. Two X- (or Z-)linked enzyme loci (Ldh and Pgd) showed allele frequency differences that delineated a hybrid zone between the subspecies in northern Pennsylvania and south-central New York. No significant linkage disequilibrium could be detected among allozyme loci within the hybrid zone samples. Genetic differentiation at Ldh and Pgd coincides with subspecies differences in diapause control and female mimicry, which are also sex-linked. Larval offspring of butterflies from 13 of the sites were then tested in the laboratory for survival and growth on Liriodendron tulipifera (Magnoliaceae), Populus deltoïdes (Salicaceae) and Betula lutea (Betulaceae). Steep clines in survival rates matched the position of the hybrid zone. Hybrid zone larvae showed intermediate survival rates and significant heterogeneity among families on all three plants, indicating presence of substantial genetic variation. The results suggest that differential host use by P. g. glaucus and P. g. canadensis is maintained primarily by independent clines in selection intensity for ability to use allopatric sets of host plants, coupled with restricted gene flow across the hybrid zone.     

A new natural intergeneric cyprinid hybrid from the Jordan River drainage, with a key to the large barbine cyprinids of the southern Levant

In Ain al-Qunaiya. an isolated source within the Jordan River drainage basin of Jordan, several fishes were found to have characters intermediate between two cyprinid species. Barbus canis Valenciennes in Cuvier & Valenciennes. 1842 and Capoeta damascina (Valenciennes in Cuvier & Valenciennes. 1842). The hypothesis of their hybrid origin is advanced. A key to the large barbine cyprinids of the Jordan River system and their hybrids is given.     

Ontogenetic patterns of enzyme locus expression in Barbus hybrids (Cypriniformes, Teleostei)

The tissue specific patterns and ontogeny of sorbitol dehydrogenase (SDH, EC 1.1.1.14). lactate dehydrogenase (LDH, EC 1.1.1.27) and isocitrate dehydrogenase (IDH, EC 1.1.1.42) are reported for Barbus tetrazona (tiger barb), B. conchonius (rosy barb), B. nigrofasciatus (black ruby barb), B. titteya (cherry barb), B. sachsi (gold barb), and in interspecific hybrids where B. tetrazona is the maternal parent. The spatial and temporal expression of SDH, LDH and IDH isozymes in Barbus is consistent with those reported for other teleosts. As the genetic distance between the parentals used in forming the hybrid increases, allelic expression proceeds from synchronous to asynchronous, with an increasing delay in embryonic gene expression. These observations are consistent with the hypothesis that parental sensor genes differ in their response to maternally controlled regulatory signals; indicative of species specific effector/activator RNA molecule concentrations and sensor/receptor gene induction thresholds.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide

This study characterizes a hybrid zone that spans a migratory divide between subspecies of the Swainson's thrush (Catharus ustulatus), a long distance migratory songbird, in the Coast Mountains of British Columbia. To assess the potential for a barrier to gene flow between the subspecies, I: (1) analyzed the shape and width of genetic and morphological clines relative to estimates of dispersal distance, (2) assessed the ratio of parental to hybrid genotypes across the hybrid zone, (3) estimated population density across the hybrid zone, and (4) compared the spatial relationship between the hybrid zone and an existing environmental gradient. The results indicate that the hybrid zone is characterized by mostly concordant character clines that are narrow relative to dispersal, the absence of a hybrid swarm, and low population density at the center of the zone. This hybrid zone and additional regions of contact between these subspecies are found on the border between coastal and interior climatic regions throughout the Pacific Northwest. An identified shift in the location, but not the width, of the mtDNA cline relative to the nuclear clines is consistent with asymmetrical hybridization. Neutral diffusion of populations following secondary contact and hybrid superiority within an ecotone are insufficient explanations for the observed patterns. The hypothesis that best fits the data is that the Swainson's thrush hybrid zone is a tension zone maintained by dispersal and ecologically mediated barriers to gene flow.     

Exogenous selection shapes germination behaviour and seedling traits of populations at different altitudes in a Senecio hybrid zone

Background and Aims

The Senecio hybrid zone on Mt Etna, Sicily, is characterized by steep altitudinal clines in quantitative traits and genetic variation. Such clines are thought to be maintained by a combination of ‘endogenous’ selection arising from genetic incompatibilities and environment-dependent ‘exogenous’ selection leading to local adaptation. Here, the hypothesis was tested that local adaptation to the altitudinal temperature gradient contributes to maintaining divergence between the parental species, S. chrysanthemifolius and S. aethnensis.

Methods

Intra- and inter-population crosses were performed between five populations from across the hybrid zone and the germination and early seedling growth of the progeny were assessed.

Key Results

Seedlings from higher-altitude populations germinated better under low temperatures (9–13 °C) than those from lower altitude populations. Seedlings from higher-altitude populations had lower survival rates under warm conditions (25/15 °C) than those from lower altitude populations, but also attained greater biomass. There was no altitudinal variation in growth or survival under cold conditions (15/5 °C). Population-level plasticity increased with altitude. Germination, growth and survival of natural hybrids and experimentally generated F₁s generally exceeded the worse-performing parent.

Conclusions

Limited evidence was found for endogenous selection against hybrids but relatively clear evidence was found for divergence in seed and seedling traits, which is probably adaptive. The combination of low-temperature germination and faster growth in warm conditions might enable high-altitude S. aethnensis to maximize its growth during a shorter growing season, while the slower growth of S. chrysanthemifolius may be an adaptation to drought stress at low altitudes. This study indicates that temperature gradients are likely to be an important environmental factor generating and maintaining adaptive divergence across the Senecio hybrid zone on Mt Etna.     

Distortion of symmetrical introgression in a hybrid zone: evidence for locus-specific selection and uni-directional range expansion

The fate of species integrity upon natural hybridization depends on the interaction between selection and dispersal. The relative significance of these processes may be studied in the initial phase of contact before selection and gene flow reach equilibrium. Here we study a hybrid zone of two salamander species, Lyciasalamandra antalyana and Lyciasalamandra billae, at the initial phase of hybridization. We quantify the degree and mode of introgression using nuclear and mtDNA markers. The hybrid zone can be characterized as an abrupt transition zone, the central hybrid zone being only c. 400 m, but introgressed genes were traced up to 3 km. Introgression was traced in both sexes but gene flow may be slightly male-biased. Indirect evidence suggests that hybrid males are less viable than females. Introgression occurred at two levels: (1) locus-specific selection led to different allelic introgression patterns independent of species, while (2) asymmetrical species-level introgression occurred predominately from L. antalyana to L. billae due to range expansion of the former. This indicates that foreign genes can be incorporated into novel genomic environments, which in turn may contribute to the great diversity of morphological variants in Lyciasalamandra.