Three Reagents for in-Solution Enrichment of Ancient Human DNA at More than a Million SNPs

In-solution enrichment for hundreds of thousands of single nucleotide polymorphisms (SNPs) has been the source of >70% of all genome-scale ancient human DNA data published to date. This approach has made it possible to generate data for one to two orders of magnitude lower cost than random shotgun sequencing, making it economical to study ancient samples with low proportions of human DNA, and increasing the rate of conversion of sampled remains into working data thereby facilitating ethical stewardship of human remains. So far, nearly all ancient DNA data obtained using in-solution enrichment has been generated using a set of bait sequences targeting about 1.24 million SNPs (the 1240k reagent). These sequences were published in 2015, but synthesis of the reagent has been cost-effective for only a few laboratories. In 2021, two companies made available reagents that target the same core set of SNPs along with supplementary content. Here, we test the properties of the three reagents on a common set of 27 ancient DNA libraries across a range of richness of DNA content and percentages of human molecules. All three reagents are highly effective at enriching many hundreds of thousands of SNPs. For all three reagents and a wide range of conditions, one round of enrichment produces data that is as useful as two rounds when tens of millions of sequences are read out as is typical for such experiments. In our testing, the Twist Ancient DNA reagent produces the highest coverages, greatest uniformity on targeted positions, and almost no bias toward enriching one allele more than another relative to shotgun sequencing. Allelic bias in 1240k enrichment has made it challenging to carry out joint analysis of these data with shotgun data, creating a situation where the ancient DNA community has been publishing two important bodes of data that cannot easily be co-analyzed by population genetic methods. To address this challenge, we introduce a subset of hundreds of thousands of SNPs for which 1240k data can be effectively co-analyzed with all other major data types.

The genetic identity of the earliest human-made hybrid animals, the kungas of Syro-Mesopotamia

Ancient DNA identifies elite equids in Bronze Age Mesopotamia as the first human-bred hybrid animals, the enigmatic kungas.

Ancient DNA Methods Improve Forensic DNA Profiling of Korean War and World War II Unknowns

The integration of massively parallel sequencing (MPS) technology into forensic casework has been of particular benefit to the identification of unknown military service members. However, highly degraded or chemically treated skeletal remains often fail to provide usable DNA profiles, even with sensitive mitochondrial (mt) DNA capture and MPS methods. In parallel, the ancient DNA field has developed workflows specifically for degraded DNA, resulting in the successful recovery of nuclear DNA and mtDNA from skeletal remains as well as sediment over 100,000 years old. In this study we use a set of disinterred skeletal remains from the Korean War and World War II to test if ancient DNA extraction and library preparation methods improve forensic DNA profiling. We identified an ancient DNA extraction protocol that resulted in the recovery of significantly more human mtDNA fragments than protocols previously used in casework. In addition, utilizing single-stranded rather than double-stranded library preparation resulted in increased attainment of reportable mtDNA profiles. This study emphasizes that the combination of ancient DNA extraction and library preparation methods evaluated here increases the success rate of DNA profiling, and likelihood of identifying historical remains.

An empirical evaluation of genotype imputation of ancient DNA

With capabilities of sequencing ancient DNA to high coverage often limited by sample quality or cost, imputation of missing genotypes presents a possibility to increase power of inference as well as cost-effectiveness for the analysis of ancient data. However, the high degree of uncertainty often associated with ancient DNA poses several methodological challenges, and performance of imputation methods in this context has not been fully explored. To gain further insights, we performed a systematic evaluation of imputation of ancient data using Beagle 4.0 and reference data from phase 3 of the 1000 Genomes project, investigating the effects of coverage, phased reference and study sample size. Making use of five ancient samples with high-coverage data available, we evaluated imputed data with respect to accuracy, reference bias and genetic affinities as captured by PCA. We obtained genotype concordance levels of over 99% for data with 1x coverage, and similar levels of accuracy and reference bias at levels as low as 0.75x. Our findings suggest that using imputed data can be a realistic option for various population genetic analyses even for data in coverage ranges below 1x. We also show that a large and varied phased reference set as well as the inclusion of low- to moderate-coverage ancient samples can increase imputation performance, particularly for rare alleles. In-depth analysis of imputed data with respect to genetic variants and allele frequencies gave further insight into the nature of errors arising during imputation, and can provide practical guidelines for post-processing and validation prior to downstream analysis.

Ancient mitochondrial genomes recovered from small vertebrate bones through minimally destructive DNA extraction: phylogeography of the New Zealand gecko genus Hoplodactylus.

Methodological and technological improvements are continually revolutionizing the field of ancient DNA. Most ancient DNA extraction methods require the partial (or complete) destruction of finite museum specimens, which disproportionately impacts small or fragmentary subfossil remains, and future analyses. We present a minimally destructive ancient DNA extraction method optimized for small vertebrate remains. We applied these methods to detect lost mainland genetic diversity in the large New Zealand diplodactylid gecko genus Hoplodactylus, which is presently restricted to predator-free island sanctuaries. We present the first mitochondrial genomes for New Zealand diplodactylid geckos, recovered from 19 modern, six historic/archival (1898 to 2011) and 16 Holocene Hoplodactylus duvaucelii sensu latu specimens, and one modern Woodworthia sp. specimen. No obvious damage was observed in post-extraction micro-CT reconstructions. All ‘large gecko’ specimens examined from extinct populations were found to be conspecific with extant Hoplodactylus species, suggesting their large relative size evolved only once in the New Zealand diplodactylid radiation. Phylogenetic analyses of Hoplodactylus samples recovered two genetically (and morphologically) distinct North and South Island clades, probably corresponding to distinct species. Finer phylogeographic structuring within Hoplodactylus spp. highlighted the impacts of Late-Cenozoic biogeographic barriers, including the opening and closure of Pliocene marine straits, fluctuations in size and suitability of glacial refugia, and eustatic sea-level change. Recent mainland extinction obscured these signals from the modern tissue derived data. These results highlight the utility of minimally destructive DNA extraction in genomic analyses of less well studied small vertebrate taxa, and the conservation of natural history collections.