Gene flow on ice: the role of sea ice and whaling in shaping Holarctic genetic diversity and population differentiation in bowhead whales (Balaena mysticetus)

Unraveling the Genetic Legacy of Commercial Whaling and Population Dynamics in Arctic Bowhead Whales and Narwhals

Assessing genetic structure and diversity in wildlife is particularly important in the context of climate change. The Arctic is rapidly warming, and endemic species must adapt quickly or face significant threats to persistence. Bowhead whales (Balaena mysticetus) and narwhals (Monodon monoceros) are two long-lived Arctic species with similar habitat requirements and are often seen together in the Canadian Arctic. Although their ranges overlap extensively, bowhead whales experienced significantly greater commercial whaling mortality than narwhals over several centuries. The similar habitat requirements but different harvest histories of these two species provide an opportunity to examine present-day genetic diversity and the demographic and genetic consequences of commercial whaling. We whole-genome resequenced contemporary Canadian Arctic bowhead whales and narwhals to delineate population structure and reconstruct demographic history. We found higher genetic diversity in bowhead whales compared to narwhals. However, bowhead whale effective population size sharply declined contemporaneously with the intense commercial whaling period. Narwhals, in contrast, exhibited recent growth in effective population size, likely reflecting exposure to limited opportunistic commercial harvest. Bowhead whales will likely continue to experience significant genetic drift in the future, leading to the erosion of genetic diversity. In contrast, narwhals do not seem to be at imminent risk of losing their current levels of genetic variation due to their long-term low effective population size and lack of evidence for a recent decline. This work highlights the importance of considering population trajectories in addition to genetic diversity when assessing the genetics of populations for conservation and management purposes.

Diversity of mitochondrial DNA in 3 species of great whales before and after modern whaling

The 20th century commercial whaling industry severely reduced populations of great whales throughout the Southern Hemisphere. The effect of this exploitation on genetic diversity and population structure remains largely undescribed. Here, we compare pre- and post-whaling diversity of mitochondrial DNA (mtDNA) control region sequences for 3 great whales in the South Atlantic, such as the blue, humpback, and fin whale. Pre-whaling diversity is described from mtDNA extracted from bones collected near abandoned whaling stations, primarily from the South Atlantic island of South Georgia. These bones are known to represent the first stage of 20th century whaling and thus pre-whaling diversity of these populations. Post-whaling diversity is described from previously published studies reporting large-scale sampling of living whales in the Southern Hemisphere. Despite relatively high levels of surviving genetic diversity in the post-whaling populations, we found evidence of a probable loss of mtDNA lineages in all 3 species. This is evidenced by the detection of a large number of haplotypes found in the pre-whaling samples that are not present in the post-whaling samples. A rarefaction analysis further supports a loss of haplotypes in the South Atlantic humpback and Antarctic blue whale populations. The bones from former whaling stations in the South Atlantic represent a remarkable molecular archive for further investigation of the decline and ongoing recovery in the great whales of the Southern Hemisphere.

High genomic diversity in the endangered East Greenland Svalbard Barents Sea stock of bowhead whales (Balaena mysticetus)

The East Greenland-Svalbard-Barents Sea (EGSB) bowhead whale stock (Balaena mysticetus) was hunted to near extinction and remains Endangered on the International Union of Conservation of Nature Red List. The intense, temporally extensive hunting pressure may have left the population vulnerable to other perturbations, such as environmental change. However, the lack of genomic baseline data renders it difficult to evaluate the impacts of various potential stressors on this stock. Twelve EGSB bowhead whales sampled in 2017/2018 were re-sequenced and mapped to a previously published draft genome. All individuals were unrelated and void of significant signs of inbreeding, with similar observed and expected homo- and heterozygosity levels. Despite the small population size, mean autosome-wide heterozygosity was 0.00102, which is higher than that of most mammals for which comparable estimates are calculated using the same parameters, and three times higher than a conspecific individual from the Eastern-Canada-West-Greenland bowhead whale stock. Demographic history analyses indicated a continual decrease of Ne from ca. 1.5 million to ca. 250,000 years ago, followed by a slight increase until ca. 100,000 years ago, followed by a rapid decrease in Ne between 50,000 and 10,000 years ago. These estimates are lower than previously suggested based on mitochondrial DNA, but suggested demographic patterns over time are similar.

Bowhead whales overwinter in the Amundsen Gulf and Eastern Beaufort Sea

