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.     

Genetic differentiation among European whitefish ecotypes based on microsatellite data

The amount of genetic differentiation at DNA microsatellite loci in European whitefish (Coregonus lavaretus) was assessed among ecotypes, populations and run-timing types. The magnitude of genetic changes potentially caused by hatchery broodstock rearing were also compared with those observed in corresponding natural populations. A total of 35 populations were studied, including 33 Coregonus lavaretus populations and two samples of Coregonus peled. Five of the six whitefish ecotypes in Finland were represented within C. lavaretus populations. Genetic diversity among C. lavaretus populations proved to be high compared to two C. peled populations. The genetic D(A) distance between these two species was as high as 0.86. The genetic differentiation among ecotypes was generally low and thus gives support for the hypothesis of one native European whitefish species in Scandinavia. Among the ecotypes the northern, large sparsely-rakered, bottom-dwelling whitefish was most unique. Thus, observed genetic differences in quantitative traits have either developed independently of phylogenetic lineages, or have mixed and later changed according to environments and selection pressures. Overall genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Populations of this area have been heavily influenced by human activities, and they also have the highest probability of mixing by natural means. In two cases, the Rivers Iijoki and Tornionjoki, statistically significant genetic differences could be observed between summer- and autumn-run spawning-time types. Wild populations had slightly higher allelic diversity than hatchery-reared populations of corresponding rivers. Although some reduction in genetic diversity during hatchery rearing is possible, it is an important aid in maintaining endangered populations.     

Postglacial recolonization patterns and genetic relationships among whitefish (Coregonus sp.) populations in Denmark, inferred from mitochondrial DNA and microsatellite markers

The genetic relationships among morphologically and geographically divergent populations of whitefish (genus: Coregonus ) from Denmark and the Baltic Sea region were studied by analysis of microsatellites and polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis of mitochondrial DNA (mtDNA) segments. The endangered North Sea houting (classified as C. oxyrhynchus ) differs morphologically and physiologically from other Danish whitefish ( C. lavaretus ). However, limited divergence of North Sea houting was observed both at the level of mtDNA and microsatellites. The implications of these results for the conservation status of North Sea houting are discussed in the light of current definitions of evolutionary significant units. Both mtDNA and microsatellite data indicated that postglacial recolonization by C. lavaretus in Denmark was less likely to have taken place from the Baltic Sea. Instead, the data suggested a recent common origin of all Danish whitefish populations, including North Sea houting, probably by recolonization via the postglacial Elbe River system. Estimates of genetic differentiation among populations based on mtDNA and microsatellites were qualitatively different. In addition, for both classes of markers analyses of genetic differentiation yielded different results, depending on whether molecular distances between alleles or haplotypes were included.     

ALLOPATRIC ORIGIN OF SYMPATRIC POPULATIONS OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS) AS REVEALED BY MITOCHONDRIAL-DNA RESTRICTION ANALYSIS

In the paper, restriction-fragment length polymorphisms in mitochondrial DNA (mtDNA) were studied to test the hypothesis that sympatric populations of lake whitefish in the Allegash basin have recently diverged through sympatric speciation. Thirteen restriction enzymes were used to analyze mtDNA of 156 specimens representing 13 populations from eastern Canada and northern Maine where normal and dwarf phenotypes of whitefish exist in sympatry and allopatry. Two monophyletic assemblages of populations that exhibit different geographic distributions were identified. One showed an eastern distribution that expands from Cape Breton to the Allegash basin and the other exhibits a more western distribution. The Allegash basin was the only area of overlap. The western assemblage exhibited the normal size phenotype in all cases, whereas the eastern assemblage exhibited the normal size phenotype in allopatric conditions and the dwarf size phenotype in sympatry. The existence of sympatric pairs in the Allegash basin result from the secondary contact of two monophyletic groups of whitefish that evolved allopatrically in separate refugia during the last glaciation events. The weak mtDNA difference of sympatric pairs suggests that speciation of lake whitefish in eastern North America was accompanied by only minor alterations of the ancestral gene pool.     

Mitochondrial DNA analysis distinguishes between British populations of the whitefish

The whitefish Coregonus lavaretus is confined to two Scottish lochs and four English and one Welsh lake in the British Isles. Restriction fragment length polymorphisms (RFLPs) in the mitochondrial (mt) DNA have been used to analyse variability within and between populations. The 14 restriction enzymes employed in this study examined approximately 3% of the mt genome and revealed 15 haplotypes among 156 fish. Haplotype diversity within populations was generally low except for the Welsh population where nine haplotypes were found among 30 fish examined. The haplotypes clustered into three distinct groups corresponding to Scottish, English and Welsh populations. No haplotypes were shared among the three groups. Possible alternative explanations for these findings are introgression from another species, stochastic lineage sorting from a polymorphic ancestral gene pool and/or separate colonizations following the last glaciation.     

