Collecting baleen whale blow samples by drone: A minimally intrusive tool for conservation genetics.

In coastal British Columbia, Canada, marine megafauna such as humpback whales (Megaptera novaeangliae) and fin whales (Balaenoptera physalus velifera) have been subject to a history of exploitation and near extirpation. While their populations have been in recovery, significant threats are posed to these vulnerable species by proposed natural resource ventures in this region, in addition to the compounding effects of anthropogenic climate change. Genetic tools play a vital role in informing conservation efforts, but the associated collection of tissue biopsy samples can be challenging for the investigators and disruptive to the ongoing behaviour of the targeted whales. Here, we evaluate a minimally intrusive approach based on collecting exhaled breath condensate, or respiratory 'blow' samples, from baleen whales using an unoccupied aerial system (UAS), within Gitga'at First Nation territory for conservation genetics. Minimal behavioural responses to the sampling technique were observed, with no response detected 87% of the time (of 112 UAS deployments). DNA from whale blow (n = 88 samples) was extracted, and DNA profiles consisting of 10 nuclear microsatellite loci, sex identification and mitochondrial (mt) DNA haplotypes were constructed. An average of 7.5 microsatellite loci per individual were successfully genotyped. The success rates for mtDNA and sex assignment were 80% and 89% respectively. Thus, this minimally intrusive sampling method can be used to describe genetic diversity and generate genetic profiles for individual identification. The results of this research demonstrate the potential of UAS-collected whale blow for conservation genetics from a remote location.

Using epigenetic clocks to investigate changes in the age structure of critically endangered Maui dolphins.

The age of an individual is an essential demographic parameter but is difficult to estimate without long-term monitoring or invasive sampling. Epigenetic approaches are increasingly used to age organisms, including nonmodel organisms such as cetaceans. Maui dolphins (Cephalorhynchus hectori maui) are a critically endangered subspecies endemic to Aotearoa New Zealand, and the age structure of this population is important for informing conservation. Here we present an epigenetic clock for aging Maui and Hector's dolphins (C. h. hectori) developed from methylation data using DNA from tooth aged individuals (n = 48). Based on this training data set, the optimal model required only eight methylation sites, provided an age correlation of .95, and had a median absolute age error of 1.54 years. A leave-one-out cross-validation analysis with the same parameters resulted in an age correlation of .87 and median absolute age error of 2.09 years. To improve age estimation, we included previously published beluga whale (Delphinapterus leucas) data to develop a joint beluga/dolphin clock, resulting in a clock with comparable performance and improved estimation of older individuals. Application of the models to DNA from skin biopsy samples of living Maui dolphins revealed a shift from a median age of 8-9 years to a younger population aged 7-8 years 10 years later. These models could be applied to other dolphin species and demonstrate the ability to construct a clock even when the number of known age samples is limited, removing this impediment to estimating demographic parameters vital to the conservation of critically endangered species.

Quantification by droplet digital PCR and species identification by metabarcoding of environmental (e)DNA from Blainville's beaked whales, with assisted localization from an acoustic array.

Detection and identification of species, subspecies or stocks of whales, dolphins and porpoises at sea remain challenging, particularly for cryptic or elusive species like beaked whales (Family: Ziphiidae). Here we investigated the potential for using an acoustically assisted sampling design to collect environmental (e)DNA from beaked whales on the U.S. Navy's Atlantic Undersea Test and Evaluation Center (AUTEC) in The Bahamas. During 12 days of August 2019, we conducted 9 small-boat surveys and collected 56 samples of seawater (paired subsamples of 1L each, including controls) using both a spatial collection design in the absence of visual confirmation of whales, and a serial collection design in the proximity of whales at the surface. There were 7 sightings of whales, including 11 Blainville's beaked whales (Mesoplodon densirostris). All whales were located initially with the assistance of information from a bottom-mounted acoustic array available on the AUTEC range. Quantification by droplet digital (dd)PCR from the four spatial design collections showed no samples of eDNA above the threshold of detection and none of these 20 samples yielded amplicons for conventional or next-generation sequencing. Quantification of the 31 samples from four serial collections identified 11 likely positive detections. eDNA barcoding by conventional sequencing and eDNA metabarcoding by next-generation sequencing confirmed species identification for 9 samples from three of the four serial collections. We further resolved five intra-specific variants (i.e., haplotypes), two of which showed an exact match to previously published haplotypes and three that have not been reported previously to the international repository, GenBank. A minimum spanning network of the five eDNA haplotypes, with all other published haplotypes of Blainville's beaked whales, suggested the potential for further resolution of differences between oceanic populations.

