Phylogenomics of the killer whale indicates ecotype divergence in sympatry

Are There Barriers Separating the Pink River Dolphin Populations (Inia boliviensis, Iniidae, Cetacea) within the Mamoré-Iténez River Basins (Bolivia)? An Analysis of Its Genetic Structure by Means of Mitochondrial and Nuclear DNA Markers

The pink river dolphin, or bufeo, is one of the dolphins which lives in the rivers of the Orinoco and Amazon basins in South America. The Bolivian bufeo population is considered a differentiated species (Inia boliviensis) from the Amazon and Orinoco species (Inia geoffrensis). Until now, no study has completed an extensive population genetics analysis of the bufeo in Bolivian rivers. We analyzed 82 bufeos from different rivers from the Mamoré and Iténez (Guaporé) river basins for the mt control region (CR), nuclear microsatellites, and DQB-1 gene sequences to determine if the inner rapids of these Bolivian river basins have some influence on the genetic structure of this species. The first relevant result was that the genetic diversity for CR, and the microsatellites were substantially lower in the Bolivian bufeos than in the dolphins studied in other areas of the Amazon and Orinoco. However, the DQB-1 gene sequences yielded similar genetic diversity to those found in other areas. The second relevant result is the existence of some significant genetic heterogeneity among the bufeo populations within Bolivia, although in a small degree, but this differentiation is independent of the inner rapids of the Bolivian rivers we sampled. The third relevant result was the existence of significant isolation by distance for the CR, but not for microsatellites and DQB-1 gene sequences. This was related to differential gene flow capacity of females (philopatric) and males (less philopatric and more migrants) and, possibly, to different selective patterns affecting the molecular markers studied. The fourth relevant result was related to diverse demographic changes of these bufeos. At least two or three bottleneck events and one or two population expansions have occurred in the Bolivian bufeo population. The major part of these events occurred during the Pleistocene.

Revised taxonomy of eastern North Pacific killer whales (Orcinus orca): Bigg's and resident ecotypes deserve species status

Killer whales (Orcinus orca) are currently recognized as a single ecologically and morphologically diverse, globally distributed species. Multiple morphotypes or ecotypes have been described, often associated with feeding specialization, and several studies have suggested taxonomic revision to include multiple subspecies or species in the genus. We review the ecological, morphological and genetic data for the well-studied 'resident' and Bigg's (aka 'transient') ecotypes in the eastern North Pacific and use quantitative taxonomic guidelines and standards to determine whether the taxonomic status of these killer whale ecotypes should be revised. Our review and new analyses indicate that species-level status is justified in both cases, and we conclude that eastern North Pacific Bigg's killer whales should be recognized as Orcinus rectipinnus (Cope in Scammon, 1869) and resident killer whales should be recognized as Orcinus ater (Cope in Scammon, 1869).

Killer whale acoustic patterns respond to prey abundance and environmental variability around the Prince Edward Islands, Southern Ocean

Killer whales are apex predators with temporally and spatially varying distributions throughout the world's oceans. Their ecology and behaviour are poorly understood in most regions due to limited research, often because of logistical challenges. Here, we used a passive acoustic monitoring device to investigate the seasonal acoustic occurrence and diel vocalizing behaviour of killer whales around the remote sub-Antarctic Prince Edward Islands (PEIs), Southern Ocean. Killer whales showed diel vocalizing patterns that varied seasonally in relation to their prey abundance and social activities. Killer whale calls were intermittently detected year-round with a high number of hours containing calls in October to December, and a secondary peak in February to May, corresponding to seal prey abundance. Random forest modelling identified wind speed as the primary predictor of the occurrence of killer whale calls (with a negative correlation) while sea surface height, chlorophyll-a and sea surface temperature were moderately important. We provide the first acoustic evidence that killer whale occurrence around the PEIs might coincide with variability in environmental conditions and prey abundance. Our results provide the first indication of diel vocalizing pattern of killer whales in the Southern Ocean. This knowledge is important for understanding killer whale ecology and adaptation to the changing oceans.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

Further evidence from common garden rearing experiments of heritable traits separating lean and siscowet lake charr (Salvelinus namaycush) ecotypes

