Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Grey wolf genomic history reveals a dual ancestry of dogs.

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived1-8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000-30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny.

The sequencing of ancient DNA has enabled the reconstruction of speciation, migration and admixture events for extinct taxa1. However, the irreversible post-mortem degradation2 of ancient DNA has so far limited its recovery-outside permafrost areas-to specimens that are not older than approximately 0.5 million years (Myr)3. By contrast, tandem mass spectrometry has enabled the sequencing of approximately 1.5-Myr-old collagen type I4, and suggested the presence of protein residues in fossils of the Cretaceous period5-although with limited phylogenetic use6. In the absence of molecular evidence, the speciation of several extinct species of the Early and Middle Pleistocene epoch remains contentious. Here we address the phylogenetic relationships of the Eurasian Rhinocerotidae of the Pleistocene epoch7-9, using the proteome of dental enamel from a Stephanorhinus tooth that is approximately 1.77-Myr old, recovered from the archaeological site of Dmanisi (South Caucasus, Georgia)10. Molecular phylogenetic analyses place this Stephanorhinus as a sister group to the clade formed by the woolly rhinoceros (Coelodonta antiquitatis) and Merck's rhinoceros (Stephanorhinus kirchbergensis). We show that Coelodonta evolved from an early Stephanorhinus lineage, and that this latter genus includes at least two distinct evolutionary lines. The genus Stephanorhinus is therefore currently paraphyletic, and its systematic revision is needed. We demonstrate that sequencing the proteome of Early Pleistocene dental enamel overcomes the limitations of phylogenetic inference based on ancient collagen or DNA. Our approach also provides additional information about the sex and taxonomic assignment of other specimens from Dmanisi. Our findings reveal that proteomic investigation of ancient dental enamel-which is the hardest tissue in vertebrates11, and is highly abundant in the fossil record-can push the reconstruction of molecular evolution further back into the Early Pleistocene epoch, beyond the currently known limits of ancient DNA preservation.

Redefining the Evolutionary History of the Rock Dove, Columba livia, Using Whole Genome Sequences.

The domestic pigeon's exceptional phenotypic diversity was key in developing Darwin's Theory of Evolution and establishing the concept of artificial selection. However, unlike its domestic counterpart, its wild progenitor, the rock dove Columba livia has received considerably less attention. Therefore, questions regarding its domestication, evolution, taxonomy, and conservation status remain unresolved. We generated whole-genome sequencing data from 65 historical rock doves that represent all currently recognized subspecies and span the species' original geographic distribution. Our dataset includes 3 specimens from Darwin's collection, and the type specimens of 5 different taxa. We characterized their population structure, genomic diversity, and gene-flow patterns. Our results show the West African subspecies C. l. gymnocyclus is basal to rock doves and domestic pigeons, and suggests gene-flow between the rock dove's sister species C. rupestris, and the ancestor of rock doves after its split from West African populations. These genomes allowed us to propose a model for the evolution of the rock dove in light of the refugia theory. We propose that rock dove genetic diversity and introgression patterns derive from a history of allopatric cycles and dispersion waves during the Quaternary glacial and interglacial periods. To explore the rock dove domestication history, we combined our new dataset with available genomes from domestic pigeons. Our results point to at least 1 domestication event in the Levant that gave rise to all domestic breeds analysed in this study. Finally, we propose a species-level taxonomic arrangement to reflect the evolutionary history of the West African rock dove populations.

Population genomics of the muskox' resilience in the near absence of genetic variation.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Modern Siberian dog ancestry was shaped by several thousand years of Eurasian-wide trade and human dispersal.

Dogs have been essential to life in the Siberian Arctic for over 9,500 y, and this tight link between people and dogs continues in Siberian communities. Although Arctic Siberian groups such as the Nenets received limited gene flow from neighboring groups, archaeological evidence suggests that metallurgy and new subsistence strategies emerged in Northwest Siberia around 2,000 y ago. It is unclear if the Siberian Arctic dog population was as continuous as the people of the region or if instead admixture occurred, possibly in relation to the influx of material culture from other parts of Eurasia. To address this question, we sequenced and analyzed the genomes of 20 ancient and historical Siberian and Eurasian Steppe dogs. Our analyses indicate that while Siberian dogs were genetically homogenous between 9,500 to 7,000 y ago, later introduction of dogs from the Eurasian Steppe and Europe led to substantial admixture. This is clearly the case in the Iamal-Nenets region (Northwestern Siberia) where dogs from the Iron Age period (∼2,000 y ago) possess substantially less ancestry related to European and Steppe dogs than dogs from the medieval period (∼1,000 y ago). Combined with findings of nonlocal materials recovered from these archaeological sites, including glass beads and metal items, these results indicate that Northwest Siberian communities were connected to a larger trade network through which they acquired genetically distinctive dogs from other regions. These exchanges were part of a series of major societal changes, including the rise of large-scale reindeer pastoralism ∼800 y ago.

