Removal of deaminated cytosines and detection of in vivo methylation in ancient DNA

Studying human and nonhuman primate evolutionary biology with powerful in vitro and in vivo functional genomics tools

In recent years, tools for functional genomic studies have become increasingly feasible for use by evolutionary anthropologists. In this review, we provide brief overviews of several exciting in vitro techniques that can be paired with "-omics" approaches (e.g., genomics, epigenomics, transcriptomics, proteomics, and metabolomics) for potentially powerful evolutionary insights. These in vitro techniques include ancestral protein resurrection, cell line experiments using primary, immortalized, and induced pluripotent stem cells, and CRISPR-Cas9 genetic manipulation. We also discuss how several of these methods can be used in vivo, for transgenic organism studies of human and nonhuman primate evolution. Throughout this review, we highlight example studies in which these approaches have already been used to inform our understanding of the evolutionary biology of modern and archaic humans and other primates while simultaneously identifying future opportunities for anthropologists to use this toolkit to help answer additional outstanding questions in evolutionary anthropology.

Targeted analysis of polymorphic loci from low-coverage shotgun sequence data allows accurate genotyping of HLA genes in historical human populations

The highly polymorphic human leukocyte antigen (HLA) plays a crucial role in adaptive immunity and is associated with various complex diseases. Accurate analysis of HLA genes using ancient DNA (aDNA) data is crucial for understanding their role in human adaptation to pathogens. Here, we describe the TARGT pipeline for targeted analysis of polymorphic loci from low-coverage shotgun sequence data. The pipeline was successfully applied to medieval aDNA samples and validated using both simulated aDNA and modern empirical sequence data from the 1000 Genomes Project. Thus the TARGT pipeline enables accurate analysis of HLA polymorphisms in historical (and modern) human populations.

Animal domestication in the era of ancient genomics

The domestication of animals led to a major shift in human subsistence patterns, from a hunter-gatherer to a sedentary agricultural lifestyle, which ultimately resulted in the development of complex societies. Over the past 15,000 years, the phenotype and genotype of multiple animal species, such as dogs, pigs, sheep, goats, cattle and horses, have been substantially altered during their adaptation to the human niche. Recent methodological innovations, such as improved ancient DNA extraction methods and next-generation sequencing, have enabled the sequencing of whole ancient genomes. These genomes have helped reconstruct the process by which animals entered into domestic relationships with humans and were subjected to novel selection pressures. Here, we discuss and update key concepts in animal domestication in light of recent contributions from ancient genomics.

Dynamic Signatures of the Epigenome: Friend or Foe?

Highly dynamic epigenetic signaling is influenced mainly by (micro)environmental stimuli and genetic factors. The exact mechanisms affecting particular epigenomic patterns differ dependently on the context. In the current review, we focus on the causes and effects of the dynamic signatures of the human epigenome as evaluated with the high-throughput profiling data and single-gene approaches. We will discuss three different aspects of phenotypic outcomes occurring as a consequence of epigenetics interplaying with genotype and environment. The first issue is related to the cases of environmental impacts on epigenetic profile, and its adverse and advantageous effects related to human health and evolutionary adaptation. The next topic will present a model of the interwoven co-evolution of genetic and epigenetic patterns exemplified with transposable elements (TEs) and their epigenetic repressors Krüppel-associated box zinc finger proteins (KRAB–ZNFs). The third aspect concentrates on the mitosis-based microevolution that takes place during carcinogenesis, leading to clonal diversity and expansion of tumor cells. The whole picture of epigenome plasticity and its role in distinct biological processes is still incomplete. However, accumulating data define epigenomic dynamics as an essential co-factor driving adaptation at the cellular and inter-species levels with a benefit or disadvantage to the host.

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.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

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.

Emergence of human-adapted Salmonella enterica is linked to the Neolithization process

It has been hypothesized that the Neolithic transition towards an agricultural and pastoralist economy facilitated the emergence of human adapted pathogens. Here, we recovered eight Salmonella enterica subsp. enterica genomes from human skeletons of transitional foragers, pastoralists, and agro-pastoralists in western Eurasia that were up to 6,500 years old. Despite the high genetic diversity of S. enterica all ancient bacterial genomes clustered in a single previously uncharacterized branch that contains S. enterica adapted to multiple mammalian species. All ancient bacterial genomes from prehistoric (agro-)pastoralists fall within a part of this branch that also includes the human-specific S. enterica Paratyphi C, illustrating the evolution of a human pathogen over a period of five thousand years. Bacterial genomic comparisons suggest that the earlier ancient strains were not host specific, differed in pathogenic potential, and experienced convergent pseudogenization that accompanied their downstream host adaptation. These observations support the concept that the emergence of human adapted S. enterica is linked to human cultural transformations.

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

More Arrows in the Ancient DNA Quiver: Use of Paleoepigenomes and Paleomicrobiomes to Investigate Animal Adaptation to Environment

