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

Coyotes can do 'puppy dog eyes' too: comparing interspecific variation in Canis facial expression muscles

Facial expressions are critical for non-verbal communication. The Canis genus epitomizes the interplay between behaviour and morphology in the evolution of non-verbal communication. Recent work suggests that the levator anguli oculi medialis (LAOM) muscle is unique to dogs (Canis familiaris) within the Canis genus and evolved due to domestication. The LAOM raises the inner eyebrows, resulting in the 'puppy dog eyes' expression. Here, we test whether the LAOM is a derived trait in dogs by (i) examining the facial expression muscles of a closely related and ancestral wild Canis species, the coyote (C. latrans) and (ii) comparing our results with other Canis and canid taxa. We discover that coyotes have a well-developed LAOM like dogs, which differs from the modified/absent LAOM in grey wolves. Our findings challenge the hypothesis that the LAOM developed due to domestication. We suggest that the LAOM is a basal trait that was lost in grey wolves. Additionally, we find inter- and intraspecific variations in the size of the muscles of the outer ear, forehead, lips and rostrum, indicating potential adaptations related to sensory perception, communication and individual-level functional variations within canids. Together, this research expands our knowledge of facial expressions, their evolution and their role in communication.

Genomic analysis of wolves from Pakistan clarifies boundaries among three divergent wolf lineages

Among the three main divergent lineages of gray wolf (Canis lupus), the Holarctic lineage is the most widespread and best studied, particularly in North America and Europe. Less is known about Tibetan (also called Himalayan) and Indian wolf lineages in southern Asia, especially in areas surrounding Pakistan where all three lineages are thought to meet. Given the endangered status of the Indian wolf in neighboring India and unclear southwestern boundary of the Tibetan wolf range, we conducted mitochondrial and genome-wide sequencing of wolves from Pakistan and Kyrgyzstan. Sequences of the mitochondrial D-loop region of 81 wolves from Pakistan indicated contact zones between Holarctic and Indian lineages across the northern and western mountains of Pakistan. Reduced-representation genome sequencing of eight wolves indicated an east-to-west cline of Indian to Holarctic ancestry, consistent with a contact zone between these two lineages in Pakistan. The western boundary of the Tibetan lineage corresponded to the Ladakh region of India's Himalayas with a narrow zone of admixture spanning this boundary from the Karakoram Mountains of northern Pakistan into Ladakh, India. Our results highlight the conservation significance of Pakistan's wolf populations, especially the remaining populations in Sindh and Southern Punjab that represent the highly endangered Indian lineage.

Conservation genomics of wolves: The global impact of RK Wayne's research

RK Wayne has arguably been the most influential geneticist of canids, famously promoting the conservation of wolves in his homeland, the United States. His influence has been felt in other countries and regions outside the contiguous United States, where he inspired others, also including former graduate students and research fellows of his, to use modern molecular techniques to examine the evolutionary biology of canids to inform the conservation and management of wolves. In this review, we focus on the implications of Wayne's work on wolves outside the United States. He envisioned a clear future for wolf conservation research, involving the study of wolves' ecological and genetic diversity, and the description of ecotypes requiring conservation. He also documented widespread hybridization among canids and introgression of DNA from domestic dogs to wolves, a process that started dozens of thousands of years ago. His work therefore calls for innovative studies, such as examining the potential fitness benefits of introgression. Inspired by his results, for example, on the purging of deleterious alleles in small populations, wolf researchers should use novel molecular tools to challenge other conservation genetics paradigms. Overall, RK Wayne's work constitutes a call for answers, which as scientists or citizens concerned with conservation matters, we are obliged to address, as we contribute to monitoring and maintaining biodiversity during our period of dramatic transformations of the biosphere.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R2>0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

In Search of the Elusive North: Evolutionary History of the Arctic Fox (Vulpes lagopus) in the Palearctic from the Late Pleistocene to the Recent Inferred from Mitogenomic Data

