Improved detection of methylation in ancient DNA.

Reconstructing premortem DNA methylation levels in ancient DNA has led to breakthrough studies such as the prediction of anatomical features of the Denisovan. These studies rely on computationally inferring methylation levels from damage signals in naturally deaminated cytosines, which requires expensive high-coverage genomes. Here, we test two methods for direct methylation measurement developed for modern DNA based on either bisulfite or enzymatic methylation treatments. Bisulfite treatment shows the least reduction in DNA yields as well as the least biases during methylation conversion, demonstrating that this method can be successfully applied to ancient DNA.

Maximizing efficiency in sedimentary ancient DNA analysis: a novel extract pooling approach.

In the last few decades, the field of ancient DNA has taken a new direction towards using sedimentary ancient DNA (sedaDNA) for studying human and mammalian population dynamics as well as past ecosystems. However, the screening of numerous sediment samples from archaeological sites remains a time-consuming and costly endeavor, particularly when targeting hominin DNA. Here, we present a novel high-throughput method that facilitates the fast and efficient analysis of sediment samples by applying a pooled testing approach. This method combines multiple extracts, enabling early parallelization of laboratory procedures and effective aDNA screening. Pooled samples with detectable aDNA signals undergo detailed analysis, while empty pools are discarded. We have successfully applied our method to multiple sediment samples from Middle and Upper Paleolithic sites in Europe, Asia, and Africa. Notably, our results reveal that an aDNA signal remains discernible even when pooled with four negative samples. We also demonstrate that the DNA yield of double-stranded libraries increases significantly when reducing the extract input, potentially mitigating the effects of inhibition. By embracing this innovative approach, researchers can analyze large numbers of sediment samples for aDNA preservation, achieving significant cost reductions of up to 70% and reducing hands-on laboratory time to one-fifth.

Link between Monkeypox Virus Genomes from Museum Specimens and 1965 Zoo Outbreak.

We used pathogen genomics to test orangutan specimens from a museum in Bonn, Germany, to identify the origin of the animals and the circumstances of their death. We found monkeypox virus genomes in the samples and determined that they represent cases from a 1965 outbreak at Rotterdam Zoo in Rotterdam, the Netherlands.

Stable population structure in Europe since the Iron Age, despite high mobility.

Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000-3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

A genetic history of continuity and mobility in the Iron Age central Mediterranean.

The Iron Age was a dynamic period in central Mediterranean history, with the expansion of Greek and Phoenician colonies and the growth of Carthage into the dominant maritime power of the Mediterranean. These events were facilitated by the ease of long-distance travel following major advances in seafaring. We know from the archaeological record that trade goods and materials were moving across great distances in unprecedented quantities, but it is unclear how these patterns correlate with human mobility. Here, to investigate population mobility and interactions directly, we sequenced the genomes of 30 ancient individuals from coastal cities around the central Mediterranean, in Tunisia, Sardinia and central Italy. We observe a meaningful contribution of autochthonous populations, as well as highly heterogeneous ancestry including many individuals with non-local ancestries from other parts of the Mediterranean region. These results highlight both the role of local populations and the extreme interconnectedness of populations in the Iron Age Mediterranean. By studying these trans-Mediterranean neighbours together, we explore the complex interplay between local continuity and mobility that shaped the Iron Age societies of the central Mediterranean.

Density separation of petrous bone powders for optimized ancient DNA yields.

