EVOLUTIONARY IMPLICATIONS OF DIVERGENT CLINES IN AN AVIAN (MANACUS: AVES) HYBRID ZONE

Abstract A previous study of the hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus) documented the unidirectional introgression of vitellinus male secondary sexual traits across the zone. Here, we examine the hybrid zone in greater genetic and morphological detail. Statistical comparisons of clines are performed using maximum‐likelihood and nonparametric bootstrap methods. Our results demonstrate that an array of six molecular and two morphometric markers agree in cline position and width. Clines for male collar and belly color are similar in width to the first eight clines, but are shifted in position by at least five cline widths. The result is that birds in intervening populations are genetically and morphometrically very like parental candei, but males have the plumage color of parental vitellinus. Neither neutral diffusion nor nonlinearity of color scales appear to be viable explanations for the large cline shifts. Genetic dominance of vitellinus plumage traits is another potential explanation that will require breeding experiments to test. Sexual selection remains a plausible explanation for the observed introgression of vitellinus color traits in these highly dimorphic, polygynous, lek‐mating birds. Two other clines, including a nondiagnostic isozyme locus, are similar in position to the main cluster of clines, but are broader in width. Thus, introgression at some loci is greater than that detected with diagnostic markers. Assuming that narrow clines are maintained by selection, variation in cline width indicates that selection is not uniform throughout the genome and that diagnostic markers are under more intense selective pressure. The traditional focus on diagnostic markers in studies of hybrid zones may therefore lead to underestimates of average introgression. This effect may be more pronounced in organisms with low levels of genetic divergence between hybridizing taxa.     

Species‐specific habitat preferences do not shape the structure of a crested newt hybrid zone (Triturus cristatus x T. carnifex)

Reproductive isolation barriers maintain the integrity of species by preventing interspecific gene flow. They involve temporal, habitat or behavioral isolation acting before fertilization, and postzygotic isolation manifested as hybrid mortality or sterility. One of the approaches of how to study reproductive isolation barriers is through the analysis of hybrid zones. In this paper, we describe the structure of a hybrid zone between two crested newt species (Triturus cristatus and T. carnifex) in the southern part of the Czech Republic using morphological, microsatellite, and mitochondrial (mtDNA) markers. Specifically, we tested the hypothesis that the structure of the hybrid zone is maintained by species‐specific habitat preferences. Comparing the genetic structure of populations with geographical and ecological parameters, we found that the hybrid zone was structured primarily geographically, with T. cristatus‐like populations occurring in the northeast and T. carnifex‐like populations in the southwest. Despite T. cristatus tending to occur in deeper ponds and T. carnifex on localities with more shading, the effect of both ecological parameters on the structure of the zone was minimal. Next, we corroborated that T. carnifex individuals and some hybrids possess mtDNA of T. dobrogicus, whose nuclear background was not detected in the studied hybrid zone. Hybridization between T. carnifex and T. dobrogicus (resulting in unidirectional mtDNA introgression) had to predate subsequent formation of the hybrid zone between T. cristatus and T. carnifex. Populations of crested newts in the southern part of the Czech Republic thus represent a genetic mosaic of nuclear and mitochondrial genomes of three species.     

PATTERNS OF HYBRIDIZATION IN THE PIRIQUETA CAROLINIANA COMPLEX IN CENTRAL FLORIDA: EVIDENCE FOR AN EXPANDING HYBRID ZONE

