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

Abstract

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...

Conservation genomics of the endangered Seychelles Magpie‐Robin (Copsychus sechellarum): a unique insight into the history of a precious endemic bird

“…Only individuals with an average genome‐wide coverage of 5× and above were included in the analysis of runs of homozygosity (ROHs) (136/141) as in Hooper et al . (2020). The largest 30 scaffolds for each of these individuals were selected for ROH analysis defined by a minimum length of 10 Mb.…”

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

“…However, a key insight from ROH studies has been that a large proportion of deleterious homozygous genotypes seem to fall into ROH (Szpiech et al, 2013), leading to the prediction that ROH detection could be one of the most powerful methods to investigate inbreeding depression (Keller et al, 2011). As genomic technologies have become more affordable, the opportunity to infer individual ROH in animals of high ecological importance or of conservation concern has therefore received increasing attention (Escoda & Castresana, 2021; Grossen et al, 2018; Humble et al, 2020; Hooer et al, 2020; Wang et al, 2021). For example, a large‐scale collared flycatcher ( Ficedula albicollis ) study on genomic inbreeding and historical demography found ROH sizes between 953 bp and 17.5 Mbp (Kardos et al, 2017), with highest ROH abundance in regions of low recombination rate.…”

Climate change and cetacean health: impacts and future directions

“…Thus, genetic diversity is also threatened under climate change. The decline in population size allows for inbreeding depression by reducing the genetic diversity and increasing recessive homozygotes created by consanguineous mating [78][79][80]. That would lead to a weakened resilience of populations against climate change effects and other anthropogenic pressures [81].…”