The Novel Evolution of the Sperm Whale Genome

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.

[Evolution of cancer resistance in the animal kingdom]

Cancer is an inevitable collateral problem inherent in the evolution of multicellular organisms, which appeared at the end of the Precambrian. Faced to this constraint, a range of diverse anticancer defenses has evolved across the animal kingdom. Today, investigating how animal organisms, especially those of large size and long lifespan, manage cancer-related issues has both fundamental and applied outcomes, as it could inspire strategies for preventing or treating human cancers. In this article, we begin by presenting the conceptual framework for understanding evolutionary theories regarding the development of anti-cancer defenses. We then present a number of examples that have been extensively studied in recent years, including naked mole rats, elephants, whales, placozoa, xenarthras (such as sloths, armadillos and anteaters) and bats. The contributions of comparative genomics to understanding evolutionary convergences are also discussed. Finally, we emphasize that natural selection has also favored anti-cancer adaptations aimed at avoiding mutagenic environments, for example by maximizing immediate reproductive efforts in the event of cancer. Exploring these adaptive solutions holds promise for identifying novel approaches to improve human health.

Molecular evolution of wound healing-related genes during cetacean secondary aquatic adaptation

The marine environment presents challenges for wound healing in cetaceans, despite their remarkable recovery abilities with minimal infections or complications. However, the molecular mechanism underlying this efficient wound healing remains underexplored. To better understand the molecular mechanisms behind wound healing in cetaceans, we investigated the evolutionary patterns of 37 wound healing-related genes in representative mammals. We found wound healing-related genes experience adaptive evolution in cetaceans: (1) Three extrinsic coagulation pathway-related genes-tissue factor (F3), coagulation factor VII (F7), and coagulation factor X (F10)-are subject to positive selection in cetaceans, which might promote efficient hemostasis after injury; positive selection in transforming growth factor-beta 2 (TGF-β2), transforming growth factor-beta 3 (TGF-β3), and platelet-derived growth factor D (PDGFD), which play immunological roles in wound healing, may help cetaceans enhance inflammatory response and tissue debridement. (2) Coagulation factor XII (F12) is the initiation factor in the intrinsic coagulation pathway. It had a premature stop codon mutation and was subjected to selective stress relaxation in cetaceans, suggesting that the early termination of F12 may help cetaceans avoid the risk of vascular blockage during diving. (3) Fibrinogen alpha chain (FGA) and FIII, which were detected to contain the specific amino acid substitutions in marine mammals, indicating similar evolutionary mechanisms might exist among marine mammals to maintain strong wound-healing ability. Thus, our research provides further impetus to study the evolution of the wound healing system in cetaceans and other marine mammals, extending knowledge of preventing coagulation disorder and atherosclerosis in humans.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Comparative Analyses Reveal the Genetic Mechanism of Ambergris Production in the Sperm Whale Based on the Chromosome-Level Genome

Sperm whales are a marine mammal famous for the aromatic substance, the ambergris, produced from its colon. Little is known about the biological processes of ambergris production, and this study aims to investigate the genetic mechanism of ambergris production in the sperm whale based on its chromosome-level genome. Comparative genomics analyses found 1207 expanded gene families and 321 positive selected genes (PSGs) in the sperm whale, and functional enrichment analyses suggested revelatory pathways and terms related to the metabolism of steroids, terpenoids, and aldosterone, as well as microbiota interaction and immune network in the intestine. Furthermore, two sperm-whale-specific missense mutations (Tyr393His and Leu567Val) were detected in the PSG LIPE, which has been reported to play vital roles in lipid and cholesterol metabolism. In total, 46 CYP genes and 22 HSD genes were annotated, and then mapped to sperm whale chromosomes. Furthermore, phylogenetic analysis of CYP genes in six mammals found that CYP2E1, CYP51A and CYP8 subfamilies exhibited relative expansion in the sperm whale. Our results could help understand the genetic mechanism of ambergris production, and further reveal the convergent evolution pattern among animals that produce similar odorants.

Current global population size, post-whaling trend and historical trajectory of sperm whales

