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

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.

Genome-wide comparison of medieval and modern Mycobacterium leprae

Leprosy was endemic in Europe until the Middle Ages. Using DNA array capture, we have obtained genome sequences of Mycobacterium leprae from skeletons of five medieval leprosy cases from the United Kingdom, Sweden, and Denmark. In one case, the DNA was so well preserved that full de novo assembly of the ancient bacterial genome could be achieved through shotgun sequencing alone. The ancient M. leprae sequences were compared with those of 11 modern strains, representing diverse genotypes and geographic origins. The comparisons revealed remarkable genomic conservation during the past 1000 years, a European origin for leprosy in the Americas, and the presence of an M. leprae genotype in medieval Europe now commonly associated with the Middle East. The exceptional preservation of M. leprae biomarkers, both DNA and mycolic acids, in ancient skeletons has major implications for palaeomicrobiology and human pathogen evolution.

The rise and fall of the Phytophthora infestans lineage that triggered the Irish potato famine

Phytophthora infestans, the cause of potato late blight, is infamous for having triggered the Irish Great Famine in the 1840s. Until the late 1970s, P. infestans diversity outside of its Mexican center of origin was low, and one scenario held that a single strain, US-1, had dominated the global population for 150 years; this was later challenged based on DNA analysis of historical herbarium specimens. We have compared the genomes of 11 herbarium and 15 modern strains. We conclude that the 19th century epidemic was caused by a unique genotype, HERB-1, that persisted for over 50 years. HERB-1 is distinct from all examined modern strains, but it is a close relative of US-1, which replaced it outside of Mexico in the 20th century. We propose that HERB-1 and US-1 emerged from a metapopulation that was established in the early 1800s outside of the species' center of diversity. DOI:http://dx.doi.org/10.7554/eLife.00731.001.

The rise and fall of the Phytophthora infestans lineage that triggered the Irish potato famine

Phytophthora infestans, the cause of potato late blight, is infamous for having triggered the Irish Great Famine in the 1840s. Until the late 1970s, P. infestans diversity outside of its Mexican center of origin was low, and one scenario held that a single strain, US-1, had dominated the global population for 150 years; this was later challenged based on DNA analysis of historical herbarium specimens. We have compared the genomes of 11 herbarium and 15 modern strains. We conclude that the 19th century epidemic was caused by a unique genotype, HERB-1, that persisted for over 50 years. HERB-1 is distinct from all examined modern strains, but it is a close relative of US-1, which replaced it outside of Mexico in the 20th century. We propose that HERB-1 and US-1 emerged from a metapopulation that was established in the early 1800s outside of the species' center of diversity. DOI:http://dx.doi.org/10.7554/eLife.00731.001.

Single-stranded DNA library preparation for the sequencing of ancient or damaged DNA

This protocol describes a method for converting short single-stranded and double-stranded DNA into libraries compatible with high-throughput sequencing using Illumina technology. This method has primarily been developed to improve sequence retrieval from ancient DNA, but it is also applicable to the sequencing of short or degraded DNA from other sources, and it can also be used for sequencing oligonucleotides. Single-stranded library preparation is performed by ligating a biotinylated adapter oligonucleotide to the 3' ends of heat-denatured DNA. The resulting strands are then immobilized on streptavidin-coated beads and copied with a polymerase. A second adapter is attached by blunt-end ligation, and library preparation is completed by PCR amplification. We estimate that intact DNA strands are recovered in the library with ∼50% efficiency. Libraries can be generated from up to 12 DNA or oligonucleotide samples in parallel within 2 d.

DNA analysis of an early modern human from Tianyuan Cave, China

Hominins with morphology similar to present-day humans appear in the fossil record across Eurasia between 40,000 and 50,000 y ago. The genetic relationships between these early modern humans and present-day human populations have not been established. We have extracted DNA from a 40,000-y-old anatomically modern human from Tianyuan Cave outside Beijing, China. Using a highly scalable hybridization enrichment strategy, we determined the DNA sequences of the mitochondrial genome, the entire nonrepetitive portion of chromosome 21 (∼30 Mbp), and over 3,000 polymorphic sites across the nuclear genome of this individual. The nuclear DNA sequences determined from this early modern human reveal that the Tianyuan individual derived from a population that was ancestral to many present-day Asians and Native Americans but postdated the divergence of Asians from Europeans. They also show that this individual carried proportions of DNA variants derived from archaic humans similar to present-day people in mainland Asia.

