Genetic evidence for different migration routes of freshwater fish into Norway revealed by analysis of current perch (Perca fluviatilis) populations in Scandinavia

Ancient DNA and osteological analyses of a unique paleo-archive reveal Early Holocene faunal expansion into the Scandinavian Arctic

Paleo-archives are essential for our understanding of species responses to climate warming, yet such archives are extremely rare in the Arctic. Here, we combine morphological analyses and bulk-bone metabarcoding to investigate a unique chronology of bone deposits sealed in the high-latitude Storsteinhola cave system (68°50' N 16°22' E) in Norway. This deposit dates to a period of climate warming from the end of the Late Glacial [~13 thousand calibrated years before the present (ka cal B.P.)] to the Holocene thermal maximum (~5.6 ka cal B.P.). Paleogenetic analyses allow us to exploit the 1000s of morphologically unidentifiable bone fragments resulting in a high-resolution sequence with 40 different taxa, including species not previously found here. Our record reveals borealization in both the marine and terrestrial environments above the Arctic Circle as a naturally recurring phenomenon in past periods of warming, providing fundamental insights into the ecosystem-wide responses that are ongoing today.

Trends of Eurasian Perch (Perca fluviatilis) mtDNA ATP6 Region Genetic Diversity within the Hydro-Systems of the Eastern Part of the Baltic Sea in the Anthropocene

The intraspecific genetic diversity of freshwater fish inhabiting hydro-systems of the macrogeographic area spreading from the Black to Baltic Seas requires comprehensive investigation from fundamental and practical perspectives. The current study focused on the involvement of the mtDNA ATP6 region in the adaptability and microevolution of Perca fluviatilis within phylogeographic and anthropogenic contexts. We sequenced a 627 bp fragment encompassing the ATP6 region and used it for genetic analysis of 193 perch caught in Latvia, Lithuania, Belarus, and Ukraine, representing natural and anthropogenically impacted populations. We evaluated patterns of intraspecific genetic diversity in the ATP6 region and phylogeographic trends within the studied area compared with previously established D-loop trends. Evaluation of ATP6 coding sequence variability revealed that among 13 newly detected haplotypes, only two were caused by non-synonymous substitutions of amino acids of the protein. PCoA revealed three genetic groups (I-III) based on the ATP6 region that encompassed four previously described genetic groups established based on the mtDNA D-loop. The two mtDNA regions (D-loop and ATP6) have microevolved at least partially independently. Prolonged anthropogenic impacts may generate new point mutations at the ATP6 locus, but this phenomenon could be mainly concealed by natural selection and reparation processes.

Shallow phylogeographic structuring of Vimba vimba across Europe suggests two distinct refugia during the last glaciation

Genetic variation and geographical structuring of vimba Vimba vimba were analysed across 26 sites (80 individuals) by means of mtDNA sequences (cyt b gene, mitochondrial control region) to localize hypothesized glacial refugia and to reconstruct postglacial recoloniation routes. Although genetic diversity among sequenced individuals was low, a combined analysis of the two sequenced fragments revealed a western (central and northern Europe: Danube, Elbe and lakes of Sweden) and an eastern clade (eastern Europe: Dnieper-South Bug, Don, Neman). Furthermore, a number of divergent ancestral haplotypes distributed around the Black and Caspian Seas became apparent. Mismatch analyses supported a sudden expansion model for the populations of the western clade between 50 and 10 000 bp. Overall, the study provides strong evidence for a northward and westward expansion of V. vimba from two refugial regions located in the Danubian drainage and the northern Pontic regions respectively.

EXTRINSIC AND INTRINSIC CONTROLS OF ZOOPLANKTON DIVERSITY IN LAKES

Pelagic crustacean zooplankton were collected from 336 Norwegian lakes covering a wide range of latitude, altitude, lake area, mean depth, production (as chlorophyll a), and fish community structure. Mean zooplankton species richness during the ice-free season was generally low at high latitudes and altitudes. Further, lower species richness was recorded in western lakes, possibly reflecting constraints on migration and dispersal. However, despite obvious spatial limitations, geographic boundaries were only weak predictors of mean zooplankton richness. Similarly, lake surface area did not contribute positively to mean richness such as seen in other ecosystem surveys. Rather, intrinsic factors such as primary production and fish community (planktivore) structure were identified by regression analysis as the major predictors of zooplankton diversity, while a positive correlation was observed between species richness and total zooplankton biomass. However, in spite of a large number of variables included in this study, the predictive power of multiple regression models was modest (<50% variance explained), pointing to a major role for within-lake properties, as yet unidentified intrinsic forces, stochasticity, or dispersal as constraints on zooplankton diversity in these lakes.

