An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers.

We report an ancient genome from the Indus Valley Civilization (IVC). The individual we sequenced fits as a mixture of people related to ancient Iranians (the largest component) and Southeast Asian hunter-gatherers, a unique profile that matches ancient DNA from 11 genetic outliers from sites in Iran and Turkmenistan in cultural communication with the IVC. These individuals had little if any Steppe pastoralist-derived ancestry, showing that it was not ubiquitous in northwest South Asia during the IVC as it is today. The Iranian-related ancestry in the IVC derives from a lineage leading to early Iranian farmers, herders, and hunter-gatherers before their ancestors separated, contradicting the hypothesis that the shared ancestry between early Iranians and South Asians reflects a large-scale spread of western Iranian farmers east. Instead, sampled ancient genomes from the Iranian plateau and IVC descend from different groups of hunter-gatherers who began farming without being connected by substantial movement of people.

Reconstructing double-stranded DNA fragments on a single-molecule level reveals patterns of degradation in ancient samples.

Extensive manipulations involved in the preparation of DNA samples for sequencing have hitherto made it impossible to determine the precise structure of double-stranded DNA fragments being sequenced, such as the presence of blunt ends, single-stranded overhangs, or single-strand breaks. We here describe MatchSeq, a method that combines single-stranded DNA library preparation from diluted DNA samples with computational sequence matching, allowing the reconstruction of double-stranded DNA fragments on a single-molecule level. The application of MatchSeq to Neanderthal DNA, a particularly complex source of degraded DNA, reveals that 1- or 2-nt overhangs and blunt ends dominate the ends of ancient DNA molecules and that short gaps exist, which are predominantly caused by the loss of individual purines. We further show that deamination of cytosine to uracil occurs in both single- and double-stranded contexts close to the ends of molecules, and that single-stranded parts of DNA fragments are enriched in pyrimidines. MatchSeq provides unprecedented resolution for interrogating the structures of fragmented double-stranded DNA and can be applied to fragmented double-stranded DNA isolated from any biological source. The method relies on well-established laboratory techniques and can easily be integrated into routine data generation. This possibility is shown by the successful reconstruction of double-stranded DNA fragments from previously published single-stranded sequence data, allowing a more comprehensive characterization of the biochemical properties not only of ancient DNA but also of cell-free DNA from human blood plasma, a clinically relevant marker for the diagnosis and monitoring of disease.

Systematic benchmark of ancient DNA read mapping.

The current standard practice for assembling individual genomes involves mapping millions of short DNA sequences (also known as DNA 'reads') against a pre-constructed reference genome. Mapping vast amounts of short reads in a timely manner is a computationally challenging task that inevitably produces artefacts, including biases against alleles not found in the reference genome. This reference bias and other mapping artefacts are expected to be exacerbated in ancient DNA (aDNA) studies, which rely on the analysis of low quantities of damaged and very short DNA fragments (~30-80 bp). Nevertheless, the current gold-standard mapping strategies for aDNA studies have effectively remained unchanged for nearly a decade, during which time new software has emerged. In this study, we used simulated aDNA reads from three different human populations to benchmark the performance of 30 distinct mapping strategies implemented across four different read mapping software-BWA-aln, BWA-mem, NovoAlign and Bowtie2-and quantified the impact of reference bias in downstream population genetic analyses. We show that specific NovoAlign, BWA-aln and BWA-mem parameterizations achieve high mapping precision with low levels of reference bias, particularly after filtering out reads with low mapping qualities. However, unbiased NovoAlign results required the use of an IUPAC reference genome. While relevant only to aDNA projects where reference population data are available, the benefit of using an IUPAC reference demonstrates the value of incorporating population genetic information into the aDNA mapping process, echoing recent results based on graph genome representations.

Integration of ancient DNA with transdisciplinary dataset finds strong support for Inca resettlement in the south Peruvian coast.