Genetic evidence of selection for complex and polygenically regulated phenotypes can easily become masked by neutral population genetic structure and phenotypic plasticity. Without direct evidence of genotype‐phenotype associations it can be difficult to conclude to what degree a phenotype is heritable or a product of environment. Common garden laboratory studies control for environmental stochasticity and help to determine the mechanism that regulate traits. Here we assess lipid content, growth, weight, and length variation in full and hybrid F₁ crosses of deep and shallow water sympatric lake charr ecotypes reared for nine years in a common garden experiment. Redundancy analysis (RDA) and quantitative‐trait‐loci (QTL) genomic scans are used to identify associations between genotypes at 19,714 single nucleotide polymorphisms (SNPs) aligned to the lake charr genome and individual phenotypes to determine the role that genetic inheritance plays in ecotype phenotypic diversity. Lipid content, growth, length, and weight differed significantly among lake charr crosses throughout the experiment suggesting that pedigree plays a large roll in lake charr development. Polygenic scores of 15 SNPs putatively associated with lipid content and/or condition factor indicated that ecotype distinguishing traits are polygenically regulated and additive. A QTL identified on chromosome 38 contained >200 genes, some of which were associated with lipid metabolism and growth, demonstrating the complex nature of ecotype diversity. The results of our common garden study further indicate that lake charr ecotypes observed in nature are predetermined at birth and that ecotypes differ fundamentally in lipid metabolism and growth.

The relevance of genetic structure in ecotype designation and conservation management

The concept of ecotypes is complex, partly because of its interdisciplinary nature, but the idea is intrinsically valuable for evolutionary biology and applied conservation. The complex nature of ecotypes has spurred some confusion and inconsistencies in the literature, thereby limiting broader theoretical development and practical application. We provide suggestions for how incorporating genetic analyses can ease confusion and help define ecotypes. We approach this by systematically reviewing 112 publications across taxa that simultaneously mention the terms ecotype, conservation and management, to examine the current use of the term in the context of conservation and management. We found that most ecotype studies involve fish, mammals and plants with a focus on habitat use, which at 60% was the most common criterion used for categorization of ecotypes. Only 53% of the studies incorporated genetic analyses, and major discrepancies in available genomic resources among taxa could have contributed to confusion about the role of genetic structure in delineating ecotypes. Our results show that the rapid advances in genetic methods, also for nonmodel organisms, can help clarify the spatiotemporal distribution of adaptive and neutral genetic variation and their relevance to ecotype designations. Genetic analyses can offer empirical support for the ecotype concept and provide a timely measure of evolutionary potential, especially in changing environmental conditions. Genetic variation that is often difficult to detect, including polygenic traits influenced by small contributions from several genes, can be vital for adaptation to rapidly changing environments. Emerging ecotypes may signal speciation in progress, and findings from genome-enabled organisms can help clarify important selective factors driving ecotype development and persistence, and thereby improve preservation of interspecific genetic diversity. Incorporation of genetic analyses in ecotype studies will help connect evolutionary biology and applied conservation, including that of problematic groups such as natural hybrid organisms and urban or anthropogenic ecotypes.

Speciation-by-depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second-generation sequence data. Here, we revisit the case of two depth-segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coral Eunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth-segregated populations (Agaricia fragilis and Pocillopora damicornis) suggest divergence leading to depth-segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence-with-gene flow may account for much of the differentiation that separates closely related, depth-segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.

Concentrations and profiles of organochlorine contaminants in North Pacific resident and transient killer whale (Orcinus orca) populations