The bowhead whale is the only baleen whale endemic to the Arctic and is well adapted to this environment. Bowheads live near the polar ice edge for much of the year and although sea ice dynamics are not the only driver of their annual migratory movements, it likely plays a key role. Given the intrinsic variability of open water and ice, one might expect bowhead migratory plasticity to be high and linked to this proximate environmental factor. Here, through a network of underwater passive acoustic recorders, we document the first known occurrence of bowheads overwintering in what is normally their summer foraging grounds in the Amundsen Gulf and eastern Beaufort Sea. The underlying question is whether this is the leading edge of a phenological shift in a species' migratory behaviour in an environment undergoing dramatic shifts due to climate change.

Mitogenomics and the genetic differentiation of contemporary Balaena mysticetus (Cetacea) from Svalbard

Full mitochondrial genomes were assembled for 12 recently sampled animals from the Svalbard bowhead whale (Balaena mysticetus) stock via high-throughput sequencing data, facilitating analysis of the demographic history of the population for the first time. The Svalbard population has retained noticeable amounts of mitochondrial genome diversity despite extreme historical harvest levels. Haplotype and nucleotide diversities were similar to those estimated earlier for other bowhead whale populations. The reconstructed demographic history was in accordance with a boom–bust scenario, combining a slight Pleistocene population growth 25 000–35 000 years ago and a Holocene decline. Employing a mutation rate of 3.418 × 10–8 substitutions per site per year, the time to the most recent common ancestor for the mitochondrial genomes of the contemporary Svalbard bowhead whales was estimated to be 68 782 (54 353–83 216) years before the present. Based on 370 bp fragments of the D-loop region, significant genetic differentiation was detected between all extant bowhead whale populations across the circumpolar Arctic. Thus, the Svalbard bowhead whales can be regarded as a population with its own genetic legacy.

Cycles of trans-Arctic dispersal and vicariance, and diversification of the amphi-boreal marine fauna

The amphi-boreal faunal element comprises closely related species and conspecific populations with vicarious distributions in the North Atlantic and North Pacific basins. It originated from an initial trans-Arctic dispersal in the Pliocene after the first opening of the Bering Strait, and subsequent inter-oceanic vicariance through the Pleistocene when the passage through the Arctic was severed by glaciations and low sea levels. Opportunities for further trans-Arctic dispersal have risen at times, however, and molecular data now expose more complex patterns of inter-oceanic affinities and dispersal histories. For a general view on the trans-Arctic dynamics and of the roles of potential dispersal-vicariance cycles in generating systematic diversity, we produced new phylogeographic data sets for amphi-boreal taxa in 21 genera of invertebrates and vertebrates, and combined them with similar published data sets of mitochondrial coding gene variation, adding up to 89 inter-oceanic comparisons involving molluscs, crustaceans, echinoderms, polychaetes, fishes and mammals. Only 39% of the cases correspond to a simple history of Pliocene divergence; in most taxonomical groups, the range of divergence estimates implies connections through the entire Pliocene-Pleistocene-Holocene time frame. Repeated inter-oceanic exchange was inferred for 23 taxa, and the latest connection was usually post-glacial. Such repeated invasions have usually led to secondary contacts and occasionally to widespread hybridization between the different invasion waves. Late- or post-glacial exchange was inferred in 46% of the taxa, stressing the importance of the relatively recent invasions to the current diversity in the North Atlantic. Individual taxa also showed complex idiosyncratic patterns and histories, and several instances of cryptic speciation were recognized. In contrast to a simple inter-oceanic vicariance scenario underlying amphi-boreal speciation, the data expose complex patterns of reinvasion and reticulation that complicate the interpretation of taxon boundaries in the region.

The endangered Spitsbergen bowhead whales' secrets revealed after hundreds of years in hiding

Spitsbergen's bowhead whales (Balaena mysticetus) were hunted to near extinction in the world's first commercial whaling enterprise; this population clearly remains threatened, but nothing is known about its distribution, making assessment unfeasible. In this study, we document range, movement patterns and habitat preferences of this population, based on tagging done from an icebreaker-based helicopter. Despite their reduced abundance, Spitsbergen's bowhead whales occupy much of their historical range, stretching across the northern Barents Region from East Greenland eastward to Franz Josef Land. Unlike larger bowhead populations to the west, they do not migrate in a classical sense, but rather disperse from wintering grounds in the northernmost parts of their range during spring, returning northward again in autumn, a pattern opposite in terms of directionality compared to other Arctic bowhead whale populations. The extreme affiliation of this population with cold, ice-filled waters is a concern given ongoing climate warming and concomitant rapid sea ice habitat loss.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean F_ST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

Effects of sea ice on Arctic biota: an emerging crisis discipline