Despite the high level of interest, the population history of arctic foxes during the Late Pleistocene and Holocene remains poorly understood. Here we aimed to fill gaps in the demographic and colonization history of the arctic fox by analyzing new ancient DNA data from fossil specimens aged from 50 to 1 thousand years from the Northern and Polar Urals, historic DNA from museum specimens from the Novaya Zemlya Archipelago and the Taymyr Peninsula and supplementing these data by previously published sequences of recent and extinct arctic foxes from other regions. This dataset was used for reconstruction of a time-calibrated phylogeny and a temporal haplotype network covering four time intervals: Late Pleistocene (ranging from 30 to 13 thousand years bp), Holocene (ranging from 4 to 1 thousand years bp), historical (approximately 150 years), and modern. Our results revealed that Late Pleistocene specimens showed no genetic similarity to either modern or historical specimens, thus supporting the earlier hypothesis on local extinction rather than habitat tracking.

Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Evolutionary history of the extinct wolf population from France in the context of global phylogeographic changes throughout the Holocene

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

The complete mitochondrial genome and phylogenetic characterization of two putative subspecies of golden jackal (Canis aureus cruesemanni and Canis aureus moreotica)

The golden jackal (Canis aureus) is a canid species found across southern Eurasia. Several subspecies of this animal have been genetically studied in regions such as Europe, the Middle East, and India. However, one subspecies that lacks current research is the Indochinese jackal (Canis aureus cruesemanni), which is primarily found in Southeast Asia. Using a genome skimming approach, we assembled the first complete mitochondrial genome for an Indochinese jackal from Thailand. To expand the number of available Canis aureus mitogenomes, we also assembled and sequenced the first complete mitochondrial genome of a golden jackal from Turkey, representing the C. a. moreotica subspecies. The mitogenomes contained 37 annotated genes and are 16,729 bps (C. a. cruesemanni) and 16,669 bps (C. a. moreotica) in length. Phylogenetic analysis with 26 additional canid mitogenomes and analyses of a cytochrome b gene-only data set together support the Indochinese jackal as a distinct and early-branching lineage among golden jackals, thereby supporting its recognition as a possible subspecies. These analyses also demonstrate that the golden jackal from Turkey is likely not a distinct lineage due to close genetic relationships with golden jackals from India and Israel.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

From high masked to high realized genetic load in inbred Scandinavian wolves

When new mutations arise at functional sites they are more likely to impair than improve fitness. If not removed by purifying selection, such deleterious mutations will generate a genetic load that can have negative fitness effects in small populations and increase the risk of extinction. This is relevant for the highly inbred Scandinavian wolf (Canis lupus) population, founded by only three wolves in the 1980s and suffering from inbreeding depression. We used functional annotation and evolutionary conservation scores to study deleterious variation in a total of 209 genomes from both the Scandinavian and neighbouring wolf populations in northern Europe. The masked load (deleterious mutations in heterozygote state) was highest in Russia and Finland with deleterious alleles segregating at lower frequency than neutral variation. Genetic drift in the Scandinavian population led to the loss of ancestral alleles, fixation of deleterious variants and a significant increase in the per-individual realized load (deleterious mutations in homozygote state; an increase by 45% in protein-coding genes) over five generations of inbreeding. Arrival of immigrants gave a temporary genetic rescue effect with ancestral alleles re-entering the population and thereby shifting deleterious alleles from homozygous into heterozygote genotypes. However, in the absence of permanent connectivity to Finnish and Russian populations, inbreeding has then again led to the exposure of deleterious mutations. These observations provide genome-wide insight into the magnitude of genetic load and genetic rescue at the molecular level, and in relation to population history. They emphasize the importance of securing gene flow in the management of endangered populations.

Pleistocene parades of carnivores into North America

The distribution and movement of species, broadly known as biogeography, is one of the fundamental subfields of ecology and evolutionary biology. However, significant mysteries remain about the processes that gave rise to the modern distribution of biodiversity across the globe. Over the last several decades, the genetic study of ancient and subfossil specimens has started to shed light on past migrations of some species, with a particular focus on humans and megafauna. In this issue of Molecular Ecology, Salis et al. (2021) use ancient mitogenomes and a new phylogeographic method to add an important new piece of evidence to the mystery of megafaunal migrations into North America during the Pleistocene. They found a striking synchronicity of brown bear (Ursus arctos) and lion (Panthera spp.) migrations across the Bering Land Bridge at several time points during the late Pleistocene, which highlights the lasting impact of sea level change on the prehistoric and modern dispersal of terrestrial carnivores across continents.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

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.

Tail Wags Dog’s SINE: Retropositional Mechanisms of Can SINE Depend on Its A-Tail Structure

