Interordinal gene capture, the phylogenetic position of Steller's sea cow based on molecular and morphological data, and the macroevolutionary history of Sirenia

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.

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

From teeth to pad: tooth loss and development of keratinous structures in sirenians

Sirenians are a well-known example of morphological adaptation to a shallow-water grazing diet characterized by a modified feeding apparatus and orofacial morphology. Such adaptations were accompanied by an anterior tooth reduction associated with the development of keratinized pads, the evolution of which remains elusive. Among sirenians, the recently extinct Steller's sea cow represents a special case for being completely toothless. Here, we used μ-CT scans of sirenian crania to understand how motor-sensor systems associated with tooth innervation responded to innovations such as keratinized pads and continuous dental replacement. In addition, we surveyed nine genes associated with dental reduction for signatures of loss of function. Our results reveal how patterns of innervation changed with modifications of the dental formula, especially continuous replacement in manatees. Both our morphological and genomic data show that dental development was not completely lost in the edentulous Steller's sea cows. By tracing the phylogenetic history of tooth innervation, we illustrate the role of development in promoting the innervation of keratinized pads, similar to the secondary use of dental canals for innervating neomorphic keratinized structures in other tetrapod groups.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

Attenuated evolution of mammals through the Cenozoic

The Cenozoic diversification of placental mammals is the archetypal adaptive radiation. Yet, discrepancies between molecular divergence estimates and the fossil record fuel ongoing debate around the timing, tempo, and drivers of this radiation. Analysis of a three-dimensional skull dataset for living and extinct placental mammals demonstrates that evolutionary rates peak early and attenuate quickly. This long-term decline in tempo is punctuated by bursts of innovation that decreased in amplitude over the past 66 million years. Social, precocial, aquatic, and herbivorous species evolve fastest, especially whales, elephants, sirenians, and extinct ungulates. Slow rates in rodents and bats indicate dissociation of taxonomic and morphological diversification. Frustratingly, highly similar ancestral shape estimates for placental mammal superorders suggest that their earliest representatives may continue to elude unequivocal identification.

Total evidence time-scaled phylogenetic and biogeographic models for the evolution of sea cows (Sirenia, Afrotheria)

Molecular phylogenetic studies that have included sirenians from the genera Trichechus, Dugong, and Hydrodamalis have resolved their interrelationships but have yielded divergence age estimates that are problematically discordant. The ages of these lineage splits have profound implications for how to interpret the sirenian fossil record-including clade membership, biogeographic patterns, and correlations with Earth history events. In an effort to address these issues, here we present a total evidence phylogenetic analysis of Sirenia that includes living and fossil species and applies Bayesian tip-dating methods to estimate their interrelationships and divergence times. In addition to extant sirenians, our dataset includes 56 fossil species from 106 dated localities and numerous afrotherian outgroup taxa. Genetic, morphological, temporal, and biogeographic data are assessed simultaneously to bring all available evidence to bear on sirenian phylogeny. The resulting time-tree is then used for Bayesian geocoordinates reconstruction analysis, which models ancestral geographic areas at splits throughout the phylogeny, thereby allowing us to infer the direction and timing of dispersals. Our results suggest that Pan-Sirenia arose in North Africa during the latest Paleocene and that the Eocene evolution of stem sirenians was primarily situated in the Tethyan realm. In the late Eocene, some lineages moved into more northern European latitudes, an area that became the source region for a key trans-Atlantic dispersal towards the Caribbean and northern-adjacent west Atlantic. This event led to the phylogenetic and biogeographic founding of crown Sirenia with the Dugongidae-Trichechidae split occurring at the Eocene-Oligocene boundary (~33.9 Ma), temporally coincident with the onset of dropping global sea levels and temperatures. This region became the nexus of sirenian diversification and supported taxonomically-rich dugongid communities until the earliest Pliocene. The Dugonginae-Hydrodamalinae split occurred near Florida during the early Miocene (~21.2 Ma) and was followed by a west-bound dispersal that gave rise to the Pacific hydrodamalines. The late middle Miocene (~12.2 Ma) split of Dugong from all other dugongines also occurred near Florida and our analyses suggest that the Indo-Pacific distribution of modern dugongs is the result of a trans-Pacific dispersal. From at least the early Miocene, trichechid evolution was based entirely in South America, presumably within the Pebas Wetlands System. We infer that the eventual establishment of Amazon drainage into the South Atlantic allowed the dispersal of Trichechus out of South America no earlier than the mid-Pliocene. Our analyses provide a new temporal and biogeographic framework for understanding major events in sirenian evolution and their possible relationships to oceanographic and climatic changes. These hypotheses can be further tested with the recovery and integration of new fossil evidence.

Simultaneous Barcode Sequencing of Diverse Museum Collection Specimens Using a Mixed RNA Bait Set

