Reconstructing ancient genomes and epigenomes

Century-old chromatin architecture revealed in formalin-fixed vertebrates

Gene expression is regulated by changes in chromatin architecture intrinsic to cellular differentiation and as an active response to environmental stimuli. Chromatin dynamics are a major driver of phenotypic diversity, regulation of development, and manifestation of disease. Remarkably, we know little about the evolutionary dynamics of chromatin reorganisation through time, data essential to characterise the impact of environmental stress during the ongoing biodiversity extinction crisis (20th-21st century). Linking the disparate fields of chromatin biology and museum science through their common use of the preservative formaldehyde (a constituent of formalin), we have generated historical chromatin profiles in museum specimens up to 117 years old. Historical chromatin profiles are reproducible, tissue-specific, sex-specific, and environmental condition-dependent in vertebrate specimens. Additionally, we show that over-fixation modulates differential chromatin accessibility to enable semi-quantitative estimates of relative gene expression in vertebrates and a yeast model. Our approach transforms formalin-fixed biological collections into an accurate, comprehensive, and global record of environmental impact on gene expression and phenotype.

Generation of genome-wide SNP markers from minimally invasive sampling in endangered animals and applications in species ecology and conservation

High-density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high-quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA.

Investigating food production-associated DNA methylation changes in paleogenomes: Lack of consistent signals beyond technical noise

The Neolithic transition introduced major diet and lifestyle changes to human populations across continents. Beyond well-documented bioarcheological and genetic effects, whether these changes also had molecular-level epigenetic repercussions in past human populations has been an open question. In fact, methylation signatures can be inferred from UDG-treated ancient DNA through postmortem damage patterns, but with low signal-to-noise ratios; it is thus unclear whether published paleogenomes would provide the necessary resolution to discover systematic effects of lifestyle and diet shifts. To address this we compiled UDG-treated shotgun genomes of 13 pre-Neolithic hunter-gatherers (HGs) and 21 Neolithic farmers (NFs) individuals from West and North Eurasia, published by six different laboratories and with coverage c.1×-58× (median = 9×). We used epiPALEOMIX and a Monte Carlo normalization scheme to estimate methylation levels per genome. Our paleomethylome dataset showed expected genome-wide methylation patterns such as CpG island hypomethylation. However, analyzing the data using various approaches did not yield any systematic signals for subsistence type, genetic sex, or tissue effects. Comparing the HG-NF methylation differences in our dataset with methylation differences between hunter-gatherers versus farmers in modern-day Central Africa also did not yield consistent results. Meanwhile, paleomethylome profiles did cluster strongly by their laboratories of origin. Using larger data volumes, minimizing technical noise and/or using alternative protocols may be necessary for capturing subtle environment-related biological signals from paleomethylomes.

Exploring the complex pre-adaptations of invasive plants to anthropogenic disturbance: a call for integration of archaeobotanical approaches

Pre-adaptation to anthropogenic disturbance is broadly considered key for plant invasion success. Nevertheless, empirical evidence remains scarce and fragmentary, given the multifaceted nature of anthropogenic disturbance itself and the complexity of other evolutionary forces shaping the (epi)-genomes of recent native and invasive plant populations. Here, we review and critically revisit the existing theory and empirical evidence in the field of evolutionary ecology and highlight novel integrative research avenues that work at the interface with archaeology to solve open questions. The approaches suggested so far focus on contemporary plant populations, although their genomes have rapidly changed since their initial introduction in response to numerous selective and stochastic forces. We elaborate that a role of pre-adaptation to anthropogenic disturbance in plant invasion success should thus additionally be validated based on the analyses of archaeobotanical remains. Such materials, in the light of detailed knowledge on past human societies could highlight fine-scale differences in the type and timing of past disturbances. We propose a combination of archaeobotanical, ancient DNA and morphometric analyses of plant macro- and microremains to assess past community composition, and species' functional traits to unravel the timing of adaptation processes, their drivers and their long-term consequences for invasive species. Although such methodologies have proven to be feasible for numerous crop plants, they have not been yet applied to wild invasive species, which opens a wide array of insights into their evolution.

Search for differentially methylated regions in ancient and modern genomes

Currently, active research is focused on investigating the mechanisms that regulate the development of various pathologies and their evolutionary dynamics. Epigenetic mechanisms, such as DNA methylation, play a significant role in evolutionary processes, as their changes have a faster impact on the phenotype compared to mutagenesis. In this study, we attempted to develop an algorithm for identifying differentially methylated regions associated with metabolic syndrome, which have undergone methylation changes in humans during the transition from a hunter-gatherer to a sedentary lifestyle. The application of existing whole-genome bisulfite sequencing methods is limited for ancient samples due to their low quality and fragmentation, and the approach to obtaining DNA methylation profiles differs significantly between ancient hunter-gatherer samples and modern tissues. In this study, we validated DamMet, an algorithm for reconstructing ancient methylomes. Application of DamMet to Neanderthal and Denisovan genomes showed a moderate level of correlation with previously published methylation profiles and demonstrated an underestimation of methylation levels in the reconstructed profiles by an average of 15-20 %. Additionally, we developed a new Python-based algorithm that allows for the comparison of methylomes in ancient and modern samples, despite the absence of methylation profiles in modern bone tissue within the context of obesity. This analysis involves a two-step data processing approach, where the first step involves the identification and filtration of tissue-specific methylation regions, and the second step focuses on the direct search for differentially methylated regions in specific areas associated with the researcher's target condition. By applying this algorithm to test data, we identified 38 differentially methylated regions associated with obesity, the majority of which were located in promoter regions. The pipeline demonstrated sufficient efficiency in detecting these regions. These results confirm the feasibility of reconstructing DNA methylation profiles in ancient samples and comparing them with modern methylomes. Furthermore, possibilities for further methodological development and the implementation of a new step for studying differentially methylated positions associated with evolutionary processes are discussed.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