Homo sapiens, Homo neanderthalensis and the Denisova specimen: New insights on their evolutionary histories using whole-genome comparisons

After a brief review of the most recent findings in the study of human evolution, an extensive comparison of the complete genomes of our nearest relative, the chimpanzee (Pan troglodytes), of extant Homo sapiens, archaic Homo neanderthalensis and the Denisova specimen were made. The focus was on non-synonymous mutations, which consequently had an impact on protein levels and these changes were classified according to degree of effect. A total of 10,447 non-synonymous substitutions were found in which the derived allele is fixed or nearly fixed in humans as compared to chimpanzee. Their most frequent location was on chromosome 21. Their presence was then searched in the two archaic genomes. Mutations in 381 genes would imply radical amino acid changes, with a fraction of these related to olfaction and other important physiological processes. Eight new alleles were identified in the Neanderthal and/or Denisova genetic pools. Four others, possibly affecting cognition, occured both in the sapiens and two other archaic genomes. The selective sweep that gave rise to Homo sapiens could, therefore, have initiated before the modern/archaic human divergence.

A recent evolutionary change affects a regulatory element in the human FOXP2 gene

The FOXP2 gene is required for normal development of speech and language. By isolating and sequencing FOXP2 genomic DNA fragments from a 49,000-year-old Iberian Neandertal and 50 present-day humans, we have identified substitutions in the gene shared by all or nearly all present-day humans but absent or polymorphic in Neandertals. One such substitution is localized in intron 8 and affects a binding site for the transcription factor POU3F2, which is highly conserved among vertebrates. We find that the derived allele of this site is less efficient than the ancestral allele in activating transcription from a reporter construct. The derived allele also binds less POU3F2 dimers than POU3F2 monomers compared with the ancestral allele. Because the substitution in the POU3F2 binding site is likely to alter the regulation of FOXP2 expression, and because it is localized in a region of the gene associated with a previously described signal of positive selection, it is a plausible candidate for having caused a recent selective sweep in the FOXP2 gene.

Archaic human genomics

For much of the 20th century, the predominant view of human evolutionary history was derived from the fossil record. Homo erectus was seen arising in Africa from an earlier member of the genus and then spreading throughout the Old World and into the Oceania. A regional continuity model of anagenetic change from H. erectus via various intermediate archaic species into the modern humans in each of the regions inhabited by H. erectus was labeled the multiregional model of human evolution (MRE). A contrasting model positing a single origin, in Africa, of anatomically modern H. sapiens with some populations later migrating out of Africa and replacing the local archaic populations throughout the world with complete replacement became known as the recent African origin (RAO) model. Proponents of both models used different interpretations of the fossil record to bolster their views for decades. In the 1980s, molecular genetic techniques began providing evidence from modern human variation that allowed not only the different models of modern human origins to be tested but also the exploration demographic history and the types of selection that different regions of the genome and even specific traits had undergone. The majority of researchers interpreted these data as strongly supporting the RAO model, especially analyses of mitochondrial DNA (mtDNA). Extrapolating backward from modern patterns of variation and using various calibration points and substitution rates, a consensus arose that saw modern humans evolving from an African population around 200,000 years ago. Much later, around 50,000 years ago, a subset of this population migrated out of Africa replacing Neanderthals in Europe and western Asia as well as archaics in eastern Asia and Oceania. mtDNA sequences from more than two-dozen Neanderthals and early modern humans re-enforced this consensus. In 2010, however, the complete draft genomes of Neanderthals and of heretofore unknown hominins from Siberia, called Denisovans, demonstrated gene flow between these archaic human species and modern Eurasians but not sub-Saharan Africans. Although the levels of gene flow may be very limited, this unexpected finding does not fit well with either the RAO model or MRE model. More thorough sampling of modern human diversity, additional fossil discoveries, and the sequencing of additional hominin fossils are necessary to throw light onto our origins and our history.

Length and GC-biases during sequencing library amplification: a comparison of various polymerase-buffer systems with ancient and modern DNA sequencing libraries