Whether and how epigenetic mechanisms and the microbiome play a role in mammalian adaptation raised considerable attention and controversy, mainly because they have the potential to add new insights into the Modern Synthesis. Recent attempts to reconcile neo-Darwinism and neo-Lamarckism in a unified theory of molecular evolution give epigenetic mechanisms and microbiome a prominent role. However, supporting empirical data are still largely missing. Because experimental studies using extant animals can hardly be done over evolutionary timescales, we propose that advances in ancient DNA techniques provide a valid alternative. In this piece, we evaluate 1) the possible roles of epigenomes and microbiomes in animal adaptation, 2) advances in the retrieval of paleoepigenome and paleomicrobiome data using ancient DNA techniques, and 3) the plasticity of either and interactions between the epigenome and the microbiome, while emphasizing that it is essential to take both into account, as well as the underlying genetic factors that may confound the findings. We propose that advanced ancient DNA techniques should be applied to a wide range of past animals, so novel dynamics in animal evolution and adaption can be revealed.

Ancient West African foragers in the context of African population history

Our knowledge of ancient human population structure in sub-Saharan Africa, particularly prior to the advent of food production, remains limited. Here we report genome-wide DNA data from four children-two of whom were buried approximately 8,000 years ago and two 3,000 years ago-from Shum Laka (Cameroon), one of the earliest known archaeological sites within the probable homeland of the Bantu language group1-11. One individual carried the deeply divergent Y chromosome haplogroup A00, which today is found almost exclusively in the same region12,13. However, the genome-wide ancestry profiles of all four individuals are most similar to those of present-day hunter-gatherers from western Central Africa, which implies that populations in western Cameroon today-as well as speakers of Bantu languages from across the continent-are not descended substantially from the population represented by these four people. We infer an Africa-wide phylogeny that features widespread admixture and three prominent radiations, including one that gave rise to at least four major lineages deep in the history of modern humans.

Distinguishing extant elephants ivory from mammoth ivory using a short sequence of cytochrome b gene

Trade in ivory from extant elephant species namely Asian elephant (Elephas maximus), African savanna elephant (Loxodonta africana) and African forest elephant (Loxodonta cyclotis) is regulated internationally, while the trade in ivory from extinct species of Elephantidae, including woolly mammoth, is unregulated. This distinction creates opportunity for laundering and trading elephant ivory as mammoth ivory. The existing morphological and molecular genetics methods do not reliably distinguish the source of ivory items that lack clear identification characteristics or for which the quality of extracted DNA cannot support amplification of large gene fragments. We present a PCR-sequencing method based on 116 bp target sequence of the cytochrome b gene to specifically amplify elephantid DNA while simultaneously excluding non-elephantid species and ivory substitutes, and while avoiding contamination by human DNA. The partial Cytochrome b gene sequence enabled accurate association of ivory samples with their species of origin for all three extant elephants and from mammoth. The detection limit of the PCR system was as low as 10 copy numbers of target DNA. The amplification and sequencing success reached 96.7% for woolly mammoth ivory and 100% for African savanna elephant and African forest elephant ivory. This is the first validated method for distinguishing elephant from mammoth ivory and it provides forensic support for investigation of ivory laundering cases.

snpAD: an ancient DNA genotype caller

Motivation

The study of ancient genomes can elucidate the evolutionary past. However, analyses are complicated by base-modifications in ancient DNA molecules that result in errors in DNA sequences. These errors are particularly common near the ends of sequences and pose a challenge for genotype calling.

Results

I describe an iterative method that estimates genotype frequencies and errors along sequences to allow for accurate genotype calling from ancient sequences. The implementation of this method, called snpAD, performs well on high-coverage ancient data, as shown by simulations and by subsampling the data of a high-coverage Neandertal genome. Although estimates for low-coverage genomes are less accurate, I am able to derive approximate estimates of heterozygosity from several low-coverage Neandertals. These estimates show that low heterozygosity, compared to modern humans, was common among Neandertals.

Availability and implementation

The C++ code of snpAD is freely available at http://bioinf.eva.mpg.de/snpAD/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Phylogeography of the second plague pandemic revealed through analysis of historical Yersinia pestis genomes

The second plague pandemic, caused by Yersinia pestis, devastated Europe and the nearby regions between the 14th and 18th centuries AD. Here we analyse human remains from ten European archaeological sites spanning this period and reconstruct 34 ancient Y. pestis genomes. Our data support an initial entry of the bacterium through eastern Europe, the absence of genetic diversity during the Black Death, and low within-outbreak diversity thereafter. Analysis of post-Black Death genomes shows the diversification of a Y. pestis lineage into multiple genetically distinct clades that may have given rise to more than one disease reservoir in, or close to, Europe. In addition, we show the loss of a genomic region that includes virulence-related genes in strains associated with late stages of the pandemic. The deletion was also identified in genomes connected with the first plague pandemic (541-750 AD), suggesting a comparable evolutionary trajectory of Y. pestis during both events.

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.

