Demographic histories and genetic diversities of Fennoscandian marine and landlocked ringed seal subspecies

Island populations are on average smaller, genetically less diverse, and at a higher risk to go extinct than mainland populations. Low genetic diversity may elevate extinction probability, but the genetic component of the risk can be affected by the mode of diversity loss, which, in turn, is connected to the demographic history of the population. Here, we examined the history of genetic erosion in three Fennoscandian ringed seal subspecies, of which one inhabits the Baltic Sea ‘mainland’ and two the ‘aquatic islands’ composed of Lake Saimaa in Finland and Lake Ladoga in Russia. Both lakes were colonized by marine seals after their formation c. 9500 years ago, but Lake Ladoga is larger and more contiguous than Lake Saimaa. All three populations suffered dramatic declines during the 20th century, but the bottleneck was particularly severe in Lake Saimaa. Data from 17 microsatellite loci and mitochondrial control‐region sequences show that Saimaa ringed seals have lost most of the genetic diversity present in their Baltic ancestors, while the Ladoga population has experienced only minor reductions. Using Approximate Bayesian computing analyses, we show that the genetic uniformity of the Saimaa subspecies derives from an extended founder event and subsequent slow erosion, rather than from the recent bottleneck. This suggests that the population has persisted for nearly 10,000 years despite having low genetic variation. The relatively high diversity of the Ladoga population appears to result from a high number of initial colonizers and a high post‐colonization population size, but possibly also by a shorter isolation period and/or occasional gene flow from the Baltic Sea.     

COLONIZATION HISTORY OF THE BALTIC HARBOR SEALS: INTEGRATING ARCHAEOLOGICAL, BEHAVIORAL, AND GENETIC DATA

Detailed knowledge about the history of colonization, population dynamics and behavior greatly enhance evaluation of genetic models of population units and migration rates in spatially structured populations. Here, the genetic uniqueness of harbor seals ( Phoca vitulinia ) in the eastern Baltic is evaluated in the light of new information on the distribution and abundance of Baltic and eastern North Sea populations during the last 11,000 yr, recent hunting statistics, and population counts. Archaeological records reveal that the Baltic population of harbor seals was founded about 8,000 yr ago. Adjacent populations in the North Sea areas were either small, or went extinct, and became significant only during the last 300 yr. This information generates the hypothesis that the Baltic population has been isolated during the last 8,000 yr, despite the lack of geographical barriers. We show that stochastic effects, isolation, and a documented recent population bottleneck can account for the low observed genetic variation in Baltic harbor seals.     

An environmental gradient dominates ecological and genetic differentiation of marine invertebrates between the North and Baltic Sea

Environmental gradients have emerged as important barriers to structuring populations and species distributions. We set out to test whether the strong salinity gradient from the marine North Sea to the brackish Baltic Sea in northern Europe represents an ecological and genetic break, and to identify life history traits that correlate with the strength of this break. We accumulated mitochondrial cytochrome oxidase subunit 1 sequence data, and data on the distribution, salinity tolerance, and life history for 28 species belonging to the Cnidaria, Crustacea, Echinodermata, Mollusca, Polychaeta, and Gastrotricha. We included seven non‐native species covering a broad range of times since introduction, in order to gain insight into the pace of adaptation and differentiation. We calculated measures of genetic diversity and differentiation across the environmental gradient, coalescent times, and migration rates between North and Baltic Sea populations, and analyzed correlations between genetic and life history data. The majority of investigated species is either genetically differentiated and/or adapted to the lower salinity conditions of the Baltic Sea. Species exhibiting population structure have a range of patterns of genetic diversity in comparison with the North Sea, from lower in the Baltic Sea to higher in the Baltic Sea, or equally diverse in North and Baltic Sea. Two of the non‐native species showed signs of genetic differentiation, their times since introduction to the Baltic Sea being about 80 and >700 years, respectively. Our results indicate that the transition from North Sea to Baltic Sea represents a genetic and ecological break: The diversity of genetic patterns points toward independent trajectories in the Baltic compared with the North Sea, and ecological differences with regard to salinity tolerance are common. The North Sea–Baltic Sea region provides a unique setting to study evolutionary adaptation during colonization processes at different stages by jointly considering native and non‐native species.     

Microsatellite and mtDNA analysis of the population structure of grey seals (Halichoerus grypus) from three breeding areas in the Baltic Sea