Ancient DNA (aDNA) analysis provides a powerful means of investigating human migration, social organization, and a plethora of other crucial questions about humanity's past. Recently, specialists have suggested that the ideal research design involving aDNA would include multiple independent lines of evidence. In this paper, we adopt a transdisciplinary approach integrating aDNA with archaeological, biogeochemical, and historical data to investigate six individuals found in two cemeteries that date to the Late Horizon (1400 to 1532 CE) and Colonial (1532 to 1825 CE) periods in the Chincha Valley of southern Peru. Genomic analyses indicate that these individuals are genetically most similar to ancient and present-day populations from the north Peruvian coast located several hundred kilometers away. These genomic data are consistent with 16th century written records as well as ceramic, textile, and isotopic data. These results provide some of the strongest evidence yet of state-sponsored resettlement in the pre-Colonial Andes. This study highlights the power of transdisciplinary research designs when using aDNA data and sets a methodological standard for investigating ancient mobility in complex societies.

Genetic contributions to variation in human stature in prehistoric Europe.

The relative contributions of genetics and environment to temporal and geographic variation in human height remain largely unknown. Ancient DNA has identified changes in genetic ancestry over time, but it is not clear whether those changes in ancestry are associated with changes in height. Here, we directly test whether changes over the past 38,000 y in European height predicted using DNA from 1,071 ancient individuals are consistent with changes observed in 1,159 skeletal remains from comparable populations. We show that the observed decrease in height between the Early Upper Paleolithic and the Mesolithic is qualitatively predicted by genetics. Similarly, both skeletal and genetic height remained constant between the Mesolithic and Neolithic and increased between the Neolithic and Bronze Age. Sitting height changes much less than standing height-consistent with genetic predictions-although genetics predicts a small post-Neolithic increase that is not observed in skeletal remains. Geographic variation in stature is also qualitatively consistent with genetic predictions, particularly with respect to latitude. Finally, we hypothesize that an observed decrease in genetic heel bone mineral density in the Neolithic reflects adaptation to the decreased mobility indicated by decreased femoral bending strength. This study provides a model for interpreting phenotypic changes predicted from ancient DNA and demonstrates how they can be combined with phenotypic measurements to understand the relative contribution of genetic and developmentally plastic responses to environmental change.

Ancient Genomes Reveal the Evolutionary History and Origin of Cashmere-Producing Goats in China.

Goats are one of the most widespread farmed animals across the world; however, their migration route to East Asia and local evolutionary history remain poorly understood. Here, we sequenced 27 ancient Chinese goat genomes dating from the Late Neolithic period to the Iron Age. We found close genetic affinities between ancient and modern Chinese goats, demonstrating their genetic continuity. We found that Chinese goats originated from the eastern regions around the Fertile Crescent, and we estimated that the ancestors of Chinese goats diverged from this population in the Chalcolithic period. Modern Chinese goats were divided into a northern and a southern group, coinciding with the most prominent climatic division in China, and two genes related to hair follicle development, FGF5 and EDA2R, were highly divergent between these populations. We identified a likely causal de novo deletion near FGF5 in northern Chinese goats that increased to high frequency over time, whereas EDA2R harbored standing variation dating to the Neolithic. Our findings add to our understanding of the genetic composition and local evolutionary process of Chinese goats.

Testing a series of modifications on genomic library preparation methods for ancient or degraded DNA.

Technological advancements have revolutionized ancient and degraded DNA analysis, moving the field to the Next Generation Sequencing era. One of the advancements, the ancient DNA-oriented high-throughput library preparation methods, enabled the sequencing of more endogenous molecules. Although fairly optimized, both single- and double-stranded library preparation methods hold the potential for further improvement. Here, we test a series of modifications made at different steps of both single- and double-stranded library preparation methods. Given all the modifications tested, we found that two of them provide further benefits, including the use of Endonuclease VIII as a pre-treatment step before preparing single-stranded libraries and the use of a modified second adapter of the single stranded-libraries as an alternative option to enable sequencing of single stranded-libraries with the standard Illumina sequencing primer instead of the custom designed as described in the single stranded library preparation method. Furthermore, we propose uracil-DNA-glycosylase (UDG) could also be considered for both single- and double-stranded library preparation methods, although additional parameters should be taken into account depending on the sequencing strategy and the sample characteristics. Further modifications were also tested and although they were not advantageous, they could be considered as equivalent to the published options.

