Microsatellites in grayling (Thymallus thymallus): comparison of two geographically remote populations from the Danubian and Adriatic river basin in Slovenia

Neotype designation for Thymallus aeliani Valenciennes, 1848 from a museum topotype specimen and its affiliation with Adriatic grayling on the basis of mitochondrial DNA

In 1848, the grayling Thymallus aeliani (Valenciennes) was described from Lake Maggiore, Italy, in the north Adriatic basin. Genetic analyses of the mitochondrial control region showed a unique evolutionary history of grayling inhabiting the rivers of northern Adriatic basin, from the upper reaches of the Po River and its left tributaries in the west to the Soča River in the east, which resulted in the designation of this phylogenetic lineage as Adriatic grayling. Consequently, the name T. aeliani was connected to the Adriatic lineage, re-establishing the validity of this taxon. However, the mitochondrial haplotypes belonging to Adriatic grayling were never compared with the type specimens of T. aeliani, as their whereabouts were unknown. In this study, a neotype for T. aeliani was designated using topotypical specimens stored at the Natural History Museum in Vienna. The neotype (NMW 68027:2 labelled as "Lago Maggiore, Bellotti, 1880") was designated pursuant to the conditions stipulated in Article 75.3 of the International Code of Zoological Nomenclature. Furthermore, the mitochondrial control region of the neotype was compared to haplotypes of the Adriatic lineage and showed high genetic similarity, which therefore connects the species name T. aeliani to the Adriatic grayling. This crucial step in fixing nomenclatural status of this species is very important for its protection and management.

Geographic hierarchical population genetic structuring in British European whitefish (Coregonus lavaretus) and its implications for conservation

The European whitefish Coregonus lavaretus complex represents one of the most diverse radiations within salmonids, with extreme morphological and genetic differentiation across its range. Such variation has led to the assignment of many populations to separate species. In Great Britain, the seven native populations of C. lavaretus (two in Scotland, four in England, one in Wales) were previously classified into three species, and recent taxonomic revision resurrected the previous nomenclature. Here we used a dataset of 15 microsatellites to: (1) investigate the genetic diversity of British populations, (2) assess the level of population structure and the relationships between British populations. Genetic diversity was highest in Welsh (HO = 0.50, AR = 5.29), intermediate in English (HO = 0.41–0.50, AR = 2.83–3.88), and lowest in Scottish populations (HO = 0.28–0.35, AR = 2.56–3.04). Population structure analyses indicated high genetic differentiation (global FST = 0.388) between all populations but for the two Scottish populations (FST = 0.063) and two English populations (FST = 0.038). Principal component analysis and molecular ANOVA revealed separation between Scottish, English, and Welsh populations, with the Scottish populations being the most diverged. We argue that the data presented here are not sufficient to support a separation of the British European whitefish populations into three separate species, but support the delineation of different ESUs for these populations.

Sex differentiation in grayling (Salmonidae) goes through an all-male stage and is delayed in genetic males who instead grow faster

Fish populations can be threatened by distorted sex ratios that arise during sex differentiation. Here we describe sex differentiation in a wild grayling (Thymallus thymallus) population that suffers from distorted sex ratios. We verified that sex determination is linked to the sex determining locus (sdY) of salmonids. This allowed us to study sex-specific gene expression and gonadal development. Sex-specific gene expression could be observed during embryogenesis and was strong around hatching. About half of the fish showed immature testes around eleven weeks after fertilization. This phenotype was mostly replaced by the "testis-to-ovary" or "ovaries" phenotypes during development. The gonads of the remaining fish stayed undifferentiated until six months after fertilization. Genetic sexing revealed that fish with undifferentiated gonads were all males, who grew larger than the genetic females during the observational period. Only 12% of the genetic males showed testicular tissue six months after fertilization. We conclude that sex differentiation starts before hatching, goes through an all-male stage for both sexes (which represents a rare case of "undifferentiated" gonochoristic species that usually go through an all-female stage), and is delayed in males. During these juvenile stages males grow faster than females instead of developing their gonads.

