Patterns of nucleotide misincorporations during enzymatic amplification and direct large-scale sequencing of ancient DNA

Unlocking the vault: next-generation museum population genomics

Natural history museum collections provide unique resources for understanding how species respond to environmental change, including the abrupt, anthropogenic climate change of the past century. Ideally, researchers would conduct genome-scale screening of museum specimens to explore the evolutionary consequences of environmental changes, but to date such analyses have been severely limited by the numerous challenges of working with the highly degraded DNA typical of historic samples. Here, we circumvent these challenges by using custom, multiplexed, exon capture to enrich and sequence ~11,000 exons (~4 Mb) from early 20th-century museum skins. We used this approach to test for changes in genomic diversity accompanying a climate-related range retraction in the alpine chipmunks (Tamias alpinus) in the high Sierra Nevada area of California, USA. We developed robust bioinformatic pipelines that rigorously detect and filter out base misincorporations in DNA derived from skins, most of which likely resulted from postmortem damage. Furthermore, to accommodate genotyping uncertainties associated with low-medium coverage data, we applied a recently developed probabilistic method to call single-nucleotide polymorphisms and estimate allele frequencies and the joint site frequency spectrum. Our results show increased genetic subdivision following range retraction, but no change in overall genetic diversity at either nonsynonymous or synonymous sites. This case study showcases the advantages of integrating emerging genomic and statistical tools in museum collection-based population genomic applications. Such technical advances greatly enhance the value of museum collections, even where a pre-existing reference is lacking and points to a broad range of potential applications in evolutionary and conservation biology.

Differentiation in neutral genes and a candidate gene in the pied flycatcher: using biological archives to track global climate change

Global climate change is one of the major driving forces for adaptive shifts in migration and breeding phenology and possibly impacts demographic changes if a species fails to adapt sufficiently. In Western Europe, pied flycatchers (Ficedula hypoleuca) have insufficiently adapted their breeding phenology to the ongoing advance of food peaks within their breeding area and consequently suffered local population declines. We address the question whether this population decline led to a loss of genetic variation, using two neutral marker sets (mitochondrial control region and microsatellites), and one potentially selectively non-neutral marker (avian Clock gene). We report temporal changes in genetic diversity in extant populations and biological archives over more than a century, using samples from sites differing in the extent of climate change. Comparing genetic differentiation over this period revealed that only the recent Dutch population, which underwent population declines, showed slightly lower genetic variation than the historic Dutch population. As that loss of variation was only moderate and not observed in all markers, current gene flow across Western and Central European populations might have compensated local loss of variation over the last decades. A comparison of genetic differentiation in neutral loci versus the Clock gene locus provided evidence for stabilizing selection. Furthermore, in all genetic markers, we found a greater genetic differentiation in space than in time. This pattern suggests that local adaptation or historic processes might have a stronger effect on the population structure and genetic variation in the pied flycatcher than recent global climate changes.

Evidence of coat color variation sheds new light on ancient canids

We have used a paleogenetics approach to investigate the genetic landscape of coat color variation in ancient Eurasian dog and wolf populations. We amplified DNA fragments of two genes controlling coat color, Mc1r (Melanocortin 1 Receptor) and CBD103 (canine-β-defensin), in respectively 15 and 19 ancient canids (dogs and wolf morphotypes) from 14 different archeological sites, throughout Asia and Europe spanning from ca. 12 000 B.P. (end of Upper Palaeolithic) to ca. 4000 B.P. (Bronze Age). We provide evidence of a new variant (R301C) of the Melanocortin 1 receptor (Mc1r) and highlight the presence of the beta-defensin melanistic mutation (CDB103-K locus) on ancient DNA from dog-and wolf-morphotype specimens. We show that the dominant K^B allele (CBD103), which causes melanism, and R301C (Mc1r), the variant that may cause light hair color, are present as early as the beginning of the Holocene, over 10 000 years ago. These results underline the genetic diversity of prehistoric dogs. This diversity may have partly stemmed not only from the wolf gene pool captured by domestication but also from mutations very likely linked to the relaxation of natural selection pressure occurring in-line with this process.

Extracting DNA from museum bird eggs, and whole genome amplification of archive DNA

We present a comprehensive protocol for extracting DNA from egg membranes and other internal debris recovered from the interior of blown museum bird eggs. A variety of commercially available DNA extraction methods were found to be applicable. DNA sequencing of polymerase chain reaction (PCR) products for a 176-bp fragment of mitochondrial DNA was successful for most egg samples (> 78%) even though the amount of DNA extracted (mean = 14.71 ± 4.55 ng/µL) was significantly less than that obtained for bird skin samples (mean = 67.88 ± 4.77 ng/µL). For PCR and sequencing of snipe (Gallinago) DNA, we provide eight new primers for the 'DNA barcode' region of COI mtDNA. In various combinations, the primers target a range of PCR products sized from 72 bp to the full 'barcode' of 751 bp. Not all possible combinations were tested with archive snipe DNA, but we found a significantly better success rate of PCR amplification for a shorter 176-bp target compared with a larger 288-bp fragment (67% vs. 39%). Finally, we explored the feasibility of whole genome amplification (WGA) for extending the use of archive DNA in PCR and sequencing applications. Of two WGA approaches, a PCR-based method was found to be able to amplify whole genomic DNA from archive skins and eggs from museum bird collections. After WGA, significantly more archive egg samples produced visible PCR products on agarose (56.9% before WGA vs. 79.0% after WGA). However, overall sequencing success did not improve significantly (78.8% compared with 83.0%).

Ancient DNA damage

Under favorable conditions DNA can survive for thousands of years in the remains of dead organisms. The DNA extracted from such remains is invariably degraded to a small average size by processes that at least partly involve depurination. It also contains large amounts of deaminated cytosine residues that are accumulated toward the ends of the molecules, as well as several other lesions that are less well characterized.

No ancient DNA damage in Actinobacteria from the Neanderthal bone

Background

