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.