Extraction of highly degraded DNA from ancient bones, teeth and sediments for high-throughput sequencing

The influence of sample quantity and lysis parameters on the success of ancient DNA extraction from skeletal remains

DNA extraction is of utmost importance in archaeobiology, as it determines the success of further DNA analyses. This study concentrates on the success of ancient DNA extraction using silica spin columns and PCR-based analysis from archaeological skeletal material and investigates the influence of sample quantity, lysis time and lysis temperature during sample preparation. The results show that lysis times ranging from 2 to 48 h are suitable, and that lysis should be carried out at a constant temperature of 56°C. Concerning sample quantity, 10 mg for mitochondrial DNA and 50 mg for chromosomal DNA are sufficient for high quality analyses. Thus invaluable sample material can be saved, and time of sample preparation can be reduced considerably.

Human mobility at Tell Atchana (Alalakh), Hatay, Turkey during the 2nd millennium BC: Integration of isotopic and genomic evidence

The Middle and Late Bronze Age, a period roughly spanning the 2nd millennium BC (ca. 2000-1200 BC) in the Near East, is frequently referred to as the first 'international age', characterized by intense and far-reaching contacts between different entities from the eastern Mediterranean to the Near East and beyond. In a large-scale tandem study of stable isotopes and ancient DNA of individuals excavated at Tell Atchana (Alalakh, located in Hatay, Turkey), we explored the role of mobility at the capital of a regional kingdom, named Mukish during the Late Bronze Age, which spanned the Amuq Valley and some areas beyond. We generated strontium and oxygen isotope data from dental enamel for 53 individuals and 77 individuals, respectively, and added ancient DNA data of 10 newly sequenced individuals to a dataset of 27 individuals published in 2020. Additionally, we improved the DNA coverage of one individual from this 2020 dataset. The DNA data revealed a very homogeneous gene pool. This picture of an overwhelmingly local ancestry was consistent with the evidence of local upbringing in most of the individuals indicated by the isotopic data, where only five were found to be non-local. High levels of contact, trade, and exchange of ideas and goods in the Middle and Late Bronze Ages, therefore, seem not to have translated into high levels of individual mobility detectable at Tell Atchana.

A method for the temperature-controlled extraction of DNA from ancient bones

Contamination with microbial and other exogenous DNA poses a significant challenge in the generation of genome-wide sequence data from ancient skeletal remains. Here we describe a method for separating ancient DNA into multiple fractions during DNA extraction by sequential temperature-controlled release of DNA into sodium phosphate buffer. An evaluation of the effectiveness of the method using a set of three ancient bones resulted in between 1.6- and 32-fold enrichment of endogenous DNA compared with regular DNA extraction. For two bones, the method outperformed previous methods of decontaminating ancient bones, including hypochlorite treatment, which resulted in near-complete destruction of DNA in the worst-preserved sample. This extraction method expands the spectrum of methods available for depleting contaminant DNA from ancient skeletal remains.

Pleistocene sediment DNA reveals hominin and faunal turnovers at Denisova Cave

Denisova Cave in southern Siberia is the type locality of the Denisovans, an archaic hominin group who were related to Neanderthals^1–4. The dozen hominin remains recovered from the deposits also include Neanderthals^5,6 and the child of a Neanderthal and a Denisovan⁷, which suggests that Denisova Cave was a contact zone between these archaic hominins. However, uncertainties persist about the order in which these groups appeared at the site, the timing and environmental context of hominin occupation, and the association of particular hominin groups with archaeological assemblages^5,8–11. Here we report the analysis of DNA from 728 sediment samples that were collected in a grid-like manner from layers dating to the Pleistocene epoch. We retrieved ancient faunal and hominin mitochondrial (mt)DNA from 685 and 175 samples, respectively. The earliest evidence for hominin mtDNA is of Denisovans, and is associated with early Middle Palaeolithic stone tools that were deposited approximately 250,000 to 170,000 years ago; Neanderthal mtDNA first appears towards the end of this period. We detect a turnover in the mtDNA of Denisovans that coincides with changes in the composition of faunal mtDNA, and evidence that Denisovans and Neanderthals occupied the site repeatedly—possibly until, or after, the onset of the Initial Upper Palaeolithic at least 45,000 years ago, when modern human mtDNA is first recorded in the sediments.

Ancient mitochondrial DNA from sediments reveals the sequence of Denisovan, Neanderthal and faunal occupation of Denisova Cave, and evidence for the appearance of modern humans at least 45,000 years ago.

A Fast and Efficient Single-stranded Genomic Library Preparation Method Optimized for Ancient DNA

We present a protocol to prepare extracted DNA for sequencing on the Illumina sequencing platform that has been optimized for ancient and degraded DNA. Our approach, the Santa Cruz Reaction or SCR, uses directional splinted ligation of Illumina’s P5 and P7 adapters to convert natively single-stranded DNA and heat denatured double-stranded DNA into sequencing libraries in a single enzymatic reaction. To demonstrate its efficacy in converting degraded DNA molecules, we prepare 5 ancient DNA extracts into sequencing libraries using the SCR and 2 of the most commonly used approaches for preparing degraded DNA for sequencing: BEST, which targets and converts double-stranded DNA, and ssDNA2.0, which targets and converts single-stranded DNA. We then compare the efficiency with which each approach recovers unique molecules, or library complexity, given a standard amount of DNA input. We find that the SCR consistently outperforms the BEST protocol in recovering unique molecules and, despite its relative simplicity to perform and low cost per library, has similar performance to ssDNA2.0 across a wide range of DNA inputs. The SCR is a cost- and time-efficient approach that minimizes the loss of unique molecules and makes accessible a taxonomically, geographically, and a temporally broader sample of preserved remains for genomic analysis.

