Genetic diversity of historical Atlantic walruses (Odobenus rosmarus rosmarus) from Bjørnøya and Håøya (Tusenøyane), Svalbard, Norway
BMC Research Notes volume 9, Article number: 112 (2016)
The population size of Atlantic walruses (Odobenus rosmarus rosmarus) is depleted relative to historical abundance levels. In Svalbard, centuries of over-exploitation brought the walrus herds to the verge of extinction, and such bottlenecks may have caused loss of genetic variation. To address this for Svalbard walruses, mitochondrial haplotypes of historical walruses from two major haul-out sites, Bjørnøya and Håøya, within the Archipelago were explored using bone samples from animals killed during the peak period of harvesting.
Using ancient DNA methodologies, the mitochondrial NADH dehydrogenase 1 (ND1) gene, the cytochrome c oxidase 1 (COI) gene, and the control region (CR) were targeted for 15 specimens from Bjørnøya (of which five were entirely negative) and 9 specimens from Håøya (of which one was entirely negative). While ND1 and COI sequences were obtained for only a few samples, the CR delivered the most comprehensive data set, and the average genetic distance among historic Svalbard samples was 0.0028 (SD = 0.0023).
The CR sequences from the historical samples appear to be nested among contemporary Atlantic walruses, and no distinct mitochondrial haplogroups were identified in the historical samples that may have been lost during the periods of extensive hunting. However, given the low sample size and poor phylogenetic resolution it cannot be excluded that such haplogroups existed.
The walrus (Odobenus rosmarus) has a circumpolar Arctic distribution. The Atlantic walrus (O. r. rosmarus) is distributed from the Canadian Arctic in the west to the Kara Sea in Russia in the east [1, 2]. The population size of Atlantic walruses is depleted relative to historical abundance levels due to over-harvesting. For example, the number of walruses in Central West Greenland (a population shared with Canada) was reduced by about 80 % between 1900 and 1960 [1, 3]. The exact number of total extant Atlantic walruses is not known but is estimated to be close to 35,000 . In Svalbard the first recorded walrus hunt took place in 1604 , but the extensive commercial hunting began around 1820 . By the middle of the 19th century the stock showed clear signs of declining, but harvesting continued. The centuries of walrus hunting in the Svalbard Archipelago brought the reportedly large herds to the verge of extinction . They were finally given total protection in 1952. It was estimated that around 1980 there were only approximately 100 individuals summering in Svalbard, most of them males . The population size of North Atlantic walruses in Svalbard has been increasing in recent decades. Lydersen et al.  estimated the number of walruses in Svalbard in 2006 to be about 2600. These authors assumed that the total population size including the individuals from Franz Josef Land to be approximately 5000. A more recent estimate from 2012 for Svalbard indicated a continued upward trend in the population size in recent years to approximately 3900 .
There is no data on the original walrus population size in Svalbard prior to hunting, but it is thought to have been very large [11, 12]. Substantial declines in population size (bottlenecks) can be accompanied by a loss of genetic variation, and this may have been the case for the walrus populations in the northern Barents Sea. Historically there were four major haul-out regions used by walruses in Svalbard (Fig. 1): (1) Bjørnøya; (2) Southeast Edgeøya including Tusenøyane; (3) Kvitøya and Nordaustland in the northeast; and (4) West- and Northwest Spitsbergen from Isfjorden to Wijdefjorden . In recent decades, walruses have started occupying former haul-out areas throughout the Archipelago (details of these observations can be found in a metadatabase at http://www.npolar.no/en/services/mms/), though sightings on Bjørnøya are still rare .
Various studies on the distribution [8, 13], and movement patterns [14–16], as well as population structure [17–19] have shown that the walruses in Svalbard and Franz Josef Land belong to the same population. In the summer months most males stay in Svalbard and most females and calves remain in the northeastern parts of Svalbard and in Franz Josef Land. Migration between the Svalbard–Franz Josef Land population and the neighboring population in East Greenland is thought to be rare. The relationship between the Svalbard–Franz Josef Land population and walruses in Novaya Zemlya and the Pechora and Kara Seas in the southeastern Barents Sea remains unknown [1, 16, 20, 21].
In this study, mitochondrial haplotype diversity of Atlantic walrus in the Svalbard area at the peak of the harvesting period was explored. Bone samples collected from historical haul-out sites on Bjørnøya (15 specimens) and Håøya (nine specimens) were included. Details regarding the samples are presented in Table 1. Unfortunately, radiocarbon dates are not available for these samples. However, the samples from Bjørnøya most likely originated from animals hunted before 1867 because hunting did not occur in this area subsequent to this date . The mandibles collected at Håøya in the Tusenøyane Archipelago are probably remains from one of the most famous slaughters of walruses in Svalbard , which occurred in 1852. The exact number of walruses present on the island at that time is not known, but it is believed that several thousands might have been hauled out, and many hundreds were killed in just 1 day . Later, Håøya was only rarely used as a haul-out site by walruses [23, 24], even into the first half of the 20th century , though haul-out at this site is more common today . It is therefore likely that the mandibles from Håøya are remains of the 1852 slaughter.
Ancient DNA methodologies used in this study followed the protocols described previously . In short, DNA was extracted from 0.1–0.5 g of finely ground bone powder as described in . Standard precautions and measures for ancient DNA analyses demonstrating authenticity of the obtained data were followed. The mitochondrial NADH dehydrogenase 1 (ND1) gene was targeted with amplification of 12, the cytochrome c oxidase 1 (COI) gene of four, and the control region (CR) of four overlapping fragments (Table 2). Herein, we report the nucleotide sequences obtained and their possible implications for the genetic structure of the Atlantic walruses in the Svalbard Archipelago.
Five samples from Bjørnøya (SBj05, SBj10, SBj12, SBj13, and SBj14) and one sample from Håøya (SH23) did not deliver any PCR products, and it was therefore concluded that no DNA survived in these bones. The SBj07 samples delivered only a short sequence for the mitochondrial control region. PCR success depended on the specific target region suggesting that some primer pairs were not optimally designed for the purpose of amplifying degraded walrus target DNA. This may be particularly the case for some ND1, but also to some extent for the targeted COI and CR fragments (Table 1), although primer pairs for the latter targets were successfully used in earlier studies .
Cytochrome c oxidase 1 (COI)
PCR amplification and sequencing of the COI was only successful for the 177 bp comprising region 3 for eight of the nine Håøya samples (Additional file 1). No sequence information was obtained for any of the 15 samples from Bjørnøya; all attempts for PCR amplification failed for DNA extracts from samples of this site. Two haplotypes were detected for the Håøya samples that differed by only one synonymous transition.
NADH dehydrogenase 1 (ND1)
Nucleotide sequence information was determined for samples from both Bjørnøya and Håøya (Additional file 2) for this gene. Again, a higher proportion of Håøya samples delivered nucleotide sequences than those from Bjørnøya. Nevertheless, for seven samples (two from Bjørnøya and five from Håøya) full sequence information was obtained. Nucleotide substitutions were detected at only six sites leading to seven distinct haplotypes for these samples.
Control region (CR)
This target region delivered sequence information for 18 samples (Additional file 3). Nucleotide substitutions were detected at eight positions. In the sequence obtained for SH20 one position remained ambiguous. A variable number of nucleotides were observed in a polyT/polyC region. This could reflect either natural variation or degradation of the old template DNA. However, similar variation was also observed in modern sequences (see respective sequences from available GenBank entries in Additional file 3), which may be taken as an indication that this reflects real natural variation. Nevertheless, if this variable polyT/polyC region is not considered there were seven distinct haplotypes in the old dataset.
The obtained CR sequences can be compared to two other recent studies that also reported mitchondrial CR sequence variation for walruses in the North Atlantic region. Lindqvist et al.  assessed the taxonomic status of the Laptev walrus (O. r. laptevi), including nine samples from the Tusenøyane (Svalbard), five samples from Franz Josef Land, and another five samples from Scoresby Sound (East Greenland). McLeod et al.  compared walrus specimens from the eastern Canadian Arctic, including the extirpated southeastern Canadian Maritimes walrus. Available CR sequences from these previous studies were downloaded from GenBank (Additional file 4: Table S1) and compared to the new data generated from the historical Håøya and Bjørnøya samples. The most comprehensive dataset was assembled for the CR, including information from 88 haplotypes (Additional file 3). Since only a few nucleotide sequences were obtained for COI and ND1, they will not be discussed further, and in the following, we focus on the discussion of the results obtained for the CR sequences. The CR dataset was subjected to maximum likelihood and Bayesian phylogenetic analyses in order to address whether the historical samples from Håøya or Bjørnøya could be assigned to specific mitochondrial lineages. The best-fit nucleotide substitution model was determined for the dataset using jModelTest v.1.1 . The HKY + I model was favored and implemented in the phylogenetic analyses. Bayesian analyses were run in BEAST v.1.8 package  with the Yule process tree prior, a lognormal relaxed clock model, for 50 million generations. The program Tracer v.1.5 was used to check for stationarity (the ESS was >200 for all the parameters sampled), and TreeAnnotator v.1.8 was used to summarize the set of post burn-in trees and their parameters (burn-in set to 500), and to produce a maximum clade credibility (MCC) phylogeny (shown in Fig. 2). A maximum likelihood tree was reconstructed using the RAxML Blackbox tool  available through the Cipres Web Portal (http://www.phylo.org). A phylogenetic network was created with Neighbor-Net  in SplitsTree v.4.13.1  using uncorrected p-distances. For genetic distance calculations the most parameterized model available in SplitsTree was chosen, with an HKY85 model, transitions: transversions weighted 2:1, gamma model of rate heterogeneity, with base frequencies estimated empirically. Nucleotide sequence distances for eastern Atlantic walruses (east Greenland, Svalbard, and Franz Josef Land) were calculated using the program DNADIST v.3.5 (copyright J. Felsenstein 1986–2008) as implemented in the program BioEdit v.7.2.5 .
The sequences obtained from the Håøya and Bjørnøya samples were very similar to sequences obtained from modern Atlantic walruses from the Svalbard Archipelago. In the Bayesian maximum clade credibility tree, the haplotypes of the Håøya and Bjørnøya samples group within a statistically well supported clade of Atlantic walruses. Within the Atlantic walrus clade there are essentially three subclades (A, B, and C), of which only clade C is statistically supported. The historical Håøya and Bjørnøya samples group within clade A, which also comprises walruses from East Greenland, Svalbard, and Franz Josef Land (Fig. 2). Clades B and C are comprised of modern walruses from West Greenland and Canada, as well as extirpated southeastern Canadian Maritimes walrus. A Maximum Likelihood (ML) analysis (Additional file 5: Figure S1) and a ML analysis on a trimmed dataset (Additional file 6: Figure S2), which was conducted in order to reduce the potential impact of missing information on the tree topology, supported the split between Atlantic and Pacific walruses. Within the Atlantic walrus only clade C received moderate bootstrap support in the analysis with the full dataset. Subclades A and B were rendered unresolved. The median-joining network also recognized the Pacific–Atlantic walrus split, and grouped the haplotypes from the historical samples from Bjørnøya and Håøya with modern walruses from East Greenland, Svalbard, and Franz Josef Land (Additional file 7: Figure S3). It is worth noting that for the most part the Håøya and Bjørnøya samples exhibited distinct haplotypes from modern walruses. Few shared haplotypes were observed between three Håøya individuals (SH17, SH19, SH24) and modern walruses from Tusenøyane in Svalbard (ATL01 and ATL03), and between one Bjørnøya sample and a modern individual from Lågøya in Svalbard (ATL04). Håøya is part of Tusenøyane, whereas Lågøya is situated in the northwest of the archipelago (in contrast to Bjørnøya which is the southernmost island in the archipelago). Based on nucleotide sequence distances between the historical walruses in this study (12 haplotypes), modern walruses from East Greenland, Svalbard, and Franz Josef Land (14 haplotypes), and modern Svalbard walruses (9 haplotypes), the average distances are slightly lower for the historical walruses (0.0028, SD = 0.0023) than for modern eastern Atlantic walruses (0.0050; SD = 0.0024) and Svalbard walruses (0.0046; SD = 0.0025).
Although there are distinct ancient haplotypes that are not recovered among the modern samples, the historical walrus haplotypes reported from Håøya and Bjørnøya appear to be nested among modern haplotypes from the Svalbard archipelago. Furthermore, higher nucleotide sequence variation was recovered from modern eastern Atlantic and Svalbard walruses as compared to the historical samples analyzed here. This might indicate that the intense hunting of the Atlantic walrus at haul-out sites on Bjørnøya and Håøya did not cause an extensive loss of mitochondrial sequence variation and possibly distinct mitochondrial haplogroups (evolutionary closely related haplotypes). However, the current study includes limited sampling and the mitochondrial sequence information provides poor phylogenetic resolution. Hence, it is possible that the relatively lower genetic variation recovered among our historical Atlantic walrus samples may be caused by limited sampling of ancient walrus genetic diversity. For example, the more remote eastern sectors of the Barents Sea may have served as a refugium for East Atlantic walruses during the height of walrus hunting, preserving most of the mitochondrial diversity found within the population today.
Bayesian information criterion
cytochrome c oxidase 1
effective sample size
maximum clade credibility
NADH dehydrogenase 1
polymerase chain reaction
Born EW, Gjertz I, Reeves RR. Population assessment of Atlantic walrus. Norsk Polarinst Meddelelser. 1995;138:1–100.
Wiig Ø, Born EW, Stewart RAE. Management of Atlantic walrus (Odobenus rosmarus rosmarus) in the arctic Atlantic. NAMMCO Scientific Publication. 2014. doi:10.7557/3.2855.
Witting L, Born EW. Population dynamics of walruses in Greenland. NAMMCO Scientific Publication. 2014. doi:10.7557/3.2612.
Laidre KL, Stern H, Kovacs KM, Lowry L, Moore SE, Regehr EV, Ferguson SH, Wiig Ø, Boveng P, Angliss AP, Born EW, Litovka D, Quakenbush L, Lydersen C, Vongraven D, Ugarte F. A circumpolar assessment of Arctic marine mammals and sea ice loss, with conservation recommendations for the 21st century. Cons Biol. 2015. doi:10.1111/cobi.12474.
Gjertz I, Wiig Ø. The number of walruses (Odobenus rosmarus) in Svalbard in summer. Polar Biol. 1995. doi:10.1007/BF00237468.
Lønø O. The catch of walrus (Odobenus rosmarus) in the areas of Svalbard, Novaja Zemlja. Frans Josef Land. Norsk Polarinst Årbok. 1970;1972:199–212.
Norderhaug M. Hvalrossens (Odobenus rosmarus) forekomst i Svalbardområdet 1960-1967. Norsk Polarinst Årbok. 1969;146-50.
Born EW. Status of the Atlantic walrus Odobenus rosmarus rosmarus in the Svalbard area. Polar Res. 1984. doi:10.1111/j.1751-8369.1984.tb00484.x.
Lydersen C, Aars J, Kovacs KM. Estimating the number of walruses in Svalbard from aerial surveys and behavioural data from satellite telemetry. Arctic. 2008;61:119–28.
Kovacs KM, Aars J, Lydersen C. Walruses recovering after 60+ years of protection in Svalbard. Norway. Polar Res. 2014. doi:10.3402/polar.v33.26034.
Reeves RR. Atlantic walrus (Odobenus rosmarus rosmarus): A literature survey and status report. Wildlife Research Report 10. Washington DC: US Fish and Wildlife Service, Department of the Interior; 1978.
Gjertz I, Wiig Ø, Øristland NA. Back calculation of original population size for walruses Odobenus rosmarus in Franz Josef Land. Wildlife Biology. 1998;4:223–30.
Gjertz I, Wiig Ø. Past and present distribution of walruses in Svalbard. Arctic. 1994;47:34–42.
Wiig Ø, Gjertz I, Griffiths D. Migration of walruses (Odobenus rosmarus) in the Svalbard and Franz Josef Land area. J Zool. 1996. doi:10.1111/j.1469-7998.1996.tb05429.x.
Freitas C, Kovacs KM, Ims RA, Fedak MA, Lydersen C. Deep into the ice: overwintering and habitat selection in Atlantic walruses. Mar Ecol Prog Ser. 2009. doi:10.3354/meps07725.
Lydersen C, Kovacs KM. Walrus Odobenus rosmarus research in Svalbard, Norway, 2000-2010. NAMMCO Scientific Publication. 2014. doi:10.7557/3.2613.
Andersen LW, Born EW, Gjertz I, Wiig Ø, Holm LE, Bendixen C. Population structure and gene flow of the Atlantic walrus (Odobenus rosmarus rosmarus) in the eastern Atlantic Arctic based on mitochondrial DNA and microsatellite variation. Mol Ecol. 1998. doi:10.1046/j.1365-294x.1998.00455.x.
Born EW, Andersen LW, Gjertz I, Wiig Ø. A review of the genetic relationships of Atlantic walrus (Odobenus rosmarus rosmarus) east and west of Greenland. Polar Biol. 2001. doi:10.1007/s003000100277.
Andersen LW, Born EW, Doidge DW, Gjertz I, Wiig Ø, Waples RS. Genetic signals of historic and recent migration between sub-populations of Atlantic walrus Odobenus rosamarus rosmarus west and east of Greenland. Endang Species Res. 2009. doi:10.3354/esr00242.
NAMMCO. Annual Report 2006–Volume I. North Atlantic Marine Mammal Commission. Tromsø, Norway: 2006; 277. http://www.nammco.no/webcronize/images/Nammco/907.pdf.
Boltunov AN, Belikov SE, Gorbunov YA, Menis DT, Semenova VS. The Atlantic walrus of the southeastern Barents Sea and adjacent regions: review of the present-day status. Marine Mammal Council: WWF-Russia; 2010.
Wiig Ø, Born EW, Gjertz I, Lydersen C, Stewart REA. Historical sex-specific distribution of Atlantic walrus (Odobenus rosmarus rosmarus) in Svalbard assessed by mandible measurements. Polar Biol. 2007. doi:10.1007/s00300-007-0334-7.
Guernsey AH. Sporting in Spitzbergen. Harper’s New Mon Mag. 1861;23:606–17.
Lamont J. Seasons with the sea-horses; or sporting adventures in the northern seas. London: Hurst and Blackett; 1861.
Lindqvist C, Bachmann L, Andersen LW, Born EW, Arnason U, Kovacs KM, Lydersen C, Abramov AV, Wiig Ø. The Laptev Sea walrus Odobenus rosmarus laptevi: an enigma revisited. Zool Scr. 2009. doi:10.1111/j.1463-6409.2008.00364.x.
Borge T, Bachmann L, Bjørnstad G, Wiig Ø. Genetic variation in Holocene bowhead whales from Svalbard. Mol Ecol. 2007. doi:10.1111/j.1365-294X.2007.03287.x.
McLeod BA, Frasier TR, Lucas Z. Assessment of the extirpated Maritimes walrus using morphological and ancient DNA analysis. PLoS One. 2014. doi:10.1371/journal.pone.0099569.
Posada D. JModelTest: phylogenetic model averaging. Mol Biol Evol. 2008. doi:10.1093/molbev/msn083.
Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007. doi:10.1186/1471-2148-7-21.
Stamatakis A, Hoover P, Rougemont J. A Rapid Bootstrap Algorithm for the RAxML Web-Servers. Syst Biol. 2008. doi:10.1080/10635150802429642.
Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006. doi:10.1093/molbev/msj030.
Huson DH. SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics. 1998. doi:10.1093/bioinformatics/14.1.68.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.
Arnason U, Adegoke JA, Bodin K, Born EW, Esa YB, Gullberg A, Nilsson M, Short RV, Xu X, Janke A. Mammalian mitogenomic relationships and the root of the eutherian tree. Proc Natl Acad Sci. 2002. doi:10.1073/pnas.102164299.
Cronin MA, Hills S, Born EW, Patton JC. Mitochondrial DNA variation in Atlantic and Pacific walruses. Can J Zool. 1994. doi:10.1139/z94-140.
All authors designed the study. CLi did the lab work, and CLi, TR, and LB conducted the analyses. ØW, KK, and CLy collected the samples. All authors contributed to the interpretation of the results and to writing of the manuscript. All authors read and approved the final version of the manuscript.
This work was financially supported by the Natural History Museum, University of Oslo and the Norwegian Polar Institute. We are grateful to Anders Skoglund, Norwegian Polar Institute, for preparing Fig. 1.
Availability of supporting data
The nucleotide sequences are available in GenBank; accession numbers KU710175–KU710220.
The authors declare that they have no competing interests.
The current study did not involve living animals. All bone samples were collected according to international guidelines.
Additional file 1. Aligned nucleotide sequences of the mitochondrial cytochrome c oxidase 1 (COI) obtained from historical walrus samples. Reference sequences from other Atlantic walrus samples obtained from GenBank are included.
Additional file 2. Aligned nucleotide sequences of the mitochondrial NADH dehydrogenase 1 (ND1) obtained from historical walrus samples. Reference sequences from other Atlantic walrus samples obtained from GenBank are included.
Additional file 3. Aligned nucleotide sequences of the mitochondrial control region (CR) obtained from historical walrus samples. Reference sequences from other Atlantic walrus samples obtained from GenBank are included.
Additional file 5: Figure S1. Maximum likelihood tree obtained from the analyses of the mitochondrial control region haplotypes from historical samples from Bjørnøya and Håøya and modern walrus samples (ATL + LAP + PAC) used in this study. The tree was reconstructed using the RAxML Blackbox tool  as implemented on the Cipres Web Portal (http://www.phylo.org). Bootstrap (BS) percentages of ≥50 are indicated.
Additional file 6: Figure S2. Maximum likelihood tree obtained from the analyses of the mitochondrial control region sequences from historical samples from Bjørnøya and Håøya and modern walrus samples (ATL + LAP + PAC) used in this study. The dataset was trimmed in order to minimize the impact of missing information on the analysis. The tree was reconstructed using the RAxML Blackbox tool  as implemented on the Cipres Web Portal (www.phylo.org). Bootstrap (BS) percentages of ≥50 are indicated.
Additional file 7: Figure S3. Neighbor network based on all mitochondrial control region sequences from historical samples from Bjørnøya and Håøya and modern walrus samples (ATL + LAP + PAC) used in this study. Nodes are represented by circles, and their sizes are proportional to the number of taxa sharing the haplotype.
About this article
Cite this article
Lindqvist, C., Roy, T., Lydersen, C. et al. Genetic diversity of historical Atlantic walruses (Odobenus rosmarus rosmarus) from Bjørnøya and Håøya (Tusenøyane), Svalbard, Norway. BMC Res Notes 9, 112 (2016). https://doi.org/10.1186/s13104-016-1907-8