A growing number of publications presenting results from sequencing natural history collection specimens reflect the importance of DNA sequence information from such samples. Ancient DNA extraction and library preparation methods in combination with target gene capture are a way of unlocking archival DNA, including from formalin-fixed wet-collection material. Here we report on an experiment, in which we used an RNA bait set containing baits from a wide taxonomic range of species for DNA hybridisation capture of nuclear and mitochondrial targets for analysing natural history collection specimens. The bait set used consists of 2,492 mitochondrial and 530 nuclear RNA baits and comprises specific barcode loci of diverse animal groups including both invertebrates and vertebrates. The baits allowed to capture DNA sequence information of target barcode loci from 84% of the 37 samples tested, with nuclear markers being captured more frequently and consensus sequences of these being more complete compared to mitochondrial markers. Samples from dry material had a higher rate of success than wet-collection specimens, although target sequence information could be captured from 50% of formalin-fixed samples. Our study illustrates how efforts to obtain barcode sequence information from natural history collection specimens may be combined and are a way of implementing barcoding inventories of scientific collection material.

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.

Genomic basis for skin phenotype and cold adaptation in the extinct Steller's sea cow

Steller's sea cow, an extinct sirenian and one of the largest Quaternary mammals, was described by Georg Steller in 1741 and eradicated by humans within 27 years. Here, we complement Steller's descriptions with paleogenomic data from 12 individuals. We identified convergent evolution between Steller's sea cow and cetaceans but not extant sirenians, suggesting a role of several genes in adaptation to cold aquatic (or marine) environments. Among these are inactivations of lipoxygenase genes, which in humans and mouse models cause ichthyosis, a skin disease characterized by a thick, hyperkeratotic epidermis that recapitulates Steller's sea cows' reportedly bark-like skin. We also found that Steller's sea cows' abundance was continuously declining for tens of thousands of years before their description, implying that environmental changes also contributed to their extinction.

Genomic evidence for the parallel regression of melatonin synthesis and signaling pathways in placental mammals

Background: The study of regressive evolution has yielded a wealth of examples where the underlying genes bear molecular signatures of trait degradation, such as pseudogenization or deletion. Typically, it appears that such disrupted genes are limited to the function of the regressed trait, whereas pleiotropic genes tend to be maintained by natural selection to support their myriad purposes. One such set of pleiotropic genes is involved in the synthesis ( AANAT, ASMT) and signaling ( MTNR1A, MTNR1B) of melatonin, a hormone secreted by the vertebrate pineal gland. Melatonin provides a signal of environmental darkness, thereby influencing the circadian and circannual rhythmicity of numerous physiological traits. Therefore, the complete loss of a pineal gland and the underlying melatonin pathway genes seems likely to be maladaptive, unless compensated by extrapineal sources of melatonin. Methods: We examined AANAT, ASMT, MTNR1A and MTNR1B in 123 vertebrate species, including pineal-less placental mammals and crocodylians. We searched for inactivating mutations and modelled selective pressures (dN/dS) to test whether the genes remain functionally intact. Results: We report that crocodylians retain intact melatonin genes and express AANAT and ASMT in their eyes, whereas all four genes have been repeatedly inactivated in the pineal-less xenarthrans, pangolins, sirenians, and whales. Furthermore, colugos have lost these genes, and several lineages of subterranean mammals have partial melatonin pathway dysfunction. These results are supported by the presence of shared inactivating mutations across clades and analyses of selection pressure based on the ratio of non-synonymous to synonymous substitutions (dN/dS), suggesting extended periods of relaxed selection on these genes. Conclusions: The losses of melatonin synthesis and signaling date to tens of millions of years ago in several lineages of placental mammals, raising questions about the evolutionary resilience of pleiotropic genes, and the causes and consequences of losing melatonin pathways in these species.

Successful application of ancient DNA extraction and library construction protocols to museum wet collection specimens

Millions of scientific specimens are housed in museum collections, a large part of which are fluid preserved. The use of formaldehyde as fixative and subsequent storage in ethanol is especially common in ichthyology and herpetology. This type of preservation damages DNA and reduces the chance of successful retrieval of genetic data. We applied ancient DNA extraction and single stranded library construction protocols to a variety of vertebrate samples obtained from wet collections and of different ages. Our results show that almost all samples tested yielded endogenous DNA. Archival DNA extraction was successful across different tissue types as well as using small amounts of tissue. Conversion of archival DNA fragments into single-stranded libraries resulted in usable data even for samples with initially undetectable DNA amounts. Subsequent target capture approaches for mitochondrial DNA using homemade baits on a subset of 30 samples resulted in almost complete mitochondrial genome sequences in several instances. Thus, application of ancient DNA methodology makes wet collection specimens, including type material as well as rare, old or extinct species, accessible for genetic and genomic analyses. Our results, accompanied by detailed step-by-step protocols, are a large step forward to open the DNA archive of museum wet collections for scientific studies.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