An ecologist's guide for studying DNA methylation variation in wild vertebrates

The field of molecular biology is advancing fast with new powerful technologies, sequencing methods and analysis software being developed constantly. Commonly used tools originally developed for research on humans and model species are now regularly used in ecological and evolutionary research. There is also a growing interest in the causes and consequences of epigenetic variation in natural populations. Studying ecological epigenetics is currently challenging, especially for vertebrate systems, because of the required technical expertise, complications with analyses and interpretation, and limitations in acquiring sufficiently high sample sizes. Importantly, neglecting the limitations of the experimental setup, technology and analyses may affect the reliability and reproducibility, and the extent to which unbiased conclusions can be drawn from these studies. Here, we provide a practical guide for researchers aiming to study DNA methylation variation in wild vertebrates. We review the technical aspects of epigenetic research, concentrating on DNA methylation using bisulfite sequencing, discuss the limitations and possible pitfalls, and how to overcome them through rigid and reproducible data analysis. This review provides a solid foundation for the proper design of epigenetic studies, a clear roadmap on the best practices for correct data analysis and a realistic view on the limitations for studying ecological epigenetics in vertebrates. This review will help researchers studying the ecological and evolutionary implications of epigenetic variation in wild populations.

A Large-Scale Study into Protist-Animal Interactions Based on Public Genomic Data Using DNA Barcodes

With the birth of next-generation sequencing (NGS) technology, genomic data in public databases have increased exponentially. Unfortunately, exogenous contamination or intracellular parasite sequences in assemblies could confuse genomic analysis. Meanwhile, they can provide a valuable resource for studies of host-microbe interactions. Here, we used a strategy based on DNA barcodes to scan protistan contamination in the GenBank WGS/TSA database. The results showed a total of 13,952 metazoan/animal assemblies in GenBank, where 17,036 contigs were found to be protistan contaminants in 1507 assemblies (10.8%), with even higher contamination rates in taxa of Cnidaria (150/281), Crustacea (237/480), and Mollusca (107/410). Taxonomic analysis of the protists derived from these contigs showed variations in abundance and evenness of protistan contamination across different metazoan taxa, reflecting host preferences of Apicomplexa, Ciliophora, Oomycota and Symbiodiniaceae for mammals and birds, Crustacea, insects, and Cnidaria, respectively. Finally, mitochondrial proteins COX1 and CYTB were predicted from these contigs, and the phylogenetic analysis corroborated the protistan origination and heterogeneous distribution of the contaminated contigs. Overall, in this study, we conducted a large-scale scan of protistan contaminant in genomic resources, and the protistan sequences detected will help uncover the protist diversity and relationships of these picoeukaryotes with Metazoa.

Evolving Perspective on the Origin and Diversification of Cellular Life and the Virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

SPIN enables high throughput species identification of archaeological bone by proteomics

Species determination based on genetic evidence is an indispensable tool in archaeology, forensics, ecology, and food authentication. Most available analytical approaches involve compromises with regard to the number of detectable species, high cost due to low throughput, or a labor-intensive manual process. Here, we introduce "Species by Proteome INvestigation" (SPIN), a shotgun proteomics workflow for analyzing archaeological bone capable of querying over 150 mammalian species by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Rapid peptide chromatography and data-independent acquisition (DIA) with throughput of 200 samples per day reduce expensive MS time, whereas streamlined sample preparation and automated data interpretation save labor costs. We confirm the successful classification of known reference bones, including domestic species and great apes, beyond the taxonomic resolution of the conventional peptide mass fingerprinting (PMF)-based Zooarchaeology by Mass Spectrometry (ZooMS) method. In a blinded study of degraded Iron-Age material from Scandinavia, SPIN produces reproducible results between replicates, which are consistent with morphological analysis. Finally, we demonstrate the high throughput capabilities of the method in a high-degradation context by analyzing more than two hundred Middle and Upper Palaeolithic bones from Southern European sites with late Neanderthal occupation. While this initial study is focused on modern and archaeological mammalian bone, SPIN will be open and expandable to other biological tissues and taxa.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Epigenomic Modifications in Modern and Ancient Genomes

