The critically endangered vaquita is not doomed to extinction by inbreeding depression.

In cases of severe wildlife population decline, a key question is whether recovery efforts will be impeded by genetic factors, such as inbreeding depression. Decades of excess mortality from gillnet fishing have driven Mexico's vaquita porpoise (Phocoena sinus) to ~10 remaining individuals. We analyzed whole-genome sequences from 20 vaquitas and integrated genomic and demographic information into stochastic, individual-based simulations to quantify the species' recovery potential. Our analysis suggests that the vaquita's historical rarity has resulted in a low burden of segregating deleterious variation, reducing the risk of inbreeding depression. Similarly, genome-informed simulations suggest that the vaquita can recover if bycatch mortality is immediately halted. This study provides hope for vaquitas and other naturally rare endangered species and highlights the utility of genomics in predicting extinction risk.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples.

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (FST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Reference genome and demographic history of the most endangered marine mammal, the vaquita.

The vaquita is the most critically endangered marine mammal, with fewer than 19 remaining in the wild. First described in 1958, the vaquita has been in rapid decline for more than 20 years resulting from inadvertent deaths due to the increasing use of large-mesh gillnets. To understand the evolutionary and demographic history of the vaquita, we used combined long-read sequencing and long-range scaffolding methods with long- and short-read RNA sequencing to generate a near error-free annotated reference genome assembly from cell lines derived from a female individual. The genome assembly consists of 99.92% of the assembled sequence contained in 21 nearly gapless chromosome-length autosome scaffolds and the X-chromosome scaffold, with a scaffold N50 of 115 Mb. Genome-wide heterozygosity is the lowest (0.01%) of any mammalian species analysed to date, but heterozygosity is evenly distributed across the chromosomes, consistent with long-term small population size at genetic equilibrium, rather than low diversity resulting from a recent population bottleneck or inbreeding. Historical demography of the vaquita indicates long-term population stability at less than 5,000 (Ne) for over 200,000 years. Together, these analyses indicate that the vaquita genome has had ample opportunity to purge highly deleterious alleles and potentially maintain diversity necessary for population health.

Genetic homogeneity in the face of morphological heterogeneity in the harbor porpoise from the Black Sea and adjacent waters (Phocoena phocoena relicta).

Absence of genetic differentiation is usually taken as an evidence of panmixia, but can also reflect other situations, including even nearly complete demographic independence among large-sized populations. Deciphering which situation applies has major practical implications (e.g., in conservation biology). The endangered harbor porpoises in the Black Sea illustrates this point well. While morphological heterogeneity suggested that population differentiation may exist between individuals from the Black and Azov seas, no genetic study provided conclusive evidence or covered the entire subspecies range. Here, we assessed the genetic structure at ten microsatellite loci and a 3904 base-pairs mitochondrial fragment in 144 porpoises across the subspecies range (i.e., Aegean, Marmara, Black, and Azov seas). Analyses of the genetic structure, including FST, Bayesian clustering, and multivariate analyses revealed a nearly complete genetic homogeneity. Power analyses rejected the possibility of underpowered analyses (power to detect FST ≥ 0.008 at microsatellite loci). Simulations under various demographic models, evaluating the evolution of FST, showed that a time-lag effect between demographic and genetic subdivision is also unlikely. With a realistic effective population size of 1000 individuals, the expected "gray zone" would be at most 20 generations under moderate levels of gene flow (≤10 migrants per generation). After excluding alternative hypotheses, panmixia remains the most likely hypothesis explaining the genetic homogeneity in the Black Sea porpoises. Morphological heterogeneity may thus reflect other processes than population subdivision (e.g., plasticity, selection). This study illustrates how combining empirical and theoretical approaches can contribute to understanding patterns of weak population structure in highly mobile marine species.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena).

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Hepatic DNA damage in harbour porpoises (Phocoena phocoena) stranded along the English and Welsh coastlines.