Different maternal lineages revealed by ancient mitochondrial genome of Camelus bactrianus from China

Domestic Bactrian camel (Camelus bactrianus) used to be one of the most important livestock species in Chinese history, as well as the major transport carrier on the ancient Silk Road. However, archeological studies on Chinese C. bactrianus are still limited, and molecular biology research on this species is mainly focused on modern specimens. In this study, we retrieved the complete mitochondrial genome from a C. bactrianus specimen, which was excavated from northwestern China and dated at 1290-1180 cal. years before present (yBP). Phylogenetic analyses using 18 mitochondrial genomes indicated that the C. bactrianus clade was divided into two maternal lineages. The majority of samples originating from Iran to Japan and Mongolia belong to subclade A1, while our sample together with two Mongolian individuals formed the much smaller subclade A2. Furthermore, the divergence time of these two maternal lineages was estimated as 165 Kya (95% credibility interval 117-222 Kya), this might indicate that several different evolutionary lineages were incorporated into the domestic gene pool during the initial domestication process. Bayesian skyline plot (BSP) analysis suggest a slow increase in female effective population size of C. bactrianus from 5000 years ago, which corresponds to the beginning of domestication of C. bactrianus. The present study also revealed that there were extensive exchanges of genetic information among C. bactrianus populations in regions along the Silk Road.

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny

Ancient DNA (aDNA) sequencing has enabled reconstruction of speciation, migration, and admixture events for extinct taxa¹. Outside the permafrost, however, irreversible aDNA post-mortem degradation² has so far limited aDNA recovery to the past ~0.5 million years (Ma)³. Contrarily, tandem mass spectrometry (MS) allowed sequencing ~1.5 million year (Ma) old collagen type I (COL1)⁴ and suggested the presence of protein residues in Cretaceous fossil remains⁵, although with limited phylogenetic use⁶. In the absence of molecular evidence, the speciation of several Early and Middle Pleistocene extinct species remain contentious. In this study, we address the phylogenetic relationships of the Eurasian Pleistocene Rhinocerotidae^7–9 using a ~1.77 Ma old dental enamel proteome of a Stephanorhinus specimen from the Dmanisi archaeological site in Georgia (South Caucasus)¹⁰. Molecular phylogenetic analyses place the Dmanisi Stephanorhinus as a sister group to the woolly (Coelodonta antiquitatis) and Merck’s rhinoceros (S. kirchbergensis) clade. We show that Coelodonta evolved from an early Stephanorhinus lineage and that the latter includes at least two distinct evolutionary lines. As such, the genus Stephanorhinus is currently paraphyletic and its systematic revision is therefore needed. We demonstrate that Early Pleistocene dental enamel proteome sequencing overcomes the limits of ancient collagen- and aDNA-based phylogenetic inference. It also provides additional information about the sex and taxonomic assignment of the specimens analysed. Dental enamel, the hardest tissue in vertebrates¹¹, is highly abundant in the fossil record. Our findings reveal that palaeoproteomic investigation of this material can push biomolecular investigation further back into the Early Pleistocene.

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 pathogen genomics as an emerging tool for infectious disease research

Over the past decade, a genomics revolution, made possible through the development of high-throughput sequencing, has triggered considerable progress in the study of ancient DNA, enabling complete genomes of past organisms to be reconstructed. A newly established branch of this field, ancient pathogen genomics, affords an in-depth view of microbial evolution by providing a molecular fossil record for a number of human-associated pathogens. Recent accomplishments include the confident identification of causative agents from past pandemics, the discovery of microbial lineages that are now extinct, the extrapolation of past emergence events on a chronological scale and the characterization of long-term evolutionary history of microorganisms that remain relevant to public health today. In this Review, we discuss methodological advancements, persistent challenges and novel revelations gained through the study of ancient pathogen genomes.

The origins and adaptation of European potatoes reconstructed from historical genomes

Potato, one of the most important staple crops, originates from the highlands of the equatorial Andes. There, potatoes propagate vegetatively via tubers under short days, constant throughout the year. After their introduction to Europe in the sixteenth century, potatoes adapted to a shorter growing season and to tuber formation under long days. Here, we traced the demographic and adaptive history of potato introduction to Europe. To this end, we sequenced 88 individuals that comprise landraces, modern cultivars and historical herbarium samples, including specimens collected by Darwin during the voyage of the Beagle. Our findings show that European potatoes collected during the period 1650-1750 were closely related to Andean landraces. After their introduction to Europe, potatoes admixed with Chilean genotypes. We identified candidate genes putatively involved in long-day pre-adaptation, and showed that the 1650-1750 European individuals were not long-day adapted through previously described allelic variants of the CYCLING DOF FACTOR1 gene. Such allelic variants were detected in Europe during the nineteenth century. Our study highlights the power of combining contemporary and historical genomes to understand the complex evolutionary history of crop adaptation to new environments.

Ancient Yersinia pestis genomes from across Western Europe reveal early diversification during the First Pandemic (541-750)