The rapid decline in Arctic sea ice (ASI) extent, area and volume during recent decades is occurring before we can understand many of the mechanisms through which ASI interacts with biological processes both at sea and on land. As a consequence, our ability to predict and manage the effects of this enormous environmental change is limited, making this a crisis discipline Here, we propose a framework to study these effects, defining direct effects as those acting on life-history events of Arctic biota, and indirect effects, where ASI acts upon biological systems through chains of events, normally involving other components of the physical system and/or biotic interactions. Given the breadth and complexity of ASI's effects on Arctic biota, Arctic research requires a truly multidisciplinary approach to address this issue. In the absence of effective global efforts to tackle anthropogenic global warming, ASI will likely continue to decrease, compromising the conservation of many ASI-related taxonomic groups and ecosystems. Mitigation actions will rely heavily on the knowledge acquired on the mechanisms and components involved with the biological effects of ASI.

Incorporating non-equilibrium dynamics into demographic history inferences of a migratory marine species

Understanding how dispersal and gene flow link geographically separated the populations over evolutionary history is challenging, particularly in migratory marine species. In southern right whales (SRWs, Eubalaena australis), patterns of genetic diversity are likely influenced by the glacial climate cycle and recent history of whaling. Here we use a dataset of mitochondrial DNA (mtDNA) sequences (n = 1327) and nuclear markers (17 microsatellite loci, n = 222) from major wintering grounds to investigate circumpolar population structure, historical demography and effective population size. Analyses of nuclear genetic variation identify two population clusters that correspond to the South Atlantic and Indo-Pacific ocean basins that have similar effective breeder estimates. In contrast, all wintering grounds show significant differentiation for mtDNA, but no sex-biased dispersal was detected using the microsatellite genotypes. An approximate Bayesian computation (ABC) approach with microsatellite markers compared the scenarios with gene flow through time, or isolation and secondary contact between ocean basins, while modelling declines in abundance linked to whaling. Secondary-contact scenarios yield the highest posterior probabilities, implying that populations in different ocean basins were largely isolated and came into secondary contact within the last 25,000 years, but the role of whaling in changes in genetic diversity and gene flow over recent generations could not be resolved. We hypothesise that these findings are driven by factors that promote isolation, such as female philopatry, and factors that could promote dispersal, such as oceanographic changes. These findings highlight the application of ABC approaches to infer the connectivity in mobile species with complex population histories and, currently, low levels of differentiation.

Melting barriers to faunal exchange across ocean basins

Accelerated loss of sea ice in the Arctic is opening routes connecting the Atlantic and Pacific Oceans for longer periods each year. These changes may increase the ease and frequency with which marine birds and mammals move between the Pacific and Atlantic Ocean basins. Indeed, recent observations of birds and mammals suggest these movements have intensified in recent decades. Reconnection of the Pacific and Atlantic Ocean basins will present both challenges to marine ecosystem conservation and an unprecedented opportunity to examine the ecological and evolutionary consequences of interoceanic faunal exchange in real time. To understand these changes and implement effective conservation of marine ecosystems, we need to further develop modeling efforts to predict the rate of dispersal and consequences of faunal exchange. These predictions can be tested by closely monitoring wildlife dispersal through the Arctic Ocean and using modern methods to explore the ecological and evolutionary consequences of these movements.

Geographic and temporal dynamics of a global radiation and diversification in the killer whale

Global climate change during the Late Pleistocene periodically encroached and then released habitat during the glacial cycles, causing range expansions and contractions in some species. These dynamics have played a major role in geographic radiations, diversification and speciation. We investigate these dynamics in the most widely distributed of marine mammals, the killer whale (Orcinus orca), using a global data set of over 450 samples. This marine top predator inhabits coastal and pelagic ecosystems ranging from the ice edge to the tropics, often exhibiting ecological, behavioural and morphological variation suggestive of local adaptation accompanied by reproductive isolation. Results suggest a rapid global radiation occurred over the last 350 000 years. Based on habitat models, we estimated there was only a 15% global contraction of core suitable habitat during the last glacial maximum, and the resources appeared to sustain a constant global effective female population size throughout the Late Pleistocene. Reconstruction of the ancestral phylogeography highlighted the high mobility of this species, identifying 22 strongly supported long-range dispersal events including interoceanic and interhemispheric movement. Despite this propensity for geographic dispersal, the increased sampling of this study uncovered very few potential examples of ancestral dispersal among ecotypes. Concordance of nuclear and mitochondrial data further confirms genetic cohesiveness, with little or no current gene flow among sympatric ecotypes. Taken as a whole, our data suggest that the glacial cycles influenced local populations in different ways, with no clear global pattern, but with secondary contact among lineages following long-range dispersal as a potential mechanism driving ecological diversification.