The sperm whale lives in most deep ice-free waters of the globe. It was targeted during two periods of whaling peaking in the 1840's and 1960's. Using a habitat suitability model, we extrapolated estimates of abundance from visual and acoustic surveys to give a global estimate of 736,053 sperm whales (CV = 0.218) in 1993. Estimates of trends in the post-whaling era suggest that: whaling, by affecting the sex ratio and/or the social cohesion of females, reduced recovery rates well after whaling ceased; preferentially-targeted adult males show the best evidence of recovery, presumably due to recruitment from breeding populations; several decades post-whaling, sperm whale populations not facing much human impact are recovering slowly, but populations may be declining in areas with substantial anthropogenic footprint. A theta-logistic population model enhanced to simulate spatial structure and the non-removal impacts of whaling indicated a pre-whaling population of 1,949,698 (CV = 0.178) in 1710 being reduced by whaling, and then then recovering a little to about 844,761 (CV = 0.209) in 2022. There is much uncertainty about these numbers and trends. A larger population estimate than produced by a similar analysis in 2002 is principally due to a better assessment of ascertainment bias.

Genomics reveal population structure, evolutionary history, and signatures of selection in the northern bottlenose whale, Hyperoodon ampullatus

Information on wildlife population structure, demographic history, and adaptations are fundamental to understanding species evolution and informing conservation strategies. To study this ecological context for a cetacean of conservation concern, we conducted the first genomic assessment of the northern bottlenose whale, Hyperoodon ampullatus, using whole-genome resequencing data (n = 37) from five regions across the North Atlantic Ocean. We found a range-wide pattern of isolation-by-distance with a genetic subdivision distinguishing three subgroups: the Scotian Shelf, western North Atlantic, and Jan Mayen regions. Signals of elevated levels of inbreeding in the Endangered Scotian Shelf population indicate this population may be more vulnerable than the other two subgroups. In addition to signatures of inbreeding, evidence of local adaptation in the Scotian Shelf was detected across the genome. We found a long-term decline in effective population size for the species, which poses risks to their genetic diversity and may be exacerbated by the isolating effects of population subdivision. Protecting important habitat and migratory corridors should be prioritized to rebuild population sizes that were diminished by commercial whaling, strengthen gene flow, and ensure animals can move across regions in response to environmental changes.

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.

The genome of the naturally evolved obesity-prone Ossabaw miniature pig

The feral pigs of Ossabaw Island (USA) have an outstanding propensity to obesity and develop complete metabolic syndrome (MetS) upon prolonged high energy dieting. We now report the first high quality genome of the Ossabaw pig with Contig N50 of ∼6.03 Mb, significantly higher than most other published pig genomes. Genomic comparison to Duroc reveals that variations including SNPs, INDELs and one ∼2 Mb inversion identified in Ossabaw pig may be related to its "thrifty" phenotype. Finally, an important positively selected gene (PSG) was found to be LEPR (leptin receptor) containing two positively selected sites which may lead to pseudogenization of this gene with possible significant effects on obesity and inflammation development. This work provides the first complete mapping of a genome representing a naturally 'feast and famine' evolved phenotype of MetS, serving as a blueprint to guide the search for new targets and new biomarkers for obesity comorbidities.

Integrated Full-Length Transcriptome and RNA-Seq to Identify Immune System Genes from the Skin of Sperm Whale (Physeter macrocephalus)

Cetaceans are a group of secondary aquatic mammals whose ancestors returned to the ocean from land, and during evolution, their immune systems adapted to the aquatic environment. Their skin, as the primary barrier to environmental pathogens, supposedly evolved to adapt to a new living environment. However, the immune system in the skin of cetaceans and the associated molecular mechanisms are still largely unknown. To better understand the immune system, we extracted RNA from the sperm whale's (Physeter macrocephalus) skin and performed PacBio full-length sequencing and RNA-seq sequencing. We obtained a total of 96,350 full-length transcripts with an average length of 1705 bp and detected 5150 genes that were associated with 21 immune-related pathways by gene annotation enrichment analysis. Moreover, we found 89 encoding genes corresponding to 33 proteins were annotated in the NOD-like receptor (NLR)-signaling pathway, including NOD1, NOD2, RIP2, and NF-kB genes, which were discussed in detail and predicted to play essential roles in the immune system of the sperm whale. Furthermore, NOD1 was highly conservative during evolution by the sequence comparison and phylogenetic tree. These results provide new information about the immune system in the skin of cetaceans, as well as the evolution of immune-related genes.

