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

Molecular Evidence for Relaxed Selection on the Enamel Genes of Toothed Whales (Odontoceti) with Degenerative Enamel Phenotypes

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development.

Guidelines and quantitative standards for improved cetacean taxonomy using full mitochondrial genomes

In many organisms, especially those of conservation concern, traditional lines of evidence for taxonomic delineation, such as morphological data, are often difficult to obtain. In these cases, genetic data are often the only source of information available for taxonomic studies. In particular, population surveys of mitochondrial genomes offer increased resolution and precision in support of taxonomic decisions relative to conventional use of the control region or other gene fragments of the mitochondrial genome. To improve quantitative guidelines for taxonomic decisions in cetaceans, we build on a previous effort targeting the control region and evaluate, for whole mitogenome sequences, a suite of divergence and diagnosability estimates for pairs of recognized cetacean populations, subspecies, and species. From this overview, we recommend new guidelines based on complete mitogenomes, combined with other types of evidence for isolation and divergence, which will improve resolution for taxonomic decisions, especially in the face of small sample sizes or low levels of genetic diversity. We further use simulated data to assist interpretations of divergence in the context of varying forms of historical demography, culture, and ecology.

Population Structure of Finless Porpoise (Neophocaena phocaenoides) Discovered off Coastal Waters, Republic of Korea

The finless porpoise (Neophocaena phocaenoides Cuvier, 1829) is distributed in the coastal waters of Asia, throughout Indonesia to the east, and as far north as the Taiwan Strait. The finless porpoise has been declared critically endangered by the WWF (World Wide Fund for Nature), and in 2017 was rated vulnerable on the IUCN Red Threatened Species List. Since this species is distributed near the coast and has many interactions with humans, effective conservation of the species requires further studies into their genetic diversity and population. In this study, 45 samples were obtained from bycatch or stranded individuals in the East, South, and West Seas, where Korean porpoises were mainly distributed from 2017–2021. We compared 473 bp mtDNA sequences from the control region. Pairwise fixation indices (F_ST) revealed that the two populations differed significantly (F_ST = 0.4557, p = 0.000). In contrast to high levels of genetic differentiation, gene flow was identified as medium levels (Nm of 0.04–0.71). Our data suggest that finless porpoises may have undergone a historic differentiation event, and that finless porpoises in the three regions could be divided into two populations: West and East/South.

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.

Vertebral morphology in extant porpoises: radiation and functional implications

Vertebral morphology has profound biomechanical implications and plays an important role in adaptation to different habitats and foraging strategies for cetaceans. Extant porpoise species (Phocoenidae) display analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. We employed 3D geometric morphometrics to study vertebral morphology in five porpoise species with contrasting habitats: the coastal Indo-Pacific finless porpoise (Neophocaena phocaenoides); the mostly coastal harbor porpoise (Phocoena phocoena) and Burmeister's porpoise (Phocoena spinipinnis); and the oceanic spectacled porpoise (Phocoena dioptrica) and Dall's porpoise (Phocoenoides dalli). We evaluated the radiation of vertebral morphology, both in size and shape, using multivariate statistics. We supplemented data with samples of an early-radiating delphinoid species, the narwhal (Monodon monoceros); and an early-radiating delphinid species, the white-beaked dolphin (Lagenorhynchus albirostris). Principal component analyses were used to map shape variation onto phylogenies, and phylogenetic constraints were investigated through permutation tests. We established links between vertebral morphology and movement patterns through biomechanical inferences from morphological presentations. We evidenced divergence in size between species with contrasting habitats, with coastal species tending to decrease in size from their estimated ancestral state, and oceanic species tending to increase in size. Regarding vertebral shape, coastal species had longer centra and shorter neural processes, but longer transverse processes, whilst oceanic species tended to have disk-shaped vertebrae with longer neural processes. Within Phocoenidae, the absence of phylogenetic constraints in vertebral morphology suggests a high level of evolutionary lability. Overall, our results are in accordance with the hypothesis of speciation within the family from a coastal ancestor, through adaptation to particular habitats. Variation in vertebral morphology in this group of small odontocetes highlights the importance of environmental complexity and particular selective pressures for the speciation process through the development of adaptations that minimize energetic costs during locomotion and prey capture. This article is protected by copyright. All rights reserved.

No leading-edge effect in North Atlantic harbor porpoises: Evolutionary and conservation implications

Understanding species responses to past environmental changes can help forecast how they will cope with ongoing climate changes. Harbor porpoises are widely distributed in the North Atlantic and were deeply impacted by the Pleistocene changes with the split of three subspecies. Despite major impacts of fisheries on natural populations, little is known about population connectivity and dispersal, how they reacted to the Pleistocene changes, and how they will evolve in the future. Here, we used phylogenetics, population genetics, and predictive habitat modeling to investigate population structure and phylogeographic history of the North Atlantic porpoises. A total of 925 porpoises were characterized at 10 microsatellite loci and one quarter of the mitogenome (mtDNA). A highly divergent mtDNA lineage was uncovered in one porpoise off Western Greenland, suggesting that a cryptic group may occur and could belong to a recently discovered mesopelagic ecotype off Greenland. Aside from it and the southern subspecies, spatial genetic variation showed that porpoises from both sides of the North Atlantic form a continuous system belonging to the same subspecies (Phocoena phocoena phocoena). Yet, we identified important departures from random mating and restricted dispersal forming a highly significant isolation by distance (IBD) at both mtDNA and nuclear markers. A ten times stronger IBD at mtDNA compared with nuclear loci supported previous evidence of female philopatry. Together with the lack of spatial trends in genetic diversity, this IBD suggests that migration-drift equilibrium has been reached, erasing any genetic signal of a leading-edge effect that accompanied the predicted recolonization of the northern habitats freed from Pleistocene ice. These results illuminate the processes shaping porpoise population structure and provide a framework for designing conservation strategies and forecasting future population evolution.

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 (F_ST  = 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.

Fresh water skin disease in dolphins: a case definition based on pathology and environmental factors in Australia

A distinct ulcerative dermatitis known as “freshwater skin disease” is an emerging clinical and pathological presentation in coastal cetaceans worldwide. In Australia, two remarkably similar mortality events enabled the creation of a case definition based on pathology and environmental factors. The first affected a community of endemic Tursiops australis in the Gippsland Lakes, Victoria, while the second occurred among T. aduncus resident in the Swan-Canning River system, Western Australia. The common features of both events were (1) an abrupt and marked decrease in salinity (from > 30ppt to < 5ppt) due to rainfall in the catchments, with hypo-salinity persisting weeks to months, and (2) dermatitis characterized grossly by patchy skin pallor that progressed to variable circular or targetoid, often raised, and centrally ulcerated lesions covering up to 70% of the body surface. The affected skin was often colonized by a variety of fungal, bacterial and algal species that imparted variable yellow, green or orange discoloration. Histologic lesions consisted of epidermal hydropic change leading to vesiculation and erosion; alternately, or in addition, the formation of intra-epithelial pustules resulting in ulceration and hypodermal necrosis. Thus, the environmental factors and characteristic pathologic lesions, are necessary components of the case definition for freshwater skin disease.

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.