We have identified a broad zone of hybridization between two morphologically and ecologically distinct herbaceous perennial taxa (morphotypes) within the Piriqueta caroliniana complex, which extends more than 300 km across the central Florida peninsula. Phylogeographic analyses indicate that the caroliniana morphotype has been present in north and central Florida since the early Pleistocene and that the viridis morphotype has immigrated into southern Florida much more recently. We examine the distribution of diagnostic morphological characters and nuclear genetic markers to assess the extent and patterns of introgression in this system. The results are consistent with the hypothesis that this hybrid zone has expanded north in recent history as viridis alleles have introgressed into regions that were previously occupied by populations of caroliniana. Genetic markers diagnostic for caroliniana have consistently high frequencies across the hybrid zone, whereas markers for viridis are extremely variable among populations with frequency reversals in adjacent populations. The latter pattern is probably the result of the combined stochastic effects of dispersal and drift on viridis alleles as they introgressed northward. Additional evidence for the recent expansion of this hybrid zone comes from patterns of variation for morphological and genetic markers. As expected for an expanding hybrid zone, within-population morphological variation was greatest toward the advancing front of introgression and levels of genetic variation for neutral diagnostic markers were greatest in the region of initial contact and lower in areas of recent expansion. The observed patterns of variation suggest that at least some hybrid genotypes have high fitnesses, which has led to the expansion of the hybrid zone via the displacement of parental genotypes in central Florida.     

Chemical communication,sexual selection,and introgression in wall lizards

Divergence in communication systems should influence the likelihood that individuals from different lineages interbreed, and consequently shape the direction and rate of hybridization. Here, we studied the role of chemical communication in hybridization, and its contribution to asymmetric and sexually selected introgression between two lineages of the common wall lizard (Podarcis muralis). Males of the two lineages differed in the chemical composition of their femoral secretions. Chemical profiles provided information regarding male secondary sexual characters, but the associations were variable and inconsistent between lineages. In experimental contact zones, chemical composition was weakly associated with male reproductive success, and did not predict the likelihood of hybridization. Consistent with these results, introgression of chemical profiles in a natural hybrid zone resembled that of neutral nuclear genetic markers overall, but one compound in particular (tocopherol methyl ether) matched closely the introgression of visual sexual characters. These results imply that associations among male chemical profiles, sexual characters, and reproductive success largely reflect transient and environmentally driven effects, and that genetic divergence in chemical composition is largely neutral. We therefore suggest that femoral secretions in wall lizards primarily provide information about residency and individual identity rather than function as sexual signals.     

Differential introgression of uniparentally inherited markers in bison populations with hybrid ancestries

Historical hybridization between Bison bison (bison) and Bos taurus (cattle) has been well documented and resulted in cattle mitochondrial DNA (mtDNA) introgression, previously identified in six different bison populations. In order to examine Y chromosome introgression, a microsatellite marker (BYM-1) with non-overlapping allele size distributions in bison and cattle was isolated from a bacterial artificial chromosome (BAC) clone, and was physically assigned to the Y chromosome by fluorescence in situ hybridization. BYM-1 genotypes for a sample of 143 male bison from 10 populations, including all six populations where cattle mtDNA haplotypes were previously identified, indicated that cattle Y chromosome introgression had not occurred in these bison populations. The differential permeability of uniparentally inherited markers to introgression is consistent with observations of sterility among first generation hybrid males and a sexual asymmetry in the direction of hybridization favouring matings between male bison and female cattle.     

Phylogeography of two European newt species--discordance between mtDNA and morphology

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5-1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long-term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.     

Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers

Evolutionary history of Muscicapidae flycatchers is inferred from nuclear and mitochondrial DNA (mtDNA) sequence comparisons and population genetic analysis of nuclear and mtDNA markers. Phylogenetic reconstruction based on sequences from the two genomes yielded similar trees with respect to the order at which the species split off. However, the genetic distances fitted a nonlinear, polynomial model reflecting diminishing divergence rate of the mtDNA sequences compared to the nuclear DNA sequences. This could be explained by Haldane's rule because genetic isolation might evolve more rapidly on the mitochondrial rather than the nuclear genome in birds. This is because hybrid sterility of the heterogametic sex (females) would predate that of the homogametic sex (males), leading to sex biased introgression of nuclear genes. Analyses of present hybrid zones of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) may indicate a slight sexual bias in rate of introgression, but the introgression rates were too low to allow proper statistical analyses. It is suggested, however, that the observed deviation from linearity can be explained by a more rapid mutational saturation of the mtDNA sequences than of the nuclear DNA sequences, as supported by analyses of third codon position transversions at two protein coding mtDNA genes. A phylogeographic scenario for the black and white flycatcher species is suggested based on interpretation of the genetic data obtained. Four species appear to have diverged from a common ancestor relatively simultaneously during the Pleistocene. After the last glaciation period, pied and collared flycatchers expanded their breeding ranges and eventually came into secondary contact in Central and Eastern Europe and on the Baltic Isles.     

