Runs of homozygosity (ROH) occur when offspring receive the same ancestral haplotype from both parents, and, accordingly, reduce individual heterozygosity. Their distribution throughout the genome contains information on the probability of inbreeding mediated by mating system and population demography. Here, we investigate variation in killer whale demographic history as reflected in genome-wide heterozygosity, using a global dataset of 26 genomes. We find an overall pattern of lower heterozygosity in genomes sampled at high latitudes, with hundreds of short ROH (< 1Mbp) reflecting high background relatedness due to coalescence of haplotypes during bottlenecks associated with founder events during post-glacial range expansions. Across most of the species' range, intermediate length ROH (1-10Mb) revealed long-term inbreeding in 22 of the 26 sampled killer whale genomes, consistent with the high social philopatry observed in all populations studied to date. Inbreeding coefficients (FROH) were comparable to those reported in other taxa with long-term low population size, such as bonobos and the Native American Karitiana of the Brazilian Amazon. The extreme outlier in this dataset, a Scottish killer whale, was homozygous over one-third of the autosomes (41.6%) with a distinct distribution of ROH length, indicating generations of inbreeding. This exceeds autozygosity in emblematic examples of long-term inbreeding, such as the Altai Neanderthal, and eastern lowland and mountain gorillas. The fate of this Scottish killer whale population, in which no calves have been born in over two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.3 Species ranges change in response to environmental oscillation [1,2], and rapid shifts can currently be observed during ongoing global warming [3,4]. Understanding how range shifts influence the genetic diversity of natural populations at the range edge is an increasingly important conservation consideration [5]. Killer whales (Orcinus orca) are comparable to humans in their global distribution, having colonised all the major oceans [6]. Killer whale occurrence is correlated with ocean productivity; highest densities are therefore at high latitudes, decreasing by 1-2 orders of magnitude from the Arctic and Antarctic to the tropics ( Figure S1) [6]. The Northern and Southern extremes of this range were covered by ice sheets during the Last Glacial Maximum ( Figure 1) and must therefore have been colonised through range expansion from lower latitudes in the last 16 Kyr [7][8][9]. Independent post-glacial range expansions in different ocean basins offer the possibility to explore the genetic outcomes of independent, parallel demographic histories. Killer whales thus represent a useful system for studying the relationship between demographic history, range expansions, and genetic diversity. Here, we report the analyses of whole genome sequences from a global dataset of killer whales (Figure 1).
RESULTS
Heterozygosity Correlates with Latitude...