Epigenetic changes have been identified as a major driver of fundamental metabolic pathways. More specifically, the importance of epigenetic regulatory mechanisms for biological processes like speciation and embryogenesis has been well documented and revealed the direct link between epigenetic modifications and various diseases. In this review, we focus on epigenetic changes in animals with special attention on human DNA methylation utilizing ancient and modern genomes. Acknowledging the latest developments in ancient DNA research, we further discuss paleoepigenomic approaches as the only means to infer epigenetic changes in the past. Investigating genome-wide methylation patterns of ancient humans may ultimately yield in a more comprehensive understanding of how our ancestors have adapted to the changing environment, and modified their lifestyles accordingly. We discuss the difficulties of working with ancient DNA in particular utilizing paleoepigenomic approaches, and assess new paleoepigenomic data, which might be helpful in future studies.

Extraction and high-throughput sequencing of oak heartwood DNA: Assessing the feasibility of genome-wide DNA methylation profiling

Tree ring features are affected by environmental factors and therefore are the basis for dendrochronological studies to reconstruct past environmental conditions. Oak wood often provides the data for these studies because of the durability of oak heartwood and hence the availability of samples spanning long time periods of the distant past. Wood formation is regulated in part by epigenetic mechanisms such as DNA methylation. Studies of the methylation state of DNA preserved in oak heartwood thus could identify epigenetic tree ring features informing on past environmental conditions. In this study, we aimed to establish protocols for the extraction of DNA, the high-throughput sequencing of whole-genome DNA libraries (WGS) and the profiling of DNA methylation by whole-genome bisulfite sequencing (WGBS) for oak (Quercus robur) heartwood drill cores taken from the trunks of living standing trees spanning the AD 1776-2014 time period. Heartwood contains little DNA, and large amounts of phenolic compounds known to hinder the preparation of high-throughput sequencing libraries. Whole-genome and DNA methylome library preparation and sequencing consistently failed for oak heartwood samples more than 100 and 50 years of age, respectively. DNA fragmentation increased with sample age and was exacerbated by the additional bisulfite treatment step during methylome library preparation. Relative coverage of the non-repetitive portion of the oak genome was sparse. These results suggest that quantitative methylome studies of oak hardwood will likely be limited to relatively recent samples and will require a high sequencing depth to achieve sufficient genome coverage.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Unraveling Human Brain Development and Evolution Using Organoid Models

Brain organoids are proving to be physiologically relevant models for studying human brain development in terms of temporal transcriptional signature recapitulation, dynamic cytoarchitectural development, and functional electrophysiological maturation. Several studies have employed brain organoid technologies to elucidate human-specific processes of brain development, gene expression, and cellular maturation by comparing human-derived brain organoids to those of non-human primates (NHPs). Brain organoids have been established from a variety of NHP pluripotent stem cell (PSC) lines and many protocols are now available for generating brain organoids capable of reproducibly representing specific brain region identities. Innumerous combinations of brain region specific organoids derived from different human and NHP PSCs, with CRISPR-Cas9 gene editing techniques and strategies to promote advanced stages of maturation, will successfully establish complex brain model systems for the accurate representation and elucidation of human brain development. Identified human-specific processes of brain development are likely vulnerable to dysregulation and could result in the identification of therapeutic targets or disease prevention strategies. Here, we discuss the potential of brain organoids to successfully model human-specific processes of brain development and explore current strategies for pinpointing these differences.

Uncovering the genomic and metagenomic research potential in old ethanol-preserved snakes

Natural history museum collections worldwide represent a tremendous resource of information on past and present biodiversity. Fish, reptiles, amphibians and many invertebrate collections have often been preserved in ethanol for decades or centuries and our knowledge on the genomic and metagenomic research potential of such material is limited. Here, we use ancient DNA protocols, combined with shotgun sequencing to test the molecular preservation in liver, skin and bone tissue from five old (1842 to 1964) museum specimens of the common garter snake (Thamnophis sirtalis). When mapping reads to a T. sirtalis reference genome, we find that the DNA molecules are highly damaged with short average sequence lengths (38-64 bp) and high C-T deamination, ranging from 9% to 21% at the first position. Despite this, the samples displayed relatively high endogenous DNA content, ranging from 26% to 56%, revealing that genome-scale analyses are indeed possible from all specimens and tissues included here. Of the three tested types of tissue, bone shows marginally but significantly higher DNA quality in these metrics. Though at least one of the snakes had been exposed to formalin, neither the concentration nor the quality of the obtained DNA was affected. Lastly, we demonstrate that these specimens display a diverse and tissue-specific microbial genetic profile, thus offering authentic metagenomic data despite being submerged in ethanol for many years. Our results emphasize that historical museum collections continue to offer an invaluable source of information in the era of genomics.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

DNA methylation profiling in mummified human remains from the eighteenth-century

Reconstruction of ancient epigenomes by DNA methylation (DNAm) can shed light into the composition of cell types, disease states, and age at death. However, such analysis is hampered by impaired DNA quality and little is known how decomposition affects DNAm. In this study, we determined if EPIC Illumina BeadChip technology is applicable for specimens from mummies of the eighteenth century CE. Overall, the signal intensity on the microarray was extremely low, but for one of two samples we were able to detect characteristic DNAm signals in a subset of CG dinucleotides (CpGs), which were selected with a stringent processing pipeline. Using only these CpGs we could train epigenetic signatures with reference DNAm profiles of multiple tissues and our predictions matched the fact that the specimen was lung tissue from a 28-year-old woman. Thus, we provide proof of principle that Illumina BeadChips are applicable for DNAm profiling in ancient samples.