Simple Summary

The genomes of higher organisms including humans are invaded by millions of repetitive elements (transposons), which can sometimes be deleterious or beneficial for hosts. Many aspects of the mechanisms underlying the expansion of transposons in the genomes remain unclear. Short retrotransposons (SINEs) are one of the most abundant classes of genomic repeats. Their amplification relies on two major processes: transcription and reverse transcription. Here, short retrotransposons of dogs and other canids called Can SINE were analyzed. Their amplification was extraordinarily active in the wolf and, particularly, dog breeds relative to other canids. We also studied a variation of their transcription mechanism involving the polyadenylation of transcripts. An analysis of specific signals involved in this process allowed us to conclude that Can SINEs could alternate amplification with and without polyadenylation in their evolution. Understanding the mechanisms of transposon replication can shed light on the mechanisms of genome function.

Abstract

SINEs, non-autonomous short retrotransposons, are widespread in mammalian genomes. Their transcripts are generated by RNA polymerase III (pol III). Transcripts of certain SINEs can be polyadenylated, which requires polyadenylation and pol III termination signals in their sequences. Our sequence analysis divided Can SINEs in canids into four subfamilies, older a1 and a2 and younger b1 and b2. Can_b2 and to a lesser extent Can_b1 remained retrotranspositionally active, while the amplification of Can_a1 and Can_a2 ceased long ago. An extraordinarily high Can amplification was revealed in different dog breeds. Functional polyadenylation signals were analyzed in Can subfamilies, particularly in fractions of recently amplified, i.e., active copies. The transcription of various Can constructs transfected into HeLa cells proposed AATAAA and (TC)_n as functional polyadenylation signals. Our analysis indicates that older Can subfamilies (a1, a2, and b1) with an active transcription terminator were amplified by the T⁺ mechanism (with polyadenylation of pol III transcripts). In the currently active Can_b2 subfamily, the amplification mechanisms with (T⁺) and without the polyadenylation of pol III transcripts (T⁻) irregularly alternate. The active transcription terminator tends to shorten, which renders it nonfunctional and favors a switch to the T⁻ retrotransposition. The activity of a truncated terminator is occasionally restored by its elongation, which rehabilitates the T⁺ retrotransposition for a particular SINE copy.

The genome and diet of a 35,000-year-old Canis lupus specimen from the Paleolithic painted cave, Chauvet-Pont d'Arc, France

The Chauvet-Pont-d'Arc Cave (Ardèche, France) contains some of the oldest Paleolithic paintings recorded to date, as well as thousands of bones of the extinct cave bear, and some remains and footprints of other animals. As part of the interdisciplinary research project devoted to this reference cave site, we analyzed a coprolite collected within the deep cave. AMS radiocarbon dating of bone fragments from the coprolite yielded an age of 30,450 ± 550 RC yr. BP (AAR-19656; 36,150-34,000 cal BP), similar to ages assigned to Paleolithic artwork and cave bear remains from the same cave sector. Using high-throughput shotgun DNA sequencing, we demonstrated a high abundance of canid DNA and lesser amounts of DNA from the extinct cave bear. We interpret the sample as feces from a canid that had consumed cave bear tissue. The high amount of canid DNA allowed us to reconstruct a complete canid mitochondrial genome sequence (average coverage: 83×) belonging to a deeply divergent clade of extinct mitochondrial wolf lineages that are most closely related to coeval (~35 ka) Belgian wolves. Analysis of the nuclear genome yielded a similar coverage for the X chromosome (2.4×) and the autosomes (range: 2.3-3.2×), indicating that the Chauvet canid was a female. Comparing the relationship of the nuclear genome of this specimen with that of a variety of canids, we found it more closely related to gray wolves' genomes than to other wild canid or dog genomes, especially wolf genomes from Europe and the Middle East. We conclude that the coprolite is feces from an animal within an extinct wolf lineage. The consumption of cave bear by this wolf likely explains its intrusion into the dark cave sectors and sheds new light on the paleoecology of a major cave site.

Demography and evolutionary history of grey wolf populations around the Bering Strait

