Influences of past climatic changes on historical population structure and demography of a cosmopolitan marine predator, the common dolphin (genus Delphinus)

Chromosome-level genome provides insights into environmental adaptability and innate immunity in the common dolphin (delphinus delphis)

The common dolphin (Delphinus delphis) is widely distributed worldwide and well adapted to various habitats. Animal genomes store clues about their pasts, and can reveal the genes underlying their evolutionary success. Here, we report the first high-quality chromosome-level genome of D. delphis. The assembled genome size was 2.56 Gb with a contig N50 of 63.85 Mb. Phylogenetically, D. delphis was close to Tursiops truncatus and T. aduncus. The genome of D. delphis exhibited 428 expanded and 1,885 contracted gene families, and 120 genes were identified as positively selected. The expansion of the HSP70 gene family suggested that D. delphis has a powerful system for buffering stress, which might be associated with its broad adaptability, longevity, and detoxification capacity. The expanded IFN-α and IFN-ω gene families, as well as the positively selected genes encoding tripartite motif-containing protein 25, peptidyl-prolyl cis-trans isomerase NIMA-interacting 1, and p38 MAP kinase, were all involved in pathways for antiviral, anti-inflammatory, and antineoplastic mechanisms. The genome data also revealed dramatic fluctuations in the effective population size during the Pleistocene. Overall, the high-quality genome assembly and annotation represent significant molecular resources for ecological and evolutionary studies of Delphinus and help support their sustainable treatment and conservation.

Ecomorphology of toothed whales (Cetacea, Odontoceti) as revealed by 3D skull geometry

Extant odontocetes (toothed whales) exhibit differences in body size and brain mass, biosonar mode, feeding strategies, and diving and habitat adaptations. Strong selective pressures associated with these factors have likely contributed to the morphological diversification of their skull. Here, we used 3D landmark geometric morphometric data from the skulls of 60 out of ~ 72 extant odontocete species and a well-supported phylogenetic tree to test whether size and shape variation are associated with ecological adaptations at an interspecific scale. Odontocete skull morphology exhibited a significant phylogenetic signal, with skull size showing stronger signal than shape. After accounting for phylogeny, significant associations were detected between skull size and biosonar mode, body length, brain and body mass, maximum and minimum prey size, and maximum peak frequency. Brain mass was also strongly correlated with skull shape together with surface temperature and average and minimum prey size. When asymmetric and symmetric components of shape were analysed separately, a significant correlation was detected between sea surface temperature and both symmetric and asymmetric components of skull shape, and between diving ecology and the asymmetric component. Skull shape variation of odontocetes was strongly influenced by evolutionary allometry but most of the associations with ecological variables were not supported after phylogenetic correction. This suggests that ecomorphological feeding adaptations vary more between, rather than within, odontocete families, and functional anatomical patterns across odontocete clades are canalised by size constraints.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

Seascape Genetics of the Atlantic Spotted Dolphin (Stenella frontalis) Based on Mitochondrial DNA

The Atlantic spotted dolphin (Stenella frontalis) is endemic to tropical, subtropical, and warm temperate waters of the Atlantic Ocean. Throughout its distribution, both geographic distance and environmental variation may contribute to population structure of the species. In this study, we follow a seascape genetics approach to investigate population differentiation of Atlantic spotted dolphins based on a large worldwide dataset and the relationship with marine environmental variables. The results revealed that the Atlantic spotted dolphin exhibits population genetic structure across its distribution based on mitochondrial DNA control region (mtDNA-CR) data. Analyses based on the contemporary landscape suggested, at both the individual and population level, that the population genetic structure is consistent with the isolation-by-distance model. However, because geography and environmental matrices were correlated, and because in some, but not all analyses, we found a significant effect for the environment, we cannot rule out the addition contribution of environmental factors in structuring genetic variation. Future analyses based on nuclear data are needed to evaluate whether local processes, such as social structure and some level of philopatry within populations, may be contributing to the associations among genetic structure, geographic, and environmental distance.

Mitochondrial genomics reveals the evolutionary history of the porpoises (Phocoenidae) across the speciation continuum

