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

Mycobacterium leprae diversity and population dynamics in medieval Europe from novel ancient genomes

BACKGROUND

Hansen's disease (leprosy), widespread in medieval Europe, is today mainly prevalent in tropical and subtropical regions with around 200,000 new cases reported annually. Despite its long history and appearance in historical records, its origins and past dissemination patterns are still widely unknown. Applying ancient DNA approaches to its major causative agent, Mycobacterium leprae, can significantly improve our understanding of the disease's complex history. Previous studies have identified a high genetic continuity of the pathogen over the last 1500 years and the existence of at least four M. leprae lineages in some parts of Europe since the Early Medieval period.

RESULTS

Here, we reconstructed 19 ancient M. leprae genomes to further investigate M. leprae's genetic variation in Europe, with a dedicated focus on bacterial genomes from previously unstudied regions (Belarus, Iberia, Russia, Scotland), from multiple sites in a single region (Cambridgeshire, England), and from two Iberian leprosaria. Overall, our data confirm the existence of similar phylogeographic patterns across Europe, including high diversity in leprosaria. Further, we identified a new genotype in Belarus. By doubling the number of complete ancient M. leprae genomes, our results improve our knowledge of the past phylogeography of M. leprae and reveal a particularly high M. leprae diversity in European medieval leprosaria.

CONCLUSIONS

Our findings allow us to detect similar patterns of strain diversity across Europe with branch 3 as the most common branch and the leprosaria as centers for high diversity. The higher resolution of our phylogeny tree also refined our understanding of the interspecies transfer between red squirrels and humans pointing to a late antique/early medieval transmission. Furthermore, with our new estimates on the past population diversity of M. leprae, we gained first insights into the disease's global history in relation to major historic events such as the Roman expansion or the beginning of the regular transatlantic long distance trade. In summary, our findings highlight how studying ancient M. leprae genomes worldwide improves our understanding of leprosy's global history and can contribute to current models of M. leprae's worldwide dissemination, including interspecies transmissions.

The origin and legacy of the Etruscans through a 2000-year archeogenomic time transect

The origin, development, and legacy of the enigmatic Etruscan civilization from the central region of the Italian peninsula known as Etruria have been debated for centuries. Here we report a genomic time transect of 82 individuals spanning almost two millennia (800 BCE to 1000 CE) across Etruria and southern Italy. During the Iron Age, we detect a component of Indo-European–associated steppe ancestry and the lack of recent Anatolian-related admixture among the putative non–Indo-European–speaking Etruscans. Despite comprising diverse individuals of central European, northern African, and Near Eastern ancestry, the local gene pool is largely maintained across the first millennium BCE. This drastically changes during the Roman Imperial period where we report an abrupt population-wide shift to ~50% admixture with eastern Mediterranean ancestry. Last, we identify northern European components appearing in central Italy during the Early Middle Ages, which thus formed the genetic landscape of present-day Italian populations.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Population genomics of the critically endangered kākāpō

The kākāpō is a flightless parrot endemic to New Zealand. Once common in the archipelago, only 201 individuals remain today, most of them descending from an isolated island population. We report the first genome-wide analyses of the species, including a high-quality genome assembly for kākāpō, one of the first chromosome-level reference genomes sequenced by the Vertebrate Genomes Project (VGP). We also sequenced and analyzed 35 modern genomes from the sole surviving island population and 14 genomes from the extinct mainland population. While theory suggests that such a small population is likely to have accumulated deleterious mutations through genetic drift, our analyses on the impact of the long-term small population size in kākāpō indicate that present-day island kākāpō have a reduced number of harmful mutations compared to mainland individuals. We hypothesize that this reduced mutational load is due to the island population having been subjected to a combination of genetic drift and purging of deleterious mutations, through increased inbreeding and purifying selection, since its isolation from the mainland ∼10,000 years ago. Our results provide evidence that small populations can survive even when isolated for hundreds of generations. This work provides key insights into kākāpō breeding and recovery and more generally into the application of genetic tools in conservation efforts for endangered species.

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.

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.

Exceptional ancient DNA preservation and fibre remains of a Sasanian saltmine sheep mummy in Chehrābād, Iran