Differing effects of size and lifestyle on bone structure in mammals

BACKGROUND

Mammals are a highly diverse group, with body mass ranging from 2 g to 170 t, and encompassing species with terrestrial, aquatic, aerial, and subterranean lifestyles. The skeleton is involved in most aspects of vertebrate life history, but while previous macroevolutionary analyses have shown that structural, phylogenetic, and functional factors influence the gross morphology of skeletal elements, their inner structure has received comparatively little attention. Here we analysed bone structure of the humerus and mid-lumbar vertebrae across mammals and their correlations with different lifestyles and body size.

RESULTS

We acquired bone structure parameters in appendicular and axial elements (humerus and mid-lumbar vertebra) from 190 species across therian mammals (placentals + marsupials). Our sample captures all transitions to aerial, fully aquatic, and subterranean lifestyles in extant therian clades. We found that mammalian bone structure is highly disparate and we show that the investigated vertebral structure parameters mostly correlate with body size, but not lifestyle, while the opposite is true for humeral parameters. The latter also show a high degree of convergence among the clades that have acquired specialised (non-terrestrial) lifestyles.

CONCLUSIONS

In light of phylogenetic, size, and functional factors, the distribution of each investigated structural parameter reveals patterns explaining the construction of appendicular and axial skeletal elements in mammalian species spanning most of the extant diversity of the clade in terms of body size and lifestyle. These patterns should be further investigated with analyses focused on specific lifestyle transitions that would ideally include key fossils.

Steller's sea cow genome suggests this species began going extinct before the arrival of Paleolithic humans

Anthropogenic activity is the top factor directly related to the extinction of several animal species. The last Steller's sea cow (Hydrodamalis gigas) population on the Commander Islands (Russia) was wiped out in the second half of the 18^th century due to sailors and fur traders hunting it for the meat and fat. However, new data suggests that the extinction process of this species began much earlier. Here, we present a nuclear de novo assembled genome of H. gigas with a 25.4× depth coverage. Our results demonstrate that the heterozygosity of the last population of this animal is low and comparable to the last woolly mammoth population that inhabited Wrangel Island 4000 years ago. Besides, as a matter of consideration, our findings also demonstrate that the extinction of this marine mammal starts along the North Pacific coastal line much earlier than the first Paleolithic humans arrived in the Bering sea region.

Redescription of the skull of the Australian flatback sea turtle, Natator depressus, provides new morphological evidence for phylogenetic relationships among sea turtles (Chelonioidea)

Chelonioidea (sea turtles) are a group where available morphological evidence for crown-group relationships are incongruent with those established using molecular data. However, morphological surveys of crown-group taxa tend to focus on a recurring subset of the extant species. The Australian flatback sea turtle, Natator depressus, is often excluded from comparisons and it is the most poorly known of the seven extant species of Chelonioidea. Previous descriptions of its skull morphology are limited and conflict. Here we describe three skulls of adult N. depressus and re-examine the phylogenetic relationships according to morphological character data. Using X-ray micro Computed Tomography we describe internal structures of the braincase and identify new phylogenetically informative characters not previously reported. Phylogenetic analysis using a Bayesian approach strongly supports a sister-group relationship between Chelonia mydas and N. depressus, a topology that was not supported by previous analyses of morphological data but one that matches the topology supported by analysis of molecular data. Our results highlight the general need to sample the morphological anatomy of crown-group taxa more thoroughly before concluding that morphological and molecular evidence are incongruous.

Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation

Big, time-scaled phylogenies are fundamental to connecting evolutionary processes to modern biodiversity patterns. Yet inferring reliable phylogenetic trees for thousands of species involves numerous trade-offs that have limited their utility to comparative biologists. To establish a robust evolutionary timescale for all approximately 6,000 living species of mammals, we developed credible sets of trees that capture root-to-tip uncertainty in topology and divergence times. Our "backbone-and-patch" approach to tree building applies a newly assembled 31-gene supermatrix to two levels of Bayesian inference: (1) backbone relationships and ages among major lineages, using fossil node or tip dating, and (2) species-level "patch" phylogenies with nonoverlapping in-groups that each correspond to one representative lineage in the backbone. Species unsampled for DNA are either excluded ("DNA-only" trees) or imputed within taxonomic constraints using branch lengths drawn from local birth-death models ("completed" trees). Joining time-scaled patches to backbones results in species-level trees of extant Mammalia with all branches estimated under the same modeling framework, thereby facilitating rate comparisons among lineages as disparate as marsupials and placentals. We compare our phylogenetic trees to previous estimates of mammal-wide phylogeny and divergence times, finding that (1) node ages are broadly concordant among studies, and (2) recent (tip-level) rates of speciation are estimated more accurately in our study than in previous "supertree" approaches, in which unresolved nodes led to branch-length artifacts. Credible sets of mammalian phylogenetic history are now available for download at http://vertlife.org/phylosubsets, enabling investigations of long-standing questions in comparative biology.