Diversity of mitochondrial DNA in 3 species of great whales before and after modern whaling.

The 20th century commercial whaling industry severely reduced populations of great whales throughout the Southern Hemisphere. The effect of this exploitation on genetic diversity and population structure remains largely undescribed. Here, we compare pre- and post-whaling diversity of mitochondrial DNA (mtDNA) control region sequences for 3 great whales in the South Atlantic, such as the blue, humpback, and fin whale. Pre-whaling diversity is described from mtDNA extracted from bones collected near abandoned whaling stations, primarily from the South Atlantic island of South Georgia. These bones are known to represent the first stage of 20th century whaling and thus pre-whaling diversity of these populations. Post-whaling diversity is described from previously published studies reporting large-scale sampling of living whales in the Southern Hemisphere. Despite relatively high levels of surviving genetic diversity in the post-whaling populations, we found evidence of a probable loss of mtDNA lineages in all 3 species. This is evidenced by the detection of a large number of haplotypes found in the pre-whaling samples that are not present in the post-whaling samples. A rarefaction analysis further supports a loss of haplotypes in the South Atlantic humpback and Antarctic blue whale populations. The bones from former whaling stations in the South Atlantic represent a remarkable molecular archive for further investigation of the decline and ongoing recovery in the great whales of the Southern Hemisphere.

A surplus no more? Variation in krill availability impacts reproductive rates of Antarctic baleen whales.

The krill surplus hypothesis of unlimited prey resources available for Antarctic predators due to commercial whaling in the 20th century has remained largely untested since the 1970s. Rapid warming of the Western Antarctic Peninsula (WAP) over the past 50 years has resulted in decreased seasonal ice cover and a reduction of krill. The latter is being exacerbated by a commercial krill fishery in the region. Despite this, humpback whale populations have increased but may be at a threshold for growth based on these human-induced changes. Understanding how climate-mediated variation in prey availability influences humpback whale population dynamics is critical for focused management and conservation actions. Using an 8-year dataset (2013-2020), we show that inter-annual humpback whale pregnancy rates, as determined from skin-blubber biopsy samples (n = 616), are positively correlated with krill availability and fluctuations in ice cover in the previous year. Pregnancy rates showed significant inter-annual variability, between 29% and 86%. Our results indicate that krill availability is in fact limiting and affecting reproductive rates, in contrast to the krill surplus hypothesis. This suggests that this population of humpback whales may be at a threshold for population growth due to prey limitations. As a result, continued warming and increased fishing along the WAP, which continue to reduce krill stocks, will likely impact this humpback whale population and other krill predators in the region. Humpback whales are sentinel species of ecosystem health, and changes in pregnancy rates can provide quantifiable signals of the impact of environmental change at the population level. Our findings must be considered paramount in developing new and more restrictive conservation and management plans for the Antarctic marine ecosystem and minimizing the negative impacts of human activities in the region.

mtDNA heteroplasmy gives rise to a new maternal lineage in North Pacific humpback whales (Megaptera novaeangliae).

Heteroplasmy in the mitochondrial genome offers a rare opportunity to track the evolution of a newly arising maternal lineage in populations of non-model species. Here, we identified a previously unreported mitochondrial DNA haplotype while assembling an integrated database of DNA profiles and photo-identification records from humpback whales in southeastern Alaska (SEAK). The haplotype, referred to as A8, was shared by only 2 individuals, a mature female with her female calf, and differed by only a single base pair from a common haplotype in the North Pacific, referred to as A-. To investigate the origins of the A8 haplotype, we reviewed n = 1,089 electropherograms (including replicate samples) of n = 710 individuals with A- haplotypes from an existing collection. From this review, we found 20 individuals with clear evidence of heteroplasmy for A-/A8 (parental/derived) haplotypes. Of these, 15 were encountered in SEAK, 4 were encountered on the Hawaiian breeding ground (the primary migratory destination for whales in SEAK), and 1 was encountered in the northern Gulf of Alaska. We used genotype exclusion and likelihood to identify one of the heteroplasmic females as the likely mother of the A8 cow and grandmother of the A8 calf, establishing the inheritance and germ-line fixation of the new haplotype from the parental heteroplasmy. The mutation leading to this heteroplasmy and the fixation of the A8 haplotype provide an opportunity to document the population dynamics and regional fidelity of a newly arising maternal lineage in a population recovering from exploitation.