Colonialism in South Africa leaves a lasting legacy of reduced genetic diversity in Cape buffalo.

The iconic Cape buffalo has experienced several documented population declines in recent history. These declines have been largely attributed to the late 19th century rinderpest pandemic. However, the effect of the rinderpest pandemic on their genetic diversity remains contentious, and other factors that have potentially affected this diversity include environmental changes during the Pleistocene, range expansions and recent human activity. Motivated by this, we present analyses of whole genome sequencing data from 59 individuals from across the Cape buffalo range to assess present-day levels of genome-wide genetic diversity and what factors have influenced these levels. We found that the Cape buffalo has high average heterozygosity overall (0.40%), with the two southernmost populations having significantly lower heterozygosity levels (0.33% and 0.29%) on par with that of the domesticated water buffalo (0.29%). Interestingly, we found that these lower levels are probably due to recent inbreeding (average fraction of runs of homozygosity 23.7% and 19.9%) rather than factors further back in time during the Pleistocene. Moreover, detailed investigations of recent demographic history show that events across the past three centuries were the main drivers of the exceptional loss of genetic diversity in the southernmost populations, coincident with the onset of colonialism in the southern extreme of the Cape buffalo range. Hence, our results add to the growing body of studies suggesting that multiple recent human-mediated impacts during the colonial period caused massive losses of large mammal abundance in southern Africa.

Genomic diversity and differentiation between island and mainland populations of white-tailed eagles (Haliaeetus albicilla).

Divergence in the face of high dispersal capabilities is a documented but poorly understood phenomenon. The white-tailed eagle (Haliaeetus albicilla) has a large geographic dispersal capability and should theoretically be able to maintain genetic homogeneity across its dispersal range. However, following analysis of the genomic variation of white-tailed eagles, from both historical and contemporary samples, clear signatures of ancient biogeographic substructure across Europe and the North-East Atlantic is observed. The greatest genomic differentiation was observed between island (Greenland and Iceland) and mainland (Denmark, Norway and Estonia) populations. The two island populations share a common ancestry from a single mainland population, distinct from the other sampled mainland populations, and despite the potential for high connectivity between Iceland and Greenland they are well separated from each other and are characterized by inbreeding and little variation. Temporal differences also highlight a pattern of regional populations persisting despite the potential for admixture. All sampled populations generally showed a decline in effective population size over time, which may have been shaped by four historical events: (1) Isolation of refugia during the last glacial period 110-115,000 years ago, (2) population divergence following the colonization of the deglaciated areas ~10,000 years ago, (3) human population expansion, which led to the settlement in Iceland ~1100 years ago, and (4) human persecution and exposure to toxic pollutants during the last two centuries.

"Type D" killer whale genomes reveal long-term small population size and low genetic diversity.

Genome sequences can reveal the extent of inbreeding in small populations. Here, we present the first genomic characterization of type D killer whales, a distinctive eco/morphotype with a circumpolar, subantarctic distribution. Effective population size is the lowest estimated from any killer whale genome and indicates a severe population bottleneck. Consequently, type D genomes show among the highest level of inbreeding reported for any mammalian species (FROH ≥ 0.65). Detected recombination cross-over events of different haplotypes are up to an order of magnitude rarer than in other killer whale genomes studied to date. Comparison of genomic data from a museum specimen of a type D killer whale that stranded in New Zealand in 1955, with 3 modern genomes from the Cape Horn area, reveals high covariance and identity-by-state of alleles, suggesting these genomic characteristics and demographic history are shared among geographically dispersed social groups within this morphotype. Limitations to the insights gained in this study stem from the nonindependence of the 3 closely related modern genomes, the recent coalescence time of most variation within the genomes, and the nonequilibrium population history which violates the assumptions of many model-based methods. Long-range linkage disequilibrium and extensive runs of homozygosity found in type D genomes provide the potential basis for both the distinctive morphology, and the coupling of genetic barriers to gene flow with other killer whale populations.

The evolutionary history of extinct and living lions.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Ancient RNA from Late Pleistocene permafrost and historical canids shows tissue-specific transcriptome survival.