Inter and intra-population phenotypic and genotypic structuring in the European whitefish Coregonus lavaretus, a rare freshwater fish in Scotland

This study revealed between-lake genetic structuring between Coregonus lavaretus collected from the only two native populations of this species in Scotland, U.K. (Lochs Eck and Lomond) evidenced by the existence of private alleles (12 in Lomond and four in Eck) and significant genetic differentiation (FST = 0·056) across 10 microsatellite markers. Juvenile C. lavaretus originating from eggs collected from the two lakes and reared in a common-garden experiment showed clear phenotypic differences in trophic morphology (i.e. head and body shape) between these populations indicating that these characteristics were, at least partly, inherited. Microsatellite analysis of adults collected from different geographic regions within Loch Lomond revealed detectable and statistically significant but relatively weak genetic structuring (FST = 0·001-0·024) and evidence of private alleles related to the basin structure of the lake. Within-lake genetic divergence patterns suggest three possibilities for this observed pattern: (1) differential selection pressures causing divergence into separate gene pools, (2) a collapse of two formerly divergent gene pools and (3) a stable state maintained by balancing selection forces resulting from spatial variation in selection and lake heterogeneity. Small estimates of effective population sizes for the populations in both lakes suggest that the capacity of both populations to adapt to future environmental change may be limited.

Genetic parentage analysis confirms a polygynandrous breeding system in the European grayling (Thymallus thymallus)

Knowing the breeding system of a species is important in order to understand individual variation in reproductive success. Large variation in reproductive success and thus reproductive skew strongly impacts on the effective number of breeders and thus the long-term effective population size (Ne). Fishes, in particular species belonging to the salmonid family, exhibit a wide diversity of breeding systems. In general, however, breeding systems are rarely studied in detail in the wild. Here we examine the breeding system of the spring-spawning European grayling Thymallus thymallus from a small Norwegian stream using parentage assignment based on the genotyping of 19 polymorphic microsatellite loci. In total 895 individual grayling fry and 154 mature grayling (57 females and 97 males) were genotyped. A total of 466 offspring were assigned a father, a mother, or a parent pair with a confidence of 90% or higher. Successfully reproducing males had on average 11.9 ± 13.3 (SD) offspring with on average 2.1 ± 1.2 partners, whereas successful females had on average 9.5 ± 12.8 offspring and 2.3 ± 1.5 partners. Parents with more partners also produced more offspring. Thus the grayling breeding system within this small stream revealed a polygynandrous breeding system, similar to what has been observed for many other salmonid fish species. The present study thus unambiguously corroborates a polygynadrous breeding system in the European grayling. This knowledge is critical for managing populations of this species, which has suffered significant local population declines throughout its range over the last several decades.

Defining conservation units in a stocking-induced genetic melting pot: unraveling native and multiple exotic genetic imprints of recent and historical secondary contact in Adriatic grayling

The definition of conservation units is crucial for the sustainable management of endangered species, though particularly challenging when recent and past anthropogenic and natural gene flow might have played a role. The conservation of the European grayling, Thymallus thymallus, is particularly complex in its southern distribution area, where the Adriatic evolutionary lineage is endangered by a long history of anthropogenic disturbance, intensive stocking and potentially widespread genetic introgression. We provide mtDNA sequence and microsatellite data of 683 grayling from 30 sites of Adriatic as well as Danubian and Atlantic origin. We apply Bayesian clustering and Approximate Bayesian Computation (ABC) to detect microgeographic population structure and to infer the demographic history of the Adriatic populations, to define appropriate conservation units. Varying frequencies of indigenous genetic signatures of the Adriatic grayling were revealed, spanning from marginal genetic introgression to the collapse of native gene pools. Genetic introgression involved multiple exotic source populations of Danubian and Atlantic origin, thus evidencing the negative impact of few decades of stocking. Within the Adige River system, a contact zone of western Adriatic and eastern Danubian populations was detected, with ABC analyses suggesting a historical anthropogenic origin of eastern Adige populations, most likely founded by medieval translocations. Substantial river-specific population substructure within the Adriatic grayling Evolutionary Significant Unit points to the definition of different conservation units. We finally propose a catalog of management measures, including the legal prohibition of stocking exotic grayling and the use of molecular markers in supportive- and captive-breeding programs.

