Social and genetic connectivity despite ecological variation in a killer whale network.

Philopatric kin-based societies encourage a narrow breadth of conservative behaviours owing to individuals primarily learning from close kin, promoting behavioural homogeneity. However, weaker social ties beyond kin, and across a behaviourally diverse social landscape, could be sufficient to induce variation and a greater ecological niche breadth. We investigated a network of 457 photo-identified killer whales from Norway (548 encounters in 2008-2021) with diet data available (46 mixed-diet individuals feeding on both fish and mammals, and 411 exclusive fish-eaters) to quantify patterns of association within and between diet groups, and to identify underlying correlates. We genotyped a subset of 106 whales to assess patterns of genetic differentiation. Our results suggested kinship as main driver of social bonds within and among cohesive social units, while diet was most likely a consequence reflective of cultural diffusion, rather than a driver. Flexible associations within and between ecologically diverse social units led to a highly connected network, reducing social and genetic differentiation between diet groups. Our study points to a role of social connectivity, in combination with individual behavioural variation, in influencing population ecology in killer whales.

Diversity and community structure of anaerobic gut fungi in the rumen of wild and domesticated herbivores.

The rumen houses a diverse community that plays a major role in the digestion process in ruminants. Anaerobic gut fungi (AGF) are key contributors to plant digestion in the rumen. Here, we present a global amplicon-based survey of the rumen AGF mycobiome by examining 206 samples from 15 animal species, 15 countries, and 6 continents. The rumen AGF mycobiome was highly diverse, with 81 out of 88 currently recognized AGF genera or candidate genera identified. However, only six genera (Neocallimastix, Orpinomyces, Caecomyces, Cyllamyces, NY9, and Piromyces) were present at >4% relative abundance. AGF diversity was higher in members of the families Antilocapridae and Cervidae compared to Bovidae. Community structure analysis identified a pattern of phylosymbiosis, where host family (10% of total variance) and species (13.5%) partially explained the rumen mycobiome composition. As well, diet composition (9%-19%), domestication (11.14%), and biogeography (14.1%) also partially explained AGF community structure; although sampling limitation, geographic range restrictions, and direct association between different factors hindered accurate elucidation of the relative contribution of each factor. Pairwise comparison of rumen and fecal samples obtained from the same subject (n = 13) demonstrated greater diversity and inter-sample variability in rumen versus fecal samples. The genera Neocallimastix and Orpinomyces were present in higher abundance in rumen samples, while Cyllamyces and Caecomyces were enriched in fecal samples. Comparative analysis of global rumen and feces data sets revealed a similar pattern. Our results provide a global view of AGF community in the rumen and identify patterns of AGF variability between rumen and feces in herbivores Gastrointestinal (GI) tract.IMPORTANCERuminants are highly successful and economically important mammalian suborder. Ruminants are herbivores that digest plant material with the aid of microorganisms residing in their GI tract. In ruminants, the rumen compartment represents the most important location where microbially mediated plant digestion occurs, and is known to house a bewildering array of microbial diversity. An important component of the rumen microbiome is the anaerobic gut fungi (AGF), members of the phylum Neocallimastigomycota. So far, studies examining AGF diversity have mostly employed fecal samples, and little is currently known regarding the identity of AGF residing in the rumen compartment, factors that impact the observed patterns of diversity and community structure of AGF in the rumen, and how AGF communities in the rumen compare to AGF communities in feces. Here, we examined the rumen AGF diversity using an amplicon-based survey targeting a wide range of wild and domesticated ruminants (n = 206, 15 different animal species) obtained from 15 different countries. Our results demonstrate that while highly diverse, no new AGF genera were identified in the rumen mycobiome samples examined. Our analysis also indicate that animal host phylogeny, diet, biogeography, and domestication status could play a role in shaping AGF community structure. Finally, we demonstrate that a greater level of diversity and higher inter-sample variability was observed in rumen compared to fecal samples, with two genera (Neocallimastix and Orpinomyces) present in higher abundance in rumen samples, and two others (Cyllamyces and Caecomyces) enriched in fecal samples. Our results provide a global view of the identity, diversity, and community structure of AGF in ruminants, elucidate factors impacting diversity and community structure of the rumen mycobiome, and identify patterns of AGF community variability between the rumen and feces in the herbivorous GI tract.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse.

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Climate change introduces threatened killer whale populations and conservation challenges to the Arctic.

