Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

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.

The origins and spread of domestic horses from the Western Eurasian steppes.

Domestication of horses fundamentally transformed long-range mobility and warfare1. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling2-4 at Botai, Central Asia around 3500 BC3. Other longstanding candidate regions for horse domestication, such as Iberia5 and Anatolia6, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 BC, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association7 between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 BC8,9 driving the spread of Indo-European languages10. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium BC Sintashta culture11,12.

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.

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.

Synchronous genetic turnovers across Western Eurasia in Late Pleistocene collared lemmings.

Recent palaeogenetic studies indicate a highly dynamic history in collared lemmings (Dicrostonyx spp.), with several demographical changes linked to climatic fluctuations that took place during the last glaciation. At the western range margin of D. torquatus, these changes were characterized by a series of local extinctions and recolonizations. However, it is unclear whether this pattern represents a local phenomenon, possibly driven by ecological edge effects, or a global phenomenon that took place across large geographical scales. To address this, we explored the palaeogenetic history of the collared lemming using a next-generation sequencing approach for pooled mitochondrial DNA amplicons. Sequences were obtained from over 300 fossil remains sampled across Eurasia and two sites in North America. We identified five mitochondrial lineages of D. torquatus that succeeded each other through time across Europe and western Russia, indicating a history of repeated population extinctions and recolonizations, most likely from eastern Russia, during the last 50 000 years. The observation of repeated extinctions across such a vast geographical range indicates large-scale changes in the steppe-tundra environment in western Eurasia during the last glaciation. All Holocene samples, from across the species' entire range, belonged to only one of the five mitochondrial lineages. Thus, extant D. torquatus populations only harbour a small fraction of the total genetic diversity that existed across different stages of the Late Pleistocene. In North American samples, haplotypes belonging to both D. groenlandicus and D. richardsoni were recovered from a Late Pleistocene site in south-western Canada. This suggests that D. groenlandicus had a more southern and D. richardsoni a more northern glacial distribution than previously thought. This study provides significant insights into the population dynamics of a small mammal at a large geographical scale and reveals a rather complex demographical history, which could have had bottom-up effects in the Late Pleistocene steppe-tundra ecosystem.

Serial population extinctions in a small mammal indicate Late Pleistocene ecosystem instability.

The Late Pleistocene global extinction of many terrestrial mammal species has been a subject of intensive scientific study for over a century, yet the relative contributions of environmental changes and the global expansion of humans remain unresolved. A defining component of these extinctions is a bias toward large species, with the majority of small-mammal taxa apparently surviving into the present. Here, we investigate the population-level history of a key tundra-specialist small mammal, the collared lemming (Dicrostonyx torquatus), to explore whether events during the Late Pleistocene had a discernible effect beyond the large mammal fauna. Using ancient DNA techniques to sample across three sites in North-West Europe, we observe a dramatic reduction in genetic diversity in this species over the last 50,000 y. We further identify a series of extinction-recolonization events, indicating a previously unrecognized instability in Late Pleistocene small-mammal populations, which we link with climatic fluctuations. Our results reveal climate-associated, repeated regional extinctions in a keystone prey species across the Late Pleistocene, a pattern likely to have had an impact on the wider steppe-tundra community, and one that is concordant with environmental change as a major force in structuring Late Pleistocene biodiversity.

On the origin of the Norwegian lemming.

The Pleistocene glacial cycles resulted in significant changes in species distributions, and it has been discussed whether this caused increased rates of population divergence and speciation. One species that is likely to have evolved during the Pleistocene is the Norwegian lemming (Lemmus lemmus). However, the origin of this species, both in terms of when and from what ancestral taxon it evolved, has been difficult to ascertain. Here, we use ancient DNA recovered from lemming remains from a series of Late Pleistocene and Holocene sites to explore the species' evolutionary history. The results revealed considerable genetic differentiation between glacial and contemporary samples. Moreover, the analyses provided strong support for a divergence time prior to the Last Glacial Maximum (LGM), therefore likely ruling out a postglacial colonization of Scandinavia. Consequently, it appears that the Norwegian lemming evolved from a small population that survived the LGM in an ice-free Scandinavian refugium.