Mining ancient microbiomes using selective enrichment of damaged DNA molecules.

BACKGROUND: The identification of bona fide microbial taxa in microbiomes derived from ancient and historical samples is complicated by the unavoidable mixture between DNA from ante- and post-mortem microbial colonizers. One possibility to distinguish between these sources of microbial DNA is querying for the presence of age-associated degradation patterns typical of ancient DNA (aDNA). The presence of uracils, resulting from cytosine deamination, has been detected ubiquitously in aDNA retrieved from diverse sources, and used as an authentication criterion. Here, we employ a library preparation method that separates molecules that carry uracils from those that do not for a set of samples that includes Neandertal remains, herbarium specimens and archaeological plant remains. RESULTS: We show that sequencing DNA libraries enriched in molecules carrying uracils effectively amplifies age associated degradation patterns in microbial mixtures of ancient and historical origin. This facilitates the discovery of authentic ancient microbial taxa in cases where degradation patterns are difficult to detect due to large sequence divergence in microbial mixtures. Additionally, the relative enrichment of taxa in the uracil enriched fraction can help to identify bona fide ancient microbial taxa that could be missed using a more targeted approach. CONCLUSIONS: Our experiments show, that in addition to its use in enriching authentic endogenous DNA of organisms of interest, the selective enrichment of damaged DNA molecules can be a valuable tool in the discovery of ancient microbial taxa.

The paradox of HBV evolution as revealed from a 16th century mummy.

Hepatitis B virus (HBV) is a ubiquitous viral pathogen associated with large-scale morbidity and mortality in humans. However, there is considerable uncertainty over the time-scale of its origin and evolution. Initial shotgun data from a mid-16th century Italian child mummy, that was previously paleopathologically identified as having been infected with Variola virus (VARV, the agent of smallpox), showed no DNA reads for VARV yet did for hepatitis B virus (HBV). Previously, electron microscopy provided evidence for the presence of VARV in this sample, although similar analyses conducted here did not reveal any VARV particles. We attempted to enrich and sequence for both VARV and HBV DNA. Although we did not recover any reads identified as VARV, we were successful in reconstructing an HBV genome at 163.8X coverage. Strikingly, both the HBV sequence and that of the associated host mitochondrial DNA displayed a nearly identical cytosine deamination pattern near the termini of DNA fragments, characteristic of an ancient origin. In contrast, phylogenetic analyses revealed a close relationship between the putative ancient virus and contemporary HBV strains (of genotype D), at first suggesting contamination. In addressing this paradox we demonstrate that HBV evolution is characterized by a marked lack of temporal structure. This confounds attempts to use molecular clock-based methods to date the origin of this virus over the time-frame sampled so far, and means that phylogenetic measures alone cannot yet be used to determine HBV sequence authenticity. If genuine, this phylogenetic pattern indicates that the genotypes of HBV diversified long before the 16th century, and enables comparison of potential pathogenic similarities between modern and ancient HBV. These results have important implications for our understanding of the emergence and evolution of this common viral pathogen.

A new model for ancient DNA decay based on paleogenomic meta-analysis.

The persistence of DNA over archaeological and paleontological timescales in diverse environments has led to a revolutionary body of paleogenomic research, yet the dynamics of DNA degradation are still poorly understood. We analyzed 185 paleogenomic datasets and compared DNA survival with environmental variables and sample ages. We find cytosine deamination follows a conventional thermal age model, but we find no correlation between DNA fragmentation and sample age over the timespans analyzed, even when controlling for environmental variables. We propose a model for ancient DNA decay wherein fragmentation rapidly reaches a threshold, then subsequently slows. The observed loss of DNA over time may be due to a bulk diffusion process in many cases, highlighting the importance of tissues and environments creating effectively closed systems for DNA preservation. This model of DNA degradation is largely based on mammal bone samples due to published genomic dataset availability. Continued refinement to the model to reflect diverse biological systems and tissue types will further improve our understanding of ancient DNA breakdown dynamics.