The Arctic is the fastest-warming region on the planet, and the lengthening ice-free season is opening Arctic waters to sub-Arctic species such as the killer whale (Orcinus orca). As apex predators, killer whales can cause significant ecosystem-scale changes. Setting conservation priorities for killer whales and their Arctic prey species requires knowledge of their evolutionary history and demographic trajectory. Using whole-genome resequencing of 24 killer whales sampled in the northwest Atlantic, we first explored the population structure and demographic history of Arctic killer whales. To better understand the broader geographic relationship of these Arctic killer whales to other populations, we compared them to a globally sampled dataset. Finally, we assessed threats to Arctic killer whales due to anthropogenic harvest by reviewing the peer-reviewed and gray literature. We found that there are two highly genetically distinct, non-interbreeding populations of killer whales using the eastern Canadian Arctic. These populations appear to be as genetically different from each other as are ecotypes described elsewhere in the killer whale range; however, our data cannot speak to ecological differences between these populations. One population is newly identified as globally genetically distinct, and the second is genetically similar to individuals sampled from Greenland. The effective sizes of both populations recently declined, and both appear vulnerable to inbreeding and reduced adaptive potential. Our survey of human-caused mortalities suggests that harvest poses an ongoing threat to both populations. The dynamic Arctic environment complicates conservation and management efforts, with killer whales adding top-down pressure on Arctic food webs crucial to northern communities' social and economic well-being. While killer whales represent a conservation priority, they also complicate decisions surrounding wildlife conservation and resource management in the Arctic amid the effects of climate change.

Low-Protein Diets Differentially Regulate Energy Balance during Thermoneutral and Heat Stress in Cobb Broiler Chicken (Gallus domesticus).

The objective was to assess whether low-protein (LP) diets regulate food intake (FI) and thermogenesis differently during thermoneutral (TN) and heat stress (HS) conditions. Two-hundred-day-old male broiler chicks were weight-matched and assigned to 36 pens with 5-6 chicks/pen. After 2 weeks of acclimation, birds were subjected into four groups (9 pens/group) including (1) a normal-protein diet under TN (ambient temperature), (2) an LP diet under TN, (3) a normal-protein diet under HS (35 °C for 7 h/day), and (4) an LP diet under HS, for 4 weeks. During HS, but not TN, LP tended to decrease FI, which might be associated with a lower mRNA abundance of duodenal ghrelin and higher GIP during HS. The LP group had a higher thermal radiation than NP under TN, but during HS, the LP group had a lower thermal radiation than NP. This was linked with higher a transcript of muscle ß1AR and AMPKα1 during TN, but not HS. Further, LP increased the gene expression of COX IV during TN but reduced COX IV and the sirtuin 1 abundance during HS. The dietary protein content differentially impacted plasma metabolome during TN and HS with divergent changes in amino acids such as tyrosine and tryptophan. Compared to NP, LP had increased abundances of p_Tenericutes, c_Mollicutes, c_Mollicutes_RF9, and f_tachnospiraceae under HS. Overall, LP diets may mitigate the negative outcome of heat stress on the survivability of birds by reducing FI and heat production. The differential effect of an LP diet on energy balance during TN and HS is likely regulated by gut and skeletal muscle and alterations in plasma metabolites and cecal microbiota.

Temporal progression of anaerobic fungal communities in dairy calves from birth to maturity.

Establishment of microbial communities in neonatal calves is vital for their growth and overall health. While this process has received considerable attention for bacteria, our knowledge on temporal progression of anaerobic gut fungi (AGF) in calves is lacking. Here, we examined AGF communities in faecal samples from six dairy cattle collected at 24 different time points during the pre-weaning (days 1-48), weaning (days 48-60), and post-weaning (days 60-360) phases. Quantitative polymerase chain reaction indicated that AGF colonisation occurs within 24 h after birth, with loads slowly increasing during pre-weaning and weaning, then drastically increasing post-weaning. Culture-independent amplicon surveys identified higher alpha diversity during pre-weaning/weaning, compared to post-weaning. AGF community structure underwent a drastic shift post-weaning, from a community enriched in genera commonly encountered in hindgut fermenters to one enriched in genera commonly encountered in adult ruminants. Comparison of AGF community between calves day 1 post-birth and their mothers suggest a major role for maternal transmission, with additional input from cohabitating subjects. This distinct pattern of AGF progression could best be understood in-light of their narrower niche preferences, metabolic specialisation, and physiological optima compared to bacteria, hence eliciting a unique response to changes in feeding pattern and associated structural GIT development during maturation.

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.

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution.

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.

