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

Population Turnover in Remote Oceania Shortly after Initial Settlement

Ancient DNA from Vanuatu and Tonga dating to about 2,900-2,600 years ago (before present, BP) has revealed that the "First Remote Oceanians" associated with the Lapita archaeological culture were directly descended from the population that, beginning around 5000 BP, spread Austronesian languages from Taiwan to the Philippines, western Melanesia, and eventually Remote Oceania. Thus, ancestors of the First Remote Oceanians must have passed by the Papuan-ancestry populations they encountered in New Guinea, the Bismarck Archipelago, and the Solomon Islands with minimal admixture [1]. However, all present-day populations in Near and Remote Oceania harbor >25% Papuan ancestry, implying that additional eastward migration must have occurred. We generated genome-wide data for 14 ancient individuals from Efate and Epi Islands in Vanuatu from 2900-150 BP, as well as 185 present-day individuals from 18 islands. We find that people of almost entirely Papuan ancestry arrived in Vanuatu by around 2300 BP, most likely reflecting migrations a few hundred years earlier at the end of the Lapita period, when there is also evidence of changes in skeletal morphology and cessation of long-distance trade between Near and Remote Oceania [2, 3]. Papuan ancestry was subsequently diluted through admixture but remains at least 80%-90% in most islands. Through a fine-grained analysis of ancestry profiles, we show that the Papuan ancestry in Vanuatu derives from the Bismarck Archipelago rather than the geographically closer Solomon Islands. However, the Papuan ancestry in Polynesia-the most remote Pacific islands-derives from different sources, documenting a third stream of migration from Near to Remote Oceania.

Ancient genomes revisit the ancestry of domestic and Przewalski’s horses

The Eneolithic Botai culture of the Central Asian steppes provides the earliest archaeological evidence for horse husbandry, ~5500 years ago, but the exact nature of early horse domestication remains controversial. We generated 42 ancient-horse genomes, including 20 from Botai. Compared to 46 published ancient- and modern-horse genomes, our data indicate that Przewalski’s horses are the feral descendants of horses herded at Botai and not truly wild horses. All domestic horses dated from ~4000 years ago to present only show ~2.7% of Botai-related ancestry. This indicates that a massive genomic turnover underpins the expansion of the horse stock that gave rise to modern domesticates, which coincides with large-scale human population expansions during the Early Bronze Age.

The genomic history of southeastern Europe

Farming was first introduced to Europe in the mid-7^th millennium BCE–associated with migrants from Anatolia who settled in the Southeast before spreading throughout Europe. To understand the dynamics of this process, we analyzed genome-wide ancient DNA data from 225 individuals who lived in southeastern Europe and surrounding regions between 12,000 and 500 BCE. We document a West-East cline of ancestry in indigenous hunter-gatherers and–in far-eastern Europe–early stages in the formation of Bronze Age Steppe ancestry. We show that the first farmers of northern and western Europe passed through southeastern Europe with limited hunter-gatherer admixture, but that some groups that remained mixed extensively, without the male-biased hunter-gatherer admixture that prevailed later in the North and West. Southeastern Europe continued to be a nexus between East and West, with intermittent genetic contact with the Steppe up to 2000 years before the migrations that replaced much of northern Europe’s population.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

gDNA extraction yield and methylation status of blood samples are affected by long-term storage conditions

Epigenetics is believed to provide great chances for a better understanding of the development and treatment of many diseases where the analysis of genomic DNA has so far failed to provide conclusive answers. Methylcytosine is a frequently used quantitative marker of epigenetic studies. Since immediate analysis of sampled material is in most cases not possible, storage time and conditions are critical aspects regarding the quality of genomic DNA and reliability of analysis. Blood is frequently used for such analyses. We, therefore, collected blood samples of ten volunteers and stored them under various conditions for ten months: -70°C, -20°C, 2-8°C and room temperature. An additional aliquot was frozen at -70°C and thawed once a week at room temperature. We then compared the DNA extraction yields and methylation status in relation to storage time and conditions. We found significantly lower DNA extraction yields (up to -97.45%; p ≤ 0.001) as well as significantly higher methylation levels after ten months of storage (up to +42.0%; p ≤ 0.001). These results suggest that storage time has an important influence on DNA analyses of blood samples for all storage conditions. This might be due to differences in stability of methylated and non-methylated DNA. Our study indicates that storage conditions and time may be a critical factor for epigenetic methylation studies and require rigorous validation. For reliable analyses we, therefore, recommend to perform epigenetic analysis directly after sample collection.

Inferring genetic origins and phenotypic traits of George Bähr, the architect of the Dresden Frauenkirche