Phylogeography of the dugong (Dugong dugon) based on historical samples identifies vulnerable Indian Ocean populations

We investigated the phylogeography of the dugong (Dugong dugon) across its original range using museum material from 14 natural history museum and university collections. The mitochondrial DNA control region was successfully amplified from samples of bone or tooth powder from 162 individuals. These samples range from 1827 to 1996 and span the historical distribution range of the dugong. We were able to successfully amplify overlapping fragments of the D-loop region of the mitochondrial DNA (mtDNA) resulting in sequences of a 355 bp fragment for 162 individuals for the final analyses. This included a new sequence (189 bp) from a previously unidentified piece of skin of the extinct Steller’s sea cow (Hydrodamalis gigas), as an outgroup. The resulting dugong sequences match those from previous studies of dugongs from Australia and Indonesia, but revealed several new and divergent mtDNA lineages in the Indian Ocean. One mtDNA lineage includes most specimens from the Western Indian Ocean, with another distinct lineage isolated to nearby Madagascar and Comores. There is little geographic structuring detectable among other populations in the Western Indian Ocean and all populations from that region appear to have historically contained comparatively low levels of genetic diversity. The genetic diversity of several Indian Ocean samples collected after 1950 was lower than that of the samples collected earlier from similar locations, a result coincident with the anecdotal reductions in population size. The new lineages and potential loss of diversity highlight the particular conservation importance and vulnerability of dugong populations in the Western Indian Ocean.

Emergence of a Chimeric Globin Pseudogene and Increased Hemoglobin Oxygen Affinity Underlie the Evolution of Aquatic Specializations in Sirenia

As limits on O₂ availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller’s sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb–O₂ affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O₂ extraction from finite lung stores and suppress tissue O₂ offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb–O₂ affinity in (sub)Arctic Steller’s sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.

The Marine Mammal Class II Major Histocompatibility Complex Organization

Sirenians share with cetaceans and pinnipeds several convergent traits selected for the aquatic lifestyle. Living in water poses new challenges not only for locomotion and feeding but also for combating new pathogens, which may render the immune system one of the best tools aquatic mammals have for dealing with aquatic microbial threats. So far, only cetaceans have had their class II Major Histocompatibility Complex (MHC) organization characterized, despite the importance of MHC genes for adaptive immune responses. This study aims to characterize the organization of the marine mammal class II MHC using publicly available genomes. We located class II sequences in the genomes of one sirenian, four pinnipeds and eight cetaceans using NCBI-BLAST and reannotated the sequences using local BLAST search with exon and intron libraries. Scaffolds containing class II sequences were compared using dotplot analysis and introns were used for phylogenetic analysis. The manatee class II region shares overall synteny with other mammals, however most DR loci were translocated from the canonical location, past the extended class II region. Detailed analysis of the genomes of closely related taxa revealed that this presumed translocation is shared with all other living afrotherians. Other presumptive chromosome rearrangements in Afrotheria are the deletion of DQ loci in Afrosoricida and deletion of DP in E. telfairi. Pinnipeds share the main features of dog MHC: lack of a functional pair of DPA/DPB genes and inverted DRB locus between DQ and DO subregions. All cetaceans share the Cetartiodactyla inversion separating class II genes into two subregions: class IIa, with DR and DQ genes, and class IIb, with non-classic genes and a DRB pseudogene. These results point to three distinct and unheralded class II MHC structures in marine mammals: one canonical organization but lacking DP genes in pinnipeds; one bearing an inversion separating IIa and IIb subregions lacking DP genes found in cetaceans; and one with a translocation separating the most diverse class II gene from the MHC found in afrotherians and presumptive functional DR, DQ, and DP genes. Future functional research will reveal how these aquatic mammals cope with pathogen pressures with these divergent MHC organizations.

Evolutionary adaptation to aquatic lifestyle in extinct sloths can lead to systemic alteration of bone structure

Through phenotypic plasticity, bones can change in structure and morphology, in response to physiological and biomechanical influences over the course of individual life. Changes in bones also occur in evolution as functional adaptations to the environment. In this study, we report on the evolution of bone mass increase (BMI) that occurred in the postcranium and skull of extinct aquatic sloths. Although non-pathological BMI in postcranial skeleton has been known in aquatic mammals, we here document general BMI in the skull for the first time. We present evidence of thickening of the nasal turbinates, nasal septum and cribriform plate, further thickening of the frontals, and infilling of sinus spaces by compact bone in the late and more aquatic species of the extinct sloth Thalassocnus Systemic bone mass increase occurred among the successively more aquatic species of Thalassocnus, as an evolutionary adaptation to the lineage's changing environment. The newly documented pachyostotic turbinates appear to have conferred little or no functional advantage and are here hypothesized as a correlation with or consequence of the systemic BMI among Thalassocnus species. This could, in turn, be consistent with a genetic accommodation of a physiological adjustment to a change of environment.

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.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1359-6) contains supplementary material, which is available to authorized users.