High-throughput sequencing technologies frequently necessitate the use of PCR for sequencing library amplification. PCR is a sometimes enigmatic process and is known to introduce biases. Here we perform a simple amplification-sequencing assay using 10 commercially available polymerase-buffer systems to amplify libraries prepared from both modern and ancient DNA. We compare the performance of the polymerases with respect to a previously uncharacterized template length bias, as well as GC-content bias, and find that simply avoiding certain polymerase can dramatically decrease the occurrence of both. For amplification of ancient DNA, we found that some commonly used polymerases strongly bias against amplification of endogenous DNA in favor of GC-rich microbial contamination, in our case reducing the fraction of endogenous sequences to almost half.

A high-coverage genome sequence from an archaic Denisovan individual

We present a DNA library preparation method that has allowed us to reconstruct a high-coverage (30×) genome sequence of a Denisovan, an extinct relative of Neandertals. The quality of this genome allows a direct estimation of Denisovan heterozygosity indicating that genetic diversity in these archaic hominins was extremely low. It also allows tentative dating of the specimen on the basis of "missing evolution" in its genome, detailed measurements of Denisovan and Neandertal admixture into present-day human populations, and the generation of a near-complete catalog of genetic changes that swept to high frequency in modern humans since their divergence from Denisovans.

[Progresses on Neandertal genomics]

Neandertal is our closest known relative and also an archaic hominid reserving the richest fossils. Whether the Neandertals exchanged their DNA with modern human or not is a matter of debate on the modern human origin. The progresses on the mitochondrial and nuclear genomes of Neandertals in recent years were reviewed in this paper. Recent study has revealed possible genetic contribution of Neandertals to the modern human to some extent, which arose the rethinking of modern human origin. The experiences gained in the research on Neandertals will benefit the study on archaic hominids, unravel the mystery of modern human origin, and enrich the relative theoretical systems in evolutionary biological field.

Improving the performance of true single molecule sequencing for ancient DNA

BACKGROUND

Second-generation sequencing technologies have revolutionized our ability to recover genetic information from the past, allowing the characterization of the first complete genomes from past individuals and extinct species. Recently, third generation Helicos sequencing platforms, which perform true Single-Molecule DNA Sequencing (tSMS), have shown great potential for sequencing DNA molecules from Pleistocene fossils. Here, we aim at improving even further the performance of tSMS for ancient DNA by testing two novel tSMS template preparation methods for Pleistocene bone fossils, namely oligonucleotide spiking and treatment with DNA phosphatase.

RESULTS

We found that a significantly larger fraction of the horse genome could be covered following oligonucleotide spiking however not reproducibly and at the cost of extra post-sequencing filtering procedures and skewed %GC content. In contrast, we showed that treating ancient DNA extracts with DNA phosphatase improved the amount of endogenous sequence information recovered per sequencing channel by up to 3.3-fold, while still providing molecular signatures of endogenous ancient DNA damage, including cytosine deamination and fragmentation by depurination. Additionally, we confirmed the existence of molecular preservation niches in large bone crystals from which DNA could be preferentially extracted.

CONCLUSIONS

We propose DNA phosphatase treatment as a mechanism to increase sequence coverage of ancient genomes when using Helicos tSMS as a sequencing platform. Together with mild denaturation temperatures that favor access to endogenous ancient templates over modern DNA contaminants, this simple preparation procedure can improve overall Helicos tSMS performance when damaged DNA templates are targeted.

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.

Next-generation sequencing offers new insights into DNA degradation

The processes underlying DNA degradation are central to various disciplines, including cancer research, forensics and archaeology. The sequencing of ancient DNA molecules on next-generation sequencing platforms provides direct measurements of cytosine deamination, depurination and fragmentation rates that previously were obtained only from extrapolations of results from in vitro kinetic experiments performed over short timescales. For example, recent next-generation sequencing of ancient DNA reveals purine bases as one of the main targets of postmortem hydrolytic damage, through base elimination and strand breakage. It also shows substantially increased rates of DNA base-loss at guanosine. In this review, we argue that the latter results from an electron resonance structure unique to guanosine rather than adenosine having an extra resonance structure over guanosine as previously suggested.

Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA

DNA that survives in museum specimens, bones and other tissues recovered by archaeologists is invariably fragmented and chemically modified. The extent to which such modifications accumulate over time is largely unknown but could potentially be used to differentiate between endogenous old DNA and present-day DNA contaminating specimens and experiments. Here we examine mitochondrial DNA sequences from tissue remains that vary in age between 18 and 60,000 years with respect to three molecular features: fragment length, base composition at strand breaks, and apparent C to T substitutions. We find that fragment length does not decrease consistently over time and that strand breaks occur preferentially before purine residues by what may be at least two different molecular mechanisms that are not yet understood. In contrast, the frequency of apparent C to T substitutions towards the 5′-ends of molecules tends to increase over time. These nucleotide misincorporations are thus a useful tool to distinguish recent from ancient DNA sources in specimens that have not been subjected to unusual or harsh treatments.

Analysis of human accelerated DNA regions using archaic hominin genomes

Several previous comparisons of the human genome with other primate and vertebrate genomes identified genomic regions that are highly conserved in vertebrate evolution but fast-evolving on the human lineage. These human accelerated regions (HARs) may be regions of past adaptive evolution in humans. Alternatively, they may be the result of non-adaptive processes, such as biased gene conversion. We captured and sequenced DNA from a collection of previously published HARs using DNA from an Iberian Neandertal. Combining these new data with shotgun sequence from the Neandertal and Denisova draft genomes, we determine at least one archaic hominin allele for 84% of all positions within HARs. We find that 8% of HAR substitutions are not observed in the archaic hominins and are thus recent in the sense that the derived allele had not come to fixation in the common ancestor of modern humans and archaic hominins. Further, we find that recent substitutions in HARs tend to have come to fixation faster than substitutions elsewhere in the genome and that substitutions in HARs tend to cluster in time, consistent with an episodic rather than a clock-like process underlying HAR evolution. Our catalog of sequence changes in HARs will help prioritize them for functional studies of genomic elements potentially responsible for modern human adaptations.

Rapid retrieval of DNA target sequences by primer extension capture

There is a widespread need for methods to enrich DNA samples for sequences of interest prior to high-throughput sequencing and to reduce the costs associated with a shotgun approach. While useful for targeting megabase-sized regions in a few samples, hybridization capture approaches such as those using microarrays currently involve bulky handling steps, long incubation times, and high per-sample costs. In contrast, the primer extension capture (PEC) method allows direct selection of small genomic regions from DNA sources within 2 h, with low costs for use with parallel samples. PEC promises useful applications in studies such as ancient DNA or forensic sequencing, taxonomic surveying of metagenomic samples, or genomic mapping of repetitive elements.

Ancient human DNA

The contribution of palaeogenetic data to the study of various aspects of hominin biology and evolution has been significant, and has the potential to increase substantially with the widespread implementation of next generation sequencing techniques. Here we discuss the present state-of-the-art of ancient human DNA analysis and the characteristics of hominin aDNA that make sequence validation particularly complex. A brief overview of the development of anthropological palaeogenetic analysis is given to illustrate the technical challenges motivating recent technological advancements.

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.

Analysis of high-throughput ancient DNA sequencing data

Advances in sequencing technologies have dramatically changed the field of ancient DNA (aDNA). It is now possible to generate an enormous quantity of aDNA sequence data both rapidly and inexpensively. As aDNA sequences are generally short in length, damaged, and at low copy number relative to coextracted environmental DNA, high-throughput approaches offer a tremendous advantage over traditional sequencing approaches in that they enable a complete characterization of an aDNA extract. However, the particular qualities of aDNA also present specific limitations that require careful consideration in data analysis. For example, results of high-throughout analyses of aDNA libraries may include chimeric sequences, sequencing error and artifacts, damage, and alignment ambiguities due to the short read lengths. Here, I describe typical primary data analysis workflows for high-throughput aDNA sequencing experiments, including (1) separation of individual samples in multiplex experiments; (2) removal of protocol-specific library artifacts; (3) trimming adapter sequences and merging paired-end sequencing data; (4) base quality score filtering or quality score propagation during data analysis; (5) identification of endogenous molecules from an environmental background; (6) quantification of contamination from other DNA sources; and (7) removal of clonal amplification products or the compilation of a consensus from clonal amplification products, and their exploitation for estimation of library complexity.

Preparation of next-generation sequencing libraries from damaged DNA