The first historically documented pandemic caused by Yersinia pestis began as the Justinianic Plague in 541 within the Roman Empire and continued as the so-called First Pandemic until 750. Although paleogenomic studies have previously identified the causative agent as Y. pestis, little is known about the bacterium's spread, diversity, and genetic history over the course of the pandemic. To elucidate the microevolution of the bacterium during this time period, we screened human remains from 21 sites in Austria, Britain, Germany, France, and Spain for Y. pestis DNA and reconstructed eight genomes. We present a methodological approach assessing single-nucleotide polymorphisms (SNPs) in ancient bacterial genomes, facilitating qualitative analyses of low coverage genomes from a metagenomic background. Phylogenetic analysis on the eight reconstructed genomes reveals the existence of previously undocumented Y. pestis diversity during the sixth to eighth centuries, and provides evidence for the presence of multiple distinct Y. pestis strains in Europe. We offer genetic evidence for the presence of the Justinianic Plague in the British Isles, previously only hypothesized from ambiguous documentary accounts, as well as the parallel occurrence of multiple derived strains in central and southern France, Spain, and southern Germany. Four of the reported strains form a polytomy similar to others seen across the Y. pestis phylogeny, associated with the Second and Third Pandemics. We identified a deletion of a 45-kb genomic region in the most recent First Pandemic strains affecting two virulence factors, intriguingly overlapping with a deletion found in 17th- to 18th-century genomes of the Second Pandemic.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

DamMet: ancient methylome mapping accounting for errors, true variants, and post-mortem DNA damage

BACKGROUND

Recent computational advances in ancient DNA research have opened access to the detection of ancient DNA methylation footprints at the genome-wide scale. The most commonly used approach infers the methylation state of a given genomic region on the basis of the amount of nucleotide mis-incorporations observed at CpG dinucleotide sites. However, this approach overlooks a number of confounding factors, including the presence of sequencing errors and true variants. The scale and distribution of the inferred methylation measurements are also variable across samples, precluding direct comparisons.

FINDINGS

Here, we present DamMet, an open-source software program retrieving maximum likelihood estimates of regional CpG methylation levels from ancient DNA sequencing data. It builds on a novel statistical model of post-mortem DNA damage for dinucleotides, accounting for sequencing errors, genotypes, and differential post-mortem cytosine deamination rates at both methylated and unmethylated sites. To validate DamMet, we extended gargammel, a sequence simulator for ancient DNA data, by introducing methylation-dependent features of post-mortem DNA decay. This new simulator provides direct validation of DamMet predictions. Additionally, the methylation levels inferred by DamMet were found to be correlated to those inferred by epiPALEOMIX and both on par and directly comparable to those measured from whole-genome bisulphite sequencing experiments of fresh tissues.

CONCLUSIONS

DamMet provides genuine estimates for local DNA methylation levels in ancient individual genomes. The returned estimates are directly cross-sample comparable, and the software is available as an open-source C++ program hosted at https://gitlab.com/KHanghoj/DamMet along with a manual and tutorial.

Assessing alignment-based taxonomic classification of ancient microbial DNA

The field of palaeomicrobiology-the study of ancient microorganisms-is rapidly growing due to recent methodological and technological advancements. It is now possible to obtain vast quantities of DNA data from ancient specimens in a high-throughput manner and use this information to investigate the dynamics and evolution of past microbial communities. However, we still know very little about how the characteristics of ancient DNA influence our ability to accurately assign microbial taxonomies (i.e. identify species) within ancient metagenomic samples. Here, we use both simulated and published metagenomic data sets to investigate how ancient DNA characteristics affect alignment-based taxonomic classification. We find that nucleotide-to-nucleotide, rather than nucleotide-to-protein, alignments are preferable when assigning taxonomies to short DNA fragment lengths routinely identified within ancient specimens (<60 bp). We determine that deamination (a form of ancient DNA damage) and random sequence substitutions corresponding to ∼100,000 years of genomic divergence minimally impact alignment-based classification. We also test four different reference databases and find that database choice can significantly bias the results of alignment-based taxonomic classification in ancient metagenomic studies. Finally, we perform a reanalysis of previously published ancient dental calculus data, increasing the number of microbial DNA sequences assigned taxonomically by an average of 64.2-fold and identifying microbial species previously unidentified in the original study. Overall, this study enhances our understanding of how ancient DNA characteristics influence alignment-based taxonomic classification of ancient microorganisms and provides recommendations for future palaeomicrobiological studies.

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.

Historical Genomes Reveal the Genomic Consequences of Recent Population Decline in Eastern Gorillas

Many endangered species have experienced severe population declines within the last centuries [1, 2]. However, despite concerns about negative fitness effects resulting from increased genetic drift and inbreeding, there is a lack of empirical data on genomic changes in conjunction with such declines [3-7]. Here, we use whole genomes recovered from century-old historical museum specimens to quantify the genomic consequences of small population size in the critically endangered Grauer's and endangered mountain gorillas. We find a reduction of genetic diversity and increase in inbreeding and genetic load in the Grauer's gorilla, which experienced severe population declines in recent decades. In contrast, the small but relatively stable mountain gorilla population has experienced little genomic change during the last century. These results suggest that species histories as well as the rate of demographic change may influence how population declines affect genome diversity.

