Species-specific responses of Late Quaternary megafauna to climate and humans.

Despite decades of research, the roles of climate and humans in driving the dramatic extinctions of large-bodied mammals during the Late Quaternary period remain contentious. Here we use ancient DNA, species distribution models and the human fossil record to elucidate how climate and humans shaped the demographic history of woolly rhinoceros, woolly mammoth, wild horse, reindeer, bison and musk ox. We show that climate has been a major driver of population change over the past 50,000 years. However, each species responds differently to the effects of climatic shifts, habitat redistribution and human encroachment. Although climate change alone can explain the extinction of some species, such as Eurasian musk ox and woolly rhinoceros, a combination of climatic and anthropogenic effects appears to be responsible for the extinction of others, including Eurasian steppe bison and wild horse. We find no genetic signature or any distinctive range dynamics distinguishing extinct from surviving species, emphasizing the challenges associated with predicting future responses of extant mammals to climate and human-mediated habitat change.

Ancient human genome sequence of an extinct Palaeo-Eskimo.

We report here the genome sequence of an ancient human. Obtained from approximately 4,000-year-old permafrost-preserved hair, the genome represents a male individual from the first known culture to settle in Greenland. Sequenced to an average depth of 20x, we recover 79% of the diploid genome, an amount close to the practical limit of current sequencing technologies. We identify 353,151 high-confidence single-nucleotide polymorphisms (SNPs), of which 6.8% have not been reported previously. We estimate raw read contamination to be no higher than 0.8%. We use functional SNP assessment to assign possible phenotypic characteristics of the individual that belonged to a culture whose location has yielded only trace human remains. We compare the high-confidence SNPs to those of contemporary populations to find the populations most closely related to the individual. This provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of that giving rise to the modern Native Americans and Inuit.

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).

Recharacterization of ancient DNA miscoding lesions: insights in the era of sequencing-by-synthesis.

Although ancient DNA (aDNA) miscoding lesions have been studied since the earliest days of the field, their nature remains a source of debate. A variety of conflicting hypotheses exist about which miscoding lesions constitute true aDNA damage as opposed to PCR polymerase amplification error. Furthermore, considerable disagreement and speculation exists on which specific damage events underlie observed miscoding lesions. The root of the problem is that it has previously been difficult to assemble sufficient data to test the hypotheses, and near-impossible to accurately determine the specific strand of origin of observed damage events. With the advent of emulsion-based clonal amplification (emPCR) and the sequencing-by-synthesis technology this has changed. In this paper we demonstrate how data produced on the Roche GS20 genome sequencer can determine miscoding lesion strands of origin, and subsequently be interpreted to enable characterization of the aDNA damage behind the observed phenotypes. Through comparative analyses on 390,965 bp of modern chloroplast and 131,474 bp of ancient woolly mammoth GS20 sequence data we conclusively demonstrate that in this sample at least, a permafrost preserved specimen, Type 2 (cytosine-->thymine/guanine-->adenine) miscoding lesions represent the overwhelming majority of damage-derived miscoding lesions. Additionally, we show that an as yet unidentified guanine-->adenine analogue modification, not the conventionally argued cytosine-->uracil deamination, underpins a significant proportion of Type 2 damage. How widespread these implications are for aDNA will become apparent as future studies analyse data recovered from a wider range of substrates.

Assessing the fidelity of ancient DNA sequences amplified from nuclear genes.