Earliest known human burial in Africa

The origin and evolution of hominin mortuary practices are topics of intense interest and debate^1-3. Human burials dated to the Middle Stone Age (MSA) are exceedingly rare in Africa and unknown in East Africa^1-6. Here we describe the partial skeleton of a roughly 2.5- to 3.0-year-old child dating to 78.3 ± 4.1 thousand years ago, which was recovered in the MSA layers of Panga ya Saidi (PYS), a cave site in the tropical upland coast of Kenya^7,8. Recent excavations have revealed a pit feature containing a child in a flexed position. Geochemical, granulometric and micromorphological analyses of the burial pit content and encasing archaeological layers indicate that the pit was deliberately excavated. Taphonomical evidence, such as the strict articulation or good anatomical association of the skeletal elements and histological evidence of putrefaction, support the in-place decomposition of the fresh body. The presence of little or no displacement of the unstable joints during decomposition points to an interment in a filled space (grave earth), making the PYS finding the oldest known human burial in Africa. The morphological assessment of the partial skeleton is consistent with its assignment to Homo sapiens, although the preservation of some primitive features in the dentition supports increasing evidence for non-gradual assembly of modern traits during the emergence of our species. The PYS burial sheds light on how MSA populations interacted with the dead.

Comparative analysis of DNA extraction protocols for ancient soft tissue museum samples

DNA studies of endangered or extinct species often rely on ancient or degraded remains. The majority of ancient DNA (aDNA) extraction protocols focus on skeletal elements, with skin and hair samples rarely explored. Similar to that found in bones and teeth, DNA extracted from historical or ancient skin and fur samples is also extremely fragmented with low endogenous content due to natural degradation processes. Thus, the development of effective DNA extraction methods is required for these materials. Here, we compared the performance of two DNA extraction protocols (commercial and custom laboratory aDNA methods) on hair and skin samples from decades-old museum specimens to Iron Age archaeological material. We found that apart from the impact sample-specific taphonomic and handling history has on the quantity and quality of DNA preservation, skin yielded more endogenous DNA than hair of the samples and protocols tested. While both methods recovered DNA from ancient soft tissue, the laboratory method performed better overall in terms of DNA yield and quality, which was primarily due to the poorer performance of the commercial binding buffer in recovering aDNA.

Unearthing Neanderthal population history using nuclear and mitochondrial DNA from cave sediments

Bones and teeth are important sources of Pleistocene hominin DNA, but are rarely recovered at archaeological sites. Mitochondrial DNA (mtDNA) has been retrieved from cave sediments but provides limited value for studying population relationships. We therefore developed methods for the enrichment and analysis of nuclear DNA from sediments and applied them to cave deposits in western Europe and southern Siberia dated to between 200,000 and 50,000 years ago. We detected a population replacement in northern Spain about 100,000 years ago, which was accompanied by a turnover of mtDNA. We also identified two radiation events in Neanderthal history during the early part of the Late Pleistocene. Our work lays the ground for studying the population history of ancient hominins from trace amounts of nuclear DNA in sediments.

Heterogeneous Hunter-Gatherer and Steppe-Related Ancestries in Late Neolithic and Bell Beaker Genomes from Present-Day France

The transition from the Late Neolithic to the Bronze Age has witnessed important population and societal changes in western Europe.¹ These include massive genomic contributions of pastoralist herders originating from the Pontic-Caspian steppes^2,3 into local populations, resulting from complex interactions between collapsing hunter-gatherers and expanding farmers of Anatolian ancestry.^4-8 This transition is documented through extensive ancient genomic data from present-day Britain,^9,10 Ireland,^11,12 Iberia,¹³ Mediterranean islands,^14,15 and Germany.⁸ It remains, however, largely overlooked in France, where most focus has been on the Middle Neolithic (n = 63),^8,9,16 with the exception of one Late Neolithic genome sequenced at 0.05× coverage.¹⁶ This leaves the key transitional period covering ∼3,400-2,700 cal. years (calibrated years) BCE genetically unsampled and thus the exact time frame of hunter-gatherer persistence and arrival of steppe migrations unknown. To remediate this, we sequenced 24 ancient human genomes from France spanning ∼3,400-1,600 cal. years BCE. This reveals Late Neolithic populations that are genetically diverse and include individuals with dark skin, hair, and eyes. We detect heterogeneous hunter-gatherer ancestries within Late Neolithic communities, reaching up to ∼63.3% in some individuals, and variable genetic contributions of steppe herders in Bell Beaker populations. We provide an estimate as late as ∼3,800 years BCE for the admixture between Neolithic and Mesolithic populations and as early as ∼2,650 years BCE for the arrival of steppe-related ancestry. The genomic heterogeneity characterized underlines the complex history of human interactions even at the local scale.

Reevaluating the timing of Neanderthal disappearance in Northwest Europe