For historic individuals, the outward appearance and other phenotypic characteristics remain often non-resolved. Unfortunately, images or detailed written sources are only scarcely available in many cases. Attempts to study historic individuals with genetic data so far focused on hypervariable regions of mitochondrial DNA and to some extent on complete mitochondrial genomes. To elucidate the potential of in-solution based genome-wide SNP capture methods - as now widely applied in population genetics - we extracted DNA from the 17th century remains of George Bähr, the architect of the Dresdner Frauenkirche. We were able to identify the remains to be of male origin, showing sufficient DNA damage, deriving from a single person and being thus likely authentic. Furthermore, we were able to show that George Bähr had light skin pigmentation and most likely brown eyes. His genomic DNA furthermore points to a Central European origin. We see this analysis as an example to demonstrate the prospects that new in-solution SNP capture methods can provide for historic cases of forensic interest, using methods well established in ancient DNA (aDNA) research and population genetics.

The rate and potential relevance of new mutations in a colonizing plant lineage

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by "natural evolution experiments" that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions-the subset of mutations that survived natural selection and drift-, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

Salmonella enterica genomes from victims of a major sixteenth-century epidemic in Mexico

Indigenous populations of the Americas experienced high mortality rates during the early contact period as a result of infectious diseases, many of which were introduced by Europeans. Most of the pathogenic agents that caused these outbreaks remain unknown. Through the introduction of a new metagenomic analysis tool called MALT, applied here to search for traces of ancient pathogen DNA, we were able to identify Salmonella enterica in individuals buried in an early contact era epidemic cemetery at Teposcolula-Yucundaa, Oaxaca in southern Mexico. This cemetery is linked, based on historical and archaeological evidence, to the 1545-1550 CE epidemic that affected large parts of Mexico. Locally, this epidemic was known as 'cocoliztli', the pathogenic cause of which has been debated for more than a century. Here, we present genome-wide data from ten individuals for Salmonella enterica subsp. enterica serovar Paratyphi C, a bacterial cause of enteric fever. We propose that S. Paratyphi C be considered a strong candidate for the epidemic population decline during the 1545 cocoliztli outbreak at Teposcolula-Yucundaa.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500–6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east–west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

The Scandinavian peninsula was the last part of Europe to be colonized after the Last Glacial Maximum. The migration routes, cultural networks, and the genetic makeup of the first Scandinavians remain elusive and several hypotheses exist based on archaeology, climate modeling, and genetics. By analyzing the genomes of early Scandinavian hunter-gatherers, we show that their migrations followed two routes: one from the south and another from the northeast along the ice-free Norwegian Atlantic coast. These groups met and mixed in Scandinavia, creating a population more diverse than contemporaneous central and western European hunter-gatherers. As northern Europe is associated with cold and low light conditions, we investigated genomic patterns of adaptation to these conditions and genes known to be involved in skin pigmentation. We demonstrate that Mesolithic Scandinavians had higher levels of light pigmentation variants compared to the respective source populations of the migrations, suggesting adaptation to low light levels and a surprising signal of genetic continuity in TMEM131, a gene that may be involved in long-term adaptation to the cold.

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

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

Human evolution: a tale from ancient genomes

The field of human ancient DNA (aDNA) has moved from mitochondrial sequencing that suffered from contamination and provided limited biological insights, to become a fully genomic discipline that is changing our conception of human history. Recent successes include the sequencing of extinct hominins, and true population genomic studies of Bronze Age populations. Among the emerging areas of aDNA research, the analysis of past epigenomes is set to provide more new insights into human adaptation and disease susceptibility through time. Starting as a mere curiosity, ancient human genetics has become a major player in the understanding of our evolutionary history.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Mitogenome evolution in the last surviving woolly mammoth population reveals neutral and functional consequences of small population size

The onset of the Holocene was associated with a global temperature increase, which led to a rise in sea levels and isolation of the last surviving population of woolly mammoths on Wrangel Island. Understanding what happened with the population's genetic diversity at the time of the isolation and during the ensuing 6000 years can help clarify the effects of bottlenecks and subsequent limited population sizes in species approaching extinction. Previous genetic studies have highlighted questions about how the Holocene Wrangel population was established and how the isolation event affected genetic diversity. Here, we generated high‐quality mitogenomes from 21 radiocarbon‐dated woolly mammoths to compare the ancestral large and genetically diverse Late Pleistocene Siberian population and the small Holocene Wrangel population. Our results indicate that mitogenome diversity was reduced to one single haplotype at the time of the isolation, and thus that the Holocene Wrangel Island population was established by a single maternal lineage. Moreover, we show that the ensuing small effective population size coincided with fixation of a nonsynonymous mutation, and a comparative analysis of mutation rates suggests that the evolutionary rate was accelerated in the Holocene population. These results suggest that isolation on Wrangel Island led to an increase in the frequency of deleterious genetic variation, and thus are consistent with the hypothesis that strong genetic drift in small populations leads to purifying selection being less effective in removing deleterious mutations.