Holarctic genetic structure and range dynamics in the woolly mammoth.

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32-34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20-15 kyr BP, when a severe population size decline occurred.

The human-initiated model of wolf domestication - An expansion based on human-dingo relations in Aboriginal Australia.

The historically known relationship of interspecies companionship between Aboriginal foraging communities in Australia and free-ranging dingoes provides a model for understanding the human-canid relations that gave rise to the first domesticated dogs. Here, we propose that a broadly similar relationship might have developed early in time between wild-living wolves and mobile groups of foragers in Late Pleistocene Eurasia, with hunter-gatherers routinely raiding wild wolf dens for pre-weaned pups, which were socialized to humans and kept in camp as tamed companions ("pets"). We outline a model in which captive wolf pups that reverted to the wild to breed when they were sexually mature established their territories in the vicinity of foraging communities - in a "liminal" ecological zone between humans and truly wild-living wolves. Many (or most) of the wolf pups humans took from the wilderness to rear in camp may have derived from these liminal dens where the breeding pairs had been under indirect human selection for tameness over many generations. This highlights the importance of the large seasonal hunting/aggregation camps associated with mammoth kill-sites in Gravettian/Epigravettian central Europe. Large numbers of foragers gathered regularly at these locations during the wild wolf birthing season. We infer that if a pattern of this kind occurred over long periods of time then there might have been a pronounced effect on genetic variation in free-ranging wolves that denned and whelped in the liminal zones in the vicinity of these human seasonal aggregation sites. The argument is not that wolves were domesticated in central Europe. Rather, it is this pattern of hunter-gatherers who caught and reared wild wolf pups gathering seasonally in large numbers that might have been the catalyst for the early changes leading to the first domesticated dogs - whether in western Eurasia or further afield.

Relationships of Late Pleistocene giant deer as revealed by Sinomegaceros mitogenomes from East Asia.

The giant deer, widespread in northern Eurasia during the Late Pleistocene, have been classified as western Megaloceros and eastern Sinomegaceros through morphological studies. While Megaloceros's evolutionary history has been unveiled through mitogenomes, Sinomegaceros remains molecularly unexplored. Herein, we generated mitogenomes of giant deer from East Asia. We find that, in contrast to the morphological differences between Megaloceros and Sinomegaceros, they are mixed in the mitochondrial phylogeny, and Siberian specimens suggest a range contact or overlap between these two groups. Meanwhile, one deep divergent clade and another surviving until 20.1 thousand years ago (ka) were detected in northeastern China, the latter implying this area as a potential refugium during the Last Glacial Maximum (LGM). Moreover, stable isotope analyses indicate correlations between climate-introduced vegetation changes and giant deer extinction. Our study demonstrates the genetic relationship between eastern and western giant deer and explores the promoters of their extirpation in northern East Asia.

Predormancy omnivory in European cave bears evidenced by a dental microwear analysis of Ursus spelaeus from Goyet, Belgium.

Previous morphological and isotopic studies indicate that Late Pleistocene cave bear (Ursus spelaeus) diet ranged from mostly vegetarian to omnivory or even carnivory. However, such analyses do not provide information on seasonal diets, and only provide an average record of diet. A dental microwear analysis of 43 young and adult individuals demonstrate that, during the predormancy period, cave bears from Goyet (Late Pleistocene, Belgium) were not strictly herbivorous, but had a mixed diet composed of hard items (e.g., possibly bone), invertebrates (e.g., insects), meat (ungulates, small vertebrates), and/or plant matter (hard mast, seeds, herbaceous vegetations, and fruits). Therefore, our results indicate that cave bears at Goyet were generalist omnivores during the predormancy period, which is consistent with current data on the dietary ecology of extant bears during this season. These data also raise questions about the ecological role and causes of the extinction of cave bears.

Morphological differences between putative Paleolithic dogs and wolves: A commentary to Janssens et al. (2021).

