Ancient DNA is typically highly degraded with appreciable cytosine deamination, and contamination with present-day DNA often complicates the identification of endogenous molecules. Together, these factors impede accurate assembly of the endogenous ancient mitochondrial genome. We present schmutzi, an iterative approach to jointly estimate present-day human contamination in ancient human DNA datasets and reconstruct the endogenous mitochondrial genome. By using sequence deamination patterns and fragment length distributions, schmutzi accurately reconstructs the endogenous mitochondrial genome sequence even when contamination exceeds 50 %. Given sufficient coverage, schmutzi also produces reliable estimates of contamination across a range of contamination rates. Availability: https://bioinf.eva.mpg.de/schmutzi/ license:GPLv3.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-015-0776-0) contains supplementary material, which is available to authorized users.

SummarySmallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1, 2, 3, 4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4, 5, 6, 7, 8, 9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4, 5, 6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18th and 19th centuries, concomitant with the development of modern vaccination.

Widespread circulation of SARS-CoV-2 in humans raises the theoretical risk of reverse zoonosis events with wildlife, reintroductions of SARS-CoV-2 into permissive nondomesticated animals. Here we report that North American deer mice (Peromyscus maniculatus) are susceptible to SARS-CoV-2 infection following intranasal exposure to a human isolate, resulting in viral replication in the upper and lower respiratory tract with little or no signs of disease. Further, shed infectious virus is detectable in nasal washes, oropharyngeal and rectal swabs, and viral RNA is detectable in feces and occasionally urine. We further show that deer mice are capable of transmitting SARS-CoV-2 to naïve deer mice through direct contact. The extent to which these observations may translate to wild deer mouse populations remains unclear, and the risk of reverse zoonosis and/or the potential for the establishment of Peromyscus rodents as a North American reservoir for SARS-CoV-2 remains unknown.

Sedimentary ancient DNA (sedaDNA) has been established as a viable biomolecular proxy for tracking taxon presence through time in a local environment, even in the total absence of surviving tissues. SedaDNA is thought to survive through mineral binding, facilitating long-term biomolecular preservation, but also challenging DNA isolation. Two common limitations in sedaDNA extraction are the carryover of other substances that inhibit enzymatic reactions, and the loss of authentic sedaDNA when attempting to reduce inhibitor co-elution. Here, we present a sedaDNA extraction procedure paired with targeted enrichment intended to maximize DNA recovery. Our procedure exhibits a 7.7–19.3x increase in on-target plant and animal sedaDNA compared to a commercial soil extraction kit, and a 1.2–59.9x increase compared to a metabarcoding approach. To illustrate the effectiveness of our cold spin extraction and PalaeoChip capture enrichment approach, we present results for the diachronic presence of plants and animals from Yukon permafrost samples dating to the Pleistocene-Holocene transition, and discuss new potential evidence for the late survival (~9700 years ago) of mammoth (Mammuthus sp.) and horse (Equus sp.) in the Klondike region of Yukon, Canada. This enrichment approach translates to a more taxonomically diverse dataset and improved on-target sequencing.

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

A widely accepted two-wave scenario of human settlement of Oceania involves the first out-of-Africa migration circa 50,000 years ago (ya), and the more recent Austronesian expansion, which reached the Bismarck Archipelago by 3,450 ya. Whereas earlier genetic studies provided evidence for extensive sex-biased admixture between the incoming and the indigenous populations, some archaeological, linguistic, and genetic evidence indicates a more complicated picture of settlement. To study regional variation in Oceania in more detail, we have compiled a genome-wide data set of 823 individuals from 72 populations (including 50 populations from Oceania) and over 620,000 autosomal single nucleotide polymorphisms (SNPs). We show that the initial dispersal of people from the Bismarck Archipelago into Remote Oceania occurred in a “leapfrog” fashion, completely by-passing the main chain of the Solomon Islands, and that the colonization of the Solomon Islands proceeded in a bidirectional manner. Our results also support a divergence between western and eastern Solomons, in agreement with the sharp linguistic divide known as the Tryon–Hackman line. We also report substantial post-Austronesian gene flow across the Solomons. In particular, Santa Cruz (in Remote Oceania) exhibits extraordinarily high levels of Papuan ancestry that cannot be explained by a simple bottleneck/founder event scenario. Finally, we use simulations to show that discrepancies between different methods for dating admixture likely reflect different sensitivities of the methods to multiple admixture events from the same (or similar) sources. Overall, this study points to the importance of fine-scale sampling to understand the complexities of human population history.

Archaeology, linguistics, and existing genetic studies indicate that Oceania was settled by two major waves of migration. The first migration took place approximately 40 thousand years ago and these migrants, Papuans, colonized much of Near Oceania. Approximately 3.5 thousand years ago, a second expansion of Austronesian-speakers arrived in Near Oceania and the descendants of these people spread to the far corners of the Pacific, colonizing Remote Oceania. To assess the female contribution of these two human expansions to modern populations and to investigate the potential impact of other migrations, we obtained 1,331 whole mitochondrial genome sequences from 34 populations spanning both Near and Remote Oceania. Our results quantify the magnitude of the Austronesian expansion and demonstrate the homogenizing effect of this expansion on almost all studied populations. With regards to Papuan influence, autochthonous haplogroups support the hypothesis of a long history in Near Oceania, with some lineages suggesting a time depth of 60 thousand years, and offer insight into historical interpopulation dynamics. Santa Cruz, a population located in Remote Oceania, is an anomaly with extreme frequencies of autochthonous haplogroups of Near Oceanian origin; simulations to investigate whether this might reflect a pre-Austronesian versus Austronesian settlement of the island failed to provide unequivocal support for either scenario.

No abstract