PHYLOGEOGRAPHIC STRUCTURE IN MITOCHONDRIAL DNA OF THE LAKE WHITEFISH (COREGONUS CLUPEAFORMIS) AND ITS RELATION TO PLEISTOCENE GLACIATIONS

Restriction-fragment length polymorphisms were employed to evaluate the phylogenetic relationships, the genetic diversity and the geographic structure in mitochondrial DNA (mtDNA) lineages of the lake whitefish, Coregonus clupeaformis. Thirteen restriction enzymes that produced 148 restriction fragments were used to assay mtDNAs of 525 specimens collected among 41 populations. The sampling covered the entire range of the species, from Alaska to Labrador. Four distinct phylogeographic assemblages were identified. The Beringian assemblage, confined to Yukon and Alaska, was phylogenetically distinct from other assemblages and exhibited the highest level of nucleotide diversity. The Acadian assemblage was confined to southeastern North America and composed of a unique mtDNA clade. The Atlantic assemblage was confined to southern Québec and the northeastern United States and was also observed among anadromous populations of northern Hudson Bay. This group was highly polymorphic and responsible for most of the mtDNA diversity observed outside Beringia. The Mississippian assemblage occupied most of the actual range of lake whitefish, from the Mackenzie delta to Labrador. Ninety-two percent of all whitefish of this proposed origin belonged to a single mtDNA haplotype. Overall, the diversity, the geographic structure and the times of divergence of mtDNA phylogenetic assemblages correlate with the Pleistocene glaciations classically assumed to have dramatically altered the genetic diversity of northern fishes in recent evolutionary times. Our results emphasize the dominant role of these catastrophic events in shaping the population genetic structure of lake whitefish.     

Morphology and mitochondrial DNA variation of the Siberian whitefish <Emphasis Type="Italic">Coregonus lavaretus pidschian</Emphasis> (Gmelin) in the upstream water bodies of the Ob and Yenisei Rivers

Combining morphological, ecological and genetic analyses, we compared patterns of diversification within and among populations of the southern Siberian whitefish species Coregonus lavaretus pidschian (Gmelin) to illuminate their evolutionary history. Using sequencing data from 1,930 bp of NADH dehydrogenase subunit 1 (ND1) mitochondrial DNA regions, we documented phylogeographic relationships among populations and developed a phylogeny of mtDNA haplotypes. We found significant differences in the perforated lateral-line scale numbers within and between some populations. Clear differences in the number of gill rakers on the first branchial arch were only exhibited between populations of C. l. pidschian and Coregonus lavaretus pravdinellus Dulkeit. Concordance between different morphological groups based on two meristic traits and mtDNA patterns was also tested.     

GENETIC EVIDENCE FOR REPRODUCTIVE ISOLATION AND MULTIPLE ORIGINS OF SYMPATRIC TROPHIC ECOTYPES OF WHITEFISH (COREGONUS)

We assessed variation in mitochondrial DNA (mtDNA) by restriction fragment length polymorphism (RFLP) analysis and in nuclear genes by allozyme analysis among sympatric pairs of limnetic and benthic ecotypes of whitefish (Coregonus) coexisting in three lakes of southern Yukon to address three evolutionary questions regarding their origins. Are sympatric low and high gill-raker count ecotypes genetically differentiated? Are they issued from monophyletic or polyphyletic evolutionary events? If they are polyphyletic in origins, did they originate from multiple allopatric speciation events or intralacustrine radiation? Our results corroborated previous genetic and ecological studies of these ecotypes, indicating that they represent genetically distinct reproductive units, and therefore refuting the hypothesis of phenotypic polymorphism within a single population. However, the amount of gene flow between ecotypes varied among lakes, correlating with the extent of morphological differentiation and the potential for premating reproductive isolation. The results indicated a polyphyletic origin of ecotypes whereby each of them have been expressed independently more than once. In the two lakes of Squanga Creek drainage, the existence of sympatric pairs was best explained by the secondary contact of two monophyletic whitefish groups that evolved in allopatry during the last glaciation events. In Dezadeash L. of Alsek R. drainage, our results could not verify either sympatric or allopatric (or microallopatric) origin of ecotypes. Regardless of the mode of speciation involved in their origins, these sympatric whitefish populations provided further evidence that Pleistocene glaciation events created conditions favoring rapid divergence and phenotypic differentiation among northern freshwater fishes.     