Glacial and interglacial periods throughout the Pleistocene have been substantial drivers of change in species distributions. Earlier analyses suggested that modern grey wolves (Canis lupus) trace their origin to a single Late Pleistocene Beringian population that expanded east and westwards, starting c. 25,000 years ago (ya). Here, we examined the demographic and phylogeographic histories of extant populations around the Bering Strait with wolves from two inland regions of the Russian Far East (RFE) and one coastal and two inland regions of North‐western North America (NNA), genotyped for 91,327 single nucleotide polymorphisms. Our results indicated that RFE and NNA wolves had a common ancestry until c. 34,400 ya, suggesting that these populations started to diverge before the previously proposed expansion out of Beringia. Coastal and inland NNA populations diverged c. 16,000 ya, concordant with the minimum proposed date for the ecological viability of the migration route along the Pacific Northwest coast. Demographic reconstructions for inland RFE and NNA populations reveal spatial and temporal synchrony, with large historical effective population sizes that declined throughout the Pleistocene, possibly reflecting the influence of broadscale climatic changes across continents. In contrast, coastal NNA wolves displayed a consistently lower effective population size than the inland populations. Differences between the demographic history of inland and coastal wolves may have been driven by multiple ecological factors, including historical gene flow patterns, natural landscape fragmentation, and more recent anthropogenic disturbance.

Overview of Current Research on Wolves in Russia

This paper provides an overview of wolf research in Russia at the beginning of the 21st century. Wolf research covered various directions, including population density estimation, management methods and minimization of human-wildlife conflicts, general and behavioral ecology, behavior, wolf population genetics and morphology, paleontology, dog domestication, helminthology and the wolves’ role in the rabies transmission. Some studies are performed with state-of-art methodology using molecular genetics, mathematical modeling, camera traps, and GPS telemetry.

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) lived^1–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.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

The Diversity in the Genus Canis Challenges Conservation Biology: A Review of Available Data on Asian Wolves

Taxa belonging to the Genus Canis can challenge taxonomists because species boundaries and distribution ranges are often gradual. Species delineation within Canis is currently not based on consistent criteria, and is hampered by geographical bias and lack of taxonomic research. But a consistent taxonomy is critical, given its importance for assigning legal protection, conservation priorities, and financial resources. We carried out a qualitative review of the major wolf lineages so far identified from Asia from historical to contemporary time and considered relevant morphological, ecological, and genetic evidence. We present full mitochondrial phylogenies and genetic distances between these lineages. This review aims to summarize the available data on contemporary Asian wolf lineages within the context of the larger phylogenetic Canis group and to work toward a taxonomy that is consistent within the Canidae. We found support for the presence and taxon eligibility of Holarctic gray, Himalayan/Tibetan, Indian, and Arabian wolves in Asia and recommend their recognition at the taxonomic levels consistent within the group.

Placing Ancient DNA Sequences into Reference Phylogenies

Joint phylogenetic analysis of ancient DNA (aDNA) with modern phylogenies is hampered by low sequence coverage and post-mortem deamination, often resulting in overconservative or incorrect assignment. We provide a new efficient likelihood-based workflow, pathPhynder, that takes advantage of all the polymorphic sites in the target sequence. This effectively evaluates the number of ancestral and derived alleles present on each branch and reports the most likely placement of an ancient sample in the phylogeny and a haplogroup assignment, together with alternatives and supporting evidence. To illustrate the application of pathPhynder, we show improved Y chromosome assignments for published aDNA sequences, using a newly compiled Y variation data set (120,908 markers from 2,014 samples) that significantly enhances Y haplogroup assignment for low coverage samples. We apply the method to all published male aDNA samples from Africa, giving new insights into ancient migrations and the relationships between ancient and modern populations. The same software can be used to place samples with large amounts of missing data into other large non-recombining phylogenies such as the mitochondrial tree.

Informing conservation strategies with museum genomics: Long-term effects of past anthropogenic persecution on the elusive European wildcat