Application of Solid-State Capture for the Retrieval of Small-to-Medium Sized Target Loci from Ancient DNA

Genetic studies that include ancient samples are often hampered by the low amount of endogenous DNA that ancient samples often contain, relative to co-extracted "contaminant" DNA from other organisms. One approach to mitigate this challenge is to perform hybridization-based capture of target genomic regions using DNA or RNA baits. Such baits are designed to have high sequence similarity to the target genomic regions and can reduce the off-target fraction in DNA sequencing libraries. Here, we present a protocol to use Agilent SureSelect microarrays to enrich ancient DNA libraries for small-to-medium-sized target loci, such as mitochondrial genomes, from ancient DNA extracts. The protocol that we present builds on previously published work by introducing improvements that improve recovery of short DNA fragments while minimizing the cost and duration of the experiment.

A phylogenomic approach to reconstruct interrelationships of main clupeocephalan lineages with a critical discussion of morphological apomorphies

BACKGROUND

Previous molecular studies on the phylogeny and classification of clupeocephalan fishes revealed numerous new taxonomic entities. For re-analysing these taxa, we perform target gene capturing and subsequent next generation sequencing of putative ortholog exons of major clupeocephalan lineages. Sequence information for the RNA bait design was derived from publicly available genomes of bony fishes. Newly acquired sequence data comprising > 800 exon sequences was subsequently used for phylogenetic reconstructions.

RESULTS

Our results support monophyletic Otomorpha comprising Alepocephaliformes. Within Ostariophysi, Gonorynchiformes are sister to a clade comprising Cypriniformes, Characiformes, Siluriformes and Gymnotiformes, where the interrelationships of Characiformes, Siluriformes and Gymnotiformes remain enigmatic. Euteleosts comprise four major clades: Lepidogalaxiiformes, Protacanthopterygii, Stomiatii, and Galaxiiformes plus Neoteleostei. The monotypic Lepidogalaxiiformes form the sister-group to all remaining euteleosts. Protacanthopterygii, comprising Argentini-, Esoci- and Salmoniformes, is sister to Stomiatii (Osmeriformes and Stomiatiformes) and Galaxiiformes plus Neoteleostei.

CONCLUSIONS

Several proposed monophyla defined by morphological apomorphies within the Clupeocephalan phylogeny are confirmed by the phylogenetic estimates presented herein. However, other morphologically described groups cannot be reconciled with molecular phylogenies. Thus, numerous morphological apomoprhies of supposed monophyla are called into question. The interpretation of suggested morphological synapomorphies of otomorph fishes is strongly affected by the inclusion of deep-sea inhabiting, and to that effect morphologically adapted Alepocephaliformes. Our revision of these potential synapomorphies, in the context that Alepocephaliformes are otomorph fishes, reveals that only a single character of the total nine characters proposed as synapomorphic for the group is clearly valid for all otomorphs. Three further characters remain possible apomorphies since their status remains unclear in the deep-sea adapted Alepocephaliformes showing developmental lag and lacking a swim bladder. Further, our analysis places Galaxiiformes as sister group to neoteleosts, which contradicts some previous molecular phylogenetic studies. This needs further investigation from a morphological perspective, as suggested synapomophies for several euteleostean lineages are challenged or still lacking. For the verification of results presented herein, a denser phylogenomic-level taxon sampling should be applied.

Molecular phylogeny of the extinct Steller's sea cow and other Sirenia species based on their complete mitochondrial genomes

The Steller's sea cow - Hydrodamalis gigas (Dugongidae: Sirenia) - is an extinct herbivorous marine mammal which inhabited the North Pacific Ocean during the Pleistocene and Holocene. H. gigas was the largest member of the Sirenia order and disappeared in the middle of the 18th century. Here, we present the complete sequence of the mitochondrial genome of this extinct animal. The Steller's sea cow mitochondrial DNA (mtDNA) is 16,872 base pairs (bp) in length and contains a set of mitochondrial genes typical for mammals. Phylogenetic analysis based on complete mitochondrial genomes of the sirenian species allows accurate assessment of the degree of their mitogenomic diversification during millions of years of evolution.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Evolutionary adaptation to aquatic lifestyle in extinct sloths can lead to systemic alteration of bone structure