Historical variation in food resources is expected to be a major driver of cetacean evolution, especially for the smallest species like porpoises. Despite major conservation issues among porpoise species (e.g., vaquita and finless), their evolutionary history remains understudied. Here, we reconstructed their evolutionary history across the speciation continuum. Phylogenetic analyses of 63 mitochondrial genomes suggest that porpoises radiated during the deep environmental changes of the Pliocene. However, all intra-specific subdivisions were shaped during the Quaternary glaciations. We observed analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. This suggests that similar mechanisms are driving species diversification in northern (harbor and Dall's) and southern species (spectacled and Burmeister's). In contrast to previous studies, spectacled and Burmeister's porpoises shared a more recent common ancestor than with the vaquita that diverged from southern species during the Pliocene. The low genetic diversity observed in the vaquita carried signatures of a very low population size since the last 5,000 years. Cryptic lineages within Dall's, spectacled and Pacific harbor porpoises suggest a richer evolutionary history than previously suspected. These results provide a new perspective on the mechanisms driving diversification in porpoises and an evolutionary framework for their conservation.

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.

Population genetic structure of the short-beaked common dolphin from the Black Sea and the Turkish Straits System

Our study aims to assess the population connectivity, evolutionary history, and conservation status of the short-beaked common dolphin in the Black Sea and Turkish Straits System (TSS). We also include DNA sequences from the Atlantic Ocean and the Mediterranean Sea to provide a regional perspective to our localized study. Analysis of 366 base pairs of mitochondrial DNA D-loop fragments from 37 samples collected from short-beaked common dolphins in the Black Sea, TSS, and Aegean Sea revealed 13 haplotypes, eight of which have not been previously reported. While analysis of samples archived on GenBank revealed 89 different haplotypes across the region. The haplotype network contains two main peripheral groups that include individuals from all locations. Haplotypes from the Atlantic Ocean are scattered across the network and no obvious population separation was detected. Some shared haplotypes potentially indicate multi-directional colonization events of the Mediterranean Sea from the eastern Atlantic Ocean. Moreover, some less widely distributed haplotypes suggest some level of more recent genetic connectivity through the Strait of Gibraltar and the TSS and point out the importance of these straits in the dispersal of short-beaked common dolphins. The haplotype and nucleotide diversity values were lower in the Black Sea, TSS, and western Mediterranean Sea when compared to the Atlantic Ocean, supporting the expansion of Atlantic populations into the Mediterranean and the Black Seas. Differentiation was observed between the Atlantic Ocean, and the Mediterranean Sea, TSS and the Black Sea based on Фst but not between Mediterranean and the Black Seas. For common dolphins, which have high dispersal potential, the protection of open seas and narrow seaways to enhance connectivity may be crucial.

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.

High genetic structure and low mitochondrial diversity in bottlenose dolphins of the Archipelago of Bocas del Toro, Panama: A population at risk?

The current conservation status of the bottlenose dolphin (Tursiops truncatus) under the IUCN is ‘least concern’. However, in the Caribbean, small and localized populations of the ‘inshore form’ may be at higher risk of extinction than the ‘worldwide distributed form’ due to a combination of factors including small population size, high site fidelity, genetic isolation, and range overlap with human activities. Here, we study the population genetic structure of bottlenose dolphins from the Archipelago of Bocas del Toro in Panama. This is a small population characterized by high site fidelity and is currently heavily-impacted by the local dolphin-watching industry. We collected skin tissue samples from 25 dolphins to study the genetic diversity and structure of this population. We amplified a portion of the mitochondrial Control Region (mtDNA-CR) and nine microsatellite loci. The mtDNA-CR analyses revealed that dolphins in Bocas del Toro belong to the ‘inshore form’, grouped with the Bahamas-Colombia-Cuba-Mexico population unit. They also possess a unique haplotype new for the Caribbean. The microsatellite data indicated that the Bocas del Toro dolphin population is highly structured, likely due to restricted movement patterns. Previous abundance estimates obtained with mark-recapture methods reported a small population of 80 dolphins (95% CI = 72–87), which is similar to the contemporary effective population size estimated in this study (Ne = 73 individuals; CI = 18.0 - ∞; 0.05). The combination of small population size, high degree of genetic isolation, and intense daily interactions with dolphin-watching boats puts the Bocas del Toro dolphin to at high risk of extinction. Despite national guidelines to regulate the dolphin-watching industry in Bocas del Toro and ongoing educational programs for tour operators, only in 2012 seven animals have died due to boat collisions. Our results suggest that the conservation status of bottlenose dolphins in Bocas del Toro should be elevated to ‘endangered’ at the national level, as a precautionary measure while population and viability estimates are conducted.

Mixing of porpoise ecotypes in southwestern UK waters revealed by genetic profiling