Postglacial colonization shows evidence for sympatric population splitting of Eurasian perch (Perca fluviatilis L.) in Lake Constance

Previous microsatellite analysis showed that two subpopulations of perch (Perca fluviatilis L) exist in Lake Constance. This raises questions of whether (i). Lake Constance was colonized by two populations that diverged in allopatry, or (ii). the two subpopulations diverged in sympatry. Sequence analysis of a 365 bp mtDNA fragment (5'-end of the D-loop) of perch from Lake Constance and adjacent waters revealed 10 haplotypes. We suggest colonization via the Danube river, based on the frequency and dispersion of haplotypes, and knowledge of the lake's palaeohydrological development. Pairwise FST-values using mitochondrial DNA sequences showed no significant population subdivision. Our study provides strong evidence that subpopulations of perch in Lake Constance have diverged in sympatry.

Phylogeography of competing sexual and parthenogenetic forms of a freshwater flatworm: patterns and explanations

BACKGROUND

Models of the maintenance of sex predict that one reproductive strategy, sexual or parthenogenetic, should outcompete the other. Distribution patterns may reflect the outcome of this competition as well as the effect of chance and historical events. We review the distribution data of sexual and parthenogenetic biotypes of the planarian Schmidtea polychroa.

RESULTS

S. polychroa lives in allopatry or sympatry across Europe except for Central and North-Western Europe, where sexual individuals have never been reported. A phylogenetic relationship between 36 populations based on a 385 bp fragment of the mitochondrial cytochrome oxidase I gene revealed that haplotypes were often similar over large geographic distances. In North Italian lakes, however, diversity was extreme, with sequence differences of up to 5% within the same lake in both sexuals and parthenogens. Mixed populations showed "endemic" parthenogenetic lineages that presumably originated from coexisting sexuals, and distantly related ones that probably result from colonization by parthenogens independent from sexuals.

CONCLUSIONS

Parthenogens originated repeatedly from sexuals, mainly in Italy, but the same may apply to other Mediterranean regions (Spain, Greece). The degree of divergence between populations suggests that S. polychroa survived the ice ages in separate ice-free areas in Central, Eastern and Southern Europe and re-colonised Europe after the retreat of the major glaciers. Combining these results with those based on nuclear markers, the data suggest that repeated hybridisation between sexuals and parthenogenetic lineages in mixed populations maintains high levels of genetic diversity in parthenogens. This can explain why parthenogens persist in populations that were originally sexual. Exclusive parthenogenesis in central and western populations suggests better colonisation capacity, possibly because of inbreeding costs as well as hybridisation of sexuals with parthenogens.

Phylogeography of the vairone (Leuciscus souffia, Risso 1826) in Central Europe

The vairone Leuciscus souffia is a cyprinid fish that inhabits river systems in and around the Alps. The complete mitochondrial DNA control region (945 bp) was sequenced in 295 vairone from 22 populations in Central Europe. A total of 51 haplotypes were identified with a maximum pairwise distance between haplotypes of 5.6%. Phylogenetic analyses revealed two major clades in L. souffia, an 'Italian' clade, and an 'Alpine' clade. Two hybrid zones exist, in the Mediterranean Alps and in the Soca basin. The position of the sister species of L. souffia, L. turskyi, to the 'Alpine' and the 'Italian' clade could not be resolved unambiguously. However, a linearized tree analysis indicated that L. turskyi represents a third lineage, that originated at the same time as the 'Alpine' and the 'Italian' clades of L. souffia. In the 'Alpine' clade two groups were resolved, a subclade with haplotypes from the Rhône and Var basins and a cluster with haplotypes from the Danube and Rhine systems. Our data suggest a long history of the vairone in Central Europe, predating Pleistocene glacial cycles. Two main refugia during glaciations must have existed, one in Italy and another one most probably in the Danube system. However, age estimates based on molecular clock calibrations suggest the survival of 'Alpine' haplotypes in several drainages during the last glaciation cycles. The Rhine system was only recently colonized.

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.

Genetic differentiation of a European caddisfly: past and present gene flow among fragmented larval habitats

We describe the genetic structure of a freshwater insect species and interpret it in terms of present-day population dynamics and possible postglacial colonization history. The sampling regime represented a large area of the species range in northwest Europe, particularly focusing on Britain, a region relatively neglected in molecular population genetic studies. Plectrocnemia conspersa generally showed low levels of genetic variation across the sampled populations (Nei's D = 0.0138) and subdivision was unrelated to the pattern of the drainage network. However, the results do suggest that populations across the region are not at equilibrium and that British populations still show effects of the recolonization of the species following the last glacial maximum. Levels of genetic diversity were lower in Britain than in mainland Europe. Two-dimensional scaling showed genetic differentiation between major regions and the pattern of genetic diversity indicates a more recent origin of populations in the north and west of the area compared with the south and east. We argue that, despite the highly fragmented larval habitat, dispersal over tens of kilometres is frequent. Over longer distances, however, P. conspersa does still show evidence of founder effects and postglacial range expansion into Britain.

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

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 mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance (AMOVA) indicates high genetic differentiation both among populations (PhiST = 0.511, P < 0.001) and among regions (PhiCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one-dimensional stepping stone models. A neighbour-joining tree recovered by (DAMBE) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.

Genetic lineages and postglacial colonization of grayling (Thymallus thymallus, Salmonidae) in Europe, as revealed by mitochondrial DNA analyses

In stark contrast to other species within the Salmonidae family, phylogeographic information on European grayling, Thymallus thymallus, is virtually nonexistent. In this paper, we utilized mitochondrial DNA polymerase chain reaction-restriction fragment length polymorphism (mtDNA PCR-RFLP) and sequence variation to infer the postglacial dispersal routes of T. thymallus into and within northern Europe, and to locate geographically, potential evolutionarily distinct populations. Mitochondrial analyses revealed a total of 27 T. thymallus haplotypes which clustered into three distinct lineages. Average pairwise interlineage divergence was four and nine times higher than average intralineage divergence for RFLP and sequence data, respectively. Two European grayling individuals from the easternmost sample in Russia exhibited haplotypes more genetically diverged from any T. thymallus haplotype than T. arcticus haplotype, and suggested that hybridization/introgression zone of these two sister species may extend much further west than previously thought. Geographic division of the lineages was generally very clear with northern Europe comprising of two genetically differentiated areas: (i) Finland, Estonia and north-western Russia; and (ii) central Germany, Poland and western Fennoscandia. Average interpopulation divergence in North European T. thymallus was 10 times higher than that observed in a recent mtDNA study of North American T. arcticus. We conclude that (i) North European T. thymallus populations have survived dramatic Pleistocene temperature oscillations and originate from ancient eastern and central European refugia; (ii) genetic divergence of population groups within northern Europe is substantial and geographically distinct; and (iii) the remainder of Europe harbours additional differentiated assemblages that likely descend from a Danubian refugium. These findings should provide useful information for developing appropriate conservation strategies for European grayling and exemplify a case with a clear need for multinational co-operation for managing and conserving biodiversity.

Phylogeography of the bullhead Cottus gobio (Pisces: Teleostei: Cottidae) suggests a pre-pleistocene origin of the major central European populations

The bullhead Cottus gobio is a small, bottom-dwelling fish consisting of populations that have not been subject to transplantations or artificial stocking. It is therefore an ideal model species for studying the colonization history of central European freshwater systems, in particular with respect to the possible influences of the Pleistocene glaciation cycles. We sampled Cottus populations across most of its distribution range, with a special emphasis on southern Germany where the major European drainage systems are in closest contact. Mitochondrial D-loop sequencing of more than 400 specimens and phylogenetic network analysis allowed us to draw a detailed picture of the colonization of Europe by C. gobio. Moreover, the molecular distances between the haplotypes enabled us to infer an approximate time frame for the origin of the various populations. The founder population of C. gobio stems apparently from the Paratethys and invaded Europe in the Pliocene. From there, the first colonization into central Europe occurred via the ancient lower Danube, with a separate colonization of the eastern European territories. During the late Pliocene, one of the central European populations must have reached the North Sea in a second step after which it then started to colonize the Atlantic drainages via coastal lines. Accordingly, we found very distinct populations in the upper and lower Rhine, which can be explained by the fact that the lower Rhine was disconnected from the upper Rhine until approximately 1 million years ago (Ma). More closely related, but still distinct, populations were found in the Elbe, the Main and the upper Danube, all presumably of Pleistocene origin. Intriguingly, they have largely maintained their population identity, despite the strong disturbance caused by the glaciation cycles in these areas. On the other hand, a mixing of populations during postglacial recolonization could be detected in the lower Rhine and its tributaries. However, the general pattern that emerges from our analysis suggests that the glaciation cycles did not have a major impact on the general population structure of C. gobio in central Europe.

Genetic divergence and phylogeographic relationships among european perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization

We used the widely distributed freshwater fish, perch (Perca fluviatilis), to investigate the postglacial colonization routes of freshwater fishes in Europe. Genetic variability within and among drainages was assessed using mitochondrial DNA (mtDNA) D-loop sequencing and RAPD markers from 55 populations all over Europe as well as one Siberian population. High level of structuring for both markers was observed among drainages and regions, while little differentiation was seen within drainages and regions. Phylogeographic relationships among European perch were determined from the distribution of 35 mtDNA haplotypes detected in the samples. In addition to a distinct southern European group, which includes a Greek and a southern Danubian population, three major groups of perch are observed: the western European drainages, the eastern European drainages including the Siberian population, and Norwegian populations from northern Norway, and western side of Oslofjord. Our data suggest that present perch populations in western and northern Europe were colonized from three main refugia, located in southeastern, northeastern and western Europe. In support of this, nested cladistic analysis of mtDNA clade and nested clade distances suggested historical range expansion as the main factor determining geographical distribution of haplotypes. The Baltic Sea has been colonized from all three refugia, and northeastern Europe harbours descendants from both eastern European refugia. In the upper part of the Danube lineages from the western European and the southern European refugia meet. The southern European refugium probably did not contribute to the recolonization of other western and northern European drainages after the last glaciation. However, phylogenetic analyses suggest that the southern European mtDNA lineage is the most ancient, and therefore likely to be the founder of all present perch lineages. The colonization routes used by perch probably also apply to other freshwater species with similar distribution patterns.

Phylogeography and postglacial dispersion of the chub (Leuciscus cephalus) in Europe

A phylogeographic analysis of mitochondrial DNA variation was performed in order to test the hypothesis of a postglacial recolonization of mid- and north-European rivers from a Danubian refuge. Over 345 chub specimens from European rivers covering most of the species' native range were investigated using 600 bp of the cytochrome b gene. Chub in European rivers belong to four highly divergent mitochondrial groups (lineages) differing by mean divergence estimates from 5.2% to 7.89%. These four lineages have a largely allopatric distribution, implying four geographical sets: two Mediterranean, and two north-European sets. This pattern provided strong evidence for: (i) the eradication of this species from most of Europe during maximum ice extent; (ii) its survival in four refugia (Adriatic side of the Balkans, eastern Greece (Aegean Rivers), southern tributaries of the Danube, and periphery of Black and Caspian Seas); (iii) a differential postglacial recolonization of mid- and northern Europe from the last two refugia only; (iv) the occurrence of this recolonization in two steps for the Danubian (western) lineage that entered western Europe (Rhine-Rhone-Loire drainages) during the Riss-Würm interglacial period and survived the last glaciation there before colonizing Garonne, UK and German drainages up to the Elbe during the Holocene; and (v) the occurrence of this recolonization in a single step for the Ponto-Caspian (eastern) lineage that entered the Baltic area as far as the Oder in the Holocene. Both lineages came into contact in the River Elbe without evident mixing.

Heteroplasmy, length and sequence variation in the mtDNA control regions of three percid fish species (Perca fluviatilis, Acerina cernua, Stizostedion lucioperca)

The nucleotide sequence of the control region and flanking tRNA genes of perch (Perca fluviatilis) mtDNA was determined. The organization of this region is similar to that of other vertebrates. A tandem array of 10-bp repeats, associated with length variation and heteroplasmy was observed in the 5' end. While the location of the array corresponds to that reported in other species, the length of the repeated unit is shorter than previously observed for tandem repeats in this region. The repeated sequence was highly similar to the Mt5 element which has been shown to specifically bind a putative D-loop DNA termination protein. Of 149 perch analyzed, 74% showed length variation heteroplasmy. Single-cell PCR on oocytes suggested that the high level of heteroplasmy is passively maintained by maternal transmission. The array was also observed in the two other percid species, ruffe (Acerina cernua) and zander (Stizostedion lucioperca). The array and the associated length variation heteroplasmy are therefore likely to be general features of percid mtDNAs. Among the perch repeats, the mutation pattern is consistent with unidirectional slippage, and statistical analyses supported the notion that the various haplotypes are associated with different levels of heteroplasmy. The variation in array length among and within species is ascribed to differences in predicted stability of secondary structures made between repeat units.