The growing number of grey seals in the Baltic Sea has led to a dramatic increase in interactions between seals and fisheries. The conflict has become such a problem that hunting was introduced in Finland in 1998 and the Swedish Environment Protection Agency recommended a cull of grey seals starting in 2001. Culling has been implemented despite the lack of data on population structure. Low levels of migration between regions would mean that intensive culling in specific geographic areas would have disproportionate effects on local population structure and genetic diversity. We used eight microsatellite loci and a 489 bp section of the mtDNA control region to examine the genetic variability and differentiation between three breeding sites in the Baltic Sea and two in the UK. We found high levels of genetic variability in all sampled Baltic groups for both the microsatellites and the control region. There were highly significant differences in microsatellite allele frequencies between all three Baltic breeding sites and between the Baltic sites and the UK sites. However, there were no significant differences in mtDNA control region haplotypes between the Baltic sites. This genetic substructure of the Baltic grey seal populations should be taken into consideration when managing the seal population to prevent the hunting regime from having an adverse effect on genetic diversity by setting hunting quotas separately for the different subpopulations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Outlier Loci Detect Intraspecific Biodiversity amongst Spring and Autumn Spawning Herring across Local Scales

Herring, Clupea harengus, is one of the ecologically and commercially most important species in European northern seas, where two distinct ecotypes have been described based on spawning time; spring and autumn. To date, it is unknown if these spring and autumn spawning herring constitute genetically distinct units. We assessed levels of genetic divergence between spring and autumn spawning herring in the Baltic Sea using two types of DNA markers, microsatellites and Single Nucleotide Polymorphisms, and compared the results with data for autumn spawning North Sea herring. Temporally replicated analyses reveal clear genetic differences between ecotypes and hence support reproductive isolation. Loci showing non-neutral behaviour, so-called outlier loci, show convergence between autumn spawning herring from demographically disjoint populations, potentially reflecting selective processes associated with autumn spawning ecotypes. The abundance and exploitation of the two ecotypes have varied strongly over space and time in the Baltic Sea, where autumn spawners have faced strong depression for decades. The results therefore have practical implications by highlighting the need for specific management of these co-occurring ecotypes to meet requirements for sustainable exploitation and ensure optimal livelihood for coastal communities.     

Temporal genetic analysis of the endangered tidewater goby: extinction–colonization dynamics or drift in isolation?

Extinction and colonization dynamics are critical to understanding the evolution and conservation of metapopulations. However, traditional field studies of extinction–colonization are potentially fraught with detection bias and have rarely been validated. Here, we provide a comparison of molecular and field‐based approaches for assessment of the extinction–colonization dynamics of tidewater goby (Eucyclogobius newberryi) in northern California. Our analysis of temporal genetic variation across 14 northern California tidewater goby populations failed to recover genetic change expected with extinction–colonization cycles. Similarly, analysis of site occupancy data from field studies (94 sites) indicated that extinction and colonization are very infrequent for our study populations. Comparison of the approaches indicated field data were subject to imperfect detection, and falsely implied extinction–colonization cycles in several instances. For northern California populations of tidewater goby, we interpret the strong genetic differentiation between populations and high degree of within‐site temporal stability as consistent with a model of drift in the absence of migration, at least over the past 20–30 years. Our findings show that tidewater goby exhibit different population structures across their geographic range (extinction–colonization dynamics in the south vs. drift in isolation in the north). For northern populations, natural dispersal is too infrequent to be considered a viable approach for recolonizing extirpated populations, suggesting that species recovery will likely depend on artificial translocation in this region. More broadly, this work illustrates that temporal genetic analysis can be used in combination with field data to strengthen inference of extinction–colonization dynamics or as a stand‐alone tool when field data are lacking.     

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.     

Phylogeographical structure, distribution and genetic variation of the green algae Ulva intestinalis and U. compressa (Chlorophyta) in the Baltic Sea area

The marine algae Ulva intestinalis and U. compressa are morphologically plastic with many overlapping characters and are therefore difficult to distinguish from each other. The present distribution of U. intestinalis and U. compressa is investigated along the salinity gradient in the Baltic Sea area through analyses of internal transcribed spacer (ITS) sequence data. Also, the amount and distribution of intraspecific genetic polymorphism in the ITS region is studied allowing inferences on the phylogeographical pattern and postglacial recolonization of the Baltic Sea area. The data show that of the two species only U. intestinalis occurs in the Baltic Sea. The distribution of U. compressa is more restricted than previously reported, and it was not found in salinities lower than 15 ppt. All of Scandinavia and the Baltic Sea were covered with ice during the last ice age and the organisms in the Baltic Sea must have colonized the area after the ice had started to melt. The genetic diversity of U. intestinalis and U. compressa in the Baltic Sea and the neighbouring area was found to be reduced compared to that in the British Isles. This reduction may be the result of either a historical reduction of diversity or an adaptation of specific clones to the northern environmental conditions.     

A PREHISTORIC BREEDING POPULATION OF HARP SEALS (PHOCA GROENLANDICA) IN THE BALTIC SEA