Organochlorine (OC) profiles have been used as chemical "fingerprints" to infer an animal's foraging area. North Pacific killer whale (Orcinus orca) populations are exposed to different levels and patterns of OCs based on their prey, distribution, and amount of time spent in a particular area. To characterize concentrations and profiles of OCs found in various populations of North Pacific killer whales, polychlorinated biphenyls (PCBs), including dioxin-like congeners, DDTs, and hexachlorobenzene (HCB), were measured in biopsy blubber samples of photo-identified resident (fish-eating) and transient (mammal-eating) killer whales collected from 1994 through 2002 from Russian Far East waters to the waters of the west coast of the United States, representing 10 populations. We compared blubber OC concentrations based on ecotype (resident vs. transient), sex and reproductive maturity, and geographic area. We also examined OC mixtures to determine if we could detect segregated geographical areas (foraging areas) among the six populations with sufficient sample sizes. Transients had significantly higher OC concentrations than residents and adult male whales had consistently higher OC levels compared to adult females, regardless of ecotype. Our OC profile findings indicate segregated foraging areas for the North Pacific killer whales, consistent with observations of their geographic distributions. Several potential health risks have also been associated with exposure to high levels of contaminants in top-level predators including reproductive impairment, immune suppression, skeletal deformities, and carcinoma. The results of this baseline study provide information on the geographic distribution of OCs found in North Pacific killer whales, results which are crucial for assessing the potential health risks associated with OC exposure in this species.

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.

Metabolic asymmetry and the global diversity of marine predators

Species richness of marine mammals and birds is highest in cold, temperate seas—a conspicuous exception to the general latitudinal gradient of decreasing diversity from the tropics to the poles. We compiled a comprehensive dataset for 998 species of sharks, fish, reptiles, mammals, and birds to identify and quantify inverse latitudinal gradients in diversity, and derived a theory to explain these patterns. We found that richness, phylogenetic diversity, and abundance of marine predators diverge systematically with thermoregulatory strategy and water temperature, reflecting metabolic differences between endotherms and ectotherms that drive trophic and competitive interactions. Spatial patterns of foraging support theoretical predictions, with total prey consumption by mammals increasing by a factor of 80 from the equator to the poles after controlling for productivity.

Killer whales (Orcinus orca) in Iceland show weak genetic structure among diverse isotopic signatures and observed movement patterns

Local adaption through ecological niche specialization can lead to genetic structure between and within populations. In the Northeast Pacific, killer whales (Orcinus orca) of the same population have uniform specialized diets that are non-overlapping with other sympatric, genetically divergent, and socially isolated killer whale ecotypes. However, killer whales in Iceland show intrapopulation variation of isotopic niches and observed movement patterns: some individuals appear to specialize on herring and follow it year-round while others feed upon herring only seasonally or opportunistically. We investigated genetic differentiation among Icelandic killer whales with different isotopic signatures and observed movement patterns. This information is key for management and conservation purposes but also for better understanding how niche specialization drives genetic differentiation. Photo-identified individuals (N = 61) were genotyped for 22 microsatellites and a 611 bp portion of the mitochondrial DNA (mtDNA) control region. Photo-identification of individuals allowed linkage of genetic data to existing data on individual isotopic niche, observed movement patterns, and social associations. Population subdivision into three genetic units was supported by a discriminant analysis of principal components (DAPC). Genetic clustering corresponded to the distribution of isotopic signatures, mtDNA haplotypes, and observed movement patterns, but genetic units were not socially segregated. Genetic differentiation was weak (F ST < 0.1), suggesting ongoing gene flow or recent separation of the genetic units. Our results show that killer whales in Iceland are not as genetically differentiated, ecologically discrete, or socially isolated as the Northeast Pacific prey-specialized killer whales. If any process of ecological divergence and niche specialization is taking place among killer whales in Iceland, it is likely at a very early stage and has not led to the patterns observed in the Northeast Pacific.

Colonizing the Wild West: Low Diversity of Complete Mitochondrial Genomes in Western North Pacific Killer Whales Suggests a Founder Effect

In the North Pacific, fish-eating R-type "resident" and mammal-eating T-type "transient" killer whales do not interbreed and differ in ecology and behavior. Full-length mitochondrial genomes (about 16.4 kbp) were sequenced and assembled for 12 R-type and 14 T-type killer whale samples from different areas of the western North Pacific. All R-type individuals had the same haplotype, previously described for R-type killer whales from both eastern and western North Pacific. However, haplotype diversity of R-type killer whales was much lower in the western North Pacific than in the Aleutian Islands and the eastern North Pacific. T-type whales had 3 different haplotypes, including one previously undescribed. Haplotype diversity of T-type killer whales in the Okhotsk Sea was also much lower than in the Aleutian Islands and the eastern North Pacific. The highest haplotype diversity for both R- and T-type killer whales was observed in the Aleutian Islands. We discuss how the environmental conditions during the last glacial period might have shaped the history of killer whale populations in the North Pacific. Our results suggest the recent colonization or re-colonization of the western North Pacific by small groups of killer whales originating from the central or eastern North Pacific, possibly due to favorable environmental changes after the Last Glacial Maximum.