Elucidating when Neanderthal populations disappeared from Eurasia is a key question in paleoanthropology, and Belgium is one of the key regions for studying the Middle to Upper Paleolithic transition. Previous radiocarbon dating placed the Spy Neanderthals among the latest surviving Neanderthals in Northwest Europe with reported dates as young as 23,880 ± 240 B.P. (OxA-8912). Questions were raised, however, regarding the reliability of these dates. Soil contamination and carbon-based conservation products are known to cause problems during the radiocarbon dating of bulk collagen samples. Employing a compound-specific approach that is today the most efficient in removing contamination and ancient genomic analysis, we demonstrate here that previous dates produced on Neanderthal specimens from Spy were inaccurately young by up to 10,000 y due to the presence of unremoved contamination. Our compound-specific radiocarbon dates on the Neanderthals from Spy and those from Engis and Fonds-de-Forêt demonstrate that they disappeared from Northwest Europe at 44,200 to 40,600 cal B.P. (at 95.4% probability), much earlier than previously suggested. Our data contribute significantly to refining models for Neanderthal disappearance in Europe and, more broadly, show that chronometric models regarding the appearance or disappearance of animal or hominin groups should be based only on radiocarbon dates obtained using robust pretreatment methods.

Genomic insights into the formation of human populations in East Asia

The deep population history of East Asia remains poorly understood owing to a lack of ancient DNA data and sparse sampling of present-day people1,2. Here we report genome-wide data from 166 East Asian individuals dating to between 6000 BC and AD 1000 and 46 present-day groups. Hunter-gatherers from Japan, the Amur River Basin, and people of Neolithic and Iron Age Taiwan and the Tibetan Plateau are linked by a deeply splitting lineage that probably reflects a coastal migration during the Late Pleistocene epoch. We also follow expansions during the subsequent Holocene epoch from four regions. First, hunter-gatherers from Mongolia and the Amur River Basin have ancestry shared by individuals who speak Mongolic and Tungusic languages, but do not carry ancestry characteristic of farmers from the West Liao River region (around 3000 BC), which contradicts theories that the expansion of these farmers spread the Mongolic and Tungusic proto-languages. Second, farmers from the Yellow River Basin (around 3000 BC) probably spread Sino-Tibetan languages, as their ancestry dispersed both to Tibet-where it forms approximately 84% of the gene pool in some groups-and to the Central Plain, where it has contributed around 59-84% to modern Han Chinese groups. Third, people from Taiwan from around 1300 BC to AD 800 derived approximately 75% of their ancestry from a lineage that is widespread in modern individuals who speak Austronesian, Tai-Kadai and Austroasiatic languages, and that we hypothesize derives from farmers of the Yangtze River Valley. Ancient people from Taiwan also derived about 25% of their ancestry from a northern lineage that is related to, but different from, farmers of the Yellow River Basin, which suggests an additional north-to-south expansion. Fourth, ancestry from Yamnaya Steppe pastoralists arrived in western Mongolia after around 3000 BC but was displaced by previously established lineages even while it persisted in western China, as would be expected if this ancestry was associated with the spread of proto-Tocharian Indo-European languages. Two later gene flows affected western Mongolia: migrants after around 2000 BC with Yamnaya and European farmer ancestry, and episodic influences of later groups with ancestry from Turan.

A minimally destructive protocol for DNA extraction from ancient teeth

Ancient DNA sampling methods-although optimized for efficient DNA extraction-are destructive, relying on drilling or cutting and powdering (parts of) bones and teeth. As the field of ancient DNA has grown, so have concerns about the impact of destructive sampling of the skeletal remains from which ancient DNA is obtained. Due to a particularly high concentration of endogenous DNA, the cementum of tooth roots is often targeted for ancient DNA sampling, but destructive sampling methods of the cementum often result in the loss of at least one entire root. Here, we present a minimally destructive method for extracting ancient DNA from dental cementum present on the surface of tooth roots. This method does not require destructive drilling or grinding, and, following extraction, the tooth remains safe to handle and suitable for most morphological studies, as well as other biochemical studies, such as radiocarbon dating. We extracted and sequenced ancient DNA from 30 teeth (and nine corresponding petrous bones) using this minimally destructive extraction method in addition to a typical tooth sampling method. We find that the minimally destructive method can provide ancient DNA that is of comparable quality to extracts produced from teeth that have undergone destructive sampling processes. Further, we find that a rigorous cleaning of the tooth surface combining diluted bleach and UV light irradiation seems sufficient to minimize external contaminants usually removed through the physical removal of a superficial layer when sampling through regular powdering methods.

A genetic history of the pre-contact Caribbean

Humans settled the Caribbean about 6,000 years ago, and ceramic use and intensified agriculture mark a shift from the Archaic to the Ceramic Age at around 2,500 years ago1-3. Here we report genome-wide data from 174 ancient individuals from The Bahamas, Haiti and the Dominican Republic (collectively, Hispaniola), Puerto Rico, Curaçao and Venezuela, which we co-analysed with 89 previously published ancient individuals. Stone-tool-using Caribbean people, who first entered the Caribbean during the Archaic Age, derive from a deeply divergent population that is closest to Central and northern South American individuals; contrary to previous work4, we find no support for ancestry contributed by a population related to North American individuals. Archaic-related lineages were >98% replaced by a genetically homogeneous ceramic-using population related to speakers of languages in the Arawak family from northeast South America; these people moved through the Lesser Antilles and into the Greater Antilles at least 1,700 years ago, introducing ancestry that is still present. Ancient Caribbean people avoided close kin unions despite limited mate pools that reflect small effective population sizes, which we estimate to be a minimum of 500-1,500 and a maximum of 1,530-8,150 individuals on the combined islands of Puerto Rico and Hispaniola in the dozens of generations before the individuals who we analysed lived. Census sizes are unlikely to be more than tenfold larger than effective population sizes, so previous pan-Caribbean estimates of hundreds of thousands of people are too large5,6. Confirming a small and interconnected Ceramic Age population7, we detect 19 pairs of cross-island cousins, close relatives buried around 75 km apart in Hispaniola and low genetic differentiation across islands. Genetic continuity across transitions in pottery styles reveals that cultural changes during the Ceramic Age were not driven by migration of genetically differentiated groups from the mainland, but instead reflected interactions within an interconnected Caribbean world1,8.

