Population genetic structure and evolutionary history of North Atlantic beluga whales (Delphinapterus leucas) from West Greenland, Svalbard and the White Sea

Population structure in many Arctic marine mammal species reflects a dynamic interplay between physical isolating mechanisms and the extent to which dispersal opportunities are met. We examined variation within mtDNA and eight microsatellite markers to investigate population structure and demographic history in beluga whales in the North Atlantic. Genetic heterogeneity was observed between Svalbard and West Greenland that reveals limited gene flow over ecological time scales. Differentiation was also recorded between Atlantic belugas and two previously studied populations in the North Pacific, the Beaufort Sea and Gulf of Alaska. However, Bayesian cluster analysis of the nDNA data identified two population clusters that did not correspond to the respective ocean basins, as predicted, but to: (1) Arctic (Svalbard-White Sea-Greenland-Beaufort Sea) and (2) Subarctic (Gulf of Alaska) regions. Similarly, the deepest phylogeographic signal was between the Arctic populations and the Gulf of Alaska. Fitting an isolation-with-migration model yielded genetic abundance estimates that match census estimates and revealed that Svalbard and the Beaufort Sea likely diverged 7,600-35,400 years ago but have experienced recurrent periods with gene flow since then, most likely via the Russian Arctic during subsequent warm periods. Considering current projections of continued sea ice losses in the Arctic, this study identified likely routes of future contact among extant beluga populations, and other mobile marine species, which have implications for genetic introgression, health, ecology and behavior.