DNA在鸟类分子系统发育研究中的应用

鸟类分子系统发育研究中常用的DNA技术有DNA杂交、RFLP和DNA序列分析等。DNA杂交技术曾在鸟类中有过大规模的应用,并由此诞生了一套新的鸟类分类系统。在鸟类的RFLP分析中,用的最多的靶序列是线粒体DNA。DNA序列分析技术被认为是进行分子系统发育研究最有效、最可靠的方法。在DNA序列分析中,线粒体基因应用最广泛,但由于其自身的一些不足,近年来,不少学者把目光投向了核基因,将线粒体基因和核基因结合起来进行系统发育研究。目前在鸟类分子系统发育中,应用较多的核基因是scnDNA,其内含子可以用于中等阶元水平的系统研究,而外显子主要用于高等阶元的系统研究。除了分子标记自身的问题之外,鸟类分子系统发育研究中还存在着方法上的问题,包括分子标记的选择,样本数量以及数据处理等。今后鸟类分子系统发育研究应该更加注重方法的标准化。     

Interspecific hybridization causes long‐term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes

Classification, phylogeography and the testing of evolutionary hypotheses rely on correct estimation of species phylogeny. Early molecular phylogenies often relied on mtDNA alone, which acts as a single linkage group with one history. Over the last decade, the use of multiple nuclear sequences has often revealed conflict among gene trees. This observation can be attributed to hybridization, lineage sorting, paralogy or selection. Here, we use 54 groups of fishes from 48 studies to estimate the degree of concordance between mitochondrial and nuclear gene trees in two ecological grades of fishes: marine and freshwater. We test the hypothesis that freshwater fish phylogenies should, on average, show more discordance because of their higher propensity for hybridization in the past. In keeping with this idea, concordance between mitochondrial and nuclear gene trees (as measured by proportion of components shared) is on average 50% higher in marine fishes. We discuss why this difference almost certainly results from introgression caused by greater historical hybridization among lineages in freshwater groups, and further emphasize the need to use multiple nuclear genes, and identify conflict among them, in estimation of species phylogeny.     

Unexpected population genetic structure of European roe deer in Poland: an invasion of the mtDNA genome from Siberian roe deer

Introgressive hybridization is a widespread evolutionary phenomenon which may lead to increased allelic variation at selective neutral loci and to transfer of fitness‐related traits to introgressed lineages. We inferred the population genetic structure of the European roe deer (Capreolus capreolus) in Poland from mitochondrial (CR and cyt b) and sex‐linked markers (ZFX, SRY, DBY4 and DBY8). Analyses of CR mtDNA sequences from 452 individuals indicated widespread introgression of Siberian roe deer (C. pygargus) mtDNA in the European roe deer genome, 2000 km from the current distribution range of C. pygargus. Introgressed individuals constituted 16.6% of the deer studied. Nearly 75% of them possessed haplotypes belonging to the group which arose 23 kyr ago and have not been detected within the natural range of Siberian roe deer, indicating that majority of present introgression has ancient origin. Unlike the mtDNA results, sex‐specific markers did not show signs of introgression. Species distribution modelling analyses suggested that C. pygargus could have extended its range as far west as Central Europe after last glacial maximum. The main hybridization event was probably associated with range expansion of the most abundant European roe deer lineage from western refugia and took place in Central Europe after the Younger Dryas (10.8–10.0 ka BP). Initially, introgressed mtDNA variants could have spread out on the wave of expansion through the mechanism of gene surfing, reaching high frequencies in European roe deer populations and leading to observed asymmetrical gene flow. Human‐mediated introductions of C. pygargus had minimal effect on the extent of mtDNA introgression.     

The phylogeny of the hominoid primates,as indicated by DNA-DNA hybridization

Summary The living hominoid primates are Man, the chimpanzees, the Gorilla, the Orangutan, and the gibbons. The cercopithecoids (Old World monkeys) are the sister group of the hominoids. The composition of the Hominoidea is not in dispute, but a consensus has not yet been reached concerning the phylogenetic branching pattern and the dating of divergence nodes. We have compared the single-copy nuclear DNA sequences of the hominoid genera using DNA-DNA hybridization to produce a complete matrix of delta T₅₀H values. The data show that the branching sequence of the lineages, from oldest to most recent, was: Old World monkeys, gibbons, Orangutan, Gorilla, chimpanzees, and Man. The calibration of the delta T₅₀H scale in absolute time needs further refinement, but the ranges of our estimates of the datings of the divergence nodes are: Cercopithecoidea, 27–33 million years ago (MYA); gibbons, 18–22 MYA; Orangutan, 13–16 MYA; Gorilla, 8–10 MYA; and chimpanzees-Man, 6.3–7.7 MYA.     