Mitochondrial diversity and inter-specific phylogeny among dolphins of the genus Stenella in the Southwest Atlantic Ocean.

The genus Stenella is comprised of five species occurring in all oceans. Despite its wide distribution, genetic diversity information on these species is still scarce especially in the Southwest Atlantic Ocean. Some features of this genus can enhance opportunities for potential introgressive hybridization, e.g. sympatric distibution along the Brazilian coast, mixed known associations among species, karyotype uniformity and genome permeability. In this study we analyzed three genes of the mitochondrial genome to investigate the genetic diversity and occurrence of genetic mixture among eighty specimens of Stenella. All species exhibited moderate to high levels of genetic diversity (h = 0.833 to h = 1.000 and π = 0.006 to π = 0.015). Specimens of S. longirostris, S. attenuata and S. frontalis were clustered into differentiated haplogroups, in contrast, haplotypes of S. coeruleoalba and S. clymene were clustered together. We detected phylogenetic structure of mixed clades for S. clymene and S. coeruleoalba specimens, in the Southwest Atlantic Ocean, and also between S. frontalis and S. attenuata in the Northeast Atlantic Ocean, and between S. frontalis and S. longirostris in the Northwest Atlantic Ocean. These specimes were morphologically identified as one species but exhibited the maternal lineage of another species, by mitochondrial DNA. Our results demonstrate that ongoing gene flow is occurring among species of the genus Stenella reinforcing that this process could be one of the reasons for the confusing taxonomy and difficulties in elucidating phylogenetic relationships within this group.

Speciation in the deep: genomics and morphology reveal a new species of beaked whale Mesoplodon eueu.

The deep sea has been described as the last major ecological frontier, as much of its biodiversity is yet to be discovered and described. Beaked whales (ziphiids) are among the most visible inhabitants of the deep sea, due to their large size and worldwide distribution, and their taxonomic diversity and much about their natural history remain poorly understood. We combine genomic and morphometric analyses to reveal a new Southern Hemisphere ziphiid species, Ramari's beaked whale, Mesoplodon eueu, whose name is linked to the Indigenous peoples of the lands from which the species holotype and paratypes were recovered. Mitogenome and ddRAD-derived phylogenies demonstrate reciprocally monophyletic divergence between M. eueu and True's beaked whale (M. mirus) from the North Atlantic, with which it was previously subsumed. Morphometric analyses of skulls also distinguish the two species. A time-calibrated mitogenome phylogeny and analysis of two nuclear genomes indicate divergence began circa 2 million years ago (Ma), with geneflow ceasing 0.35-0.55 Ma. This is an example of how deep sea biodiversity can be unravelled through increasing international collaboration and genome sequencing of archival specimens. Our consultation and involvement with Indigenous peoples offers a model for broadening the cultural scope of the scientific naming process.

Corrigendum to: Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales.

[This corrects the article DOI: 10.1093/conphys/cow050.].

Cultural Transmission of Fine-Scale Fidelity to Feeding Sites May Shape Humpback Whale Genetic Diversity in Russian Pacific Waters.

Mitochondrial DNA (mtDNA) differences between humpback whales on different feeding grounds can reflect the cultural transmission of migration destinations over generations, and therefore represent one of the very few cases of gene-culture coevolution identified in the animal kingdom. In Russian Pacific waters, photo-identification (photo-ID) studies have shown minimal interchange between whales feeding off the Commander Islands and those feeding in the Karaginsky Gulf, regions that are separated by only 500 km and have previously been lumped together as a single Russian feeding ground. Here, we assessed whether genetic differentiation exists between these 2 groups of humpback whales. We discovered a strong mtDNA differentiation between the 2 feeding sites (FST = 0.18, ΦST = 0.14, P < 0.001). In contrast, nuclear DNA (nuDNA) polymorphisms, determined at 8 microsatellite loci, did not reveal any differentiation. Comparing our mtDNA results with those from a previous ocean-basin study reinforced the differences between the 2 feeding sites. Humpback whales from the Commanders appeared most similar to those of the western Gulf of Alaska and the Aleutian feeding grounds, whereas Karaginsky differed from all other North Pacific feeding grounds. Comparison to breeding grounds suggests mixed origins for the 2 feeding sites; there are likely connections between Karaginsky and the Philippines and to a lesser extent to Okinawa, Japan, whereas the Commanders are linked to the Mexican breeding grounds. The mtDNA differentiation between the Commander Islands and Karaginsky Gulf suggests a case of gene-culture coevolution, correlated to fidelity to a specific feeding site within a particular feeding ground. From a conservation perspective, our findings emphasize the importance of considering these 2 feeding sites as separate management units.