Contact zones between ecotypes are windows for understanding how species may react to climate changes. Here, we analysed the fine-scale genetic and morphological variation in harbour porpoises (Phocoena phocoena) around the UK by genotyping 591 stranded animals at nine microsatellite loci. The data were integrated with a prior study to map at high resolution the contact zone between two previously identified ecotypes meeting in the northern Bay of Biscay. Clustering and spatial analyses revealed that UK porpoises are derived from two genetic pools with porpoises from the southwestern UK being genetically differentiated, and having larger body sizes compared to those of other UK areas. Southwestern UK porpoises showed admixed ancestry between southern and northern ecotypes with a contact zone extending from the northern Bay of Biscay to the Celtic Sea and Channel. Around the UK, ancestry blends from one genetic group to the other along a southwest--northeast axis, correlating with body size variation, consistent with previously reported morphological differences between the two ecotypes. We also detected isolation by distance among juveniles but not in adults, suggesting that stranded juveniles display reduced intergenerational dispersal. The fine-scale structure of this admixture zone raises the question of how it will respond to future climate change and provides a reference point for further study.

Post-glacial habitat release and incipient speciation in the genus Delphinus

The role of ecological and changing environmental factors in the radiation of species diversity is a fundamental question in evolutionary biology. Of particular interest is the potential for these factors to determine the boundary between what we would consider differentiation among populations and incipient speciation. Dolphins in the genus Delphinus provide a useful test case, exhibiting morphological variation in beak length, coloration and body size across their wide geographic distribution, and in particular among coastal and more pelagic habitats. Two species have been proposed, D. delphis and D. capensis, but morphologically similar allopatric populations are not monophyletic, indicating that the mostly coastal 'long-beaked' D. capensis form is not a single globally distributed species. However, the sympatric populations in the Eastern North Pacific currently designated as these two species are both morphologically and genetically differentiated. Here we use microsatellite DNA and mitochondrial DNA markers to investigate the evolutionary mechanisms that led to this incipient speciation event. We used coalescent and assignment methods to investigate the timing and extent of reproductive isolation. Our data indicate that although there is some level of on-going gene flow, the putative species found in the Eastern North Pacific are reciprocally monophyletic. The timing of isolation appears to be associated with regional changes in paleoceanographic conditions within the Holocene timeframe.

Molecular and Morphological Differentiation of Common Dolphins (Delphinus sp.) in the Southwestern Atlantic: Testing the Two Species Hypothesis in Sympatry

The taxonomy of common dolphins (Delphinus sp.) has always been controversial, with over twenty described species since the original description of the type species of the genus (Delphinus delphis Linnaeus, 1758). Two species and four subspecies are currently accepted, but recent molecular data have challenged this view. In this study we investigated the molecular taxonomy of common dolphins through analyses of cytochrome b sequences of 297 individuals from most of their distribution. We included 37 novel sequences from the Southwestern Atlantic Ocean, a region where the short- and long-beaked morphotypes occur in sympatry, but which had not been well sampled before. Skulls of individuals from the Southwestern Atlantic were measured to test the validity of the rostral index as a diagnostic character and confirmed the presence of the two morphotypes in our genetic sample. Our genetic results show that all common dolphins in the Atlantic Ocean belong to a single species, Delphinus delphis. According to genetic data, the species Delphinus capensis is invalid. Long-beaked common dolphins from the Northeastern Pacific Ocean may constitute a different species. Our conclusions prompt the need for revision of currently accepted common dolphin species and subspecies and of Delphinus delphis distribution.

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.

Understanding the mechanisms of antitropical divergence in the seabird White-faced Storm-petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White-faced Storm-petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial ΦST = 0.89 (P < 0.01) and global nuclear ΦST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long-term antitropical divergence in isolation. Coalescent-based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000-1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000-600 000), both probably associated with the profound palaeo-oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long-distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.

Phylogenomics of the killer whale indicates ecotype divergence in sympatry

For many highly mobile species, the marine environment presents few obvious barriers to gene flow. Even so, there is considerable diversity within and among species, referred to by some as the 'marine speciation paradox'. The recent and diverse radiation of delphinid cetaceans (dolphins) represents a good example of this. Delphinids are capable of extensive dispersion and yet many show fine-scale genetic differentiation among populations. Proposed mechanisms include the division and isolation of populations based on habitat dependence and resource specializations, and habitat release or changing dispersal corridors during glacial cycles. Here we use a phylogenomic approach to investigate the origin of differentiated sympatric populations of killer whales (Orcinus orca). Killer whales show strong specialization on prey choice in populations of stable matrifocal social groups (ecotypes), associated with genetic and phenotypic differentiation. Our data suggest evolution in sympatry among populations of resource specialists.