Like many carnivore species, European wildcats (Felis silvestris) have suffered severe anthropogenic population declines in the past, resulting in a strong population bottleneck at the beginning of the 20th century. In Germany, the species has managed to survive its near extinction in small isolated areas and is currently recolonizing former habitats owing to legal protection and concerted conservation efforts. Here, we SNP-genotyped and mtDNA-sequenced 56 historical and 650 contemporary samples to assess the impact of massive persecution on genetic diversity, population structure, and hybridization dynamics of wildcats. Spatiotemporal analyses suggest that the presumed postglacial differentiation between two genetically distinct metapopulations in Germany is in fact the result of the anthropogenic bottleneck followed by re-expansion from few secluded refugia. We found that, despite the bottleneck, populations experienced no severe genetic erosion, nor suffered from elevated inbreeding or showed signs of increased hybridization with domestic cats. Our findings have significant implications for current wildcat conservation strategies, as the data analyses show that the two presently recognized wildcat population clusters should be treated as a single conservation unit. Although current populations appear under no imminent threat from genetic factors, fostering connectivity through the implementation of forest corridors will facilitate the preservation of genetic diversity and promote long-term viability. The present study documents how museum collections can be used as essential resource for assessing long-term anthropogenic effects on natural populations, for example, regarding population structure and the delineation of appropriate conservation units, potentially informing todays' species conservation.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Disentangling the admixed trails of grey wolf evolution

The predominant phylogenetic patterns within a genome do not always reflect correctly the history of evolutionary divergence and speciation, and the true phylogenetic signal tends to be concentrated within low-recombination regions of the genome. In this issue of Molecular Ecology, Hennelly et al. show that this is also the case for intraspecific relationships that are characterized by considerable gene flow between lineages. The study reconstructs the phylogenetic relationships of Indian and Tibetan wolves with other grey wolf (Canis lupus) populations worldwide, and demonstrates that these two populations represent phylogenetically distinct lineages. This inference was supported by using low-recombination regions of autosomal chromosomes and the X chromosome, which proved to be essential for correct inference of the lineage splitting order. Their study illustrates the power of analytical approaches that implement knowledge of genome evolution patterns to reconstruct complex intraspecific evolutionary relationships. The study also provides a compelling example of the application of modern phylogenomic approaches in the identification of evolutionarily significant units for the purpose of species conservation.

A multi-proxy approach to exploring Homo sapiens' arrival, environments and adaptations in Southeast Asia

The capability of Pleistocene hominins to successfully adapt to different types of tropical forested environments has long been debated. In order to investigate environmental changes in Southeast Asia during a critical period for the turnover of hominin species, we analysed palaeoenvironmental proxies from five late Middle to Late Pleistocene faunas. Human teeth discoveries have been reported at Duoi U'Oi, Vietnam (70-60 ka) and Nam Lot, Laos (86-72 ka). However, the use of palaeoproteomics allowed us to discard the latter, and, to date, no human remains older than ~ 70 ka are documented in the area. Our findings indicate that tropical rainforests were highly sensitive to climatic changes over that period, with significant fluctuations of the canopy forests. Locally, large-bodied faunas were resilient to these fluctuations until the cooling period of the Marine Isotope Stage 4 (MIS 4; 74-59 ka) that transformed the overall biotope. Then, under strong selective pressures, populations with new phenotypic characteristics emerged while some other species disappeared. We argue that this climate-driven shift offered new foraging opportunities for hominins in a novel rainforest environment and was most likely a key factor in the settlement and dispersal of our species during MIS 4 in SE Asia.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

A mummified Pleistocene gray wolf pup

In July 2016, a mummified carcass of an ancient wolf (Canis lupus) pup (specimen YG 648.1) was discovered in thawing permafrost in the Klondike goldfields, near Dawson City, Yukon, Canada (Figure 1A). The wolf pup mummy was recovered along a small tributary of Last Chance Creek during hydraulic thawing that exposed the permafrost sediment in which it was preserved. This mummified wolf pup is important to the local Tr'ondëk Hwëch'in people, who named it Zhùr, meaning 'wolf' in the Hän language of their community. Here, we report detailed morphometric, isotopic, and genetic analyses of Zhùr that reveal details of her appearance, evolutionary relationships to other wolves and short life-history and ecology. Zhùr is the most complete wolf mummy known. She lived approximately 57,000 years ago and died in her den during a collapse of the sediments. During her short life, she ate aquatic resources, and is related to ancient Beringian and Russian gray wolves and her clade is basal to all living gray wolves. VIDEO ABSTRACT.

Ancient divergence of Indian and Tibetan wolves revealed by recombination-aware phylogenomics

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000-700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome-wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high-recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

Considering Pleistocene North American wolves and coyotes in the eastern Canis origin story

