Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae).

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].

Positive Selection and Inactivation in the Vision and Hearing Genes of Cetaceans.

The transition to an aquatic lifestyle in cetaceans (whales and dolphins) resulted in a radical transformation in their sensory systems. Toothed whales acquired specialized high-frequency hearing tied to the evolution of echolocation, whereas baleen whales evolved low-frequency hearing. More generally, all cetaceans show adaptations for hearing and seeing underwater. To determine the extent to which these phenotypic changes have been driven by molecular adaptation, we performed large-scale targeted sequence capture of 179 sensory genes across the Cetacea, incorporating up to 54 cetacean species from all major clades as well as their closest relatives, the hippopotamuses. We screened for positive selection in 167 loci related to vision and hearing and found that the diversification of cetaceans has been accompanied by pervasive molecular adaptations in both sets of genes, including several loci implicated in nonsyndromic hearing loss. Despite these findings, however, we found no direct evidence of positive selection at the base of odontocetes coinciding with the origin of echolocation, as found in studies examining fewer taxa. By using contingency tables incorporating taxon- and gene-based controls, we show that, although numbers of positively selected hearing and nonsyndromic hearing loss genes are disproportionately high in cetaceans, counts of vision genes do not differ significantly from expected values. Alongside these adaptive changes, we find increased evidence of pseudogenization of genes involved in cone-mediated vision in mysticetes and deep-diving odontocetes.

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.

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 $\sim $3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from $>$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.].

Ancestral resurrection of anthropoid estrogen receptor ß demonstrates functional consequences of positive selection.

Anthropoid primates arose during the Eocene approximately 55 million years ago (mya), and extant anthropoids share a most recent common ancestor ∼40mya. Paleontology has been very successful at describing the morphological phenotypes of extinct anthropoids. Less well understood is the molecular biology of these extinct species as well as the phenotypic consequences of evolutionary variation in their genomes. Here we resurrect the most recent common ancestral anthropoid estrogen receptor ß gene (ESR2) and demonstrate that the function of this ancestral estrogen receptor has been maintained during human descent but was altered during early New World monkey (NWM) evolution by becoming a more potent transcriptional activator. We tested hypotheses of adaptive evolution in the protein coding sequences of ESR2, and determined that ESR2 evolved via episodic positive selection on the NWM stem lineage. We separately co-transfected ESR2 constructs for human, NWM, and the anthropoid ancestor along with reporter gene vectors and performed hormone binding dose response experiments that measure transactivation activity. We found the transactivation potentials of the ancestral and human sequences to be significantly lower (p<0.0001 in each comparison) than that of the NWM when treated with estradiol, the most prevalent estrogen. We conclude the difference in fold activation is due to positive selection in the NWM ERß ligand binding domain. Our study validates inferential methods for detecting adaptive evolution that predict functional consequences of nucleotide substitutions and points a way toward examining the functional consequences of positive Darwinian selection.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals.

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

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.

Drivers of morphological evolution in the toothed whale jaw.

Toothed whales (odontocetes) emit high-frequency underwater sounds (echolocate)-an extreme and unique innovation allowing them to sense their prey and environment. Their highly specialized mandible (lower jaw) allows high-frequency sounds to be transmitted back to the inner ear. Echolocation is evident in the earliest toothed whales, but little research has focused on the evolution of mandibular form regarding this unique adaptation. Here, we use a high-density, three-dimensional geometric morphometric analysis of 100 living and extinct cetacean species spanning their ∼50-million-year evolutionary history. Our analyses demonstrate that most shape variation is found in the relative length of the jaw and the mandibular symphysis. The greatest morphological diversity was obtained during two periods of rapid evolution: the initial evolution of archaeocetes (stem whales) in the early to mid-Eocene as they adapted to an aquatic lifestyle, representing one of the most extreme adaptive transitions known, and later on in the mid-Oligocene odontocetes as they became increasingly specialized for a range of diets facilitated by increasingly refined echolocation. Low disparity in the posterior mandible suggests the shape of the acoustic window, which receives sound, has remained conservative since the advent of directional hearing in the aquatic archaeocetes, even as the earliest odontocetes began to receive sounds from echolocation. Diet, echolocation, feeding method, and dentition type strongly influence mandible shape. Unlike in the toothed whale cranium, we found no significant asymmetry in the mandible. We suggest that a combination of refined echolocation and associated dietary specializations have driven morphology and disparity in the toothed whale mandible.

Relaxed clocks and inferences of heterogeneous patterns of nucleotide substitution and divergence time estimates across whales and dolphins (Mammalia: Cetacea).

Various nucleotide substitution models have been developed to accommodate among lineage rate heterogeneity, thereby relaxing the assumptions of the strict molecular clock. Recently developed "uncorrelated relaxed clock" and "random local clock" (RLC) models allow decoupling of nucleotide substitution rates between descendant lineages and are thus predicted to perform better in the presence of lineage-specific rate heterogeneity. However, it is uncertain how these models perform in the presence of punctuated shifts in substitution rate, especially between closely related clades. Using cetaceans (whales and dolphins) as a case study, we test the performance of these two substitution models in estimating both molecular rates and divergence times in the presence of substantial lineage-specific rate heterogeneity. Our RLC analyses of whole mitochondrial genome alignments find evidence for up to ten clade-specific nucleotide substitution rate shifts in cetaceans. We provide evidence that in the uncorrelated relaxed clock framework, a punctuated shift in the rate of molecular evolution within a subclade results in posterior rate estimates that are either misled or intermediate between the disparate rate classes present in baleen and toothed whales. Using simulations, we demonstrate abrupt changes in rate isolated to one or a few lineages in the phylogeny can mislead rate and age estimation, even when the node of interest is calibrated. We further demonstrate how increasing prior age uncertainty can bias rate and age estimates, even while the 95% highest posterior density around age estimates decreases; in other words, increased precision for an inaccurate estimate. We interpret the use of external calibrations in divergence time studies in light of these results, suggesting that rate shifts at deep time scales may mislead inferences of absolute molecular rates and ages.

A phylogenetic blueprint for a modern whale.

The emergence of Cetacea in the Paleogene represents one of the most profound macroevolutionary transitions within Mammalia. The move from a terrestrial habitat to a committed aquatic lifestyle engendered wholesale changes in anatomy, physiology, and behavior. The results of this remarkable transformation are extant whales that include the largest, biggest brained, fastest swimming, loudest, deepest diving mammals, some of which can detect prey with a sophisticated echolocation system (Odontoceti - toothed whales), and others that batch feed using racks of baleen (Mysticeti - baleen whales). A broad-scale reconstruction of the evolutionary remodeling that culminated in extant cetaceans has not yet been based on integration of genomic and paleontological information. Here, we first place Cetacea relative to extant mammalian diversity, and assess the distribution of support among molecular datasets for relationships within Artiodactyla (even-toed ungulates, including Cetacea). We then merge trees derived from three large concatenations of molecular and fossil data to yield a composite hypothesis that encompasses many critical events in the evolutionary history of Cetacea. By combining diverse evidence, we infer a phylogenetic blueprint that outlines the stepwise evolutionary development of modern whales. This hypothesis represents a starting point for more detailed, comprehensive phylogenetic reconstructions in the future, and also highlights the synergistic interaction between modern (genomic) and traditional (morphological+paleontological) approaches that ultimately must be exploited to provide a rich understanding of evolutionary history across the entire tree of Life.