To date, the field of ancient DNA has relied almost exclusively on mitochondrial DNA (mtDNA) sequences. However, a number of recent studies have reported the successful recovery of ancient nuclear DNA (nuDNA) sequences, thereby allowing the characterization of genetic loci directly involved in phenotypic traits of extinct taxa. It is well documented that postmortem damage in ancient mtDNA can lead to the generation of artifactual sequences. However, as yet no one has thoroughly investigated the damage spectrum in ancient nuDNA. By comparing clone sequences from 23 fossil specimens, recovered from environments ranging from permafrost to desert, we demonstrate the presence of miscoding lesion damage in both the mtDNA and nuDNA, resulting in insertion of erroneous bases during amplification. Interestingly, no significant differences in the frequency of miscoding lesion damage are recorded between mtDNA and nuDNA despite great differences in cellular copy numbers. For both mtDNA and nuDNA, we find significant positive correlations between total sequence heterogeneity and the rates of type 1 transitions (adenine --> guanine and thymine --> cytosine) and type 2 transitions (cytosine --> thymine and guanine --> adenine), respectively. Type 2 transitions are by far the most dominant and increase relative to those of type 1 with damage load. The results suggest that the deamination of cytosine (and 5-methyl cytosine) to uracil (and thymine) is the main cause of miscoding lesions in both ancient mtDNA and nuDNA sequences. We argue that the problems presented by postmortem damage, as well as problems with contamination from exogenous sources of conserved nuclear genes, allelic variation, and the reliance on single nucleotide polymorphisms, call for great caution in studies relying on ancient nuDNA sequences.

Crosslinks rather than strand breaks determine access to ancient DNA sequences from frozen sediments.

Diagenesis was studied in DNA obtained from Siberian permafrost (permanently frozen soil) ranging from 10,000 to 400,000 years in age. Despite optimal preservation conditions, we found the sedimentary DNA to be severely modified by interstrand crosslinks; single- and double-stranded breaks; and freely exposed sugar, phosphate, and hydroxyl groups. Intriguingly, interstrand crosslinks were found to accumulate approximately 100 times faster than single-stranded breaks, suggesting that crosslinking rather than depurination is the primary limiting factor for ancient DNA amplification under frozen conditions. The results question the reliability of the commonly used models relying on depurination kinetics for predicting the long-term survival of DNA under permafrost conditions and suggest that new strategies for repair of ancient DNA must be considered if the yield of amplifiable DNA from permafrost sediments is to be significantly increased. Using the obtained rate constant for interstrand crosslinks the maximal survival time of amplifiable 120-bp fragments of bacterial 16S ribosomal DNA was estimated to be approximately 400,000 years. Additionally, a clear relationship was found between DNA damage and sample age, contradicting previously raised concerns about the possible leaching of free DNA molecules between permafrost layers.

The last Viking King: a royal maternity case solved by ancient DNA analysis.

The last of the Danish Viking Kings, Sven Estridsen, died in a.d. 1074 and is entombed in Roskilde Cathedral with other Danish kings and queens. Sven's mother, Estrid, is entombed in a pillar across the chancel. However, while there is no reasonable doubt about the identity of Sven, there have been doubts among historians whether the woman entombed was indeed Estrid. To shed light on this problem, we have extracted and analysed mitochondrial DNA (mtDNA) from pulp of teeth from each of the two royals. Four overlapping DNA-fragments covering about 400bp of hypervariable region 1 (HVR-1) of the D-loop were PCR amplified, cloned and a number of clones with each segment were sequenced. Also a segment containing the H/non-H specific nucleotide 7028 was sequenced. Consensus sequences were determined and D-loop results were replicated in an independent laboratory. This allowed the assignment of King Sven Estridsen to haplogroup H; Estrid's sequence differed from that of Sven at two positions in HVR-1, 16093T-->C and 16304T-->C, indicating that she belongs to subgroup H5a. Given the maternal inheritance of mtDNA, offspring will have the same mtDNA sequence as their mother with the exception of rare cases where the sequence has been altered by a germ line mutation. Therefore, the observation of two sequence differences makes it highly unlikely that the entombed woman was the mother of Sven. In addition, physical examination of the skeleton and the teeth strongly indicated that this woman was much younger (approximately 35 years) at the time of death than the 70 years history records tell. Although the entombed woman cannot be the Estrid, she may well be one of Sven's two daughters-in-law who were also called Estrid and who both became queens.

DNA from pre-Clovis human coprolites in Oregon, North America.