Killer whale nuclear genome and mtDNA reveal widespread population bottleneck during the last glacial maximum

Ecosystem function and resilience is determined by the interactions and independent contributions of individual species. Apex predators play a disproportionately determinant role through their influence and dependence on the dynamics of prey species. Their demographic fluctuations are thus likely to reflect changes in their respective ecological communities and habitat. Here, we investigate the historical population dynamics of the killer whale based on draft nuclear genome data for the Northern Hemisphere and mtDNA data worldwide. We infer a relatively stable population size throughout most of the Pleistocene, followed by an order of magnitude decline and bottleneck during the Weichselian glacial period. Global mtDNA data indicate that while most populations declined, at least one population retained diversity in a stable, productive ecosystem off southern Africa. We conclude that environmental changes during the last glacial period promoted the decline of a top ocean predator, that these events contributed to the pattern of diversity among extant populations, and that the relatively high diversity of a population currently in productive, stable habitat off South Africa suggests a role for ocean productivity in the widespread decline.

Hybrid speciation in a marine mammal: the clymene dolphin (Stenella clymene)

Natural hybridization may result in the exchange of genetic material between divergent lineages and even the formation of new taxa. Many of the Neo-Darwinian architects argued that, particularly for animal clades, natural hybridization was maladaptive. Recent evidence, however, has falsified this hypothesis, instead indicating that this process may lead to increased biodiversity through the formation of new species. Although such cases of hybrid speciation have been described in plants, fish and insects, they are considered exceptionally rare in mammals. Here we present evidence for a marine mammal, Stenella clymene, arising through natural hybridization. We found phylogenetic discordance between mitochondrial and nuclear markers, which, coupled with a pattern of transgressive segregation seen in the morphometric variation of some characters, support a case of hybrid speciation. S. clymene is currently genetically differentiated from its putative parental species, Stenella coerueloalba and Stenella longisrostris, although low levels of introgressive hybridization may be occurring. Although non-reticulate forms of evolution, such as incomplete lineage sorting, could explain our genetic results, we consider that the genetic and morphological evidence taken together argue more convincingly towards a case of hybrid speciation. We anticipate that our study will bring attention to this important aspect of reticulate evolution in non-model mammal species. The study of speciation through hybridization is an excellent opportunity to understand the mechanisms leading to speciation in the context of gene flow.

De novo assembly of the Indo-Pacific humpback dolphin leucocyte transcriptome to identify putative genes involved in the aquatic adaptation and immune response

Background

The Indo-Pacific humpback dolphin (Sousa chinensis), a marine mammal species inhabited in the waters of Southeast Asia, South Africa and Australia, has attracted much attention because of the dramatic decline in population size in the past decades, which raises the concern of extinction. So far, this species is poorly characterized at molecular level due to little sequence information available in public databases. Recent advances in large-scale RNA sequencing provide an efficient approach to generate abundant sequences for functional genomic analyses in the species with un-sequenced genomes.

Principal Findings

We performed a de novo assembly of the Indo-Pacific humpback dolphin leucocyte transcriptome by Illumina sequencing. 108,751 high quality sequences from 47,840,388 paired-end reads were generated, and 48,868 and 46,587 unigenes were functionally annotated by BLAST search against the NCBI non-redundant and Swiss-Prot protein databases (E-value<10⁻⁵), respectively. In total, 16,467 unigenes were clustered into 25 functional categories by searching against the COG database, and BLAST2GO search assigned 37,976 unigenes to 61 GO terms. In addition, 36,345 unigenes were grouped into 258 KEGG pathways. We also identified 9,906 simple sequence repeats and 3,681 putative single nucleotide polymorphisms as potential molecular markers in our assembled sequences. A large number of unigenes were predicted to be involved in immune response, and many genes were predicted to be relevant to adaptive evolution and cetacean-specific traits.

Conclusion

This study represented the first transcriptome analysis of the Indo-Pacific humpback dolphin, an endangered species. The de novo transcriptome analysis of the unique transcripts will provide valuable sequence information for discovery of new genes, characterization of gene expression, investigation of various pathways and adaptive evolution, as well as identification of genetic markers.