Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Faunal record identifies Bering isthmus conditions as constraint to end-Pleistocene migration to the New World.

Human colonization of the New World is generally believed to have entailed migrations from Siberia across the Bering isthmus. However, the limited archaeological record of these migrations means that details of the timing, cause and rate remain cryptic. Here, we have used a combination of ancient DNA, 14C dating, hydrogen and oxygen isotopes, and collagen sequencing to explore the colonization history of one of the few other large mammals to have successfully migrated into the Americas at this time: the North American elk (Cervus elaphus canadensis), also known as wapiti. We identify a long-term occupation of northeast Siberia, far beyond the species's current Old World distribution. Migration into North America occurred at the end of the last glaciation, while the northeast Siberian source population became extinct only within the last 500 years. This finding is congruent with a similar proposed delay in human colonization, inferred from modern human mitochondrial DNA, and suggestions that the Bering isthmus was not traversable during parts of the Late Pleistocene. Our data imply a fundamental constraint in crossing Beringia, placing limits on the age and mode of human settlement in the Americas, and further establish the utility of ancient DNA in palaeontological investigations of species histories.

Sequencing the nuclear genome of the extinct woolly mammoth.

In 1994, two independent groups extracted DNA from several Pleistocene epoch mammoths and noted differences among individual specimens. Subsequently, DNA sequences have been published for a number of extinct species. However, such ancient DNA is often fragmented and damaged, and studies to date have typically focused on short mitochondrial sequences, never yielding more than a fraction of a per cent of any nuclear genome. Here we describe 4.17 billion bases (Gb) of sequence from several mammoth specimens, 3.3 billion (80%) of which are from the woolly mammoth (Mammuthus primigenius) genome and thus comprise an extensive set of genome-wide sequence from an extinct species. Our data support earlier reports that elephantid genomes exceed 4 Gb. The estimated divergence rate between mammoth and African elephant is half of that between human and chimpanzee. The observed number of nucleotide differences between two particular mammoths was approximately one-eighth of that between one of them and the African elephant, corresponding to a separation between the mammoths of 1.5-2.0 Myr. The estimated probability that orthologous elephant and mammoth amino acids differ is 0.002, corresponding to about one residue per protein. Differences were discovered between mammoth and African elephant in amino-acid positions that are otherwise invariant over several billion years of combined mammalian evolution. This study shows that nuclear genome sequencing of extinct species can reveal population differences not evident from the fossil record, and perhaps even discover genetic factors that affect extinction.

Holarctic genetic structure and range dynamics in the woolly mammoth.

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

Ancient DNA analyses exclude humans as the driving force behind late Pleistocene musk ox (Ovibos moschatus) population dynamics.

The causes of the late Pleistocene megafaunal extinctions are poorly understood. Different lines of evidence point to climate change, the arrival of humans, or a combination of these events as the trigger. Although many species went extinct, others, such as caribou and bison, survived to the present. The musk ox has an intermediate story: relatively abundant during the Pleistocene, it is now restricted to Greenland and the Arctic Archipelago. In this study, we use ancient DNA sequences, temporally unbiased summary statistics, and Bayesian analytical techniques to infer musk ox population dynamics throughout the late Pleistocene and Holocene. Our results reveal that musk ox genetic diversity was much higher during the Pleistocene than at present, and has undergone several expansions and contractions over the past 60,000 years. Northeast Siberia was of key importance, as it was the geographic origin of all samples studied and held a large diverse population until local extinction at approximately 45,000 radiocarbon years before present ((14)C YBP). Subsequently, musk ox genetic diversity reincreased at ca. 30,000 (14)C YBP, recontracted at ca. 18,000 (14)C YBP, and finally recovered in the middle Holocene. The arrival of humans into relevant areas of the musk ox range did not affect their mitochondrial diversity, and both musk ox and humans expanded into Greenland concomitantly. Thus, their population dynamics are better explained by a nonanthropogenic cause (for example, environmental change), a hypothesis supported by historic observations on the sensitivity of the species to both climatic warming and fluctuations.

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

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

Genetic structure and extinction of the woolly mammoth, Mammuthus primigenius.

The interval since circa 50 Ka has been a period of significant species extinctions among the large mammal fauna. However, the relative roles of an increasing human presence and a synchronous series of complex environmental changes in these extinctions have yet to be fully resolved. Recent analyses of fossil material from Beringia have clarified our understanding of the spatiotemporal pattern of Late Pleistocene extinctions, identifying periods of population turnover well before the last glacial maximum (LGM: circa 21 Ka) or subsequent human expansion. To examine the role of pre-LGM population changes in shaping the genetic structure of an extinct species, we analyzed the mitochondrial DNA of woolly mammoths in western Beringia and across its range. We identify genetic signatures of a range expansion of mammoths, from eastern to western Beringia, after the last interglacial (circa 125 Ka), and then an extended period during which demographic inference indicates no population-size increase. The most marked change in diversity at this time is the loss of one of two major mitochondrial lineages.