Mummified remains have long attracted interest as a potential source of ancient DNA. However, mummification is a rare process that requires an anhydrous environment to rapidly dehydrate and preserve tissue before complete decomposition occurs. We present the whole-genome sequences (3.94 X) of an approximately 1600-year-old naturally mummified sheep recovered from Chehrābād, a salt mine in northwestern Iran. Comparative analyses of published ancient sequences revealed the remarkable DNA integrity of this mummy. Hallmarks of postmortem damage, fragmentation and hydrolytic deamination are substantially reduced, likely owing to the high salinity of this taphonomic environment. Metagenomic analyses reflect the profound influence of high-salt content on decomposition; its microbial profile is predominated by halophilic archaea and bacteria, possibly contributing to the remarkable preservation of the sample. Applying population genomic analyses, we find clustering of this sheep with Southwest Asian modern breeds, suggesting ancestry continuity. Genotyping of a locus influencing the woolly phenotype showed the presence of an ancestral 'hairy' allele, consistent with hair fibre imaging. This, along with derived alleles associated with the fat-tail phenotype, provides genetic evidence that Sasanian-period Iranians maintained specialized sheep flocks for different uses, with the 'hairy', 'fat-tailed'-genotyped sheep likely kept by the rural community of Chehrābād's miners.

Recovery and analysis of ancient beetle DNA from subfossil packrat middens using high-throughput sequencing

The study of ancient DNA is revolutionizing our understanding of paleo-ecology and the evolutionary history of species. Insects are essential components in many ecosystems and constitute the most diverse group of animals. Yet they are largely neglected in ancient DNA studies. We report the results of the first targeted investigation of insect ancient DNA to positively identify subfossil insects to species, which includes the recovery of endogenous content from samples as old as ~ 34,355 ybp. Potential inhibitors currently limiting widespread research on insect ancient DNA are discussed, including the lack of closely related genomic reference sequences (decreased mapping efficiency) and the need for more extensive collaborations with insect taxonomists. The advantages of insect-based studies are also highlighted, especially in the context of understanding past climate change. In this regard, insect remains from ancient packrat middens are a rich and largely uninvestigated resource for exploring paleo-ecology and species dynamics over time.

Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains

The Aceramic Neolithic (∼9600 to 7000 cal BC) period in the Zagros Mountains, western Iran, provides some of the earliest archaeological evidence of goat (Capra hircus) management and husbandry by circa 8200 cal BC, with detectable morphological change appearing ∼1,000 y later. To examine the genomic imprint of initial management and its implications for the goat domestication process, we analyzed 14 novel nuclear genomes (mean coverage 1.13X) and 32 mitochondrial (mtDNA) genomes (mean coverage 143X) from two such sites, Ganj Dareh and Tepe Abdul Hosein. These genomes show two distinct clusters: those with domestic affinity and a minority group with stronger wild affinity, indicating that managed goats were genetically distinct from wild goats at this early horizon. This genetic duality, the presence of long runs of homozygosity, shared ancestry with later Neolithic populations, a sex bias in archaeozoological remains, and demographic profiles from across all layers of Ganj Dareh support management of genetically domestic goat by circa 8200 cal BC, and represent the oldest to-this-date reported livestock genomes. In these sites a combination of high autosomal and mtDNA diversity, contrasting limited Y chromosomal lineage diversity, an absence of reported selection signatures for pigmentation, and the wild morphology of bone remains illustrates domestication as an extended process lacking a strong initial bottleneck, beginning with spatial control, demographic manipulation via biased male culling, captive breeding, and subsequently phenotypic and genomic selection.

Progress in forensic bone DNA analysis: Lessons learned from ancient DNA

Research on ancient and forensic DNA is related in many ways, and the two fields must deal with similar obstacles. Therefore, communication between these two communities has the potential to improve results in both research fields. Here, we present the insights gained in the ancient DNA community with regard to analyzing DNA from aged skeletal material and the potential use of the developed protocols in forensic work. We discuss the various steps, from choosing samples for DNA extraction to deciding between classical PCR amplification and massively parallel sequencing approaches. Based on the progress made in ancient DNA analyses combined with the requirements of forensic work, we suggest that there is substantial potential for incorporating ancient DNA approaches into forensic protocols, a process that has already begun to a considerable extent. However, taking full advantage of the experiences gained from ancient DNA work will require comparative studies by the forensic DNA community to tailor the methods developed for ancient samples to the specific needs of forensic studies and case work. If successful, in our view, the benefits for both communities would be considerable.

Successful application of ancient DNA extraction and library construction protocols to museum wet collection specimens