A high-coverage Neandertal genome from Vindija Cave in Croatia

To date, the only Neandertal genome that has been sequenced to high quality is from an individual found in Southern Siberia. We sequenced the genome of a female Neandertal from ~50,000 years ago from Vindija Cave, Croatia, to ~30-fold genomic coverage. She carried 1.6 differences per 10,000 base pairs between the two copies of her genome, fewer than present-day humans, suggesting that Neandertal populations were of small size. Our analyses indicate that she was more closely related to the Neandertals that mixed with the ancestors of present-day humans living outside of sub-Saharan Africa than the previously sequenced Neandertal from Siberia, allowing 10 to 20% more Neandertal DNA to be identified in present-day humans, including variants involved in low-density lipoprotein cholesterol concentrations, schizophrenia, and other diseases.

Inferring Past Environments from Ancient Epigenomes

Analyzing the conditions in which past individuals lived is key to understanding the environments and cultural transitions to which humans had to adapt. Here, we suggest a methodology to probe into past environments, using reconstructed premortem DNA methylation maps of ancient individuals. We review a large body of research showing that differential DNA methylation is associated with changes in various external and internal factors, and propose that loci whose DNA methylation level is environmentally responsive could serve as markers to infer about ancient daily life, diseases, nutrition, exposure to toxins, and more. We demonstrate this approach by showing that hunger-related DNA methylation changes are found in ancient hunter-gatherers. The strategy we present here opens a window to reconstruct previously inaccessible aspects of the lives of past individuals.

Reconstructing Prehistoric African Population Structure

We assembled genome-wide data from 16 prehistoric Africans. We show that the anciently divergent lineage that comprises the primary ancestry of the southern African San had a wider distribution in the past, contributing approximately two-thirds of the ancestry of Malawi hunter-gatherers ∼8,100-2,500 years ago and approximately one-third of the ancestry of Tanzanian hunter-gatherers ∼1,400 years ago. We document how the spread of farmers from western Africa involved complete replacement of local hunter-gatherers in some regions, and we track the spread of herders by showing that the population of a ∼3,100-year-old pastoralist from Tanzania contributed ancestry to people from northeastern to southern Africa, including a ∼1,200-year-old southern African pastoralist. The deepest diversifications of African lineages were complex, involving either repeated gene flow among geographically disparate groups or a lineage more deeply diverging than that of the San contributing more to some western African populations than to others. We finally leverage ancient genomes to document episodes of natural selection in southern African populations. PAPERCLIP.

DNA damage is a pervasive cause of sequencing errors, directly confounding variant identification

Mutations in somatic cells generate a heterogeneous genomic population and may result in serious medical conditions. Although cancer is typically associated with somatic variations, advances in DNA sequencing indicate that cell-specific variants affect a number of phenotypes and pathologies. Here, we show that mutagenic damage accounts for the majority of the erroneous identification of variants with low to moderate (1 to 5%) frequency. More important, we found signatures of damage in most sequencing data sets in widely used resources, including the 1000 Genomes Project and The Cancer Genome Atlas, establishing damage as a pervasive cause of sequencing errors. The extent of this damage directly confounds the determination of somatic variants in these data sets.

Reducing Artifactual EGFR T790M Mutations in DNA from Formalin-Fixed Paraffin-Embedded Tissue by Use of Thymine-DNA Glycosylase

BACKGROUND

False-positive EGFR T790M mutations have been reported in formalin-fixed lung tumors, but the cause of the false positives has not been identified. The T790M mutation results from a C>T change at the cytosine of a CpG dinucleotide. The presence or absence of methylation at this cytosine has different consequences following deamination, resulting in a thymine or uracil, respectively, both of which however result in an artifactual change. Uracil-DNA glycosylase (UDG) can be used to eliminate DNA templates with uracil residues but is not active against artifactual thymines. We therefore investigated the use of thymine-DNA glycosylase (TDG) to reduce artifactual T790M mutations.

METHODS

Formalin-fixed normal lung tissues and lung squamous cell carcinomas were tested to measure the frequency of false-positive EGFR mutations by use of droplet digital PCR before and after treatment with either UDG or TDG. Methylation at the cytosine at EGFR T790 was assessed by pyrosequencing and by analysis of public databases.

RESULTS

