No distinct local cuisines among humpback whales: A population diet comparison in the Southern Hemisphere.

Southern hemisphere humpback whale (Megaptera novaeangliae, SHHW) breeding populations follow a high-fidelity Antarctic krill (Euphausia superba) diet while feeding in distinct sectors of the Southern Ocean. Their capital breeding life history requires predictable ecosystem productivity to fuel migration and migration-related behaviours. It is therefore postulated that populations feeding in areas subject to the strongest climate change impacts are more likely to show the first signs of a departure from a high-fidelity krill diet. We tested this hypothesis by investigating blubber fatty acid profiles and skin stable isotopes obtained from five SHHW populations in 2019, and comparing them to Antarctic krill stable isotopes sampled in three SHHW feeding areas in the Southern Ocean in 2019. Fatty acid profiles and δ13C and δ15N varied significantly among all five populations, however, calculated trophic positions did not (2.7 to 3.1). Similarly, fatty acid ratios, 16:1ω7c/16:0 and 20:5ω3/22:6ω3 were above 1, showing that whales from all five populations are secondary heterotrophs following an omnivorous diet with a diatom-origin. Thus, evidence for a potential departure from a high-fidelity Antarctic krill diet was not seen in any population. δ13C of all populations were similar to δ13C of krill sampled in productive upwelling areas or the marginal sea-ice zone. Consistency in trophic position and diet origin but significant fatty acid and stable isotope differences demonstrate that the observed variability arises at lower trophic levels. Our results indicate that, at present, there is no evidence of a divergence from a high-fidelity krill diet. Nevertheless, the characteristic isotopic signal of whales feeding in productive upwelling areas, or in the marginal sea-ice zone, implies that future cryosphere reductions could impact their feeding ecology.

Ancestry testing of "Old Tom," a killer whale central to mutualistic interactions with human whalers.

Cooperative hunting between humans and killer whales (Orcinus orca) targeting baleen whales was reported in Eden, New South Wales, Australia, for almost a century. By 1928, whaling operations had ceased, and local killer whale sightings became scarce. A killer whale from the group, known as "Old Tom," washed up dead in 1930 and his skeleton was preserved. How these killer whales from Eden relate to other populations globally and whether their genetic descendants persist today remains unknown. We extracted and sequenced DNA from Old Tom using ancient DNA techniques. Genomic sequences were then compared with a global dataset of mitochondrial and nuclear genomes. Old Tom shared a most recent common ancestor with killer whales from Australasia, the North Atlantic, and the North Pacific, having the highest genetic similarity with contemporary New Zealand killer whales. However, much of the variation found in Old Tom's genome was not shared with these widespread populations, suggesting ancestral rather than ongoing gene flow. Our genetic comparisons also failed to find any clear descendants of Tom, raising the possibility of local extinction of this group. We integrated Traditional Custodian knowledge to recapture the events in Eden and recognize that Indigenous Australians initiated the relationship with the killer whales before European colonization and the advent of commercial whaling locally. This study rectifies discrepancies in local records and provides new insight into the origins of the killer whales in Eden and the history of Australasian killer whales.

"Type D" killer whale genomes reveal long-term small population size and low genetic diversity.

Genome sequences can reveal the extent of inbreeding in small populations. Here, we present the first genomic characterization of type D killer whales, a distinctive eco/morphotype with a circumpolar, subantarctic distribution. Effective population size is the lowest estimated from any killer whale genome and indicates a severe population bottleneck. Consequently, type D genomes show among the highest level of inbreeding reported for any mammalian species (FROH ≥ 0.65). Detected recombination cross-over events of different haplotypes are up to an order of magnitude rarer than in other killer whale genomes studied to date. Comparison of genomic data from a museum specimen of a type D killer whale that stranded in New Zealand in 1955, with 3 modern genomes from the Cape Horn area, reveals high covariance and identity-by-state of alleles, suggesting these genomic characteristics and demographic history are shared among geographically dispersed social groups within this morphotype. Limitations to the insights gained in this study stem from the nonindependence of the 3 closely related modern genomes, the recent coalescence time of most variation within the genomes, and the nonequilibrium population history which violates the assumptions of many model-based methods. Long-range linkage disequilibrium and extensive runs of homozygosity found in type D genomes provide the potential basis for both the distinctive morphology, and the coupling of genetic barriers to gene flow with other killer whale populations.

Runs of homozygosity in killer whale genomes provide a global record of demographic histories.

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically, the probability of inbreeding mediated by mating system and/or population demography. Here, we investigated whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global data set of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstructed demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We found a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1 Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5 Mb) were found in low latitude populations, and populations of known conservation concern. These include a Scottish killer whale, for which 37.8% of the autosomes were comprised of ROH >1.5 Mb in length. The fate of this population, in which only two adult males have been sighted in the past five years, and zero fecundity over the last two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.

Killer whale genomes reveal a complex history of recurrent admixture and vicariance.

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.