Millions of scientific specimens are housed in museum collections, a large part of which are fluid preserved. The use of formaldehyde as fixative and subsequent storage in ethanol is especially common in ichthyology and herpetology. This type of preservation damages DNA and reduces the chance of successful retrieval of genetic data. We applied ancient DNA extraction and single stranded library construction protocols to a variety of vertebrate samples obtained from wet collections and of different ages. Our results show that almost all samples tested yielded endogenous DNA. Archival DNA extraction was successful across different tissue types as well as using small amounts of tissue. Conversion of archival DNA fragments into single-stranded libraries resulted in usable data even for samples with initially undetectable DNA amounts. Subsequent target capture approaches for mitochondrial DNA using homemade baits on a subset of 30 samples resulted in almost complete mitochondrial genome sequences in several instances. Thus, application of ancient DNA methodology makes wet collection specimens, including type material as well as rare, old or extinct species, accessible for genetic and genomic analyses. Our results, accompanied by detailed step-by-step protocols, are a large step forward to open the DNA archive of museum wet collections for scientific studies.

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

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

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.

Reconstruction of ancient microbial genomes from the human gut

Loss of gut microbial diversity^1–6 in industrial populations is associated with chronic diseases⁷, underscoring the importance of studying our ancestral gut microbiome. However, relatively little is known about the composition of pre-industrial gut microbiomes. Here we performed a large-scale de novo assembly of microbial genomes from palaeofaeces. From eight authenticated human palaeofaeces samples (1,000–2,000 years old) with well-preserved DNA from southwestern USA and Mexico, we reconstructed 498 medium- and high-quality microbial genomes. Among the 181 genomes with the strongest evidence of being ancient and of human gut origin, 39% represent previously undescribed species-level genome bins. Tip dating suggests an approximate diversification timeline for the key human symbiont Methanobrevibacter smithii. In comparison to 789 present-day human gut microbiome samples from eight countries, the palaeofaeces samples are more similar to non-industrialized than industrialized human gut microbiomes. Functional profiling of the palaeofaeces samples reveals a markedly lower abundance of antibiotic-resistance and mucin-degrading genes, as well as enrichment of mobile genetic elements relative to industrial gut microbiomes. This study facilitates the discovery and characterization of previously undescribed gut microorganisms from ancient microbiomes and the investigation of the evolutionary history of the human gut microbiota through genome reconstruction from palaeofaeces.

Ancient microbiomes from palaeofaeces are more similar to non-industrialized than industrialized human gut microbiomes regardless of geography, but 39% of their de novo reconstructed genomes represent previously undescribed microbial species.

The evolution and changing ecology of the African hominid oral microbiome

The microbiome plays key roles in human health, but little is known about its evolution. We investigate the evolutionary history of the African hominid oral microbiome by analyzing dental biofilms of humans and Neanderthals spanning the past 100,000 years and comparing them with those of chimpanzees, gorillas, and howler monkeys. We identify 10 core bacterial genera that have been maintained within the human lineage and play key biofilm structural roles. However, many remain understudied and unnamed. We find major taxonomic and functional differences between the oral microbiomes of Homo and chimpanzees but a high degree of similarity between Neanderthals and modern humans, including an apparent Homo-specific acquisition of starch digestion capability in oral streptococci, suggesting microbial coadaptation with host diet.

The oral microbiome plays key roles in human biology, health, and disease, but little is known about the global diversity, variation, or evolution of this microbial community. To better understand the evolution and changing ecology of the human oral microbiome, we analyzed 124 dental biofilm metagenomes from humans, including Neanderthals and Late Pleistocene to present-day modern humans, chimpanzees, and gorillas, as well as New World howler monkeys for comparison. We find that a core microbiome of primarily biofilm structural taxa has been maintained throughout African hominid evolution, and these microbial groups are also shared with howler monkeys, suggesting that they have been important oral members since before the catarrhine–platyrrhine split ca. 40 Mya. However, community structure and individual microbial phylogenies do not closely reflect host relationships, and the dental biofilms of Homo and chimpanzees are distinguished by major taxonomic and functional differences. Reconstructing oral metagenomes from up to 100 thousand years ago, we show that the microbial profiles of both Neanderthals and modern humans are highly similar, sharing functional adaptations in nutrient metabolism. These include an apparent Homo-specific acquisition of salivary amylase-binding capability by oral streptococci, suggesting microbial coadaptation with host diet. We additionally find evidence of shared genetic diversity in the oral bacteria of Neanderthal and Upper Paleolithic modern humans that is not observed in later modern human populations. Differences in the oral microbiomes of African hominids provide insights into human evolution, the ancestral state of the human microbiome, and a temporal framework for understanding microbial health and disease.