Speciation reversal in European whitefish (Coregonus lavaretus (L.)) caused by competitor invasion

Invasion of exotic species has caused the loss of biodiversity and imparts evolutionary and ecological changes in the introduced systems. In northern Fennoscandia, European whitefish (Coregonus lavaretus (L.)) is a highly polymorphic species displaying adaptive radiations into partially reproductively isolated and thus genetically differentiated sympatric morphs utilizing the planktivorous and benthivorous food niche in many lakes. In 1993, Lake Skrukkebukta was invaded by vendace (Coregonus albula (L.)) which is a zooplanktivorous specialist. The vendace displaced the densely rakered whitefish from its preferred pelagic niche to the benthic habitat harbouring the large sparsely rakered whitefish. In this study, we investigate the potential influence of the vendace invasion on the breakdown of reproductive isolation between the two whitefish morphs. We inferred the genotypic and phenotypic differentiation between the two morphs collected at the arrival (1993) and 15 years after (2008) the vendace invasion using 16 microsatellite loci and gill raker numbers, the most distinctive adaptive phenotypic trait between them. The comparison of gill raker number distributions revealed two modes growing closer over 15 years following the invasion. Bayesian analyses of genotypes revealed that the two genetically distinct whitefish morphs that existed in 1993 had collapsed into a single population in 2008. The decline in association between the gill raker numbers and admixture values over 15 years corroborates the findings from the Bayesian analysis. Our study thus suggests an apparent decrease of reproductive isolation in a morph-pair of European whitefish within 15 years (≃ 3 generations) following the invasion of a superior trophic competitor (vendace) in a subarctic lake, reflecting a situation of "speciation in reverse".

Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral, pelagic, and profundal lake habitats

Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatric European whitefish morphs from three northern Fennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.

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.

Genetic differentiation of European grayling (Thymallus thymallus) populations in Serbia, based on mitochondrial and nuclear DNA analyses

Background

The structure and diversity of grayling (Thymallus thymallus) populations have been well studied in most of its native habitat; however the southernmost populations of the Balkan Peninsula remain largely unexplored. The purpose of this study was to assess the genetic diversity of Serbian grayling populations, detect the impact of stocking and provide guidelines for conservation and management.

Methods