Artifactual EGFR T790M mutations were detected in all of the archival formalin-fixed normal lung and lung squamous cell carcinomas at mutant allele frequencies of 1% or lower. The cytosine at EGFR T790 showed high levels of methylation in all lung cancer samples and normal tissues. Pretreatment of the formalin-fixed DNA with either UDG or TDG reduced the false EGFR T790M mutations, but a greater reduction was seen with the TDG treatment.

CONCLUSIONS

Both U:G and T:G lesions in formalin-fixed tissue are sources of false-positive EGFR T790M mutations. This is the first report of the use of TDG to reduce sequence artifacts in formalin-fixed DNA and is applicable to the accurate detection of mutations arising at methylated cytosines.

Preservation of the metaproteome: variability of protein preservation in ancient dental calculus

Proteomic analysis of dental calculus is emerging as a powerful tool for disease and dietary characterisation of archaeological populations. To better understand the variability in protein results from dental calculus, we analysed 21 samples from three Roman-period populations to compare: 1) the quantity of extracted protein; 2) the number of mass spectral queries; and 3) the number of peptide spectral matches and protein identifications. We found little correlation between the quantity of calculus analysed and total protein identifications, as well as no systematic trends between site location and protein preservation. We identified a wide range of individual variability, which may be associated with the mechanisms of calculus formation and/or post-depositional contamination, in addition to taphonomic factors. Our results suggest dental calculus is indeed a stable, long-term reservoir of proteins as previously reported, but further systematic studies are needed to identify mechanisms associated with protein entrapment and survival in dental calculus.

Mining Metagenomic Data Sets for Ancient DNA: Recommended Protocols for Authentication

While a comparatively young area of research, investigations relying on ancient DNA data have been highly valuable in revealing snapshots of genetic variation in both the recent and the not-so-recent past. Born out of a tradition of single-locus PCR-based approaches that often target individual species, stringent criteria for both data acquisition and analysis were introduced early to establish high standards of data quality. Today, the immense volume of data made available through next-generation sequencing has significantly increased the analytical resolution offered by processing ancient tissues and permits parallel analyses of host and microbial communities. The adoption of this new approach to data acquisition, however, requires an accompanying update on methods of DNA authentication, especially given that ancient molecules are expected to exist in low proportions in archaeological material, where an environmental signal is likely to dominate. In this review, we provide a summary of recent data authentication approaches that have been successfully used to distinguish between endogenous and nonendogenous DNA sequences in metagenomic data sets. While our discussion mostly centers on the detection of ancient human and ancient bacterial pathogen DNA, their applicability is far wider.

Ancient Egyptian mummy genomes suggest an increase of Sub-Saharan African ancestry in post-Roman periods

Egypt, located on the isthmus of Africa, is an ideal region to study historical population dynamics due to its geographic location and documented interactions with ancient civilizations in Africa, Asia and Europe. Particularly, in the first millennium BCE Egypt endured foreign domination leading to growing numbers of foreigners living within its borders possibly contributing genetically to the local population. Here we present 90 mitochondrial genomes as well as genome-wide data sets from three individuals obtained from Egyptian mummies. The samples recovered from Middle Egypt span around 1,300 years of ancient Egyptian history from the New Kingdom to the Roman Period. Our analyses reveal that ancient Egyptians shared more ancestry with Near Easterners than present-day Egyptians, who received additional sub-Saharan admixture in more recent times. This analysis establishes ancient Egyptian mummies as a genetic source to study ancient human history and offers the perspective of deciphering Egypt's past at a genome-wide level.

Archaeological and historical records had shown ancient Egypt before and after Ptolemaic and Roman periods to be a hub of human migration and exchange. Here, Schuenemann and colleagues analyse ancient mitochondrial and nuclear DNA to investigate the genetic history of Egypt.

Ancient genomic changes associated with domestication of the horse

The domestication of the horse was a seminal event in human cultural evolution. Libradoet al.obtained genome sequences from 14 horses from the Bronze and Iron Ages, about 2000 to 4000 years ago, soon after domestication. They identified variants determining coat color and genes selected during the domestication process. They could also see evidence of admixture with archaic horses and the demography of the domestication process, which included the accumulation of deleterious variants. The horse appears to have undergone a different type of domestication process than animals that were domesticated simply for food.

Science, this issue p.442

Improving ancient DNA genome assembly

Most reconstruction methods for genomes of ancient origin that are used today require a closely related reference. In order to identify genomic rearrangements or the deletion of whole genes, de novo assembly has to be used. However, because of inherent problems with ancient DNA, its de novo assembly is highly complicated. In order to tackle the diversity in the length of the input reads, we propose a two-layer approach, where multiple assemblies are generated in the first layer, which are then combined in the second layer. We used this two-layer assembly to generate assemblies for two different ancient samples and compared the results to current de novo assembly approaches. We are able to improve the assembly with respect to the length of the contigs and can resolve more repetitive regions.