Systematic benchmark of ancient DNA read mapping

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

Million-year-old DNA sheds light on the genomic history of mammoths

Temporal genomic data hold great potential for studying evolutionary processes, including speciation. However, sampling across speciation events would in many cases require genomic time series that stretch well into the Early Pleistocene (>1 million years). Although theoretical models suggest that DNA should survive on this timescale¹, the oldest genomic data recovered so far is from a 560-780 ka old horse specimen². Here we report the recovery of genome-wide data from three Early and Middle Pleistocene mammoth specimens, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these gave rise to the woolly mammoth, whereas the other represents a previously unrecognised lineage that was ancestral to the first mammoths to colonise North America. Our analyses reveal that the North American Columbian mammoth traces its ancestry to a Middle Pleistocene hybridisation between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were present already a million years ago. These findings highlight the potential of deep time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

Later Stone Age human hair from Vaalkrans Shelter, Cape Floristic Region of South Africa, reveals genetic affinity to Khoe groups

Previous studies show that the indigenous people of the southern Cape of South Africa were dramatically impacted by the arrival of European colonists starting ~400 years ago and their descendants are today mixed with Europeans and Asians. To gain insight on the occupants of the Vaalkrans Shelter located at the southernmost tip of Africa, we investigated the genetic make-up of an individual who lived there about 200 years ago. We further contextualize the genetic ancestry of this individual among prehistoric and current groups. From a hair sample excavated at the shelter, which was indirectly dated to about 200 years old, we sequenced the genome (1.01 times coverage) of a Later Stone Age individual. We analyzed the Vaalkrans genome together with genetic data from 10 ancient (pre-colonial) individuals from southern Africa spanning the last 2000 years. We show that the individual from Vaalkrans was a man who traced ~80% of his ancestry to local southern San hunter-gatherers and ~20% to a mixed East African-Eurasian source. This genetic make-up is similar to modern-day Khoekhoe individuals from the Northern Cape Province (South Africa) and Namibia, but in the southern Cape, the Vaalkrans man's descendants have likely been assimilated into mixed-ancestry "Coloured" groups. The Vaalkrans man's genome reveals that Khoekhoe pastoralist groups/individuals lived in the southern Cape as late as 200 years ago, without mixing with non-African colonists or Bantu-speaking farmers. Our findings are also consistent with the model of a Holocene pastoralist migration, originating in Eastern Africa, shaping the genomic landscape of historic and current southern African populations.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10^-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

A genetic history of the pre-contact Caribbean

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

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

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

Moose genomes reveal past glacial demography and the origin of modern lineages

Background

Numerous megafauna species from northern latitudes went extinct during the Pleistocene/Holocene transition as a result of climate-induced habitat changes. However, several ungulate species managed to successfully track their habitats during this period to eventually flourish and recolonise the holarctic regions. So far, the genomic impacts of these climate fluctuations on ungulates from high latitudes have been little explored. Here, we assemble a de-novo genome for the European moose (Alces alces) and analyse it together with re-sequenced nuclear genomes and ancient and modern mitogenomes from across the moose range in Eurasia and North America.

Results

We found that moose demographic history was greatly influenced by glacial cycles, with demographic responses to the Pleistocene/Holocene transition similar to other temperate ungulates. Our results further support that modern moose lineages trace their origin back to populations that inhabited distinct glacial refugia during the Last Glacial Maximum (LGM). Finally, we found that present day moose in Europe and North America show low to moderate inbreeding levels resulting from post-glacial bottlenecks and founder effects, but no evidence for recent inbreeding resulting from human-induced population declines.

Conclusions

Taken together, our results highlight the dynamic recent evolutionary history of the moose and provide an important resource for further genomic studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07208-3.

Competitive mapping allows for the identification and exclusion of human DNA contamination in ancient faunal genomic datasets

BACKGROUND

After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data.

RESULTS

Here we investigate whether human contaminating DNA can be found in ancient faunal sequencing datasets. We identify variable levels of human contamination, which persists even after the sequence reads have been mapped to the faunal reference genomes. This contamination has the potential to affect a range of downstream analyses.

CONCLUSIONS

We propose a fast and simple method, based on competitive mapping, which allows identifying and removing human contamination from ancient faunal DNA datasets with limited losses of true ancient data. This method could represent an important tool for the ancient DNA field.

The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics

The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.

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.

Were domestic camelids present on the prehispanic South American agricultural frontier? An ancient DNA study