Eighty grayling individuals were collected from four rivers (Ibar, Lim, Drina and Rzav). The mitochondrial DNA control region (CR; 595 bp of the 3'end and 74 bp of flanking tRNA) and the ATP6 gene (630 bp fragment) were sequenced for 20 individuals (five from each locality). In addition, all individuals were genotyped with 12 microsatellite loci. The diversity and structure of the populations as well as the recent and ancient population declines were studied using specialized software.

Results

We detected three new haplotypes in the mtDNA CR and four haplotypes in the ATP6 gene of which three had not been described before. Previously, one CR haplotype and two ATP6 gene haplotypes had been identified as allochthonous, originating from Slovenia. Reconstruction of phylogenetic relations placed the remaining two CR haplotypes from the River Danube drainage of Serbia into a new clade, which is related to the previously described sister Slovenian clade. These two clades form a new Balkan clade. Microsatellite marker analysis showed that all four populations are genetically distinct from each other without any sign of intra-population structure, although stocking of the most diverse population (Drina River) was confirmed by mtDNA analysis. Recent and historical population declines of Serbian grayling do not differ from those of other European populations.

Conclusions

Our study shows that (1) the Ibar, Lim and Drina Rivers grayling populations are genetically distinct from populations outside of Serbia and thus should be managed as native populations in spite of some introgression in the Drina River population and (2) the Rzav River population is not appropriate for further stocking activities since it originates from stocked Slovenian grayling. However, the Rzav River population does not represent an immediate danger to other populations because it is physically isolated from these.

Unanticipated population structure of European grayling in its northern distribution: implications for conservation prioritization

BACKGROUND

The European grayling (Thymallus thymallus) is a salmonid fish native to Europe, with a distribution ranging from England and France to the Ural Mountains of north-western Russia. The majority of grayling populations inhabit freshwater rivers and lakes but some populations also occupy brackish water in northern parts of the Baltic Sea. Previous population genetic studies have demonstrated that grayling populations in Finland, Estonia and Russia belong to a single mitochondrial lineage and exhibit high levels of differentiation even at a small geographic scale. As a result, we predicted that grayling populations should not cluster regionally. Despite the extensive amount of genetic research that has been carried out on grayling, comprehensive national-level information on population structure of grayling in Northern Europe is still lacking. Yet this is the level at which populations are currently managed.

RESULTS

We found unanticipated population structure of grayling clustering into three groups largely corresponding to the northern, Baltic and south-eastern geographic areas of Finland using 13 microsatellite loci. We also found a high level of genetic differentiation among the groups and moderate to high differentiation within the groups. This combined with low variability strongly indicates that genetic drift and limited migration have a major impact on grayling population structure. An allele size permutation test indicated that mutations at microsatellite loci have not significantly contributed to genetic differentiation among the three Finnish groups. However, at the European scale, mutations had significantly contributed to population differentiation.

CONCLUSION

This research provides novel genetic information on European grayling in its northern distribution range and has clear implications for supporting country-scale conservation efforts. Specifically, the strong between population divergence observed indicates that single populations should generally be recognized as separate management units. We also introduced an alternative prioritization strategy for population conservation based on the evaluation of the relative roles of different evolutionary forces shaping the gene pools. We envision that the proposed approach to categorize populations for conservation will be a useful tool for wildlife researchers and conservationists working on a diverse range of organisms.

Reproductive isolation, evolutionary distinctiveness and setting conservation priorities: the case of European lake whitefish and the endangered North Sea houting (Coregonus spp.)

Background

Adaptive radiation within fishes of the Coregonus lavaretus complex has created numerous morphs, posing significant challenges for taxonomy and conservation priorities. The highly endangered North Sea houting (C. oxyrhynchus; abbreviated NSH) has been considered a separate species from European lake whitefish (C. lavaretus; abbreviated ELW) due to morphological divergence and adaptation to oceanic salinities. However, its evolutionary and taxonomic status is controversial. We analysed microsatellite DNA polymorphism in nine populations from the Jutland Peninsula and the Baltic Sea, representing NSH (three populations, two of which are reintroduced) and ELW (six populations). The objectives were to: 1) analyse postglacial recolonization of whitefish in the region; 2) assess the evolutionary distinctiveness of NSH, and 3) apply several approaches for defining conservation units towards setting conservation priorities for NSH.

Results

Bayesian cluster analyses of genetic differentiation identified four major groups, corresponding to NSH and three groups of ELW (Western Jutland, Central Jutland, Baltic Sea). Estimates of historical migration rates indicated recolonization in a north-eastern direction, suggesting that all except the Baltic Sea population predominantly represent postglacial recolonization via the ancient Elbe River. Contemporary gene flow has not occurred between NSH and ELW, with a divergence time within the last 4,000 years suggested from coalescence methods. NSH showed interbreeding with ELW when brought into contact by stocking. Thus, reproductive isolation of NSH was not absolute, although possible interbreeding beyond the F1 level could not be resolved.

Conclusion

Fishes of the C. lavaretus complex in the Jutland Peninsula originate from the same recolonization event. NSH has evolved recently and its species status may be questioned due to incomplete reproductive isolation from ELW, but it was shown to merit consideration as an independent conservation unit. Yet, application of several approaches for defining conservation units generated mixed outcomes regarding its conservation priority. Within the total species complex, it remains one among many recently evolved unique forms. Its uniqueness and high conservation priority is more evident at a local geographical scale, where conservation efforts will also benefit populations of a number of other endangered species.

Phylogenetic analysis of the genus Thymallus (grayling) based on mtDNA control region and ATPase 6 genes, with inferences on control region constraints and broad-scale Eurasian phylogeography