The pelagic and gregarious, low Arctic harp seal ( Phoca groenlandica ) is the most common seal species in most refuse faunas from coastal hunter-gatherer sites dating from the late Atlantic to the early Subboreal period ( ca. 4000-2000 cal B. C.) in the Baltic Sea. Our main objective was to examine the migration contra breeding population hypotheses regarding the Baltic harp seals. Analyses of epiphyseal fusion data and osteometry of archeological harp seal remains from 25 dwelling-sites suggest that a local breeding population established itself in the early Subboreal period. In the Middle Neolithic the rookery possibly was situated in the Baltic proper, south of Aland and west of Gotland. The mean adult size of the Baltic harp seals decreased, suggesting minimal genetic exchange with the north Atlantic Ocean population. Genetic drift, interspecific competition, and over-hunting by humans are all factors likely to have contributed to the eventual extinction of harp seals in the Baltic Sea.     

Adaptive and Neutral Genetic Variation and Colonization History of Atlantic Salmon, Salmo salar

Synopsis I combined neutral microsatellite markers with the major histocompatibility complex (MHC) class IIB to study genetic differentiation and colonization history in Atlantic salmon, Salmo salar, in the Baltic Sea and in the north-eastern Atlantic. Baltic salmon populations have lower levels of microsatellite genetic variation, in terms of heterozygosity and allelic richness than Atlantic populations, confirming earlier findings with other genetic markers, suggesting that the Baltic Sea populations have been exposed to genetic bottlenecks, most likely at a founding event. On the other hand, the level of MHC variation was similar in the Baltic and in the north-eastern Atlantic, indicating that positive balancing selection has increased the level of MHC-variation. Both microsatellite and MHC class IIB genetic variation give strong support to the hypothesis that the Baltic salmon are of a biphyletic origin, the southern population in this study is strongly differentiated from both the northern Baltic salmon populations and from the north-eastern Atlantic populations. Salmon may have colonized the northern Baltic Sea either from the south, via the so called “N?rke strait” or from the north, via a proposed historical connection between the White Sea and the northern Baltic. At microsatellites, no significant isolation-by distance was found at either colonization route. At the MHC, populations were significantly isolated by distance when assuming that colonization occurred via the “N?rke strait”.     

Further expansion of the genus Cercopagis (Crustacea,Branchiopoda, Onychopoda) in the Baltic Sea,with notes on the taxa present and their ecology

The onychopod cladoceran Cercopagis that recently invaded the Baltic Sea is reported from new zones of the northern Baltic proper. Because of successful survival and an expanding distribution range, the addition of Cercopagis to the Baltic fauna is considered to be permanent. What has previously been cited as Cercopagis pengoi encompasses the morphology of several other species, subspecies and forms. Either a number of morphologically similar species is present, or there is a number of spurious species in Cercopagis. The last hypothesis is favoured. The spatial distribution pattern of Cercopagis, as well as that of total zooplankton, was correlated with depth. Deep (>100 m) and shallow (<10 m) stations had significantly lower abundance than stations of intermediate depth (<100 m). An overview of the distribution of C. pengoi group in fresh and brackish waters suggests a high tolerance to environmental factors, but with differences among taxa. Due to this ecological flexibility, the colonization of the Baltic is not unexpected. Increasing salinity may restrict dispersal of cercopagids to the southern areas of the Baltic and to the North Sea, but inland lakes (e.g. in Sweden) present an ecological profile suitable for colonization.     

The 1988 and 2002 phocine distemper virus epidemics in European harbour seals

We present new and revised data for the phocine distemper virus (PDV) epidemics that resulted in the deaths of more than 23 000 harbour seals Phoca vitulina in 1988 and 30,000 in 2002. On both occasions the epidemics started at the Danish island of Anholt in central Kattegat, and subsequently spread to adjacent colonies in a stepwise fashion. However, this pattern was not maintained throughout the epidemics and new centres of infection appeared far from infected populations on some occasions: in 1988 early positive cases were observed in the Irish Sea, and in 2002 the epidemic appeared in the Dutch Wadden Sea, 6 wk after the initiation of the outbreak at Anholt Island. Since the harbour seal is a rather sedentary species, such 'jumps' in the spread among colonies suggest that another vector species could have been involved. We discussed the role of sympatric species as disease vectors, and suggested that grey seal populations could act as reservoirs for PDV if infection rates in sympatric species are lower than in harbour seals. Alternatively, grey seals could act as subclinical infected carriers of the virus between Arctic and North Sea seal populations. Mixed colonies of grey and harbour seal colonies are found at all locations where the jumps occurred. It seems likely that grey seals, which show long-distance movements, contributed to the spread among regions. The harbour seal populations along the Norwegian coast and in the Baltic escaped both epidemics, which could be due either to genetic differences among harbour seal populations or to immunity. Catastrophic events such as repeated epidemics should be accounted for in future models and management strategies of wildlife populations.     

Phylogeographic evidence for the postglacial colonization of the North and Baltic Sea coasts from inland glacial refugia by Triglochin maritima L.