Cryptic diversity and unexpected evolutionary patterns in the meadow lizard,Darevskia praticola (Eversmann, 1834)

Darevskia praticola differs from the other species of the genus in having a large but disjunct distribution, covering the Balkan and the Caucasus regions. Furthermore, most Darevskia species occupy saxicolous habitats, whereas D. praticola inhabits meadows and forest environments. Here we determine the phylogeographic and phylogenetic relationships of Darevskia praticola sensu lato and evaluate the current, morphology-based taxonomy. We sequenced two mtDNA genes (Cyt-b and ND4) and two nuclear loci (MC1R and RELN) for samples collected across the species range. Because our sequences amplified with the Cyt-b primers appear to represent a nuclear pseudogene we excluded this marker from the final analysis. Our results support monophyly of D. praticola and show its division into three clades. The first divergence, dated to the Late Pliocene, is between the Balkans and the Caucasus. The Caucasus lineage is further subdivided in a western Greater Caucasus and a Transcaucasia clade, likely due to subsequent differentiation during the Pleistocene. Our findings do not support the current taxonomic arrangement within D. praticola. The main geographic divergence likely happened due to a vicariance event associated with Plio-Pleistocene climatic and vegetation oscillations.     

Introgressive hybridization and natural selection in Darwin's finches

Introgressive hybridization, i.e. hybridization with backcrossing, can lead to the fusion of two species, but it can also lead to evolution of a new trajectory through an enhancement of genetic variation in a new or changed ecological environment. On Daphne Major Island in the Galápagos archipelago, ~1–2% of Geospiza fortis finches breed with the resident G. scandens and with the rare immigrant species G. fuliginosa in each breeding season. Previous research has demonstrated morphological convergence of G. fortis and G. scandens over a 30‐year period as a result of bidirectional introgression. Here we examine the role of hybridization with G. fuliginosa in the evolutionary trajectory of G. fortis. Geospiza fuliginosa (~12 g) is smaller and has a more pointed beak than G. fortis (~17 g). Genetic variation of the G. fortis population was increased by receiving genes more frequently from G. fuliginosa than from G. scandens (~21 g). A severe drought in 2003–2005 resulted in heavy and selective mortality of G. fortis with large beaks, and they became almost indistinguishable morphologically from G. fuliginosa. This was followed by continuing hybridization, a further decrease in beak size and a potential morphological fusion of the two species under entirely natural conditions.     

The biogeography of mitochondrial and nuclear discordance in animals

Combining nuclear (nuDNA) and mitochondrial DNA (mtDNA) markers has improved the power of molecular data to test phylogenetic and phylogeographic hypotheses and has highlighted the limitations of studies using only mtDNA markers. In fact, in the past decade, many conflicting geographic patterns between mitochondrial and nuclear genetic markers have been identified (i.e. mito-nuclear discordance). Our goals in this synthesis are to: (i) review known cases of mito-nuclear discordance in animal systems, (ii) to summarize the biogeographic patterns in each instance and (iii) to identify common drivers of discordance in various groups. In total, we identified 126 cases in animal systems with strong evidence of discordance between the biogeographic patterns obtained from mitochondrial DNA and those observed in the nuclear genome. In most cases, these patterns are attributed to adaptive introgression of mtDNA, demographic disparities and sex-biased asymmetries, with some studies also implicating hybrid zone movement, human introductions and Wolbachia infection in insects. We also discuss situations where divergent mtDNA clades seem to have arisen in the absence of geographic isolation. For those cases where foreign mtDNA haplotypes are found deep within the range of a second taxon, data suggest that those mtDNA haplotypes are more likely to be at a high frequency and are commonly driven by sex-biased asymmetries and/or adaptive introgression. In addition, we discuss the problems with inferring the processes causing discordance from biogeographic patterns that are common in many studies. In many cases, authors presented more than one explanation for discordant patterns in a given system, which indicates that likely more data are required. Ideally, to resolve this issue, we see important future work shifting focus from documenting the prevalence of mito-nuclear discordance towards testing hypotheses regarding the drivers of discordance. Indeed, there is great potential for certain cases of mitochondrial introgression to become important natural systems within which to test the effect of different mitochondrial genotypes on whole-animal phenotypes.     

Natural hybridization and gene flow between twotridentiger gobies in lake hinuma (ibaraki prefecture, japan)