Through phenotypic plasticity, bones can change in structure and morphology, in response to physiological and biomechanical influences over the course of individual life. Changes in bones also occur in evolution as functional adaptations to the environment. In this study, we report on the evolution of bone mass increase (BMI) that occurred in the postcranium and skull of extinct aquatic sloths. Although non-pathological BMI in postcranial skeleton has been known in aquatic mammals, we here document general BMI in the skull for the first time. We present evidence of thickening of the nasal turbinates, nasal septum and cribriform plate, further thickening of the frontals, and infilling of sinus spaces by compact bone in the late and more aquatic species of the extinct sloth Thalassocnus Systemic bone mass increase occurred among the successively more aquatic species of Thalassocnus, as an evolutionary adaptation to the lineage's changing environment. The newly documented pachyostotic turbinates appear to have conferred little or no functional advantage and are here hypothesized as a correlation with or consequence of the systemic BMI among Thalassocnus species. This could, in turn, be consistent with a genetic accommodation of a physiological adjustment to a change of environment.

Evolutionary adaptation to aquatic lifestyle in extinct sloths can lead to systemic alteration of bone structure

Through phenotypic plasticity, bones can change in structure and morphology, in response to physiological and biomechanical influences over the course of individual life. Changes in bones also occur in evolution as functional adaptations to the environment. In this study, we report on the evolution of bone mass increase (BMI) that occurred in the postcranium and skull of extinct aquatic sloths. Although non-pathological BMI in postcranial skeleton has been known in aquatic mammals, we here document general BMI in the skull for the first time. We present evidence of thickening of the nasal turbinates, nasal septum and cribriform plate, further thickening of the frontals, and infilling of sinus spaces by compact bone in the late and more aquatic species of the extinct sloth Thalassocnus Systemic bone mass increase occurred among the successively more aquatic species of Thalassocnus, as an evolutionary adaptation to the lineage's changing environment. The newly documented pachyostotic turbinates appear to have conferred little or no functional advantage and are here hypothesized as a correlation with or consequence of the systemic BMI among Thalassocnus species. This could, in turn, be consistent with a genetic accommodation of a physiological adjustment to a change of environment.

First adequately-known quadrupedal sirenian from Eurasia (Eocene, Bay of Biscay, Huesca, northeastern Spain)

Sirenians are the only extant herbivorous mammals fully adapted to an aquatic lifestyle. They originated in Africa during the Paleocene from an undetermined clade of afrotherian mammals, and by the end of the Eocene they were widely distributed across the tropical latitudes. Here we introduce Sobrarbesiren cardieli gen. et sp. nov. It is the first adequately-known quadrupedal sirenian from Eurasia and the oldest record of this clade from western Europe. Fossils have been recovered from the middle Lutetian (SBZ15) site of Castejón de Sobrarbe-41 (Huesca, Spain), and comprise many cranial and postcranial remains, including pelvic girdle and hind limb bones, from at least six sirenian individuals of different ontogenetic stages. Sobrarbesiren shows a suite of characters previously considered synapomorphies of different clades of derived sirenians, such as the presence of the processus retroversus of the squamosal and the pterygoid fossa, combined with ancestral characters such as the presence of an alisphenoid canal, a permanent P5, at least two sacral vertebrae, a primitive pelvis and functional femora and fibulae. Sobrarbesiren is recovered as the sister taxon of Dugongidae and represents a transitional stage of adaptation to aquatic life between the amphibious quadrupedal prorastomids and the aquatic quadrupedal protosirenids.

Their loss is our gain: regressive evolution in vertebrates provides genomic models for uncovering human disease loci

Throughout Earth's history, evolution's numerous natural 'experiments' have resulted in a diverse range of phenotypes. Though de novo phenotypes receive widespread attention, degeneration of traits inherited from an ancestor is a very common, yet frequently neglected, evolutionary path. The latter phenomenon, known as regressive evolution, often results in vertebrates with phenotypes that mimic inherited disease states in humans. Regressive evolution of anatomical and/or physiological traits is typically accompanied by inactivating mutations underlying these traits, which frequently occur at loci identical to those implicated in human diseases. Here we discuss the potential utility of examining the genomes of vertebrates that have experienced regressive evolution to inform human medical genetics. This approach is low cost and high throughput, giving it the potential to rapidly improve knowledge of disease genetics. We discuss two well-described examples, rod monochromacy (congenital achromatopsia) and amelogenesis imperfecta, to demonstrate the utility of this approach, and then suggest methods to equip non-experts with the ability to corroborate candidate genes and uncover new disease loci.

Evolution of UCP1 Transcriptional Regulatory Elements Across the Mammalian Phylogeny