Different growth trends of whitefish (Coregonus lavaretus) forms in the northern Baltic Sea

The sea growth of two whitefish forms, anadromous (Coregonus lavaretus lavaretus) and sea‐spawning (Coregonus lavaretus widegreni), was analysed using samples collected from the commercial sea catch in the Gulf of Bothnia (GoB) in the northern Baltic Sea during 1998–2014. In the GoB area, these two forms are possible to identify because the gill‐raker number and size at maturity vary between forms. The growth rate of the forms is linked to their feeding area. Sea‐spawning whitefish, which has a feeding migration near its home site, was shorter in the northern GoB (66°N–64°N) at the ages of 3–11 than those in the southern GoB (64°N–60°30′N). In the data, most whitefish were caught with gill nets in the GoB. The mesh sizes of gill nets capturing the anadromous form were mostly 35–45 mm, while those capturing the sea‐spawning form were <35 mm in the northern GoB. It is likely that the different growth trends for small and large whitefish were connected with differences in their recruitment for fishing. The length of anadromous females at the age of four sea years increased significantly, but the length of six‐year‐old anadromous female whitefish decreased over the catch years from 1998–2014. In contrast, the length of slow‐growing sea‐spawning whitefish of six years or older increased significantly in relation to the catch year in the gill‐net catch. The increase in the growth of young age groups in both forms was probably associated with the increasing temperature and the low fishing pressure on small fish. The decreasing age at capture for both forms and the depression of the mean size of old anadromous whitefish are signs of high fishing pressure with a high gill‐net effort that selectively removes the largest and oldest individuals of both forms.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

MULTIPLE MODES OF SPECIATION INVOLVED IN THE PARALLEL EVOLUTION OF SYMPATRIC MORPHOTYPES OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS,SALMONIDAE)

We performed a phylogenetic analysis of mtDNA variation among seven sympatric pairs of dwarf and normal morphotypes of whitefish from northern Québec and the St. John River drainage to address three questions relevant to understanding their radiation. Are all sympatric pairs reproductively isolated? Do phylogenetic analyses confirm that sympatric whitefish morphotypes found in eastern North America represent the outcome of polyphyletic evolutionary events? If so, did all sympatric pairs from the St. John River drainage originate from the same scenario of allopatric divergence and secondary contact? The hypothesis of genetic differentiation was supported for all sympatric pairs from the St. John River drainage, whereas lack of mtDNA diversity precluded any test of reproductive isolation for northern Québec populations. Patterns of mtDNA variation confirmed that dwarf and normal morphotypes evolved in parallel among independent, yet closely related, lineages, thus providing indirect evidence for the role of natural selection in promoting phenotypic radiation in whitefish. Patterns of mtDNA diversity among sympatric pairs of the St. John River indicated a complex picture of whitefish evolution that implied sympatric divergence and multiple allopatric divergence/secondary contact events on a small geographic scale. These results suggests that ecological opportunities, namely trophic niche availability, may promote population divergence in whitefish.     

Late glacial history of the cold-adapted freshwater fish Cottus gobio,revealed by microsatellites

The distribution of genetic diversity at 10 highly polymorphic microsatellite loci within the European freshwater fish, Cottus gobio, L. was examined. The sampling range comprised a large geographical scale including lineages known to be highly divergent at both mitochondrial DNA (mtDNA) and allozymes. An analysis of genetic variability within populations showed that expected heterozygosity and allelic richness could be explained largely by current effective population sizes. Evidence was found, however, that historical processes predating the last major glaciation affected allelic richness. In addition to confirming the large-scale patterns from earlier studies, the microsatellite data revealed new insights into recent processes by analysing genetic structure within ancient lineages defined by mtDNA data. Stepwise mutation model (SMM) and nonSMM-based methods demonstrated a clear genetic structuring within the Northwestern European lineage comprising populations from Britain and Belgium, and within the Central European lineage populations from the rivers Danube, Elbe and Main. Supported by an analysis of genetic variability within populations these results showed that the bullhead populations most probably persisted throughout the last major glaciation within the British Isles and within the drainages of the rivers Elbe and Main. Such observations provide the first genetic evidence for a glacial refugium in such close proximity to the European glacial margins.     

Phylogeography of Varicorhinus barbatulus (Cyprinidae) in Taiwan based on nucleotide variation of mtDNA and allozymes