Of 851 specimens ofTridentiger obscurus andt. brevispinis collected from Lake Hinuma (Ibaraki Prefecture, Japan) from July 1996 to February 1998, 49 (5.8%) comprised F₁ hybrids and backcross progeny of the two species. Since the mitochondrial DNA haplotypes of the F₁ hybrids reflected those ofT. brevispinis, most instances of hybridization are thought to have occurred between maleT. obscurus and femaleT. brevispinis. Although both allozyme and mitochondrial DNA analyses indicated backcrossing and introgression of mitochondrial DNA, the frequency of backcross progeny was relatively low, suggesting the existence of a natural selection to backrossing.     

Interspecific differentiation and intraspecific substructure in two closely related clupeids with extensive hybridization, Alosa alosa and Alosa fallax

In this study, the evolutionary relationships within and among populations of European shads, Alosa alosa and Alosa fallax , was investigated. Screening of allelic variation across eight allozyme loci and sequencing 448 bp of the mtDNA cytochrome b gene in 14 rivers throughout the range of the species supported that the two taxa were independent lineages (1·3% net nucleotide divergence) despite extensive hybridization. Genetic diversity and structure was considerably higher in A. fallax than A. alosa and the former species revealed evidence of distinct lineages in the Mediterranean and Atlantic basins. A Bayesian clustering approach combined with gill raker counts verified that individuals of the two species could be assigned to their parent group with relatively high confidence. Evaluation of hybridization in the Lima and Mondego Rivers in Portugal provided evidence that introgression is extensive but is not currently obscuring (through hybrid swarming) the diagnosability of the two species.     

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.     

A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes

Phylogenetic relationships of members of the subfamily Poeciliinae (Cyprinodontiformes) are investigated to test alternate hypotheses of diversification resulting from the assembly of the Central America and the Caribbean from the Cretaceous period onwards. We use 4333 aligned base pairs of mitochondrial DNA and 1549 aligned base pairs of nuclear DNA from 55 samples representing 48 ingroup and seven outgroup species to test this hypothesis. Mitochondrial genes analyzed include those encoding the 12S and 16S ribosomal RNAs; transfer RNAs coding for valine, leucine, isoleucine, glutamine, methionine, tryptophan, alanine, asparagine, cysteine and tyrosine; and complete cytochrome b and NADH dehydrogenase subunit I and II; nuclear gene analyzed included the third exon of the recombination activation gene 1 (RAG1). Analyses of combined mtDNA and nuclear DNA data sets result in a well-supported phylogenetic hypothesis. This hypothesis is in conflict with the classical taxonomic assignment of genera into tribes and phylogenetic hypotheses based on the taxonomy; however, the molecular hypothesis defines nine clades that are geographically restricted and consistent with the geological evolution of Central America and the Caribbean. Our analyses support multiple colonization events of Middle America followed by a mix of vicariance and dispersal events.     

Extensive trans-species mitochondrial polymorphisms in the carabid beetles Carabus subgenus Ohomopterus caused by repeated introgressive hybridization

To study the potential importance of introgressive hybridization to the evolutionary diversification of a carabid beetle lineage, we studied intraspecific and trans-species polymorphisms in the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequence (1083 bp) in four species of the subgenus Ohomopterus (genus Carabus) in central and eastern Honshu, Japan. Of the four species, C. insulicola is parapatric with the other three, and can hybridize naturally with at least two. This species possesses two haplotypes of remote lineages. We classified ND5 haplotypes using polymerase chain reaction-restriction fragment length polymorphism with TaqI endonuclease for 524 specimens, and sequenced 143 samples. Analysis revealed that each species was polyphyletic in its mitochondrial DNA phylogeny, representing a marked case of trans-species polymorphism. Recent one-way introgression of mitochondria from C. arrowianus nakamurai to C. insulicola, and from C. insulicola to C. esakii, was inferred from the frequency of identical sequences between these species and from direct evidence of hybridization in their contact zones. Other intraspecific polymorphisms in the four species may be due to undetected introgressive hybridization (e.g. C. insulicola to C. maiyasanus) or from stochastic lineage sorting of ancestral polymorphisms. This beetle group has a genital lock-and-key system, with species-specific or subspecies-specific genital morphology that may act as a barrier to hybridization. However, our results demonstrate that introgressive hybridization has occurred multiple times, at least for mitochondria, despite differences among, and stability within, morphological characters that distinguish local populations. Thus, hybridization and introgression could have been key processes in the evolutionary diversification of Ohomopterus.     