High pregnancy rates in humpback whales (Megaptera novaeangliae) around the Western Antarctic Peninsula, evidence of a rapidly growing population.

Antarctic humpback whales are recovering from near extirpation from commercial whaling. To understand the dynamics of this recovery and establish a baseline to monitor impacts of a rapidly changing environment, we investigated sex ratios and pregnancy rates of females within the Western Antarctic Peninsula (WAP) feeding population. DNA profiling of 577 tissue samples (2010-2016) identified 239 males and 268 females. Blubber progesterone levels indicated 63.5% of the females biopsied were pregnant. This proportion varied significantly across years, from 36% in 2010 to 86% in 2014. A comparison of samples collected in summer versus fall showed significant increases in the proportion of females present (50% to 59%) and pregnant (59% to 72%), consistent with demographic variation in migratory timing. We also found evidence of annual reproduction among females; 54.5% of females accompanied by a calf were pregnant. These high pregnancy rates are consistent with a population recovering from past exploitation, but appear inconsistent with recent estimates of WAP humpback population growth. Thus, our results will help to better understand population growth potential and set a current baseline from which to determine the impact of climate change and variability on fecundity and reproductive rates.

Erratum: Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales.

[This corrects the article DOI: 10.1093/conphys/cow050.].

Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales.

Understanding reproductive rates of wild animal populations is crucially important for management and conservation. Assessing pregnancy status of free-ranging cetaceans has historically been difficult; however, recent advances in analytical techniques have allowed the diagnosis of pregnancy from small samples of blubber tissue. The primary objectives of this study were as follows: (i) to test the efficacy of blubber progesterone assays as a tool for diagnosing pregnancy in humpback whales (Megaptera novaeangliae); (ii) to estimate the pregnancy rate of humpback whales in Monterey Bay, California; and (iii) to investigate the relationship between stable isotopes and reproductive status of these whales. Progesterone concentrations of female whales fell into two distinct groups, allowing for diagnostic separation of pregnant and non-pregnant individuals. Pregnancy rate varied between years of the study (48.4%% in 2011 and 18.5% in 2012), but fell within the range of other estimates of reproductive success for this population. Stable carbon and nitrogen isotope ratios were examined to investigate the impacts of pregnancy on these values. Neither δ15N nor δ13C varied in a consistent way among animals of different sex or reproductive status. The relationship between δ15N and δ13C was strongly positive for male and non-pregnant female humpbacks; however, no relationship existed for pregnant whales. This difference may be indicative of the effects of pregnancy on δ15N, resulting from tissue synthesis and reduced excretion of nitrogenous waste, as well as on δ13C through increased mobilization of lipid stores to meet the energetic demands of pregnancy. Ultimately, our results support the use of blubber progesterone assays for diagnosing pregnancy in humpback whales and indicate that, when paired with other approaches (e.g. stable isotope analysis), pregnancy status can be an informative tool for addressing questions about animal physiology, ecology and population biology. This information will provide for more effective management and conservation efforts in a rapidly changing world.

What influences the worldwide genetic structure of sperm whales (Physeter macrocephalus)?

The interplay of natural selection and genetic drift, influenced by geographic isolation, mating systems and population size, determines patterns of genetic diversity within species. The sperm whale provides an interesting example of a long-lived species with few geographic barriers to dispersal. Worldwide mtDNA diversity is relatively low, but highly structured among geographic regions and social groups, attributed to female philopatry. However, it is unclear whether this female philopatry is due to geographic regions or social groups, or how this might vary on a worldwide scale. To answer these questions, we combined mtDNA information for 1091 previously published samples with 542 newly obtained DNA profiles (394-bp mtDNA, sex, 13 microsatellites) including the previously unsampled Indian Ocean, and social group information for 541 individuals. We found low mtDNA diversity (π = 0.430%) reflecting an expansion event <80 000 years bp, but strong differentiation by ocean, among regions within some oceans, and among social groups. In comparison, microsatellite differentiation was low at all levels, presumably due to male-mediated gene flow. A hierarchical amova showed that regions were important for explaining mtDNA variance in the Indian Ocean, but not Pacific, with social group sampling in the Atlantic too limited to include in analyses. Social groups were important in partitioning mtDNA and microsatellite variance within both oceans. Therefore, both geographic philopatry and social philopatry influence genetic structure in the sperm whale, but their relative importance differs by sex and ocean, reflecting breeding behaviour, geographic features and perhaps a more recent origin of sperm whales in the Pacific. By investigating the interplay of evolutionary forces operating at different temporal and geographic scales, we show that sperm whales are perhaps a unique example of a worldwide population expansion followed by rapid assortment due to female social organization.