The Neanderthal genome was recently sequenced using DNA extracted from a 38,000-year-old fossil. At the start of the project, the fraction of mammalian and bacterial DNA in the sample was estimated to be <6% and 9%, respectively. Treatment with restriction enzymes prior to sequencing increased the relative proportion of mammalian DNA to 15%, but the large majority of sequences remain uncharacterized.

Principal Findings

Our taxonomic profiling of 3.95 Gb of Neanderthal DNA isolated from the Vindija Neanderthal Vi33.16 fossil showed that 90% of about 50,000 rRNA gene sequence reads were of bacterial origin, of which Actinobacteria accounted for more than 75%. Actinobacteria also represented more than 80% of the PCR-amplified 16S rRNA gene sequences from a cave sediment sample taken from the same G layer as the Neanderthal bone. However, phylogenetic analyses did not identify any sediment clones that were closely related to the bone-derived sequences. We analysed the patterns of nucleotide differences in the individual sequence reads compared to the assembled consensus sequences of the rRNA gene sequences. The typical ancient nucleotide substitution pattern with a majority of C to T changes indicative of DNA damage was observed for the Neanderthal rRNA gene sequences, but not for the Streptomyces-like rRNA gene sequences.

Conclusions/Significance

Our analyses suggest that the Actinobacteria, and especially members of the Streptomycetales, contribute the majority of sequences in the DNA extracted from the Neanderthal fossil Vi33.16. The bacterial DNA showed no signs of damage, and we hypothesize that it was derived from bacteria that have been enriched inside the bone. The bioinformatic approach used here paves the way for future studies of microbial compositions and patterns of DNA damage in bacteria from archaeological bones. Such studies can help identify targeted measures to increase the relative amount of endogenous DNA in the sample.

Mitochondrial lineage sorting in action--historical biogeography of the Hyles euphorbiae complex (Sphingidae, Lepidoptera) in Italy

Background

Mitochondrial genes are among the most commonly used markers in studies of species’ phylogeography and to draw conclusions about taxonomy. The Hyles euphorbiae complex (HEC) comprises six distinct mitochondrial lineages in the Mediterranean region, of which one exhibits a cryptic disjunct distribution. The predominant mitochondrial lineage in most of Europe, euphorbiae, is also present on Malta; however, it is nowadays strangely absent from Southern Italy and Sicily, where it is replaced by 'italica'. A separate biological entity in Italy is further corroborated by larval colour patterns with a congruent, confined suture zone along the Northern Apennines. By means of historic DNA extracted from museum specimens, we aimed to investigate the evolution of the mitochondrial demographic structure of the HEC in Italy and Malta throughout the Twentieth Century.

Results

At the beginning of the Twentieth Century, the European mainland lineages were also present at a moderate frequency in Southern Italy and Sicily. The proportion of 'italica' then steadily increased in this area from below 60 percent to near fixation in about 120 years. Thus, geographical sorting of mitochondrial lineages in the HEC was not as complete then as the current demography suggests. The pattern of an integral 'italica' core region and a disjunct euphorbiae distribution evolved very recently. To explain these strong demographic changes, we propose genetic drift due to anthropogenic habitat loss and fragmentation in combination with an impact from recent climate warming that favoured the spreading of the potentially better adapted 'italica' populations.

Conclusions

The pattern of geographically separated mitochondrial lineages is commonly interpreted as representing long term separated entities. However, our results indicate that such a pattern can emerge surprisingly quickly, even in a widespread and rather common taxon. We thus caution against drawing hasty taxonomic conclusions from biogeographical patterns of mitochondrial markers derived from modern sampling alone.

Ancient DNA analysis affirms the canid from Altai as a primitive dog

The origin of domestic dogs remains controversial, with genetic data indicating a separation between modern dogs and wolves in the Late Pleistocene. However, only a few dog-like fossils are found prior to the Last Glacial Maximum, and it is widely accepted that the dog domestication predates the beginning of agriculture about 10,000 years ago. In order to evaluate the genetic relationship of one of the oldest dogs, we have isolated ancient DNA from the recently described putative 33,000-year old Pleistocene dog from Altai and analysed 413 nucleotides of the mitochondrial control region. Our analyses reveal that the unique haplotype of the Altai dog is more closely related to modern dogs and prehistoric New World canids than it is to contemporary wolves. Further genetic analyses of ancient canids may reveal a more exact date and centre of domestication.

On the accuracy of short read mapping

The development of high-throughput sequencing technologies has revolutionized the way we study genomes and gene regulation. In a single experiment, millions of reads are produced. To gain knowledge from these experiments the first thing to be done is finding the genomic origin of the reads, i.e., mapping the reads to a reference genome. In this new situation, conventional alignment tools are obsolete, as they cannot handle this huge amount of data in a reasonable amount of time. Thus, new mapping algorithms have been developed, which are fast at the expense of a small decrease in accuracy. In this chapter we discuss the current problems in short read mapping and show that mapping reads correctly is a nontrivial task. Through simple experiments with both real and synthetic data, we demonstrate that different mappers can give different results depending on the type of data, and that a considerable fraction of uniquely mapped reads is potentially mapped to an incorrect location. Furthermore, we provide simple statistical results on the expected number of random matches in a genome (E-value) and the probability of a random match as a function of read length. Finally, we show that quality scores contain valuable information for mapping and why mapping quality should be evaluated in a probabilistic manner. In the end, we discuss the potential of improving the performance of current methods by considering these quality scores in a probabilistic mapping program.

Molecular and phylogenetic approaches for assessing sources of Cryptosporidium contamination in water