Global diversity and oceanic divergence of humpback whales (Megaptera novaeangliae)

Humpback whales (Megaptera novaeangliae) annually undertake the longest migrations between seasonal feeding and breeding grounds of any mammal. Despite this dispersal potential, discontinuous seasonal distributions and migratory patterns suggest that humpbacks form discrete regional populations within each ocean. To better understand the worldwide population history of humpbacks, and the interplay of this species with the oceanic environment through geological time, we assembled mitochondrial DNA control region sequences representing approximately 2700 individuals (465 bp, 219 haplotypes) and eight nuclear intronic sequences representing approximately 70 individuals (3700 bp, 140 alleles) from the North Pacific, North Atlantic and Southern Hemisphere. Bayesian divergence time reconstructions date the origin of humpback mtDNA lineages to the Pleistocene (880 ka, 95% posterior intervals 550-1320 ka) and estimate radiation of current Northern Hemisphere lineages between 50 and 200 ka, indicating colonization of the northern oceans prior to the Last Glacial Maximum. Coalescent analyses reveal restricted gene flow between ocean basins, with long-term migration rates (individual migrants per generation) of less than 3.3 for mtDNA and less than 2 for nuclear genomic DNA. Genetic evidence suggests that humpbacks in the North Pacific, North Atlantic and Southern Hemisphere are on independent evolutionary trajectories, supporting taxonomic revision of M. novaeangliae to three subspecies.

High genetic diversity in a small population: the case of Chilean blue whales

It is generally assumed that species with low population sizes have lower genetic diversities than larger populations and vice versa. However, this would not be the case for long-lived species with long generation times, and which populations have declined due to anthropogenic effects, such as the blue whale (Balaenoptera musculus). This species was intensively decimated globally to near extinction during the 20th century. Along the Chilean coast, it is estimated that at least 4288 blue whales were hunted from an apparently pre-exploitation population size (k) of a maximum of 6200 individuals (Southeastern Pacific). Thus, here, we describe the mtDNA (control region) and nDNA (microsatellites) diversities of the Chilean blue whale aggregation site in order to verify the expectation of low genetic diversity in small populations. We then compare our findings with other blue whale aggregations in the Southern Hemisphere. Interestingly, although the estimated population size is small compared with the pre-whaling era, there is still considerable genetic diversity, even after the population crash, both in mitochondrial (N = 46) and nuclear (N = 52) markers (Hd = 0.890 and Ho = 0.692, respectively). Our results suggest that this diversity could be a consequence of the long generation times and the relatively short period of time elapsed since the end of whaling, which has been observed in other heavily-exploited whale populations. The genetic variability of blue whales on their southern Chile feeding grounds was similar to that found in other Southern Hemisphere blue whale feeding grounds. Our phylogenetic analysis of mtDNA haplotypes does not show extensive differentiation of populations among Southern Hemisphere blue whale feeding grounds. The present study suggests that although levels of genetic diversity are frequently used as estimators of population health, these parameters depend on the biology of the species and should be taken into account in a monitoring framework study to obtain a more complete picture of the conservation status of a population.

Molecular insights into the historic demography of bowhead whales: understanding the evolutionary basis of contemporary management practices

Patterns of genetic variation observed within species reflect evolutionary histories that include signatures of past demography. Understanding the demographic component of species' history is fundamental to informed management because changes in effective population size affect response to environmental change and evolvability, the strength of genetic drift, and maintenance of genetic variability. Species experiencing anthropogenic population reductions provide valuable case studies for understanding the genetic response to demographic change because historic changes in the census size are often well documented. A classic example is the bowhead whale, Balaena mysticetus, which experienced dramatic population depletion due to commercial whaling in the late 19th and early 20th centuries. Consequently, we analyzed a large multi-marker dataset of bowhead whales using a variety of analytical methods, including extended Bayesian skyline analysis and approximate Bayesian computation, to characterize genetic signatures of both ancient and contemporary demographic histories. No genetic signature of recent population depletion was recovered through any analysis incorporating realistic mutation assumptions, probably due to the combined influences of long generation time, short bottleneck duration, and the magnitude of population depletion. In contrast, a robust signal of population expansion was detected around 70,000 years ago, followed by a population decline around 15,000 years ago. The timing of these events coincides to a historic glacial period and the onset of warming at the end of the last glacial maximum, respectively. By implication, climate driven long-term variation in Arctic Ocean productivity, rather than recent anthropogenic disturbance, appears to have been the primary driver of historic bowhead whale demography.