Why are whales big? Genes behind ocean giants.

Gigantism is prevalent in animals, but it has never reached more extreme levels than in aquatic mammals such as whales, dolphins, and porpoises. A new study by Silva et al. has uncovered five genes underlying this gigantism, a phenotype with important connections to aging and cancer suppression in long-lived animals.

Ancient whale rhodopsin reconstructs dim-light vision over a major evolutionary transition: Implications for ancestral diving behavior.

Cetaceans are fully aquatic mammals that descended from terrestrial ancestors, an iconic evolutionary transition characterized by adaptations for underwater foraging via breath-hold diving. Although the evolutionary history of this specialized behavior is challenging to reconstruct, coevolving sensory systems may offer valuable clues. The dim-light visual pigment, rhodopsin, which initiates phototransduction in the rod photoreceptors of the eye, has provided insight into the visual ecology of depth in several aquatic vertebrate lineages. Here, we use ancestral sequence reconstruction and protein resurrection experiments to quantify light-activation metrics in rhodopsin pigments from ancestors bracketing the cetacean terrestrial-to-aquatic transition. By comparing multiple reconstruction methods on a broadly sampled cetartiodactyl species tree, we generated highly robust ancestral sequence estimates. Our experimental results provide direct support for a blue-shift in spectral sensitivity along the branch separating cetaceans from terrestrial relatives. This blue-shift was 14 nm, resulting in a deep-sea signature (λmax = 486 nm) similar to many mesopelagic-dwelling fish. We also discovered that the decay rates of light-activated rhodopsin increased in ancestral cetaceans, which may indicate an accelerated dark adaptation response typical of deeper-diving mammals. Because slow decay rates are thought to help sequester cytotoxic photoproducts, this surprising result could reflect an ecological trade-off between rod photoprotection and dark adaptation. Taken together, these ancestral shifts in rhodopsin function suggest that some of the first fully aquatic cetaceans could dive into the mesopelagic zone (>200 m). Moreover, our reconstructions indicate that this behavior arose before the divergence of toothed and baleen whales.

Conservation and selective pressures shaping baleen whale olfactory receptor genes supports their use of olfaction in the marine environment.