Ancient DNA analysis of rare genetic bone disorders

OBJECTIVE

Review of the current advancements in the field of paleogenetics that provide new opportunities in studying the evolution of rare genetic bone diseases.

MATERIAL AND METHODS

Based on cases from the literature, the genetics of rare bone diseases will be introduced and the main methodological issues will be addressed, focusing on the opportunities presented by the application of aDNA analyses in the field of paleopathology.

RESULTS

Medical literature provides large datasets on the genes responsible for rare bone disorders. These genes, subdivided in functional categories, display important future targets when analyzing rare genetic bone disorders in ancient human remains.

CONCLUSIONS

Knowledge on both phenotype and genotype is required to study rare diseases in ancient human remains.

SIGNIFICANCE

The proposed interdisciplinary research will provide new insight into the occurrence and spread of genetic risk factors in the past and will help in the diagnostics of these rare and often neglected diseases.

LIMITATIONS

The current limitations in ancient DNA research and targeting the disease-causing specific mutations (e.g., somatic or germline).

SUGGESTIONS FOR FURTHER RESEARCH

Methodological advancements and candidate gene lists provide the optimal basis for future interdisciplinary studies of rare genetic bone disorders in ancient human remains.

Genomic reconstruction of fossil and living microorganisms in ancient Siberian permafrost

BACKGROUND

Total DNA (intracellular, iDNA and extracellular, eDNA) from ancient permafrost records the mixed genetic repository of the past and present microbial populations through geological time. Given the exceptional preservation of eDNA under perennial frozen conditions, typical metagenomic sequencing of total DNA precludes the discrimination between fossil and living microorganisms in ancient cryogenic environments. DNA repair protocols were combined with high throughput sequencing (HTS) of separate iDNA and eDNA fraction to reconstruct metagenome-assembled genomes (MAGs) from ancient microbial DNA entrapped in Siberian coastal permafrost.

RESULTS

Despite the severe DNA damage in ancient permafrost, the coupling of DNA repair and HTS resulted in a total of 52 MAGs from sediments across a chronosequence (26-120 kyr). These MAGs were compared with those derived from the same samples but without utilizing DNA repair protocols. The MAGs from the youngest stratum showed minimal DNA damage and thus likely originated from viable, active microbial species. Many MAGs from the older and deeper sediment appear related to past aerobic microbial populations that had died upon freezing. MAGs from anaerobic lineages, including Asgard archaea, however exhibited minimal DNA damage and likely represent extant living microorganisms that have become adapted to the cryogenic and anoxic environments. The integration of aspartic acid racemization modeling and metaproteomics further constrained the metabolic status of the living microbial populations. Collectively, combining DNA repair protocols with HTS unveiled the adaptive strategies of microbes to long-term survivability in ancient permafrost.

CONCLUSIONS

Our results indicated that coupling of DNA repair protocols with simultaneous sequencing of iDNA and eDNA fractions enabled the assembly of MAGs from past and living microorganisms in ancient permafrost. The genomic reconstruction from the past and extant microbial populations expanded our understanding about the microbial successions and biogeochemical alterations from the past paleoenvironment to the present-day frozen state. Furthermore, we provided genomic insights into long-term survival mechanisms of microorganisms under cryogenic conditions through geological time. The combined strategies in this study can be extrapolated to examine other ancient non-permafrost environments and constrain the search for past and extant extraterrestrial life in permafrost and ice deposits on Mars. Video abstract.

The genomic history of the Aegean palatial civilizations

The Cycladic, the Minoan, and the Helladic (Mycenaean) cultures define the Bronze Age (BA) of Greece. Urbanism, complex social structures, craft and agricultural specialization, and the earliest forms of writing characterize this iconic period. We sequenced six Early to Middle BA whole genomes, along with 11 mitochondrial genomes, sampled from the three BA cultures of the Aegean Sea. The Early BA (EBA) genomes are homogeneous and derive most of their ancestry from Neolithic Aegeans, contrary to earlier hypotheses that the Neolithic-EBA cultural transition was due to massive population turnover. EBA Aegeans were shaped by relatively small-scale migration from East of the Aegean, as evidenced by the Caucasus-related ancestry also detected in Anatolians. In contrast, Middle BA (MBA) individuals of northern Greece differ from EBA populations in showing ∼50% Pontic-Caspian Steppe-related ancestry, dated at ca. 2,600-2,000 BCE. Such gene flow events during the MBA contributed toward shaping present-day Greek genomes.

Systematic benchmark of ancient DNA read mapping

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

Environmental palaeogenomic reconstruction of an Ice Age algal population

Palaeogenomics has greatly increased our knowledge of past evolutionary and ecological change, but has been restricted to the study of species that preserve either as or within fossils. Here we show the potential of shotgun metagenomics to reveal population genomic information for a taxon that does not preserve in the body fossil record, the algae Nannochloropsis. We shotgun sequenced two lake sediment samples dated to the Last Glacial Maximum and reconstructed full chloroplast and mitochondrial genomes to explore within-lake population genomic variation. This revealed two major haplogroups for each organellar genome, which could be assigned to known varieties of N. limnetica, although we show that at least three haplotypes were present using our minimum haplotype diversity estimation method. These approaches demonstrate the utility of lake sedimentary ancient DNA (sedaDNA) for population genomic analysis, thereby opening the door to environmental palaeogenomics, which will unlock the full potential of sedaDNA.