The high sequence diversity and heterogeneity observed within species or genotypes of Cryptosporidium requires phylogenetic approaches for the identification of novel sequences obtained from the environment. A long-term study on Cryptosporidium in the agriculturally-intensive South Nation River watershed in Ontario, Canada was undertaken, in which 60 sequence types were detected. Of these sequence types 33 were considered novel with no identical matches in GenBank. Detailed phylogenetic analysis identified that most sequences belonged to 17 previously described species: Cryptosporidium andersoni, Cryptosporidium baileyi, Cryptosporidium hominis, Cryptosporidium parvum, Cryptosporidium ubiquitum, Cryptosporidium meleagridis, muskrat I, muskrat II, deer mouse II, fox, vole, skunk, shrew, W12, W18, W19 and W25 genotypes. In addition, two new genotypes were identified, W27 and W28. C. andersoni and the muskrat II genotype were most frequently detected in the water samples. Species associated with livestock made up 39% of the total molecular detections, while wildlife associated species and genotypes accounted for 55% of the Cryptosporidium identified. The human pathogenic species C. hominis and C. parvum had an overall prevalence of 1.6% in the environment, indicating a small risk to humans from the Cryptosporidium present in the watershed. Phylogenetic analysis and knowledge of host-parasite relationships are fundamental in using Cryptosporidium as a source-tracking or human health risk assessment tool.

Detection of ultra-rare mutations by next-generation sequencing

Next-generation DNA sequencing promises to revolutionize clinical medicine and basic research. However, while this technology has the capacity to generate hundreds of billions of nucleotides of DNA sequence in a single experiment, the error rate of ~1% results in hundreds of millions of sequencing mistakes. These scattered errors can be tolerated in some applications but become extremely problematic when "deep sequencing" genetically heterogeneous mixtures, such as tumors or mixed microbial populations. To overcome limitations in sequencing accuracy, we have developed a method termed Duplex Sequencing. This approach greatly reduces errors by independently tagging and sequencing each of the two strands of a DNA duplex. As the two strands are complementary, true mutations are found at the same position in both strands. In contrast, PCR or sequencing errors result in mutations in only one strand and can thus be discounted as technical error. We determine that Duplex Sequencing has a theoretical background error rate of less than one artifactual mutation per billion nucleotides sequenced. In addition, we establish that detection of mutations present in only one of the two strands of duplex DNA can be used to identify sites of DNA damage. We apply the method to directly assess the frequency and pattern of random mutations in mitochondrial DNA from human cells.

Improving ancient DNA read mapping against modern reference genomes

Background

Next-Generation Sequencing has revolutionized our approach to ancient DNA (aDNA) research, by providing complete genomic sequences of ancient individuals and extinct species. However, the recovery of genetic material from long-dead organisms is still complicated by a number of issues, including post-mortem DNA damage and high levels of environmental contamination. Together with error profiles specific to the type of sequencing platforms used, these specificities could limit our ability to map sequencing reads against modern reference genomes and therefore limit our ability to identify endogenous ancient reads, reducing the efficiency of shotgun sequencing aDNA.

Results

In this study, we compare different computational methods for improving the accuracy and sensitivity of aDNA sequence identification, based on shotgun sequencing reads recovered from Pleistocene horse extracts using Illumina GAIIx and Helicos Heliscope platforms. We show that the performance of the Burrows Wheeler Aligner (BWA), that has been developed for mapping of undamaged sequencing reads using platforms with low rates of indel-types of sequencing errors, can be employed at acceptable run-times by modifying default parameters in a platform-specific manner. We also examine if trimming likely damaged positions at read ends can increase the recovery of genuine aDNA fragments and if accurate identification of human contamination can be achieved using a strategy previously suggested based on best hit filtering. We show that combining our different mapping and filtering approaches can increase the number of high-quality endogenous hits recovered by up to 33%.

Conclusions

We have shown that Illumina and Helicos sequences recovered from aDNA extracts could not be aligned to modern reference genomes with the same efficiency unless mapping parameters are optimized for the specific types of errors generated by these platforms and by post-mortem DNA damage. Our findings have important implications for future aDNA research, as we define mapping guidelines that improve our ability to identify genuine aDNA sequences, which in turn could improve the genotyping accuracy of ancient specimens. Our framework provides a significant improvement to the standard procedures used for characterizing ancient genomes, which is challenged by contamination and often low amounts of DNA material.

Coprolites as a source of information on the genome and diet of the cave hyena

We performed high-throughput sequencing of DNA from fossilized faeces to evaluate this material as a source of information on the genome and diet of Pleistocene carnivores. We analysed coprolites derived from the extinct cave hyena (Crocuta crocuta spelaea), and sequenced 90 million DNA fragments from two specimens. The DNA reads enabled a reconstruction of the cave hyena mitochondrial genome with up to a 158-fold coverage. This genome, and those sequenced from extant spotted (Crocuta crocuta) and striped (Hyaena hyaena) hyena specimens, allows for the establishment of a robust phylogeny that supports a close relationship between the cave and the spotted hyena. We also demonstrate that high-throughput sequencing yields data for cave hyena multi-copy and single-copy nuclear genes, and that about 50 per cent of the coprolite DNA can be ascribed to this species. Analysing the data for additional species to indicate the cave hyena diet, we retrieved abundant sequences for the red deer (Cervus elaphus), and characterized its mitochondrial genome with up to a 3.8-fold coverage. In conclusion, we have demonstrated the presence of abundant ancient DNA in the coprolites surveyed. Shotgun sequencing of this material yielded a wealth of DNA sequences for a Pleistocene carnivore and allowed unbiased identification of diet.

Ancient DNA assessment of tiger salamander population in Yellowstone National Park

Recent data indicates that blotched tiger salamanders (Ambystoma tigrinum melanostictum) in northern regions of Yellowstone National Park are declining due to climate-related habitat changes. In this study, we used ancient and modern mitochondrial haplotype diversity to model the effective size of this amphibian population through recent geological time and to assess past responses to climatic changes in the region. Using subfossils collected from a cave in northern Yellowstone, we analyzed >700 base pairs of mitochondrial sequence from 16 samples ranging in age from 100 to 3300 years old and found that all shared an identical haplotype. Although mitochondrial diversity was extremely low within the living population, we still were able to detect geographic subdivision within the local area. Using serial coalescent modelling with Bayesian priors from both modern and ancient genetic data we simulated a range of probable population sizes and mutation rates through time. Our simulations suggest that regional mitochondrial diversity has remained relatively constant even through climatic fluctuations of recent millennia.