Next-generation sequencing (NGS) has revolutionized ancient DNA research, especially when combined with high-throughput target enrichment methods. However, attaining high sequencing depth and accuracy from samples often remains problematic due to the damaged state of ancient DNA, in particular the extremely low copy number of ancient DNA and the abundance of uracil residues derived from cytosine deamination that lead to miscoding errors. It is therefore critical to use a highly efficient procedure for conversion of a raw DNA extract into an adaptor-ligated sequencing library, and equally important to reduce errors from uracil residues. We present a protocol for NGS library preparation that allows highly efficient conversion of DNA fragments into an adaptor-ligated form. The protocol incorporates an option to remove the vast majority of uracil miscoding lesions as part of the library preparation process. The procedure requires only two spin column purification steps and no gel purification or bead handling. Starting from an aliquot of DNA extract, a finished, highly amplified library can be generated in 5 h, or under 3 h if uracil removal is not required.

Capillary electrophoresis of human mtDNA control region sequences from highly degraded samples using short mtDNA amplicons

The forensic applications of mtDNA sequencing center primarily on samples that are either highly degraded or contain little or no nuclear DNA, since the testing of these sample types is often unsuccessful with more widely used nuclear STR profiling assays. In these cases, sequence data from the noncoding mtDNA control region are targeted due to its high variability. However, the ease of authentic DNA recovery and the strategy used for recovery depend strictly on the quality of the sample. In this chapter, we will cover mitochondrial DNA sequencing procedures for short mtDNA amplicons which range in size from 100 to 350 bp. Generally speaking, amplicons of this size are required only for the most degraded specimens, and the protocols described here have been specifically developed for recalcitrant human skeletal remains encountered during the course of a large-scale missing persons' identification effort. DNA templates from these types of specimens tend to exhibit various forms of intrastrand damage that, in turn, manifest as artifacts in the sequence data. Because these artifacts are not generally observed among sequence data from pristine templates, we address the particular data idiosyncrasies that warrant additional scrutiny. Although this chapter will primarily highlight this particular application, the basic experimental parameters and data considerations should easily extend to other applications and/or sample types. The protocols described here have been deliberately designed to produce raw sequence electropherograms and final mtDNA profiles that adhere to the strictest forensic guidelines in terms of overall data quality.

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.

Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform

Due to the increasing throughput of current DNA sequencing instruments, sample multiplexing is necessary for making economical use of available sequencing capacities. A widely used multiplexing strategy for the Illumina Genome Analyzer utilizes sample-specific indexes, which are embedded in one of the library adapters. However, this and similar multiplex approaches come with a risk of sample misidentification. By introducing indexes into both library adapters (double indexing), we have developed a method that reveals the rate of sample misidentification within current multiplex sequencing experiments. With ~0.3% these rates are orders of magnitude higher than expected and may severely confound applications in cancer genomics and other fields requiring accurate detection of rare variants. We identified the occurrence of mixed clusters on the flow as the predominant source of error. The accuracy of sample identification is further impaired if indexed oligonucleotides are cross-contaminated or if indexed libraries are amplified in bulk. Double-indexing eliminates these problems and increases both the scope and accuracy of multiplex sequencing on the Illumina platform.

Learning about human population history from ancient and modern genomes

Genome-wide data, both from SNP arrays and from complete genome sequencing, are becoming increasingly abundant and are now even available from extinct hominins. These data are providing new insights into population history; in particular, when combined with model-based analytical approaches, genome-wide data allow direct testing of hypotheses about population history. For example, genome-wide data from both contemporary populations and extinct hominins strongly support a single dispersal of modern humans from Africa, followed by two archaic admixture events: one with Neanderthals somewhere outside Africa and a second with Denisovans that (so far) has only been detected in New Guinea. These new developments promise to reveal new stories about human population history, without having to resort to storytelling.

Addressing challenges in the production and analysis of illumina sequencing data