Cytochrome C Oxidase Subunit 1, Internal Transcribed Spacer 1, Nicotinamide Adenine Dinucleotide Hydrogen Dehydrogenase Subunits 2 and 5 of Clonorchis sinensis Ancient DNA Retrieved from Joseon Dynasty Mummy Specimens.

We analyzed Clonorchis sinensis ancient DNA (aDNA) acquired from the specimens of the Joseon mummies. The target regions were cytochrome C oxidase subunit 1 (CO1), internal transcribed spacer 1 (ITS1), nicotinamide adenine dinucleotide hydrogen (NADH) dehydrogenase subunits 2 (NAD2) and 5 (NAD5). The sequences of C. sinensis aDNA was completely or almost identical to modern C. sinensis sequences in GenBank. We also found that ITS1, NAD2 and NAD5 could be good markers for molecular diagnosis between C. sinensis and the other trematode parasite species. The current result could improve our knowledge about genetic history of C. sinensis.

gargammel: a sequence simulator for ancient DNA.

Summary: Ancient DNA has emerged as a remarkable tool to infer the history of extinct species and past populations. However, many of its characteristics, such as extensive fragmentation, damage and contamination, can influence downstream analyses. To help investigators measure how these could impact their analyses in silico , we have developed gargammel, a package that simulates ancient DNA fragments given a set of known reference genomes. Our package simulates the entire molecular process from post-mortem DNA fragmentation and DNA damage to experimental sequencing errors, and reproduces most common bias observed in ancient DNA datasets. Availability and Implementation: The package is publicly available on github: https://grenaud.github.io/gargammel/ and released under the GPL. Contact: gabriel.renaud@snm.ku.dk. Supplementary information: Supplementary data are available at Bioinformatics online.

High-Throughput DNA sequencing of ancient wood.

Reconstructing the colonization and demographic dynamics that gave rise to extant forests is essential to forecasts of forest responses to environmental changes. Classical approaches to map how population of trees changed through space and time largely rely on pollen distribution patterns, with only a limited number of studies exploiting DNA molecules preserved in wooden tree archaeological and subfossil remains. Here, we advance such analyses by applying high-throughput (HTS) DNA sequencing to wood archaeological and subfossil material for the first time, using a comprehensive sample of 167 European white oak waterlogged remains spanning a large temporal (from 550 to 9,800 years) and geographical range across Europe. The successful characterization of the endogenous DNA and exogenous microbial DNA of 140 (~83%) samples helped the identification of environmental conditions favouring long-term DNA preservation in wood remains, and started to unveil the first trends in the DNA decay process in wood material. Additionally, the maternally inherited chloroplast haplotypes of 21 samples from three periods of forest human-induced use (Neolithic, Bronze Age and Middle Ages) were found to be consistent with those of modern populations growing in the same geographic areas. Our work paves the way for further studies aiming at using ancient DNA preserved in wood to reconstruct the micro-evolutionary response of trees to climate change and human forest management.

ProtASR: An Evolutionary Framework for Ancestral Protein Reconstruction with Selection on Folding Stability.

