Yakutia, Sakha Republic, in the Siberian Far East, represents one of the coldest places on Earth, with winter record temperatures dropping below −70°C. Nevertheless, Yakutian horses survive all year round in the open air due to striking phenotypic adaptations, including compact body conformations, extremely hairy winter coats, and acute seasonal differences in metabolic activities. The evolutionary origins of Yakutian horses and the genetic basis of their adaptations remain, however, contentious. Here, we present the complete genomes of nine present-day Yakutian horses and two ancient specimens dating from the early 19th century and ∼5,200 y ago. By comparing these genomes with the genomes of two Late Pleistocene, 27 domesticated, and three wild Przewalski's horses, we find that contemporary Yakutian horses do not descend from the native horses that populated the region until the mid-Holocene, but were most likely introduced following the migration of the Yakut people a few centuries ago. Thus, they represent one of the fastest cases of adaptation to the extreme temperatures of the Arctic. We find cisregulatory mutations to have contributed more than nonsynonymous changes to their adaptation, likely due to the comparatively limited standing variation within gene bodies at the time the population was founded. Genes involved in hair development, body size, and metabolic and hormone signaling pathways represent an essential part of the Yakutian horse adaptive genetic toolkit. Finally, we find evidence for convergent evolution with native human populations and woolly mammoths, suggesting that only a few evolutionary strategies are compatible with survival in extremely cold environments. The Yakutian horse is the most northerly distributed horse on the planet and certainly the most resistant to cold. In contrast to SignificanceYakutia is among the coldest regions in the Northern Hemisphere, showing ∼40% of its territory above the Arctic Circle. Native horses are particularly adapted to this environment, with body sizes and thick winter coats minimizing heat loss. We sequenced complete genomes of two ancient and nine present-day Yakutian horses to elucidate their evolutionary origins. We find that the contemporary population descends from domestic livestock, likely brought by early horse-riders who settled in the region a few centuries ago. The metabolic, anatomical, and physiological adaptations of these horses therefore emerged on very short evolutionary time scales. We show the relative importance of regulatory changes in the adaptive process and identify genes independently selected in cold-adapted human populations and woolly mammoths.

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.

A megaslump at Batagaika, in northern Yakutia, exposes a remarkable stratigraphic sequence of permafrost deposits 50-80 m thick. To determine their potential for answering key questions about Quaternary environmental and climatic change in northeast Siberia, we carried out a reconnaissance study of their cryostratigraphy and paleoecology, supported by four rangefinder 14 C ages. The sequence includes two ice complexes separated by a unit of fine sand containing narrow syngenetic ice wedges and multiple paleosols. Overall, the sequence developed as permafrost grew syngenetically through an eolian sand sheet aggrading on a hillslope. Wood remains occur in two forest beds, each associated with a reddened weathering horizon. The lower bed contains high amounts of Larix pollen (>20%), plus small amounts of Picea and Pinus pumila, and is attributed to interglacial conditions. Pollen from the overlying sequence is dominated by herbaceous taxa (~70%-80%) attributed to an open tundra landscape during interstadial climatic conditions. Of three hypothetical age schemes considered, we tentatively attribute much of the Batagaika sequence to Marine Oxygen Isotope Stage (MIS) 3. The upper and lower forest beds may represent a mid-MIS 3 optimum and MIS 5, respectively, although we cannot discount alternative attributions to MIS 5 and 7.

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.

The present phylogeographic pattern of red deer in Eurasia is not only a result of the contraction of their distribution range into glacial refugia and postglacial expansion, but probably also an effect of replacement of some red deer s.l. mtDNA lineages by others during the last 50 000 years. To better recognize this process, we analysed 501 sequences of mtDNA cytochrome b, including 194 ancient and 75 contemporary samples newly obtained for this study. The inclusion of 161 radiocarbon-dated samples enabled us to study the phylogeny in a temporal context and conduct divergence-time estimation and molecular dating. Depending on methodology, our estimate of divergence between Cervus elaphus and Cervus canadensis varied considerably (370 000 or 1.37 million years BP, respectively). The divergence times of genetic lineages and haplogroups corresponded to large environmental changes associated with stadials and interstadials of the Late Pleistocene. Due to the climatic oscillations, the distribution of C. elaphus and C. canadensis fluctuated in north–south and east–west directions. Some haplotypes dated to pre-Last Glacial Maximum periods were not detected afterwards, representing possibly extinct populations. We indicated with a high probability the presence of red deer sensu lato in south-eastern Europe and western Asia during the Last Glacial Maximum.