The evolution of factor XI and the kallikrein-kinin system

Factor XI (FXI) is the zymogen of a plasma protease (FXIa) that contributes to hemostasis by activating factor IX (FIX). In the original cascade model of coagulation, FXI is converted to FXIa by factor XIIa (FXIIa), a component, along with prekallikrein and high-molecular-weight kininogen (HK), of the plasma kallikrein-kinin system (KKS). More recent coagulation models emphasize thrombin as a FXI activator, bypassing the need for FXIIa and the KKS. We took an evolutionary approach to better understand the relationship of FXI to the KKS and thrombin generation. BLAST searches were conducted for FXI, FXII, prekallikrein, and HK using genomes for multiple vertebrate species. The analysis shows the KKS appeared in lobe-finned fish, the ancestors of all land vertebrates. FXI arose later from a duplication of the prekallikrein gene early in mammalian evolution. Features of FXI that facilitate efficient FIX activation are present in all living mammals, including primitive egg-laying monotremes, and may represent enhancement of FIX-activating activity inherent in prekallikrein. FXI activation by thrombin is a more recent acquisition, appearing in placental mammals. These findings suggest FXI activation by FXIIa may be more important to hemostasis in primitive mammals than in placental mammals. FXI activation by thrombin places FXI partially under control of the vitamin K-dependent coagulation mechanism, reducing the importance of the KKS in blood coagulation. This would explain why humans with FXI deficiency have a bleeding abnormality, whereas those lacking components of the KKS do not.

Genome sequencing of deep-sea hydrothermal vent snails reveals adaptions to extreme environments

BACKGROUND

The scaly-foot snail (Chrysomallon squamiferum) is highly adapted to deep-sea hydrothermal vents and has drawn much interest since its discovery. However, the limited information on its genome has impeded further related research and understanding of its adaptation to deep-sea hydrothermal vents.

FINDINGS

Here, we report the whole-genome sequencing and assembly of the scaly-foot snail and another snail (Gigantopelta aegis), which inhabits similar environments. Using Oxford Nanopore Technology, 10X Genomics, and Hi-C technologies, we obtained a chromosome-level genome of C. squamiferum with an N50 size of 20.71 Mb. By constructing a phylogenetic tree, we found that these 2 deep-sea snails evolved independently of other snails. Their divergence from each other occurred ∼66.3 million years ago. Comparative genomic analysis showed that different snails have diverse genome sizes and repeat contents. Deep-sea snails have more DNA transposons and long terminal repeats but fewer long interspersed nuclear elements than other snails. Gene family analysis revealed that deep-sea snails experienced stronger selective pressures than freshwater snails, and gene families related to the nervous system, immune system, metabolism, DNA stability, antioxidation, and biomineralization were significantly expanded in scaly-foot snails. We also found 251 H-2 Class II histocompatibility antigen, A-U α chain-like (H2-Aal) genes, which exist uniquely in the Gigantopelta aegis genome. This finding is important for investigating the evolution of major histocompatibility complex (MHC) genes.

CONCLUSION

Our study provides new insights into deep-sea snail genomes and valuable resources for further studies.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.]

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

•

Deducing chromosome evolution from ancestral Cetruminantia and ancestral Odontoceti

•

Reconstructing the demographic history of Sousa chinensis

•

Implying high risk of extinction and strong conservation requirement for S. chinensis

Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta)

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

BlobToolKit - Interactive Quality Assessment of Genome Assemblies

Reconstruction of target genomes from sequence data produced by instruments that are agnostic as to the species-of-origin may be confounded by contaminant DNA. Whether introduced during sample processing or through co-extraction alongside the target DNA, if insufficient care is taken during the assembly process, the final assembled genome may be a mixture of data from several species. Such assemblies can confound sequence-based biological inference and, when deposited in public databases, may be included in downstream analyses by users unaware of underlying problems. We present BlobToolKit, a software suite to aid researchers in identifying and isolating non-target data in draft and publicly available genome assemblies. BlobToolKit can be used to process assembly, read and analysis files for fully reproducible interactive exploration in the browser-based Viewer. BlobToolKit can be used during assembly to filter non-target DNA, helping researchers produce assemblies with high biological credibility. We have been running an automated BlobToolKit pipeline on eukaryotic assemblies publicly available in the International Nucleotide Sequence Data Collaboration and are making the results available through a public instance of the Viewer at https://blobtoolkit.genomehubs.org/view. We aim to complete analysis of all publicly available genomes and then maintain currency with the flow of new genomes. We have worked to embed these views into the presentation of genome assemblies at the European Nucleotide Archive, providing an indication of assembly quality alongside the public record with links out to allow full exploration in the Viewer.