Reinforcing plant evolutionary genomics using ancient DNA

Improved understanding of ancient DNA (aDNA) biochemical properties coupled with application of next generation sequencing (NGS) methods enabled sequencing and authenticating genomes of historical samples. This advancement ignited a revolution in plant evolutionary genomics by allowing direct observations of past molecular diversity. Analyses of genomes sequenced from temporally distributed samples of Gossypium sp., Phytophthora infestans and Arabidopsis thaliana improved our understanding of the evolutionary rates and time scales at which genome remodeling takes place. Comparison of historical samples of barley (Hordeum vulgare) and maize (Zea mays ssp. mays) with their present-day counterparts enabled assessment of selection during different stages of domestication. These examples show how aDNA already improved our evolutionary inferences. Increasing quality and amount of sequencing data retrieved from historical plants will further advance our understanding of plant evolution.

Comparing Ancient DNA Preservation in Petrous Bone and Tooth Cementum

Large-scale genomic analyses of ancient human populations have become feasible partly due to refined sampling methods. The inner part of petrous bones and the cementum layer in teeth roots are currently recognized as the best substrates for such research. We present a comparative analysis of DNA preservation in these two substrates obtained from the same human skulls, across a range of different ages and preservation environments. Both substrates display significantly higher endogenous DNA content (average of 16.4% and 40.0% for teeth and petrous bones, respectively) than parietal skull bone (average of 2.2%). Despite sample-to-sample variation, petrous bone overall performs better than tooth cementum (p = 0.001). This difference, however, is driven largely by a cluster of viking skeletons from one particular locality, showing relatively poor molecular tooth preservation (<10% endogenous DNA). In the remaining skeletons there is no systematic difference between the two substrates. A crude preservation (good/bad) applied to each sample prior to DNA-extraction predicted the above/below 10% endogenous DNA threshold in 80% of the cases. Interestingly, we observe signficantly higher levels of cytosine to thymine deamination damage and lower proportions of mitochondrial/nuclear DNA in petrous bone compared to tooth cementum. Lastly, we show that petrous bones from ancient cremated individuals contain no measurable levels of authentic human DNA. Based on these findings we discuss the pros and cons of sampling the different elements.

A High-Coverage Yersinia pestis Genome from a Sixth-Century Justinianic Plague Victim

The Justinianic Plague, which started in the sixth century and lasted to the mid eighth century, is thought to be the first of three historically documented plague pandemics causing massive casualties. Historical accounts and molecular data suggest the bacterium Yersinia pestis as its etiological agent. Here we present a new high-coverage (17.9-fold) Y. pestis genome obtained from a sixth-century skeleton recovered from a southern German burial site close to Munich. The reconstructed genome enabled the detection of 30 unique substitutions as well as structural differences that have not been previously described. We report indels affecting a lacl family transcription regulator gene as well as nonsynonymous substitutions in the nrdE, fadJ, and pcp genes, that have been suggested as plague virulence determinants or have been shown to be upregulated in different models of plague infection. In addition, we identify 19 false positive substitutions in a previously published lower-coverage Y. pestis genome from another archaeological site of the same time period and geographical region that is otherwise genetically identical to the high-coverage genome sequence reported here, suggesting low-genetic diversity of the plague during the sixth century in rural southern Germany.

Native American Genomics and Population Histories

Studies of Native American genetic diversity and population history have been transformed over the last decade by important developments in anthropological genetics. During this time, researchers have adopted new DNA technologies and computational approaches for analyzing genomic data, and they have become increasingly sensitive to the views of research participants and communities. As new methods are applied to long-standing questions, and as more research is conducted in collaboration with indigenous communities, we are gaining new insights into the history and diversity of indigenous populations. This review discusses the recent methodological advances and genetic studies that have improved our understanding of Native American genomics and population histories. We synthesize current knowledge about Native American genomic variation and build a model of population history in the Americas. We also discuss the broader implications of this research for anthropology and related disciplines, and we highlight challenges and other considerations for future research.

Fast, Accurate and Automatic Ancient Nucleosome and Methylation Maps with epiPALEOMIX