The southern boundary of prehispanic farming in South America occurs in central Mendoza Province, Argentina at approximately 34 degrees south latitude. Archaeological evidence of farming includes the recovery of macrobotanical remains of cultigens and isotopic chemistry of human bone. Since the 1990s, archaeologists have also hypothesized that the llama (Lama glama), a domesticated South American camelid, was also herded near the southern boundary of prehispanic farming. The remains of a wild congeneric camelid, the guanaco (Lama guanicoe), however, are common in archaeological sites throughout Mendoza Province. It is difficult to distinguish bones of the domestic llama from wild guanaco in terms of osteological morphology, and therefore, claims that llama were in geographic areas where guanaco were also present based on osteometric analysis alone remain equivocal. A recent study, for example, claimed that twenty-five percent of the camelid remains from the high elevation Andes site of Laguna del Diamante S4 were identified based on osteometric evidence as domestic llama, but guanaco are also a likely candidate since the two species overlap in size. We test the hypothesis that domesticated camelids occurred in prehispanic, southern Mendoza through analysis of ancient DNA. We generated whole mitochondrial genome datasets from 41 samples from southern Mendoza late Holocene archaeological sites, located between 450 and 3400 meters above sea level (masl). All camelid samples from those sites were identified as guanaco; thus, we have no evidence to support the hypothesis that the domestic llama occurred in prehispanic southern Mendoza.

Origins and genetic legacy of prehistoric dogs

Dogs were the first domestic animal, but little is known about their population history and to what extent it was linked to humans. We sequenced 27 ancient dog genomes and found that all dogs share a common ancestry distinct from present-day wolves, with limited gene flow from wolves since domestication but substantial dog-to-wolf gene flow. By 11,000 years ago, at least five major ancestry lineages had diversified, demonstrating a deep genetic history of dogs during the Paleolithic. Coanalysis with human genomes reveals aspects of dog population history that mirror humans, including Levant-related ancestry in Africa and early agricultural Europe. Other aspects differ, including the impacts of steppe pastoralist expansions in West and East Eurasia and a near-complete turnover of Neolithic European dog ancestry.

A systematic investigation of human DNA preservation in medieval skeletons

Ancient DNA (aDNA) analyses necessitate the destructive sampling of archaeological material. Currently, the cochlea, part of the osseous inner ear located inside the petrous pyramid, is the most sought after skeletal element for molecular analyses of ancient humans as it has been shown to yield high amounts of endogenous DNA. However, destructive sampling of the petrous pyramid may not always be possible, particularly in cases where preservation of skeletal morphology is of top priority. To investigate alternatives, we present a survey of human aDNA preservation for each of ten skeletal elements in a skeletal collection from Medieval Germany. Through comparison of human DNA content and quality we confirm best performance of the petrous pyramid and identify seven additional sampling locations across four skeletal elements that yield adequate aDNA for most applications in human palaeogenetics. Our study provides a better perspective on DNA preservation across the human skeleton and takes a further step toward the more responsible use of ancient materials in human aDNA studies.

Alteration of methylation status in archival DNA samples: A qualitative assessment by methylation specific polymerase chain reaction

DNA methylation is an epigenetic mechanism that regulates gene expression, which also facilitates genomic imprinting. Genomic imprinting is responsible for differential expression of genes based on parent of origin. Altered methylation of parental alleles results in imprinting disorders diagnosed by methylation specific polymerase chain reaction (MS-PCR) technique. With increasing evidence of genes under epigenetic influence, methylation studies are extensively performed on archival samples. To evaluate effect of storage and storage conditions on DNA methylation, a systematic MS-PCR based analysis was planned on an imprinted gene, SNRPN, located on chromosome 15q11.2. It was assessed by MS-PCR on fresh, 4 -20, and -80°C stored DNA samples for different time periods for systematic evaluation of methylation status. Technical factors like type of sample processing, method of DNA isolation, primer region polymorphism, sample heterogeneity were also evaluated. DNA methylation was observed to be altered for SNRPN gene after storage at -80°C from 2 months onwards. Long-term storage of DNA at -80°C results in altered DNA methylation status. This may lead to false MS-PCR diagnosis of imprinting disorders. Our proof of concept study should be followed with quantitative validation since the findings have critical implications in the present era of biobanking.