Lammers et al. use sedimentary ancient DNA to reconstruct palaeogenomes of Nannochloropsis. This study demonstrates the value of sedaDNA for palaeogenomic reconstructions and population genomic analysis.

Formation and repair of unavoidable, endogenous interstrand cross-links in cellular DNA

Genome integrity is essential for life and, as a result, DNA repair systems evolved to remove unavoidable DNA lesions from cellular DNA. Many forms of life possess the capacity to remove interstrand DNA cross-links (ICLs) from their genome but the identity of the naturally-occurring, endogenous substrates that drove the evolution and retention of these DNA repair systems across a wide range of life forms remains uncertain. In this review, we describe more than a dozen chemical processes by which endogenous ICLs plausibly can be introduced into cellular DNA. The majority involve DNA degradation processes that introduce aldehyde residues into the double helix or reactions of DNA with endogenous low molecular weight aldehyde metabolites. A smaller number of the cross-linking processes involve reactions of DNA radicals generated by oxidation.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

Filling Important Gaps in the Genomic History of Southwest Asia

Many crucial developments in human prehistory occurred in Southwest Asia, including the transition to agriculture as well as the emergence of writing and of the earliest civilization. Two new studies in this issue of Cell map the genetic composition of human groups inhabiting the region during the sixth and first millennia Before Common Era (BCE) and uncover periods of significant population turnover.

Cost-effective straightforward method for captured whole mitogenome sequencing of ancient DNA

Working with mitochondrial DNA from highly degraded samples is challenging. We present a whole mitogenome Illumina-based sequencing method suitable for highly degraded samples. The method makes use of double-stranded library preparation with hybridization-based target enrichment. The aim of the study was to implement a new user-friendly method for analysing many ancient DNA samples at low cost. The method combines the Swift 2S™ Turbo library preparation kit and xGen® panel for mitogenome enrichment. Swift allows to use low input of aDNA and own adapters and primers, handles inhibitors well, and has only two purification steps. xGen is straightforward to use and is able to leverage already pooled libraries. Given the ancient DNA is more challenging to work with, the protocol was developed with several improvements, especially multiplying DNA input in case of low concentration DNA extractions followed by AMPure® beads size selection and real-time pre-capture PCR monitoring in order to avoid cycle-optimization step. Nine out of eleven analysed samples successfully retrieved mitogenomes. Hence, our method provides an effective analysis of whole mtDNA, and has proven to be fast, cost-effective, straightforward, with utilisation in population-wide research of burial sites.

Using in silico predicted ancestral genomes to improve the efficiency of paleogenome reconstruction

Paleogenomics is the nascent discipline concerned with sequencing and analysis of genome-scale information from historic, ancient, and even extinct samples. While once inconceivable due to the challenges of DNA damage, contamination, and the technical limitations of PCR-based Sanger sequencing, following the dawn of the second-generation sequencing revolution, it has rapidly become a reality. However, a significant challenge facing ancient DNA studies on extinct species is the lack of closely related reference genomes against which to map the sequencing reads from ancient samples. Although bioinformatic efforts to improve the assemblies have focused mainly in mapping algorithms, in this article we explore the potential of an alternative approach, namely using reconstructed ancestral genome as reference for mapping DNA sequences of ancient samples. Specifically, we present a preliminary proof of concept for a general framework and demonstrate how under certain evolutionary divergence thresholds, considerable mapping improvements can be easily obtained.

Isolation, Library Preparation, and Bioinformatic Analysis of Historical and Ancient Plant DNA

The ability to sequence DNA retrieved from ancient and historical material plays a crucial role in reinforcing evolutionary and anthropological inference. While the focus of the field is largely on analyzing DNA from ancient hominids and other animals, we have also learned from plant ancient DNA (aDNA), in particular, about human farming practices, crop domestication, environment management, species invasion, and adaptation to various environmental conditions. In the following protocols, we outline best practices for plant aDNA isolation, preparation for sequencing, bioinformatic processing, and authentication. We describe the process all the way from processing of archaeological or historical plant material to characterizing and authenticating sequencing reads. In alternative protocols, we include modifications to this process that are tailored to strongly degraded DNA. Throughout, we stress the importance of precautionary measures to successfully analyze aDNA. Finally, we discuss the evolution of the archaeogenomics field and the development of new methods, which both shaped this protocol. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Isolation of aDNA Alternate Protocol 1: Isolation of ultra-short DNA (Dabney modification) Support Protocol 1: Preparation of PTB-based mix Support Protocol 2: Preparation of binding buffer Basic Protocol 2: Preparation of genomic libraries Alternate Protocol 2: Preparation of genomic libraries with uracil removal Basic Protocol 3: Bioinformatic processing and authentication of aDNA.

Genomics and the Evolutionary History of Equids