Origins and wanderings of the Finnish hunting spitzes

Deducing the evolutionary histories of dog breeds can be challenging due to convergent traits and frequent admixture. In this report, we have explored the relationships of indigenous Finnish hunting spitz breeds among other northern Eurasian hunting breeds using commercially available SNP analysis (the MyDogDNA panel test). We find that Nordic hunting breeds Finnish Spitz, Nordic Spitz and the Karelian Bear Dog, as well as the reindeer herding Lapphund and Lapponian herder are all closely related and have common origins with the northeastern Eurasian Laika breeds, rather than with other Scandinavian Spitz breeds, such as Elkhounds and Swedish Vallhund. By tracing admixture events and direction of gene flow, we also elucidate the complex interactions between the breeds and provide new insight into the history of Swedish Elkhound and Russian-European Laika. The findings, together with an analysis of genetic differentiation between the populations, not only help to understand the origins of the breeds but also provide interesting possibilities to revive genetic diversity, lost during the breeding history, by backcrossing breeds to their hypothetical ancestry.

Eco-Evolutionary Processes Generating Diversity Among Bottlenose Dolphin, Tursiops truncatus, Populations off Baja California, Mexico

For highly mobile species that nevertheless show fine-scale patterns of population genetic structure, the relevant evolutionary mechanisms determining structure remain poorly understood. The bottlenose dolphin (Tursiops truncatus) is one such species, exhibiting complex patterns of genetic structure associated with local habitat dependence in various geographic regions. Here we studied bottlenose dolphin populations in the Gulf of California and Pacific Ocean off Baja California where habitat is highly structured to test associations between ecology, habitat dependence and genetic differentiation. We investigated population structure at a fine geographic scale using both stable isotope analysis (to assess feeding ecology) and molecular genetic markers (to assess population structure). Our results show that there are at least two factors affecting population structure for both genetics and feeding ecology (as indicated by stable isotope profiles). On the one hand there is a signal for the differentiation of individuals by ecotype, one foraging more offshore than the other. At the same time, there is differentiation between the Gulf of California and the west coast of Baja California, meaning that for example, nearshore ecotypes were both genetically and isotopically differentiated either side of the peninsula. We discuss these data in the context of similar studies showing fine-scale population structure for delphinid species in coastal waters, and consider possible evolutionary mechanisms.

Gene-culture coevolution in whales and dolphins

Whales and dolphins (Cetacea) have excellent social learning skills as well as a long and strong mother-calf bond. These features produce stable cultures, and, in some species, sympatric groups with different cultures. There is evidence and speculation that this cultural transmission of behavior has affected gene distributions. Culture seems to have driven killer whales into distinct ecotypes, which may be incipient species or subspecies. There are ecotype-specific signals of selection in functional genes that correspond to cultural foraging behavior and habitat use by the different ecotypes. The five species of whale with matrilineal social systems have remarkably low diversity of mtDNA. Cultural hitchhiking, the transmission of functionally neutral genes in parallel with selective cultural traits, is a plausible hypothesis for this low diversity, especially in sperm whales. In killer whales the ecotype divisions, together with founding bottlenecks, selection, and cultural hitchhiking, likely explain the low mtDNA diversity. Several cetacean species show habitat-specific distributions of mtDNA haplotypes, probably the result of mother-offspring cultural transmission of migration routes or destinations. In bottlenose dolphins, remarkable small-scale differences in haplotype distribution result from maternal cultural transmission of foraging methods, and large-scale redistributions of sperm whale cultural clans in the Pacific have likely changed mitochondrial genetic geography. With the acceleration of genomics new results should come fast, but understanding gene-culture coevolution will be hampered by the measured pace of research on the socio-cultural side of cetacean biology.

Intrapopulation variability in the timing of ontogenetic habitat shifts in sea turtles revealed using δ15 N values from bone growth rings