The evolutionary origins and hybridization patterns of Canis species in North America have been hotly debated for the past 30 years. Disentangling ancestry and timing of hybridization in Great Lakes wolves, eastern Canadian wolves, red wolves, and eastern coyotes are most often partitioned into a 2-species model that assigns all ancestry to gray wolves and/or coyotes, and a 3-species model that includes a third, North American evolved eastern wolf genome. The proposed models address recent or sometimes late Holocene hybridization events but have largely ignored potential Pleistocene era progenitors and opportunities for hybridization that may have impacted the current mixed genomes in eastern Canada and the United States. Here, we re-analyze contemporary and ancient mitochondrial DNA genomes with Bayesian phylogenetic analyses to more accurately estimate divergence dates among lineages. We combine that with a review of the literature on Late Pleistocene Canis distributions to: (a) identify potential Pleistocene progenitors to southern North American gray wolves and eastern wolves; and (b) illuminate opportunities for ancient hybridization events. Specifically, we propose that Beringian gray wolves (C. lupus) and extinct large wolf-like coyotes (C. latrans orcutti) are likely progenitors to Mexican and Plains gray wolves and eastern wolves, respectively, and may represent a potentially unrecognized source of introgressed genomic variation within contemporary Canis genomes. These events speak to the potential origins of contemporary genomes and provide a new perspective on Canis ancestry, but do not negate current conservation priorities of dwindling wolf populations with unique genomic signatures and key ecologically critical roles.

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.

Dog domestication and the dual dispersal of people and dogs into the Americas

Advances in the isolation and sequencing of ancient DNA have begun to reveal the population histories of both people and dogs. Over the last 10,000 y, the genetic signatures of ancient dog remains have been linked with known human dispersals in regions such as the Arctic and the remote Pacific. It is suspected, however, that this relationship has a much deeper antiquity, and that the tandem movement of people and dogs may have begun soon after the domestication of the dog from a gray wolf ancestor in the late Pleistocene. Here, by comparing population genetic results of humans and dogs from Siberia, Beringia, and North America, we show that there is a close correlation in the movement and divergences of their respective lineages. This evidence places constraints on when and where dog domestication took place. Most significantly, it suggests that dogs were domesticated in Siberia by ∼23,000 y ago, possibly while both people and wolves were isolated during the harsh climate of the Last Glacial Maximum. Dogs then accompanied the first people into the Americas and traveled with them as humans rapidly dispersed into the continent beginning ∼15,000 y ago.

Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge

The Bering Land Bridge (BLB) last connected Eurasia and North America during the Late Pleistocene. Although the BLB would have enabled transfers of terrestrial biota in both directions, it also acted as an ecological filter whose permeability varied considerably over time. Here we explore the possible impacts of this ecological corridor on genetic diversity within, and connectivity among, populations of a once wide-ranging group, the caballine horses (Equus spp.). Using a panel of 187 mitochondrial and eight nuclear genomes recovered from present-day and extinct caballine horses sampled across the Holarctic, we found that Eurasian horse populations initially diverged from those in North America, their ancestral continent, around 1.0-0.8 million years ago. Subsequent to this split our mitochondrial DNA analysis identified two bidirectional long-range dispersals across the BLB ~875-625 and ~200-50 thousand years ago, during the Middle and Late Pleistocene. Whole genome analysis indicated low levels of gene flow between North American and Eurasian horse populations, which probably occurred as a result of these inferred dispersals. Nonetheless, mitochondrial and nuclear diversity of caballine horse populations retained strong phylogeographical structuring. Our results suggest that barriers to gene flow, currently unidentified but possibly related to habitat distribution across Beringia or ongoing evolutionary divergence, played an important role in shaping the early genetic history of caballine horses, including the ancestors of living horses within Equus ferus.

A refined proposal for the origin of dogs: the case study of Gnirshöhle, a Magdalenian cave site

