The genetic relationship between extirpated and contemporary Atlantic salmon Salmo salar L. lines from the southern Baltic Sea
© Bernaś et al. 2016
Received: 11 December 2015
Accepted: 22 March 2016
Published: 1 April 2016
The genetic relationship between original Atlantic salmon populations that are now extinct in the southern Baltic Sea and the present-day populations has long been controversial. To investigate and clarify this issue, we successfully genotyped individuals of the historical populations from the Oder and Vistula Rivers using DNA extracted from dried scales with the Atlantic salmon single nucleotide polymorphism array.
Our results showed a global F ST of 0.2515 for all pairs of loci, which indicates a high level of genetic differentiation among the groups analyzed in this study. Pairwise F ST values were significant for all comparisons and the highest values were found between present-day reintroduced Slupia River salmon and extinct Vistula River Atlantic salmon. Bayesian analysis of genetic structure revealed the existence of substructures in the extirpated Polish populations and three main clades among studied stocks.
The historical salmon population from the Oder River was genetically closer to present-day salmon from the Neman River than to the historical salmon from the Vistula River. Vistula salmon clearly separated from all other analyzed salmon stocks. It is likely that the origins of the Atlantic salmon population from the Morrum River and the Polish historical native populations are different.
Currently, the European distribution of the Atlantic salmon Salmo salar L. ranges from northern Portugal to the Pechora River in northwest Russia and encompasses Iceland, the British Isles and the Baltic Sea . Baltic salmon populations reproduce in about 45 river systems of which at least 29 hold native salmon populations or partly mixed following stocking practices . The native Atlantic salmon in Poland became extinct in the middle of the twentieth century. Salmon disappeared from the upper Vistula River in the 1950s, from its lower course after the Włocławek Dam was opened in the 1960s and from most Pomeranian rivers during the same period. Salmon also disappeared from the Drawa River (Oder basin) at the end of the 1980s.
Historically, the largest population of salmon ascended the Vistula River and migrated to the tributaries of the upper Vistula River . In the Oder River, salmon migrated upstream to the mountain tributaries, but due to dam construction on these rivers during the second half of the twentieth century, salmon spawning grounds were observed mainly in the Drawa River . The present Polish wild salmon population in the Slupia River originated from the Latvian Daugava River, from which individuals were collected for a recovery program that was launched in the 1990s . To date, the genetic make-up of the Polish salmon breeding population is still very similar to that of the Daugava stock .
Ongoing discussions about the post-glacial colonization routes that were followed by the Baltic salmon (Salmo salar) have, to date, resulted in three hypotheses. The first one argues that Baltic salmon may have derived from a refugial population that survived in eastern pre-glacial lakes and colonized the Baltic basin before a marine strait connected the Baltic Sea to the North Sea [5, 6]. The second hypothesis suggests that Baltic salmon have a western Atlantic origin , whereas the third one proposes that colonization spread from both east and west areas . Knowledge on how the southern Baltic Sea was colonized by Atlantic salmon populations is still scarce and the only published results concern populations from the southern Swedish rivers, Morrum and Eman  and the southeast Neman River . Until now, extinction of Polish populations has prevented analysis of these lines. However, preserved salmon scales and tissue fragments that were collected before 1970 can be used to extract DNA samples.
Basic statistics of five salmon populations from the southern Baltic Sea
Genetic diversity indices for the investigated salmon populations
During the last decade, many new fish passes have been constructed on Pomeranian rivers in the southern Baltic area (e.g. seven in the Slupia River basin) and a previous fish pass on the Włocławek dam in the Vistula River was rebuilt in 2015. These systems allow the passage of Atlantic salmons and have increased their migration but the recovery program that still continues is only moderately successful . In the Baltic Sea area, the percentage of wild stray salmon is relatively low and they only occasionally contribute to spawning nonnative rivers. However due to stocking activity and especially to the release of juveniles in river-mouth areas, the percentage of strayers can increase by up to dozens of percent .
Our findings demonstrate that preserved fish scales can be a very valuable source of genomic DNA in population genetic studies. However, analyses based on markers from non-coding regions such as microsatellites often lead to a large amount of allele dropout , whereas with SNPs this phenomenon is much less common  although archival samples can be affected by a reduction in the number of heterozygous alleles. Johnston et al.  reported that genotyping DNA from preserved Atlantic salmon scales was highly successful for samples up to 24 years old. Although the degree of DNA degradation increases with age, we were able to successfully genotype DNA from scales that were older than 50 years. Furthermore, as Johnston et al.  indicated, we observed that, with such samples, incorrectly aligned reads occur mainly in the case of multi-site variants, thus analysis of most of this type of genetic variation was not possible.
In our study, the lowest level of heterozygosity, most departures from HWE and the lowest allelic richness were found for the salmon individuals that were sampled from the Vistula River. These results can be explained by a reduction in population size that was inferred from a drop in commercial and research catches prior to complete extinction in the 1960 s .
Säisä et al.  suggested that the southern Baltic Sea was colonized by salmon from a southern refugium during the Baltic Ice Lake stage since deglaciation occurred very early in the Baltic Sea history. Pre-glacial ice-dammed lakes are also known to have existed in the Neman, Vistula and Oder basins and they may have served as glacial salmon refugia [16, 17]. Our results partially support this theory. Among the rivers studied here, the Oder River was the first deglaciated area about 14.5 ka BP, followed by the Vistula and Neman Rivers about 13.8 to 13.2 ka BP, and finally the Morrum and Daugava rivers between 12.8 and 12.5 ka BP . The Bayesian analysis showed that the populations from the Oder, Vistula and Neman Rivers were separated by relatively large genetic distances, in particular, between those from the Vistula and Oder Rivers, which supports the hypothesis that they originated from geographically close but separated refugia.
Our data showed that the historical salmon population from the Oder River was genetically closer to salmon from the Neman River than from the Vistula River. Salmon from the Vistula River were clearly separated from all other analyzed stocks. Although a mitochondrial DNA study reported that salmon populations from the Neman River clustered with populations from the eastern Baltic Sea including the Daugava River , our data did not confirm this result, which may be due to mitochondrial and nuclear markers having different discriminating potential .
Results from the STRUCTURE analysis provided evidence of population substructures in the samples from the Vistula and Oder Rivers, which originate from the presence of migrant fish since among the samples from both of these extinct populations, potential migrant fish from the Neman River were detected (PCoA and assignment test). This may be explained by the presence of strayers from the Neman River that could have been caught in Polish rivers or the effect of stocking activities. In our analysis, the salmon population from the Morrum River and the Polish present-day line introduced from the Daugava River stock each formed a cluster, which indicates that the post-glacial origin of the Atlantic salmon from southern Sweden and the origins of the Polish populations from the Vistula and Oder Rivers differ. Since the formation of the southern Baltic rivers occurred much earlier than that of rivers from the Baltic Main Basin, this may explain the differences among Polish stocks and populations from the south Scandinavian Peninsula.
Native Atlantic salmon populations from the southern Baltic Sea originated from diverse sources after the retreat of the glacial ice sheet. Although the polish historical salmon populations from the Oder and Vistula Rivers clustered together, they were separated by a significant genetic distance. The population from the Oder River was more closely related to the present-day Neman River stock than that from the Vistula River and genetic substructure was revealed in these historical populations. The Polish native populations had a different origin and clustered separately from the Morrum and Daugava River salmon.
RB collected fish samples. RB and APK conducted the molecular work. RB and PD participated in design of the study. RW conceived the study, and participated in its design and coordination. RB and APK analyzed the data. RB, RW and PD wrote the manuscript. All authors read and approved the final manuscript.
This study was partially funded by the Ministry of Agricultural and Rural Development Board, Poland, operating under the “Stocking of Polish Marine Areas” program, statutory topic S-25 in the IFI Olsztyn and by statutory topic IV.1 in the IO PAS. Bayesian computations were performed on clusters from PL-GRID framework. Final Illumina scoring for historical samples was carried out in the Laboratory of High Throughput Technologies, Institute of Molecular Biology and Biotechnology, Adam Mickiewicz University in Poznan.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Parrish DL, Behnke RJ, Gephard SR, McCormick SD, Reeves GH. Why aren’t there more Atlantic salmon (Salmo salar)? Can J Fish Aquat Sci. 1998;55:281–7.View ArticleGoogle Scholar
- Mannerla M, Andersson M, Birzaks J, Debowski P, Degerman E, Huhmarniemi A, et al. Salmon and sea trout populations and rivers in the Baltic Sea: HELCOM assessment of salmon (Salmo salar) and sea trout (Salmo trutta) populations and habitats in rivers flowing to the Baltic Sea. Baltic Sea Environ Proc. 2011;126A:1–80.Google Scholar
- Bartel R. The restoration of Atlantic salmon (Salmo salar L.) in Poland. Arch Pol Fish. 2001;9:219–28.Google Scholar
- Popović D. Genetic polymorphism of populations of four anadromous fish species restituted in polish waters. PhD Thesis, Warsaw University. 2008.Google Scholar
- Kazakov RV, Titov SF. Geographical patterns in the population genetics of Atlantic salmon, Salmo salar L., on U.S.S.R. territory, as evidence for colonization routes. J Fish Biol. 1991;39:1–6.View ArticleGoogle Scholar
- Nilsson J, Gross R, Asplund T, Dove O, Jansson H, Kelloniemi J, et al. Matrilinear phylogeography of Atlantic salmon (Salmo salar L.) in Europe and postglacial colonization of the Baltic Sea area. Mol Ecol. 2001;10:89–102.View ArticlePubMedGoogle Scholar
- Verspoor E, McCarthy EM, Knox D, Bourke EA, Cross TF. The phylogeography of European Atlantic salmon (Salmo salar L.) based on RFLP analysis of the ND1/16sRNA region of the mtDNA. Biol J Linnean Soc. 1999;68:129–46.View ArticleGoogle Scholar
- Koljonen ML, Jansson H, Paaver T, Vasin O, Koskiniemi J. Phylogeographic lineages and differentiation pattern of Atlantic salmon (Salmo salar) in the Baltic Sea with management implications. Can J Fish Aquat Sci. 1999;56:1766–80.View ArticleGoogle Scholar
- Säisä M, Koljonen ML, Gross R, Nilsson J, Tähtinen J, Koskiniemi J, Vasemägi A. Population genetic structure and postglacial colonization of Atlantic salmon (Salmo salar) in the Baltic Sea area based on microsatellite DNA variation. Can J Fish Aquat Sci. 2005;62:1887–904.View ArticleGoogle Scholar
- Leliūna E, Virbickas J. Phylogeographic characteristics of the Atlantic salmon (Salmo salar L) population of the Nemunas river. Acta Zool Litu. 2006;16:229–34.View ArticleGoogle Scholar
- Poćwierz-Kotus A, Bernaś R, Kent MP, Lien S, Leliűna E, Dębowski P, et al. Restitution and genetic differentiation of salmon populations in the southern Baltic genotyped with the Atlantic salmon 7 K SNP array. Genet Sel Evol. 2015;47:39.View ArticlePubMedPubMed CentralGoogle Scholar
- ICES. Report of the Baltic salmon and trout assessment working group (WGBAST): 3–12 April 2013. Tallinn. ICES CM 2013/ACOM:08. 2013:334 pp.Google Scholar
- Nielsen EE, Hansen MM, Loeschcke V. Analysis of applications DNA from old scale samples: technical aspects, and perspectives for conservation. Hereditas. 1999;130:265–76.View ArticleGoogle Scholar
- Speller CF, Hauser L, Lepofsky D, Moore J, Rodrigues AT, Moss ML, et al. High potential for using DNA from ancient herring bones to inform modern fisheries management and conservation. PLoS One. 2012;7:e51122.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnston SE, Lindqvist M, Niemelä E, Orell P, Erkinaro J, Kent MP, et al. Fish scales and SNP chips: SNP genotyping and allele frequency estimation in individual and pooled DNA from historical samples of Atlantic salmon (Salmo salar). BMC Genomics. 2013;14:439.View ArticlePubMedPubMed CentralGoogle Scholar
- Marks L. Last glacial maximum in Poland. Quat Sci Rev. 2002;21:103–10.View ArticleGoogle Scholar
- Tonteri A, Titov S, Veselov A, Zubchenko A, Koskinen MT, Lesbarreres D, et al. Phylogeography of anadromous and non-anadromous Atlantic salmon (Salmo salar) from northern Europe. Ann Zool Fennici. 2005;42:1–22.Google Scholar
- Uścinowicz S. Southern Baltic area during the last deglaciation. Geol Q. 1999;43:137–48.Google Scholar
- Crespi BJ, Fulton MJ. Molecular systematics of Salmonidae: combined nuclear data yields a robust phylogeny. Mol Phylogenet Evol. 2004;31:658–79.View ArticlePubMedGoogle Scholar
- Bourret V, Kent MP, Primmer CR, Vasemägi A, Karlsson S, Hindar K, et al. SNP-array reveals genome-wide patterns of geographical and potential adaptive divergence across the natural range of Atlantic salmon (Salmo salar). Mol Ecol. 2013;22:532–51.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2007;1:47–50.PubMedPubMed CentralGoogle Scholar
- Kalinowski ST. HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes. 2005;5:187–9.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23:1801–6.View ArticlePubMedGoogle Scholar
- Peakall R, Smouse P. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Piry S, Alapetite A, Cornuet J-M, Paetkau D, Baudouin L, Estoup A. GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered. 2004;95:536–9.View ArticlePubMedGoogle Scholar
- Paetkau D, Calvert W, Stirling I, Strobeck C. Microsatellite analysis of population structure in Canadian polar bears. Mol Ecol. 1995;4:347–54.View ArticlePubMedGoogle Scholar