High-resolution analysis of cytosine methylation in ancient DNA

Epigenetic changes to gene expression can result in heritable phenotypic characteristics that are not encoded in the DNA itself, but rather by biochemical modifications to the DNA or associated chromatin proteins. Interposed between genes and environment, these epigenetic modifications can be influenced by environmental factors to affect phenotype for multiple generations. This raises the possibility that epigenetic states provide a substrate for natural selection, with the potential to participate in the rapid adaptation of species to changes in environment. Any direct test of this hypothesis would require the ability to measure epigenetic states over evolutionary timescales. Here we describe the first single-base resolution of cytosine methylation patterns in an ancient mammalian genome, by bisulphite allelic sequencing of loci from late Pleistocene Bison priscus remains. Retrotransposons and the differentially methylated regions of imprinted loci displayed methylation patterns identical to those derived from fresh bovine tissue, indicating that methylation patterns are preserved in the ancient DNA. Our findings establish the biochemical stability of methylated cytosines over extensive time frames, and provide the first direct evidence that cytosine methylation patterns are retained in DNA from ancient specimens. The ability to resolve cytosine methylation in ancient DNA provides a powerful means to study the role of epigenetics in evolution.

Case study: ancient DNA recovered from pleistocene-age remains of a Florida armadillo

Warm, humid regions are not ideal for long-term DNA preservation. Consequently, little ancient DNA research has been carried out involving taxa that lived in, for example, tropical and subtropical regions. Those studies that have isolated ancient DNA from warm environments have mostly been limited to the most recent several thousand years. Here, we discuss an ancient DNA experiment in which we attempt to amplify mitochondrial DNA from remains of armadillo, glyptodont, and pampathere from sites in Florida, USA, all believed to be around 10,000-12,000 years old. We were successful in recovering DNA from only one of these samples. However, based on the amount and distribution of DNA damage, the ancient DNA recovered was well-preserved despite the age and preservation environment. In this case study chapter, we discuss the experimental procedure we used to characterize the DNA from the Floridian samples, focusing on challenges of working with ancient specimens from warm environments and steps taken to confirm the authenticity of the recovered sequence.

Target enrichment via DNA hybridization capture

Recent advances in high-throughput DNA sequencing technologies have allowed entire nuclear genomes to be shotgun sequenced from ancient DNA (aDNA) extracts. Nonetheless, targeted analyses of specific genomic loci will remain an important tool for future aDNA studies. DNA capture via hybridization allows the efficient exploitation of current high-throughput sequencing for population genetic analyses using aDNA samples. Specifically, hybridization capture allows larger data sets to be generated for multiple target loci as well as for multiple samples in parallel. "Bait" molecules are used to select target regions from DNA libraries for sequencing. Here we present a brief overview of the currently available hybridization capture protocols using either an in-solution or a solid-phase (immobilized) approach. While it is possible to purchase ready-made kits for this purpose, I present a protocol that allows users to generate their own custom bait to be used for hybridization capture.

DNA damage in plant herbarium tissue

Dried plant herbarium specimens are potentially a valuable source of DNA. Efforts to obtain genetic information from this source are often hindered by an inability to obtain amplifiable DNA as herbarium DNA is typically highly degraded. DNA post-mortem damage may not only reduce the number of amplifiable template molecules, but may also lead to the generation of erroneous sequence information. A qualitative and quantitative assessment of DNA post-mortem damage is essential to determine the accuracy of molecular data from herbarium specimens. In this study we present an assessment of DNA damage as miscoding lesions in herbarium specimens using 454-sequencing of amplicons derived from plastid, mitochondrial, and nuclear DNA. In addition, we assess DNA degradation as a result of strand breaks and other types of polymerase non-bypassable damage by quantitative real-time PCR. Comparing four pairs of fresh and herbarium specimens of the same individuals we quantitatively assess post-mortem DNA damage, directly after specimen preparation, as well as after long-term herbarium storage. After specimen preparation we estimate the proportion of gene copy numbers of plastid, mitochondrial, and nuclear DNA to be 2.4–3.8% of fresh control DNA and 1.0–1.3% after long-term herbarium storage, indicating that nearly all DNA damage occurs on specimen preparation. In addition, there is no evidence of preferential degradation of organelle versus nuclear genomes. Increased levels of C→T/G→A transitions were observed in old herbarium plastid DNA, representing 21.8% of observed miscoding lesions. We interpret this type of post-mortem DNA damage-derived modification to have arisen from the hydrolytic deamination of cytosine during long-term herbarium storage. Our results suggest that reliable sequence data can be obtained from herbarium specimens.

Evaluating the impact of post-mortem damage in ancient DNA: a theoretical approach

The growth of ancient DNA research has offered exceptional opportunities and raised great expectations, but has also presented some considerable challenges. One of the ongoing issues is the impact of post-mortem damage in DNA molecules. Nucleotide alterations and DNA strand breakages lead to a significant decrease in the quantity of DNA molecules of useful length in a sample and to errors in the final DNA sequences obtained. We present a model of age-dependent DNA damage and quantify the influence of that damage on subsequent steps in the sequencing process, including the polymerase chain reaction and cloning. Calculations using our model show that deposition conditions, rather than the age of a sample, have the greatest influence on the level of DNA damage. In turn, this affects the probability of interpreting an erroneous (possessing damage-derived mutations) sequence as being authentic. We also evaluated the effect of post-mortem damage on real data sets using a Bayesian phylogenetic approach. According to our study, damage-derived sequence alterations appear to have little impact on the final DNA sequences. This indicates the effectiveness of current methods for sequence authentication and validation.

Metastatic colorectal cancer KRAS genotyping in routine practice: results and pitfalls

