Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments.

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.1 In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.2 Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia3-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,6-12 in paleodietary reconstructions,13-19 and facilitated genomic reconstructions from cave contexts.8,20-22 Murchie et al.6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.24 Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 40,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

Collapse of the mammoth-steppe in central Yukon as revealed by ancient environmental DNA.

The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

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.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears.

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

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.

Genomic Adaptations and Evolutionary History of the Extinct Scimitar-Toothed Cat, Homotherium latidens.

Homotherium was a genus of large-bodied scimitar-toothed cats, morphologically distinct from any extant felid species, that went extinct at the end of the Pleistocene [1-4]. They possessed large, saber-form serrated canine teeth, powerful forelimbs, a sloping back, and an enlarged optic bulb, all of which were key characteristics for predation on Pleistocene megafauna [5]. Previous mitochondrial DNA phylogenies suggested that it was a highly divergent sister lineage to all extant cat species [6-8]. However, mitochondrial phylogenies can be misled by hybridization [9], incomplete lineage sorting (ILS), or sex-biased dispersal patterns [10], which might be especially relevant for Homotherium since widespread mito-nuclear discrepancies have been uncovered in modern cats [10]. To examine the evolutionary history of Homotherium, we generated a ∼7x nuclear genome and a ∼38x exome from H. latidens using shotgun and target-capture sequencing approaches. Phylogenetic analyses reveal Homotherium as highly divergent (∼22.5 Ma) from living cat species, with no detectable signs of gene flow. Comparative genomic analyses found signatures of positive selection in several genes, including those involved in vision, cognitive function, and energy consumption, putatively consistent with diurnal activity, well-developed social behavior, and cursorial hunting [5]. Finally, we uncover relatively high levels of genetic diversity, suggesting that Homotherium may have been more abundant than the limited fossil record suggests [3, 4, 11-14]. Our findings complement and extend previous inferences from both the fossil record and initial molecular studies, enhancing our understanding of the evolution and ecology of this remarkable lineage.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations.

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

The evolutionary history of extinct and living lions.

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

Author Correction: Giant beaver palaeoecology inferred from stable isotopes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Giant beaver palaeoecology inferred from stable isotopes.

This is a multi-individual (n = 11), stable carbon and nitrogen isotope study of bone collagen (δ13Ccol and δ15Ncol) from the giant beaver (genus Castoroides). The now-extinct giant beaver was once one of the most widespread Pleistocene megafauna in North America. We confirm that Castoroides consumed a diet of predominantly submerged aquatic macrophytes. These dietary preferences rendered the giant beaver highly dependent on wetland habitat for survival. Castoroides' δ13Ccol and δ15Ncol do not support the hypothesis that the giant beaver consumed trees or woody plants, which suggests that it did not share the same behaviours as Castor (i.e., tree-cutting and harvesting). The onset of warmer, more arid conditions likely contributed to the extinction of Castoroides. Six new radiocarbon dates help establish the chronology of the northward dispersal of the giant beaver in Beringia, indicating a correlation with ice sheet retreat.

Nitrogen isotopes suggest a change in nitrogen dynamics between the Late Pleistocene and modern time in Yukon, Canada.

A magnificent repository of Late Pleistocene terrestrial megafauna fossils is contained in ice-rich loess deposits of Alaska and Yukon, collectively eastern Beringia. The stable carbon (δ13C) and nitrogen (δ15N) isotope compositions of bone collagen from these fossils are routinely used to determine paleodiet and reconstruct the paleoecosystem. This approach requires consideration of changes in C- and N-isotope dynamics over time and their effects on the terrestrial vegetation isotopic baseline. To test for such changes between the Late Pleistocene and modern time, we compared δ13C and δ15N for vegetation and bone collagen and structural carbonate of some modern, Yukon, arctic ground squirrels with vegetation and bones from Late Pleistocene fossil arctic ground squirrel nests preserved in Yukon loess deposits. The isotopic discrimination between arctic ground squirrel bone collagen and their diet was measured using modern samples, as were isotopic changes during plant decomposition; Over-wintering decomposition of typical vegetation following senescence resulted in a minor change (~0-1 ‰) in δ13C of modern Yukon grasses. A major change (~2-10 ‰) in δ15N was measured for decomposing Yukon grasses thinly covered by loess. As expected, the collagen-diet C-isotope discrimination measured for modern samples confirms that modern vegetation δ13C is a suitable proxy for the Late Pleistocene vegetation in Yukon Territory, after correction for the Suess effect. The N-isotope composition of vegetation from the fossil arctic ground squirrel nests, however, is determined to be ~2.8 ‰ higher than modern grasslands in the region, after correction for decomposition effects. This result suggests a change in N dynamics in this region between the Late Pleistocene and modern time.