One level at which persistent organic pollutants (POPs) and polycyclic aromatic hydrocarbons PAHs) can exert damage is by causing DNA strand-breaks or nucleotide base modifications, which, if unrepaired, can lead to embryonic mutations, abnormal development and cancer. In marine ecosystems, genotoxicity is expected to be particularly strong in long-lived apex predators due to pollutant bioaccumulation. We conducted 32 P-postlabeling analyses optimized for the detection and quantification of aromatic/hydrophobic DNA adducts in the livers of 40 sexually-mature North Atlantic harbour porpoises (Phocoena phocoena) stranded along the English and Welsh coastlines. We examined hepatic tissue to search for inflammatory and preneoplastic lesions and examine their association with adduct levels. Adducts were found in all porpoises (mean: 17.56 ± 11.95 per 108 nucleotides), and were higher than levels reported for marine vertebrates from polluted sites. The pollutants causing the induced DNA adducts could not be further characterized. Hepatic DNA damage did not correlate with levels of blubber POP concentrations (including total polychlorinated biphenyl [PCBs], dichlorodiphenyltrichloroethane [DDT] and dieldrin); PAH concentrations were not available for the present study. However, DNA damage predicted occurrence of inflammatory and preneoplastic lesions. Further, our data showed a reduction in hepatic DNA adduct levels with age in the 40 animals examined while POP concentrations, particularly PCBs, increased with age. Using a different dataset of 145 mature male harbour porpoises confirmed that higher contaminant levels (total PCBs, DDT and dieldrin) are found in older animals. The reduction in hepatic DNA adduct levels in older animals was in accordance with other studies which show that suppression of hepatic CYP1A enzyme activity at high PCB concentrations might impact on CYP1A-mediated DNA adduct formation of PAHs which are ubiquitous environmental pollutants and readily metabolized by CYP1A to species binding to DNA. In summary, our study shows that pollutant-induced DNA damage is prevalent in harbour porpoises from UK waters and may lead to detectable sub-lethal hepatic damage. Environ. Mol. Mutagen. 59:613-624, 2018. © 2018 The Authors Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Candidate gene molecular markers as tools for analyzing genetic susceptibility to morbillivirus infection in stranded Cetaceans.

Morbilliviruses, such as Cetacean morbillivirus (CeMV) or Phocine distemper virus (PDV), represent a growing threat for marine mammals on both hemispheres. Because free-ranging animal populations strongly rely on natural resistance mechanisms, innate immunity-related genes and virus cell entry receptor genes may represent key factors involved in susceptibility to CeMV in Cetaceans. Using the next generation sequencing technology, we have sequenced 11 candidate genes in two model species, Stenella coeruleoalba and Phocoena phocoena. Suitable single nucleotide polymorphism markers of potential functional importance, located in genes coding for basigin (BSG, CD147), the signaling lymphocyte activating molecule (SLAMF1), the poliovirus-related receptor-4 (NECTIN4, PVRL4), toll-like receptors 3, 7, 8 (TLR3, TLR7, TLR8), natural resistance-associated macrophage protein (SLC11A1) and natural cytotoxicity triggering receptor 1 (NCR1), were identified in each model species, along with MHC-DQB haplotypes unique for each species. This set of molecular markers represents a potentially useful tool for studying host genetic variation and susceptibility to morbillivirus infection in Cetaceans as well as for studying functionally important genetic diversity of selected Cetacean populations.

Is there a distinct harbor porpoise subpopulation in the Marmara Sea?

Genetic population structure of geographically isolated endangered Black Sea harbor porpoise (Phocoena phocoena relicta) is little known in Turkish waters, especially in the Turkish Straits System (TSS- Marmara Sea, Bosphorus and Dardanelles), which connects the Black Sea and the Aegean Sea. Mitochondrial DNA sequences of 70 new individuals sampled in the Turkish Black Sea, TSS and Aegean Sea, revealed five new haplotypes from the Black Sea. The findings support the idea that harbor porpoises from the Black Sea dispersed into the Aegean through the TSS. Considering signatures of population expansion, all subpopulations showed a signature of population expansion. The network data and the Фst calculations indicated that the Marmara Sea subpopulation was significantly differentiated from all of the other subpopulations, and supports the notion of its isolated. The finding of a potential management unit (MU) within an already heavily impacted subpopulation as a whole suggests that the individuals of P. p. relicta inhabiting the Marmara Sea require a very rigorous conservation strategy to ensure the survival of this subpopulation, represented by its unique haplotype.