KRAS genotyping is mandatory before anti-epidermal growth factor receptor monoclonal antibody therapy in metastatic colorectal cancer, which is the second leading cause of cancer-related death in the United States and in Europe. Thus, large-scale KRAS mutation screening is needed for efficient patient management and in the future metastatic colorectal cancer genotyping might also include the detection of the BRAF V600E mutation, which is a very strong negative prognostic factor in colorectal cancer. We report our experience of routine KRAS/BRAF mutation screening practice performed on 1130 formalin-fixed paraffin-embedded tumor samples from 992 colorectal cancer patients. DNA was extracted from macrodissected tumor areas highlighted by a pathologist, KRAS codons 12/13 and BRAF V600E mutations were assessed in a single SNaPshot® multiplex assay and each mutation was confirmed by an independent analysis. KRAS and BRAF mutations were, respectively, present in 41.8 and 6.5% of the tumor samples. If KRAS and BRAF mutations were mutually exclusive, four samples presented two concomitant KRAS mutations. Genotyping of paired primary tumors and metastases from 44 patients indicated that 5 patients (11.4%) presented discordant KRAS mutational status. KRAS genotype heterogeneity was also observed within primary tumor sites in seven cases. Non-reproducible KRAS artefactual mutations were detected in 53 samples (4.7%). We found that the prominent mechanism leading to these artefactual mutations was the fragmentation of DNA occurring during tissue processing. Routine KRAS genotyping performed on formalin-fixed paraffin-embedded tissues requires, therefore, the development of quality control scheme for molecular pathology, especially because of DNA damages induced by formalin fixation. The tumor heterogeneity observed in some patients indicates that it should be more appropriate to perform KRAS genotyping on metastases if sample is available.

True single-molecule DNA sequencing of a pleistocene horse bone

Second-generation sequencing platforms have revolutionized the field of ancient DNA, opening access to complete genomes of past individuals and extinct species. However, these platforms are dependent on library construction and amplification steps that may result in sequences that do not reflect the original DNA template composition. This is particularly true for ancient DNA, where templates have undergone extensive damage post-mortem. Here, we report the results of the first "true single molecule sequencing" of ancient DNA. We generated 115.9 Mb and 76.9 Mb of DNA sequences from a permafrost-preserved Pleistocene horse bone using the Helicos HeliScope and Illumina GAIIx platforms, respectively. We find that the percentage of endogenous DNA sequences derived from the horse is higher among the Helicos data than Illumina data. This result indicates that the molecular biology tools used to generate sequencing libraries of ancient DNA molecules, as required for second-generation sequencing, introduce biases into the data that reduce the efficiency of the sequencing process and limit our ability to fully explore the molecular complexity of ancient DNA extracts. We demonstrate that simple modifications to the standard Helicos DNA template preparation protocol further increase the proportion of horse DNA for this sample by threefold. Comparison of Helicos-specific biases and sequence errors in modern DNA with those in ancient DNA also reveals extensive cytosine deamination damage at the 3' ends of ancient templates, indicating the presence of 3'-sequence overhangs. Our results suggest that paleogenomes could be sequenced in an unprecedented manner by combining current second- and third-generation sequencing approaches.

A 33,000-year-old incipient dog from the Altai Mountains of Siberia: evidence of the earliest domestication disrupted by the Last Glacial Maximum

Background

Virtually all well-documented remains of early domestic dog (Canis familiaris) come from the late Glacial and early Holocene periods (ca. 14,000–9000 calendar years ago, cal BP), with few putative dogs found prior to the Last Glacial Maximum (LGM, ca. 26,500–19,000 cal BP). The dearth of pre-LGM dog-like canids and incomplete state of their preservation has until now prevented an understanding of the morphological features of transitional forms between wild wolves and domesticated dogs in temporal perspective.

Methodology/Principal Finding

We describe the well-preserved remains of a dog-like canid from the Razboinichya Cave (Altai Mountains of southern Siberia). Because of the extraordinary preservation of the material, including skull, mandibles (both sides) and teeth, it was possible to conduct a complete morphological description and comparison with representative examples of pre-LGM wild wolves, modern wolves, prehistoric domesticated dogs, and early dog-like canids, using morphological criteria to distinguish between wolves and dogs. It was found that the Razboinichya Cave individual is most similar to fully domesticated dogs from Greenland (about 1000 years old), and unlike ancient and modern wolves, and putative dogs from Eliseevichi I site in central Russia. Direct AMS radiocarbon dating of the skull and mandible of the Razboinichya canid conducted in three independent laboratories resulted in highly compatible ages, with average value of ca. 33,000 cal BP.

Conclusions/Significance

The Razboinichya Cave specimen appears to be an incipient dog that did not give rise to late Glacial – early Holocene lineages and probably represents wolf domestication disrupted by the climatic and cultural changes associated with the LGM. The two earliest incipient dogs from Western Europe (Goyet, Belguim) and Siberia (Razboinichya), separated by thousands of kilometers, show that dog domestication was multiregional, and thus had no single place of origin (as some DNA data have suggested) and subsequent spread.

Molecular identification of bacteria by total sequence screening: determining the cause of death in ancient human subjects

Research of ancient pathogens in ancient human skeletons has been mainly carried out on the basis of one essential historical or archaeological observation, permitting specific pathogens to be targeted. Detection of ancient human pathogens without such evidence is more difficult, since the quantity and quality of ancient DNA, as well as the environmental bacteria potentially present in the sample, limit the analyses possible. Using human lung tissue and/or teeth samples from burials in eastern Siberia, dating from the end of 17(th) to the 19(th) century, we propose a methodology that includes the: 1) amplification of all 16S rDNA gene sequences present in each sample; 2) identification of all bacterial DNA sequences with a degree of identity ≥ 95%, according to quality criteria; 3) identification and confirmation of bacterial pathogens by the amplification of the rpoB gene; and 4) establishment of authenticity criteria for ancient DNA. This study demonstrates that from teeth samples originating from ancient human subjects, we can realise: 1) the correct identification of bacterial molecular sequence signatures by quality criteria; 2) the separation of environmental and pathogenic bacterial 16S rDNA sequences; 3) the distribution of bacterial species for each subject and for each burial; and 4) the characterisation of bacteria specific to the permafrost. Moreover, we identified three pathogens in different teeth samples by 16S rDNA sequence amplification: Bordetella sp., Streptococcus pneumoniae and Shigella dysenteriae. We tested for the presence of these pathogens by amplifying the rpoB gene. For the first time, we confirmed sequences from Bordetella pertussis in the lungs of an ancient male Siberian subject, whose grave dated from the end of the 17(th) century to the early 18(th) century.