Uncoupling protein 1 (UCP1) permits non-shivering thermogenesis (NST) when highly expressed in brown adipose tissue (BAT) mitochondria. Exclusive to placental mammals, BAT has commonly been regarded to be advantageous for thermoregulation in hibernators, small-bodied species, and the neonates of larger species. While numerous regulatory control motifs associated with UCP1 transcription have been proposed for murid rodents, it remains unclear whether these are conserved across the eutherian mammal phylogeny and hence essential for UCP1 expression. To address this shortcoming, we conducted a broad comparative survey of putative UCP1 transcriptional regulatory elements in 139 mammals (135 eutherians). We find no evidence for presence of a UCP1 enhancer in monotremes and marsupials, supporting the hypothesis that this control region evolved in a stem eutherian ancestor. We additionally reveal that several putative promoter elements (e.g., CRE-4, CCAAT) identified in murid rodents are not conserved among BAT-expressing eutherians, and together with the putative regulatory region (PRR) and CpG island do not appear to be crucial for UCP1 expression. The specificity and importance of the upTRE, dnTRE, URE1, CRE-2, RARE-2, NBRE, BRE-1, and BRE-2 enhancer elements first described from rats and mice are moreover uncertain as these motifs differ substantially—but generally remain highly conserved—in other BAT-expressing eutherians. Other UCP1 enhancer motifs (CRE-3, PPRE, and RARE-3) as well as the TATA box are also highly conserved in nearly all eutherian lineages with an intact UCP1. While these transcriptional regulatory motifs are generally also maintained in species where this gene is pseudogenized, the loss or degeneration of key basal promoter (e.g., TATA box) and enhancer elements in other UCP1-lacking lineages make it unlikely that the enhancer region is pleiotropic (i.e., co-regulates additional genes). Importantly, differential losses of (or mutations within) putative regulatory elements among the eutherian lineages with an intact UCP1 suggests that the transcriptional control of gene expression is not highly conserved in this mammalian clade.

Inactivation of thermogenic UCP1 as a historical contingency in multiple placental mammal clades