The computational reconstruction of ancestral proteins provides information on past biological events and has practical implications for biomedicine and biotechnology. Currently available tools for ancestral sequence reconstruction (ASR) are often based on empirical amino acid substitution models that assume that all sites evolve at the same rate and under the same process. However, this assumption is frequently violated because protein evolution is highly heterogeneous due to different selective constraints among sites. Here, we present ProtASR, a new evolutionary framework to infer ancestral protein sequences accounting for selection on protein stability. First, ProtASR generates site-specific substitution matrices through the structurally constrained mean-field (MF) substitution model, which considers both unfolding and misfolding stability. We previously showed that MF models outperform empirical amino acid substitution models, as well as other structurally constrained substitution models, both in terms of likelihood and correctly inferring amino acid distributions across sites. In the second step, ProtASR adapts a well-established maximum-likelihood (ML) ASR procedure to infer ancestral proteins under MF models. A known bias of ML ASR methods is that they tend to overestimate the stability of ancestral proteins by underestimating the frequency of deleterious mutations. We compared ProtASR under MF to two empirical substitution models (JTT and CAT), reconstructing the ancestral sequences of simulated proteins. ProtASR yields reconstructed proteins with less biased stabilities, which are significantly closer to those of the simulated proteins. Analysis of extant protein families suggests that folding stability evolves through time across protein families, potentially reflecting neutral fluctuation. Some families exhibit a more constant protein folding stability, while others are more variable. ProtASR is freely available from https://github.com/miguelarenas/protasr and includes detailed documentation and ready-to-use examples. It runs in seconds/minutes depending on protein length and alignment size. [Ancestral sequence reconstruction; folding stability; molecular adaptation; phylogenetics; protein evolution; protein structure.].

Ancient plant DNA in lake sediments.

Contents 924 I. 925 II. 925 III. 927 IV. 929 V. 930 VI. 930 VII. 931 VIII. 933 IX. 935 X. 936 XI. 938 938 References 938 SUMMARY: Recent advances in sequencing technologies now permit the analyses of plant DNA from fossil samples (ancient plant DNA, plant aDNA), and thus enable the molecular reconstruction of palaeofloras. Hitherto, ancient frozen soils have proved excellent in preserving DNA molecules, and have thus been the most commonly used source of plant aDNA. However, DNA from soil mainly represents taxa growing a few metres from the sampling point. Lakes have larger catchment areas and recent studies have suggested that plant aDNA from lake sediments is a more powerful tool for palaeofloristic reconstruction. Furthermore, lakes can be found globally in nearly all environments, and are therefore not limited to perennially frozen areas. Here, we review the latest approaches and methods for the study of plant aDNA from lake sediments and discuss the progress made up to the present. We argue that aDNA analyses add new and additional perspectives for the study of ancient plant populations and, in time, will provide higher taxonomic resolution and more precise estimation of abundance. Despite this, key questions and challenges remain for such plant aDNA studies. Finally, we provide guidelines on technical issues, including lake selection, and we suggest directions for future research on plant aDNA studies in lake sediments.

Ancient mitochondrial genome reveals trace of prehistoric migration in the east Pamir by pastoralists.

The complete mitochondrial genome of one 700-year-old individual found in Tashkurgan, Xinjiang was target enriched and sequenced in order to shed light on the population history of Tashkurgan and determine the phylogenetic relationship of haplogroup U5a. The ancient sample was assigned to a subclade of haplogroup U5a2a1, which is defined by two rare and stable transversions at 16114A and 13928C. Phylogenetic analysis shows a distribution pattern for U5a2a that is indicative of an origin in the Volga-Ural region and exhibits a clear eastward geographical expansion that correlates with the pastoral culture also entering the Eurasian steppe. The haplogroup U5a2a present in the ancient Tashkurgan individual reveals prehistoric migration in the East Pamir by pastoralists. This study shows that studying an ancient mitochondrial genome is a useful approach for studying the evolutionary process and population history of Eastern Pamir.

Ancient DNA in historical parchments - identifying a procedure for extraction and amplification of genetic material.

Historical parchments in the form of documents, manuscripts, books, or letters, make up a large portion of cultural heritage collections. Their priceless historical value is associated with not only their content, but also the information hidden in the DNA deposited on them. Analyses of ancient DNA (aDNA) retrieved from parchments can be used in various investigations, including, but not limited to, studying their authentication, tracing the development of the culture, diplomacy, and technology, as well as obtaining information on the usage and domestication of animals. This article proposes and verifies a procedure for aDNA recovery from historical parchments and its appropriate preparation for further analyses. This study involved experimental selection of an aDNA extraction method with the highest efficiency and quality of extracted genetic material, from among the multi-stage phenol-chloroform extraction methods, and the modern, column-based techniques that use selective DNA-binding membranes. Moreover, current techniques to amplify entire genetic material were questioned, and the possibility of using mitochondrial DNA for species identification was analyzed. The usefulness of the proposed procedure was successfully confirmed in identification tests of historical parchments dating back to the 13-16th century AD.