Primers for identification of type and other archived specimens of Aphrodes leafhoppers (Hemiptera, Cicadellidae)

Primers were developed for leafhoppers of the genus Aphrodes amplifying 84-244 bp fragments of the mitochondrial cytochrome oxidase subunit I gene. DNA was extracted from legs of over 100-year-old archived museum specimens, amplified and sequenced. The fragments contain sufficient variation to unequivocally identify the different species. The majority of the analysed museum specimens, including three specimens of the syntype series for the UK endemic species A. aestuarina (Edwards 1908), were found to have been assigned to the wrong species. This work clearly underlines the need to validate museum specimens using molecular methods where identity is in doubt, based on reliable standards for species discrimination.

Influence of climate warming on arctic mammals? New insights from ancient DNA studies of the collared lemming Dicrostonyx torquatus

Background

Global temperature increased by approximately half a degree (Celsius) within the last 150 years. Even this moderate warming had major impacts on Earth's ecological and biological systems, especially in the Arctic where the magnitude of abiotic changes even exceeds those in temperate and tropical biomes. Therefore, understanding the biological consequences of climate change on high latitudes is of critical importance for future conservation of the species living in this habitat. The past 25,000 years can be used as a model for such changes, as they were marked by prominent climatic changes that influenced geographical distribution, demographic history and pattern of genetic variation of many extant species. We sequenced ancient and modern DNA of the collared lemming (Dicrostonyx torquatus), which is a key species of the arctic biota, from a single site (Pymva Shor, Northern Pre Urals, Russia) to see if climate warming events after the Last Glacial Maximum had detectable effects on the genetic variation of this arctic rodent species, which is strongly associated with a cold and dry climate.

Results

Using three dimensional network reconstructions we found a dramatic decline in genetic diversity following the LGM. Model-based approaches such as Approximate Bayesian Computation and Markov Chain Monte Carlo based Bayesian inference show that there is evidence for a population decline in the collared lemming following the LGM, with the population size dropping to a minimum during the Greenland Interstadial 1 (Bølling/Allerød) warming phase at 14.5 kyrs BP.

Conclusion

Our results show that previous climate warming events had a strong influence on genetic diversity and population size of collared lemmings. Due to its already severely compromised genetic diversity a similar population reduction as a result of the predicted future climate change could completely abolish the remaining genetic diversity in this population. Local population extinctions of collared lemmings would have severe effects on the arctic ecosystem, as collared lemmings are a key species in the trophic interactions and ecosystem processes in the Arctic.

A draft sequence of the Neandertal genome

Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

Computational challenges in the analysis of ancient DNA

A new method of next-generation sequencing analysis is presented which takes into account the biases characteristic of ancient, including Neandertal, DNA samples.

High-throughput sequencing technologies have opened up a new avenue for studying extinct organisms. Here we identify and quantify biases introduced by particular characteristics of ancient DNA samples. These analyses demonstrate the importance of closely related genomic sequence for correctly identifying and classifying bona fide endogenous DNA fragments. We show that more accurate genome divergence estimates from ancient DNA sequence can be attained using at least two outgroup genomes and appropriate filtering.

Phylogeographic history of grey wolves in Europe

Background

While it is generally accepted that patterns of intra-specific genetic differentiation are substantially affected by glacial history, population genetic processes occurring during Pleistocene glaciations are still poorly understood. In this study, we address the question of the genetic consequences of Pleistocene glaciations for European grey wolves. Combining our data with data from published studies, we analysed phylogenetic relationships and geographic distribution of mitochondrial DNA haplotypes for 947 contemporary European wolves. We also compared the contemporary wolf sequences with published sequences of 24 ancient European wolves.

Results

We found that haplotypes representing two haplogroups, 1 and 2, overlap geographically, but substantially differ in frequency between populations from south-western and eastern Europe. A comparison between haplotypes from Europe and other continents showed that both haplogroups are spread throughout Eurasia, while only haplogroup 1 occurs in contemporary North American wolves. All ancient wolf samples from western Europe that dated from between 44,000 and 1,200 years B.P. belonged to haplogroup 2, suggesting the long-term predominance of this haplogroup in this region. Moreover, a comparison of current and past frequencies and distributions of the two haplogroups in Europe suggested that haplogroup 2 became outnumbered by haplogroup 1 during the last several thousand years.

Conclusions

Parallel haplogroup replacement, with haplogroup 2 being totally replaced by haplogroup 1, has been reported for North American grey wolves. Taking into account the similarity of diets reported for the late Pleistocene wolves from Europe and North America, the correspondence between these haplogroup frequency changes may suggest that they were associated with ecological changes occurring after the Last Glacial Maximum.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius)

BACKGROUND

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

METHODOLOGY

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

CONCLUSIONS

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

Characterization of nucleotide misincorporation patterns in the iceman's mitochondrial DNA

Background

The degradation of DNA represents one of the main issues in the genetic analysis of archeological specimens. In the recent years, a particular kind of post-mortem DNA modification giving rise to nucleotide misincorporation (“miscoding lesions”) has been the object of extensive investigations.

Methodology/Principal Findings

To improve our knowledge regarding the nature and incidence of ancient DNA nucleotide misincorporations, we have utilized 6,859 (629,975 bp) mitochondrial (mt) DNA sequences obtained from the 5,350–5,100-years-old, freeze-desiccated human mummy popularly known as the Tyrolean Iceman or Ötzi. To generate the sequences, we have applied a mixed PCR/pyrosequencing procedure allowing one to obtain a particularly high sequence coverage. As a control, we have produced further 8,982 (805,155 bp) mtDNA sequences from a contemporary specimen using the same system and starting from the same template copy number of the ancient sample. From the analysis of the nucleotide misincorporation rate in ancient, modern, and putative contaminant sequences, we observed that the rate of misincorporation is significantly lower in modern and putative contaminant sequence datasets than in ancient sequences. In contrast, type 2 transitions represent the vast majority (85%) of the observed nucleotide misincorporations in ancient sequences.