The phylogeographical patterns and population genetic structures of Varicorhinus barbatulus in Taiwan were investigated based on genetic diversity of 34 allozyme loci and nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mtDNA. Allozyme and mtDNA analyses revealing evident geographical structuring suggest limited gene flow between populations (F(ST)=0.511 and 0.791, respectively). Low genetic variability within populations (P=5.56%; He=0.018) based on allozymes and significantly negative Tajima's D statistics based on mtDNA suggest that most populations in Taiwan may have originated from a small number of founders followed by demographic expansion. The gene genealogy of mtDNA identified six lineages corresponding to major drainages that were separated by the geological barriers due to vicariant events. A minimum spanning network based on nucleotide substitutions reflects divergence from populations of the Miao-li Plateau to northern and southern regions of the island. In contrast to a previous hypothesis that suggests an early invasion to eastern part of Taiwan prior to the lifting of central mountain range some one million years ago, the mtDNA genealogy and molecular dating reveal very recent colonization of the eastern population. Nested clade analyses revealing significant associations between genetic structure and geographical division identify past fragmentation and range expansion as major phylogeographical events that shaped the geographical distribution of this species in Taiwan.     

Contemporary factors influencing genetic diversity in the Alaska humpback whitefish Coregonus clupeaformis complex

Thirteen microsatellite loci were used to address three hypotheses regarding genetic diversity in the humpback whitefish Coregonus clupeaformis complex in Alaska. The test results provided further insight into the factors influencing C. clupeaformis complex population structure and level of genetic variation. First, the microsatellite data did not provide evidence of two spatially distinct Beringian and Eurasian refugial groups as revealed in previous phylogeographic analyses of mitochondrial DNA variation. Rather, the population structure inferred from the microsatellite variation appears to reveal the influence of factors on a more recent time scale, including gene flow among the refugial groups and isolation of some anadromous and freshwater resident populations. Second, anadromous C. clupeaformis complex collections exhibited higher intra‐population genetic diversity than freshwater resident collections. This outcome is consistent with previous meta analyses suggesting that freshwater resident populations probably have smaller historical effective population sizes and less conspecific gene flow because the habitat tends to be smaller and supports fewer and smaller populations. Finally, the analysis of contemporary immigration rates was consistent with, but did not provide statistical support for, the hypothesis that gene flow among anadromous C. clupeaformis complex populations along coastal Alaska is influenced by the Alaska Coastal Current. Further studies are needed to evaluate gene flow among coastal Alaska C. clupeaformis complex populations.     

Gene flow and effective population size in two life-history types of broad whitefish Coregonus nasus from the Canadian Arctic

In this study, the magnitude and direction of gene flow and estimates of effective population sizes (N(e) ) were quantified among two life-history types (lacustrine and anadromous) of broad whitefish Coregonus nasus in the lower Mackenzie River system. The data suggest that dispersal and subsequent gene flow occurs between these groups, with the former appearing to be asymmetrical. Gene flow may potentially be directionally biased as well, a result attributed to source-sink population dynamics and the ongoing process of post-glacial colonization and contemporary range expansion. Additionally, average N(e) estimates were consistently lower for lacustrine populations of C. nasus although confidence intervals for both contemporary and historical estimates broadly overlapped. The lower average estimates of N(e) for lacustrine populations was suggested to be the result of more recent founding events following post-glacial dispersal. This study provides one of the first assessments of gene flow and N(e) in an Arctic coregonine, results that may be relevant to other freshwater and anadromous Arctic species persisting in systems near the periphery of their range.     

AFLP utility for population assignment studies: analytical investigation and empirical comparison with microsatellites

Individual-based population assignment tests have thus far mainly relied on the use of microsatellite loci. However, the logistic difficulty of screening large numbers of loci required to reach sufficient statistical power hampers the usefulness of microsatellites in situations of weak population structuring. Amplified fragment length polymorphisms (AFLP) represents an alternative for overcoming this logistical issue as the technique allows the user to characterize a much larger number of loci with a comparable analytical effort. In this study, an assignment test based on maximum likelihood for dominant markers was used to investigate the potential usefulness of AFLP for population assignment. We also compared assignment success achieved with AFLP with that obtained using microsatellites in a case study of low population differentiation involving whitefish (Coregonus clupeaformis) sympatric ecotypes. The analytical investigation showed that the minimum number of AFLP loci required to reach an assignment success of 95% stood within values that are easily achievable in many situations. This also showed how assignment success varied according to the number of AFLP loci used, their absolute frequency and their frequency differential and sampling errors, as well as the number of putative source populations. The case study showed that given a comparable analytical effort in the laboratory, AFLP were much more efficient than the microsatellite loci in discriminating the source of an individual among putative populations. AFLP resulted in higher assignment success at all levels of stringency and the log-likelihood differences between populations obtained with AFLP for each individual were much larger than those obtained with microsatellites. These results indicate that research involving individual-based population assignment methods should benefit importantly from the use of AFLP markers, especially in systems characterized by weak population structuring.     