We present first insights into the molecular phylogeny of the grayling genus Thymallus (Salmonidae) using sequences from the mitochondrial control region and ATPase6 genes. A suite of analytical approaches were applied for each gene separately and for the combined data. The ATPase6 gene is shown to have a mean divergence rate across the genus of 2.46 times faster than the complete control region. Based on the combined data, four major (internal) clades, presumably originating in the Pliocene, were resolved with high support in all analyses and represented two distinct lineages in the Amur basin, one lineage in all remaining Siberian and Mongolian drainages, and one lineage corresponding to European grayling Thymallus thymallus. The resolution of multiple lineages, from both additional internal and terminal clades, within each major drainage basin underscores the complexity and effects that Pleistocene hydrological dynamics have had on the distribution of biodiversity in Siberia.

Experience from Lipizzan horse and salmonid species endemic to the Adriatic river system. Examples for the application of molecular markers for preservation of biodiversity and management of animal genetic resources

Management of breeding- and free-living populations, traditionally based on phenotypic traits, relays more and more on availability of reliable information about the basic population genetic parameters as heterozygosity, mean number of alleles per locus, percentage of polymorphic loci, population structuring, genetic distances and others. Therefore, the application of molecular markers, revealing a great deal of phenotypically hidden information, becomes inevitable for population analysis. Conservation geneticists use this information for implementation of appropriate management policies. Application of molecular markers in Lipizzan horse breed, which is an example for a pedigreed breeding population, and in two endangered salmonid fish populations in Slovenia, are presented. In the Lipizzan horse breed, an insight in the population structure, overall heterozygosity, relationship between population parameters and phenotypic traits and reliability of pedigree data was gained by using molecular markers. Marble trout and Adriatic grayling were selected as examples for free-living populations, seriously endangered by human activity in the past. Development of informative molecular markers, their application and suggestions for appropriate conservation actions are described.

Genetic introgression between wild and stocked salmonids and the prospects for using molecular markers in population rehabilitation: the case of the Adriatic grayling (Thymallus thymallus L. 1785)

In the north Adriatic basin, a morphologically and genetically distinct lineage of grayling is found, designated as the Adriatic grayling. In Slovenia, the Adriatic grayling is restricted to the Soca river system, where it is critically endangered. The most pertinent threat is stocking with non-native, highly divergent Sava (Danubian) drainage stock, and this activity has been going on for more than four decades. The present study was designed to characterise the genetic structure of the Adriatic grayling in Slovenia, with particular emphasis on estimating the degree of introgression with non-indigenous stocked grayling. We analysed polymorphism at 154 microsatellite loci in samples representing grayling from the Adriatic and Danubian drainage stock. A relatively high number (12) of alleles, diagnostic for the Adriatic grayling, were identified. However, a correspondence analysis based on individual multilocus genotypes also revealed that there is no distinctive Adriatic group but rather a dispersed multitude of individuals that cannot be unambiguously distinguished from the more homogenous Danubian population. A Bayesian analysis of individual admixture coefficients confirmed this pattern and revealed extensive introgression between the Adriatic grayling and stocked grayling of Danubian origin. Average individual admixture coefficients showed that only between 50 and 60% of the original gene pools remained, and only few non-introgressed indigenous individuals could be identified. Microsatellite-based individual admixture analysis appear to be an important tool for identifying remaining non-introgressed indigenous individuals that could be used for restoring the original populations.

Identification of reproductively isolated lineages of Amur grayling (Thymallus grubii Dybowski 1869): concordance between phenotypic and genetic variation