Density separation is a process routinely used to segregate minerals, organic matter, and even microplastics, from soils and sediments. Here we apply density separation to archaeological bone powders before DNA extraction to increase endogenous DNA recovery relative to a standard control extraction of the same powders. Using nontoxic heavy liquid solutions, we separated powders from the petrous bones of 10 individuals of similar archaeological preservation into eight density intervals (2.15 to 2.45 g/cm3, in 0.05 increments). We found that the 2.30 to 2.35 g/cm3 and 2.35 to 2.40 g/cm3 intervals yielded up to 5.28-fold more endogenous unique DNA than the corresponding standard extraction (and up to 8.53-fold before duplicate read removal), while maintaining signals of ancient DNA authenticity and not reducing library complexity. Although small 0.05 g/cm3 intervals may maximally optimize yields, a single separation to remove materials with a density above 2.40 g/cm3 yielded up to 2.57-fold more endogenous DNA on average, which enables the simultaneous separation of samples that vary in preservation or in the type of material analyzed. While requiring no new ancient DNA laboratory equipment and fewer than 30 min of extra laboratory work, the implementation of density separation before DNA extraction can substantially boost endogenous DNA yields without decreasing library complexity. Although subsequent studies are required, we present theoretical and practical foundations that may prove useful when applied to other ancient DNA substrates such as teeth, other bones, and sediments.

Northeastern Asian and Jomon-related genetic structure in the Three Kingdoms period of Gimhae, Korea.

The genetic history of prehistoric and protohistoric Korean populations is not well understood because only a small number of ancient genomes are available. Here, we report the first paleogenomic data from the Korean Three Kingdoms period, a crucial point in the cultural and historic formation of Korea. These data comprise eight shotgun-sequenced genomes from ancient Korea (0.7×-6.1× coverage). They were derived from two archeological sites in Gimhae: the Yuha-ri shell mound and the Daesung-dong tumuli, the latter being the most important funerary complex of the Gaya confederacy. All individuals are from between the 4th and 5th century CE and are best modeled as an admixture between a northern China Bronze Age genetic source and a source of Jomon-related ancestry that shares similarities with the present-day genomes from Japan. The observed substructure and proportion of Jomon-related ancestry suggest the presence of two genetic groups within the population and diversity among the Gaya population. We could not correlate the genomic differences between these two groups with either social status or sex. All the ancient individuals' genomic profiles, including phenotypically relevant SNPs associated with hair and eye color, facial morphology, and myopia, imply strong genetic and phenotypic continuity with modern Koreans for the last 1,700 years.

Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment.

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

Ancient Rome: A genetic crossroads of Europe and the Mediterranean.

Ancient Rome was the capital of an empire of ~70 million inhabitants, but little is known about the genetics of ancient Romans. Here we present 127 genomes from 29 archaeological sites in and around Rome, spanning the past 12,000 years. We observe two major prehistoric ancestry transitions: one with the introduction of farming and another prior to the Iron Age. By the founding of Rome, the genetic composition of the region approximated that of modern Mediterranean populations. During the Imperial period, Rome's population received net immigration from the Near East, followed by an increase in genetic contributions from Europe. These ancestry shifts mirrored the geopolitical affiliations of Rome and were accompanied by marked interindividual diversity, reflecting gene flow from across the Mediterranean, Europe, and North Africa.

Authentication and Assessment of Contamination in Ancient DNA.

Contamination from both present-day humans and postmortem microbial sources is a common challenge in ancient DNA studies. Here we present a suite of tools to assist in the assessment of contamination in ancient DNA data sets. These tools perform standard tests of authenticity of ancient DNA data including detecting the presence of postmortem damage signatures in sequence alignments and quantifying the amount of present-day human contamination.

A fourth Denisovan individual.

The presence of Neandertals in Europe and Western Eurasia before the arrival of anatomically modern humans is well supported by archaeological and paleontological data. In contrast, fossil evidence for Denisovans, a sister group of Neandertals recently identified on the basis of DNA sequences, is limited to three specimens, all of which originate from Denisova Cave in the Altai Mountains (Siberia, Russia). We report the retrieval of DNA from a deciduous lower second molar (Denisova 2), discovered in a deep stratigraphic layer in Denisova Cave, and show that this tooth comes from a female Denisovan individual. On the basis of the number of "missing substitutions" in the mitochondrial DNA determined from the specimen, we find that Denisova 2 is substantially older than two of the other Denisovans, reinforcing the view that Denisovans were likely to have been present in the vicinity of Denisova Cave over an extended time period. We show that the level of nuclear DNA sequence diversity found among Denisovans is within the lower range of that of present-day human populations.