Phylogenetic evidence of mitochondrial DNA introgression among pocket gophers in New Mexico (family Geomyidae)

Mitochondrial DNA (mtDNA) variation in the cytochrome b gene was determined for two divergent taxa of pocket gophers, Thomomys bottae actuosus and T. b. ruidosae . These two taxa hybridize in a narrow contact zone, but introgression of nuclear markers such as allozymes or chromosomes does not extend much beyond the hybrid zone (Patton et al. 1979). We found that despite their distinctness, the two subspecies shared very similar mtDNA haplotypes. By a comparison of phylogenetic histories derived from nuclear markers (allozymes) and from mtDNA haplotypes sampled in different populations of T. bottae from New Mexico, we show that apparent similarity is due to an introgression of T. b. ruidosae mtDNA into T. b. actuosus nuclear background. Evidence of introgression is not limited to the present-day contact zone between these two taxa, but extends at least 75 km away from it. The actuosus haplotype coexists along with the ruidosae mtDNA in the Gallinas Mts., which are inhabited by otherwise pure T. b. actuosus , while further north only typical actuosus haplotypes were detected. Of several potential mechanisms which could lead to such a geographical pattern of variation, we argue that a combination of range shifts due to climatic fluctuations, and genetic drift are most likely. Horizontal gene transfers due to hybridization are historical events which seem rather common among pocket gophers. Although they can be identified with careful phylogenetic study using independent data sets, the potential for misinterpreting a gene tree as an organismal tree is great in this and other groups of animals.     

Patterns of introgression across an expanding hybrid zone: analysing historical patterns of gene flow using nonequilibrium approaches

Previous studies suggest that the hybrid zone between two taxa in the Piriqueta caroliniana complex in central Florida (south-eastern North America) has recently expanded with hybrids replacing parental genotypes across a broad region of the Florida peninsula. Here I use patterns of genetic disequilibria and levels of differentiation among populations to infer historical patterns of introgression across this broad hybrid zone. There were strong positive associations among taxon-specific alleles at the southern extent of hybridization, but disequilibria values were close to zero across the central and northern portions of the hybrid zone. Levels of among-population differentiation in the central portion of the hybrid zone were relatively low, and increased towards the northern extent of hybridization. The high levels of disequilibria are coincident with the sharp clines at the southern end of the hybrid zone, suggesting that there is a tension zone in this region that is maintained by selection against hybrids and dispersal from parental regions. The levels of disequilibria within populations and patterns of differentiation among populations are consistent with historical introgression and northward expansion of this hybrid zone, which may have slowed or ceased in recent generations.     

Morphological and phylogenetic evidence for hybridization and introgression in a sea star secondary contact zone

Abstract. Glacial cycles and other climatic events have been widely invoked as factors promoting divergence, secondary contact, and hybridization between populations of terrestrial organisms, but the origin and fate of secondary contact in the sea is much less well understood. We studied the distribution of morphological and genetic variation in a northwest Atlantic zone of secondary contact between congeneric sea stars of Asterias that probably separated after the Pliocene as part of the trans-Arctic interchange. These species have similar reproductive biology and can hybridize in the laboratory. However, multivariate analysis of morphological traits scored from sea stars inside and outside the zone of secondary contact clearly indicated two clusters of phenotypes that corresponded to the two taxonomic species. A quantitative analysis of this clustering pattern did not support the hypothesis of a third grouping that might correspond to intermediate hybrid phenotypes. Known F₁ hybrids from laboratory matings grouped with one of the two taxonomic species. However, a survey of mtDNA sequence variation clearly indicated that ∼13% of individuals of one species ( Asterias forbesi ) are descendants of hybridization events that resulted in introgression of haplotypes of Asterias rubens into populations of A. forbesi . We conclude that morphological phenotypes are inadequate to identify hybrids of Asterias and their descendants, and that hybridization and introgression might be common in this secondary contact zone.     