Evaluation of DNA Extraction Methods Developed for Forensic and Ancient DNA Applications Using Bone Samples of Different Age

The efficient extraction of DNA from challenging samples, such as bones, is critical for the success of downstream genotyping analysis in molecular genetic disciplines. Even though the ancient DNA community has developed several protocols targeting small DNA fragments that are typically present in decomposed or old specimens, only recently forensic geneticists have started to adopt those protocols. Here, we compare an ancient DNA extraction protocol (Dabney) with a bone extraction method (Loreille) typically used in forensics. Real-time quantitative PCR and forensically representative typing methods including fragment size analysis and sequencing were used to assess protocol performance. We used four bone samples of different age in replicates to study the effects of both extraction methods. Our results confirm Loreille’s overall increased gain of DNA when enough tissue is available and Dabney’s improved efficiency for retrieving shorter DNA fragments that is beneficial when highly degraded DNA is present. The results suggest that the choice of extraction method needs to be based on available sample, degradation state, and targeted genotyping method. We modified the Dabney protocol by pooling parallel lysates prior to purification to study gain and performance in single tube typing assays and found that up to six parallel lysates lead to an almost linear gain of extracted DNA. These data are promising for further forensic investigations as the adapted Dabney protocol combines increased sensitivity for degraded DNA with necessary total DNA amount for forensic applications.

Ancient DNA from Guam and the peopling of the Pacific

Humans reached the Mariana Islands in the western Pacific by ∼3,500 y ago, contemporaneous with or even earlier than the initial peopling of Polynesia. They crossed more than 2,000 km of open ocean to get there, whereas voyages of similar length did not occur anywhere else until more than 2,000 y later. Yet, the settlement of Polynesia has received far more attention than the settlement of the Marianas. There is uncertainty over both the origin of the first colonizers of the Marianas (with different lines of evidence suggesting variously the Philippines, Indonesia, New Guinea, or the Bismarck Archipelago) as well as what, if any, relationship they might have had with the first colonizers of Polynesia. To address these questions, we obtained ancient DNA data from two skeletons from the Ritidian Beach Cave Site in northern Guam, dating to ∼2,200 y ago. Analyses of complete mitochondrial DNA genome sequences and genome-wide SNP data strongly support ancestry from the Philippines, in agreement with some interpretations of the linguistic and archaeological evidence, but in contradiction to results based on computer simulations of sea voyaging. We also find a close link between the ancient Guam skeletons and early Lapita individuals from Vanuatu and Tonga, suggesting that the Marianas and Polynesia were colonized from the same source population, and raising the possibility that the Marianas played a role in the eventual settlement of Polynesia.

The elusive parasite: comparing macroscopic, immunological, and genomic approaches to identifying malaria in human skeletal remains from Sayala, Egypt (third to sixth centuries AD)

Although malaria is one of the oldest and most widely distributed diseases affecting humans, identifying and characterizing its presence in ancient human remains continue to challenge researchers. We attempted to establish a reliable approach to detecting malaria in human skeletons using multiple avenues of analysis: macroscopic observations, rapid diagnostic tests, and shotgun-capture sequencing techniques, to identify pathological changes, Plasmodium antigens, and Plasmodium DNA, respectively. Bone and tooth samples from ten individuals who displayed skeletal lesions associated with anaemia, from a site in southern Egypt (third to sixth centuries AD), were selected. Plasmodium antigens were detected in five of the ten bone samples, and traces of Plasmodium aDNA were detected in six of the twenty bone and tooth samples. There was relatively good synchronicity between the biomolecular findings, despite not being able to authenticate the results. This study highlights the complexity and limitations in the conclusive identification of the Plasmodium parasite in ancient human skeletons. Limitations regarding antigen and aDNA preservation and the importance of sample selection are at the forefront of the search for malaria in the past. We confirm that, currently, palaeopathological changes such as cribra orbitalia are not enough to be certain of the presence of malaria. While biomolecular methods are likely the best chance for conclusive identification, we were unable to obtain results which correspond to the current authentication criteria of biomolecules. This study represents an important contribution in the refinement of biomolecular techniques used; also, it raises new insight regarding the consistency of combining several approaches in the identification of malaria in past populations.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12520-021-01350-z.

An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole-genome sequencing data