Advances in DNA sequencing technologies have made it possible to generate large amounts of sequence data very rapidly and at substantially lower cost than capillary sequencing. These new technologies have specific characteristics and limitations that require either consideration during project design, or which must be addressed during data analysis. Specialist skills, both at the laboratory and the computational stages of project design and analysis, are crucial to the generation of high quality data from these new platforms. The Illumina sequencers (including the Genome Analyzers I/II/IIe/IIx and the new HiScan and HiSeq) represent a widely used platform providing parallel readout of several hundred million immobilized sequences using fluorescent-dye reversible-terminator chemistry. Sequencing library quality, sample handling, instrument settings and sequencing chemistry have a strong impact on sequencing run quality. The presence of adapter chimeras and adapter sequences at the end of short-insert molecules, as well as increased error rates and short read lengths complicate many computational analyses. We discuss here some of the factors that influence the frequency and severity of these problems and provide solutions for circumventing these. Further, we present a set of general principles for good analysis practice that enable problems with sequencing runs to be identified and dealt with.

New insights in the removal of the hydantoins, oxidation product of pyrimidines, via the base excision and nucleotide incision repair pathways

Background

Oxidative damage to DNA, if not repaired, can be both miscoding and blocking. These genetic alterations can lead to mutations and/or cell death, which in turn cause cancer and aging. Oxidized DNA bases are substrates for two overlapping repair pathways: base excision (BER) and nucleotide incision repair (NIR). Hydantoin derivatives such as 5-hydroxyhydantoin (5OH-Hyd) and 5-methyl-5-hydroxyhydantoin (5OH-5Me-Hyd), major products of cytosine and thymine oxidative degradation pathways, respectively, have been detected in cancer cells and ancient DNA. Hydantoins are blocking lesions for DNA polymerases and excised by bacterial and yeast DNA glycosylases in the BER pathway. However little is known about repair of pyrimidine-derived hydantoins in human cells.

Methodology/Principal Findings

Here, using both denaturing PAGE and MALDI-TOF MS analyses we report that the bacterial, yeast and human AP endonucleases can incise duplex DNA 5′ next to 5OH-Hyd and 5OH-5Me-Hyd thus initiating the NIR pathway. We have fully reconstituted the NIR pathway for these lesions in vitro using purified human proteins. Depletion of Nfo in E. coli and APE1 in HeLa cells abolishes the NIR activity in cell-free extracts. Importantly, a number of redundant DNA glycosylase activities can excise hydantoin residues, including human NTH1, NEIL1 and NEIL2 and the former protein being a major DNA glycosylase activity in HeLa cells extracts.

Conclusions/Significance

This study demonstrates that both BER and NIR pathways can compete and/or back-up each other to remove hydantoin DNA lesions in vivo.

Genetic history of an archaic hominin group from Denisova Cave in Siberia

Using DNA extracted from a finger bone found in Denisova Cave in southern Siberia, we have sequenced the genome of an archaic hominin to about 1.9-fold coverage. This individual is from a group that shares a common origin with Neanderthals. This population was not involved in the putative gene flow from Neanderthals into Eurasians; however, the data suggest that it contributed 4-6% of its genetic material to the genomes of present-day Melanesians. We designate this hominin population 'Denisovans' and suggest that it may have been widespread in Asia during the Late Pleistocene epoch. A tooth found in Denisova Cave carries a mitochondrial genome highly similar to that of the finger bone. This tooth shares no derived morphological features with Neanderthals or modern humans, further indicating that Denisovans have an evolutionary history distinct from Neanderthals and modern humans.

Genomic DNA sequences from mastodon and woolly mammoth reveal deep speciation of forest and savanna elephants

To elucidate the history of living and extinct elephantids, we generated 39,763 bp of aligned nuclear DNA sequence across 375 loci for African savanna elephant, African forest elephant, Asian elephant, the extinct American mastodon, and the woolly mammoth. Our data establish that the Asian elephant is the closest living relative of the extinct mammoth in the nuclear genome, extending previous findings from mitochondrial DNA analyses. We also find that savanna and forest elephants, which some have argued are the same species, are as or more divergent in the nuclear genome as mammoths and Asian elephants, which are considered to be distinct genera, thus resolving a long-standing debate about the appropriate taxonomic classification of the African elephants. Finally, we document a much larger effective population size in forest elephants compared with the other elephantid taxa, likely reflecting species differences in ancient geographic structure and range and differences in life history traits such as variance in male reproductive success.

An efficient multistrategy DNA decontamination procedure of PCR reagents for hypersensitive PCR applications

Background