Multiproxy evidence highlights a complex evolutionary legacy of maize in South America

Domesticated maize evolved from wild teosinte under human influences in Mexico beginning around 9000 years before the present (yr B.P.), traversed Central America by ~7500 yr B.P., and spread into South America by ~6500 yr B.P. Landrace and archaeological maize genomes from South America suggest that the ancestral population to South American maize was brought out of the domestication center in Mexico and became isolated from the wild teosinte gene pool before traits of domesticated maize were fixed. Deeply structured lineages then evolved within South America out of this partially domesticated progenitor population. Genomic, linguistic, archaeological, and paleoecological data suggest that the southwestern Amazon was a secondary improvement center for partially domesticated maize. Multiple waves of human-mediated dispersal are responsible for the diversity and biogeography of modern South American maize.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.

Reconstructing the Deep Population History of Central and South America

We report genome-wide ancient DNA from 49 individuals forming four parallel time transects in Belize, Brazil, the Central Andes, and the Southern Cone, each dating to at least ∼9,000 years ago. The common ancestral population radiated rapidly from just one of the two early branches that contributed to Native Americans today. We document two previously unappreciated streams of gene flow between North and South America. One affected the Central Andes by ∼4,200 years ago, while the other explains an affinity between the oldest North American genome associated with the Clovis culture and the oldest Central and South Americans from Chile, Brazil, and Belize. However, this was not the primary source for later South Americans, as the other ancient individuals derive from lineages without specific affinity to the Clovis-associated genome, suggesting a population replacement that began at least 9,000 years ago and was followed by substantial population continuity in multiple regions.

Archival May-Grünwald-Giemsa-Stained Bone Marrow Smears Are an Eligible Source for Molecular DNA Research

Biobanking is increasingly important in studying complex heterogeneous diseases. Therefore, it is essential to ensure the sample quality after long-term storage for reliable downstream analyses. The Clinical Biobank of the Jessa Hospital and the University Biobank Limburg (UBiLim) hold a continuously growing collection of hematological samples, including May-Grünwald-Giemsa (MGG)- and Perls' Prussian Blue (PPB)-stained bone marrow (BM) smears, stored at room temperature (RT) for up to 20 years. In this study, we investigated the effect of short- and long-term storage on the quality of DNA and RNA extracted from these BM smears to assess their fitness-for-purpose in downstream molecular applications, including agarose gel electrophoresis, bio-analyzer analysis, quantitative polymerase chain reaction (qPCR), and targeted next-generation sequencing (NGS). The RNA quality was very low for all samples, independent of storage time or staining method. The DNA from PPB-stained BM smears was already degraded after 1 year of storage and correspondingly could not be used for reliable downstream molecular analysis. In contrast, DNA extracted from MGG-stained BM smears stored for up to 10 years was able to generate high-quality data in qPCR and targeted NGS analyses. Longer storage periods (>15 years) of these samples revealed a high degree of degradation and a significant amount of DNA transitions and transversions. In conclusion, the DNA extracted from archival MGG-stained BM smears with a storage time up to at least 10 years was qualitatively good and fit for downstream analysis, including targeted NGS. This indicates that these samples are an eligible source for molecular DNA research and for studying complex diseases.

Early human dispersals within the Americas

Studies of the peopling of the Americas have focused on the timing and number of initial migrations. Less attention has been paid to the subsequent spread of people within the Americas. We sequenced 15 ancient human genomes spanning from Alaska to Patagonia; six are ≥10,000 years old (up to ~18× coverage). All are most closely related to Native Americans, including those from an Ancient Beringian individual and two morphologically distinct “Paleoamericans.” We found evidence of rapid dispersal and early diversification that included previously unknown groups as people moved south. This resulted in multiple independent, geographically uneven migrations, including one that provides clues of a Late Pleistocene Australasian genetic signal, as well as a later Mesoamerican-related expansion. These led to complex and dynamic population histories from North to South America.

Intragenus (Homo) variation in a chemokine receptor gene (CCR5)