Identification of palearctic coregonid fish species using mtDNA and allozyme genetic markers

Two types of molecular genetic markers were used for genetic identification of species and local stocks of palearctic coregonids (Coregonidae, Salmoniformes, Teleostei). Seven nominate species of whitefishes and ciscoes Coregonus , spp. of Eurasia Arctic Sea basin and inconnu Stenodus leucichthys nelma , represented by specimens from North America were studied. Using restriction analysis of PCR-amplified products of the ND-1 gene of mitochondrial DNA (mtDNA) and allelic composition at several allozyme loci discrimination was successful between C. lavaretus pidschian , Siberian whitefish, C. nasus , broad whitefish, C. autumnalis , Arctic cisco, C. migratorius , Baikal omul, C. peled , peled, and C. sardinella , least cisco. Muksun C. muksun , was indistinguishable from Siberian whitefish. Creatine kinase (CK) isozyme patterns and Rsa , I restriction patterns of ND-1 were the most effective markers allowing discrimination among species. Intra-specific differentiation in mtDNA was found in all species but was much less pronounced than inter-species variation. In several specimens composite haplotypes typical of another species were found that reflect probable gene introgression by hybridization. A combination of mtDNA and nuclear genetic markers is suggested for reliable identification of both typical species representatives and hybrids.     

How normal is abnormal? Discrimination between deformations and natural variation in gonad morphology of European whitefish Coregonus lavaretus

The gonad morphology of whitefish Coregonus lavaretus collected in Lake Thun, Switzerland, and two neighbouring lakes was assessed in order to differentiate between 'normal' and 'abnormal' character states of gonad morphology, which had been previously described in C. lavaretus from Lake Thun (constrictions, asymmetries, aplasia, compartmentations, fusions and hermaphroditism). In total, 4668 fish were collected and analysed using two complementary sampling schemes: (1) monthly samples of catches by the commercial fishermen and (2) samples of ripe spawners of all known 33 spawning sites of the three lakes. Considerable variation in gonad morphology in C. lavaretus populations of all lakes was found. Notably, all deviation types were observed in fish of all three lakes. Asymmetries and constrictions were frequent in all three lakes and showed systematic differences in frequency between the two sampling strategies. This indicates that asymmetries and constrictions represent to a large extent natural variation in gonad morphology of C. lavaretus and are also prone to considerable measurement error. In contrast, aplasia, fusions, compartmentations and hermaphroditism occurred predominantly in one C. lavaretus form of Lake Thun and in particular in populations spawning at great depths. This suggests that these deviation types are probably reliable indicators for gonad deformations and supports the interpretation that Lake Thun harbours a unique case of deformed gonads in C. lavaretus of yet unknown origin.     

Species introduction promotes hybridization and introgression in Coregonus: is there sign of selection against hybrids?

Species introductions are considered one of the major drivers of biodiversity loss via ecological interactions and genetic admixture with local fauna. We examined two well-recognized fish species, native whitefish (Coregonus lavaretus) and introduced vendace (Coregonus albula), as well as their morphological hybrids in a single lake to test for selection against hybrids and backcrosses in the wild. A representative random subsample of 693 individuals (27.8%) was taken from the total catch of coregonids. This subsample was examined with the aim to select c. 50 individuals of pure whitefish (n = 52), pure vendace (n = 55) and putative hybrid (n = 19) for genetic analyses. The subsequent microsatellites and mitochondrial (mt) DNA analyses provided compelling evidence of hybridization and introgression. Of the 126 fish examined, four were found to be F(1) , 14 backcrosses to whitefish and seven backcrosses to vendace. The estimates of historical gene flow suggested higher rates from introduced vendace into native whitefish than vice versa, whereas estimates of contemporary gene flow were equal. Mitochondrial introgression was skewed, with 18 backcrosses having vendace mtDNA and only three with whitefish mtDNA. Hybrids and backcrosses had intermediate morphology and niche utilization compared with parental species. No evidence of selection against hybrids or backcrosses was apparent, as both hybrid and backcross growth rates and fecundities were high. Hybrids (F(1) ) were only detected in 2 year-classes, suggesting temporal variability in mating between vendace and whitefish. However, our data show that hybrids reached sexual maturity and reproduced actively, with backcrosses recorded from six consecutive year-classes, whereas no F(2) individuals were found. The results indicate widespread introgression, as 10.8% of coregonids were estimated to be backcrosses.