The first epigenomes from archaic hominins (AH) and ancient anatomically modern humans (AMH) have recently been characterized, based, however, on a limited number of samples. The extent to which ancient genome-wide epigenetic landscapes can be reconstructed thus remains contentious. Here, we present epiPALEOMIX, an open-source and user-friendly pipeline that exploits post-mortem DNA degradation patterns to reconstruct ancient methylomes and nucleosome maps from shotgun and/or capture-enrichment data. Applying epiPALEOMIX to the sequence data underlying 35 ancient genomes including AMH, AH, equids and aurochs, we investigate the temporal, geographical and preservation range of ancient epigenetic signatures. We first assess the quality of inferred ancient epigenetic signatures within well-characterized genomic regions. We find that tissue-specific methylation signatures can be obtained across a wider range of DNA preparation types than previously thought, including when no particular experimental procedures have been used to remove deaminated cytosines prior to sequencing. We identify a large subset of samples for which DNA associated with nucleosomes is protected from post-mortem degradation, and nucleosome positioning patterns can be reconstructed. Finally, we describe parameters and conditions such as DNA damage levels and sequencing depth that limit the preservation of epigenetic signatures in ancient samples. When such conditions are met, we propose that epigenetic profiles of CTCF binding regions can be used to help data authentication. Our work, including epiPALEOMIX, opens for further investigations of ancient epigenomes through time especially aimed at tracking possible epigenetic changes during major evolutionary, environmental, socioeconomic, and cultural shifts.

DNA Sequencing in Cultural Heritage

During the last three decades, DNA analysis on degraded samples revealed itself as an important research tool in anthropology, archaeozoology, molecular evolution, and population genetics. Application on topics such as determination of species origin of prehistoric and historic objects, individual identification of famous personalities, characterization of particular samples important for historical, archeological, or evolutionary reconstructions, confers to the paleogenetics an important role also for the enhancement of cultural heritage. A really fast improvement in methodologies in recent years led to a revolution that permitted recovering even complete genomes from highly degraded samples with the possibility to go back in time 400,000 years for samples from temperate regions and 700,000 years for permafrozen remains and to analyze even more recent material that has been subjected to hard biochemical treatments. Here we propose a review on the different methodological approaches used so far for the molecular analysis of degraded samples and their application on some case studies.

Systematic relationships of five newly sequenced cervid species

Cervid phylogenetics has been puzzling researchers for over 150 years. In recent decades, molecular systematics has provided new input for both the support and revision of the previous results from comparative anatomy but has led to only partial consensus. Despite all of the efforts to reach taxon-wide species sampling over the last two decades, a number of cervid species still lack molecular data because they are difficult to access in the wild. By extracting ancient DNA from museum specimens, in this study, we obtained partial mitochondrial cytochrome b gene sequences for Mazama bricenii, Mazama chunyi, Muntiacus atherodes, Pudu mephistophiles, and Rusa marianna, including three holotypes. These new sequences were used to enrich the existing mitochondrial DNA alignments and yielded the most taxonomically complete data set for cervids to date. Phylogenetic analyses provide new insights into the evolutionary history of these five species. However, systematic uncertainties within Muntiacus persist and resolving phylogenetic relationships within Pudu and Mazama remain challenging.

Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel

The performance of hybridization capture combined with next‐generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient‐domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187‐fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient‐domestic dromedaries with 17–95% length coverage and 1.27–47.1‐fold read depths for the covered regions. Using whole‐genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1–1.06‐fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.

Artifactual mutations resulting from DNA lesions limit detection levels in ultrasensitive sequencing applications

The need in cancer research or evolutionary biology to detect rare mutations or variants present at very low frequencies (<10⁻⁵) poses an increasing demand on lowering the detection limits of available methods. Here we demonstrated that amplifiable DNA lesions introduce important error sources in ultrasensitive technologies such as single molecule PCR (smPCR) applications (e.g. droplet-digital PCR), or next-generation sequencing (NGS) based methods. Using templates with known amplifiable lesions (8-oxoguanine, deaminated 5-methylcytosine, uracil, and DNA heteroduplexes), we assessed with smPCR and duplex sequencing that templates with these lesions were amplified very efficiently by proofreading polymerases (except uracil), leading to G->T, and to a lesser extent, to unreported G->C substitutions at 8-oxoguanine lesions, and C->T transitions in amplified uracil containing templates. Long heat incubations common in many DNA extraction protocols significantly increased the number of G->T substitutions. Moreover, in ∼50-80% smPCR reactions we observed the random amplification preference of only one of both DNA strands explaining the known ‘PCR jackpot effect’, with the result that a lesion became indistinguishable from a true mutation or variant. Finally, we showed that artifactual mutations derived from uracil and 8-oxoguanine could be significantly reduced by DNA repair enzymes.

Opsin gene repertoires in northern archaic hominids