Ancient DNA research has developed rapidly over the past few decades due to improvements in PCR and next-generation sequencing (NGS) technologies, but challenges still exist. One major challenge in relation to ancient DNA research is to recover genuine endogenous ancient DNA sequences from raw sequencing data. This is often difficult due to degradation of ancient DNA and high levels of contamination, especially homologous contamination that has extremely similar genetic background with that of the real ancient DNA. In this study, we collected whole-genome sequencing (WGS) data from 6 ancient samples to compare different mapping algorithms. To further explore more effective methods to separate endogenous DNA from homologous contaminations, we attempted to recover reads based on ancient DNA specific characteristics of deamination, depurination, and DNA fragmentation with different parameters. We propose a quick and improved pipeline for separating endogenous ancient DNA while simultaneously decreasing homologous contaminations to very low proportions. Our goal in this research was to develop useful recommendations for ancient DNA mapping and for separation of endogenous DNA to facilitate future studies of ancient DNA.

Three Phases of Ancient Migration Shaped the Ancestry of Human Populations in Vanuatu

The archipelago of Vanuatu has been at the crossroads of human population movements in the Pacific for the past three millennia. To help address several open questions regarding the history of these movements, we generated genome-wide data for 11 ancient individuals from the island of Efate dating from its earliest settlement to the recent past, including five associated with the Chief Roi Mata's Domain World Heritage Area, and analyzed them in conjunction with 34 published ancient individuals from Vanuatu and elsewhere in Oceania, as well as present-day populations. Our results outline three distinct periods of population transformations. First, the four earliest individuals, from the Lapita-period site of Teouma, are concordant with eight previously described Lapita-associated individuals from Vanuatu and Tonga in having almost all of their ancestry from a "First Remote Oceanian" source related to East and Southeast Asians. Second, both the Papuan ancestry predominating in Vanuatu for the past 2,500 years and the smaller component of Papuan ancestry found in Polynesians can be modeled as deriving from a single source most likely originating in New Britain, suggesting that the movement of people carrying this ancestry to Remote Oceania closely followed that of the First Remote Oceanians in time and space. Third, the Chief Roi Mata's Domain individuals descend from a mixture of Vanuatu- and Polynesian-derived ancestry and are related to Polynesian-influenced communities today in central, but not southern, Vanuatu, demonstrating Polynesian genetic input in multiple groups with independent histories.

The evolutionary history of Neanderthal and Denisovan Y chromosomes

Ancient DNA has provided new insights into many aspects of human history. However, we lack comprehensive studies of the Y chromosomes of Denisovans and Neanderthals because the majority of specimens that have been sequenced to sufficient coverage are female. Sequencing Y chromosomes from two Denisovans and three Neanderthals shows that the Y chromosomes of Denisovans split around 700 thousand years ago from a lineage shared by Neanderthals and modern human Y chromosomes, which diverged from each other around 370 thousand years ago. The phylogenetic relationships of archaic and modern human Y chromosomes differ from the population relationships inferred from the autosomal genomes and mirror mitochondrial DNA phylogenies, indicating replacement of both the mitochondrial and Y chromosomal gene pools in late Neanderthals. This replacement is plausible if the low effective population size of Neanderthals resulted in an increased genetic load in Neanderthals relative to modern humans.

Human DNA extraction from highly degraded skeletal remains: How to find a suitable method?

Retrieving DNA from highly degraded human skeletal remains is still a challenge due to low concentration and fragmentation, which makes it difficult to extract and purify. Recent works showed that silica-based methods allow better DNA recovery and this fact may be attributed to the type of bones and the quality of the preserved tissue. However, more systematic studies are needed to evaluate the efficiency of the different silica-based extraction methods considering the type of bones. The main goal of the present study is to establish the best extraction method and the type of bone that can maximize the recovery of PCR-amplifiable DNA and the subsequent retrieval of mitochondrial and nuclear genetic information. Five individuals were selected from an archaeological site located in Catalonia-Spain dating from 5th to 11th centuries AD. For each individual, five samples from different skeletal regions were collected: petrous bone, pulp cavity and cementum of tooth, and rib and limb bones. Four extraction methods were tested, three silica-based (silica in-suspension, HE column and XS plasma column) and the classical method based on phenol-chloroform. Silica in-suspension method from petrous bone and pulp cavity showed the best results. However, the remains preservation will ultimately be the key to the molecular result success.

New perspectives on Neanderthal dispersal and turnover from Stajnia Cave (Poland)

The Micoquian is the broadest and longest enduring cultural facies of the Late Middle Palaeolithic that spread across the periglacial and boreal environments of Europe between Eastern France, Poland, and Northern Caucasus. Here, we present new data from the archaeological record of Stajnia Cave (Poland) and the paleogenetic analysis of a Neanderthal molar S5000, found in a Micoquian context. Our results demonstrate that the mtDNA genome of Stajnia S5000 dates to MIS 5a making the tooth the oldest Neanderthal specimen from Central-Eastern Europe. Furthermore, S5000 mtDNA has the fewest number of differences to mtDNA of Mezmaiskaya 1 Neanderthal from Northern Caucasus, and is more distant from almost contemporaneous Neanderthals of Scladina and Hohlenstein-Stadel. This observation and the technological affinity between Poland and the Northern Caucasus could be the result of increased mobility of Neanderthals that changed their subsistence strategy for coping with the new low biomass environments and the increased foraging radius of gregarious animals. The Prut and Dniester rivers were probably used as the main corridors of dispersal. The persistence of the Micoquian techno-complex in South-Eastern Europe infers that this axis of mobility was also used at the beginning of MIS 3 when a Neanderthal population turnover occurred in the Northern Caucasus.

Genetic Diversity of Historical and Modern Populations of Russian Cattle Breeds Revealed by Microsatellite Analysis