The equid family contains only one single extant genus, Equus, including seven living species grouped into horses on the one hand and zebras and asses on the other. In contrast, the equine fossil record shows that an extraordinarily richer diversity existed in the past and provides multiple examples of a highly dynamic evolution punctuated by several waves of explosive radiations and extinctions, cross-continental migrations, and local adaptations. In recent years, genomic technologies have provided new analytical solutions that have enhanced our understanding of equine evolution, including the species radiation within Equus; the extinction dynamics of several lineages; and the domestication history of two individual species, the horse and the donkey. Here, we provide an overview of these recent developments and suggest areas for further research.

Ancient DNA Studies in Pre-Columbian Mesoamerica

Mesoamerica is a historically and culturally defined geographic area comprising current central and south Mexico, Belize, Guatemala, El Salvador, and border regions of Honduras, western Nicaragua, and northwestern Costa Rica. The permanent settling of Mesoamerica was accompanied by the development of agriculture and pottery manufacturing (2500 BCE-150 CE), which led to the rise of several cultures connected by commerce and farming. Hence, Mesoamericans probably carried an invaluable genetic diversity partly lost during the Spanish conquest and the subsequent colonial period. Mesoamerican ancient DNA (aDNA) research has mainly focused on the study of mitochondrial DNA in the Basin of Mexico and the Yucatán Peninsula and its nearby territories, particularly during the Postclassic period (900-1519 CE). Despite limitations associated with the poor preservation of samples in tropical areas, recent methodological improvements pave the way for a deeper analysis of Mesoamerica. Here, we review how aDNA research has helped discern population dynamics patterns in the pre-Columbian Mesoamerican context, how it supports archaeological, linguistic, and anthropological conclusions, and finally, how it offers new working hypotheses.

The Recovery, Interpretation and Use of Ancient Pathogen Genomes

The ability to sequence genomes from ancient biological material has provided a rich source of information for evolutionary biology and engaged considerable public interest. Although most studies of ancient genomes have focused on vertebrates, particularly archaic humans, newer technologies allow the capture of microbial pathogens and microbiomes from ancient and historical human and non-human remains. This coming of age has been made possible by techniques that allow the preferential capture and amplification of discrete genomes from a background of predominantly host and environmental DNA. There are now near-complete ancient genome sequences for three pathogens of considerable historical interest - pre-modern bubonic plague (Yersinia pestis), smallpox (Variola virus) and cholera (Vibrio cholerae) - and for three equally important endemic human disease agents - Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy) and Treponema pallidum pallidum (syphilis). Genomic data from these pathogens have extended earlier work by paleopathologists. There have been efforts to sequence the genomes of additional ancient pathogens, with the potential to broaden our understanding of the infectious disease burden common to past populations from the Bronze Age to the early 20th century. In this review we describe the state-of-the-art of this rapidly developing field, highlight the contributions of ancient pathogen genomics to multidisciplinary endeavors and describe some of the limitations in resolving questions about the emergence and long-term evolution of pathogens.

Population genetics of wild Macaca fascicularis with low-coverage shotgun sequencing of museum specimens

OBJECTIVES

Long-tailed macaques (Macaca fascicularis) are widely distributed throughout the mainland and islands of Southeast Asia, making them a useful model for understanding the complex biogeographical history resulting from drastic changes in sea levels throughout the Pleistocene. Past studies based on mitochondrial genomes (mitogenomes) of long-tailed macaque museum specimens have traced their colonization patterns throughout the archipelago, but mitogenomes trace only the maternal history. Here, our objectives were to trace phylogeographic patterns of long-tailed macaques using low-coverage nuclear DNA (nDNA) data from museum specimens.

METHODS

We performed population genetic analyses and phylogenetic reconstruction on nuclear single nucleotide polymorphisms (SNPs) from shotgun sequencing of 75 long-tailed macaque museum specimens from localities throughout Southeast Asia.

RESULTS

We show that shotgun sequencing of museum specimens yields sufficient genome coverage (average ~1.7%) for reconstructing population relationships using SNP data. Contrary to expectations of divergent results between nuclear and mitochondrial genomes for a female philopatric species, phylogeographical patterns based on nuclear SNPs proved to be closely similar to those found using mitogenomes. In particular, population genetic analyses and phylogenetic reconstruction from the nDNA identify two major clades within M. fascicularis: Clade A includes all individuals from the mainland along with individuals from northern Sumatra, while Clade B consists of the remaining island-living individuals, including those from southern Sumatra.

CONCLUSIONS

Overall, we demonstrate that low-coverage sequencing of nDNA from museum specimens provides enough data for examining broad phylogeographic patterns, although greater genome coverage and sequencing depth would be needed to distinguish between very closely related populations, such as those throughout the Philippines.

The landscape of persistent human DNA viruses in femoral bone