Determining location and timing of ontogenetic shifts in the habitat use of highly migratory species, along with possible intrapopulation variation in these shifts, is essential for understanding mechanisms driving alternate life histories and assessing overall population trends. Measuring variations in multi-year habitat-use patterns is especially difficult for remote oceanic species. To investigate the potential for differential habitat use among migratory marine vertebrates, we measured the naturally occurring stable nitrogen isotope (δ¹⁵ N) patterns that differentiate distinct ocean regions to create a 'regional isotope characterization', analysed the δ¹⁵ N values from annual bone growth layer rings from dead-stranded animals, and then combined the bone and regional isotope data to track individual animal movement patterns over multiple years. We used humeri from juvenile North Pacific loggerhead turtles (Caretta caretta), animals that undergo long migrations across the North Pacific Ocean (NPO), using multiple discrete regions as they develop to adulthood. Typical of many migratory marine species, ontogenetic changes in habitat use throughout their decades-long juvenile stage is poorly understood, but each potential habitat has unique foraging opportunities and spatially explicit natural and anthropogenic threats that could affect key life-history parameters. We found a bimodal size/age distribution in the timing that juveniles underwent an ontogenetic habitat shift from the oceanic central North Pacific (CNP) to the neritic east Pacific region near the Baja California Peninsula (BCP) (42·7 ± 7·2 vs. 68·3 ± 3·4 cm carapace length, 7·5 ± 2·7 vs. 15·6 ± 1·7 years). Important to the survival of this population, these disparate habitats differ considerably in their food availability, energy requirements and threats, and these differences can influence life-history parameters such as growth, survival and future fecundity. This is the first evidence of alternative ontogenetic shifts and habitat-use patterns for juveniles foraging in the eastern NPO. We combine two techniques, skeletochronology and stable isotope analysis, to reconstruct multi-year habitat-use patterns of a remote migratory species, linked to estimated ages and body sizes of individuals, to reveal variable ontogeny during the juvenile life stage that could drive alternate life histories and that has the potential to illuminate the migration patterns for other species with accretionary tissues.

High gene flow in the American badger overrides habitat preferences and limits broadscale genetic structure

Habitat associations are a function of habitat preferences and dispersal capabilities, both of which can influence how species responded to Quaternary climatic changes and contemporary habitat heterogeneity. Predicting resultant genetic structure is not always straightforward, especially in species where high dispersal potential and habitat preferences yield opposing predictions. The American badger has high dispersal capabilities that predict widespread panmixia, but avoids closed-canopy forests and clay soils, which could restrict gene flow and create ecologically based population genetic structure. We used mitochondrial sequence and microsatellite data sets to characterize how these opposing forces contribute to genetic structure in badgers at a continent-wide scale. Our data revealed an overall lack of ecologically based population genetic structure, suggesting that high dispersal capabilities were sufficiently realized to overcome most habitat-based genetic structure. At a broadscale, badger gene flow is limited only by geographic distance (isolation by distance) and large water barriers (Lake Michigan and the Mississippi River). The absence of genetic structure in a species with strong avoidance of unsuitable habitats advances our understanding of when and how genetic structure emerges in widespread, highly mobile species.

Genome-wide SNP data suggest complex ancestry of sympatric North Pacific killer whale ecotypes

Three ecotypes of killer whale occur in partial sympatry in the North Pacific. Individuals assortatively mate within the same ecotype, resulting in correlated ecological and genetic differentiation. A key question is whether this pattern of evolutionary divergence is an example of incipient sympatric speciation from a single panmictic ancestral population, or whether sympatry could have resulted from multiple colonisations of the North Pacific and secondary contact between ecotypes. Here, we infer multilocus coalescent trees from >1000 nuclear single-nucleotide polymorphisms (SNPs) and find evidence of incomplete lineage sorting so that the genealogies of SNPs do not all conform to a single topology. To disentangle whether uncertainty in the phylogenetic inference of the relationships among ecotypes could also result from ancestral admixture events we reconstructed the relationship among the ecotypes as an admixture graph and estimated f₄-statistics using TreeMix. The results were consistent with episodes of admixture between two of the North Pacific ecotypes and the two outgroups (populations from the Southern Ocean and the North Atlantic). Gene flow may have occurred via unsampled 'ghost' populations rather than directly between the populations sampled here. Our results indicate that because of ancestral admixture events and incomplete lineage sorting, a single bifurcating tree does not fully describe the relationship among these populations. The data are therefore most consistent with the genomic variation among North Pacific killer whale ecotypes resulting from multiple colonisation events, and secondary contact may have facilitated evolutionary divergence. Thus, the present-day populations of North Pacific killer whale ecotypes have a complex ancestry, confounding the tree-based inference of ancestral geography.