Assessing natural and disturbed population structure in European grayling Thymallus thymallus: melding phylogeographic,population genetic and jurisdictional perspectives for conservation planning

Natural and human‐mediated population structure of European grayling Thymallus thymallus was assessed in primarily Austrian reaches of the Danube River catchment. Data on phylogeographic structure based on mitochondrial DNA (mtDNA) were combined with variation found across 13 microsatellite loci to assess introgression stemming from stocking activities. Populations north of the Alps generally showed no signs of long‐term introgression from allochthonous stocks, although one population comprised two distinct genotypic groups, one of which may stem from recently stocked material. South of the Alps, levels of introgression with stocks from the north were extensive and only one of six populations is believed to represent a reference state. Allelic diversity and expected heterozygosity were positively correlated with mtDNA admixture supporting more recent introgression and that there is little or no reproductive barrier between these two major lineages north and south of the Alps. A third unrelated mtDNA lineage is described from the Lafnitz, a tributary of the Raab drainage. The high lineage diversity in the upper Austrian Danube is not concordant with an existing model of T. thymallus evolutionarily significant units in Europe. Caveats in naming such units or following a strict hierarchical conservation unit structure for broadly distributed species with complex phylogeographic distributions stretching over various jurisdictions are discussed. The necessity of using both phylogeographic and population genetic approaches in evaluating the history and conservation value of populations in a conservation context is additionally highlighted.     

Ecology and genetics of interspecific hybridization in the swallowtails, Papilio hospiton Gene and P. machaon L., in Corsica (Lepidoptera: Papiliomdae)

The Corsican swallowtail butterfly, Papilio hospiton , is endemic to Corsica and Sardinia (France and Italy) and included in the list of endangered species by the Washington Convention, It is spread all over Corsica in scattered populations linked to diverse habitats. A study by enzyme electrophoresis showed that the genetic diversity of the species is of the same order of magnitude as that of Papilio machaon from continental France and Corsica. The differentiation between populations is rather low, which is consistent with the high vagility of the adults. Natural hybridization between P. hospiton and P. machaon is frequent, and laboratory crosses show that the hybrids are not sterile. However, developmental perturbations impair the viability of further hybrid progenies. Although limited introgression between the two species is likely to take place, enzyme electrophoresis and PCR-RFLP analysis of mitochondrial DNA show that their gene pools remain distincfTtienetic assimilation by P. machaon therefore does not seem to be a threat for P. hospiton.     

Phylogeography and introgressive hybridization: chipmunks (genus Tamias) in the northern Rocky Mountains

Abstract. If phylogeographic studies are to be broadly used for assessing population-level processes relevant to speciation and systematics, the ability to identify and incorporate instances of hybridization into the analytical framework is essential. Here, we examine the evolutionary history of two chipmunk species, Tamias ruficaudus and Tamias amoenus , in the northern Rocky Mountains by integrating multivariate morphometrics of bacular (os penis) variation, phylogenetic estimation, and nested clade analysis with regional biogeography. Our results indicate multiple examples of mitochondrial DNA introgression layered within the evolutionary history of these nonsister species. Three of these events are most consistent with recent and/or ongoing asymmetric introgression of mitochondrial DNA across morphologically defined secondary contact zones. In addition, we find preliminary evidence where a fourth instance of nonconcordant characters may represent complete fixation of introgressed mitochondrial DNA via a more ancient hybridization event, although alternative explanations of convergence or incomplete sorting of ancestral polymorphisms cannot be dismissed with these data. The demonstration of hybridization among chipmunks with strongly differentiated bacular morphology contradicts long-standing assumptions that variation within this character is diagnostic of complete reproductive isolation within Tamias . Our results illustrate the utility of phylogeographic analyses for detecting instances of reticulate evolution and for incorporating this and other information in the inference of the evolutionary history of species.     

Rare and asymmetrical hybridization of the endemic Barbus carpathicus with its widespread congener Barbus barbus