A new genus of horse from Pleistocene North America.

The extinct 'New World stilt-legged', or NWSL, equids constitute a perplexing group of Pleistocene horses endemic to North America. Their slender distal limb bones resemble those of Asiatic asses, such as the Persian onager. Previous palaeogenetic studies, however, have suggested a closer relationship to caballine horses than to Asiatic asses. Here, we report complete mitochondrial and partial nuclear genomes from NWSL equids from across their geographic range. Although multiple NWSL equid species have been named, our palaeogenomic and morphometric analyses support the idea that there was only a single species of middle to late Pleistocene NWSL equid, and demonstrate that it falls outside of crown group Equus. We therefore propose a new genus, Haringtonhippus, for the sole species H. francisci. Our combined genomic and phenomic approach to resolving the systematics of extinct megafauna will allow for an improved understanding of the full extent of the terminal Pleistocene extinction event.

Evolutionary History of Saber-Toothed Cats Based on Ancient Mitogenomics.

Saber-toothed cats (Machairodontinae) are among the most widely recognized representatives of the now largely extinct Pleistocene megafauna. However, many aspects of their ecology, evolution, and extinction remain uncertain. Although ancient-DNA studies have led to huge advances in our knowledge of these aspects of many other megafauna species (e.g., mammoths and cave bears), relatively few ancient-DNA studies have focused on saber-toothed cats [1-3], and they have been restricted to short fragments of mitochondrial DNA. Here we investigate the evolutionary history of two lineages of saber-toothed cats (Smilodon and Homotherium) in relation to living carnivores and find that the Machairodontinae form a well-supported clade that is distinct from all living felids. We present partial mitochondrial genomes from one S. populator sample and three Homotherium sp. samples, including the only Late Pleistocene Homotherium sample from Eurasia [4]. We confirm the identification of the unique Late Pleistocene European fossil through ancient-DNA analyses, thus strengthening the evidence that Homotherium occurred in Europe over 200,000 years later than previously believed. This in turn forces a re-evaluation of its demography and extinction dynamics. Within the Machairodontinae, we find a deep divergence between Smilodon and Homotherium (∼18 million years) but limited diversity between the American and European Homotherium specimens. The genetic data support the hypothesis that all Late Pleistocene (or post-Villafrancian) Homotherium should be considered a single species, H. latidens, which was previously proposed based on morphological data [5, 6].

Nitrogen and carbon isotopic dynamics of subarctic soils and plants in southern Yukon Territory and its implications for paleoecological and paleodietary studies.

We examine here the carbon and nitrogen isotopic compositions of bulk soils (8 topsoil and 7 subsoils, including two soil profiles) and five different plant parts of 79 C3 plants from two main functional groups: herbs and shrubs/subshrubs, from 18 different locations in grasslands of southern Yukon Territory, Canada (eastern shoreline of Kluane Lake and Whitehorse area). The Kluane Lake region in particular has been identified previously as an analogue for Late Pleistocene eastern Beringia. All topsoils have higher average total nitrogen δ15N and organic carbon δ13C than plants from the same sites with a positive shift occurring with depth in two soil profiles analyzed. All plants analyzed have an average whole plant δ13C of -27.5 ± 1.2 ‰ and foliar δ13C of -28.0 ± 1.3 ‰, and average whole plant δ15N of -0.3 ± 2.2 ‰ and foliar δ15N of -0.6 ± 2.7 ‰. Plants analyzed here showed relatively smaller variability in δ13C than δ15N. Their average δ13C after suitable corrections for the Suess effect should be suitable as baseline for interpreting diets of Late Pleistocene herbivores that lived in eastern Beringia. Water availability, nitrogen availability, spacial differences and intra-plant variability are important controls on δ15N of herbaceous plants in the study area. The wider range of δ15N, the more numerous factors that affect nitrogen isotopic composition and their likely differences in the past, however, limit use of the modern N isotopic baseline for vegetation in paleodietary models for such ecosystems. That said, the positive correlation between foliar δ15N and N content shown for the modern plants could support use of plant δ15N as an index for plant N content and therefore forage quality. The modern N isotopic baseline cannot be applied directly to the past, but it is prerequisite to future efforts to detect shifts in N cycling and forage quality since the Late Pleistocene through comparison with fossil plants from the same region.