Janssens et al. (2021, doi: 10.1002/ar.24624) recently commented on our article (Galeta et al., 2021, doi: 10.1002/ar.24500) regarding the morphological differences between putative Paleolithic dog and Pleistocene wolf crania. The authors argued that these differences reflect the normal population variation of wolves, that some of the cranial measurements used do not reflect morphological changes during domestication, and that our canid dataset was small because we inexplicably omitted several specimens we analyzed in our previous publications. In this commentary, we briefly address the issue of within and between morpho-population variability. The results based on our canid sample suggest that the magnitude of morphological differences between distinct morpho-populations (i.e., recent northern dogs and wolves) is at least twice as large as that observed within morpho-populations (between two groups of recent northern wolves segregated by cluster analysis). The morphological differences between putative Paleolithic dogs and Pleistocene wolves are relatively large, which may indicate that they did not likely represent a single Late Pleistocene morpho-population. Finally, we clarified the rationale behind the composition of our 2021 dataset to show that we did not adjust the list of the analyzed specimens. Although the sample size was small, the randomization analysis published in 2021 confirmed that the unbalanced composition of the reference sample did not affect the reliability of the morphological segregation of putative Paleolithic dogs and Pleistocene wolves.

Population dynamics and demographic history of Eurasian collared lemmings.

BACKGROUND: Ancient DNA studies suggest that Late Pleistocene climatic changes had a significant effect on population dynamics in Arctic species. The Eurasian collared lemming (Dicrostonyx torquatus) is a keystone species in the Arctic ecosystem. Earlier studies have indicated that past climatic fluctuations were important drivers of past population dynamics in this species. RESULTS: Here, we analysed 59 ancient and 54 modern mitogenomes from across Eurasia, along with one modern nuclear genome. Our results suggest population growth and genetic diversification during the early Late Pleistocene, implying that collared lemmings may have experienced a genetic bottleneck during the warm Eemian interglacial. Furthermore, we find multiple temporally structured mitogenome clades during the Late Pleistocene, consistent with earlier results suggesting a dynamic late glacial population history. Finally, we identify a population in northeastern Siberia that maintained genetic diversity and a constant population size at the end of the Pleistocene, suggesting suitable conditions for collared lemmings in this region during the increasing temperatures associated with the onset of the Holocene. CONCLUSIONS: This study highlights an influence of past warming, in particular the Eemian interglacial, on the evolutionary history of the collared lemming, along with spatiotemporal population structuring throughout the Late Pleistocene.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFß signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes.

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.

Morphological evidence for early dog domestication in the European Pleistocene: New evidence from a randomization approach to group differences.

The antiquity of the wolf/dog domestication has been recently pushed back in time from the Late Upper Paleolithic (~14,000 years ago) to the Early Upper Paleolithic (EUP; ~36,000 years ago). Some authors questioned this early dog domestication claiming that the putative (EUP) Paleolithic dogs fall within the morphological range of recent wolves. In this study, we reanalyzed a data set of large canid skulls using unbalanced- and balanced-randomized discriminant analyses to assess whether the putative Paleolithic dogs are morphologically unique or whether they represent a subsample of the wolf morpho-population. We evaluated morphological differences between 96 specimens of the 4 a priori reference groups (8 putative Paleolithic dogs, 41 recent northern dogs, 7 Pleistocene wolves, and 40 recent northern wolves) using discriminant analysis based on 5 ln-transformed raw and allometrically size-adjusted cranial measurements. Putative Paleolithic dogs are classified with high accuracies (87.5 and 100.0%, cross-validated) and randomization experiment suggests that these classification rates cannot be exclusively explained by the small and uneven sample sizes of reference groups. It indicates that putative Upper Paleolithic dogs may represent a discrete canid group with morphological signs of domestication (a relatively shorter skull and wider palate and braincase) that distinguish them from sympatric Pleistocene wolves. The present results add evidence to the view that these specimens could represent incipient Paleolithic dogs that were involved in daily activities of European Upper Paleolithic forager groups.