The timing of the first human migration into the Americas and its relation to the appearance of the Clovis technological complex in North America at about 11,000 to 10,800 radiocarbon years before the present (14C years B.P.) remains contentious. We establish that humans were present at Paisley 5 Mile Point Caves, in south-central Oregon, by 12,300 14C years B.P., through the recovery of human mitochondrial DNA (mtDNA) from coprolites, directly dated by accelerator mass spectrometry. The mtDNA corresponds to Native American founding haplogroups A2 and B2. The dates of the coprolites are >1000 14C years earlier than currently accepted dates for the Clovis complex.

800,000 year old mammoth DNA, modern elephant DNA or PCR artefact?

Poulakakis and colleagues (Poulakakis et al. 2006: Biol. Lett. 2, 451-454), report the recovery of 'authentic' mammoth DNA from an 800,000-year-old fragment of bone excavated on the island of Crete. In light of results from other ancient DNA studies that indicate how DNA survival is unlikely in samples, which are recovered from warm environments and are relatively old (e.g. more than 100,000 years), these findings come as a great surprise. Here, we show that problems exist with the methodological approaches used in the study. First, the nested PCR technique as reported is nonsensical--one of the second round 'nested' primers falls outside the amplicon of the first round PCR. More worryingly, the binding region of one of the first round primers (Elcytb320R) falls within the short 43 base pair reported mammoth sequence, specifically covering two of the three reportedly diagnostic Elephas polymorphisms. Finally, we demonstrate using a simple BLAST search in GenBank that the claimed 'uniquely derived character state' for mammoths is in fact also found within modern elephants.

The use of coded PCR primers enables high-throughput sequencing of multiple homolog amplification products by 454 parallel sequencing.

BACKGROUND: The invention of the Genome Sequence 20 DNA Sequencing System (454 parallel sequencing platform) has enabled the rapid and high-volume production of sequence data. Until now, however, individual emulsion PCR (emPCR) reactions and subsequent sequencing runs have been unable to combine template DNA from multiple individuals, as homologous sequences cannot be subsequently assigned to their original sources. METHODOLOGY: We use conventional PCR with 5'-nucleotide tagged primers to generate homologous DNA amplification products from multiple specimens, followed by sequencing through the high-throughput Genome Sequence 20 DNA Sequencing System (GS20, Roche/454 Life Sciences). Each DNA sequence is subsequently traced back to its individual source through 5'tag-analysis. CONCLUSIONS: We demonstrate that this new approach enables the assignment of virtually all the generated DNA sequences to the correct source once sequencing anomalies are accounted for (miss-assignment rate<0.4%). Therefore, the method enables accurate sequencing and assignment of homologous DNA sequences from multiple sources in single high-throughput GS20 run. We observe a bias in the distribution of the differently tagged primers that is dependent on the 5' nucleotide of the tag. In particular, primers 5' labelled with a cytosine are heavily overrepresented among the final sequences, while those 5' labelled with a thymine are strongly underrepresented. A weaker bias also exists with regards to the distribution of the sequences as sorted by the second nucleotide of the dinucleotide tags. As the results are based on a single GS20 run, the general applicability of the approach requires confirmation. However, our experiments demonstrate that 5'primer tagging is a useful method in which the sequencing power of the GS20 can be applied to PCR-based assays of multiple homologous PCR products. The new approach will be of value to a broad range of research areas, such as those of comparative genomics, complete mitochondrial analyses, population genetics, and phylogenetics.

Diverse plant and animal genetic records from Holocene and Pleistocene sediments.

Genetic analyses of permafrost and temperate sediments reveal that plant and animal DNA may be preserved for long periods, even in the absence of obvious macrofossils. In Siberia, five permafrost cores ranging from 400,000 to 10,000 years old contained at least 19 different plant taxa, including the oldest authenticated ancient DNA sequences known, and megafaunal sequences including mammoth, bison, and horse. The genetic data record a number of dramatic changes in the taxonomic diversity and composition of Beringian vegetation and fauna. Temperate cave sediments in New Zealand also yielded DNA sequences of extinct biota, including two species of ratite moa, and 29 plant taxa characteristic of the prehuman environment. Therefore, many sedimentary deposits may contain unique, and widespread, genetic records of paleoenvironments.

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.