In endemic species that co-occur with widespread congeners, hybridization can lead to an influx of novel and beneficial genetic variation, but high rates of introgression may cause genetic swamping of the endemic species and have detrimental effects on its survival potential. This study examines hybridization between sympatric populations of the Carpathian barbel Barbus carpathicus , a recently discovered cryptic species with a restricted range, and the widespread common barbel Barbus barbus . Based on six diagnostic allozyme loci, a microsatellite locus and mtDNA, hybrids were found to be present at multiple localities within the Vistula River drainage (Baltic Sea) as well as in the Tisza River system of the Danube River drainage (Black Sea). However, the numbers of hybrids were very low; four individuals of 230 fish sampled from the Vistula drainage. Bayesian assessment of their nuclear genotypes suggested that two hybrids in the Vistula drainage and nine in the Tisza system were F₁ generation, and one in the Vistula drainage and one in the Tisza system were backcrosses (BC) to B. barbus , while no F₂ or BC to B. carpathicus were detected. No hybrid carried B. carpathicus mtDNA and cytonuclear linkage disequilibria showed significant positive associations between hybrid genotypes and B. barbus mtDNA, suggesting unidirectionality in the interspecific mating with a disproportionate contribution of B. barbus mothers. Despite geographically broad occurrence of hybrids, these data provide evidence of strong constraints on hybridization in the native breeding habitats and the lack of introgression towards B. carpathicus .     

Evidence of natural hybridization among homothallic members of the basidiomycete Armillaria mellea sensu stricto

Populations of Armillaria mellea (Basidiomycota, Agaricales) across much of its range are heterothallic; homothallic populations occur only in Africa (A. mellea ssp. africana), China (China Biological Species CBS G), and Japan (A. mellea ssp. nipponica). Monosporous isolates of heterothallic A. mellea are haploid and their mating behaviour is consistent with the requirement of two different alleles at two mating-type loci (tetrapolar mating system) to create a diploid individual. In contrast, monosporous isolates of homothallic A. mellea are putatively diploid; they bypass the haploid phase by undergoing karyogamy in the basidium (a unique type of secondary homothallism/pseudohomothallism). In order to determine the genetic origin of this homothallism, we analyzed genetic variation of 47 heterothallic isolates from China, Europe, and North America, and 14 homothallic isolates from Africa, China, and Japan. Gene trees and mutational networks were constructed for partial mitochondrial gene ATP synthase subunit 6 (ATP6) and for the following nuclear genes: actin (ACTIN), elongation factor subunit 1-alpha (EFA), glyceraldehyde 3-phosphate dehydrogenase (GPD), and the RNA polymerase subunit II (RPB2). Homothallic isolates from Africa and Japan shared a common mitochondrial ATP6 haplotype with homothallic isolates from China, and are likely introductions. Homothallic isolates from China that shared a common mitochondrial haplotype with all European isolates did not share European nuclear haplotypes, as revealed by median-joining networks, but instead clustered with haplotypes from China or were intermediate between those of China and Europe. Such mitochondrial-nuclear discordance in homothallic isolates from China is indicative of hybridization between lineages originating from China and Europe.     

THE EXTENT OF INTROGRESSION OUTSIDE THE CONTACT ZONE BETWEEN NOTROPIS CORNUTUS AND NOTROPIS CHRYSOCEPHALUS (TELEOSTEI: CYPRINIDAE)

The cyprinid fishes, Notropis cornutus and N. chrysocephalus, hybridize in a long, narrow zone in the midwestern United States. To quantify the extent of introgression of genetic markers outside of this zone, samples were collected along transects starting near the region of contact (as defined by morphological characters), followed by samples progressively more distant. Diagnostic allozymic and mitochondrial DNA (mtDNA) restriction site markers were used to estimate the extent of introgression outside of the zone, while polymorphic allozyme and mtDNA markers were used to evaluate the potential for gene flow among populations within transects. Analysis of populations from the northern transect provided evidence for differentiation of populations for some of the markers; however, on average, enough gene flow has occurred to overcome substantial differentiation. Introgressed mtDNA and allozyme haplotypes were rare and found only in the population closest to the contact zone. The rarity of introgressed alleles in the more northern populations is consistent with the recent origin of these populations after the Wisconsin glaciation (less than 12,000 years bp) and/or selection maintaining the northern boundary of the contact zone. Analysis of populations from the southern transect revealed evidence for population subdivision but no evidence for introgression at the diagnostic allozyme loci; however, nearly all individuals from this transect possessed introgressed mtDNA haplotypes, with samples furthest from the contact zone exhibiting the highest frequencies of introgression. Patterns of variation for one of the polymorphic allozyme markers (Est-A) and introgressed mtDNAs were highly correlated, suggesting that allozymic heterogeneity at this locus is also the result of introgression. The most likely explanation for these data is that these introgressed haplotypes are indicators of a more southern position of the contact zone during the Pleistocene, with the contact zone shifting northward with the recession of the glacial front. Such movement implicates selection in the maintenance of distributional limits of these species, and hence, the width and position of the contact zone.     

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.