We investigated the geographical distribution of genetic variation in 67 individuals of Triglochin maritima from 38 localities across Europe using AFLP markers. Analysis of genetic variation resulted in the recognition of two major genetic groups. Apart from few geographical outliers, these are distributed (1) along the Atlantic coasts of Portugal, Spain and France and (2) in the North Sea area, the Baltic Sea area, at central European inland localities, the northern Adriatic Sea coast and the Mediterranean coast of southwest France. Considering possible range shifts of T. maritima in reaction to Quaternary climatic changes as deduced from the present-day northern temperature limit of the species, Quaternary changes of coastline in the North Sea area and the very recent origin of the Baltic Sea, we conclude that the coastal populations of T. maritima in the North Sea and Baltic Sea areas originated from inland populations.     

Combined Genetic and Telemetry Data Reveal High Rates of Gene Flow,Migration, and Long-Distance Dispersal Potential in Arctic Ringed Seals (Pusa hispida)

Ringed seals (Pusa hispida) are broadly distributed in seasonally ice covered seas, and their survival and reproductive success is intricately linked to sea ice and snow. Climatic warming is diminishing Arctic snow and sea ice and threatens to endanger ringed seals in the foreseeable future. We investigated the population structure and connectedness within and among three subspecies: Arctic (P. hispida hispida), Baltic (P. hispida botnica), and Lake Saimaa (P. hispida saimensis) ringed seals to assess their capacity to respond to rapid environmental changes. We consider (a) the geographical scale of migration, (b) use of sea ice, and (c) the amount of gene flow between subspecies. Seasonal movements and use of sea ice were determined for 27 seals tracked via satellite telemetry. Additionally, population genetic analyses were conducted using 354 seals representative of each subspecies and 11 breeding sites. Genetic analyses included sequences from two mitochondrial regions and genotypes of 9 microsatellite loci. We found that ringed seals disperse on a pan-Arctic scale and both males and females may migrate long distances during the summer months when sea ice extent is minimal. Gene flow among Arctic breeding sites and between the Arctic and the Baltic Sea subspecies was high; these two subspecies are interconnected as are breeding sites within the Arctic subspecies.     

Colonization and/or mitochondrial selective sweeps across the North Atlantic intertidal assemblage revealed by multi‐taxa approximate Bayesian computation

Intertidal and subtidal communities of the western and eastern coasts of the North Atlantic Ocean were greatly affected by Pleistocene glaciations, with some taxa persisting on both coasts, and others recolonizing after being extirpated on one coast during the Last Glacial Maximum. In the original spirit of comparative phylogeography, we conducted a comparative analysis using mtDNA sequence data and a hierarchical approximate Bayesian computation (ABC) approach for testing these two scenarios across 12 intertidal and subtidal coastal invertebrates spanning the North Atlantic to determine the temporal dynamics of species membership of these two ephemeral communities. Conditioning on a low gene‐flow model, our results suggested that a colonization or mitochondrial selective sweep history was predominant across all taxa, with only the bivalve mollusc Mytilus edulis showing a history of trans‐Atlantic persistence. Conditioning on a high gene‐flow model weakened the support for this assemblage‐level demographic history. The predominance of a colonization‐type history also highlights concerns about analyses based on single‐locus data where genetic hitchhiking may be incorrectly inferred as colonization. In conclusion, driving factors in shifting species range distributions and membership of ephemeral coastal communities could be species‐specific environmental tolerances, species interactions, and/or stochastic demographic extinction. Through a re‐examination of a long‐standing question of North Atlantic phylogeography, we highlight the flexibility and statistical honesty of using a model‐based ABC approach.     

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene‐linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high‐resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.     

Local extinction and recolonization, species effective population size, and modern human origins

A primary objection from a population genetics perspective to a multiregional model of modern human origins is that the model posits a large census size, whereas genetic data suggest a small effective population size. The relationship between census size and effective size is complex, but arguments based on an island model of migration show that if the effective population size reflects the number of breeding individuals and the effects of population subdivision, then an effective population size of 10,000 is inconsistent with the census size of 500,000 to 1,000,000 that has been suggested by archeological evidence. However, these models have ignored the effects of population extinction and recolonization, which increase the expected variance among demes and reduce the inbreeding effective population size. Using models developed for population extinction and recolonization, we show that a large census size consistent with the multiregional model can be reconciled with an effective population size of 10,000, but genetic variation among demes must be high, reflecting low interdeme migration rates and a colonization process that involves a small number of colonists or kin-structured colonization. Ethnographic and archeological evidence is insufficient to determine whether such demographic conditions existed among Pleistocene human populations, and further work needs to be done. More realistic models that incorporate isolation by distance and heterogeneity in extinction rates and effective deme sizes also need to be developed. However, if true, a process of population extinction and recolonization has interesting implications for human demographic history.