Differential mechanisms of tolerance to extreme environmental conditions in tardigrades

Tardigrades, also known as water bears, are small aquatic animals that inhabit marine, fresh water or limno-terrestrial environments. While all tardigrades require surrounding water to grow and reproduce, species living in limno-terrestrial environments (e.g. Ramazzottius varieornatus) are able to undergo almost complete dehydration by entering an arrested state known as anhydrobiosis, which allows them to tolerate ionic radiation, extreme temperatures and intense pressure. Previous studies based on comparison of the genomes of R. varieornatus and Hypsibius dujardini - a less tolerant tardigrade - have pointed to potential mechanisms that may partially contribute to their remarkable ability to resist extreme physical conditions. In this work, we have further annotated the genomes of both tardigrades using a guided approach in search for novel mechanisms underlying the extremotolerance of R. varieornatus. We have found specific amplifications of several genes, including MRE11 and XPC, and numerous missense variants exclusive of R. varieornatus in CHEK1, POLK, UNG and TERT, all of them involved in important pathways for DNA repair and telomere maintenance. Taken collectively, these results point to genomic features that may contribute to the enhanced ability to resist extreme environmental conditions shown by R. varieornatus.

Genes lost during the transition from land to water in cetaceans highlight genomic changes associated with aquatic adaptations

The transition from land to water in whales and dolphins (cetaceans) was accompanied by remarkable adaptations. To reveal genomic changes that occurred during this transition, we screened for protein-coding genes that were inactivated in the ancestral cetacean lineage. We found 85 gene losses. Some of these were likely beneficial for cetaceans, for example, by reducing the risk of thrombus formation during diving (F12 and KLKB1), erroneous DNA damage repair (POLM), and oxidative stress-induced lung inflammation (MAP3K19). Additional gene losses may reflect other diving-related adaptations, such as enhanced vasoconstriction during the diving response (mediated by SLC6A18) and altered pulmonary surfactant composition (SEC14L3), while loss of SLC4A9 relates to a reduced need for saliva. Last, loss of melatonin synthesis and receptor genes (AANAT, ASMT, and MTNR1A/B) may have been a precondition for adopting unihemispheric sleep. Our findings suggest that some genes lost in ancestral cetaceans were likely involved in adapting to a fully aquatic lifestyle.

Whole Genome Sequencing of Chinese White Dolphin (Sousa chinensis) for High-Throughput Screening of Antihypertensive Peptides

Chinese white dolphin (Sousa chinensis), also known as the Indo-Pacific humpback dolphin, has been classified as “Vulnerable” on the IUCN Red List of Threatened Species. It is a special cetacean species that lives in tropical and subtropical nearshore waters, with significant differences from other cetaceans. Here, we sequenced and assembled a draft genome of the Chinese white dolphin with a total length of 2.3 Gb and annotation of 18,387 protein-coding genes. Genes from certain expanded families are potentially involved in DNA replication and repairing, suggesting that they may be related to adaptation of this marine mammal to nearshore environments. We also discovered that its historical population had undergone a remarkable bottleneck incident before the Mindel glaciation. In addition, a comparative genomic survey on antihypertensive peptides (AHTPs) among five representative mammals with various residential habitats (such as remarkable differences in exogenous ion concentrations and sea depth) revealed that these small bioactive peptides were highly conserved among these examined mammals, and they had the most abundant hits in collagen subunit proteins, especially for two putative AHTP peptides Gly-Leu-Pro (GLP) and Leu-Gly-Pro (LGP). Our genome assembly will be a valuable resource for further genetic researches on adaptive ecology and conservation biology of cetaceans, and for in-depth investigations into bioactive peptides in aquatic and terrestrial mammals for development of peptide-based drugs to treat various human cardiovascular diseases.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae)

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

A survey of the sperm whale (Physeter catodon) commensal microbiome

Background