Humans have a comparatively higher rate of more polymorphisms in regulatory regions of the primate CCR5 gene, an immune system gene with both general and specific functions. This has been interpreted as allowing flexibility and diversity of gene expression in response to varying disease loads. A broad expression repertoire is useful to humans-the only globally distributed primate-due to our unique adaptive pattern that increased pathogen exposure and disease loads (e.g., sedentism, subsistence practices). The main objective of the study was to determine if the previously observed human pattern of increased variation extended to other members of our genus, Homo. The data for this study are mined from the published genomes of extinct hominins (four Neandertals and two Denisovans), an ancient human (Ust'-Ishim), and modern humans (1000 Genomes). An average of 15 polymorphisms per individual were found in human populations (with a total of 262 polymorphisms). There were 94 polymorphisms identified across extinct Homo (an average of 13 per individual) with 41 previously observed in modern humans and 53 novel polymorphisms (32 in Denisova and 21 in Neandertal). Neither the frequency nor distribution of polymorphisms across gene regions exhibit significant differences within the genus Homo. Thus, humans are not unique with regards to the increased frequency of regulatory polymorphisms and the evolution of variation patterns across CCR5 gene appears to have originated within the genus. A broader evolutionary perspective on regulatory flexibility may be that it provided an advantage during the transition to confrontational foraging (and later hunting) that altered human-environment interaction as well as during migration to Eurasia and encounters with novel pathogens.

Ancient goat genomes reveal mosaic domestication in the Fertile Crescent

Current genetic data are equivocal as to whether goat domestication occurred multiple times or was a singular process. We generated genomic data from 83 ancient goats (51 with genome-wide coverage) from Paleolithic to Medieval contexts throughout the Near East. Our findings demonstrate that multiple divergent ancient wild goat sources were domesticated in a dispersed process that resulted in genetically and geographically distinct Neolithic goat populations, echoing contemporaneous human divergence across the region. These early goat populations contributed differently to modern goats in Asia, Africa, and Europe. We also detect early selection for pigmentation, stature, reproduction, milking, and response to dietary change, providing 8000-year-old evidence for human agency in molding genome variation within a partner species.

Application of a novel molecular method to age free-living wild Bechstein's bats

The age profile of populations fundamentally affects their conservation status. Yet, age is frequently difficult to assess in wild animals. Here, we assessed the use of DNA methylation of homologous genes to establish the age structure of a rare and elusive wild mammal: the Bechstein's bat (Myotis bechsteinii). We collected 62 wing punches from individuals whose ages were known as a result of a long-term banding study. DNA methylation was measured at seven CpG sites from three genes, which have previously shown age-associated changes in humans and laboratory mice. All CpG sites from the tested genes showed a significant relationship between DNA methylation and age, both individually and in combination (multiple linear regression R²  = 0.58, p < 0.001). Despite slight approximation around estimates, the approach is sufficiently precise to place animals into practically useful age cohorts. This method is of considerable practical benefit as it can reliably age individual bats. It is also much faster than traditional capture-mark-recapture techniques, with the potential to collect information on the age structure of an entire colony from a single sampling session to better inform conservation actions for Bechstein's bats. By identifying three genes where DNA methylation correlates with age across distantly related species, this study also suggests that the technique can potentially be applied across a wide range of mammals.

The prehistoric peopling of Southeast Asia

The human occupation history of Southeast Asia (SEA) remains heavily debated. Current evidence suggests that SEA was occupied by Hòabìnhian hunter-gatherers until ~4000 years ago, when farming economies developed and expanded, restricting foraging groups to remote habitats. Some argue that agricultural development was indigenous; others favor the “two-layer” hypothesis that posits a southward expansion of farmers giving rise to present-day Southeast Asian genetic diversity. By sequencing 26 ancient human genomes (25 from SEA, 1 Japanese Jōmon), we show that neither interpretation fits the complexity of Southeast Asian history: Both Hòabìnhian hunter-gatherers and East Asian farmers contributed to current Southeast Asian diversity, with further migrations affecting island SEA and Vietnam. Our results help resolve one of the long-standing controversies in Southeast Asian prehistory.

Analysis of 3800-year-old Yersinia pestis genomes suggests Bronze Age origin for bubonic plague

The origin of Yersinia pestis and the early stages of its evolution are fundamental subjects of investigation given its high virulence and mortality that resulted from past pandemics. Although the earliest evidence of Y. pestis infections in humans has been identified in Late Neolithic/Bronze Age Eurasia (LNBA 5000–3500y BP), these strains lack key genetic components required for flea adaptation, thus making their mode of transmission and disease presentation in humans unclear. Here, we reconstruct ancient Y. pestis genomes from individuals associated with the Late Bronze Age period (~3800 BP) in the Samara region of modern-day Russia. We show clear distinctions between our new strains and the LNBA lineage, and suggest that the full ability for flea-mediated transmission causing bubonic plague evolved more than 1000 years earlier than previously suggested. Finally, we propose that several Y. pestis lineages were established during the Bronze Age, some of which persist to the present day.

Yersinia pestis has caused infections (plague) in humans since the Early Bronze Age (5000 years ago). Here, Spyrou et al. reconstruct Y. pestis genomes from Late Bronze Age individuals, and find genomic evidence compatible with flea-mediated transmission causing bubonic plague.

Historic Treponema pallidum genomes from Colonial Mexico retrieved from archaeological remains

Treponema pallidum infections occur worldwide causing, among other diseases, syphilis and yaws. In particular sexually transmitted syphilis is regarded as a re-emerging infectious disease with millions of new infections annually. Here we present three historic T. pallidum genomes (two from T. pallidum ssp. pallidum and one from T. pallidum ssp. pertenue) that have been reconstructed from skeletons recovered from the Convent of Santa Isabel in Mexico City, operational between the 17th and 19th century. Our analyses indicate that different T. pallidum subspecies caused similar diagnostic presentations that are normally associated with syphilis in infants, and potential evidence of a congenital infection of T. pallidum ssp. pertenue, the causative agent of yaws. This first reconstruction of T. pallidum genomes from archaeological material opens the possibility of studying its evolutionary history at a resolution previously assumed to be out of reach.

Ancient genomes document multiple waves of migration in Southeast Asian prehistory

Southeast Asia is home to rich human genetic and linguistic diversity, but the details of past population movements in the region are not well known. Here, we report genome-wide ancient DNA data from 18 Southeast Asian individuals spanning from the Neolithic period through the Iron Age (4100 to 1700 years ago). Early farmers from Man Bac in Vietnam exhibit a mixture of East Asian (southern Chinese agriculturalist) and deeply diverged eastern Eurasian (hunter-gatherer) ancestry characteristic of Austroasiatic speakers, with similar ancestry as far south as Indonesia providing evidence for an expansive initial spread of Austroasiatic languages. By the Bronze Age, in a parallel pattern to Europe, sites in Vietnam and Myanmar show close connections to present-day majority groups, reflecting substantial additional influxes of migrants.

Recently evolved human-specific methylated regions are enriched in schizophrenia signals

Background

One explanation for the persistence of schizophrenia despite the reduced fertility of patients is that it is a by-product of recent human evolution. This hypothesis is supported by evidence suggesting that recently-evolved genomic regions in humans are involved in the genetic risk for schizophrenia. Using summary statistics from genome-wide association studies (GWAS) of schizophrenia and 11 other phenotypes, we tested for enrichment of association with GWAS traits in regions that have undergone methylation changes in the human lineage compared to Neanderthals and Denisovans, i.e. human-specific differentially methylated regions (DMRs). We used analytical tools that evaluate polygenic enrichment of a subset of genomic variants against all variants.

Results

Schizophrenia was the only trait in which DMR SNPs showed clear enrichment of association that passed the genome-wide significance threshold. The enrichment was not observed for Neanderthal or Denisovan DMRs. The enrichment seen in human DMRs is comparable to that for genomic regions tagged by Neanderthal Selective Sweep markers, and stronger than that for Human Accelerated Regions. The enrichment survives multiple testing performed through permutation (n = 10,000) and bootstrapping (n = 5000) in INRICH (p < 0.01). Some enrichment of association with height was observed at the gene level.

Conclusions

Regions where DNA methylation modifications have changed during recent human evolution show enrichment of association with schizophrenia and possibly with height. Our study further supports the hypothesis that genetic variants conferring risk of schizophrenia co-occur in genomic regions that have changed as the human species evolved. Since methylation is an epigenetic mark, potentially mediated by environmental changes, our results also suggest that interaction with the environment might have contributed to that association.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1177-2) contains supplementary material, which is available to authorized users.