Megafaunal isotopes reveal role of increased moisture on rangeland during late Pleistocene extinctions.

The role of environmental change in the late Pleistocene megafaunal extinctions remains a key question, owing in part to uncertainty about landscape changes at continental scales. We investigated the influence of environmental changes on megaherbivores using bone collagen nitrogen isotopes (n = 684, 63 new) as a proxy for moisture levels in the rangelands that sustained late Pleistocene grazers. An increase in landscape moisture in Europe, Siberia and the Americas during the Last Glacial-Interglacial Transition (LGIT; ~25-10 kyr bp) directly affected megaherbivore ecology on four continents, and was associated with a key period of population decline and extinction. In all regions, the period of greatest moisture coincided with regional deglaciation and preceded the widespread formation of wetland environments. Moisture-driven environmental changes appear to have played an important part in the late Quaternary megafaunal extinctions through alteration of environments such as rangelands, which supported a large biomass of specialist grazers. On a continental scale, LGIT moisture changes manifested differently according to regional climate and geography, and the stable presence of grasslands surrounding the central forested belt of Africa during this period helps to explain why proportionally fewer African megafauna became extinct during the late Pleistocene.

Fossil and genomic evidence constrains the timing of bison arrival in North America.

The arrival of bison in North America marks one of the most successful large-mammal dispersals from Asia within the last million years, yet the timing and nature of this event remain poorly determined. Here, we used a combined paleontological and paleogenomic approach to provide a robust timeline for the entry and subsequent evolution of bison within North America. We characterized two fossil-rich localities in Canada's Yukon and identified the oldest well-constrained bison fossil in North America, a 130,000-y-old steppe bison, Bison cf. priscus We extracted and sequenced mitochondrial genomes from both this bison and from the remains of a recently discovered, ∼120,000-y-old giant long-horned bison, Bison latifrons, from Snowmass, Colorado. We analyzed these and 44 other bison mitogenomes with ages that span the Late Pleistocene, and identified two waves of bison dispersal into North America from Asia, the earliest of which occurred ∼195-135 thousand y ago and preceded the morphological diversification of North American bison, and the second of which occurred during the Late Pleistocene, ∼45-21 thousand y ago. This chronological arc establishes that bison first entered North America during the sea level lowstand accompanying marine isotope stage 6, rejecting earlier records of bison in North America. After their invasion, bison rapidly colonized North America during the last interglaciation, spreading from Alaska through continental North America; they have been continuously resident since then.

Bison phylogeography constrains dispersal and viability of the Ice Free Corridor in western Canada.

The Ice Free Corridor has been invoked as a route for Pleistocene human and animal dispersals between eastern Beringia and more southerly areas of North America. Despite the significance of the corridor, there are limited data for when and how this corridor was used. Hypothetical uses of the corridor include: the first expansion of humans from Beringia into the Americas, northward postglacial expansions of fluted point technologies into Beringia, and continued use of the corridor as a contact route between the north and south. Here, we use radiocarbon dates and ancient mitochondrial DNA from late Pleistocene bison fossils to determine the chronology for when the corridor was open and viable for biotic dispersals. The corridor was closed after ∼23,000 until 13,400 calendar years ago (cal y BP), after which we find the first evidence, to our knowledge, that bison used this route to disperse from the south, and by 13,000 y from the north. Our chronology supports a habitable and traversable corridor by at least 13,000 cal y BP, just before the first appearance of Clovis technology in interior North America, and indicates that the corridor would not have been available for significantly earlier southward human dispersal. Following the opening of the corridor, multiple dispersals of human groups between Beringia and interior North America may have continued throughout the latest Pleistocene and early Holocene. Our results highlight the utility of phylogeographic analyses to test hypotheses about paleoecological history and the viability of dispersal routes over time.