The relative importance of various sensory modalities can shift in response to evolutionary transitions, resulting in changes to underlying gene families encoding their reception systems. The rapid birth-and-death process underlying the evolution of the large olfactory receptor (OR) gene family has accelerated genomic-level change for the sense of smell in particular. The transition from the land to sea in marine mammals is an attractive model for understanding the influence of habitat shifts on sensory systems, with the retained OR repertoire of baleen whales contrasting with its loss in toothed whales. In this study, we examine to what extent the transition from a terrestrial to a marine environment has influenced the evolution of baleen whale OR repertoires. We developed Gene Mining Pipeline (GMPipe) (https://github.com/AprilJauhal/GMPipe), which can accurately identify large numbers of candidate OR genes. GMPipe identified 707 OR sequences from eight baleen whale species. These repertoires exhibited distinct family count distributions compared to terrestrial mammals, including signs of relative expansion in families OR10, OR11 and OR13. While many receptors have been lost or show signs of random drift in baleen whales, others exhibit signs of evolving under purifying or positive selection. Over 85% of OR genes could be sorted into orthologous groups of sequences containing at least four homologous sequences. Many of these groups, particularly from family OR10, presented signs of relative expansion and purifying selective pressure. Overall, our results suggest that the relatively small size of baleen whale OR repertoires result from specialisation to novel olfactory landscapes, as opposed to random drift.

Molecular identification of whales remains from the Keller Peninsula, Admiralty Bay, King George Island, Antarctica.

At the beginning of the 20th century, intense whaling activity took place in the South Shetland Islands, which is represented today in the form of ruins and numerous whale bones scattered along several Antarctic beaches. Despite being exposed to a harsh environment throughout the last decades, the present manuscript tried to answer if these bone remains still have viable DNA to allow species' identification using molecular methods. Several individuals were collected from the shores of Keller Peninsula, Admiralty Bay, Antarctica, and submitted to DNA extraction, amplification and Sanger sequencing. The challenging identification of these bone fragments proved to be still feasible. Mitochondrial DNA was successfully extracted, amplified and sequenced. A database with 43 sequences including previously published and newly determined sequences were built and enabled the precise identification to species level for some of the collected samples, therefore shedding light on the whales species that inhabited the region and how their overexploitation seems to have affected modern day presence of these species within the study area.

Genomic Impact of Whaling in North Atlantic Fin Whales.

It is generally recognized that large-scale whaling in the 19th and 20th century led to a substantial reduction of the size of many cetacean populations, particularly those of the baleen whales (Mysticeti). The impact of these operations on genomic diversity of one of the most hunted whales, the fin whale (Balaenoptera physalus), has remained largely unaddressed because of the paucity of adequate samples and the limitation of applicable techniques. Here, we have examined the effect of whaling on the North Atlantic fin whale based on genomes of 51 individuals from Icelandic waters, representing three temporally separated intervals, 1989, 2009 and 2018 and provide a reference genome for the species. Demographic models suggest a noticeable drop of the effective population size of the North Atlantic fin whale around a century ago. The present results suggest that the genome-wide heterozygosity is not markedly reduced and has remained comparable with other baleen whale species. Similarly, there are no signs of apparent inbreeding, as measured by the proportion of long runs of homozygosity, or of a distinctively increased mutational load, as measured by the amount of putative deleterious mutations. Compared with other baleen whales, the North Atlantic fin whale appears to be less affected by anthropogenic influences than other whales such as the North Atlantic right whale, consistent with the presence of long runs of homozygosity and higher levels of mutational load in an otherwise more heterozygous genome. Thus, genome-wide assessments of other species and populations are essential for future, more specific, conservation efforts.

The multicomparative 2-n-way genome suite.

To effectively analyze the increasing amounts of available genomic data, improved comparative analytical tools that are accessible to and applicable by a broad scientific community are essential. We built the "2-n-way" software suite to provide a fundamental and innovative processing framework for revealing and comparing inserted elements among various genomes. The suite comprises two user-friendly web-based modules. The 2-way module generates pairwise whole-genome alignments of target and query species. The resulting genome coordinates of blocks (matching sequences) and gaps (missing sequences) from multiple 2-ways are then transferred to the n-way module and sorted into projects, in which user-defined coordinates from reference species are projected to the block/gap coordinates of orthologous loci in query species to provide comparative information about presence (blocks) or absence (gaps) patterns of targeted elements over many entire genomes and phylogroups. Thus, the 2-n-way software suite is ideal for performing multidirectional, non-ascertainment-biased screenings to extract all possible presence/absence data of user-relevant elements in orthologous sequences. To highlight its applicability and versatility, we used 2-n-way to expose approximately 100 lost introns in vertebrates, analyzed thousands of potential phylogenetically informative bat and whale retrotransposons, and novel human exons as well as thousands of human polymorphic retrotransposons.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea.

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

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

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

Ancestry testing of "Old Tom," a killer whale central to mutualistic interactions with human whalers.

Cooperative hunting between humans and killer whales (Orcinus orca) targeting baleen whales was reported in Eden, New South Wales, Australia, for almost a century. By 1928, whaling operations had ceased, and local killer whale sightings became scarce. A killer whale from the group, known as "Old Tom," washed up dead in 1930 and his skeleton was preserved. How these killer whales from Eden relate to other populations globally and whether their genetic descendants persist today remains unknown. We extracted and sequenced DNA from Old Tom using ancient DNA techniques. Genomic sequences were then compared with a global dataset of mitochondrial and nuclear genomes. Old Tom shared a most recent common ancestor with killer whales from Australasia, the North Atlantic, and the North Pacific, having the highest genetic similarity with contemporary New Zealand killer whales. However, much of the variation found in Old Tom's genome was not shared with these widespread populations, suggesting ancestral rather than ongoing gene flow. Our genetic comparisons also failed to find any clear descendants of Tom, raising the possibility of local extinction of this group. We integrated Traditional Custodian knowledge to recapture the events in Eden and recognize that Indigenous Australians initiated the relationship with the killer whales before European colonization and the advent of commercial whaling locally. This study rectifies discrepancies in local records and provides new insight into the origins of the killer whales in Eden and the history of Australasian killer whales.

Differential MC5R loss in whales and manatees reveals convergent evolution to the marine environment.

Melanocortin 5 receptor (MC5R), which is expressed in the terminally differentiated sebaceous gland, is a G protein-coupled receptor (GPCR). MC5R exists mostly in mammals but is completely lost in whales; only the relic of MC5R can be detected in manatees, and phenotypically, they have lost sebaceous glands. Interestingly, whales and manatees are both aquatic mammals but have no immediate common ancestors. The loss of MC5R and sebaceous glands in whales and manatees is likely to be a result of convergent evolution. Here, we find that MC5R in whales and manatees are lost by two different mechanisms. Homologous recombination of MC5R in manatees and the insertion of reverse transcriptase in whales lead to the gene loss, respectively. On one hand, in manatees, there are two "TTATC" sequences flanking MC5R, and homologous recombination of the segments between the two "TTATC" sequences resulted in the partial loss of the sequence of MC5R. On the other hand, in whales, reverse transcriptase inserts between MC2R and RNMT on the chromosome led to the loss of MC5R. Based on these two different mechanisms for gene loss in whales and manatees, we finally concluded that MC5R loss might be the result of convergent evolution to the marine environment, and we explored the impact on biological function that is significant to environmental adaptation.

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.

Molecular evolutionary analyses of tooth genes support sequential loss of enamel and teeth in baleen whales (Mysticeti).

The loss of teeth and evolution of baleen racks in Mysticeti was a profound transformation that permitted baleen whales to radiate and diversify into a previously underutilized ecological niche of bulk filter-feeding on zooplankton and other small prey. Ancestral state reconstructions suggest that postnatal teeth were lost in the common ancestor of crown Mysticeti. Genomic studies provide some support for this hypothesis and suggest that the genetic toolkit for enamel production was inactivated in the common ancestor of living baleen whales. However, molecular studies to date have not provided direct evidence for the complete loss of teeth, including their dentin component, on the stem mysticete branch. Given these results, several questions remain unanswered: (1) Were teeth lost in a single step or did enamel loss precede dentin loss? (2) Was enamel lost early or late on the stem mysticete branch? (3) If enamel and dentin/tooth loss were decoupled in the ancestry of baleen whales, did dentin loss occur on the stem mysticete branch or independently in different crown mysticete lineages? To address these outstanding questions, we compiled and analyzed complete protein-coding sequences for nine tooth-related genes from cetaceans with available genome data. Seven of these genes are associated with enamel formation (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) whereas two other genes are either dentin-specific (DSPP) or tooth-specific (ODAPH) but not enamel-specific. Molecular evolutionary analyses indicate that all seven enamel-specific genes have inactivating mutations that are scattered across branches of the mysticete tree. Three of the enamel genes (ACP4, KLK4, MMP20) have inactivating mutations that are shared by all mysticetes. The two genes that are dentin-specific (DSPP) or tooth-specific (ODAPH) do not have any inactivating mutations that are shared by all mysticetes, but there are shared mutations in Balaenidae as well as in Plicogulae (Neobalaenidae + Balaenopteroidea). These shared mutations suggest that teeth were lost at most two times. Shared inactivating mutations and dN/dS analyses, in combination with cetacean divergence times, were used to estimate inactivation times of genes and by proxy enamel and tooth phenotypes at ancestral nodes. The results of these analyses are most compatible with a two-step model for the loss of teeth in the ancestry of living baleen whales: enamel was lost very early on the stem Mysticeti branch followed by the independent loss of dentin (and teeth) in the common ancestors of Balaenidae and Plicogulae, respectively. These results imply that some stem mysticetes, and even early crown mysticetes, may have had vestigial teeth comprised of dentin with no enamel. Our results also demonstrate that all odontocete species (in our study) with absent or degenerative enamel have inactivating mutations in one or more of their enamel genes.

A functional enrichment test for molecular convergent evolution finds a clear protein-coding signal in echolocating bats and whales.

Distantly related species entering similar biological niches often adapt by evolving similar morphological and physiological characters. How much genomic molecular convergence (particularly of highly constrained coding sequence) contributes to convergent phenotypic evolution, such as echolocation in bats and whales, is a long-standing fundamental question. Like others, we find that convergent amino acid substitutions are not more abundant in echolocating mammals compared to their outgroups. However, we also ask a more informative question about the genomic distribution of convergent substitutions by devising a test to determine which, if any, of more than 4,000 tissue-affecting gene sets is most statistically enriched with convergent substitutions. We find that the gene set most overrepresented (q-value = 2.2e-3) with convergent substitutions in echolocators, affecting 18 genes, regulates development of the cochlear ganglion, a structure with empirically supported relevance to echolocation. Conversely, when comparing to nonecholocating outgroups, no significant gene set enrichment exists. For aquatic and high-altitude mammals, our analysis highlights 15 and 16 genes from the gene sets most affected by molecular convergence which regulate skin and lung physiology, respectively. Importantly, our test requires that the most convergence-enriched set cannot also be enriched for divergent substitutions, such as in the pattern produced by inactivated vision genes in subterranean mammals. Showing a clear role for adaptive protein-coding molecular convergence, we discover nearly 2,600 convergent positions, highlight 77 of them in 3 organs, and provide code to investigate other clades across the tree of life.

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

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

The Karyotype of Blainville's Beaked Whale, Mesoplodon densirostris.

The karyotype of the Odontocete whale, Mesoplodon densirostris, has not been previously reported. The chromosome number is determined to be 2n = 42, and the karyotype is presented using G-, C-, and nucleolar organizer region (NOR) banding. The findings include NOR regions on 2 chromosomes, regions of heterochromatic variation, a large block of heterochromatin on the X chromosome, and a relatively large Y chromosome. The karyotype is compared to published karyograms of 2 other species of Mesoplodon.

Linking gene expression and phenotypic changes in the developmental and evolutionary origins of osteosclerosis in the ribs of bowhead whales (Balaena mysticetus).

Bowhead whales are among the longest-lived mammals with an extreme lifespan of about 211 years. During the first 25 years of their lives, rib bones increase in mineral density and the medulla transitions from compact to trabecular bone. Molecular drivers associated with these phenotypic changes in bone remain unknown. This study assessed expression levels of osteogenic genes from samples of rib bones of bowheads. Samples were harvested from prenatal to 86-year-old whales, representing the first third of the bowhead lifespan. Fetal to 2-year-old bowheads showed expression levels consistent with the rapid deposition of the bone extracellular matrix. Sexually mature animals showed expression levels associated with low rates of osteogenesis and increased osteoclastogenesis. After the first 25 years of life, declines in osteogenesis corresponded with increased expression of EZH2, an epigenetic regulator of osteogenesis. These findings suggest EZH2 may be at least one epigenetic modifier that contributes to the age-related changes in the rib bone phenotype along with the transition from compact to trabecular bone. Ancient cetaceans and their fossil relatives also display these phenotypes, suggesting EZH2 may have shaped the skeleton of whales in evolutionary history.

DNA preserved in jetsam whale ambergris.

Jetsam ambergris, found on beaches worldwide, has always been assumed to originate as a natural product of sperm whales (Physeteroidea). However, only indirect evidence has ever been produced for this, such as the presence of whale prey remains in ambergris. Here, we extracted and analysed DNA sequences from jetsam ambergris from beaches in New Zealand and Sri Lanka, and sequences from ambergris of a sperm whale beached in The Netherlands. The lipid-rich composition of ambergris facilitated high preservation-quality of endogenous DNA, upon which we performed shotgun Illumina sequencing. Alignment of mitochondrial and nuclear genome sequences with open-access reference data for multiple whale species confirms that all three jetsam samples derived originally from sperm whales (Physeter macrocephalus). Shotgun sequencing here also provides implications for metagenomic insights into ambergris-preserved DNA. These results demonstrate significant implications for elucidating the origins of jetsam ambergris as a prized natural product, and also for the understanding of sperm whale metabolism and diet, and the ecological mechanisms underlying these coproliths.

Genetic Diversity and Connectivity of Southern Right Whales (Eubalaena australis) Found in the Brazil and Chile-Peru Wintering Grounds and the South Georgia (Islas Georgias del Sur) Feeding Ground.

As species recover from exploitation, continued assessments of connectivity and population structure are warranted to provide information for conservation and management. This is particularly true in species with high dispersal capacity, such as migratory whales, where patterns of connectivity could change rapidly. Here we build on a previous long-term, large-scale collaboration on southern right whales (Eubalaena australis) to combine new (nnew) and published (npub) mitochondrial (mtDNA) and microsatellite genetic data from all major wintering grounds and, uniquely, the South Georgia (Islas Georgias del Sur: SG) feeding grounds. Specifically, we include data from Argentina (npub mtDNA/microsatellite = 208/46), Brazil (nnew mtDNA/microsatellite = 50/50), South Africa (nnew mtDNA/microsatellite = 66/77, npub mtDNA/microsatellite = 350/47), Chile-Peru (nnew mtDNA/microsatellite = 1/1), the Indo-Pacific (npub mtDNA/microsatellite = 769/126), and SG (npub mtDNA/microsatellite = 8/0, nnew mtDNA/microsatellite = 3/11) to investigate the position of previously unstudied habitats in the migratory network: Brazil, SG, and Chile-Peru. These new genetic data show connectivity between Brazil and Argentina, exemplified by weak genetic differentiation and the movement of 1 genetically identified individual between the South American grounds. The single sample from Chile-Peru had an mtDNA haplotype previously only observed in the Indo-Pacific and had a nuclear genotype that appeared admixed between the Indo-Pacific and South Atlantic, based on genetic clustering and assignment algorithms. The SG samples were clearly South Atlantic and were more similar to the South American than the South African wintering grounds. This study highlights how international collaborations are critical to provide context for emerging or recovering regions, like the SG feeding ground, as well as those that remain critically endangered, such as Chile-Peru.

Odontogenic ameloblast-associated (ODAM) is inactivated in toothless/enamelless placental mammals and toothed whales.

BACKGROUND: The gene for odontogenic ameloblast-associated (ODAM) is a member of the secretory calcium-binding phosphoprotein gene family. ODAM is primarily expressed in dental tissues including the enamel organ and the junctional epithelium, and may also have pleiotropic functions that are unrelated to teeth. Here, we leverage the power of natural selection to test competing hypotheses that ODAM is tooth-specific versus pleiotropic. Specifically, we compiled and screened complete protein-coding sequences, plus sequences for flanking intronic regions, for ODAM in 165 placental mammals to determine if this gene contains inactivating mutations in lineages that either lack teeth (baleen whales, pangolins, anteaters) or lack enamel on their teeth (aardvarks, sloths, armadillos), as would be expected if the only essential functions of ODAM are related to tooth development and the adhesion of the gingival junctional epithelium to the enamel tooth surface. RESULTS: We discovered inactivating mutations in all species of placental mammals that either lack teeth or lack enamel on their teeth. A surprising result is that ODAM is also inactivated in a few additional lineages including all toothed whales that were examined. We hypothesize that ODAM inactivation is related to the simplified outer enamel surface of toothed whales. An alternate hypothesis is that ODAM inactivation in toothed whales may be related to altered antimicrobial functions of the junctional epithelium in aquatic habitats. Selection analyses on ODAM sequences revealed that the composite dN/dS value for pseudogenic branches is close to 1.0 as expected for a neutrally evolving pseudogene. DN/dS values on transitional branches were used to estimate ODAM inactivation times. In the case of pangolins, ODAM was inactivated ~ 65 million years ago, which is older than the oldest pangolin fossil (Eomanis, 47 Ma) and suggests an even more ancient loss or simplification of teeth in this lineage. CONCLUSION: Our results validate the hypothesis that the only essential functions of ODAM that are maintained by natural selection are related to tooth development and/or the maintenance of a healthy junctional epithelium that attaches to the enamel surface of teeth.