Genomes of Pleistocene Siberian Wolves Uncover Multiple Extinct Wolf Lineages.

Extant Canis lupus genetic diversity can be grouped into three phylogenetically distinct clades: Eurasian and American wolves and domestic dogs.1 Genetic studies have suggested these groups trace their origins to a wolf population that expanded during the last glacial maximum (LGM)1-3 and replaced local wolf populations.4 Moreover, ancient genomes from the Yana basin and the Taimyr peninsula provided evidence of at least one extinct wolf lineage that dwelled in Siberia during the Pleistocene.35 Previous studies have suggested that Pleistocene Siberian canids can be classified into two groups based on cranial morphology. Wolves in the first group are most similar to present-day populations, although those in the second group possess intermediate features between dogs and wolves.67 However, whether this morphological classification represents distinct genetic groups remains unknown. To investigate this question and the relationships between Pleistocene canids, present-day wolves, and dogs, we resequenced the genomes of four Pleistocene canids from Northeast Siberia dated between >50 and 14 ka old, including samples from the two morphological categories. We found these specimens cluster with the two previously sequenced Pleistocene wolves, which are genetically more similar to Eurasian wolves. Our results show that, though the four specimens represent extinct wolf lineages, they do not form a monophyletic group. Instead, each Pleistocene Siberian canid branched off the lineage that gave rise to present-day wolves and dogs. Finally, our results suggest the two previously described morphological groups could represent independent lineages similarly related to present-day wolves and dogs.

Were ancient foxes far more carnivorous than recent ones?-Carnassial morphological evidence.

Crown shape variation of the first lower molar in the arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) was analyzed using five groups of morphotypes. Carnassial morphologies were compared between the species and between spatially and temporally distant populations: one Late Pleistocene (n = 45) and seven modern populations of the arctic fox (n = 259), and one Late Pleistocene (n = 35) and eight modern populations of the red fox (n = 606). The dentition of Holocene red foxes had larger morphotype variability than that of arctic foxes. The lower carnassials of the red fox kept have some primitive characters (additional cusps and stylids, complex shape of transverse cristid), whereas the first lower molars of the arctic fox have undergone crown shape simplification, with the occlusal part of the tooth undergoing a more pronounced adaptation to a more carnivorous diet. From the Late Pleistocene of Belgium to the present days, the arctic fox's crown shape has been simplified and some primitive characters have disappeared. In the red fox chronological changes in the morphology of the lower carnassials were not clearly identified. The phyletic tree based on morphotype carnassial characteristics indicated the distinctiveness of both foxes: in the arctic fox line, the ancient population from Belgium and recent Greenland made separate branches, whereas in the red foxes the ancient population from Belgium was most similar to modern red foxes from Belgium and Italy.

Early Pleistocene origin and extensive intra-species diversity of the extinct cave lion.

The cave lion is an extinct felid that was widespread across the Holarctic throughout the Late Pleistocene. Its closest extant relative is the lion (Panthera leo), but the timing of the divergence between these two taxa, as well as their taxonomic ranking are contentious. In this study we analyse 31 mitochondrial genome sequences from cave lion individuals that, through a combination of 14C and genetic tip dating, are estimated to be from dates extending well into the mid-Pleistocene. We identified two deeply diverged and well-supported reciprocally monophyletic mitogenome clades in the cave lion, and an additional third distinct lineage represented by a single individual. One of these clades was restricted to Beringia while the other was prevalent across western Eurasia. These observed clade distributions are in line with previous observations that Beringian and European cave lions were morphologically distinct. The divergence dates for these lineages are estimated to be far older than those between extant lions subspecies. By combining our radiocarbon tip-dates with a split time prior that takes into account the most up-to-date fossil stem calibrations, we estimated the mitochondrial DNA divergence between cave lions and lions to be 1.85 Million ya (95% 0.52- 2.91 Mya). Taken together, these results support previous hypotheses that cave lions existed as at least two subspecies during the Pleistocene, and that lions and cave lions were distinct species.