Analysis of ancient and historical DNA has great potential to trace the genetic diversity of local cattle populations during their centuries-long development. Forty-nine specimens representing five cattle breeds (Kholmogor, Yaroslavl, Great Russian, Novgorod, and Holland), dated from the end of the 19th century to the first half of the 20th century, were genotyped for nine polymorphic microsatellite loci. Using a multiple-tube approach, we determined the consensus genotypes of all samples/loci analysed. Amplification errors, including allelic drop-out (ADO) and false alleles (FA), occurred with an average frequency of 2.35% and 0.79%, respectively. A significant effect of allelic length on ADO rate (r² = 0.620, p = 0.05) was shown. We did not observe significant differences in genetic diversity among historical samples and modern representatives of Kholmogor and Yaroslavl breeds. The unbiased expected heterozygosity values were 0.726–0.774 and 0.708–0.739; the allelic richness values were 2.716–2.893 and 2.661–2.758 for the historical and modern samples, respectively. Analyses of F_ST and Jost’s D genetic distances, and the results of STRUCTURE clustering, showed the maintenance of a part of historical components in the modern populations of Kholmogor and Yaroslavl cattle. Our study contributes to the conservation of biodiversity in the local Russian genetic resources of cattle.

Ancient genomes in South Patagonia reveal population movements associated with technological shifts and geography

Archaeological research documents major technological shifts among people who have lived in the southern tip of South America (South Patagonia) during the last thirteen millennia, including the development of marine-based economies and changes in tools and raw materials. It has been proposed that movements of people spreading culture and technology propelled some of these shifts, but these hypotheses have not been tested with ancient DNA. Here we report genome-wide data from 20 ancient individuals, and co-analyze it with previously reported data. We reveal that immigration does not explain the appearance of marine adaptations in South Patagonia. We describe partial genetic continuity since ~6600 BP and two later gene flows correlated with technological changes: one between 4700-2000 BP that affected primarily marine-based groups, and a later one impacting all <2000 BP groups. From ~2200-1200 BP, mixture among neighbors resulted in a cline correlated to geographic ordering along the coast.

Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA

It has been shown that highly fragmented DNA is most efficiently converted into DNA libraries for sequencing if both strands of the DNA fragments are processed independently. We present an updated protocol for library preparation from single-stranded DNA, which is based on the splinted ligation of an adapter oligonucleotide to the 3' ends of single DNA strands, the synthesis of a complementary strand using a DNA polymerase and the addition of a 5' adapter via blunt-end ligation. The efficiency of library preparation is determined individually for each sample using a spike-in oligonucleotide. The whole workflow, including library preparation, quantification and amplification, requires two work days for up to 16 libraries. Alternatively, we provide documentation and electronic protocols enabling automated library preparation of 96 samples in parallel on a Bravo NGS Workstation (Agilent Technologies). After library preparation, molecules with uninformative short inserts (shorter than ~30-35 base pairs) can be removed by polyacrylamide gel electrophoresis if desired.

Environmental DNA and RNA as Records of Human Exposome, Including Biotic/Abiotic Exposures and Its Implications in the Assessment of the Role of Environment in Chronic Diseases

Most of environment-related diseases often result from multiple exposures of abiotic and/or biotic stressors across various life stages. The application of environmental DNA/RNA (eDNA/eRNA) to advance ecological understanding has been very successfully used. However, the eminent extension of eDNA/eRNA-based approaches to estimate human exposure to biotic and/or abiotic environmental stressors to understand the environmental causes of chronic diseases has yet to start. Here, we introduce the potential of eDNA/eRNA for bio-monitoring of human exposome and health effects in the real environmental or occupational settings. This review is the first of its kind to discuss how eDNA/eRNA-based approaches can be applied for assessing the human exposome. eDNA-based exposome assessment is expected to rely on our ability to capture the genome- and epigenome-wide signatures left behind by individuals in the indoor and outdoor physical spaces through shedding, excreting, etc. Records of eDNA/eRNA exposome may reflect the early appearance, persistence, and presence of biotic and/or abiotic-exposure-mediated modifications in these nucleic acid molecules. Functional genome- and epigenome-wide mapping of eDNA offer great promise to help elucidate the human exposome. Assessment of longitudinal exposure to physical, biological, and chemical agents present in the environment through eDNA/eRNA may enable the building of an integrative causal dynamic stochastic model to estimate environmental causes of human health deficits. This model is expected to incorporate key biological pathways and gene networks linking individuals, their geographic locations, and random multi-hits of environmental factors. Development and validation of monitoring of eDNA/eRNA exposome should seriously be considered to introduce into safety and risk assessment and as surrogates of chronic exposure to environmental stressors. Here we highlight that eDNA/eRNA reflecting longitudinal exposure of both biotic and abiotic environmental stressors may serve as records of human exposome and discuss its application as molecular tools for understanding the toxicogenomics basis of environment-related health deficits.

Present-Day DNA Contamination in Ancient DNA Datasets

Present-day contamination can lead to false conclusions in ancient DNA studies. A number of methods are available to estimate contamination, which use a variety of signals and are appropriate for different types of data. Here an overview of currently available methods highlighting their strengths and weaknesses is provided, and a classification based on the signals used to estimate contamination is proposed. This overview aims at enabling researchers to choose the most appropriate methods for their dataset. Based on this classification, potential avenues for the further development of methods are discussed.