The life aquatic: advances in marine vertebrate genomics.

The ocean is hypothesized to be where life on earth originated, and subsequent evolutionary transitions between marine and terrestrial environments have been key events in the origin of contemporary biodiversity. Here, we review how comparative genomic approaches are an increasingly important aspect of understanding evolutionary processes, such as physiological and morphological adaptation to the diverse habitats within the marine environment. In addition, we highlight how population genomics has provided unprecedented resolution for population structuring, speciation and adaptation in marine environments, which can have a low cost of dispersal and few physical barriers to gene flow, and can thus support large populations. Building upon this work, we outline the applications of genomics tools to conservation and their relevance to assessing the wide-ranging impact of fisheries and climate change on marine species.

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.

Admixture between Ancient Lineages, Selection, and the Formation of Sympatric Stickleback Species-Pairs.

Ecological speciation has become a popular model for the development and maintenance of reproductive isolation in closely related sympatric pairs of species or ecotypes. An implicit assumption has been that such pairs originate (possibly with gene flow) from a recent, genetically homogeneous ancestor. However, recent genomic data have revealed that currently sympatric taxa are often a result of secondary contact between ancestrally allopatric lineages. This has sparked an interest in the importance of initial hybridization upon secondary contact, with genomic reanalysis of classic examples of ecological speciation often implicating admixture in speciation. We describe a novel occurrence of unusually well-developed reproductive isolation in a model system for ecological speciation: the three-spined stickleback (Gasterosteus aculeatus), breeding sympatrically in multiple lagoons on the Scottish island of North Uist. Using morphological data, targeted genotyping, and genome-wide single-nucleotide polymorphism data, we show that lagoon resident and anadromous ecotypes are strongly reproductively isolated with an estimated hybridization rate of only ∼1%. We use palaeoecological and genetic data to test three hypotheses to explain the existence of these species-pairs. Our results suggest that recent, purely ecological speciation from a genetically homogeneous ancestor is probably not solely responsible for the evolution of species-pairs. Instead, we reveal a complex colonization history with multiple ancestral lineages contributing to the genetic composition of species-pairs, alongside strong disruptive selection. Our results imply a role for admixture upon secondary contact and are consistent with the recent suggestion that the genomic underpinning of ecological speciation often has an older, allopatric origin.

Phylogenomics and species delimitation for effective conservation of manta and devil rays.

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Genomics and the origin of species.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Host-derived population genomics data provides insights into bacterial and diatom composition of the killer whale skin.

Recent exploration into the interactions and relationship between hosts and their microbiota has revealed a connection between many aspects of the host's biology, health and associated micro-organisms. Whereas amplicon sequencing has traditionally been used to characterize the microbiome, the increasing number of published population genomics data sets offers an underexploited opportunity to study microbial profiles from the host shotgun sequencing data. Here, we use sequence data originally generated from killer whale Orcinus orca skin biopsies for population genomics, to characterize the skin microbiome and investigate how host social and geographical factors influence the microbial community composition. Having identified 845 microbial taxa from 2.4 million reads that did not map to the killer whale reference genome, we found that both ecotypic and geographical factors influence community composition of killer whale skin microbiomes. Furthermore, we uncovered key taxa that drive the microbiome community composition and showed that they are embedded in unique networks, one of which is tentatively linked to diatom presence and poor skin condition. Community composition differed between Antarctic killer whales with and without diatom coverage, suggesting that the previously reported episodic migrations of Antarctic killer whales to warmer waters associated with skin turnover may control the effects of potentially pathogenic bacteria such as Tenacibaculum dicentrarchi. Our work demonstrates the feasibility of microbiome studies from host shotgun sequencing data and highlights the importance of metagenomics in understanding the relationship between host and microbial ecology.

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.

Postglacial Colonization of Northern Coastal Habitat by Bottlenose Dolphins: A Marine Leading-Edge Expansion?