The imprints left by persistent DNA viruses in the tissues can testify to the changes driving virus evolution as well as provide clues on the provenance of modern and ancient humans. However, the history hidden in skeletal remains is practically unknown, as only parvovirus B19 and hepatitis B virus DNA have been detected in hard tissues so far. Here, we investigated the DNA prevalences of 38 viruses in femoral bone of recently deceased individuals. To this end, we used quantitative PCRs and a custom viral targeted enrichment followed by next-generation sequencing. The data was analyzed with a tailor-made bioinformatics pipeline. Our findings revealed bone to be a much richer source of persistent DNA viruses than earlier perceived, discovering ten additional ones, including several members of the herpes- and polyomavirus families, as well as human papillomavirus 31 and torque teno virus. Remarkably, many of the viruses found have oncogenic potential and/or may reactivate in the elderly and immunosuppressed individuals. Thus, their persistence warrants careful evaluation of their clinical significance and impact on bone biology. Our findings open new frontiers for the study of virus evolution from ancient relics as well as provide new tools for the investigation of human skeletal remains in forensic and archaeological contexts.

Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification

Genetic tools are increasingly used to identify and discriminate between species. One key transition in this process was the recognition of the potential of the ca 658bp fragment of the organelle cytochrome c oxidase I (COI) as a barcode region, which revolutionized animal bioidentification and lead, among others, to the instigation of the Barcode of Life Database (BOLD), containing currently barcodes from >7.9 million specimens. Following this discovery, suggestions for other organellar regions and markers, and the primers with which to amplify them, have been continuously proposed. Most recently, the field has taken the leap from PCR-based generation of DNA references into shotgun sequencing-based "genome skimming" alternatives, with the ultimate goal of assembling organellar reference genomes. Unfortunately, in genome skimming approaches, much of the nuclear genome (as much as 99% of the sequence data) is discarded, which is not only wasteful, but can also limit the power of discrimination at, or below, the species level. Here, we advocate that the full shotgun sequence data can be used to assign an identity (that we term for convenience its "DNA-mark") for both voucher and query samples, without requiring any computationally intensive pretreatment (e.g. assembly) of reads. We argue that if reference databases are populated with such "DNA-marks," it will enable future DNA-based taxonomic identification to complement, or even replace PCR of barcodes with genome skimming, and we discuss how such methodology ultimately could enable identification to population, or even individual, level.

SPLATCHE3: simulation of serial genetic data under spatially explicit evolutionary scenarios including long-distance dispersal

SUMMARY

SPLATCHE3 simulates genetic data under a variety of spatially explicit evolutionary scenarios, extending previous versions of the framework. The new capabilities include long-distance migration, spatially and temporally heterogeneous short-scale migrations, alternative hybridization models, simulation of serial samples of genetic data and a large variety of DNA mutation models. These implementations have been applied independently to various studies, but grouped together in the current version.

AVAILABILITY AND IMPLEMENTATION

SPLATCHE3 is written in C++ and is freely available for non-commercial use from the website http://www.splatche.com/splatche3. It includes console versions for Linux, MacOs and Windows and a user-friendly GUI for Windows, as well as detailed documentation and ready-to-use examples.

Craniodental and Postcranial Characters of Non-Avian Dinosauria Often Imply Different Trees

Despite the increasing importance of molecular sequence data, morphology still makes an important contribution to resolving the phylogeny of many groups, and is the only source of data for most fossils. Most systematists sample morphological characters as broadly as possible on the principle of total evidence. However, it is not uncommon for sampling to be focused on particular aspects of anatomy, either because characters therein are believed to be more informative, or because preservation biases restrict what is available. Empirically, the optimal trees from partitions of morphological data sets often represent significantly different hypotheses of relationships. Previous work on hard-part versus soft-part characters across animal phyla revealed significant differences in about a half of sampled studies. Similarly, studies of the craniodental versus postcranial characters of vertebrates revealed significantly different trees in about one-third of cases, with the highest rates observed in non-avian dinosaurs. We test whether this is a generality here with a much larger sample of 81 published data matrices across all major dinosaur groups. Using the incongruence length difference test and two variants of the incongruence relationship difference test, we found significant incongruence in about 50% of cases. Incongruence is not uniformly distributed across major dinosaur clades, being highest (63%) in Theropoda and lowest (25%) in Thyreophora. As in previous studies, our partition tests show some sensitivity to matrix dimensions and the amount and distribution of missing entries. Levels of homoplasy and retained synapomorphy are similar between partitions, such that incongruence must partly reflect differences in patterns of homoplasy between partitions, which may itself be a function of modularity and mosaic evolution. Finally, we implement new tests to determine which partition yields trees most similar to those from the entire matrix. Despite no bias across dinosaurs overall, there are striking differences between major groups. The craniodental characters of Ornithischia and the postcranial characters of Saurischia yield trees most similar to the "total evidence" trees derived from the entire matrix. Trees from these same character partitions also tend to be most stratigraphically congruent: a mutual consilience suggesting that those partitions yield more accurate trees. [Dinosauria; homoplasy; partition homogeneity.].

Mining ancient microbiomes using selective enrichment of damaged DNA molecules

Background

The identification of bona fide microbial taxa in microbiomes derived from ancient and historical samples is complicated by the unavoidable mixture between DNA from ante- and post-mortem microbial colonizers. One possibility to distinguish between these sources of microbial DNA is querying for the presence of age-associated degradation patterns typical of ancient DNA (aDNA). The presence of uracils, resulting from cytosine deamination, has been detected ubiquitously in aDNA retrieved from diverse sources, and used as an authentication criterion. Here, we employ a library preparation method that separates molecules that carry uracils from those that do not for a set of samples that includes Neandertal remains, herbarium specimens and archaeological plant remains.