Harbour Porpoises, Phocoena phocoena, in the Mediterranean Sea and Adjacent Regions: Biogeographic Relicts of the Last Glacial Period.

The harbour porpoise, Phocoena phocoena, is one of the best studied cetacean species owing to its common distribution along the coastal waters of the Northern Hemisphere. In European waters, strandings are common and bycatch mortalities in commercial fisheries reach alarming numbers. Lethal interactions resulting from human activities together with ongoing environmental changes raise serious concerns about population viability throughout the species' range. These concerns foster the need to fill critical gaps in knowledge of harbour porpoise biology, including population structure, feeding ecology, habitat preference and evolutionary history, that are critical information for planning effective management and conservation efforts. While the species is distributed fairly continuously in the North Atlantic Ocean, it becomes fragmented in the south-eastern part with isolated populations occurring along the Atlantic coasts of the Iberian Peninsula, Northwest Africa and the Black Sea. The latter population is separated from Atlantic populations by the Mediterranean Sea, where the species is almost entirely absent. Understanding the evolutionary history of these populations occurring in marginal habitats holds the potential to reveal fundamental aspects of the species' biology such as the factors determining its distribution, ecological niche, and how past and recent environmental variation have shaped the current population structure. This information can be critical for understanding the future evolution of the species in consideration of ongoing environmental changes. This chapter summarizes the recent advances in our knowledge regarding the populations bordering the Mediterranean Sea with a special emphasis on their ecological and evolutionary history, which has recently been reconstructed from genetic analyses.

Spatially Explicit Analysis of Genome-Wide SNPs Detects Subtle Population Structure in a Mobile Marine Mammal, the Harbor Porpoise.

The population structure of the highly mobile marine mammal, the harbor porpoise (Phocoena phocoena), in the Atlantic shelf waters follows a pattern of significant isolation-by-distance. The population structure of harbor porpoises from the Baltic Sea, which is connected with the North Sea through a series of basins separated by shallow underwater ridges, however, is more complex. Here, we investigated the population differentiation of harbor porpoises in European Seas with a special focus on the Baltic Sea and adjacent waters, using a population genomics approach. We used 2872 single nucleotide polymorphisms (SNPs), derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), as well as 13 microsatellite loci and mitochondrial haplotypes for the same set of individuals. Spatial principal components analysis (sPCA), and Bayesian clustering on a subset of SNPs suggest three main groupings at the level of all studied regions: the Black Sea, the North Atlantic, and the Baltic Sea. Furthermore, we observed a distinct separation of the North Sea harbor porpoises from the Baltic Sea populations, and identified splits between porpoise populations within the Baltic Sea. We observed a notable distinction between the Belt Sea and the Inner Baltic Sea sub-regions. Improved delineation of harbor porpoise population assignments for the Baltic based on genomic evidence is important for conservation management of this endangered cetacean in threatened habitats, particularly in the Baltic Sea proper. In addition, we show that SNPs outperform microsatellite markers and demonstrate the utility of RAD-tags from a relatively small, opportunistically sampled cetacean sample set for population diversity and divergence analysis.

Postglacial climate changes and rise of three ecotypes of harbour porpoises, Phocoena phocoena, in western Palearctic waters.

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base-pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23-19 kilo-years ago, kyrbp). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp) than either does with the European continental shelf ecotype (c. 28 kyrbp), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post-glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.

A European melting pot of harbour porpoise in the French Atlantic coasts inferred from mitochondrial and nuclear data.