Mammalian commensal microbiota play important roles in the health of its host. In comparison to terrestrial mammals, commensal microbiota of marine mammals is mainly focused on the composition and function of skin and gut microbiota, with less attention paid to the health impact of bacteria and viruses. Previous studies on sperm whales (Physeter catodon) have affirmed their important phylogenetic position; however, studies on their commensal microbiota have not been published, due to difficulty in sample collection.

Methods

Here, we sequenced the metagenomes of blood, muscle and fecal samples from a stranded sperm whale using the BGISEQ-500 platform. We compared the diversity and abundance of microbiomes from three different tissues and tried to search pathogenic bacterial and virulence genes probably related to the health of the sperm whale. We also performed 16S rDNA sequencing of the fecal sample to compare to published gut metagenome data from other marine mammals.

Results

Our results demonstrated notable differences in species richness and abundance in the three samples. Extensive bacteria, including Enterococcus faecium, Fusobacterium nucleatum, Pseudomonas aeruginosa, Streptococcus anginosus, Streptococcus pneumoniae, and Streptococcus suis, and five toxigenic Clostridium species usually associated with infection, were found in the three samples. We also found the taxa composition of sperm whale gut microbiota was similar to that of other whales, suggesting co-evolution with its host. This study is the first report of the sperm whale gut microbiome, and provides a foundation for the pathogen detection and health assessment of the sperm whale.

Differential Evolution of the Epidermal Keratin Cytoskeleton in Terrestrial and Aquatic Mammals

Keratins are the main intermediate filament proteins of epithelial cells. In keratinocytes of the mammalian epidermis they form a cytoskeleton that resists mechanical stress and thereby are essential for the function of the skin as a barrier against the environment. Here, we performed a comparative genomics study of epidermal keratin genes in terrestrial and fully aquatic mammals to determine adaptations of the epidermal keratin cytoskeleton to different environments. We show that keratins K5 and K14 of the innermost (basal), proliferation-competent layer of the epidermis are conserved in all mammals investigated. In contrast, K1 and K10, which form the main part of the cytoskeleton in the outer (suprabasal) layers of the epidermis of terrestrial mammals, have been lost in whales and dolphins (cetaceans) and in the manatee. Whereas in terrestrial mammalian epidermis K6 and K17 are expressed only upon stress-induced epidermal thickening, high levels of K6 and K17 are consistently present in dolphin skin, indicating constitutive expression and substitution of K1 and K10. K2 and K9, which are expressed in a body site-restricted manner in human and mouse suprabasal epidermis, have been lost not only in cetaceans and manatee but also in some terrestrial mammals. The evolution of alternative splicing of K10 and differentiation-dependent upregulation of K23 have increased the complexity of keratin expression in the epidermis of terrestrial mammals. Taken together, these results reveal evolutionary diversification of the epidermal cytoskeleton in mammals and suggest a complete replacement of the quantitatively predominant epidermal proteins of terrestrial mammals by originally stress-inducible keratins in cetaceans.

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.

Transition to an Aquatic Habitat Permitted the Repeated Loss of the Pleiotropic KLK8 Gene in Mammals

Kallikrein related peptidase 8 (KLK8; also called neuropsin) is a serine protease that plays distinct roles in the skin and hippocampus. In the skin, KLK8 influences keratinocyte proliferation and desquamation, and activates antimicrobial peptides in sweat. In the hippocampus, KLK8 affects memory acquisition. Here, we examined the evolution of KLK8 in mammals and discovered that, out of 70 placental mammals, KLK8 is exclusively lost in three independent fully-aquatic lineages, comprising dolphin, killer whale, minke whale, and manatee. In addition, while the sperm whale has an intact KLK8 reading frame, the gene evolves neutrally in this species. We suggest that the distinct functions of KLK8 likely became obsolete in the aquatic environment, leading to the subsequent loss of KLK8 in several fully-aquatic mammalian lineages. First, the cetacean and manatee skin lacks sweat glands as an adaptation to the aquatic environment, which likely made the epidermal function of KLK8 obsolete. Second, cetaceans and manatees exhibit a proportionally small hippocampus, which may have rendered the hippocampal functions of KLK8 obsolete. Together, our results shed light on the genomic changes that correlate with skin and neuroanatomical differences of aquatic mammals, and show that even pleiotropic genes can be lost during evolution if an environmental change nullifies the need for the different functions of such genes.