Damage patterns observed in mtDNA control region MPS data for a range of template concentrations and when using different amplification approaches

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows practitioners the ability to fully resolve heteroplasmic sites. In forensic DNA analysis, identifying heteroplasmy (a naturally occurring mixture of two mtDNA profiles) can provide additional mtDNA profile information which can lead to an increase in the discrimination potential of an mtDNA match between an evidentiary sample and reference source. Forensic samples such as hair and skeletal remains, especially older, more compromised samples, can often exhibit DNA damage. Because both damage and heteroplasmy can manifest as a mixture of two nucleotides, it is important to differentiate between the two conditions when interpreting mtDNA MPS data. In this study, DNA damage was applied under controlled conditions to samples containing a range of template concentrations, including some with identified heteroplasmy. Damage was applied via storage in water at room temperature on samples diluted before or after storage to mimic low template scenarios. Damage was assessed with respect to the following areas: mtDNA quantification and degradation ratios, MPS read depth, MPS profile results, overall damage rates, and the interpretation of heteroplasmy. Datasets were generated to assess and compare two different amplification and library preparation strategies: the Promega PowerSeq™ CRM Nested System kit and a 1.16 kb target amplicon of the entire mtDNA control region followed by a Nextera® XT library preparation. The results of this study provide an evaluation of the Promega 10-plex MPS procedure as an improved process to mitigate the impact of mtDNA damage on low template samples. Some of the negative effects of damage observed in this study were a decrease in mtDNA yield by 20-30% and lower quality MPS sequencing results. These effects were observed more frequently when samples were diluted prior to inducing damage, illustrating that low template samples are more susceptible to damage. The findings of this study will assist forensic laboratories in differentiating between damage and heteroplasmy, which is essential when developing robust mtDNA MPS interpretation guidelines such as setting appropriate reporting thresholds.

Removing reference bias and improving indel calling in ancient DNA data analysis by mapping to a sequence variation graph

BACKGROUND

During the last decade, the analysis of ancient DNA (aDNA) sequence has become a powerful tool for the study of past human populations. However, the degraded nature of aDNA means that aDNA molecules are short and frequently mutated by post-mortem chemical modifications. These features decrease read mapping accuracy and increase reference bias, in which reads containing non-reference alleles are less likely to be mapped than those containing reference alleles. Alternative approaches have been developed to replace the linear reference with a variation graph which includes known alternative variants at each genetic locus. Here, we evaluate the use of variation graph software vg to avoid reference bias for aDNA and compare with existing methods.

RESULTS

We use vg to align simulated and real aDNA samples to a variation graph containing 1000 Genome Project variants and compare with the same data aligned with bwa to the human linear reference genome. Using vg leads to a balanced allelic representation at polymorphic sites, effectively removing reference bias, and more sensitive variant detection in comparison with bwa, especially for insertions and deletions (indels). Alternative approaches that use relaxed bwa parameter settings or filter bwa alignments can also reduce bias but can have lower sensitivity than vg, particularly for indels.

CONCLUSIONS

Our findings demonstrate that aligning aDNA sequences to variation graphs effectively mitigates the impact of reference bias when analyzing aDNA, while retaining mapping sensitivity and allowing detection of variation, in particular indel variation, that was previously missed.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.

Digging into the Genomic Past of Swiss Honey Bees by Whole-Genome Sequencing Museum Specimens

Historical specimens in museum collections provide opportunities to gain insights into the genomic past. For the Western honey bee, Apis mellifera L., this is particularly important because its populations are currently under threat worldwide and have experienced many changes in management and environment over the last century. Using Swiss Apis mellifera mellifera as a case study, our research provides important insights into the genetic diversity of native honey bees prior to the industrial-scale introductions and trade of non-native stocks during the 20th century—the onset of intensive commercial breeding and the decline of wild honey bees following the arrival of Varroa destructor. We sequenced whole-genomes of 22 honey bees from the Natural History Museum in Bern collected in Switzerland, including the oldest A. mellifera sample ever sequenced. We identify both, a historic and a recent migrant, natural or human-mediated, which corroborates with the population history of honey bees in Switzerland. Contrary to what we expected, we find no evidence for a significant genetic bottleneck in Swiss honey bees, and find that genetic diversity is not only maintained, but even slightly increased, most probably due to modern apicultural practices. Finally, we identify signals of selection between historic and modern honey bee populations associated with genes enriched in functions linked to xenobiotics, suggesting a possible selective pressure from the increasing use and diversity of chemicals used in agriculture and apiculture over the last century.

2000-year-old pathogen genomes reconstructed from metagenomic analysis of Egyptian mummified individuals

