Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Matti T. Heino, Tommi Nyman, Jukka U. Palo, Jenni Harmoinen, Mia Valtonen, Małgorzata Pilot, Sanni Översti, Elina Salmela, Mervi Kunnasranta, Risto Väinölä, A. Rus Hoelzel & Jouni Aspi Title: Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections Year: 2023 Version: Published version Copyright: The Author(s) 2023 Rights: CC BY 4.0 Rights url: http://creativecommons.org/licenses/by/4.0/ Please cite the original version: Heino, M. T., Nyman, T., Palo, J. U., Harmoinen, J., Valtonen, M., Pilot, M., Översti, S., Salmela, E., Kunnasranta, M., Väinölä, R., Hoelzel, A. R., & Aspi, J. (2023). Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections. Ecology and Evolution, 13, e9720. https://doi.org/10.1002/ece3.9720 Ecology and Evolution. 2023;13:e9720.  | 1 of 16 https://doi.org/10.1002/ece3.9720 www.ecolevol.org Received: 28 July 2022  | Revised: 14 December 2022  | Accepted: 19 December 2022 DOI: 10.1002/ece3.9720 R E S E A R C H A R T I C L E Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections Matti T. Heino1,2  | Tommi Nyman3  | Jukka U. Palo2,4  | Jenni Harmoinen1,5  | Mia Valtonen5,6,7  | Małgorzata Pilot8,9,10  | Sanni Översti11,12  | Elina Salmela12,13  | Mervi Kunnasranta14,15  | Risto Väinölä16  | A. Rus Hoelzel17  | Jouni Aspi1 1Ecology and Genetics Research Unit, University of Oulu, Oulu, Finland 2Department of Forensic Medicine, University of Helsinki, Helsinki, Finland 3Department of Ecosystems in the Barents Region, Svanhovd Research Station, Norwegian Institute of Bioeconomy Research, Svanvik, Norway 4Forensic Chemistry Unit/Forensic Genetics, Finnish Institute for Health and Welfare, Helsinki, Finland 5Wildlife Ecology Group, Natural Resources Institute Finland, Helsinki, Finland 6Department of Environmental and Biological Sciences, University of Eastern Finland, Joensuu, Finland 7Institute of Biotechnology, University of Helsinki, Helsinki, Finland 8School of Biological and Biomedical Sciences, Durham University, Durham, UK 9Museum and Institute of Zoology, Polish Academy of Sciences, Gdańsk, Poland 10Faculty of Biology, University of Gdańsk, Gdańsk, Poland 11Transmission, Infection, Diversification and Evolution Group, Max- Planck Institute for the Science of Human History, Jena, Germany 12Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland 13Department of Biology, Faculty of Science, University of Turku, Turku, Finland 14University of Eastern Finland, Joensuu, Finland 15Natural Resources Institute Finland, Joensuu, Finland 16Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland 17Department of Biosciences, Durham University, Durham, UK This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. © 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. Correspondence Matti T. Heino, Department of Forensic Medicine, University of Helsinki, Helsinki, Finland. Email: matti.heino@helsinki.fi Funding information Ella and Georg Ehrnrooth Foundation; Emil Aaltonen Foundation; Marie Skłodowska Curie Intra- European Fellowship from the European Commission, Grant/Award Number: PIEF- GA- 2009- 235978; Jane and Aatos Erkko Foundation; Jenny and Antti Wihuri Foundation; Kuopio Naturalists' Society; Maj and Tor Nessling Foundation; Polish National Agency for Academic Exchange (NAWA); Polish Returns Fellowship, Grant/Award Number: PPN/ PPO/2018/1/00037; Nestori Foundation; University of Oulu Scholarship Foundation; Raija and Ossi Tuuliainen Foundation Abstract The Saimaa ringed seal (Pusa hispida saimensis) is endemic to Lake Saimaa in Finland. The subspecies is thought to have originated when parts of the ringed seal popula- tion of the Baltic region were trapped in lakes emerging due to postglacial bedrock rebound around 9000 years ago. During the 20th century, the population experi- enced a drastic human- induced bottleneck. Today encompassing a little over 400 seals with extremely low genetic diversity, it is classified as endangered. We se- quenced sections of the mitochondrial control region from 60 up to 125- years- old museum specimens of the Saimaa ringed seal. The generated dataset was combined with publicly available sequences. We studied how genetic variation has changed through time in this subspecies and how it is phylogenetically related to other ringed seal populations from the Baltic Sea, Lake Ladoga, North America, Svalbard, and the White Sea. We observed temporal fluctuations in haplotype frequencies and loss of haplotypes accompanied by a recent reduction in female effective population 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 2 of 16  |     HEINO et al. 1  |  INTRODUC TION Modern DNA sequencing technologies and analytical methods, especially simulation- based approaches, enable detailed molecu- lar inferences on population- level processes and demographic his- tory in human and wildlife populations (Foote et al., 2021; Hoban et al., 2012; Taylor et al., 2021; Vilaça et al., 2021). However, there are limits as to what can be deduced from present- day genetic data. This is largely since the contemporary genetic composition of any popula- tion has been shaped by the actions of all evolutionary forces (drift, mutation, migration, and selection) throughout the whole history of the population. Understanding the population dynamics at a certain period would require disentangling the effects of these forces, which is often impossible. The confounding signals arising from the actions of the different forces at different times can be to some extent min- imized, for example, by marker choice (ignoring selection, assuming steady mutation rates), but drift and gene flow are harder to partition. These limitations of contemporary variation can nevertheless be circumvented by directly obtaining data on past diversity. Temporal sampling of varied genetic data, for example, from subfossil, ar- cheological, or museum specimens, can reveal the developments leading to the palimpsest of current genetic diversity and provide a substantially more accurate view of past evolutionary history. In this context, zoological museum collections have proven an invalu- able source of direct knowledge on temporal changes in DNA diver- sity and, ultimately, on the demography of various wildlife species (Hofreiter et al., 2015; Huynen et al., 2012; Nakahama, 2021; Rizzi et al., 2012). The Saimaa ringed seal (Pusa hispida saimensis) offers an excel- lent model population for exploring the effects of small size and isolation on genetic diversity. Endemic to Lake Saimaa in Finland (Figure 1), its origin, and timing of isolation are deemed as well es- tablished. It is considered a subspecies of the Holarctic ringed seal, which, numbering in millions, is the most abundant Arctic pinniped (Reeves, 1998). Based on the understanding of the geological his- tory of the Baltic Sea region, the population originates from marine ringed seals that colonized the Baltic basin during the deglaciation of the Scandinavian Ice Sheet, c. 10,000 years ago (Ukkonen, 2002; size. In apparent contrast with the traditionally held view of the Baltic origin of the population, the Saimaa ringed seal mtDNA variation also shows affinities to North American ringed seals. Our results suggest that the Saimaa ringed seal has experi- enced recent genetic drift associated with small population size. The results further suggest that extant Baltic ringed seal is not representative of the ancestral popula- tion of the Saimaa ringed seal, which calls for re- evaluation of the deep history of this subspecies. K E Y W O R D S freshwater pinniped, genetic diversity, genetic drift, mitochondrial DNA, museum specimens, Saimaa ringed seal T A X O N O M Y C L A S S I F I C A T I O N Biogeography, Conservation genetics, Demography, Genetics, Phylogenetics, Population genetics F I G U R E 1 Maps showing the areas from which (a) marine and (b) freshwater ringed seals were sampled for this study (note that a permanent population is not present in Lake Pielinen). (a) (b) Lake Saimaa Lake Ladoga Lake Pielinen Baltic Sea White Sea Svalbard North America 70°N 60°N 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  3 of 16HEINO et al. Ukkonen et al., 2014). As a consequence of postglacial bedrock rebound, parts of the Baltic ringed seal population (currently P. h. botnica) were trapped in emerging freshwater lakes, including lakes Saimaa and Ladoga (subspecies P. h. ladogensis). Hence, the Saimaa ringed seal has lived in complete isolation for c. 860 generations (Palo et al., 2003), during which it has evolved into a morpholog- ically, behaviorally, and genetically distinct subspecies (Berta & Churchill, 2012; Hyvärinen & Nieminen, 1990; Kunnasranta, 2001). Today, the Saimaa ringed seal population encompasses c. 420– 430 seals (Metsähallitus, 2020) and is classified as endangered, both nationally (Liukko et al., 2016) and internationally (Kunnasranta et al., 2021). During the 20th century, the population experienced a drastic human- induced bottleneck from hunting and later due to environmental pollutants. Due to its limited lacustrine habitat, the population has proba- bly never been large; based on the size and resources of the lake, the maximum number of seals has been estimated as being few thousands (Hyvärinen et al., 1999). In the late 1800s, there were possibly as many as 1000 seals remaining in Lake Saimaa, as sug- gested by bounty statistics (Kokko et al., 1999) and local knowledge (Hyvärinen et al., 2004). The dwindling of the population started in the early 1900s and continued through most of the century de- spite the protection of the subspecies in 1955. In the 1980s, the population reached its ultimate low of fewer than 150 seals (Sipilä et al., 1990). Since then, intensified conservation efforts (aiming at N > 400) set the population to a slow increase in the 1990s, and the positive trend has now continued for more than 20 years (Kunnasranta et al., 2021). In comparison with other ringed seal subspecies— and almost any mammal population— genetic diversity of the Saimaa ringed seal is extremely low (Martinez- Bakker et al., 2013; Palo et al., 2003; Valtonen et al., 2012, 2014). Roughly 55% of the autosomal hetero- zygosity and 90% of the mitochondrial nucleotide diversity present today in the marine population has been lost during the founder event and isolation. The diversity is also substantially lower than in Ladoga ringed seals (Nyman et al., 2014; Valtonen et al., 2012, 2014). Although the Saimaa ringed seal is known to have suffered a substantial reduction in numbers during the 20th century, it is not known how much of the total reduction of genetic diversity can be attributed to these recent events. Previous genetic studies have es- timated a population trajectory spanning back to the isolation time- point (Nyman et al., 2014; Palo et al., 2003; Valtonen et al., 2012) but have suggested a minor role for the 20th century bottleneck in the total diversity loss. Also, Nyman et al. (2014), exploring the plau- sible demographic histories leading to the current diversity using Approximate Bayesian Computation (ABC) analyses, could not find conclusive support for models including the recent bottleneck, over models including only an ancient bottleneck associated with the col- onization of Lake Saimaa. Using mitochondrial DNA (mtDNA) control- region sequences obtained from 60 up to 125- years- old museum specimens, as well as previously published contemporary data, we explored changes in Saimaa ringed seal mtDNA diversity. Our main aim was to evaluate whether the loss of genetic variation of the Saimaa ringed seal can be primarily attributed to slow erosion during its isolation after a founder event, as previously suggested (Nyman et al., 2014; Palo et al., 2003; Valtonen et al., 2012), or to the 20th- century bottleneck (Peart et al., 2020; Stoffel et al., 2018). We also evaluate the history of the Saimaa population by positioning the diversity in Saimaa on the cladistics canvas of mtDNA variation found in ringed seals from the Baltic Sea, Lake Ladoga, North America (Chukchi and Beaufort Sea), Svalbard, and the White Sea. 2  |  MATERIAL S AND METHODS 2.1  |  DNA extraction, amplification, and sequencing We obtained tissue samples (bones, claws, or skin) of 63 historical Saimaa ringed seals from the collections of the Finnish Museum of Natural History (Table S1). DNA extractions, amplifications, and se- quencing were performed either at the University of Oulu, Finland, or at Durham University, UK. At the University of Oulu, DNA was extracted in the ancient- DNA laboratory of the Centre for Material Analysis, following the protocol by Rohland et al. (2010). For bones and claws, a small por- tion of the surface of each sample was removed, and circa 250 mg of powder was obtained by drilling into the sample; the powder was then subjected to DNA extraction. For hides, a circa 250- mg piece was cut into small pieces with a sterile scalpel and then sub- jected to extraction. We targeted a circa 704- bp section of the mitochondrial control region (with some length variation observed among the samples). The whole region was amplified and se- quenced in short variable- sized overlapping fragments depending on sample quality (for primers, see Table S2). We first attempted to amplify each sample using primer pairs H16305/R328, F322/ R471, and F425/L224. If this failed, amplification of shorter frag- ments was performed using primer pairs H16305/R109, F108/ R328, F322/R471, F425/R605, and F603/L224. Each PCR in- cluded 1X AmpliTaq Gold® 360 Buffer, 2 mM of MgCl2, 1 mg/ml of BSA, 200 μM of each dNTP, 0.4 μM of each primer, and 0.6 U of AmpliTaq Gold® 360 DNA Polymerase. The PCR profile con- sisted of initial denaturation in 95°C for 10 min, then 50 cycles of 95°C for 30 s, 50°C for 30 s, and 72°C for 30 s, and finally 72°C for 7 min. To identify possible mis- incorporated bases caused by postmortem changes, we PCR- amplified and sequenced every re- gion at least twice. As a control, some samples were also extracted twice. Negative controls were included in each extraction and PCR batch. PCR products were purified using Exonuclease I (Thermo Fisher Scientific) and Shrimp Alkaline Phosphatase (Thermo Fisher Scientific). Sequencing reactions were then done using BigDye Terminator v1.1 Cycle Sequencing Kit (Thermo Fisher Scientific), and the products were run on ABI 3730 DNA Analyzer (Applied Biosystems). DNA sequences were edited and assembled into con- tigs using CodonCode Aligner v4.0.4 (CodonCode Corporation). 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 4 of 16  |     HEINO et al. We also included previously unpublished sequence data from addi- tional modern samples from Lake Saimaa (N = 5) and the Baltic Sea (N = 2) (for methods, see Valtonen et al. (2012)). At the Durham University, DNA extraction was carried out fol- lowing the method described in de Bruyn et al. (2009). The mtDNA control region was amplified in three overlapping fragments of 185– 210 bp, from which 490 bp sequence was reconstructed (for primers, see Table S2). PCRs were performed in a volume of 16 μl containing 1X Qiagen Multiplex PCR Master Mix, 0.2 μM of each primer, and 5 μl of the five times diluted extract. The thermal cycling profile was as follows: 95°C for 15 min, 45 cycles of 94°C for 30 s, 50°C for 90 s, and 72°C for 60 s, followed by a final extension at 72°C for 10 min. Negative controls were used for all PCR runs. PCR products were purified using QIAquick PCR Purification Kits (Qiagen) according to the manufacturer's instructions. Purified PCR products were se- quenced at the sequencing facility of Durham University. As in Oulu, each sample was sequenced in both directions, and the chromato- grams of each sequence were checked by eye. We obtained the full target region of circa 704 bp from 56 out of 63 museum specimens, 49 of which were fully replicated. A partial target region was obtained from four additional speci- mens, three of which were fully replicated. DNA in ancient and historical samples is often damaged, which results in nucleotide misincorporations (Hansen et al., 2001; Hofreiter et al., 2001; Pääbo et al., 2004). All ambiguous bases that remained after am- plification and sequencing a position in a sample for two or more times were given an IUPAC ambiguity code (e.g., a position that could be either a cytosine (C) or a thymine (T) was coded as Y). We included samples that had unreplicated sequence areas in the analyses, after replacing any unique SNPs in their unreplicated parts with IUPAC ambiguity codes. For example, if the unrepli- cated part had a T in a position in which all other Saimaa ringed seals had a C, the T was changed to Y. We additionally included in the analyses five previously unreported modern sequences from Saimaa and two from the Baltic Sea (Table S1). All newly reported sequences have been deposited in GenBank under accession num- bers ON841430- ON841496. 2.2  |  Haplotype network A haplotype network for seals from the Baltic area was con- structed based on a combined dataset including sequences of museum and modern samples. For the museum samples, only the 56 samples from which the full target region was obtained were included. The modern data included sequences from 208 Saimaa, 21 Baltic, and 16 Ladoga ringed seals sampled during the years 1980– 2008, 2006– 2009, and 1991– 2000, respectively (Valtonen et al., 2012), as well as 57 sequences obtained from Saimaa ringed seal placentas collected from breeding sites during the years 2009– 2011 (Valtonen et al., 2015). Sequences were aligned using the FFT- NS- i method within the online version of MAFFT v7 (https://mafft.cbrc.jp/alignment/server/; Katoh et al., 2002, 2019; Katoh & Standley, 2013). A total of 24 bases from the begin- ning and 20 bases from the end of the alignment were excluded to accommodate the seven previously unpublished modern se- quences. The final alignment included 663 bp from 358 samples. We first converted the Fasta- formatted alignment file into a TCS input file using FaBox (Villesen, 2007) and then built a haplotype network using TCS (Clement et al., 2000), treating gaps as 5th state. The network was visualized with tcsBU (Múrias Dos Santos et al., 2016). 2.3  |  Temporal genetic variation To study temporal genetic variation in the Saimaa ringed seal popu- lation, the full- length sequences from Saimaa (N = 321) were divided into five temporal groups (TG) based on the history of the popula- tion outlined in Kunnasranta et al. (2021): TG1 pre- bottleneck 1894– 1939 (N = 23), TG2 bottleneck 1960– 1979 (N = 33), TG3 bottleneck 1980– 1989 (N = 81), TG4 recovery 1990– 1999 (N = 56), and TG5 increase 2000– 2011 (N = 128). We decided to treat the latter group as modern as the largest changes in population size have happened before this period. Temporal changes in the haplotype frequencies in the Saimaa ringed seal population were visualized as a haplotype network di- vided into the aforementioned five temporal groups. The figure was drawn by manually editing networks built with TCS and tcsBU based on templates constructed in TempNet (Prost & Anderson, 2011). As there was large variation in samples sizes from different temporal groups, we also built a temporal network with circle sizes reflect- ing the relative proportion of samples belonging to each haplotype. We further explored how sampling bias could affect the estimates of haplotype diversity by performing sample size based rarefaction for the haplotypes observed in each temporal period using iNEXT (iNterpolation and EXTrapolation) Online (Chao et al., 2014, 2016) using 200 bootstraps. As the sample sizes per temporal group are relatively small, the significance of the haplotype frequency differences between the contemporary population (TG5) and earlier temporal groups (TG1– TG4) was tested using a resampling approach. For each haplo- type in temporal groups TG1– TG4, N observations were drawn with a probability of success equal to the frequency of that particular haplotype in the reference population (TG5) and N equal to the sam- ple size in the temporal group in question. For the six haplotypes ab- sent from the reference population TG5, resampling was performed using a frequency of 0.023132 instead of zero, as a haplotype of that frequency would have a 5% chance of not showing up in a sample of N = 128. The resampling was repeated 100,000 times per hap- lotype and temporal group, and the resultant empirical distribution of the expected number of occurrences of each haplotype was re- corded and compared with the observed number. The underlying idea is that a nonsignificant haplotype frequency difference may be due to sampling variation, whereas a significant difference could be caused by other factors, such as genetic drift or natural selection. 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  5 of 16HEINO et al. The resampling was performed with a script written for R v3.6.0. by ES (https://github.com/esalm ela/Conti Unity). We tested for the existence of differences in the spatial dis- tribution of samples across the different temporal groups using chi- square tests. For this, we first excluded the individual from Lake Pielinen, currently uninhabited by seals, as well as nine his- torical specimens for which the collection site within Lake Saimaa was unknown. To improve concordance to the test's assumptions, we pooled time periods TG1 and TG2 and combined Kolovesi and Haukivesi into a single area, resulting in a contingency table of four time periods by four areas. As over 20% of the expected fre- quencies were still below 5, we conducted a second test using a five time periods by two areas dataset, in which sampling loca- tions were coded as “northern” (Northern Saimaa, Kolovesi, and Haukivesi basins) or “southern” (Pihlajavesi and Southern Saimaa basins). 2.4  |  Past effective female population size We studied changes in the past effective female population size of the Saimaa ringed seal by constructing Bayesian skyline plots (Drummond et al., 2005) in BEAST v1.10.4 (Suchard et al., 2018) based on the 321 full- length sequences. The best- fitting substi- tution model was inferred using model averaging with bModel- Test (Bouckaert & Drummond, 2017) in BEAST v2.6.3 (Bouckaert et al., 2019) and analyzing the results with Tracer v.1.7.1 (Rambaut et al., 2018). The resulting optimal substitution model was TN93 (Tamura & Nei, 1993) accounting for proportion of invariant sites and rate variation across sites modeled with four gamma categories (Yang, 1994). The data were fitted into two molecular clock mod- els: strict and relaxed uncorrelated lognormal clock. The most suit- able clock model was determined with Akaike's information criteria through MCMC (AICM) (Baele et al., 2012) implemented in Tracer v1.6 (Rambaut et al., 2018). A strict- clock model was preferred over a relaxed clock (AICM = 11.86 with 100 bootstrap replicates) and was used in the subsequent analyses. The sampling years of the sequences were used as a source for external calibration of the molecular clock (i.e., “tip- calibration”). We performed three paral- lel analyses, each with 500 million steps. Every 50,000th step was sampled, and the first 10% of the steps were discarded as a burn- in. Consistency between runs was evaluated in Tracer, and chains were combined with LogCombiner included in the BEAST software package. Tracer v1.7.2 (Rambaut et al., 2018) was used to con- firm that an adequate effective sample size (ESS > 200) had been reached for each parameter and for visualization of the Bayesian skyline plot. Furthermore, although it has been shown that tip- calibration yields more consistent results than internal node calibration (Rieux et al., 2014), on an evolutionary scale, samples in this study can be considered recent. Thus, to assess if the temporal structure of the data used is sufficient to calibrate the molecular clock, we conducted a date- randomization test. We performed 20 random shufflings of tip dates with the R package TIPDATINGBEAST v1.1– 0 (Rieux & Khatchikian, 2017) and ran the data- randomized datasets with BEAST v1.10.4 as described above. 2.5  |  Phylogenetic relationships between ringed seal populations Finally, we studied broader phylogenetic connections among ringed seals from Lake Saimaa and other representative subspe- cies and populations in the Arctic and Baltic regions. As an addi- tion to the data described above, sequences from Palo et al. (2003) and Martinez- Bakker et al. (2013) were included in this analysis. As for the Baltic region dataset mentioned above, sequences were aligned using the FFT- NS- i method in MAFFT. We used only the se- quence region that overlapped in all the datasets, which resulted in an alignment of 348 bp. Next, sequences that had any ambiguous bases in polymorphic sites were removed, bringing the number of sequences down to 551. In order to make the data more manageable for the analysis, sequences were collapsed into 211 unique haplo- types using the DNA to haplotype collapser and converter in FaBox (Villesen, 2007). A Bayesian phylogenetic analysis of the sequence dataset was run in MrBayes v3.2.7a (Ronquist et al., 2012) on the CIPRES Server (Miller et al., 2010). To retain information included in alignment gaps, gaps were recoded to binary presence/absence characters following the method of Simmons and Ochoterena (2000) as implemented in FastGap v1.2 (Borchsenius, 2009). In the run, we used a mixed substitution model with gamma- categorized rate vari- ation for nucleotide sites, while the recoded gaps were treated as a separate partition of variable, gamma- corrected “restriction” sites. Ratepr, statefreqpr, shape, revmat, and tratio were left as default and unlinked across partitions, but the branch length prior was set to unconstrained:Exp(100.0). Two runs with four incrementally heated chains each were run for 10 million generations, while sampling trees every 1000 generations. Run convergence was assessed in Tracer and, after removing the first 10% of the sampled trees as a burn- in, a midpoint- rooted summary tree showing all compatible groupings was calculated in MrBayes. 3  |  RESULTS 3.1  |  Haplotype network Our haplotype network for the Baltic area ringed seal mtDNA vari- ation (Figure 2) shows similar haplotype assemblage as in Valtonen et al. (2012, 2015), but the museum samples carried four additional Saimaa- specific haplotypes (H9, H10, H12, and H14) while lacking four previously observed haplotypes (H4, H6, H7, and H8). In total, the “Saimaa clade” of the network consists of 14 distinct haplotypes (H1– H14), with one haplotype from the Baltic Sea and the single se- quence from Lake Pielinen (H11) being placed in the center of this assemblage. We emphasize that Lake Pielinen has not harbored a 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 6 of 16  |     HEINO et al. permanent ringed seal population in historical times, so the sample from this lake clearly represents a migrant individual that has strayed to Lake Pielinen from Lake Saimaa via the 67 km long River Pielisjoki. As in the results of Valtonen et al. (2012), the network shows all the Lake Saimaa haplotypes clustering relatively close to each other, with two haplotypes (H1 and H3) dominating through time. Haplotype H14 encountered in Lake Saimaa has a 64- bp deletion but also ad- ditional SNPs, which explain its outlier position in the network. 3.2  |  Temporal genetic variation Division of the samples into five temporal groups (TG1– TG5) al- lowed exploring changes through time in the Lake Saimaa mtDNA haplotype pool. The temporal haplotype networks (Figures 3, S1) and the haplotype table (Table 1) show substantial fluctuation in haplotype frequencies through time. The most common haplo- type today, H3, seems to have been dominant also during the two first periods (TG1: 1894– 1939 and TG2: 1960– 1979). The second most common extant haplotype, H1, however, seems to have be- come common only after TG1: 1894– 1939, and it was even more common than H3 during TG3: 1980– 1989 and TG4: 1990– 1999. The third most common haplotype (H2) is observed only from pe- riod TG2: 1960– 1979 onwards, with a slight increase in frequency towards the present. Despite the relatively small temporal group- specific sample sizes, our resampling tests (Figure 4) showed that many of the observed dif- ferences in haplotype frequencies are unlikely to be caused by mere sampling variation: Five haplotypes in TG1– TG4 showed a nominally significant difference from TG5 within at least one time period (seven significant differences in all). While this result is not corrected for mul- tiple testing (i.e., for the 4 × 14 = 56 haplotype- wise tests), only 2.8 significant results would be expected by chance at a significance level of α = 0.05 in 56 independent tests, and the probability of obtaining seven or more significant results by chance would be 0.021. Taking into account also the Lake Pielinen specimen, eight dif- ferent haplotypes were observed in the earliest time interval TG1 (1894– 1939). Three of these (H10, H11, and H12) were not observed in later time. The resampling analysis suggested that the overrepre- sentation of H9 and H12 in this time interval cannot be explained by small sample size alone (Figure 4). H12 is only observed in TG1, with a frequency of ≈13%. The fact that three unique haplotypes F I G U R E 2 Haplotype network of ringed seals in the Baltic Sea and lakes Saimaa, Pielinen, and Ladoga. The sizes of the circles correspond to the number of sampled individuals with the haplotype, section colors denote sampling site and collection period within Lake Saimaa (see map and legend in inset). Small white circles denote unobserved haplotypes. Note that haplotype H14 is present in a single individual that has a 64- bp deletion (each bp in the deletion is counted as one change) as well as unique substitutions. Lake Saimaa 1980–2011 Lake Saimaa 1894–1976 Lake Pielinen Spatial and temporal groups Lake Ladoga Baltic Sea H14 H8 H1 H5 H10 H11 H12 H4 H6 H2 H3 H7 H13 H9 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  7 of 16HEINO et al. were observed only in the smallest sample of the oldest time period shows that some variation has been permanently lost and also sug- gests that some of the past variation remains unsampled. These con- clusions are also supported by the period- specific rarefaction curves, although the 95% confidence intervals overlap broadly across periods (Figure S2). In TG2 (1960– 1979), only six haplotypes were observed despite a larger sample size than in TG1, suggesting substantial loss of haplotype diversity between 1894 and 1960. TG3 (1980– 1989) demonstrates a step towards the modern haplotype distribution, with H1 and H3 dominating. Nevertheless, compared to modern times, H1 is nominally significantly overrepresented and H3 underrepresented both in the 1980s and 1990s samples (TG3 and TG4). Furthermore, H7, which has a frequency of ≈5% in the modern TG5 (2000– 2011), appears for the first time in TG4 (1990– 1999). Haplotype H8 is only observed in TG5 (2000– 2011) despite the larger summed sample size of the earlier time periods (N = 128 vs. N = 193). Our tests for differential sampling across time periods indi- cated differences in the four temporal groups by four areas dataset (χ2 = 12.7, df = 9, P = 0.0005), but in this case the assumptions of the test are slightly violated (over 20% of the expected frequencies were below 5). The test assumptions were, however, met by the five tem- poral groups by two areas dataset, which also resulted in statistically significant differences in sampling (χ2 = 12.744, df = 4, P = 0.013). 3.3  |  Past effective female population size The Bayesian skyline plot (Figure 5) suggests that there was a slow decline in the female effective population size (Nef) of Saimaa ringed seals from the 12th century. Following this, there was a sharp de- crease in Nef from the middle of the 19th century onwards. In our date- randomization test, the 95% highest posterior densities for the Clock.rate and TreeModel.rootHeight parameters overlapped partially between the actual and randomized datasets (Figure S3). F I G U R E 3 Haplotype network of Saimaa ringed seals during five different time intervals spanning periods from 1894 to the present. In each network, circle size is proportional to the number of samples representing each haplotype, and sector colors denote the proportion of samples belonging to each of the sampling areas indicated in the map to the right of the networks. 2000–2011 TG5 n=128 1990–1999 TG4 n=56 1980–1989 TG3 n=81 1960–1979 TG2 n=33 1894–1939 TG1 n=23 H5 H1 H8 H14 H10 H11 H13 H3 H12 H6 H4 H2 H7 H9 (72 steps) 101 3 1 1 3 1 3 1412 1 4 1 2148 2 4 3 2 1 19 1 28 1 4 1 1 63 7 4 40 10 2 1 1 1 1 Southern Saimaa Pihlajavesi area Haukivesi area Northern Saimaa Kolovesi Lake Pielinen = Unknown TA B L E 1 Number of each haplotype observed in each temporal period Period TG1 TG2 TG3 TG4 TG5 Haplotype H1 3 12 48 28 40 H2 0 4 4 4 10 H3 10 14 21 19 63 H4 0 0 3 1 2 H5 1 0 2 1 0 H6 0 0 1 1 1 H7 0 0 0 1 7 H8 0 0 0 0 4 H9 3 1 0 0 0 H10 1 0 0 0 0 H11 1 0 0 0 0 H12 3 0 0 0 0 H13 1 1 2 1 1 H14 0 1 0 0 0 Total 23 33 81 56 128 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 8 of 16  |     HEINO et al. 3.4  |  Phylogenetic relationships among ringed seal populations In our ringed seal mtDNA phylogeny (Figure 6), the Saimaa ringed seal clade is unexpectedly nested within the Arctic, and more spe- cifically North American, ringed seal diversity rather than the Baltic ringed seal diversity. The ringed seals from the White Sea, which are geographically the closest Arctic population to Saimaa, are however not closely related to Saimaa. This is also the case for the geographically adjacent Ladoga ringed seal population, as individu- als from Lake Ladoga are widely scattered throughout the phylog- eny. Interestingly, the aforementioned Saimaa ringed seal individual 5688 (haplotype H14) with the 64- bp deletion is grouped with 1.0 posterior support with two North American individuals that share the same deletion as well as SNPs. Five Baltic ringed seal haplotypes are nested inside the Saimaa ringed seal clade. F I G U R E 4 Haplotype frequency differences between the modern population and earlier temporal groups (TG1– TG4). Dark gray bars depict the expected distribution of the number of observations per haplotype in each temporal group based on 100,000 rounds of resampling from the modern distribution (TG5, N = 128). Red bars indicate the observed numbers within each temporal group, and red font denotes haplotypes with a nominally significant frequency difference as compared the modern population (i.e., with an observed number of occurrences outside 95% of the simulated probability mass). 0 3 6 9 13 17 21 25 29 33 37 41 45 49 53 (a) TG4 (1990–1999; N = 56) H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 0 4 8 13 19 25 31 37 43 49 55 61 67 73 79 (b) TG3 (1980–1989; N = 81) H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 0 2 4 6 8 10 13 16 19 22 25 28 31 (c) TG2 (1960–1979; N = 33) H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 0 2 4 6 8 10 12 14 16 18 20 22 (d) TG1 (1894–1939; N = 23) H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 Fr eq ue nc y of h ap lo ty pe c ou nt s Fr eq ue nc y of h ap lo ty pe c ou nt s Fr eq ue nc y of h ap lo ty pe c ou nt s Fr eq ue nc y of h ap lo ty pe c ou nt s Haplotype count in simulated and observed datasets Haplotype count in simulated and observed datasets 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  9 of 16HEINO et al. 4  |  DISCUSSION Improvements in molecular- genetic methods have increasingly made museum collections an invaluable source for gaining insights into the past genetic composition of animal and plant species (Burrell et al., 2015; Nakahama, 2021; Raxworthy & Smith, 2021). Museum material allows inferences of long- term changes in the frequencies of genetic variants (Bi et al., 2019) and can be used to improve estimates of demographic trajectories (Hung et al., 2014). Especially for endan- gered species, genotyping museum samples can allow direct dem- onstration of loss of genetic diversity over time due to diminished population size (Dussex et al., 2021; van der Valk et al., 2019; von Seth et al., 2021). Here, we sequenced mitochondrial control- region sequences from 60 museum- preserved samples of the endangered Saimaa ringed seal and then combined these data with previously published contemporary sequences from Lake Saimaa. We used this combined dataset spanning 117 years to infer the genetic impacts of the 20th- century population bottleneck, caused mainly by relentless persecution of Saimaa ringed seals during the first half of the cen- tury. We also combined our Saimaa dataset with existing data from other Holarctic ringed seal subspecies, in order to investigate the origin of the lake- endemic Saimaa ringed seal population. 4.1  |  The genetic signature of the 20th- century bottleneck The Saimaa ringed seal ranks among the genetically least diverse animal subspecies on the Earth, both in terms of mitochondrial (Martinez- Bakker et al., 2013; Nyman et al., 2014; Palo, 2003; Valtonen et al., 2012, 2014) and nuclear (Martinez- Bakker et al., 2013; Nyman et al., 2014; Palo, 2003; Palo et al., 2003; Stoffel et al., 2018; Valtonen et al., 2014, 2015) genetic diversity. The genetic uniform- ity of the subspecies is even more striking when considering the fact that ringed seals in general tend to be highly diverse (Martinez- Bakker et al., 2013; O'Corry- Crowe, 2008; Olsen et al., 2011; Palo et al., 2001). As estimated by Nyman et al. (2014), Palo et al. (2003), and Valtonen et al. (2012), the Saimaa ringed seal has lost 33.9% and 89.4% of its mitochondrial haplotype and nucleotide diversity, re- spectively, as well as 55– 69% of its overall heterozygosity in relation to its relatives in the Arctic Ocean and the Baltic Sea. As such, the anthropogenic population collapse of the Saimaa ringed seal in the 20th century is well known and thoroughly docu- mented (Hyvärinen et al., 1998; Kokko et al., 1998; Ranta et al., 1996; Sipilä et al., 1990). However, whether this bottleneck fully explains the extremely low genetic diversity of the population is less evident. The Saimaa ringed seal is long- lived and has an estimated generation time of 11 years (Palo et al., 2001). Hence, Nyman et al. (2014), Palo et al. (2003), and Valtonen et al. (2012) inferred that the main por- tion of diversity loss occurred well prior to the recent bottleneck, that is, during the long postglacial isolation of the Saimaa ringed seal. Nyman et al. (2014) used Approximate Bayesian Computation infer- ence based on combined nuclear and mitochondrial data and pos- tulated that incorporating the recent bottleneck does not improve model fit in relation to a demographic model that includes only a bottleneck during colonization of the Lake Saimaa system nearly 10,000 years ago. However, it should be emphasized that abrupt or slow postglacial genetic erosion does not preclude the possibility that the recent bot- tleneck would have led to further marked loss of diversity. Indeed, our Bayesian skyline plot (Figure 5) showed a sharp reduction in female effective population size towards the recent. The com- mencement of the dip in the estimated trajectory closely parallels the recent population collapse, which started in the early 19th cen- tury and intensified during the first half of the 20th century, largely as a result of a growing human population and increased access to modern firearms in the region. James and Eyre- Walker (2020) have provided quantification of the relationship between mitochondrial diversity and the census population size. Because mitochondria are haploid and have higher mutation rate than nuclear DNA, it is ex- pected that changes in effective population size are more readily vis- ible in mitochondrial than in nuclear DNA. Our results are however also consistent with the recent results of Stoffel et al. (2018) and Peart et al. (2020), who found statistically significant signatures of the recent bottleneck based on nuclear microsatellites and ddRAD- based SNPs, respectively. The effects of the still- continuing popu- lation recovery, which started in the 1980s, were not manifested in our most recent temporal samples of 2011. The patterns we observed in the Bayesian skyline plot need to be interpreted with caution, as the 95% posterior density interval of the trajectory is wide, and population substructure and unbal- anced sampling may confound skyline analyses (Heller et al., 2013). Moreover, in our date- randomization test, the 95% highest posterior densities for the clockRate and TreeHeight parameters overlapped somewhat between the actual and randomized datasets (Figure S3). This suggests that the temporal timeframe of the sampled sequences might not be wide enough to produce entirely reliable estimates of F I G U R E 5 Bayesian skyline plot depicting the estimated historical effective population size of Saimaa ringed seal females on a logarithmic scale. Calendar years are shown on the X- axis. The thick line shows the median and the shaded area represents the 95% highest posterior density interval. 100 1000 10000 100000 Year 1100 1200 1300 1400 1500 1600 1700 1800 1900 20001000900 E ffe ct iv e po pu la tio n si ze 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 10 of 16  |     HEINO et al. ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲▲ ▲ ▲ ▲ ▲▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲▲ ▲ ▲ ▲ ▲▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ●● ● ●● ● ●● ● ● ● ● ●● ●● ■ ■■■■ ■ ■ ■ ■ ■ ■■ ■ ■ ■ ■■■ ■ ■ ■ ■■ ■■ ■■ ■ ■ ■■ ■■■■■ ■■■■ ■■ ■ ■ ■■■■ ■ ■ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ▲ ● ● ●● ● ● ● ●● ●● ● ●● ● ●● ● ● ● ▲ ▲ ▲ ● ■ ● ● ● 0.04 0.83 0.74 0.51 0.54 0.75 0.56 1.00 0.62 0.87 0.53 0.74 0.85 0.73 0.78 0.71 0.94 0.67 0.88 1.00 0.51 0.99 0.5 0.95 0.93 0.97 0.73 0.82 0.56 0.72 1.00 0.57 0.96 1.00 0.58 0.89 1.00 0.86 1.00 0.93 0.51 1.00 0.86 0.57 0.79 0.5 0.63 0.84 1.00 0.58 0.90 0.96 0.67 1.00 0.59 0.73 0.74 0.71 0.6 0.57 0.75 0.5 1.00 1.00 0.78 0.82 0.90 0.6 0.94 0.93 0.5 0.7 0.6 0.71 0.68 0.82 0.58 North America Svalbard White Sea Baltic Sea Lake Saimaa Lake Ladoga Lake Pielinen Baltic Sea White Sea Svalbard North America Lake Saimaa Lake Pielinen Lake Ladoga H14 (KN-5688) 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  11 of 16HEINO et al. evolutionary rate and time. Future analyses should aim at strength- ening the temporal signal by including older samples, longer se- quences, and/or internal node calibrations, as suggested by Rieux and Khatchikian (2017). 4.2  |  Fluctuations in haplotype frequencies and recent loss of diversity We observed substantial fluctuation in mtDNA haplotype fre- quencies through time in the Saimaa ringed seal population. While Valtonen et al. (2012, 2014) demonstrated significant differentia- tion in haplotype frequencies between the years 1980 and 2008, our analyses incorporating museum material extend the pattern nearly a 100 years back in time. As shown by our temporal networks (Figures 3, S1), rarefaction curves (Figure S2), and resampling analy- ses (Figure 4), several haplotypes have probably been completely lost, while some haplotypes that are dominant today have only recently increased in frequency. The observed differences in hap- lotype frequencies in the different temporal groups are likely real and not caused by sampling bias. The observation of three unique haplotypes in the sample of 23 individuals that lived between 1894 and 1939, but not found among the nearly 300 individuals from later times, provides concrete evidence for diversity loss. Given the small sample sizes in the oldest temporal groups, it is also likely that much of the past variation remains unsampled. The observed haplotype loss and frequency changes suggest that, during the last 125 years, the reduction in effective female population size in the Saimaa ringed seal has been dramatic. In this, the fate of the Saimaa ringed seal resembles another severely persecuted mammal, the gray wolf (Canis lupus). The wolves in the past have harbored high mtDNA di- versity that has collapsed in the turn of the 20th century especially in Western Europe (Dufresnes et al., 2018). When interpreting the results above, a potential complication arises because of the previously documented genetic substructure within Lake Saimaa (Valtonen et al., 2014), which could bias infer- ences if the main regions of the lake are differentially represented in the temporal samples. Statistically significant differential sampling indeed seemed to be present in the dataset. The fact that samples from Southern Saimaa are only present in the three most recent temporal groups (TG3, TG4, and TG5) might explain the absence of some Southern Saimaa- specific haplotypes in the earlier time inter- vals and would at the same time bias estimated diversity of the most recent time periods upwards. However, regionally biased sampling does not seem to explain the loss of haplotypes H9, H10, H12, and H14 in the most recent time intervals, as these haplotypes were ob- served in lake regions that are well represented in temporal groups TG3– TG5. 4.3  |  Origin of the Saimaa ringed seal According to the commonly held view, the Saimaa ringed seal popu- lation originated from Baltic ringed seals trapped in an inland lake system that gradually emerged due to postglacial land uplift around 9000 years ago (Hyvärinen & Nieminen, 1990; Nyman et al., 2014). However, the mtDNA phylogeny presented here (Figure 6) suggests that the modern Baltic ringed seal population is a poor proxy for the ancestral Saimaa population. All the Saimaa haplotypes form a cluster that is not nested within the Baltic haplotypes. In fact, none of the geographically adjacent populations in the Baltic Sea, Lake Ladoga, and the White Sea are close to the Saimaa ringed seal clus- ter. Rather, the Saimaa ringed seal clade is nested within the Arctic, especially North American, ringed seal diversity, which suggests that the history of the Saimaa ringed seal population is more com- plex than previously thought. The fact that some Baltic ringed seal haplotypes are nested inside the Saimaa ringed seal mtDNA clade (Figures 2 and 6) might conceivably suggest returning gene flow from Lake Saimaa into the Baltic population. However, the current outlet of Lake Saimaa is River Vuoksi, which runs into Lake Ladoga, from which the route to the Baltic Sea continues via River Neva. The absence of Saimaa- clade haplotypes from Lake Ladoga therefore suggests that their presence in the Baltic Sea reflects incomplete lineage sorting. As already noted by Palo et al. (2003) and Valtonen et al. (2012), particularly problematic for the postglacial “Out of Baltic” scenario for the Saimaa ringed seal is the presence of multiple closely related, but still diverse, haplotypes in Lake Saimaa. It is unlikely that a ran- dom sampling (through lineage sorting either during the initial colo- nization phase or later) from the Baltic haplotype pool would retain such a closely related set of haplotypes. The presence of a single mtDNA clade within Lake Saimaa could theoretically be explained by an extended colonization bottleneck (Nyman et al., 2014) and differentiation of a few founding lineages during the isolation. The current diversity within Lake Saimaa is, however, too high to have emerged postglacially. As pointed out by Valtonen et al. (2012), the accumulation of the diversity observed in the Saimaa clade would re- quire at least 95,000 years (c. 8600 generations), or Saimaa- specific mutation rates that are roughly 10x faster than normally estimated for mammalian mtDNA. We note that most of the posterior sup- port values linking the Saimaa ringed seal clade together with North American ringed seals are not particularly strong in our phylogenetic tree and could possibly lead to spurious population connections. We however also note that similar connection has been previously re- ported by Palo et al. (2003). The observation of one Saimaa ringed seal individual with a characteristic 64- bp deletion is interesting. This individual not only shared the same deletion with two North American ringed seals, but also SNPs. This suggests that the deletion F I G U R E 6 Bayesian phylogenetic tree of 211 unique mitochondrial control- region haplotypes found in freshwater, Baltic, and Arctic ringed seal populations. Symbols and colors at the tips of the tree indicate the areas where each haplotype is present (see inset map and legend). Numbers above branches are posterior probabilities; only values over 0.50 are shown. The location of the Saimaa ringed seal individual with a 64- bp deletion (haplotype H14) is indicated by an arrow. 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 12 of 16  |     HEINO et al. more likely stems from shared ancestry rather than from convergent deletion events of the same sequence region. As this group of three seals is part of the Saimaa ringed seal clade, it further suggests that the diversity observed within Saimaa does not solely stem from postglacial diversification of a single lineage after the colonization of Lake Saimaa but is at least partly of ancient origin. The observed pattern could have emerged if the Baltic Sea area has experienced multiple colonization waves of ringed seals after the last glacial period (Schmölcke, 2008), so that the Saimaa pop- ulation represents survivors of an early wave, and the extant Baltic and Ladoga populations derive from a later wave that replaced the earlier Baltic population. Based on ancient specimens, Bro- Jørgensen (2021) inferred that such a replacement has happened in the Baltic Sea gray seal (Halichoerus grypus). However, according to Ukkonen et al. (2014), there is no paleontological evidence suggest- ing that the Baltic ringed seal population would have experienced extinctions and recolonizations after the initial immigration into the basin after the last glacial period, though additional immigration from the Arctic may have occurred. On the other hand, the pat- tern could be explained if Saimaa ringed seals originated from an ancestral population that inhabited proglacial lake systems situated on the edge of the Fennoscandian Ice Sheet. Current models sup- port the existence of such refugia for other aquatic species, such as the Atlantic salmon (Salmo salar L.) (Asplund et al., 2004; Nilsson et al., 2001; Tonteri et al., 2005, 2007), grayling (Thymallus thymallus) (Koskinen et al., 2000), and European perch (Perca fluviatilis) (Nesbø et al., 1999). 5  |  CONCLUSIONS The use of museum specimens allowed us to directly investigate genetic patterns in the endangered Saimaa ringed seal popula- tion through more than a 100 years. Although erosion of the ini- tial genetic diversity has continued throughout the isolation, we observed 20th- century loss of haplotypes and relatively drastic fluctuations in haplotype frequencies, demonstrating a genetic ef- fect of the human- induced population collapse. Combining newly generated and already published data from multiple ringed seal populations additionally allowed us to investigate broad phylo- geographic patterns in ringed seals. In apparent contrast with the traditionally held view of the Baltic origin of the population, the Saimaa ringed seal mtDNA variation shows enigmatic affinities to North American ringed seals. These results add to the grow- ing body of evidence which calls for a re- evaluation of the deep history of the Saimaa ringed seal population. Future data on still- unsampled populations, for example, in the Arctic Ocean, as well as genomic data and ancient DNA could provide keys to under- standing the origin and demographic history of the Saimaa ringed seal. Lake Saimaa may harbor a ringed seal population that is even more unique than previously thought, which calls to strengthen the conservation efforts of this population further. AUTHOR CONTRIBUTIONS Matti T. Heino: Conceptualization (equal); data curation (equal); formal analysis (equal); funding acquisition (equal); investiga- tion (equal); methodology (equal); project administration (equal); validation (equal); visualization (supporting); writing – original draft (equal); writing – review and editing (equal). Tommi Nyman: Conceptualization (equal); data curation (supporting); formal analy- sis (equal); funding acquisition (equal); investigation (equal); meth- odology (equal); project administration (equal); validation (equal); visualization (lead); writing – original draft (equal); writing – re- view and editing (equal). Jukka U. Palo: Conceptualization (equal); formal analysis (equal); funding acquisition (equal); investigation (equal); methodology (equal); resources (equal); validation (equal); writing – original draft (equal); writing – review and editing (equal). Jenni Harmoinen: Data curation (equal); formal analysis (equal); investigation (equal); methodology (equal); validation (equal); writ- ing – review and editing (equal). Mia Valtonen: Conceptualization (equal); data curation (equal); formal analysis (equal); funding acquisi- tion (equal); investigation (equal); validation (equal); writing – origi- nal draft (equal); writing – review and editing (equal). Małgorzata Pilot: Conceptualization (equal); data curation (equal); formal analy- sis (equal); funding acquisition (equal); investigation (equal); meth- odology (equal); project administration (equal); resources (equal); validation (equal); writing – original draft (equal); writing – review and editing (equal). Sanni Översti: Data curation (equal); formal anal- ysis (equal); investigation (equal); methodology (equal); validation (equal); visualization (supporting); writing – original draft (equal); writing – review and editing (equal). Elina Salmela: Formal analysis (equal); investigation (equal); methodology (equal); software (lead); writing – original draft (equal); writing – review and editing (equal). Mervi Kunnasranta: Conceptualization (equal); data curation (equal); funding acquisition (equal); investigation (equal); resources (equal); writing – review and editing (equal). Risto Väinölä: Data curation (equal); investigation (equal); resources (equal); writing – review and editing (equal). A. Rus Hoelzel: Conceptualization (equal); formal analysis (equal); investigation (equal); supervision (equal); writing – review and editing (equal). Jouni Aspi: Conceptualization (equal); funding acquisition (equal); investigation (equal); project administra- tion (equal); resources (equal); supervision (equal); writing – original draft (equal); writing – review and editing (equal). ACKNOWLEDG MENTS This work was carried out with the support of the Centre for Material Analysis, University of Oulu, Finland. The authors wish to acknowledge CSC— IT Center for Science, Finland, for compu- tational resources and the Finnish Museum of Natural History for access to their collections. MTH acknowledges funding from the Emil Aaltonen Foundation and the University of Oulu Scholarship Foundation. MV was funded by the Maj and Tor Nessling Foundation and Jane and Aatos Erkko Foundation. MP was sup- ported by Marie Skłodowska Curie Intra- European Fellowship from the European Commission (PIEF- GA- 2009- 235978) and the 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  13 of 16HEINO et al. Polish National Agency for Academic Exchange (NAWA; Polish Returns Fellowship PPN/PPO/2018/1/00037). ES acknowledges funding from the Jenny and Antti Wihuri Foundation and the Ella and Georg Ehrnrooth Foundation. The Raija and Ossi Tuuliainen Foundation, Kuopio Naturalists' Society, and Nestori Foundation are also acknowledged for funding the research. We like to thank Jaakko Lumme and Laura Kvist from the University Oulu and Love Dalén from the Centre for Palaeogenetics for discussion. We further thank Pirkko Ukkonen for her help with the material collection from museum specimens. The work is dedicated to our deceased colleague Minna Ruokonen, who had an instrumental role in conceiving, designing, and supervising the research and in acquiring funding. FUNDING INFORMATION Emil Aaltonen Foundation; University of Oulu Scholarship Foundation; Maj and Tor Nessling Foundation; Jane and Aatos Erkko Foundation; Marie Skłodowska Curie Intra- European Fellowship from the European Commission (PIEF- GA- 2009- 235978); Polish National Agency for Academic Exchange (NAWA; Polish Returns Fellowship PPN/PPO/2018/1/00037); Raija and Ossi Tuuliainen Foundation; Kuopio Naturalists' Society; Nestori Foundation; Jenny and Antti Wihuri Foundation; Ella and Georg Ehrnrooth Foundation. CONFLIC T OF INTERE S T The authors declare no conflicts of interest. OPEN RE SE ARCH BADG E S This article has earned an Open Data badge for making publicly available the digitally- shareable data necessary to reproduce the reported results. The data is available at https://github.com/esalm ela/Conti Unity. DATA AVAIL ABILIT Y S TATEMENT The sequences reported in this study have been deposited in GenBank under accession numbers ON841430- ON841496. The script that was used for haplotype resampling is available in https:// github.com/esalm ela/Conti Unity. ORCID Matti T. Heino https://orcid.org/0000-0003-3812-3543 Tommi Nyman https://orcid.org/0000-0003-2061-0570 Jukka U. Palo https://orcid.org/0000-0002-1955-7963 Jenni Harmoinen https://orcid.org/0000-0001-5402-7127 Mia Valtonen https://orcid.org/0000-0003-2034-2019 Małgorzata Pilot https://orcid.org/0000-0002-6057-5768 Sanni Översti https://orcid.org/0000-0001-7118-1281 Elina Salmela https://orcid.org/0000-0003-1326-4462 Mervi Kunnasranta https://orcid.org/0000-0003-3612-8842 Risto Väinölä https://orcid.org/0000-0003-0056-5741 A. Rus Hoelzel https://orcid.org/0000-0002-7265-4180 Jouni Aspi https://orcid.org/0000-0002-2451-3201 R E FE R E N C E S Asplund, T., Veselov, A., Primmer, C. R., Bakhmet, I., Potutkin, A., Titov, S., & Lumme, J. (2004). Geographical structure and postglacial history of mtDNA haplotype variation in Atlantic salmon (Salmo salar L.) among rivers of the white and Barents Sea basins. Annales Zoologici Fennici, 41(3), 465– 475. Baele, G., Lemey, P., Bedford, T., Rambaut, A., Suchard, M. A., & Alekseyenko, A. V. (2012). Improving the accuracy of demographic and molecular clock model comparison while accommodating phy- logenetic uncertainty. Molecular Biology and Evolution, 29(9), 2157– 2167. https://doi.org/10.1093/molbe v/mss084 Berta, A., & Churchill, M. (2012). Pinniped taxonomy: Review of cur- rently recognized species and subspecies, and evidence used for their description. Mammal Review, 42(3), 207– 234. https://doi. org/10.1111/j.1365- 2907.2011.00193.x Bi, K., Linderoth, T., Singhal, S., Vanderpool, D., Patton, J. L., Nielsen, R., Moritz, C., & Good, J. M. (2019). Temporal genomic contrasts reveal rapid evolutionary responses in an alpine mammal during recent climate change. PLoS Genetics, 15(5), e1008119. https://doi. org/10.1371/journ al.pgen.1008119 Borchsenius, F. (2009). FastGap 1.2. http://www.aubot.dk/FastG ap_ home.htm Bouckaert, R., Vaughan, T. G., Barido- Sottani, J., Duchêne, S., Fourment, M., Gavryushkina, A., Heled, J., Jones, G., Kühnert, D., de Maio, N., Matschiner, M., Mendes, F. K., Müller, N. F., Ogilvie, H. A., du Plessis, L., Popinga, A., Rambaut, A., Rasmussen, D., Siveroni, I., … Drummond, A. J. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Computational Biology, 15(4), e1006650. https://doi.org/10.1371/ journ al.pcbi.1006650 Bouckaert, R. R., & Drummond, A. J. (2017). bModelTest: Bayesian phylogenetic site model averaging and model comparison. BMC Evolutionary Biology, 17(1), 1– 11. https://doi.org/10.1186/s1286 2- 017- 0890- 6 Bro- Jørgensen, M. H. (2021). Ancient genomics of Baltic seals: Insights on the past Baltic grey seal and harp seal populations (PhD dissertation, Stockholm University, Department of Archaeology and Classical Studies). http://urn.kb.se/resol ve?urn=urn:nbn:se:su:diva- 189293 Burrell, A. S., Disotell, T. R., & Bergey, C. M. (2015). The use of mu- seum specimens with high- throughput DNA sequencers. Journal of Human Evolution, 79, 35– 44. https://doi.org/10.1016/j. jhevol.2014.10.015 Chao, A., Gotelli, N. J., Hsieh, T. C., Sander, E. L., Ma, K. H., Colwell, R. K., & Ellison, A. M. (2014). Rarefaction and extrapolation with hill numbers: A framework for sampling and estimation in species di- versity studies. Ecological Monographs, 84(1), 45– 67. https://doi. org/10.1890/13- 0133.1 Chao, A., Ma, K. H., & Hsieh, T. C. (2016). iNEXT (iNterpolation and EXTrapolation) online: Software for interpolation and extrapola- tion of species diversity. Program and User's Guide published at http://chao.stat.nthu.edu.tw/wordp ress/softw are_downl oad/ inext - onlin e/ Clement, M., Posada, D., & Crandall, K. A. (2000). TCS: a computer pro- gram to estimate gene genealogies. Molecular Ecology, 9(10), 1657– 1659. https://doi.org/10.1046/j.1365- 294x.2000.01020.x de Bruyn, M., Hall, B. L., Chauke, L. F., Baroni, C., Koch, P. L., & Hoelzel, A. R. (2009). Rapid response of a marine mammal species to Holocene climate and habitat change. PLoS Genetics, 5(7), e1000554. https:// doi.org/10.1371/journ al.pgen.1000554 Drummond, A. J., Rambaut, A., Shapiro, B., & Pybus, O. G. (2005). Bayesian coalescent inference of past population dynamics from 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 14 of 16  |     HEINO et al. molecular sequences. Molecular Biology and Evolution, 22(5), 1185– 1192. https://doi.org/10.1093/molbe v/msi103 Dufresnes, C., Miquel, C., Remollino, N., Biollaz, F., Salamin, N., Taberlet, P., & Fumagalli, L. (2018). Howling from the past: Historical phylogeog- raphy and diversity losses in European grey wolves. Proceedings of the Royal Society B, 285(1884). https://doi.org/10.1098/RSPB.2018.1148 Dussex, N., van der Valk, T., Morales, H. E., Wheat, C. W., Díez- del- Molino, D., von Seth, J., & Dalén, L. (2021). Population genomics of the critically endangered kākāpō. Cell Genomics, 1(1), 100002. https://doi.org/10.1016/j.xgen.2021.100002 Foote, A. D., Hooper, R., Alexander, A., Baird, R. W., Baker, C. S., Ballance, L., Barlow, J., Brownlow, A., Collins, T., Constantine, R., Dalla Rosa, L., Davison, N. J., Durban, J. W., Esteban, R., Excoffier, L., Martin, S. L. F., Forney, K. A., Gerrodette, T., Gilbert, M. T. P., … Morin, P. A. (2021). Runs of homozygosity in killer whale genomes provide a global record of demographic histories. Molecular Ecology, 30(23), 6162– 6177. https://doi.org/10.1111/mec.16137 Hansen, A., Willerslev, E., Wiuf, C., Mourier, T., & Arctander, P. (2001). Statistical evidence for miscoding lesions in ancient DNA tem- plates. Molecular Biology and Evolution, 18(2), 262– 265. https://doi. org/10.1093/oxfor djour nals.molbev.a003800 Heller, R., Chikhi, L., & Siegismund, H. R. (2013). The confounding ef- fect of population structure on Bayesian skyline plot inferences of demographic history. PLoS One, 8(5), e62992. https://doi. org/10.1371/journ al.pone.0062992 Hoban, S., Bertorelle, G., & Gaggiotti, O. E. (2012). Computer simulations: Tools for population and evolutionary genetics. Nature Reviews. Genetics, 13(2), 110– 122. https://doi.org/10.1038/nrg3130 Hofreiter, M., Jaenicke, V., Serre, D., Haeseler, A., & Pääbo, S. (2001). DNA sequences from multiple amplifications reveal artifacts in- duced by cytosine deamination in ancient DNA. Nucleic Acids Research, 29(23), 4793– 4799. Hofreiter, M., Paijmans, J. L. A., Goodchild, H., Speller, C. F., Barlow, A., Fortes, G. G., Thomas, J. A., Ludwig, A., & Collins, M. J. (2015). The future of ancient DNA: Technical advances and conceptual shifts. BioEssays: News and Reviews in Molecular, Cellular and Developmental Biology, 37(3), 284– 293. https://doi.org/10.1002/bies.20140 0160 Hung, C. M., Shaner, P. J. L., Zink, R. M., Liu, W. C., Chu, T. C., Huang, W. S., & Li, S. H. (2014). Drastic population fluctuations explain the rapid extinction of the passenger pigeon. Proceedings of the National Academy of Sciences of the United States of America, 111(29), 10636– 10641. https://doi.org/10.1073/pnas.14015 26111 Huynen, L., Millar, C. D., & Lambert, D. M. (2012). Resurrecting ancient animal genomes: The extinct moa and more. BioEssays: News and Reviews in Molecular, Cellular and Developmental Biology, 34(8), 661– 669. https://doi.org/10.1002/bies.20120 0040 Hyvärinen, H., Kunnasranta, M., Nieminen, P., & Taskinen, J. (2004). Hyle: Saimaan oma norppa. Tammi. Hyvärinen, H., & Nieminen, M. (1990). Differentiation of the ringed seal in the Baltic Sea, Lake Ladoga and Lake Saimaa. Finnish Game Research, 47, 21– 27. Hyvärinen, H., Sipilä, T., Koskela, J., & Kunnasranta, M. (1999). The Saimaa ringed seal. In E. Kuusisto (Ed.), Saimaa, a Living Lake (pp. 126– 136). Tammi. Hyvärinen, H., Sipila, T., Kunnasranta, M., & Koskela, J. T. (1998). Mercury pollution and the Saimaa ringed seal (Phoca hispida saimensis). Marine Pollution Bulletin, 36(1), 76– 81. https://doi.org/10.1016/ s0025 - 326x(98)90037 - 6 James, J., & Eyre- Walker, A. (2020). Mitochondrial DNA sequence diver- sity in mammals: A correlation between the effective and census population sizes. Genome Biology and Evolution, 12(12), 2441– 2449. doi:10.1093/gbe/evaa222 Katoh, K., Misawa, K., Kuma, K. I., & Miyata, T. (2002). MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research, 39(14), 3059– 3066. https://doi. org/10.1093/nar/gkf436 Katoh, K., Rozewicki, J., & Yamada, K. D. (2019). MAFFT online service: Multiple sequence alignment, interactive sequence choice and vi- sualization. Briefings in Bioinformatics, 20(4), 1160– 1166. https:// doi.org/10.1093/bib/bbx108 Katoh, K., & Standley, D. M. (2013). MAFFT multiple sequence align- ment software version 7: Improvements in performance and us- ability. Molecular Biology and Evolution, 30(4), 772– 780. https://doi. org/10.1093/molbe v/mst010 Kokko, H., Helle, E., Lindström, J., Ranta, E., Sipilä, T., & Courchamp, F. (1999). Backcasting population sizes of ringed and grey seals in the Baltic and Lake Saimaa during the 20th century. Annales Zoologici Fennici, 36(2), 65– 73. Kokko, H., Lindström, J., Ranta, E., Sipilä, T., & Koskela, J. (1998). Estimating the demographic effective population size of the Saimaa ringed seal (Phoca hispida saimensis Nordq.). Animal Conservation Forum, 1(1), 47– 54. https://doi.org/10.1111/j.1469- 1795.1998. tb002 25.x Koskinen, M. T., Ranta, E., Piironen, J., Veselov, A., Titov, S., Haugen, T. O., Nilsson, J., Carlstein, M., & Primmer, C. R. (2000). Genetic lineages and postglacial colonization of grayling (Thymallus thy- mallus, Salmonidae) in Europe, as revealed by mitochondrial DNA analyses. Molecular Ecology, 9(10), 1609– 1624. https://doi. org/10.1046/j.1365- 294x.2000.01065.x Kunnasranta, M. (2001). Behavioural biology of two ringed seal (Phoca hispida) subspecies in the large European lakes Saimaa and Ladoga (PhD dissertation, University of Eastern Finland, Department of Biology). Kunnasranta, M., Niemi, M., Auttila, M., Valtonen, M., Kammonen, J., & Nyman, T. (2021). Sealed in a lake— Biology and conservation of the endangered Saimaa ringed seal: A review. Biological Conservation, 253, 108908. https://doi.org/10.1016/j.biocon.2020.108908 Liukko, U.- M., Henttonen, H., Hanski, I. K., Kauhala, K., Kojola, I., Kyheröinen, E.- M., & Pitkänen, J. (2016). The 2015 red list of Finnish mammal species. Ympäristöministeriö and Suomen Ympäristökeskus. Martinez- Bakker, M. E., Sell, S. K., Swanson, B. J., Kelly, B. P., & Tallmon, D. A. (2013). Combined genetic and telemetry data reveal high rates of gene flow, migration, and long- distance dispersal poten- tial in Arctic ringed seals (Pusa hispida). PLoS One, 8(10), e77125. https://doi.org/10.1371/journ al.pone.0077125 Metsähallitus. (2020). Saimaannorppa. https://www.metsa.fi/luont o- ja- kultt uurip erint o/lajie n- suoje lu/saima annor ppa/ Miller, M. A., Pfeiffer, W., & Schwartz, T. (2010). Creating the CIPRES Science Gateway for inference of large phylogenetic trees. 2010 Gateway computing environments workshop, GCE 2010. https:// doi.org/10.1109/gce.2010.5676129 Múrias Dos Santos, A., Cabezas, M. P., Tavares, A. I., Xavier, R., & Branco, M. (2016). tcsBU: A tool to extend TCS network layout and visu- alization. Bioinformatics, 32(4), 627– 628. https://doi.org/10.1093/ bioin forma tics/btv636 Nakahama, N. (2021). Museum specimens: An overlooked and valuable material for conservation genetics. Ecological Research, 36(1), 13– 23. https://doi.org/10.1111/1440- 1703.12181 Nesbø, C. L., Fossheim, T., Vøllestad, L. A., & Jakobsen, K. S. (1999). Genetic divergence and phylogeographic relationships among European perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization. Molecular Ecology, 8(9), 1387– 1404. https://doi.org/10.1046/j.1365- 294x.1999.00699.x Nilsson, J., Gross, R., Asplund, T., Dove, O., Jansson, H., Kelloniemi, J., Kohlmann, K., Löytynoja, A., Nielsen, E. E., Paaver, T., Primmer, C. R., Titov, S., Vasemägi, A., Veselov, A., Ost, T., & Lumme, J. (2001). Matrilinear phylogeography of Atlantic salmon (Salmo salar L.) in Europe and postglacial colonization of the Baltic Sea area. Molecular Ecology, 10(1), 89– 102. https://doi. org/10.1046/j.1365- 294x.2001.01168.x Nyman, T., Valtonen, M., Aspi, J., Ruokonen, M., Kunnasranta, M., & Palo, J. U. (2014). Demographic histories and genetic diversities 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License     |  15 of 16HEINO et al. of Fennoscandian marine and landlocked ringed seal subspecies. Ecology and Evolution, 4(17), 3420– 3434. https://doi.org/10.1002/ ece3.1193 O'Corry- Crowe, G. (2008). Climate change and the molecular ecology of Arctic marine mammals. Ecological Applications, 18, S56– S76. https://doi.org/10.1890/06- 0795.1 Olsen, M. T., Volny, V. H., Bérubé, M., Dietz, R., Lydersen, C., Kovacs, K. M., Dodd, R. S., & Palsbøll, P. J. (2011). A simple route to single- nucleotide polymorphisms in a nonmodel species: Identification and characterization of SNPs in the artic ringed seal (Pusa his- pida hispida). Molecular Ecology Resources, 11, 9– 19. https://doi. org/10.1111/j.1755- 0998.2010.02941.x Pääbo, S., Poinar, H., Serre, D., Jaenicke- Despres, V., Hebler, J., Rohland, N., Kuch, M., Krause, J., Vigilant, L., & Hofreiter, M. (2004). Genetic analyses from ancient DNA. Annual Review of Genetics, 38, 645– 679. https://doi.org/10.1146/annur ev.genet.37.110801.143214 Palo, J. (2003). Genetic diversity and phylogeography of landlocked seals (PhD dissertation, University of Helsinki, Department of Ecology and Systematics). Palo, J. U., Hyvärinen, H., Helle, E., Mäkinen, H. S., & Väinölä, R. (2003). Postglacial loss of microsatellite variation in the landlocked Lake Saimaa ringed seal. Conservation Genetics, 4(2), 117– 128. https:// doi.org/10.1023/a:10233 03109701 Palo, J. U., Mäkinen, H. S., Helle, E., Stenman, O., & Väinölä, R. (2001). Microsatellite variation in ringed seals (Phoca hispida): Genetic structure and history of the Baltic Sea population. Heredity, 86(5), 609– 617. https://doi. org/10.1046/j.1365- 2540.2001.00859.x Peart, C. R., Tusso, S., Pophaly, S. D., Botero- Castro, F., Wu, C. C., Aurioles- Gamboa, D., Baird, A. B., Bickham, J. W., Forcada, J., Galimberti, F., Gemmell, N. J., Hoffman, J. I., Kovacs, K. M., Kunnasranta, M., Lydersen, C., Nyman, T., de Oliveira, L. R., Orr, A. J., Sanvito, S., … Wolf, J. B. W. (2020). Determinants of genetic variation across eco- evolutionary scales in pinnipeds. Nature Ecology & Evolution, 4(8), 1095– 1104. https://doi.org/10.1038/s4155 9- 020- 1215- 5 Prost, S., & Anderson, C. N. K. (2011). TempNet: A method to display statistical parsimony networks for heterochronous DNA sequence data. Methods in Ecology and Evolution, 2(6), 663– 667. https://doi. org/10.1111/j.2041- 210x.2011.00129.x Rambaut, A., Drummond, A. J., Xie, D., Baele, G., & Suchard, M. A. (2018). Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Systematic Biology, 67(5), 901– 904. https://doi.org/10.1093/ sysbi o/syy032 Ranta, E., Lindström, J., & Kokko, H. (1996). Ecological risk analysis: The case of the Saimaa ringed seal. Ambio, 25(5), 363– 365. Raxworthy, C. J., & Smith, B. T. (2021). Mining museums for historical DNA: Advances and challenges in museomics. Trends in Ecology & Evolution, 36(11), 1049– 1060. https://doi.org/10.1016/j. tree.2021.07.009 Reeves, R. R. (1998). Distribution, abundance and biology of ringed seals (Phoca hispida): An overview. NAMMCO Scientific Publications, 1, 9– 45. https://doi.org/10.7557/3.2979 Rieux, A., Eriksson, A., Li, M., Sobkowiak, B., Weinert, L. A., Warmuth, V., Ruiz- Linares, A., Manica, A., & Balloux, F. (2014). Improved cal- ibration of the human mitochondrial clock using ancient genomes. Molecular Biology and Evolution, 31(10), 2780– 2792. https://doi. org/10.1093/molbe v/msu222 Rieux, A., & Khatchikian, C. E. (2017). Tipdatingbeast: An r package to assist the implementation of phylogenetic tip- dating tests using Beast. Molecular Ecology Resources, 17(4), 608– 613. https://doi. org/10.1111/1755- 0998.12603 Rizzi, E., Lari, M., Gigli, E., de Bellis, G., & Caramelli, D. (2012). Ancient DNA studies: New perspectives on old samples. Genetics Selection Evolution, 44(1), 1– 19. https://doi.org/10.1186/1297- 9686- 44- 21 Rohland, N., Siedel, H., & Hofreiter, M. (2010). A rapid column- based ancient DNA extraction method for increased sample through- put. Molecular Ecology Resources, 10(4), 677– 683. https://doi. org/10.1111/j.1755- 0998.2009.02824.x Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D. L., Darling, A., Höhna, S., Larget, B., Liu, L., Suchard, M. A., & Huelsenbeck, J. P. (2012). MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic Biology, 61(3), 539– 542. https://doi.org/10.1093/sysbi o/sys029 Schmölcke, U. (2008). Holocene environmental changes and the seal (Phocidae) fauna of the Baltic Sea: Coming, going and staying. Mammal Review, 38(4), 231– 246. https://doi. org/10.1111/j.1365- 2907.2008.00131.x Simmons, M. P., & Ochoterena, H. (2000). Gaps as characters in sequence- based phylogenetic analyses. Systematic Biology, 49(2), 369– 381. https://doi.org/10.1093/sysbi o/49.2.369 Sipilä, T., Helle, E., & Hyvärinen, H. (1990). Distribution, population size and reproductivity of the Saimaa ringed seal (Phoca hispida saimen- sis Nordq.) in Finland, 1980– 1984. Finnish Game Research, 47, 3– 10. Stoffel, M. A., Humble, E., Paijmans, A. J., Acevedo- Whitehouse, K., Chilvers, B. L., Dickerson, B., Galimberti, F., Gemmell, N. J., Goldsworthy, S. D., Nichols, H. J., Krüger, O., Negro, S., Osborne, A., Pastor, T., Robertson, B. C., Sanvito, S., Schultz, J. K., Shafer, A. B. A., Wolf, J. B. W., & Hoffman, J. I. (2018). Demographic histo- ries and genetic diversity across pinnipeds are shaped by human exploitation, ecology and life- history. Nature Communications, 9(1), 1– 12. https://doi.org/10.1038/s4146 7- 018- 06695 - z Suchard, M. A., Lemey, P., Baele, G., Ayres, D. L., Drummond, A. J., & Rambaut, A. (2018). Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evolution, 4(1). https://doi. org/10.1093/ve/vey016 Tamura, K., & Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution, 10(3), 512– 526. https://doi.org/10.1093/oxfor djour nals.molbev.a040023 Taylor, R. S., Jensen, E. L., Coltman, D. W., Foote, A. D., & Lamichhaney, S. (2021). Seeing the whole picture: What molecular ecology is gain- ing from whole genomes. Molecular Ecology, 30(23), 5917– 5922. https://doi.org/10.1111/mec.16282 Tonteri, A., Veselov, A. J., Titov, S., Lumme, J., & Primmer, C. R. (2007). The effect of migratory behaviour on genetic diversity and pop- ulation divergence: A comparison of anadromous and freshwater Atlantic salmon Salmo salar. Journal of Fish Biology, 70, 381– 398. https://doi.org/10.1111/j.1095- 8649.2007.01519.x Tonteri, A., Titov, S., Veselov, A., Zubchenko, A., Koskinen, M. T., Lesbarrères, D., & Primmer, C. R. (2005). Phylogeography of anad- romous and non- anadromous Atlantic salmon (Salmo salar) from northern Europe. Annales Zoologici Fennici, 42(1), 1– 22. Ukkonen, P. (2002). The early history of seals in the northern Baltic. Annales Zoologici Fennici, 39, 187– 207. Ukkonen, P., Aaris- Sørensen, K., Arppe, L., Daugnora, L., Halkka, A., Lõugas, L., & Storå, J. (2014). An Arctic seal in temperate waters: History of the ringed seal (Pusa hispida) in the Baltic Sea and its ad- aptation to the changing environment. The Holocene, 24(12), 1694– 1706. https://doi.org/10.1177/09596 83614 551226 Valtonen, M., Heino, M., Aspi, J., Buuri, H., Kokkonen, T., Kunnasranta, M., & Nyman, T. (2015). Genetic monitoring of a critically- endangered seal population based on field- collected placentas. Annales Zoologici Fennici, 52(1– 2), 51– 65. https://doi.org/10.5735/086.052.0205 Valtonen, M., Palo, J. U., Aspi, J., Ruokonen, M., Kunnasranta, M., & Nyman, T. (2014). Causes and consequences of fine- scale pop- ulation structure in a critically endangered freshwater seal. BMC Ecology, 14(1), 1– 15. https://doi.org/10.1186/1472- 6785- 14- 22 Valtonen, M., Palo, J. U., Ruokonen, M., Kunnasranta, M., & Nyman, T. (2012). Spatial and temporal variation in genetic diversity of an 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License 16 of 16  |     HEINO et al. endangered freshwater seal. Conservation Genetics, 13(5), 1231– 1245. https://doi.org/10.1007/s1059 2- 012- 0367- 5 van der Valk, T., Díez- del- Molino, D., Marques- Bonet, T., Guschanski, K., & Dalén, L. (2019). Historical genomes reveal the genomic consequences of recent population decline in eastern gorillas. Current Biology, 29(1), 165– 170.e6. https://doi.org/10.1016/j.cub.2018.11.055 Vilaça, S. T., Piccinno, R., Rota- Stabelli, O., Gabrielli, M., Benazzo, A., Matschiner, M., Soares, L. S., Bolten, A. B., Bjorndal, K. A., & Bertorelle, G. (2021). Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species. Molecular Ecology, 30(23), 6178– 6192. https:// doi.org/10.1111/mec.16113 Villesen, P. (2007). FaBox: An online toolbox for FASTA se- quences. Molecular Ecology Notes, 7(6), 965– 968. https://doi. org/10.1111/j.1471- 8286.2007.01821.x von Seth, J., Dussex, N., Díez- del- Molino, D., van der Valk, T., Kutschera, V. E., Kierczak, M., Steiner, C. C., Liu, S., Gilbert, M. T. P., Sinding, M. S., Prost, S., Guschanski, K., Nathan, S. K. S. S., Brace, S., Chan, Y. L., Wheat, C. W., Skoglund, P., Ryder, O. A., Goossens, B., … Dalén, L. (2021). Genomic insights into the conservation status of the world's last remaining Sumatran rhinoceros populations. Nature Communications, 12(1), 1– 11. https://doi.org/10.1038/s4146 7- 021- 22386 - 8 Yang, Z. (1994). Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: Approximate meth- ods. Journal of Molecular Evolution, 39(3), 306– 314. https://doi. org/10.1007/bf001 60154 SUPPORTING INFORMATION Additional supporting information can be found online in the Supporting Information section at the end of this article. How to cite this article: Heino, M. T., Nyman, T., Palo, J. U., Harmoinen, J., Valtonen, M., Pilot, M., Översti, S., Salmela, E., Kunnasranta, M., Väinölä, R., Hoelzel, A. R., & Aspi, J. (2023). Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections. Ecology and Evolution, 13, e9720. https://doi. org/10.1002/ece3.9720 20457758, 2023, 1, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1002/ece3.9720 by Luonnonvarakeskus, W iley O nline Library on [28/02/2023]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on W iley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License