The Neanderthals' northern distribution, hunting techniques, and orbit breadths suggest that they were more active in dim light than modern humans. We surveyed visual opsin genes from four Neanderthals and two other archaic hominids to see if they provided additional support for this hypothesis. This analysis was motivated by the observation that alleles responsible for anomalous trichromacy in humans are more common in northern latitudes, by data suggesting that these variants might enhance vision in mesopic conditions, and by the observation that dim light active species often have fewer opsin genes than diurnal relatives. We also looked for evidence of convergent amino acid substitutions in Neanderthal opsins and orthologs from crepuscular or nocturnal species. The Altai Neanderthal, the Denisovan, and the Ust'-Ishim early modern human had opsin genes that encoded proteins identical to orthologs in the human reference genome. Opsins from the Vindija Cave Neanderthals (three females) had many nonsynonymous substitutions, including several predicted to influence colour vision (e.g., stop codons). However, the functional implications of these observations were difficult to assess, given that "control" loci, where no substitutions were expected, differed from humans to the same extent. This left unresolved the test for colour vision deficiencies in Vindija Cave Neanderthals.

Genomic insights into the origin of farming in the ancient Near East

We report genome-wide ancient DNA from 44 ancient Near Easterners ranging in time between ~12,000 and 1,400 bc, from Natufian hunter-gatherers to Bronze Age farmers. We show that the earliest populations of the Near East derived around half their ancestry from a 'Basal Eurasian' lineage that had little if any Neanderthal admixture and that separated from other non-African lineages before their separation from each other. The first farmers of the southern Levant (Israel and Jordan) and Zagros Mountains (Iran) were strongly genetically differentiated, and each descended from local hunter-gatherers. By the time of the Bronze Age, these two populations and Anatolian-related farmers had mixed with each other and with the hunter-gatherers of Europe to greatly reduce genetic differentiation. The impact of the Near Eastern farmers extended beyond the Near East: farmers related to those of Anatolia spread westward into Europe; farmers related to those of the Levant spread southward into East Africa; farmers related to those of Iran spread northward into the Eurasian steppe; and people related to both the early farmers of Iran and to the pastoralists of the Eurasian steppe spread eastward into South Asia.

Toward high-resolution population genomics using archaeological samples

The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

Genomic analysis of 6,000-year-old cultivated grain illuminates the domestication history of barley

The cereal grass barley was domesticated about 10,000 years before the present in the Fertile Crescent and became a founder crop of Neolithic agriculture. Here we report the genome sequences of five 6,000-year-old barley grains excavated at a cave in the Judean Desert close to the Dead Sea. Comparison to whole-exome sequence data from a diversity panel of present-day barley accessions showed the close affinity of ancient samples to extant landraces from the Southern Levant and Egypt, consistent with a proposed origin of domesticated barley in the Upper Jordan Valley. Our findings suggest that barley landraces grown in present-day Israel have not experienced major lineage turnover over the past six millennia, although there is evidence for gene flow between cultivated and sympatric wild populations. We demonstrate the usefulness of ancient genomes from desiccated archaeobotanical remains in informing research into the origin, early domestication and subsequent migration of crop species.

Epigenetics of Ancient DNA

Initially, the study of DNA isolated from ancient specimens had been based on the analysis of the primary nucleotide sequence. This approach has allowed researchers to study the evolutionary changes that occur in different populations and determine the influence of the environment on genetic selection. However, the improvement of methodological approaches to genome-wide analysis has opened up new possibilities in the search for the epigenetic mechanisms involved in the regulation of gene expression. It was discovered recently that the methylation status of the regulatory elements of the HOXD cluster and MEIS1 gene changed during human evolution. Epigenetic changes in these genes played a key role in the evolution of the limbs of modern humans. Recent works have demonstrated that it is possible to determine the transcriptional activity of genes in ancient DNA samples by combining information on DNA methylation and the DNAaseI hypersensitive sequences located at the transcription start sites of genes. In the nearest future, if a preserved fossils brain is found, it will be possible to identify the evolutionary changes in the higher nervous system associated with epigenetic differences.

Temporal patterns of damage and decay kinetics of DNA retrieved from plant herbarium specimens

Herbaria archive a record of changes of worldwide plant biodiversity harbouring millions of specimens that contain DNA suitable for genome sequencing. To profit from this resource, it is fundamental to understand in detail the process of DNA degradation in herbarium specimens. We investigated patterns of DNA fragmentation and nucleotide misincorporation by analysing 86 herbarium samples spanning the last 300 years using Illumina shotgun sequencing. We found an exponential decay relationship between DNA fragmentation and time, and estimated a per nucleotide fragmentation rate of 1.66 × 10(-4) per year, which is six times faster than the rate estimated for ancient bones. Additionally, we found that strand breaks occur specially before purines, and that depurination-driven DNA breakage occurs constantly through time and can to a great extent explain decreasing fragment length over time. Similar to what has been found analysing ancient DNA from bones, we found a strong correlation between the deamination-driven accumulation of cytosine to thymine substitutions and time, which reinforces the importance of substitution patterns to authenticate the ancient/historical nature of DNA fragments. Accurate estimations of DNA degradation through time will allow informed decisions about laboratory and computational procedures to take advantage of the vast collection of worldwide herbarium specimens.