Global diversity and oceanic divergence of humpback whales (Megaptera novaeangliae).

Humpback whales (Megaptera novaeangliae) annually undertake the longest migrations between seasonal feeding and breeding grounds of any mammal. Despite this dispersal potential, discontinuous seasonal distributions and migratory patterns suggest that humpbacks form discrete regional populations within each ocean. To better understand the worldwide population history of humpbacks, and the interplay of this species with the oceanic environment through geological time, we assembled mitochondrial DNA control region sequences representing approximately 2700 individuals (465 bp, 219 haplotypes) and eight nuclear intronic sequences representing approximately 70 individuals (3700 bp, 140 alleles) from the North Pacific, North Atlantic and Southern Hemisphere. Bayesian divergence time reconstructions date the origin of humpback mtDNA lineages to the Pleistocene (880 ka, 95% posterior intervals 550-1320 ka) and estimate radiation of current Northern Hemisphere lineages between 50 and 200 ka, indicating colonization of the northern oceans prior to the Last Glacial Maximum. Coalescent analyses reveal restricted gene flow between ocean basins, with long-term migration rates (individual migrants per generation) of less than 3.3 for mtDNA and less than 2 for nuclear genomic DNA. Genetic evidence suggests that humpbacks in the North Pacific, North Atlantic and Southern Hemisphere are on independent evolutionary trajectories, supporting taxonomic revision of M. novaeangliae to three subspecies.

Low diversity in the mitogenome of sperm whales revealed by next-generation sequencing.

Large population sizes and global distributions generally associate with high mitochondrial DNA control region (CR) diversity. The sperm whale (Physeter macrocephalus) is an exception, showing low CR diversity relative to other cetaceans; however, diversity levels throughout the remainder of the sperm whale mitogenome are unknown. We sequenced 20 mitogenomes from 17 sperm whales representative of worldwide diversity using Next Generation Sequencing (NGS) technologies (Illumina GAIIx, Roche 454 GS Junior). Resequencing of three individuals with both NGS platforms and partial Sanger sequencing showed low discrepancy rates (454-Illumina: 0.0071%; Sanger-Illumina: 0.0034%; and Sanger-454: 0.0023%) confirming suitability of both NGS platforms for investigating low mitogenomic diversity. Using the 17 sperm whale mitogenomes in a phylogenetic reconstruction with 41 other species, including 11 new dolphin mitogenomes, we tested two hypotheses for the low CR diversity. First, the hypothesis that CR-specific constraints have reduced diversity solely in the CR was rejected as diversity was low throughout the mitogenome, not just in the CR (overall diversity π = 0.096%; protein-coding 3rd codon = 0.22%; CR = 0.35%), and CR phylogenetic signal was congruent with protein-coding regions. Second, the hypothesis that slow substitution rates reduced diversity throughout the sperm whale mitogenome was rejected as sperm whales had significantly higher rates of CR evolution and no evidence of slow coding region evolution relative to other cetaceans. The estimated time to most recent common ancestor for sperm whale mitogenomes was 72,800 to 137,400 years ago (95% highest probability density interval), consistent with previous hypotheses of a bottleneck or selective sweep as likely causes of low mitogenome diversity.

Paternity assignment and demographic closure in the New Zealand southern right whale.