Advantages of an easy-to-use DNA extraction method for minimal-destructive analysis of collection specimens

Here we present and justify an approach for minimal-destructive DNA extraction from historic insect specimens for next generation sequencing applications. An increasing number of studies use insects from museum collections for biodiversity research. However, the availability of specimens for molecular analyses has been limited by the degraded nature of the DNA gained from century-old museum material and the consumptive nature of most DNA extraction procedures. The method described in this manuscript enabled us to successfully extract DNA from specimens as old as 241 years using a minimal-destructive approach. The direct comparison of the DNeasy extraction Kit and the Monarch® PCR & DNA Clean-up Kit showed a significant increase of 17.3-fold higher DNA yield extracted with the Monarch Oligo protocol on average. By using an extraction protocol originally designed for oligonucleotide clean-up, we were able to combine overcoming the restrictions by target fragment size and strand state, with minimising time consumption and labour-intensity. The type specimens used for the minimal-destructive DNA extraction exhibited no significant external change or post-extraction damage, while sufficient DNA was retrieved for analyses.

Latest trends in archaeogenetic research of west Eurasians

During the past ten years, archaeogenetic research has exponentially grown to study the genetic history of human populations, using genome-wide data from large numbers of ancient individuals. Of the entire globe, Europe and the Near East are the regions where ancient DNA data is by far most abundant with over 2500 genomes published at present. In this review, we focus on archaeological contexts that have received less attention in the literature, specifically those associated with west Eurasian hunter-gatherers as well as populations from the Iron Age and later historical periods. In addition, we emphasize a recent shift from continent-wide to regional and even site-specific studies, which is starting to provide novel insights into sociocultural aspects of past societies.

An optimized method for the extraction of ancient eukaryote DNA from marine sediments

Marine sedimentary ancient DNA (sedaDNA) provides a powerful means to reconstruct marine palaeo-communities across the food web. However, currently there are few optimized sedaDNA extraction protocols available to maximize the yield of small DNA fragments typical of ancient DNA (aDNA) across a broad diversity of eukaryotes. We compared seven combinations of sedaDNA extraction treatments and sequencing library preparations using marine sediments collected at a water depth of 104 m off Maria Island, Tasmania, in 2018. These seven methods contrasted frozen versus refrigerated sediment, bead-beating induced cell lysis versus ethylenediaminetetraacetic acid (EDTA) incubation, DNA binding in silica spin columns versus in silica-solution, diluted versus undiluted DNA in shotgun library preparations to test potential inhibition issues during amplification steps, and size-selection of low molecular-weight (LMW) DNA to increase the extraction efficiency of sedaDNA. Maximum efficiency was obtained from frozen sediments subjected to a combination of EDTA incubation and bead-beating, DNA binding in silica-solution, and undiluted DNA in shotgun libraries, across 45 marine eukaryotic taxa. We present an optimized extraction protocol integrating these steps, with an optional post-library LMW size-selection step to retain DNA fragments of ≤500 base pairs. We also describe a stringent bioinformatic filtering approach for metagenomic data and provide a comprehensive list of contaminants as a reference for future sedaDNA studies. The new extraction and data-processing protocol should improve quantitative paleo-monitoring of eukaryotes from marine sediments, as well as other studies relying on the detection of highly fragmented and degraded eukaryote DNA in sediments.

Ancient Plant Genomics in Archaeology, Herbaria, and the Environment

The ancient DNA revolution of the past 35 years has driven an explosion in the breadth, nuance, and diversity of questions that are approachable using ancient biomolecules, and plant research has been a constant, indispensable facet of these developments. Using archaeological, paleontological, and herbarium plant tissues, researchers have probed plant domestication and dispersal, plant evolution and ecology, paleoenvironmental composition and dynamics, and other topics across related disciplines. Here, we review the development of the ancient DNA discipline and the role of plant research in its progress and refinement. We summarize our understanding of long-term plant DNA preservation and the characteristics of degraded DNA. In addition, we discuss challenges in ancient DNA recovery and analysis and the laboratory and bioinformatic strategies used to mitigate them. Finally, we review recent applications of ancient plant genomic research.

Large-scale identification of human bone remains via SNP microarray analysis with reference SNP database

Single nucleotide polymorphisms (SNPs) are valuable markers complementary to conventional forensic short tandem repeat (STR) markers in genetic typing, with potential advantages in challenging forensic casework. With the advent of high-throughput technologies, such as microarrays and massively parallel sequencing, the use of SNP typing has now expanded to large-scale forensic applications. Herein, a forensic case is presented to demonstrate the usefulness of SNP typing in identifying large-scale human bone remains with reference database construction. A total of 402 bone remains were recovered from an island in the Jeju Province of Korea where a massive disaster occurred in 1948. The first phase of the identification process was accomplished via conventional DNA typing methods including autosomal and Y-chromosomal STR typing, and mitochondrial DNA sequencing, which resulted in the identification of 74 of 402 remains. The second phase of the identification involved the remaining 327 unidentified remains using SNP typing as a supplementary tool based on Affymetrix resequencing array. The SNP markers of 782 family members were also analyzed and a reference database was constructed for comparison. An additional 51 bone remains were identified in the second phase. SNP data obtained from the supplementary genotyping yielded additional genetic information as well as contributed to kinship testing to determine the second degrees of relationship. In addition SNPs are useful in discriminating ambiguous relationship when only STR data are available. A software program developed for SNP typing system enabled efficient kinship analysis for large-scale forensic identification. The results and the casework are described and discussed.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

