Pulling out the 1%: Whole-Genome Capture for the Targeted Enrichment of Ancient DNA Sequencing Libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062–147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217–73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.     

Ancient DNA: would the real Neandertal please stand up?

Mitochondrial DNA sequences recovered from eight Neandertal specimens cannot be detected in either early fossil Europeans or in modern populations. This indicates that, if Neandertals made any genetic contribution at all to modern humans, it must have been limited, though the extent of the contribution cannot be resolved at present.     

Distribution patterns of postmortem damage in human mitochondrial DNA

The distribution of postmortem damage in mitochondrial DNA retrieved from 37 ancient human DNA samples was analyzed by cloning and was compared with a selection of published animal data. A relative rate of damage (rho(v)) was calculated for nucleotide positions within the human hypervariable region 1 (HVR1) and cytochrome oxidase subunit III genes. A comparison of damaged sites within and between the regions reveals that damage hotspots exist and that, in the HVR1, these correlate with sites known to have high in vivo mutation rates. Conversely, HVR1 subregions with known structural function, such as MT5, have lower in vivo mutation rates and lower postmortem-damage rates. The postmortem data also identify a possible functional subregion of the HVR1, termed "low-diversity 1," through the lack of sequence damage. The amount of postmortem damage observed in mitochondrial coding regions was significantly lower than in the HVR1, and, although hotspots were noted, these did not correlate with codon position. Finally, a simple method for the identification of incorrect archaeological haplogroup designations is introduced, on the basis of the observed spectrum of postmortem damage.     

The genetic origins of the Andaman Islanders

Mitochondrial sequences were retrieved from museum specimens of the enigmatic Andaman Islanders to analyze their evolutionary history. D-loop and protein-coding data reveal that phenotypic similarities with African pygmoid groups are convergent. Genetic and epigenetic data are interpreted as favoring the long-term isolation of the Andamanese, extensive population substructure, and/or two temporally distinct settlements. An early colonization featured populations bearing mtDNA lineage M2, and this lineage is hypothesized to represent the phylogenetic signal of an early southern movement of humans through Asia. The results demonstrate that Victorian anthropological collections can be used to study extinct, or seriously admixed populations, to provide new data about early human origins.     

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

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

Molecular diet analysis of two african free-tailed bats (molossidae) using high throughput sequencing

Given the diversity of prey consumed by insectivorous bats, it is difficult to discern the composition of their diet using morphological or conventional PCR-based analyses of their faeces. We demonstrate the use of a powerful alternate tool, the use of the Roche FLX sequencing platform to deep-sequence uniquely 5' tagged insect-generic barcode cytochrome c oxidase I (COI) fragments, that were PCR amplified from faecal pellets of two free-tailed bat species Chaerephon pumilus and Mops condylurus (family: Molossidae). Although the analyses were challenged by the paucity of southern African insect COI sequences in the GenBank and BOLD databases, similarity to existing collections allowed the preliminary identification of 25 prey families from six orders of insects within the diet of C. pumilus, and 24 families from seven orders within the diet of M. condylurus. Insects identified to families within the orders Lepidoptera and Diptera were widely present among the faecal samples analysed. The two families that were observed most frequently were Noctuidae and Nymphalidae (Lepidoptera). Species-level analysis of the data was accomplished using novel bioinformatics techniques for the identification of molecular operational taxonomic units (MOTU). Based on these analyses, our data provide little evidence of resource partitioning between sympatric M. condylurus and C. pumilus in the Simunye region of Swaziland at the time of year when the samples were collected, although as more complete databases against which to compare the sequences are generated this may have to be re-evaluated.     

Juvenile songbirds compensate for displacement to oceanic islands during autumn migration

To what degree juvenile migrant birds are able to correct for orientation errors or wind drift is still largely unknown. We studied the orientation of passerines on the Faroe Islands far off the normal migration routes of European migrants. The ability to compensate for displacement was tested in naturally occurring vagrants presumably displaced by wind and in birds experimentally displaced 1100 km from Denmark to the Faroes. The orientation was studied in orientation cages as well as in the free-flying birds after release by tracking departures using small radio transmitters. Both the naturally displaced and the experimentally displaced birds oriented in more easterly directions on the Faroes than was observed in Denmark prior to displacement. This pattern was even more pronounced in departure directions, perhaps because of wind influence. The clear directional compensation found even in experimentally displaced birds indicates that first-year birds can also possess the ability to correct for displacement in some circumstances, possibly involving either some primitive form of true navigation, or 'sign posts', but the cues used for this are highly speculative. We also found some indications of differences between species in the reaction to displacement. Such differences might be involved in the diversity of results reported in displacement studies so far.     

Beringian paleoecology inferred from permafrost-preserved fungal DNA

The diversity of fungi in permanently frozen soil from northeastern Siberia was studied by culture-independent PCR amplification of diverse environmental 18S rRNA genes. Elaborate protocols to avoid contamination during drilling, sampling, and amplification were used. A broad diversity of eukaryotic DNA sequences that were 510 bp long, including sequences of various fungi, plants, and invertebrates, could be obtained reproducibly from samples that were up to 300,000 to 400,000 years old. The sequences revealed that ancient fungal communities included a diversity of cold-adapted yeasts, dark-pigmented fungi, plant-parasitic fungi, and lichen mycobionts. DNA traces of tree-associated macrofungi in a modern tundra sample indicated that there was a shift in fungal diversity following the last ice age and supported recent results showing that there was a severe change in the plant composition in northeastern Siberia during this period. Interestingly, DNA sequences with high homology to sequences of coprophilic and keratinophilic fungi indicated that feces, hair, skin, and nails could have been sources of ancient megafauna DNA recently reported to be present in small amounts of Siberian permafrost sediments.     

Geologically ancient DNA: fact or artefact?

Studies continue to report ancient DNA sequences and viable microbial cells that are many millions of years old. In this paper we evaluate some of the most extravagant claims of geologically ancient DNA. We conclude that although exciting, the reports suffer from inadequate experimental setup and insufficient authentication of results. Consequently, it remains doubtful whether amplifiable DNA sequences and viable bacteria can survive over geological timescales. To enhance the credibility of future studies and assist in discarding false-positive results, we propose a rigorous set of authentication criteria for work with geologically ancient DNA.