Dogs are known to be the oldest animals domesticated by humans. Although many studies have examined wolf domestication, the geographic and temporal origin of this process is still being debated. To address this issue, our study sheds new light on the early stages of wolf domestication during the Magdalenian period (16–14 ka cal BP) in the Hegau Jura region (Southwestern Germany and Switzerland). By combining morphology, genetics, and isotopes, our multidisciplinary approach helps to evaluate alternate processes driving the early phases of domestication. The isotope analysis uncovered a restricted, low δ¹⁵N protein diet for all analyzed Gnirshöhle specimens, while morphological examinations and phylogenetic relationships did not unequivocally assign them to one or the other canid lineage. Intriguingly, the newly generated mitochondrial canid genomes span the entire genetic diversity of modern dogs and wolves. Such high mitochondrial diversity could imply that Magdalenian people tamed and reared animals originating from different wolf lineages. We discuss our results in light of three ecological hypotheses and conclude that both domestication and the existence of a specialized wolf ecomorph are highly probable. However, due to their proximity to humans and a restricted diet, we propose domestication as the most likely scenario explaining the patterns observed herein.

Extended survival of Pleistocene Siberian wolves into the early 20th century on the island of Honshū

The Japanese or Honshū wolf was one the most distinct gray wolf subspecies due to its small stature and endemicity to the islands of Honshū, Shikoku, and Kyūshū. Long revered as a guardian of farmers and travellers, it was persecuted from the 17th century following a rabies epidemic, which led to its extinction in the early 20th century. To better understand its evolutionary history, we sequenced the nuclear genome of a 19th century Honshū wolf specimen to an average depth of coverage of 3.7✕. We find Honshū wolves were closely related to a lineage of Siberian wolves that were previously believed to have gone extinct in the Late Pleistocene, thereby extending the survival of this ancient lineage until the early 20th century. We also detected significant gene flow between Japanese dogs and the Honshū wolf, corroborating previous reports on Honshū wolf dog interbreeding.

•

Generated 3.7✕ nuclear genome of the extinct Honshu wolf

•

The Honshū wolf belonged to the lineage of Siberian Pleistocene wolves

•

There was gene flow between Honshū wolves and Japanese dogs

Examining Natural History through the Lens of Palaeogenomics

The many high-resolution tools that are uniquely applicable to specimens from the Quaternary period (the past ~2.5 Ma) provide an opportunity to cross-validate data and test hypotheses based on the morphology and distribution of fossils. Among these tools is palaeogenomics - the genome-scale sequencing of genetic material from ancient specimens - that can provide direct insight into ecology and evolution, potentially improving the accuracy of inferences about past ecological communities over longer timescales. Palaeogenomics has revealed instances of over- and underestimation of extinct diversity, detected cryptic faunal migration and turnover, allowed quantification of widespread sex biases and sexual dimorphism in the fossil record, revealed past hybridisation events and hybrid individuals, and has highlighted previously unrecognised routes of zoonotic disease transfer.

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.

Dog colour patterns explained by modular promoters of ancient canid origin

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records. This article is part of the theme issue 'Cross-disciplinary approaches to prehistoric demography'.

Population Dynamics in Italian Canids between the Late Pleistocene and Bronze Age

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a ~24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

Pre-extinction Demographic Stability and Genomic Signatures of Adaptation in the Woolly Rhinoceros

Ancient DNA has significantly improved our understanding of the evolution and population history of extinct megafauna. However, few studies have used complete ancient genomes to examine species responses to climate change prior to extinction. The woolly rhinoceros (Coelodonta antiquitatis) was a cold-adapted megaherbivore widely distributed across northern Eurasia during the Late Pleistocene and became extinct approximately 14 thousand years before present (ka BP). While humans and climate change have been proposed as potential causes of extinction [1-3], knowledge is limited on how the woolly rhinoceros was impacted by human arrival and climatic fluctuations [2]. Here, we use one complete nuclear genome and 14 mitogenomes to investigate the demographic history of woolly rhinoceros leading up to its extinction. Unlike other northern megafauna, the effective population size of woolly rhinoceros likely increased at 29.7 ka BP and subsequently remained stable until close to the species' extinction. Analysis of the nuclear genome from a ∼18.5-ka-old specimen did not indicate any increased inbreeding or reduced genetic diversity, suggesting that the population size remained steady for more than 13 ka following the arrival of humans [4]. The population contraction leading to extinction of the woolly rhinoceros may have thus been sudden and mostly driven by rapid warming in the Bølling-Allerød interstadial. Furthermore, we identify woolly rhinoceros-specific adaptations to arctic climate, similar to those of the woolly mammoth. This study highlights how species respond differently to climatic fluctuations and further illustrates the potential of palaeogenomics to study the evolutionary history of extinct species.