PCR amplification of minute quantities of degraded DNA for ancient DNA research, forensic analyses, wildlife studies and ultrasensitive diagnostics is often hampered by contamination problems. The extent of these problems is inversely related to DNA concentration and target fragment size and concern (i) sample contamination, (ii) laboratory surface contamination, (iii) carry-over contamination, and (iv) contamination of reagents.

Methodology/Principal Findings

Here we performed a quantitative evaluation of current decontamination methods for these last three sources of contamination, and developed a new procedure to eliminate contaminating DNA contained in PCR reagents. We observed that most current decontamination methods are either not efficient enough to degrade short contaminating DNA molecules, rendered inefficient by the reagents themselves, or interfere with the PCR when used at doses high enough to eliminate these molecules. We also show that efficient reagent decontamination can be achieved by using a combination of treatments adapted to different reagent categories. Our procedure involves γ- and UV-irradiation and treatment with a mutant recombinant heat-labile double-strand specific DNase from the Antarctic shrimp Pandalus borealis. Optimal performance of these treatments is achieved in narrow experimental conditions that have been precisely analyzed and defined herein.

Conclusions/Significance

There is not a single decontamination method valid for all possible contamination sources occurring in PCR reagents and in the molecular biology laboratory and most common decontamination methods are not efficient enough to decontaminate short DNA fragments of low concentration. We developed a versatile multistrategy decontamination procedure for PCR reagents. We demonstrate that this procedure allows efficient reagent decontamination while preserving the efficiency of PCR amplification of minute quantities of DNA.

Road blocks on paleogenomes--polymerase extension profiling reveals the frequency of blocking lesions in ancient DNA

Although the last few years have seen great progress in DNA sequence retrieval from fossil specimens, some of the characteristics of ancient DNA remain poorly understood. This is particularly true for blocking lesions, i.e. chemical alterations that cannot be bypassed by DNA polymerases and thus prevent amplification and subsequent sequencing of affected molecules. Some studies have concluded that the vast majority of ancient DNA molecules carry blocking lesions, suggesting that the removal, repair or bypass of blocking lesions might dramatically increase both the time depth and geographical range of specimens available for ancient DNA analysis. However, previous studies used very indirect detection methods that did not provide conclusive estimates on the frequency of blocking lesions in endogenous ancient DNA. We developed a new method, polymerase extension profiling (PEP), that directly reveals occurrences of polymerase stalling on DNA templates. By sequencing thousands of single primer extension products using PEP methodology, we have for the first time directly identified blocking lesions in ancient DNA on a single molecule level. Although we found clear evidence for blocking lesions in three out of four ancient samples, no more than 40% of the molecules were affected in any of the samples, indicating that such modifications are far less frequent in ancient DNA than previously thought.

Targeted investigation of the Neandertal genome by array-based sequence capture

It is now possible to perform whole-genome shotgun sequencing as well as capture of specific genomic regions for extinct organisms. However, targeted resequencing of large parts of nuclear genomes has yet to be demonstrated for ancient DNA. Here we show that hybridization capture on microarrays can successfully recover more than a megabase of target regions from Neandertal DNA even in the presence of approximately 99.8% microbial DNA. Using this approach, we have sequenced approximately 14,000 protein-coding positions inferred to have changed on the human lineage since the last common ancestor shared with chimpanzees. By generating the sequence of one Neandertal and 50 present-day humans at these positions, we have identified 88 amino acid substitutions that have become fixed in humans since our divergence from the Neandertals.

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.

The complete mitochondrial DNA genome of an unknown hominin from southern Siberia

With the exception of Neanderthals, from which DNA sequences of numerous individuals have now been determined, the number and genetic relationships of other hominin lineages are largely unknown. Here we report a complete mitochondrial (mt) DNA sequence retrieved from a bone excavated in 2008 in Denisova Cave in the Altai Mountains in southern Siberia. It represents a hitherto unknown type of hominin mtDNA that shares a common ancestor with anatomically modern human and Neanderthal mtDNAs about 1.0 million years ago. This indicates that it derives from a hominin migration out of Africa distinct from that of the ancestors of Neanderthals and of modern humans. The stratigraphy of the cave where the bone was found suggests that the Denisova hominin lived close in time and space with Neanderthals as well as with modern humans.