Results

We show that sequencing DNA libraries enriched in molecules carrying uracils effectively amplifies age associated degradation patterns in microbial mixtures of ancient and historical origin. This facilitates the discovery of authentic ancient microbial taxa in cases where degradation patterns are difficult to detect due to large sequence divergence in microbial mixtures. Additionally, the relative enrichment of taxa in the uracil enriched fraction can help to identify bona fide ancient microbial taxa that could be missed using a more targeted approach.

Conclusions

Our experiments show, that in addition to its use in enriching authentic endogenous DNA of organisms of interest, the selective enrichment of damaged DNA molecules can be a valuable tool in the discovery of ancient microbial taxa.

Museomics identifies genetic erosion in two butterfly species across the 20th century in Finland

Erosion of biodiversity generated by anthropogenic activities has been studied for decades and in many areas at the species level, using taxa monitoring. In contrast, genetic erosion within species has rarely been tracked, and is often studied by inferring past population dynamics from contemporaneous estimators. An alternative to such inferences is the direct examination of past genes, by analysing museum collection specimens. While providing direct access to genetic variation over time, historical DNA is usually not optimally preserved, and it is necessary to apply genotyping methods based on hybridization-capture to unravel past genetic variation. In this study, we apply such a method (i.e., HyRAD), to large time series of two butterfly species in Finland, and present a new bioinformatic pipeline, namely PopHyRAD, that standardizes and optimizes the analysis of HyRAD data at the within-species level. In the localities for which the data retrieved have sufficient power to accurately examine genetic dynamics through time, we show that genetic erosion has increased across the last 100 years, as revealed by signatures of allele extinctions and heterozygosity decreases, despite local variations. In one of the two butterflies (Erebia embla), isolation by distance also increased through time, revealing the effect of greater habitat fragmentation over time.

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.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

Quaternary DNA: A Multidisciplinary Research Field

The purpose of this Milankovitch review is to explain the significance of Quaternary DNA studies and the importance of the recent methodological advances that have enabled the study of late Quaternary remains in more detail, and the testing of new assumptions in evolutionary biology and phylogeography to reconstruct the past. The topic is wide, and this review is not intended to be an exhaustive account of all the aDNA work performed in the last three decades on late-Quaternary remains. Instead, it is a selection of relevant studies aimed at illustrating how aDNA has been used to reconstruct not only environments of the past, but also the history of many species including our own.

Helicobacter pylori in ancient human remains

The bacterium Helicobacter pylori (H. pylori) infects the stomachs of approximately 50% of all humans. With its universal occurrence, high infectivity and virulence properties it is considered as one of the most severe global burdens of modern humankind. It has accompanied humans for many thousands of years, and due to its high genetic variability and vertical transmission, its population genetics reflects the history of human migrations. However, especially complex demographic events such as the colonisation of Europe cannot be resolved with population genetic analysis of modern H. pylori strains alone. This is best exemplified with the reconstruction of the 5300-year-old H. pylori genome of the Iceman, a European Copper Age mummy. Our analysis provided precious insights into the ancestry and evolution of the pathogen and underlined the high complexity of ancient European population history. In this review we will provide an overview on the molecular analysis of H. pylori in mummified human remains that were done so far and we will outline methodological advancements in the field of ancient DNA research that support the reconstruction and authentication of ancient H. pylori genome sequences.

The Prevotella copri Complex Comprises Four Distinct Clades Underrepresented in Westernized Populations

Prevotella copri is a common human gut microbe that has been both positively and negatively associated with host health. In a cross-continent meta-analysis exploiting >6,500 metagenomes, we obtained >1,000 genomes and explored the genetic and population structure of P. copri. P. copri encompasses four distinct clades (>10% inter-clade genetic divergence) that we propose constitute the P. copri complex, and all clades were confirmed by isolate sequencing. These clades are nearly ubiquitous and co-present in non-Westernized populations. Genomic analysis showed substantial functional diversity in the complex with notable differences in carbohydrate metabolism, suggesting that multi-generational dietary modifications may be driving reduced prevalence in Westernized populations. Analysis of ancient metagenomes highlighted patterns of P. copri presence consistent with modern non-Westernized populations and a clade delineation time pre-dating human migratory waves out of Africa. These findings reveal that P. copri exhibits a high diversity that is underrepresented in Western-lifestyle populations.

Resolving a clinical tuberculosis outbreak using palaeogenomic genome reconstruction methodologies

This study describes the analysis of DNA from heat-killed (boilate) isolates of Mycobacterium tuberculosis from two UK outbreaks where DNA was of sub-optimal quality for the standard methodologies routinely used in microbial genomics. An Illumina library construction method developed for sequencing ancient DNA was successfully used to obtain whole genome sequences, allowing analysis of the outbreak by gene-by-gene MLST, SNP mapping and phylogenetic analysis. All cases were spoligotyped to the same Haarlem H1 sub-lineage. This is the first described application of ancient DNA library construction protocols to allow whole genome sequencing of a clinical tuberculosis outbreak. Using this method it is possible to obtain epidemiologically meaningful data even when DNA is of insufficient quality for standard methods.