Ancient DNA sequences point to a large loss of mitochondrial genetic diversity in the saiga antelope (Saiga tatarica) since the Pleistocene.

Prior to the Holocene, the range of the saiga antelope (Saiga tatarica) spanned from France to the Northwest Territories of Canada. Although its distribution subsequently contracted to the steppes of Central Asia, historical records indicate that it remained extremely abundant until the end of the Soviet Union, after which its populations were reduced by over 95%. We have analysed the mitochondrial control region sequence variation of 27 ancient and 38 modern specimens, to assay how the species' genetic diversity has changed since the Pleistocene. Phylogenetic analyses reveal the existence of two well-supported, and clearly distinct, clades of saiga. The first, spanning a time range from >49,500 (14) C ybp to the present, comprises all the modern specimens and ancient samples from the Northern Urals, Middle Urals and Northeast Yakutia. The second clade is exclusive to the Northern Urals and includes samples dating from between 40,400 to 10,250 (14) C ybp. Current genetic diversity is much lower than that present during the Pleistocene, an observation that data modelling using serial coalescent indicates cannot be explained by genetic drift in a population of constant size. Approximate Bayesian Computation analyses show the observed data is more compatible with a drastic population size reduction (c. 66-77%) following either a demographic bottleneck in the course of the Holocene or late Pleistocene, or a geographic fragmentation (followed by local extinction of one subpopulation) at the Holocene/Pleistocene transition.

Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution.

BACKGROUND: The scientific literature contains many examples where DNA sequence analyses have been used to provide definitive answers to phylogenetic problems that traditional (non-DNA based) approaches alone have failed to resolve. One notable example concerns the rhinoceroses, a group for which several contradictory phylogenies were proposed on the basis of morphology, then apparently resolved using mitochondrial DNA fragments. RESULTS: In this study we report the first complete mitochondrial genome sequences of the extinct ice-age woolly rhinoceros (Coelodonta antiquitatis), and the threatened Javan (Rhinoceros sondaicus), Sumatran (Dicerorhinus sumatrensis), and black (Diceros bicornis) rhinoceroses. In combination with the previously published mitochondrial genomes of the white (Ceratotherium simum) and Indian (Rhinoceros unicornis) rhinoceroses, this data set putatively enables reconstruction of the rhinoceros phylogeny. While the six species cluster into three strongly supported sister-pairings: (i) The black/white, (ii) the woolly/Sumatran, and (iii) the Javan/Indian, resolution of the higher-level relationships has no statistical support. The phylogenetic signal from individual genes is highly diffuse, with mixed topological support from different genes. Furthermore, the choice of outgroup (horse vs tapir) has considerable effect on reconstruction of the phylogeny. The lack of resolution is suggestive of a hard polytomy at the base of crown-group Rhinocerotidae, and this is supported by an investigation of the relative branch lengths. CONCLUSION: Satisfactory resolution of the rhinoceros phylogeny may not be achievable without additional analyses of substantial amounts of nuclear DNA. This study provides a compelling demonstration that, in spite of substantial sequence length, there are significant limitations with single-locus phylogenetics. We expect further examples of this to appear as next-generation, large-scale sequencing of complete mitochondrial genomes becomes commonplace in evolutionary studies. "The human factor in classification is nowhere more evident than in dealing with this superfamily (Rhinocerotoidea)." G. G. Simpson (1945).

First DNA sequences from Asian cave bear fossils reveal deep divergences and complex phylogeographic patterns.

Until recently, cave bears were believed to have only inhabited Europe. However, recent morphological evidence suggests that cave bears' geographic range extended as far east as Transbaikalia, Eastern Siberia. These Asian cave bears were morphologically distinct from European cave bears. However, how they related to European lineages remains unclear, stressing the need to assess the phylogenetic and phylogeographic relationship between Asian cave bears and their European relatives. In this work, we address this issue using a 227 base-pair fragment of the mitochondrial control region obtained from nine fossil bone samples from eight sites from the Urals, Caucasus, Altai Mountains, Ukraine and Yana River region in Eastern Siberia. Results of the phylogenetic analyses indicate that (i) the cave bear from the Yana River is most closely related to cave bears from the Caucasus region; (ii) the Caucasus/Yana group of bears is genetically very distinct from both European cave bears and brown bears, suggesting that these bears could represent an independent species; and (iii) the Western European cave bear lineage reached at least temporarily to the Altai Mountains, 7000 km east of their known centre of distribution. These results suggest that the diversity of cave bears was greater than previously believed, and that they could survive in a much wider range of ecological conditions than previously assumed. They also agree with recent studies on other extinct and extant species, such as wolves, hyenas and steppe bison, which have also revealed higher genetic and ecological diversity in Pleistocene populations than previously known.

Phylogeography of lions (Panthera leo ssp.) reveals three distinct taxa and a late Pleistocene reduction in genetic diversity.