Nuclear and mitochondrial DNA sequences from two Denisovan individuals.

Denisovans, a sister group of Neandertals, have been described on the basis of a nuclear genome sequence from a finger phalanx (Denisova 3) found in Denisova Cave in the Altai Mountains. The only other Denisovan specimen described to date is a molar (Denisova 4) found at the same site. This tooth carries a mtDNA sequence similar to that of Denisova 3. Here we present nuclear DNA sequences from Denisova 4 and a morphological description, as well as mitochondrial and nuclear DNA sequence data, from another molar (Denisova 8) found in Denisova Cave in 2010. This new molar is similar to Denisova 4 in being very large and lacking traits typical of Neandertals and modern humans. Nuclear DNA sequences from the two molars form a clade with Denisova 3. The mtDNA of Denisova 8 is more diverged and has accumulated fewer substitutions than the mtDNAs of the other two specimens, suggesting Denisovans were present in the region over an extended period. The nuclear DNA sequence diversity among the three Denisovans is comparable to that among six Neandertals, but lower than that among present-day humans.

Ancient human genomes suggest three ancestral populations for present-day Europeans.

We sequenced the genomes of a ∼7,000-year-old farmer from Germany and eight ∼8,000-year-old hunter-gatherers from Luxembourg and Sweden. We analysed these and other ancient genomes with 2,345 contemporary humans to show that most present-day Europeans derive from at least three highly differentiated populations: west European hunter-gatherers, who contributed ancestry to all Europeans but not to Near Easterners; ancient north Eurasians related to Upper Palaeolithic Siberians, who contributed to both Europeans and Near Easterners; and early European farmers, who were mainly of Near Eastern origin but also harboured west European hunter-gatherer related ancestry. We model these populations' deep relationships and show that early European farmers had ∼44% ancestry from a 'basal Eurasian' population that split before the diversification of other non-African lineages.

Patterns of coding variation in the complete exomes of three Neandertals.

We present the DNA sequence of 17,367 protein-coding genes in two Neandertals from Spain and Croatia and analyze them together with the genome sequence recently determined from a Neandertal from southern Siberia. Comparisons with present-day humans from Africa, Europe, and Asia reveal that genetic diversity among Neandertals was remarkably low, and that they carried a higher proportion of amino acid-changing (nonsynonymous) alleles inferred to alter protein structure or function than present-day humans. Thus, Neandertals across Eurasia had a smaller long-term effective population than present-day humans. We also identify amino acid substitutions in Neandertals and present-day humans that may underlie phenotypic differences between the two groups. We find that genes involved in skeletal morphology have changed more in the lineage leading to Neandertals than in the ancestral lineage common to archaic and modern humans, whereas genes involved in behavior and pigmentation have changed more on the modern human lineage.

The complete genome sequence of a Neanderthal from the Altai Mountains.

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.

Next-generation museomics disentangles one of the largest primate radiations.

Guenons (tribe Cercopithecini) are one of the most diverse groups of primates. They occupy all of sub-Saharan Africa and show great variation in ecology, behavior, and morphology. This variation led to the description of over 60 species and subspecies. Here, using next-generation DNA sequencing (NGS) in combination with targeted DNA capture, we sequenced 92 mitochondrial genomes from museum-preserved specimens as old as 117 years. We infer evolutionary relationships and estimate divergence times of almost all guenon taxa based on mitochondrial genome sequences. Using this phylogenetic framework, we infer divergence dates and reconstruct ancestral geographic ranges. We conclude that the extraordinary radiation of guenons has been a complex process driven by, among other factors, localized fluctuations of African forest cover. We find incongruences between phylogenetic trees reconstructed from mitochondrial and nuclear DNA sequences, which can be explained by either incomplete lineage sorting or hybridization. Furthermore, having produced the largest mitochondrial DNA data set from museum specimens, we document how NGS technologies can "unlock" museum collections, thereby helping to unravel the tree-of-life.

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.

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.