Strong selection on male plumage in a hybrid zone between a hybrid bird species and one of its parents

Homoploid hybrid speciation (HHS) requires reproductive barriers between hybrid and parent species, despite incomplete reproductive isolation (RI) between the parents. Novel secondary sexual trait values in hybrids may cause prezygotic isolation from both parents, whereas signals inherited by the hybrid from one parent species may cause prezygotic isolation with the other. Here we investigate whether differences in male plumage function as a premating barrier between the hybrid Italian sparrow and one of its parent species, the house sparrow, in a narrow Alpine hybrid zone. Italian sparrow male plumage is a composite mosaic of the parental traits, with its head plumage most similar to its other parent, the Spanish sparrow. We use geographical cline analysis to examine selection on three plumage traits, 75 nuclear single nucleotide polymorphisms (SNPs) and hybrid indices based on these SNPs. Several SNPs showed evidence of restricted introgression in the Alps, supporting earlier findings. Crown colour exhibited the narrowest plumage cline, representing a 37% (range 4–65%) drop in fitness. The cline was too narrow to be due to neutral introgression. Only crown colour was significantly bimodal in the hybrid zone. Bimodality may be due to RI or a major QTL, although fitness estimates suggest that selection contributes to the pattern. We discuss the implications with respect to HHS and the species status of the Italian sparrow.     

Populations composed entirely of hybrid colonies: bidirectional hybridization and polyandry in harvester ants

In eusocial Hymenoptera, haplodiploid life cycles, obligate sterile castes, and polyandry may facilitate selection for hybridization. We analyzed a broad hybrid zone between the ecologically distinct seed‐harvester ants Pogonomyrmex occidentalis (Cresson) and Pogonomyrmex maricopa (Wheeler) using mitochondrial (mt)DNA sequence data, eight morphological markers, and 14 random amplified polymorphic DNA (RAPD) markers. Average mtDNA sequence divergence among parental species was 11.34%, indicating secondary contact. RAPD markers were significantly correlated with morphological variation, confirming the interspecific hybrid origin of all morphologically putative hybrid colonies. A morphological hybrid index indicates an abundance of both F₁ hybrids and parental morphotypes within colonies. Individual character frequencies plotted against distance show coincident and concordant clines, suggesting little to no introgression. The structure of the hybrid zone is two‐fold. Within the western region, stark reversals in character frequencies coincide with overt soil differences, indicating a mosaic hybrid zone structure. The eastern region is a riparian habitat where four adjacent populations were composed entirely of hybrid colonies. These habitat associations suggest that hybrid worker genomes permit dispersal into intermediate environments that select against one or both parental species. The present study suggests that, in addition to retaining reproductive compatibility, ecologically distinct species of ants may generate hybrid colonies maintained by environmental selection. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 320–336.     

Inferences from a rapidly moving hybrid zone

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e.. all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the 'wake" of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline movement is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter pure fatima genotypes. This movement can be explained if random dispersal of 7-19 km x gen(1/2) is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.     

PATTERNS OF VARIATION AND LINKAGE DISEQUILIBRIUM IN A FIELD CRICKET HYBRID ZONE

The distribution of multilocus genotypes found within a natural hybrid zone is determined by the sample of genotypes present when the hybrid zone first formed, by subsequent patterns of genetic exchange between the hybridizing taxa, and by drift and selection within each of the hybrid zone populations. We have used anonymous nuclear DNA restriction fragment polymorphisms (RFLPs) to characterize the array of multilocus genotypes present within a well-studied hybrid zone between two eastern North American field crickets, Gryllus pennsylvanicus and Gryllus firmus. These crickets hybridize along a zone of contact that extends from New England to Virginia. Previous studies have shown that both premating and postmating barriers exist between the two cricket species, but the absence of diagnostic morphological and allozyme markers has made it difficult to assess the consequences of these barriers for genetic exchange. Analyses based on four diagnostic anonymous nuclear markers indicate that hybrid zone populations in Connecticut contain few F₁ hybrids, and that nonrandom associations persist among nuclear gene markers, between nuclear and cytoplasmic markers, and between molecular markers and morphology. Field cricket populations within the hybrid zone are not “hybrid swarms” but consist primarily of crickets that are very much like one or the other of the parental species. Despite ample opportunity for genetic exchange and evidence for introgression at some loci, the two species remain quite distinct. Such a pattern appears to be characteristic of many natural hybrid zones.     