Lions were the most widespread carnivores in the late Pleistocene, ranging from southern Africa to the southern USA, but little is known about the evolutionary relationships among these Pleistocene populations or the dynamics that led to their extinction. Using ancient DNA techniques, we obtained mitochondrial sequences from 52 individuals sampled across the present and former range of lions. Phylogenetic analysis revealed three distinct clusters: (i) modern lions, Panthera leo; (ii) extinct Pleistocene cave lions, which formed a homogeneous population extending from Europe across Beringia (Siberia, Alaska and western Canada); and (iii) extinct American lions, which formed a separate population south of the Pleistocene ice sheets. The American lion appears to have become genetically isolated around 340 000 years ago, despite the apparent lack of significant barriers to gene flow with Beringian populations through much of the late Pleistocene. We found potential evidence of a severe population bottleneck in the cave lion during the previous interstadial, sometime after 48 000 years, adding to evidence from bison, mammoths, horses and brown bears that megafaunal populations underwent major genetic alterations throughout the last interstadial, potentially presaging the processes involved in the subsequent end-Pleistocene mass extinctions.

Evolution, systematics, and phylogeography of pleistocene horses in the new world: a molecular perspective.

The rich fossil record of horses has made them a classic example of evolutionary processes. However, while the overall picture of equid evolution is well known, the details are surprisingly poorly understood, especially for the later Pliocene and Pleistocene, c. 3 million to 0.01 million years (Ma) ago, and nowhere more so than in the Americas. There is no consensus on the number of equid species or even the number of lineages that existed in these continents. Likewise, the origin of the endemic South American genus Hippidion is unresolved, as is the phylogenetic position of the "stilt-legged" horses of North America. Using ancient DNA sequences, we show that, in contrast to current models based on morphology and a recent genetic study, Hippidion was phylogenetically close to the caballine (true) horses, with origins considerably more recent than the currently accepted date of c. 10 Ma. Furthermore, we show that stilt-legged horses, commonly regarded as Old World migrants related to the hemionid asses of Asia, were in fact an endemic North American lineage. Finally, our data suggest that there were fewer horse species in late Pleistocene North America than have been named on morphological grounds. Both caballine and stilt-legged lineages may each have comprised a single, wide-ranging species.

Non-destructive sampling of ancient insect DNA.

BACKGROUND: A major challenge for ancient DNA (aDNA) studies on insect remains is that sampling procedures involve at least partial destruction of the specimens. A recent extraction protocol reveals the possibility of obtaining DNA from past insect remains without causing visual morphological damage. We test the applicability of this protocol on historic museum beetle specimens dating back to AD 1820 and on ancient beetle chitin remains from permafrost (permanently frozen soil) dating back more than 47,000 years. Finally, we test the possibility of obtaining ancient insect DNA directly from non-frozen sediments deposited 3280-1800 years ago -- an alternative approach that also does not involve destruction of valuable material. METHODOLOGY/PRINCIPAL FINDINGS: The success of the methodological approaches are tested by PCR and sequencing of COI and 16S mitochondrial DNA (mtDNA) fragments of 77-204 base pairs (-bp) in size using species-specific and general insect primers. CONCLUSION/SIGNIFICANCE: The applied non-destructive DNA extraction method shows promising potential on insect museum specimens of historical age as far back as AD 1820, but less so on the ancient permafrost-preserved insect fossil remains tested, where DNA was obtained from samples up to ca. 26,000 years old. The non-frozen sediment DNA approach appears to have great potential for recording the former presence of insect taxa not normally preserved as macrofossils and opens new frontiers in research on ancient biodiversity.

Whole-genome shotgun sequencing of mitochondria from ancient hair shafts.

Although the application of sequencing-by-synthesis techniques to DNA extracted from bones has revolutionized the study of ancient DNA, it has been plagued by large fractions of contaminating environmental DNA. The genetic analyses of hair shafts could be a solution: We present 10 previously unexamined Siberian mammoth (Mammuthus primigenius) mitochondrial genomes, sequenced with up to 48-fold coverage. The observed levels of damage-derived sequencing errors were lower than those observed in previously published frozen bone samples, even though one of the specimens was >50,000 14C years old and another had been stored for 200 years at room temperature. The method therefore sets the stage for molecular-genetic analysis of museum collections.

Rise and fall of the Beringian steppe bison.

The widespread extinctions of large mammals at the end of the Pleistocene epoch have often been attributed to the depredations of humans; here we present genetic evidence that questions this assumption. We used ancient DNA and Bayesian techniques to reconstruct a detailed genetic history of bison throughout the late Pleistocene and Holocene epochs. Our analyses depict a large diverse population living throughout Beringia until around 37,000 years before the present, when the population's genetic diversity began to decline dramatically. The timing of this decline correlates with environmental changes associated with the onset of the last glacial cycle, whereas archaeological evidence does not support the presence of large populations of humans in Eastern Beringia until more than 15,000 years later.