While sequencing ancient DNA (aDNA) from archaeological material is now commonplace, very few attempts to sequence ancient transcriptomes have been made, even from typically stable deposition environments such as permafrost. This is presumably due to assumptions that RNA completely degrades relatively quickly, particularly when dealing with autolytic, nuclease-rich mammalian tissues. However, given the recent successes in sequencing ancient RNA (aRNA) from various sources including plants and animals, we suspect that these assumptions may be incorrect or exaggerated. To challenge the underlying dogma, we generated shotgun RNA data from sources that might normally be dismissed for such study. Here, we present aRNA data generated from two historical wolf skins, and permafrost-preserved liver tissue of a 14,300-year-old Pleistocene canid. Not only is the latter the oldest RNA ever to be sequenced, but it also shows evidence of biologically relevant tissue specificity and close similarity to equivalent data derived from modern-day control tissue. Other hallmarks of RNA sequencing (RNA-seq) data such as exon-exon junction presence and high endogenous ribosomal RNA (rRNA) content confirms our data's authenticity. By performing independent technical library replicates using two high-throughput sequencing platforms, we show not only that aRNA can survive for extended periods in mammalian tissues but also that it has potential for tissue identification. aRNA also has possible further potential, such as identifying in vivo genome activity and adaptation, when sequenced using this technology.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Population genomics of grey wolves and wolf-like canids in North America.

North America is currently home to a number of grey wolf (Canis lupus) and wolf-like canid populations, including the coyote (Canis latrans) and the taxonomically controversial red, Eastern timber and Great Lakes wolves. We explored their population structure and regional gene flow using a dataset of 40 full genome sequences that represent the extant diversity of North American wolves and wolf-like canid populations. This included 15 new genomes (13 North American grey wolves, 1 red wolf and 1 Eastern timber/Great Lakes wolf), ranging from 0.4 to 15x coverage. In addition to providing full genome support for the previously proposed coyote-wolf admixture origin for the taxonomically controversial red, Eastern timber and Great Lakes wolves, the discriminatory power offered by our dataset suggests all North American grey wolves, including the Mexican form, are monophyletic, and thus share a common ancestor to the exclusion of all other wolves. Furthermore, we identify three distinct populations in the high arctic, one being a previously unidentified "Polar wolf" population endemic to Ellesmere Island and Greenland. Genetic diversity analyses reveal particularly high inbreeding and low heterozygosity in these Polar wolves, consistent with long-term isolation from the other North American wolves.

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

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

Ancient DNA suggests modern wolves trace their origin to a Late Pleistocene expansion from Beringia.

Grey wolves (Canis lupus) are one of the few large terrestrial carnivores that have maintained a wide geographical distribution across the Northern Hemisphere throughout the Pleistocene and Holocene. Recent genetic studies have suggested that, despite this continuous presence, major demographic changes occurred in wolf populations between the Late Pleistocene and early Holocene, and that extant wolves trace their ancestry to a single Late Pleistocene population. Both the geographical origin of this ancestral population and how it became widespread remain unknown. Here, we used a spatially and temporally explicit modelling framework to analyse a data set of 90 modern and 45 ancient mitochondrial wolf genomes from across the Northern Hemisphere, spanning the last 50,000 years. Our results suggest that contemporary wolf populations trace their ancestry to an expansion from Beringia at the end of the Last Glacial Maximum, and that this process was most likely driven by Late Pleistocene ecological fluctuations that occurred across the Northern Hemisphere. This study provides direct ancient genetic evidence that long-range migration has played an important role in the population history of a large carnivore, and provides insight into how wolves survived the wave of megafaunal extinctions at the end of the last glaciation. Moreover, because Late Pleistocene grey wolves were the likely source from which all modern dogs trace their origins, the demographic history described in this study has fundamental implications for understanding the geographical origin of the dog.

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.

Whole-Genome Sequencing of African Dogs Provides Insights into Adaptations against Tropical Parasites.

Natural selection in domestic dogs is of great interest in evolutionary biology since dogs have migrated to every inhabited continent of the world alongside humans, and adapted to diverse environments. Here, we explored their demographic history and genetic basis of adaptation to the tropical African environment using whole genome analyses of 19 African indigenous dogs from Nigeria. Demographic analysis suggests that the ancestors of these dogs migrated into Africa from Eurasia 14,000 years ago and underwent a severe founder effect before population expansion. Admixture analysis further reveals that African dog genomes contain about 1.88-3.50% introgression from African golden wolves (Canis anthus). Population genetic analysis identifies 50 positively selected genes linked with immunity, angiogenesis, ultraviolet protection, as well as insulin secretion and sensitivity that may contribute to adaptation to tropical conditions. One of the positively selected genes, adhesion G protein-coupled receptor E1 (ADGRE1), has also been found to be association with severe malaria resistance in African human populations. Functional assessments showed that ADGRE1 provides protective host defense against Plasmodium infections. This result, together with the fact that the inflammatory response to canine babesiosis is similar to complicated falciparum malaria in humans, support the dogs as a model for the study of malaria control and treatment.