Genome-culture coevolution promotes rapid divergence of killer whale ecotypes

Analysing population genomic data from killer whale ecotypes, which we estimate have globally radiated within less than 250,000 years, we show that genetic structuring including the segregation of potentially functional alleles is associated with socially inherited ecological niche. Reconstruction of ancestral demographic history revealed bottlenecks during founder events, likely promoting ecological divergence and genetic drift resulting in a wide range of genome-wide differentiation between pairs of allopatric and sympatric ecotypes. Functional enrichment analyses provided evidence for regional genomic divergence associated with habitat, dietary preferences and post-zygotic reproductive isolation. Our findings are consistent with expansion of small founder groups into novel niches by an initial plastic behavioural response, perpetuated by social learning imposing an altered natural selection regime. The study constitutes an important step towards an understanding of the complex interaction between demographic history, culture, ecological adaptation and evolution at the genomic level.

Killer whales have evolved into specialized ecotypes based on hunting strategies and ecological niches. Here, Andrew Foote and colleagues sequenced the whole genome of individual killer whales representing 5 different ecotypes from North Pacific and Antarctic, and show expansion of small founder groups to adapt to specific ecological niches.

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.

Killer whale call frequency is similar across the oceans, but varies across sympatric ecotypes

Killer whale populations may differ in genetics, morphology, ecology, and behavior. In the North Pacific, two sympatric populations ("resident" and "transient") specialize on different prey (fish and marine mammals) and retain reproductive isolation. In the eastern North Atlantic, whales from the same populations have been observed feeding on both fish and marine mammals. Fish-eating North Pacific "residents" are more genetically related to eastern North Atlantic killer whales than to sympatric mammal-eating "transients." In this paper, a comparison of frequency variables in killer whale calls recorded from four North Pacific resident, two North Pacific transient, and two eastern North Atlantic populations is reported to assess which factors drive the large-scale changes in call structure. Both low-frequency and high-frequency components of North Pacific transient killer whale calls have significantly lower frequencies than those of the North Pacific resident and North Atlantic populations. The difference in frequencies could be related to ecological specialization or to the phylogenetic history of these populations. North Pacific transient killer whales may have genetically inherited predisposition toward lower frequencies that may shape their learned repertoires.

Population genomics of the killer whale indicates ecotype evolution in sympatry involving both selection and drift

The evolution of diversity in the marine ecosystem is poorly understood, given the relatively high potential for connectivity, especially for highly mobile species such as whales and dolphins. The killer whale (Orcinus orca) has a worldwide distribution, and individual social groups travel over a wide geographic range. Even so, regional populations have been shown to be genetically differentiated, including among different foraging specialists (ecotypes) in sympatry. Given the strong matrifocal social structure of this species together with strong resource specializations, understanding the process of differentiation will require an understanding of the relative importance of both genetic drift and local adaptation. Here we provide a high-resolution analysis based on nuclear single-nucleotide polymorphic markers and inference about differentiation at both neutral loci and those potentially under selection. We find that all population comparisons, within or among foraging ecotypes, show significant differentiation, including populations in parapatry and sympatry. Loci putatively under selection show a different pattern of structure compared to neutral loci and are associated with gene ontology terms reflecting physiologically relevant functions (e.g. related to digestion). The pattern of differentiation for one ecotype in the North Pacific suggests local adaptation and shows some fixed differences among sympatric ecotypes. We suggest that differential habitat use and resource specializations have promoted sufficient isolation to allow differential evolution at neutral and functional loci, but that the process is recent and dependent on both selection and drift.