Ancient genomes reveal a high diversity of Mycobacterium leprae in medieval Europe

Studying ancient DNA allows us to retrace the evolutionary history of human pathogens, such as Mycobacterium leprae, the main causative agent of leprosy. Leprosy is one of the oldest recorded and most stigmatizing diseases in human history. The disease was prevalent in Europe until the 16th century and is still endemic in many countries with over 200,000 new cases reported annually. Previous worldwide studies on modern and European medieval M. leprae genomes revealed that they cluster into several distinct branches of which two were present in medieval Northwestern Europe. In this study, we analyzed 10 new medieval M. leprae genomes including the so far oldest M. leprae genome from one of the earliest known cases of leprosy in the United Kingdom-a skeleton from the Great Chesterford cemetery with a calibrated age of 415-545 C.E. This dataset provides a genetic time transect of M. leprae diversity in Europe over the past 1500 years. We find M. leprae strains from four distinct branches to be present in the Early Medieval Period, and strains from three different branches were detected within a single cemetery from the High Medieval Period. Altogether these findings suggest a higher genetic diversity of M. leprae strains in medieval Europe at various time points than previously assumed. The resulting more complex picture of the past phylogeography of leprosy in Europe impacts current phylogeographical models of M. leprae dissemination. It suggests alternative models for the past spread of leprosy such as a wide spread prevalence of strains from different branches in Eurasia already in Antiquity or maybe even an origin in Western Eurasia. Furthermore, these results highlight how studying ancient M. leprae strains improves understanding the history of leprosy worldwide.

Ancient hepatitis B viruses from the Bronze Age to the Medieval period