Conclusions/Significance

This study provides a further contribution to the knowledge of nucleotide misincorporation patterns in DNA sequences obtained from freeze-preserved archeological specimens. In the Iceman system, ancient sequences can be clearly distinguished from contaminants on the basis of nucleotide misincorporation rates. This observation confirms a previous identification of the ancient mummy sequences made on a purely phylogenetical basis. The present investigation provides further indication that the majority of ancient DNA damage is reflected by type 2 (cytosine→thymine/guanine→adenine) transitions and that type 1 transitions are essentially PCR artifacts.

Removal of deaminated cytosines and detection of in vivo methylation in ancient DNA

DNA sequences determined from ancient organisms have high error rates, primarily due to uracil bases created by cytosine deamination. We use synthetic oligonucleotides, as well as DNA extracted from mammoth and Neandertal remains, to show that treatment with uracil-DNA-glycosylase and endonuclease VIII removes uracil residues from ancient DNA and repairs most of the resulting abasic sites, leaving undamaged parts of the DNA fragments intact. Neandertal DNA sequences determined with this protocol have greatly increased accuracy. In addition, our results demonstrate that Neandertal DNA retains in vivo patterns of CpG methylation, potentially allowing future studies of gene inactivation and imprinting in ancient organisms.

Preferential access to genetic information from endogenous hominin ancient DNA and accurate quantitative SNP-typing via SPEX

The analysis of targeted genetic loci from ancient, forensic and clinical samples is usually built upon polymerase chain reaction (PCR)-generated sequence data. However, many studies have shown that PCR amplification from poor-quality DNA templates can create sequence artefacts at significant levels. With hominin (human and other hominid) samples, the pervasive presence of highly PCR-amplifiable human DNA contaminants in the vast majority of samples can lead to the creation of recombinant hybrids and other non-authentic artefacts. The resulting PCR-generated sequences can then be difficult, if not impossible, to authenticate. In contrast, single primer extension (SPEX)-based approaches can genotype single nucleotide polymorphisms from ancient fragments of DNA as accurately as modern DNA. A single SPEX-type assay can amplify just one of the duplex DNA strands at target loci and generate a multi-fold depth-of-coverage, with non-authentic recombinant hybrids reduced to undetectable levels. Crucially, SPEX-type approaches can preferentially access genetic information from damaged and degraded endogenous ancient DNA templates over modern human DNA contaminants. The development of SPEX-type assays offers the potential for highly accurate, quantitative genotyping from ancient hominin samples.

The Neandertal genome and ancient DNA authenticity

Recent advances in high-thoughput DNA sequencing have made genome-scale analyses of genomes of extinct organisms possible. With these new opportunities come new difficulties in assessing the authenticity of the DNA sequences retrieved. We discuss how these difficulties can be addressed, particularly with regard to analyses of the Neandertal genome. We argue that only direct assays of DNA sequence positions in which Neandertals differ from all contemporary humans can serve as a reliable means to estimate human contamination. Indirect measures, such as the extent of DNA fragmentation, nucleotide misincorporations, or comparison of derived allele frequencies in different fragment size classes, are unreliable. Fortunately, interim approaches based on mtDNA differences between Neandertals and current humans, detection of male contamination through Y chromosomal sequences, and repeated sequencing from the same fossil to detect autosomal contamination allow initial large-scale sequencing of Neandertal genomes. This will result in the discovery of fixed differences in the nuclear genome between Neandertals and current humans that can serve as future direct assays for contamination. For analyses of other fossil hominins, which may become possible in the future, we suggest a similar ‘boot-strap' approach in which interim approaches are applied until sufficient data for more definitive direct assays are acquired.

Archaeogenetic evidence of ancient nubian barley evolution from six to two-row indicates local adaptation

Background

Archaeobotanical samples of barley (Hordeum vulgare L.) found at Qasr Ibrim display a two-row phenotype that is unique to the region of archaeological sites upriver of the first cataract of the Nile, characterised by the development of distinctive lateral bracts. The phenotype occurs throughout all strata at Qasr Ibrim, which range in age from 3000 to a few hundred years.

Methodology and Findings

We extracted ancient DNA from barley samples from the entire range of occupancy of the site, and studied the Vrs1 gene responsible for row number in extant barley. Surprisingly, we found a discord between the genotype and phenotype in all samples; all the barley had a genotype consistent with the six-row condition. These results indicate a six-row ancestry for the Qasr Ibrim barley, followed by a reassertion of the two-row condition. Modelling demonstrates that this sequence of evolutionary events requires a strong selection pressure.

Conclusions

The two-row phenotype at Qasr Ibrim is caused by a different mechanism to that in extant barley. The strength of selection required for this mechanism to prevail indicates that the barley became locally adapted in the region in response to a local selection pressure. The consistency of the genotype/phenotype discord over time supports a scenario of adoption of this barley type by successive cultures, rather than the importation of new barley varieties associated with individual cultures.

Paleogenomics in a temperate environment: shotgun sequencing from an extinct Mediterranean caprine

Background

Numerous endemic mammals, including dwarf elephants, goats, hippos and deers, evolved in isolation in the Mediterranean islands during the Pliocene and Pleistocene. Most of them subsequently became extinct during the Holocene. Recently developed high-throughput sequencing technologies could provide a unique tool for retrieving genomic data from these extinct species, making it possible to study their evolutionary history and the genetic bases underlying their particular, sometimes unique, adaptations.

Methodology/Principals Findings