The establishment of a hybrid zone between red and sika deer (genus Cervus)

Japanese sika deer ( Cervus nippon nippon ) were introduced to Scotland around 80 years (20 generations) ago. The sika phenotype is expanding its range and hybridizing extensively with native red deer ( Cervus elaphus ) leading to the establishment of a hybrid zone. This zone is currently moving and cannot be considered to be at equilibrium. Cervid genotypes and mitochondrial haplotypes were mapped across the sika phenotype range, using diagnostic protein isozymes, microsatellite nuclear DNA markers and RFLPs in mtDNA. These were analysed to estimate heterozygote deficits and nuclear linkage disequilibria and cytonuclear disequilibria in relation to gene frequencies and time since contact. Introgression was found in both taxa and strong linkage disequilibria and heterozygote deficits characterize the populations longest exposed to hybridization. Populations further from the introduction site, where hybridization is facilitated by the dispersal of sika-like stags, show low values for linkage disequilibria and heterozygote deficit. The observed patterns in genotype are explained in terms of assortative mating and a selective advantage of the sika genotype. The genetic integrity of the Scottish mainland red deer is shown to be at risk from the invasion of sika.     

The modality of nine Triturus newt hybrid zones assessed with nuclear,mitochondrial and morphological data

The shape of hybrid zones provides insight into genetic isolation between the evolutionary lineages involved: the greater the limitation to gene flow, the farther along the unimodal to bimodal continuum. We study hybrid zone modality in a group of closely related species at a variety of levels of genetic divergence. We explore the degree and the geography of hybridization in 12 transects of Triturus newts, for nine of the ten secondary contact zones in Europe, using nuclear and mitochondrial encoded genetic data and morphological characteristics. Species status is verified by examining the correlation between the three marker systems in a geographical context. At 40 allozyme loci, two transects classify as unimodal, three as bimodal and two as intermediate, while the signal for hybridization is weak or absent in five transects. One zone studied in duplicate was classified as intermediate in one region and as bimodal in another region. mtDNA introgression is frequent and extends beyond nuclear introgression in two transects. Morphology provides additional evidence for hybridity, including one transect for which the signal of nuclear gene flow is weak. Compared to simulations allowing panmixia, the observed allozyme transitions at contact zones show a deficit of backcrossing to various degrees. Over all transects, there is a weak negative relationship between the level of hybridization and allozyme genetic distance for species pairs, consistent with Bateson–Dobzhansky–Muller effects. This observation, based upon highly comparable data for a single genus, supports conclusions derived from analyses over a wide variety of other taxa. © 2014 The Authors. Biological Journal of the Linnean Society published by John Wiley & Sons Ltd on behalf of The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 604–622.     

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.     

Morphological and genetic divergence of intralacustrine stickleback morphs in Iceland: a case for selective differentiation?

The evolutionary processes involved in population divergence and local adaptation are poorly understood. Theory predicts that divergence of adjacent populations is possible but depends on several factors including gene flow, divergent selection, population size and the number of genes involved in divergence and their distribution on the genome. We analyse variation in neutral markers, markers linked to putative quantitative trait loci and morphological traits in a recent (<10000 years) zone of primary divergence between stickleback morphs in Lake Thingvallavatn, Iceland. Environmental factors, especially predation, are clearly implicated in reducing gene flow between morphs. There is continuous morphological and genetic variation between habitats with a zone centre similar to secondary contact zones. Individual microsatellite loci are implicated as being linked to adaptive variation by direct tests as well as by differences in cline shape. Patterns of linkage disequilibria indicate that the morphs have diverged at several loci. This divergence shows parallels and differences with the well-studied limnetic-benthic stickleback morphs, both in phenotypic divergence and at the genomic level.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.