BACKGROUND

Recent advances in sequencing have facilitated large-scale analyses of the metagenomic composition of different samples, including the environmental microbiome of air, water, and soil, as well as the microbiome of living humans and other animals. Analyses of the microbiome of ancient human samples may provide insights into human health and disease, as well as pathogen evolution, but the field is still in its very early stages and considered highly challenging.

RESULTS

The metagenomic and pathogen content of Egyptian mummified individuals from different time periods was investigated via genetic analysis of the microbial composition of various tissues. The analysis of the dental calculus' microbiome identified Red Complex bacteria, which are correlated with periodontal diseases. From bone and soft tissue, genomes of two ancient pathogens, a 2200-year-old Mycobacterium leprae strain and a 2000-year-old human hepatitis B virus, were successfully reconstructed.

CONCLUSIONS

The results show the reliability of metagenomic studies on Egyptian mummified individuals and the potential to use them as a source for the extraction of ancient pathogen DNA.

Ancient Bacterial Genomes Reveal a High Diversity of Treponema pallidum Strains in Early Modern Europe

Syphilis is a globally re-emerging disease, which has marked European history with a devastating epidemic at the end of the 15^th century. Together with non-venereal treponemal diseases, like bejel and yaws, which are found today in subtropical and tropical regions, it currently poses a substantial health threat worldwide. The origins and spread of treponemal diseases remain unresolved, including syphilis' potential introduction into Europe from the Americas. Here, we present the first genetic data from archaeological human remains reflecting a high diversity of Treponema pallidum in early modern Europe. Our study demonstrates that a variety of strains related to both venereal syphilis and yaws-causing T. pallidum subspecies were already present in Northern Europe in the early modern period. We also discovered a previously unknown T. pallidum lineage recovered as a sister group to yaws- and bejel-causing lineages. These findings imply a more complex pattern of geographical distribution and etiology of early treponemal epidemics than previously understood.

Pre-extinction Demographic Stability and Genomic Signatures of Adaptation in the Woolly Rhinoceros

Ancient DNA has significantly improved our understanding of the evolution and population history of extinct megafauna. However, few studies have used complete ancient genomes to examine species responses to climate change prior to extinction. The woolly rhinoceros (Coelodonta antiquitatis) was a cold-adapted megaherbivore widely distributed across northern Eurasia during the Late Pleistocene and became extinct approximately 14 thousand years before present (ka BP). While humans and climate change have been proposed as potential causes of extinction [1-3], knowledge is limited on how the woolly rhinoceros was impacted by human arrival and climatic fluctuations [2]. Here, we use one complete nuclear genome and 14 mitogenomes to investigate the demographic history of woolly rhinoceros leading up to its extinction. Unlike other northern megafauna, the effective population size of woolly rhinoceros likely increased at 29.7 ka BP and subsequently remained stable until close to the species' extinction. Analysis of the nuclear genome from a ∼18.5-ka-old specimen did not indicate any increased inbreeding or reduced genetic diversity, suggesting that the population size remained steady for more than 13 ka following the arrival of humans [4]. The population contraction leading to extinction of the woolly rhinoceros may have thus been sudden and mostly driven by rapid warming in the Bølling-Allerød interstadial. Furthermore, we identify woolly rhinoceros-specific adaptations to arctic climate, similar to those of the woolly mammoth. This study highlights how species respond differently to climatic fluctuations and further illustrates the potential of palaeogenomics to study the evolutionary history of extinct species.

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

BACKGROUND

Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679).

RESULTS

A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent.

CONCLUSIONS

The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.

Present-Day DNA Contamination in Ancient DNA Datasets

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

Robbery in progress: Historical museum collections bring to light a mitochondrial capture within a bird species widespread across southern Australia, the Copperback Quail-thrush Cinclosoma clarum

