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.

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe’s largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

The extant diversity of horses (family Equidae) represents a small fraction of that occurring over their evolutionary history. One such lost lineage is the subgenus Sussemionus , which is thought to have become extinct during the Middle Pleistocene. However, recent molecular studies and morphological analysis have revealed that one of their representatives, E . ovodovi , did exist in Siberia during the Late Pleistocene. Fossil materials of E . ovodovi have thus far only been found in Russia. In this study, we extracted DNA from three equid fossil specimens excavated from northeastern China dated at 12,770–12,596, 29,525–28,887 and 40,201–38,848 cal. yBP, respectively, and retrieved three near-complete mitochondrial genomes from the specimens. Phylogenetic analyses cluster the Chinese haplotypes together with previously published Russian E . ovodovi , strongly supporting the assignment of these samples to this taxon. The molecular identification of E . ovodovi in northeastern China extends the known geographical range of this fossil species by several thousand kilometers to the east. The estimated coalescence time of all E . ovodovi haplotypes is approximately 199 Kya, with the Chinese haplotypes coalescing approximately 130 Kya. With a radiocarbon age of 12,770–12,596 cal. yBP, the youngest sample in this study represents the first E . ovodovi sample dating to the terminal Pleistocene, moving the extinction date of this species forwards considerably compared to previously documented fossils. Overall, comparison of our three mitochondrial genomes with the two published ones suggests a genetic diversity similar to several extant species of the genus Equus .

High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present ‘TAPAS’, (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material.

No abstract