Solving the woolly mammoth conundrum: amino acid ¹5N-enrichment suggests a distinct forage or habitat.

Understanding woolly mammoth ecology is key to understanding Pleistocene community dynamics and evaluating the roles of human hunting and climate change in late Quaternary megafaunal extinctions. Previous isotopic studies of mammoths' diet and physiology have been hampered by the 'mammoth conundrum': woolly mammoths have anomalously high collagen δ(15)N values, which are more similar to coeval carnivores than herbivores, and which could imply a distinct diet and (or) habitat, or a physiological adaptation. We analyzed individual amino acids from collagen of adult woolly mammoths and coeval species, and discovered greater (15)N enrichment in source amino acids of woolly mammoths than in most other herbivores or carnivores. Woolly mammoths consumed an isotopically distinct food source, reflective of extreme aridity, dung fertilization, and (or) plant selection. This dietary signal suggests that woolly mammoths occupied a distinct habitat or forage niche relative to other Pleistocene herbivores.

Genomic Data from Extinct North American Camelops Revise Camel Evolutionary History.

Recent advances in paleogenomic technologies have enabled an increasingly detailed understanding of the evolutionary relationships of now-extinct mammalian taxa. However, a number of enigmatic Quaternary species have never been characterized with molecular data, often because available fossils are rare or are found in environments that are not optimal for DNA preservation. Here, we analyze paleogenomic data extracted from bones attributed to the late Pleistocene western camel, Camelops cf. hesternus, a species that was distributed across central and western North America until its extinction approximately 13,000 years ago. Despite a modal sequence length of only around 35 base pairs, we reconstructed high-coverage complete mitochondrial genomes and low-coverage partial nuclear genomes for each specimen. We find that Camelops is sister to African and Asian bactrian and dromedary camels, to the exclusion of South American camelids (llamas, guanacos, alpacas, and vicuñas). These results contradict previous morphology-based phylogenetic models for Camelops, which suggest instead a closer relationship between Camelops and the South American camelids. The molecular data imply a Late Miocene divergence of the Camelops clade from lineages that separately gave rise to the extant camels of Eurasia. Our results demonstrate the increasing capacity of modern paleogenomic methods to resolve evolutionary relationships among distantly related lineages.

American mastodon extirpation in the Arctic and Subarctic predates human colonization and terminal Pleistocene climate change.

Existing radiocarbon ((14)C) dates on American mastodon (Mammut americanum) fossils from eastern Beringia (Alaska and Yukon) have been interpreted as evidence they inhabited the Arctic and Subarctic during Pleistocene full-glacial times (∼ 18,000 (14)C years B.P.). However, this chronology is inconsistent with inferred habitat preferences of mastodons and correlative paleoecological evidence. To establish a last appearance date (LAD) for M. americanum regionally, we obtained 53 new (14)C dates on 36 fossils, including specimens with previously published dates. Using collagen ultrafiltration and single amino acid (hydroxyproline) methods, these specimens consistently date to beyond or near the ∼ 50,000 y B.P. limit of (14)C dating. Some erroneously "young" (14)C dates are due to contamination by exogenous carbon from natural sources and conservation treatments used in museums. We suggest mastodons inhabited the high latitudes only during warm intervals, particularly the Last Interglacial [Marine Isotope Stage (MIS) 5] when boreal forests existed regionally. Our (14)C dataset suggests that mastodons were extirpated from eastern Beringia during the MIS 4 glacial interval (∼ 75,000 y ago), following the ecological shift from boreal forest to steppe tundra. Mastodons thereafter became restricted to areas south of the continental ice sheets, where they suffered complete extinction ∼ 10,000 (14)C years B.P. Mastodons were already absent from eastern Beringia several tens of millennia before the first humans crossed the Bering Isthmus or the onset of climate changes during the terminal Pleistocene. Local extirpations of mastodons and other megafaunal populations in eastern Beringia were asynchrononous and independent of their final extinction south of the continental ice sheets.