Ancient DNA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean.

Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since at least the 16th century. Here we build upon previous work on ancient DNA of BFT in the Mediterranean by comparing contemporary (2009-2012) specimens with archival (1911-1926) and archaeological (2nd century BCE-15th century CE) specimens that represent population states prior to these two major periods of exploitation, respectively. We successfully genotyped and analysed 259 contemporary and 123 historical (91 archival and 32 archaeological) specimens at 92 SNP loci that were selected for their ability to differentiate contemporary populations or their association with core biological functions. We found no evidence of genetic bottlenecks, inbreeding or population restructuring between temporal sample groups that might explain what has driven catch fluctuations since the 16th century. We also detected a putative adaptive response, involving the cytoskeletal protein synemin which may be related to muscle stress. However, these results require further investigation with more extensive genome-wide data to rule out demographic changes due to overfishing, and other natural and anthropogenic factors, in addition to elucidating the adaptive drivers related to these.

Bile acids and oxo-metabolites as markers of human faecal input in the ancient Pompeii ruins.

Small organic molecules, lipids, proteins, and DNA fragments can remain stable over centuries. Powerful and sensitive chemical analysis can therefore be used to characterize ancient remains for classical archaeological studies. This bio-ecological dimension of archaeology can contribute knowledge about several aspects of ancient life, including social organization, daily habits, nutrition, and food storage. Faecal remains (i.e. coprolites) are particularly interesting in this regard, with scientists seeking to identify new faecal markers. Here, we report the analysis of faecal samples from modern-day humans and faecal samples from a discharge pit on the site of the ruins of ancient Pompeii. We propose that bile acids and their gut microbiota oxo-metabolites are the most specific steroid markers for detecting faecal inputs. This is due to their extreme chemical stability and their exclusive occurrence in vertebrate faeces, compared to other ubiquitous sterols and steroids.

Chronic active non-lethal human-type tuberculosis in a high royal Bavarian officer of Napoleonic times-a mummy study.

In paleopathology, morphological and molecular evidence for infection by mycobacteria of the M. tuberculosis complex (MTC) is frequently associated with early death. In the present report, we describe a multidisciplinary study of a well-preserved mummy from Napoleonic times with a long-standing tuberculous infection by M. tuberculosis senso stricto who died at the age of 88 years of focal and non-MTB related bronchopneumonia. The well-preserved natural mummy of the Royal Bavarian General, Count Heinrich LII Reuss-Köstritz (1763-1851 CE), was extensively investigated by macro- and histomorphology, whole body CT scans and organ radiography, various molecular tissue analyses, including stable isotope analysis and molecular genetic tests. We identified signs for a long-standing, but terminally inactive pulmonary tuberculosis, tuberculous destruction of the second lumbar vertebral body, and a large tuberculous abscess in the right (retroperitoneal) psoas region (a cold abscess). This cold abscess harboured an active tuberculous infection as evidenced by histological and molecular tests. Radiological and histological analysis further revealed extensive arteriosclerosis with (non-obliterating) coronary and significant carotid arteriosclerosis, healthy bone tissue without evidence of age-related osteopenia, evidence for diffuse idiopathic skeletal hyperostosis and mild osteoarthrosis of few joints. This suggests excellent living conditions correlating well with his diet indicated by stable isotope results and literary evidence. Despite the clear evidence of a tuberculous cold abscess with bacterioscopic and molecular proof for a persisting MTC infection of a human-type M. tuberculosis strain, we can exclude the chronic MTC infection as cause of death. The detection of MTC in historic individuals should therefore be interpreted with great caution and include further data, such as their nutritional status.