The spread of steppe and Iranian-related ancestry in the islands of the western Mediterranean

Steppe-pastoralist-related ancestry reached Central Europe by at least 2500 BC, whereas Iranian farmer-related ancestry was present in Aegean Europe by at least 1900 BC. However, the spread of these ancestries into the western Mediterranean, where they have contributed to many populations that live today, remains poorly understood. Here, we generated genome-wide ancient-DNA data from the Balearic Islands, Sicily and Sardinia, increasing the number of individuals with reported data from 5 to 66. The oldest individual from the Balearic Islands (~2400 BC) carried ancestry from steppe pastoralists that probably derived from west-to-east migration from Iberia, although two later Balearic individuals had less ancestry from steppe pastoralists. In Sicily, steppe pastoralist ancestry arrived by ~2200 BC, in part from Iberia; Iranian-related ancestry arrived by the mid-second millennium BC, contemporary to its previously documented spread to the Aegean; and there was large-scale population replacement after the Bronze Age. In Sardinia, nearly all ancestry derived from the island's early farmers until the first millennium BC, with the exception of an outlier from the third millennium BC, who had primarily North African ancestry and who-along with an approximately contemporary Iberian-documents widespread Africa-to-Europe gene flow in the Chalcolithic. Major immigration into Sardinia began in the first millennium BC and, at present, no more than 56-62% of Sardinian ancestry is from its first farmers. This value is lower than previous estimates, highlighting that Sardinia, similar to every other region in Europe, has been a stage for major movement and mixtures of people.

The formation of human populations in South and Central Asia

By sequencing 523 ancient humans, we show that the primary source of ancestry in modern South Asians is a prehistoric genetic gradient between people related to early hunter-gatherers of Iran and Southeast Asia. After the Indus Valley Civilization's decline, its people mixed with individuals in the southeast to form one of the two main ancestral populations of South Asia, whose direct descendants live in southern India. Simultaneously, they mixed with descendants of Steppe pastoralists who, starting around 4000 years ago, spread via Central Asia to form the other main ancestral population. The Steppe ancestry in South Asia has the same profile as that in Bronze Age Eastern Europe, tracking a movement of people that affected both regions and that likely spread the distinctive features shared between Indo-Iranian and Balto-Slavic languages.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

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.

More than skin and bones: Comparing extraction methods and alternative sources of DNA from avian museum specimens

Next-generation sequencing has greatly expanded the utility and value of museum collections by revealing specimens as genomic resources. As the field of museum genomics grows, so does the need for extraction methods that maximize DNA yields. For avian museum specimens, the established method of extracting DNA from toe pads works well for most specimens. However, for some specimens, especially those of birds that are very small or very large, toe pads can be a poor source of DNA. In this study, we apply two DNA extraction methods (phenol-chloroform and silica column) to three different sources of DNA (toe pad, skin punch and bone) from 10 historical avian museum specimens. We show that a modified phenol-chloroform protocol yielded significantly more DNA than a silica column protocol (e.g., Qiagen DNeasy Blood & Tissue Kit) across all tissue types. However, extractions using the silica column protocol contained longer fragments on average than those using the phenol-chloroform protocol, probably as a result of loss of small fragments through the silica column. While toe pads yielded more DNA than skin punches and bone fragments, skin punches proved to be a reliable alternative source of DNA and might be especially appealing when toe pad extractions are impractical. Overall, we found that historical bird museum specimens contain substantial amounts of DNA for genomic studies under most extraction scenarios, but that a phenol-chloroform protocol consistently provides the high quantities of DNA required for most current genomic protocols.

The genomic history of the Iberian Peninsula over the past 8000 years

We assembled genome-wide data from 271 ancient Iberians, of whom 176 are from the largely unsampled period after 2000 BCE, thereby providing a high-resolution time transect of the Iberian Peninsula. We document high genetic substructure between northwestern and southeastern hunter-gatherers before the spread of farming. We reveal sporadic contacts between Iberia and North Africa by ~2500 BCE and, by ~2000 BCE, the replacement of 40% of Iberia's ancestry and nearly 100% of its Y-chromosomes by people with Steppe ancestry. We show that, in the Iron Age, Steppe ancestry had spread not only into Indo-European-speaking regions but also into non-Indo-European-speaking ones, and we reveal that present-day Basques are best described as a typical Iron Age population without the admixture events that later affected the rest of Iberia. Additionally, we document how, beginning at least in the Roman period, the ancestry of the peninsula was transformed by gene flow from North Africa and the eastern Mediterranean.

Ancient DNA reveals a multistep spread of the first herders into sub-Saharan Africa

How food production first entered eastern Africa ~5000 years ago and the extent to which people moved with livestock is unclear. We present genome-wide data from 41 individuals associated with Later Stone Age, Pastoral Neolithic (PN), and Iron Age contexts in what are now Kenya and Tanzania to examine the genetic impacts of the spreads of herding and farming. Our results support a multiphase model in which admixture between northeastern African-related peoples and eastern African foragers formed multiple pastoralist groups, including a genetically homogeneous PN cluster. Additional admixture with northeastern and western African-related groups occurred by the Iron Age. These findings support several movements of food producers while rejecting models of minimal admixture with foragers and of genetic differentiation between makers of distinct PN artifacts.