We analysed variation at maternally (mitochondrial DNA control region sequences) and bi-parentally (10 microsatellites) inherited genetic markers, as well as across 12 meristic characters in 7 populations of Amur grayling, Thymallus grubii, from eastern Siberia. All three data sets were concordant in supporting the existence of three diagnosable, reciprocally monophyletic, and most probably reproductively isolated, lineages of grayling within the Amur drainage. There was a significant correlation between genetic and phenotypic divergence, both within and among lineages. Two phenotypically distinct forms (with and without an orange spot on the posterior portion of the dorsal fin), found in sympatry in the lower Amur, most likely result from secondary contact, as they demonstrate 4.6% sequence divergence at the mitochondrial DNA control region. This divergence, together with the existence of at least one nearby population of orange spot grayling outside the Amur drainage (0.8% divergence) underscore the palaeo-hydrological complexity of the system, which presumably promoted genetic divergence in a shifting allopatric framework throughout the Pleistocene. Grayling from the upper Amur, corresponding to the type locality for the species, formed a sister group (1.4-1.6% divergent) to the orange spot lineage perhaps diverging in the early Pleistocene (1.4-1.6 Ma).

Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus, based on complete sequencing of the mitochondrial DNA control region

The complete mitochondrial DNA (mtDNA) control region (1043 base pairs) and 162-bp of flanking transfer RNA genes were sequenced in 316 European grayling, Thymallus thymallus, from 44 populations throughout the Western European range of the species. A total of 58 haplotypes were revealed with pairwise divergence ranging from 0.001 to 0.038. An inferred intraspecific phylogenetic tree revealed two well-supported clades within the Danube basin, one highly divergent clade in the Adriatic basin, and one large, diverse group representing most other populations. A deeply divergent haplotype fixed in the Loire basin in central France, more groups of haplotypes from distinct Danubian tributaries, and a relatively ancestral haplotype fixed in former tributaries of the Elbe in Denmark all suggest a complex pattern of interglacial and postglacial expansions originating from disjunct refugia throughout central Europe. Despite some evidence of human-mediated stock transfers, parsimony-network-based nested-clade analysis (NCA) supported specific inferences relating to corridors of postglacial expansion such as the lower Rhine (Moselle) and Elbe systems (Danish populations) serving as sources for expansion into the Baltic to the north as well as the upper Rhine and Danube to the south; and specific Rhine populations (Doller, Orbe and Reuss) serving as sources for colonization of the Rhone. The multiple divergent clades representing populations in the upper Danube, as well as the deeply divergent haplotypes found in the Adriatic and Loire basins (> 5% divergence from Asian outgroups) support the theory that European grayling have had a long history in Western Europe, pre-dating Pleistocene glacial cycles. The patterns of mtDNA divergence shown here support a perspective of rich inter- and intrabasin genetic diversity that should be protected from current trends to translocate brood stocks for rearing and release in response to declining populations, especially in southern European basins.

Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae

The phylogeography of an endangered salmonid, European grayling (Thymallus thymallus), was studied based on analysis of 17 nuclear microsatellite DNA loci. In agreement with earlier mitochondrial DNA (mtDNA) studies, phylogenetic relationships of the populations suggested that northern Europe was colonized from two distinct Pleistocene refugia. Furthermore, microsatellites revealed highly supported grouping of mainland Swedish, Norwegian, Danish, German and Slovenian populations, suggesting that grayling from the northwestern and central Europe have descended from their southern conspecifics. The level of divergence between populations was substantial, even across short geographical distances. Although this was in part due to postglacial colonization patterns and contemporary barriers for gene flow, the high divergence estimates between hydrologically connected sampling locations implied efficient interpopulation reproductive isolation. Microsatellites revealed that the populations exhibited, on average, only 3.5 (+/-2.2) alleles per locus, indicating that T. thymallus has strikingly low levels of intrapopulation genetic diversity as compared with other freshwater fish species. Accordingly, as indicated by analysis of molecular variance (AMOVA), only 49.1-58.0% of the total grayling microsatellite diversity resided within populations. A latitudinal genetic diversity gradient, potentially resulting from glaciation-mediated founder events, was not evident. Alternatively, it is possible that grayling display limited dispersal behaviour/capability, leading to low long-term effective population sizes and, consequently, depauperate intrapopulation polymorphism. These findings have implications for conservation of T. thymallus. Importantly, they exemplify that microsatellites can be highly informative for intraspecific phylogeography studies dealing with substantial divergence scales.