Hepatitis B virus (HBV) is a major cause of human hepatitis. There is considerable uncertainty about the timescale of its evolution and its association with humans. Here we present 12 full or partial ancient HBV genomes that are between approximately 0.8 and 4.5 thousand years old. The ancient sequences group either within or in a sister relationship with extant human or other ape HBV clades. Generally, the genome properties follow those of modern HBV. The root of the HBV tree is projected to between 8.6 and 20.9 thousand years ago, and we estimate a substitution rate of 8.04 × 10^-6-1.51 × 10^-5 nucleotide substitutions per site per year. In several cases, the geographical locations of the ancient genotypes do not match present-day distributions. Genotypes that today are typical of Africa and Asia, and a subgenotype from India, are shown to have an early Eurasian presence. The geographical and temporal patterns that we observe in ancient and modern HBV genotypes are compatible with well-documented human migrations during the Bronze and Iron Ages^1,2. We provide evidence for the creation of HBV genotype A via recombination, and for a long-term association of modern HBV genotypes with humans, including the discovery of a human genotype that is now extinct. These data expose a complexity of HBV evolution that is not evident when considering modern sequences alone.

Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing

Ancient DNA research has been revolutionized following development of next-generation sequencing platforms. Although a number of such platforms have been applied to ancient DNA samples, the Illumina series are the dominant choice today, mainly because of high production capacities and short read production. Recently a potentially attractive alternative platform for palaeogenomic data generation has been developed, the BGISEQ-500, whose sequence output are comparable with the Illumina series. In this study, we modified the standard BGISEQ-500 library preparation specifically for use on degraded DNA, then directly compared the sequencing performance and data quality of the BGISEQ-500 to the Illumina HiSeq2500 platform on DNA extracted from 8 historic and ancient dog and wolf samples. The data generated were largely comparable between sequencing platforms, with no statistically significant difference observed for parameters including level (P = 0.371) and average sequence length (P = 0718) of endogenous nuclear DNA, sequence GC content (P = 0.311), double-stranded DNA damage rate (v. 0.309), and sequence clonality (P = 0.093). Small significant differences were found in single-strand DNA damage rate (δS; slightly lower for the BGISEQ-500, P = 0.011) and the background rate of difference from the reference genome (θ; slightly higher for BGISEQ-500, P = 0.012). This may result from the differences in amplification cycles used to polymerase chain reaction–amplify the libraries. A significant difference was also observed in the mitochondrial DNA percentages recovered (P = 0.018), although we believe this is likely a stochastic effect relating to the extremely low levels of mitochondria that were sequenced from 3 of the samples with overall very low levels of endogenous DNA. Although we acknowledge that our analyses were limited to animal material, our observations suggest that the BGISEQ-500 holds the potential to represent a valid and potentially valuable alternative platform for palaeogenomic data generation that is worthy of future exploration by those interested in the sequencing and analysis of degraded DNA.

Estimating genetic kin relationships in prehistoric populations

Archaeogenomic research has proven to be a valuable tool to trace migrations of historic and prehistoric individuals and groups, whereas relationships within a group or burial site have not been investigated to a large extent. Knowing the genetic kinship of historic and prehistoric individuals would give important insights into social structures of ancient and historic cultures. Most archaeogenetic research concerning kinship has been restricted to uniparental markers, while studies using genome-wide information were mainly focused on comparisons between populations. Applications which infer the degree of relationship based on modern-day DNA information typically require diploid genotype data. Low concentration of endogenous DNA, fragmentation and other post-mortem damage to ancient DNA (aDNA) makes the application of such tools unfeasible for most archaeological samples. To infer family relationships for degraded samples, we developed the software READ (Relationship Estimation from Ancient DNA). We show that our heuristic approach can successfully infer up to second degree relationships with as little as 0.1x shotgun coverage per genome for pairs of individuals. We uncover previously unknown relationships among prehistoric individuals by applying READ to published aDNA data from several human remains excavated from different cultural contexts. In particular, we find a group of five closely related males from the same Corded Ware culture site in modern-day Germany, suggesting patrilocality, which highlights the possibility to uncover social structures of ancient populations by applying READ to genome-wide aDNA data. READ is publicly available from https://bitbucket.org/tguenther/read.

The Beaker phenomenon and the genomic transformation of northwest Europe

From around 2750 to 2500 bc, Bell Beaker pottery became widespread across western and central Europe, before it disappeared between 2200 and 1800 bc. The forces that propelled its expansion are a matter of long-standing debate, and there is support for both cultural diffusion and migration having a role in this process. Here we present genome-wide data from 400 Neolithic, Copper Age and Bronze Age Europeans, including 226 individuals associated with Beaker-complex artefacts. We detected limited genetic affinity between Beaker-complex-associated individuals from Iberia and central Europe, and thus exclude migration as an important mechanism of spread between these two regions. However, migration had a key role in the further dissemination of the Beaker complex. We document this phenomenon most clearly in Britain, where the spread of the Beaker complex introduced high levels of steppe-related ancestry and was associated with the replacement of approximately 90% of Britain's gene pool within a few hundred years, continuing the east-to-west expansion that had brought steppe-related ancestry into central and northern Europe over the previous centuries.