Oscillations in the Earth's temperature and the subsequent retreating and advancing of ice-sheets around the polar regions are thought to have played an important role in shaping the distribution and genetic structuring of contemporary high-latitude populations. After the Last Glacial Maximum (LGM), retreating of the ice-sheets would have enabled early colonizers to rapidly occupy suitable niches to the exclusion of other conspecifics, thereby reducing genetic diversity at the leading-edge. Bottlenose dolphins (genus Tursiops) form distinct coastal and pelagic ecotypes, with finer-scale genetic structuring observed within each ecotype. We reconstruct the postglacial colonization of the Northeast Atlantic (NEA) by bottlenose dolphins using habitat modeling and phylogenetics. The AquaMaps model hindcasted suitable habitat for the LGM in the Atlantic lower latitude waters and parts of the Mediterranean Sea. The time-calibrated phylogeny, constructed with 86 complete mitochondrial genomes including 30 generated for this study and created using a multispecies coalescent model, suggests that the expansion to the available coastal habitat in the NEA happened via founder events starting ~15 000 years ago (95% highest posterior density interval: 4 900-26 400). The founders of the 2 distinct coastal NEA populations comprised as few as 2 maternal lineages that originated from the pelagic population. The low effective population size and genetic diversity estimated for the shared ancestral coastal population subsequent to divergence from the pelagic source population are consistent with leading-edge expansion. These findings highlight the legacy of the Late Pleistocene glacial cycles on the genetic structuring and diversity of contemporary populations.

Energy balance affects pulsatile secretion of luteinizing hormone from the adenohypophesis and expression of neurokinin B in the hypothalamus of ovariectomized gilts.

The pubertal transition of gonadotropin secretion in pigs is metabolically gated. Kisspeptin (KISS1) and neurokinin B (NKB) are coexpressed in neurons within the arcuate nucleus of the hypothalamus (ARC) and are thought to play an important role in the integration of nutrition and metabolic state with the reproductive neuroendocrine axis. The hypothesis that circulating concentrations of luteinizing hormone (LH) and expression of KISS1 and tachykinin 3(TAC3, encodes NKB) in the ARC of female pigs are reduced with negative energy balance was tested using ovariectomized, prepubertal gilts fed to either gain or lose body weight. Restricted feeding of ovariectomized gilts caused a rapid and sustained metabolic response characterized by reduced concentrations of plasma urea nitrogen, insulin, leptin, and insulin-like growth factor-1 and elevated concentrations of free fatty acids. The secretory pattern of LH shifted from one of low amplitude to one of high amplitude, which caused overall circulating concentrations of LH to be greater in restricted gilts. Nutrient-restricted gilts had greater expression of follicle-stimulating hormone and gonadotropin-releasing hormone receptor, but not LH in the anterior pituitary gland. Expression of KISS1 in the ARC was not affected by dietary treatment, but expression of TAC3 was greater in restricted gilts. These data are consistent with the idea that hypothalamic expression of KISS1 is correlated with the number of LH pulse in pig, and further indicate that amplitude of LH pulses may be regulated by NKB in the gilt.

Demography or selection on linked cultural traits or genes? Investigating the driver of low mtDNA diversity in the sperm whale using complementary mitochondrial and nuclear genome analyses.

Mitochondrial DNA has been heavily utilized in phylogeography studies for several decades. However, underlying patterns of demography and phylogeography may be misrepresented due to coalescence stochasticity, selection, variation in mutation rates and cultural hitchhiking (linkage of genetic variation to culturally-transmitted traits affecting fitness). Cultural hitchhiking has been suggested as an explanation for low genetic diversity in species with strong social structures, counteracting even high mobility, abundance and limited barriers to dispersal. One such species is the sperm whale, which shows very limited phylogeographic structure and low mtDNA diversity despite a worldwide distribution and large population. Here, we use analyses of 175 globally distributed mitogenomes and three nuclear genomes to evaluate hypotheses of a population bottleneck/expansion vs. a selective sweep due to cultural hitchhiking or selection on mtDNA as the mechanism contributing to low worldwide mitochondrial diversity in sperm whales. In contrast to mtDNA control region (CR) data, mitogenome haplotypes are largely ocean-specific, with only one of 80 shared between the Atlantic and Pacific. Demographic analyses of nuclear genomes suggest low mtDNA diversity is consistent with a global reduction in population size that ended approximately 125,000 years ago, correlated with the Eemian interglacial. Phylogeographic analysis suggests that extant sperm whales descend from maternal lineages endemic to the Pacific during the period of reduced abundance and have subsequently colonized the Atlantic several times. Results highlight the apparent impact of past climate change, and suggest selection and hitchhiking are not the sole processes responsible for low mtDNA diversity in this highly social species.