A DNA extraction of a ∼6,000 year-old bone sample from an extinct caprine (Myotragus balearicus) from the Balearic Islands in the Western Mediterranean, has been subjected to shotgun sequencing with the GS FLX 454 platform. Only 0.27% of the resulting sequences, identified from alignments with the cow genome and comprising 15,832 nucleotides, with an average length of 60 nucleotides, proved to be endogenous.

Conclusions

A phylogenetic tree generated with Myotragus sequences and those from other artiodactyls displays an identical topology to that generated from mitochondrial DNA data. Despite being in an unfavourable thermal environment, which explains the low yield of endogenous sequences, our study demonstrates that it is possible to obtain genomic data from extinct species from temperate regions.

New insights from old bones: DNA preservation and degradation in permafrost preserved mammoth remains

Despite being plagued by heavily degraded DNA in palaeontological remains, most studies addressing the state of DNA degradation have been limited to types of damage which do not pose a hindrance to Taq polymerase during PCR. Application of serial qPCR to the two fractions obtained during extraction (demineralization and protein digest) from six permafrost mammoth bones and one partially degraded modern elephant bone has enabled further insight into the changes which endogenous DNA is subjected to during diagenesis. We show here that both fractions exhibit individual qualities in terms of the prevailing type of DNA (i.e. mitochondrial versus nuclear DNA) as well as the extent of damage, and in addition observed a highly variable ratio of mitochondrial to nuclear DNA among the six mammoth samples. While there is evidence suggesting that mitochondrial DNA is better preserved than nuclear DNA in ancient permafrost samples, we find the initial DNA concentration in the bone tissue to be as relevant for the total accessible mitochondrial DNA as the extent of DNA degradation post-mortem. We also evaluate the general applicability of indirect measures of preservation such as amino-acid racemization, bone crystallinity index and thermal age to these exceptionally well-preserved samples.

Sequences from first settlers reveal rapid evolution in Icelandic mtDNA pool

A major task in human genetics is to understand the nature of the evolutionary processes that have shaped the gene pools of contemporary populations. Ancient DNA studies have great potential to shed light on the evolution of populations because they provide the opportunity to sample from the same population at different points in time. Here, we show that a sample of mitochondrial DNA (mtDNA) control region sequences from 68 early medieval Icelandic skeletal remains is more closely related to sequences from contemporary inhabitants of Scotland, Ireland, and Scandinavia than to those from the modern Icelandic population. Due to a faster rate of genetic drift in the Icelandic mtDNA pool during the last 1,100 years, the sequences carried by the first settlers were better preserved in their ancestral gene pools than among their descendants in Iceland. These results demonstrate the inferential power gained in ancient DNA studies through the application of population genetics analyses to relatively large samples.

Ancient DNA studies have great potential to shed light on the evolution of populations because they provide the opportunity to sample from the same population at different points in time. However, ancient DNA studies are often based on DNA extracted from only one or a few individuals and, therefore, do not lend themselves to statistical inference. Here, we describe the analysis of a sample of mitochondrial DNA (mtDNA) control region sequences from 68 Icelandic skeletal remains that are about 1,000 years old, from the time that Iceland was first settled. We show that the ancient Icelandic mtDNA sequences are more closely related to sequences from contemporary inhabitants of Scotland, Ireland, and Scandinavia (and several other European populations) than to those from the modern Icelandic population. It appears that the array of sequences carried by the first generations of Icelanders was better preserved in the gene pools of their ancestors than among their modern descendants because of a faster rate of evolution due to genetic drift in the Icelandic mtDNA pool during the last 1,100 years. These results demonstrate the inferential power that can be gained from studies by applying the methods of population genetics to samples of ancient DNA sequences.

The development and impact of 454 sequencing

The 454 Sequencer has dramatically increased the volume of sequencing conducted by the scientific community and expanded the range of problems that can be addressed by the direct readouts of DNA sequence. Key breakthroughs in the development of the 454 sequencing platform included higher throughput, simplified all in vitro sample preparation and the miniaturization of sequencing chemistries, enabling massively parallel sequencing reactions to be carried out at a scale and cost not previously possible. Together with other recently released next-generation technologies, the 454 platform has started to democratize sequencing, providing individual laboratories with access to capacities that rival those previously found only at a handful of large sequencing centers. Over the past 18 months, 454 sequencing has led to a better understanding of the structure of the human genome, allowed the first non-Sanger sequence of an individual human and opened up new approaches to identify small RNAs. To make next-generation technologies more widely accessible, they must become easier to use and less costly. In the longer term, the principles established by 454 sequencing might reduce cost further, potentially enabling personalized genomics.

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.

Postmortem miscoding lesions in sequence analysis of human ancient mitochondrial DNA

Genetic miscoding lesions can cause inaccuracies during the interpretation of ancient DNA sequence data. In this study, genetic miscoding lesions were identified and assessed by cloning and direct sequencing of degraded, amplified mitochondrial DNA (mtDNA) extracted from human remains. Forty-two individuals, comprising nine collections from five geographic locations, were analyzed for the presence of DNA damage that can affect the generation of a correct mtDNA profile. In agreement with previous studies, high levels (56.5% of all damage sites) of proposed hydrolytic damage products were observed. Among these, type 2 transitions (cytosine --> thymine or guanine --> adenine), which are highly indicative of hydrolytic deamination, were observed in 50% of all misincorporations that occurred. In addition to hydrolytic damage products, oxidative damage products were also observed in this study and were responsible for approximately 43.5% of all misincorporations. This level of misincorporation is in contrast to previous studies characterizing miscoding lesions from the analysis of bone and teeth, where few to no oxidative damage products were observed. Of all the oxidative damage products found in this study, type 2 transversions (cytosine --> adenine/guanine --> thymine or cytosine --> guanine/guanine --> cytosine), which are commonly formed through the generation of 8-hydroxyguanine, accounted for 30.3% of all genetic miscoding lesions observed. This study identifies the previously unreported presence of oxidative DNA damage and proposes that damage to degraded DNA templates is highly specific in type, correlating with the geographic location and the taphonomic conditions of the depositional environment from which the remains are recovered.

A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.