Field surveys have reported a global shift in harbour porpoise distribution in European waters during the last 15 years, including a return to the Atlantic coasts of France. In this study, we analyzed genetic polymorphisms at a fragment of the mitochondrial control region (mtDNA CR) and 7 nuclear microsatellite loci, for 52 animals stranded and by-caught between 2000 and 2010 along the Atlantic coasts of France. The analysis of nuclear and mitochondrial loci provided contrasting results. The mtDNA revealed two genetically distinct groups, one closely related to the Iberian and African harbour porpoises, and the second related to individuals from the more northern waters of Europe. In contrast, nuclear polymorphisms did not display such a distinction. Nuclear markers suggested that harbour porpoises behaved as a randomly mating population along the Atlantic coasts of France. The difference between the two kinds of markers can be explained by differences in their mode of inheritance, the mtDNA being maternally inherited in contrast to nuclear loci that are bi-parentally inherited. Our results provide evidence that a major proportion of the animals we sampled are admixed individuals from the two genetically distinct populations previously identified along the Iberian coasts and in the North East Atlantic. The French Atlantic coasts are clearly the place where these two previously separated populations of harbour porpoises are now admixing. The present shifts in distribution of harbour porpoises along this coast is likely caused by habitat changes that will need to be further studied.

Phenotypic and genetic divergence among harbour porpoise populations associated with habitat regions in the North Sea and adjacent seas.

Determining the mechanisms that generate population structure is essential to the understanding of speciation and the evolution of biodiversity. Here, we investigate a geographical range that transects two habitat gradients, the North Sea to North Atlantic transition, and the temperate to subpolar regions. We studied the harbour porpoise (Phocoena phocoena), a small odontocete inhabiting both subpolar and temperate waters. To assess differentiation among putative populations, we measured morphological variation at cranial traits (N = 462 individuals) and variation at eight microsatellite loci for 338 of the same individuals from Norwegian, British and Danish waters. Significant morphological differentiation reflected the size of the buccal cavity. Porpoises forage in relatively shallow waters preying mainly on benthic species in British and Danish waters, and on mesopelagic and pelagic fish off the coast of Norway. We suggest that the observed differentiation may be explained by resource specialization and either adaptation or developmental responses to different local habitats.

Genomic characterization of novel dolphin papillomaviruses provides indications for recombination within the Papillomaviridae.

Phylogenetic analysis of novel dolphin (Tursiops truncatus) papillomavirus sequences, TtPV1, -2, and -3, indicates that the early and late protein coding regions of their genomes differ in evolutionary history. Sliding window bootscan analysis showed a significant a change in phylogenetic clustering, in which the grouped sequences of TtPV1 and -3 move from a cluster with the Phocoena spinipinnis PsPV1 in the early region to a cluster with TtPV2 in the late region. This provides indications for a possible recombination event near the end of E2/beginning of L2. A second possible recombination site could be located near the end of L1, in the upstream regulatory region. Selection analysis by using maximum likelihood models of codon substitutions ruled out the possibility of intense selective pressure, acting asymmetrically on the viral genomes, as an alternative explanation for the observed difference in evolutionary history between the early and late genomic regions of these cetacean papillomaviruses.

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus).

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.

Rise of oceanographic barriers in continuous populations of a cetacean: the genetic structure of harbour porpoises in Old World waters.

BACKGROUND: Understanding the role of seascape in shaping genetic and demographic population structure is highly challenging for marine pelagic species such as cetaceans for which there is generally little evidence of what could effectively restrict their dispersal. In the present work, we applied a combination of recent individual-based landscape genetic approaches to investigate the population genetic structure of a highly mobile extensive range cetacean, the harbour porpoise in the eastern North Atlantic, with regards to oceanographic characteristics that could constrain its dispersal. RESULTS: Analyses of 10 microsatellite loci for 752 individuals revealed that most of the sampled range in the eastern North Atlantic behaves as a 'continuous' population that widely extends over thousands of kilometres with significant isolation by distance (IBD). However, strong barriers to gene flow were detected in the south-eastern part of the range. These barriers coincided with profound changes in environmental characteristics and isolated, on a relatively small scale, porpoises from Iberian waters and on a larger scale porpoises from the Black Sea. CONCLUSION: The presence of these barriers to gene flow that coincide with profound changes in oceanographic features, together with the spatial variation in IBD strength, provide for the first time strong evidence that physical processes have a major impact on the demographic and genetic structure of a cetacean. This genetic pattern further suggests habitat-related fragmentation of the porpoise range that is likely to intensify with predicted surface ocean warming.