We surveyed mitochondrial, autosomal, and Z chromosome diversity within and between the Copperback Quail-thrush Cinclosoma clarum and Chestnut Quail-thrush C. castanotum, which together span the arid and semi-arid zones of southern Australia, and primarily from specimens held in museum collections. We affirm the recent taxonomic separation of the two species and then focus on diversity within the more widespread of the two species, C. clarum. To guide further study of the system and what it offers to understanding the genomics of the differentiation and speciation processes, we develop and present a hypothesis to explain mitonuclear discordance that emerged in ourdata. Following a period of historical allopatry, secondary contact has resulted in an eastern mitochondrial genome replacing the western mitochondrial genome in western populations. This is predicted under a population-level invasion in the opposite direction, that of the western population invading the range of the eastern one. Mitochondrial captures can be driven by neutral, demographic processes, or adaptive mechanisms, and we favor the hypothesized capture being driven by neutral means. We cannot fully reject the adaptive process but suggest how these alternatives may be further tested. We acknowledge an alternative hypothesis, which finds some support in phenotypic data published elsewhere, namely that outcomes of secondary contact have been more complex than our current genomic data suggest. Discriminating and reconciling these two alternative hypotheses, which may not be mutually exclusive, could be tested with closer sampling at levels of population, individual, and nucleotide than has so far been possible. This would be further aided by knowledge of the genetic basis to phenotypic variation described elsewhere.

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.

Harnessing epigenetics to study human evolution

Recent advances in ancient DNA extraction and high-throughput sequencing technologies enabled the high-quality sequencing of archaic genomes, including the Neanderthal and the Denisovan. While comparisons with modern humans revealed both archaic-specific and human-specific sequence changes, in the absence of gene expression information, understanding the functional implications of such genetic variations remains a major challenge. To study gene regulation in archaic humans, epigenetic research comes to our aid. DNA methylation, which is highly correlated with transcription, can be directly measured in modern samples, as well as reconstructed in ancient samples. This puts DNA methylation as a natural basis for comparative epigenetics between modern humans, archaic humans and nonhuman primates.

A treponemal genome from an historic plague victim supports a recent emergence of yaws and its presence in 15th century Europe

Developments in techniques for identification of pathogen DNA in archaeological samples can expand our resolution of disease detection. Our application of a non-targeted molecular screening tool for the parallel detection of pathogens in historical plague victims from post-medieval Lithuania revealed the presence of more than one active disease in one individual. In addition to Yersinia pestis, we detected and genomically characterized a septic infection of Treponema pallidum pertenue, a subtype of the treponemal disease family recognised as the cause of the tropical disease yaws. Our finding in northern Europe of a disease that is currently restricted to equatorial regions is interpreted within an historical framework of intercontinental trade and potential disease movements. Through this we offer an alternative hypothesis for the history and evolution of the treponemal diseases, and posit that yaws be considered an important contributor to the sudden epidemic of late 15^th century Europe that is widely ascribed to syphilis.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Advances in the molecular detection of tuberculosis in pre-contact Andean South America

Andean paleopathological research has significantly enhanced knowledge about the geographical distribution and evolution of tuberculosis (TB) in pre-Columbian South America. In this paper, we review the history and progress of research on ancient tuberculosis (TB) in the Andean region, focusing on the strengths and limitations of current approaches for the molecular detection of ancient pathogens, with special attention to TB. As a case study, we describe a molecular screening approach for the detection of ancient Mycobacterium tuberculosis in individuals from Late Intermediate Period (1000-1400 CE) contexts at the site of Huari, Peru. We evaluate 34 commingled human vertebrae and combine morphological assessments of pathology with high throughput sequencing and a non-selective approach to ancient pathogen DNA screening. Our method enabled the simultaneous detection of ancient M. tuberculosis DNA and an evaluation of the environmental microbial composition of each sample. Our results show that despite the dominance of environmental DNA, molecular signatures of M. tuberculosis were identified in eight vertebrae, six of which had no observable skeletal pathology classically associated tuberculosis infection. This screening approach will assist in the identification of candidate samples for downstream genomic analyses. The method permits higher resolution disease identification in cases where pathology may be absent, or where the archaeological context may necessitate a broad differential diagnosis based on morphology alone.

A Paleogenomic Reconstruction of the Deep Population History of the Andes

There are many unanswered questions about the population history of the Central and South Central Andes, particularly regarding the impact of large-scale societies, such as the Moche, Wari, Tiwanaku, and Inca. We assembled genome-wide data on 89 individuals dating from ∼9,000-500 years ago (BP), with a particular focus on the period of the rise and fall of state societies. Today's genetic structure began to develop by 5,800 BP, followed by bi-directional gene flow between the North and South Highlands, and between the Highlands and Coast. We detect minimal admixture among neighboring groups between ∼2,000-500 BP, although we do detect cosmopolitanism (people of diverse ancestries living side-by-side) in the heartlands of the Tiwanaku and Inca polities. We also highlight cases of long-range mobility connecting the Andes to Argentina and the Northwest Andes to the Amazon Basin. VIDEO ABSTRACT.

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

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

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

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