The identification and characterization of reproductively isolated subpopulations or 'stocks' are essential for effective conservation and management decisions. This can be difficult in vagile marine species like marine mammals. We used paternity assignment and 'gametic recapture' to examine the reproductive autonomy of southern right whales (Eubalaena australis) on their New Zealand (NZ) calving grounds. We derived DNA profiles for 34 mother-calf pairs from skin biopsy samples, using sex-specific markers, 13 microsatellite loci and mtDNA haplotypes. We constructed DNA profiles for 314 adult males, representing 30% of the census male abundance of the NZ stock, previously estimated from genotypic mark-recapture modelling to be 1085 (95% CL 855, 1416). Under the hypothesis of demographic closure and the assumption of equal reproductive success among males, we predict: (i) the proportion of paternities assigned will reflect the proportion of the male population sampled and (ii) the gametic mark-recapture (GMR) estimate of male abundance will be equivalent to the census male estimate for the NZ stock. Consistent with these predictions, we found that the proportion of assigned paternities equalled the proportion of the census male population size sampled. Using the sample of males as the initial capture, and paternity assignment as the recapture, the GMR estimate of male abundance was 1001 (95% CL 542, 1469), similar to the male census estimate. These findings suggest that right whales returning to the NZ calving ground are reproductively autonomous on a generational timescale, as well as isolated by maternal fidelity on an evolutionary timescale, from others in the Indo-Pacific region.

Genetic evidence of illegal trade in protected whales links Japan with the US and South Korea.

We report on genetic identification of 'whale meat' purchased in sushi restaurants in Los Angeles, CA (USA) in October 2009 and in Seoul, South Korea in June and September 2009. Phylogenetic analyses of mtDNA cytochrome b sequences confirmed that the products included three species of whale currently killed in the controversial scientific whaling programme of Japan, but which are protected from international trade: the fin, sei and Antarctic minke. The DNA profile of the fin whale sold in Seoul established a match to products purchased previously in Japan in September 2007, confirming unauthorized trade between these two countries. Following species identification, these products were handed over to the appropriate national or local authorities for further investigation. The illegal trade of products from protected species of whales, presumably taken under a national permit for scientific research, is a timely reminder of the need for independent, transparent and robust monitoring of any future whaling.

Applied conservation genetics and the need for quality control and reporting of genetic data used in fisheries and wildlife management.

Genetic data are often critical for defining populations for management purposes (e.g., identifying geographic boundaries or diagnostic characters for genetically discrete subunits) but can be called into question by both scientific and legal review. This can result in reversed or delayed implementation of management actions. We discuss methods for data quality control and quality analysis and describe examples of steps applied to 2 of the most common types of genetic data, mitochondrial DNA sequences, and microsatellite genotypes. These steps can serve both as guides to conservation geneticists and as an initial protocol for managers to determine whether genetic data will hold up against legal and scientific challenges. In addition, we suggest types of data and quality measures that should be reported as supplementary materials to published reports. These supplementary data serve to reduce the occurrence of legal and conservation controversies and improve reproducibility over time in population genetics studies where genetic monitoring is likely to play an increasing role.

Phylogeny of the beaked whale genus Mesoplodon (Ziphiidae: Cetacea) revealed by nuclear introns: implications for the evolution of male tusks.

With 14 species currently recognized, the beaked whale genus Mesoplodon (family Ziphiidae) is the most speciose in the order Cetacea. Beaked whales are widely distributed but are rarely seen at sea due to their oceanic distribution, deep-diving capacity, and apparent low abundance. Morphological differentiation among Mesoplodon species is relatively limited, with the exception of tooth form in adult males. Based on scarring patterns, males appear to use their tusk-like teeth as weapons in aggressive encounters with other males. Females are effectively toothless. We used sequences from seven nuclear introns (3348 base pairs) to construct a robust and highly resolved phylogeny, which was then used as a framework to test predictions from four hypotheses seeking to explain patterns of Mesoplodon tusk morphology and/or the processes that have driven the diversification of this genus: (1) linear progression of tusk form; (2) allopatric speciation through isolation in adjacent deep-sea canyons; (3) sympatric speciation through sexual selection on tusks; and (4) selection for species-recognition cues. Maximum-likelihood and Bayesian reconstructions confirmed the monophyly of the genus and revealed that what were considered ancestral and derived tusk forms have in fact arisen independently on several occasions, contrary to predictions from the linear-progression hypothesis. Further, none of the three well-supported species clades was confined to a single ocean basin, as might have been expected from the deep-sea canyon-isolation or sexual-selection hypotheses, and some species with similar tusks have overlapping distributions, contrary to predictions from the species-recognition hypothesis. However, the divergent tusk forms and sympatric distributions of three of the four sister-species pairs identified suggest that sexual selection on male tusks has likely played an important role in this unique radiation, although other forces are clearly also involved. To our knowledge, this is the first time that sexual selection has been explicitly implicated in the radiation of a mammalian group outside terrestrial ungulates.