Mitochondrial uncoupling protein 1 (UCP1) is essential for nonshivering thermogenesis in brown adipose tissue and is widely accepted to have played a key thermoregulatory role in small-bodied and neonatal placental mammals that enabled the exploitation of cold environments. We map ucp1 sequences from 133 mammals onto a species tree constructed from a ~51-kb sequence alignment and show that inactivating mutations have occurred in at least 8 of the 18 traditional placental orders, thereby challenging the physiological importance of UCP1 across Placentalia. Selection and timetree analyses further reveal that ucp1 inactivations temporally correspond with strong secondary reductions in metabolic intensity in xenarthrans and pangolins, or in six other lineages coincided with a ~30 million-year episode of global cooling in the Paleogene that promoted sharp increases in body mass and cladogenesis evident in the fossil record. Our findings also demonstrate that members of various lineages (for example, cetaceans, horses, woolly mammoths, Steller's sea cows) evolved extreme cold hardiness in the absence of UCP1-mediated thermogenesis. Finally, we identify ucp1 inactivation as a historical contingency that is linked to the current low species diversity of clades lacking functional UCP1, thus providing the first evidence for species selection related to the presence or absence of a single gene product.

Waking the undead: Implications of a soft explosive model for the timing of placental mammal diversification

The explosive, long fuse, and short fuse models represent competing hypotheses for the timing of placental mammal diversification. Support for the explosive model, which posits both interordinal and intraordinal diversification after the KPg mass extinction, derives from morphological cladistic studies that place Cretaceous eutherians outside of crown Placentalia. By contrast, most molecular studies favor the long fuse model wherein interordinal cladogenesis occurred in the Cretaceous followed by intraordinal cladogenesis after the KPg boundary. Phillips (2016) proposed a soft explosive model that allows for the emergence of a few lineages (Xenarthra, Afrotheria, Euarchontoglires, Laurasiatheria) in the Cretaceous, but otherwise agrees with the explosive model in positing the majority of interordinal diversification after the KPg mass extinction. Phillips (2016) argues that rate transference errors associated with large body size and long lifespan have inflated previous estimates of interordinal divergence times, and further suggests that most interordinal divergences are positioned after the KPg boundary when rate transference errors are avoided through the elimination of calibrations in large-bodied and/or long lifespan clades. Here, we show that rate transference errors can also occur in the opposite direction and drag forward estimated divergence dates when calibrations in large-bodied/long lifespan clades are omitted. This dragging forward effect results in the occurrence of more than half a billion years of 'zombie lineages' on Phillips' preferred timetree. By contrast with ghost lineages, which are a logical byproduct of an incomplete fossil record, zombie lineages occur when estimated divergence dates are younger than the minimum age of the oldest crown fossils. We also present the results of new timetree analyses that address the rate transference problem highlighted by Phillips (2016) by deleting taxa that exceed thresholds for body size and lifespan. These analyses recover all interordinal divergence times in the Cretaceous and are consistent with the long fuse model of placental diversification. Finally, we outline potential problems with morphological cladistic analyses of higher-level relationships among placental mammals that may account for the perceived discrepancies between molecular and paleontological estimates of placental divergence times.

Barcoding the largest animals on Earth: ongoing challenges and molecular solutions in the taxonomic identification of ancient cetaceans

Over the last few centuries, many cetacean species have witnessed dramatic global declines due to industrial overharvesting and other anthropogenic influences, and thus are key targets for conservation. Whale bones recovered from archaeological and palaeontological contexts can provide essential baseline information on the past geographical distribution and abundance of species required for developing informed conservation policies. Here we review the challenges with identifying whale bones through traditional anatomical methods, as well as the opportunities provided by new molecular analyses. Through a case study focused on the North Sea, we demonstrate how the utility of this (pre)historic data is currently limited by a lack of accurate taxonomic information for the majority of ancient cetacean remains. We then discuss current opportunities presented by molecular identification methods such as DNA barcoding and collagen peptide mass fingerprinting (zooarchaeology by mass spectrometry), and highlight the importance of molecular identifications in assessing ancient species' distributions through a case study focused on the Mediterranean. We conclude by considering high-throughput molecular approaches such as hybridization capture followed by next-generation sequencing as cost-effective approaches for enhancing the ecological informativeness of these ancient sample sets.This article is part of the themed issue 'From DNA barcodes to biomes'.

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.

Inactivation of C4orf26 in toothless placental mammals

Previous studies have reported inactivated copies of six enamel-related genes (AMBN, AMEL, AMTN, ENAM, KLK4, MMP20) and one dentin-related gene (DSPP) in one or more toothless vertebrates and/or vertebrates with enamelless teeth, thereby providing evidence that these genes are enamel or tooth-specific with respect to their critical functions that are maintained by natural selection. Here, we employ available genome sequences for edentulous and enamelless mammals to evaluate the enamel specificity of four genes (WDR72, SLC24A4, FAM83H, C4orf26) that have been implicated in amelogenesis imperfecta, a condition in which proper enamel formation is abrogated during tooth development. Coding sequences for WDR72, SCL24A4, and FAM83H are intact in four edentulous taxa (Chinese pangolin, three baleen whales) and three taxa (aardvark, nine-banded armadillo, Hoffmann's two-toed sloth) with enamelless teeth, suggesting that these genes have critical functions beyond their involvement in tooth development. By contrast, genomic data for C4orf26 reveal inactivating mutations in pangolin and bowhead whale as well as evidence for deletion of this gene in two minke whale species. Hybridization capture of exonic regions and PCR screens provide evidence for inactivation of C4orf26 in eight additional baleen whale species. However, C4orf26 is intact in all three species with enamelless teeth that were surveyed, as well as in 95 additional mammalian species with enamel-capped teeth. Estimates of selection intensity suggest that dN/dS ratios on branches leading to taxa with enamelless teeth are similar to the dN/dS ratio on branches leading to taxa with enamel-capped teeth. Based on these results, we conclude that C4orf26 is tooth-specific, but not enamel-specific, with respect to its essential functions that are maintained by natural selection. A caveat is that an alternative splice site variant, which translates exon 3 in a different reading frame, is putatively functional in Catarrhini and may have evolved an additional role in this primate clade.

Resurrecting phenotypes from ancient DNA sequences: promises and perspectives

Anatomical changes in extinct mammalian lineages over evolutionary time, such as the loss of fingers and teeth and the rapid increase in body size that accompanied the late Miocene dispersal of the progenitors of Steller’s sea cows (Hydrodamalis gigas (Zimmermann, 1780)) into North Pacific waters and the convergent development of a thick pelage and accompanying reductions in ear and tail surface area of woolly mammoths (Mammuthus primigenius (Blumenbach, 1799)) and woolly rhinoceros (Coelodonta antiquitatis (Blumenbach, 1799)), are prime examples of adaptive evolution underlying the exploitation of new habitats. It is likely, however, that biochemical specializations adopted during these evolutionary transitions were of similar or even greater biological importance. As these “living” processes do not fossilize, direct information regarding the physiological attributes of extinct species has largely remained beyond the range of scientific inquiry. However, the ability to retrieve genomic sequences from ancient DNA samples, combined with ectopic expression systems, now permit the evolutionary origins and structural and functional properties of authentic prehistoric proteins to be examined in great detail. Exponential technical advances in ancient DNA retrieval, enrichment, and sequencing will soon permit targeted generation of complete genomes from hundreds of extinct species across the last one million years that, in combination with emerging in vitro expression, genome engineering, and cell differentiation techniques, promises to herald an exciting new trajectory of evolutionary research at the interface of biochemistry, genomics, palaeontology, and cell biology.