Joint Estimation of Contamination, Error and Demography for Nuclear DNA from Ancient Humans

When sequencing an ancient DNA sample from a hominin fossil, DNA from present-day humans involved in excavation and extraction will be sequenced along with the endogenous material. This type of contamination is problematic for downstream analyses as it will introduce a bias towards the population of the contaminating individual(s). Quantifying the extent of contamination is a crucial step as it allows researchers to account for possible biases that may arise in downstream genetic analyses. Here, we present an MCMC algorithm to co-estimate the contamination rate, sequencing error rate and demographic parameters—including drift times and admixture rates—for an ancient nuclear genome obtained from human remains, when the putative contaminating DNA comes from present-day humans. We assume we have a large panel representing the putative contaminant population (e.g. European, East Asian or African). The method is implemented in a C++ program called ‘Demographic Inference with Contamination and Error’ (DICE). We applied it to simulations and genome data from ancient Neanderthals and modern humans. With reasonable levels of genome sequence coverage (>3X), we find we can recover accurate estimates of all these parameters, even when the contamination rate is as high as 50%.

When extracting and sequencing ancient DNA from human remains, a recurrent problem is the presence of DNA from the paleontologists, archaeologists or geneticists that may have handled the fossil. If a DNA library is highly contaminated, this will introduce biases in downstream analyses, so it is important to determine the amount of extraneous DNA. Different methods exist for this purpose, but few are applicable to the nuclear genome, and none of them can extract reliable genomic information from highly contaminated samples. Thus, samples with high rates of contamination are usually discarded. Here, we present a method to jointly estimate contamination and error rates, along with demographic parameters, like drift times and admixture rates. Our method can serve to uncover important details about the evolutionary history of archaic and early modern humans from ancient DNA samples, even if those samples are highly contaminated.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

The major histocompatibility complex in Old World camelids and low polymorphism of its class II genes

BACKGROUND

The Major Histocompatibility Complex (MHC) is a genomic region containing genes with crucial roles in immune responses. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. To counteract the high variability of pathogens, the MHC evolved into a region of considerable heterogeneity in its organization, number and extent of polymorphism. Studies of MHCs in different model species contribute to our understanding of mechanisms of immunity, diseases and their evolution. Camels are economically important domestic animals and interesting biomodels. Three species of Old World camels have been recognized: the dromedary (Camelus dromedarius), Bactrian camel (Camelus bactrianus) and the wild camel (Camelus ferus). Despite their importance, little is known about the MHC genomic region, its organization and diversity in camels. The objectives of this study were to identify, map and characterize the MHC region of Old World camelids, with special attention to genetic variation at selected class MHC II loci.

RESULTS

Physical mapping located the MHC region to the chromosome 20 in Camelus dromedarius. Cytogenetic and comparative analyses of whole genome sequences showed that the order of the three major sub-regions is "Centromere - Class II - Class III - Class I". DRA, DRB, DQA and DQB exon 2 sequences encoding the antigen binding site of the corresponding class II antigen presenting molecules showed high degree of sequence similarity and extensive allele sharing across the three species. Unexpectedly low extent of polymorphism with low numbers of alleles and haplotypes was observed in all species, despite different geographic origins of the camels analyzed. The DRA locus was found to be polymorphic, with three alleles shared by all three species. DRA and DQA sequences retrieved from ancient DNA samples of Camelus dromedarius suggested that additional polymorphism might exist.

CONCLUSIONS

This study provided evidence that camels possess an MHC comparable to other mammalian species in terms of its genomic localization, organization and sequence similarity. We described ancient variation at the DRA locus, monomorphic in most species. The extent of molecular diversity of MHC class II genes seems to be substantially lower in Old World camels than in other mammalian species.

Ancient gene flow from early modern humans into Eastern Neanderthals

It has been shown that Neanderthals contributed genetically to modern humans outside Africa 47,000-65,000 years ago. Here we analyse the genomes of a Neanderthal and a Denisovan from the Altai Mountains in Siberia together with the sequences of chromosome 21 of two Neanderthals from Spain and Croatia. We find that a population that diverged early from other modern humans in Africa contributed genetically to the ancestors of Neanderthals from the Altai Mountains roughly 100,000 years ago. By contrast, we do not detect such a genetic contribution in the Denisovan or the two European Neanderthals. We conclude that in addition to later interbreeding events, the ancestors of Neanderthals from the Altai Mountains and early modern humans met and interbred, possibly in the Near East, many thousands of years earlier than previously thought.