From AD 1347 to AD 1353, the Black Death killed tens of millions of people in Europe, leaving misery and devastation in its wake, with successive epidemics ravaging the continent until the 18th century. The etiology of this disease has remained highly controversial, ranging from claims based on genetics and the historical descriptions of symptoms that it was caused by Yersinia pestis to conclusions that it must have been caused by other pathogens. It has also been disputed whether plague had the same etiology in northern and southern Europe. Here we identified DNA and protein signatures specific for Y. pestis in human skeletons from mass graves in northern, central and southern Europe that were associated archaeologically with the Black Death and subsequent resurgences. We confirm that Y. pestis caused the Black Death and later epidemics on the entire European continent over the course of four centuries. Furthermore, on the basis of 17 single nucleotide polymorphisms plus the absence of a deletion in glpD gene, our aDNA results identified two previously unknown but related clades of Y. pestis associated with distinct medieval mass graves. These findings suggest that plague was imported to Europe on two or more occasions, each following a distinct route. These two clades are ancestral to modern isolates of Y. pestis biovars Orientalis and Medievalis. Our results clarify the etiology of the Black Death and provide a paradigm for a detailed historical reconstruction of the infection routes followed by this disease.

To date, characterization of ancient oral (dental calculus) and gut (coprolite) microbiota has been primarily accomplished through a metataxonomic approach involving targeted amplification of one or more variable regions in the 16S rRNA gene. Specifically, the V3 region (E. coli 341–534) of this gene has been suggested as an excellent candidate for ancient DNA amplification and microbial community reconstruction. However, in practice this metataxonomic approach often produces highly skewed taxonomic frequency data. In this study, we use non-targeted (shotgun metagenomics) sequencing methods to better understand skewed microbial profiles observed in four ancient dental calculus specimens previously analyzed by amplicon sequencing. Through comparisons of microbial taxonomic counts from paired amplicon (V3 U341F/534R) and shotgun sequencing datasets, we demonstrate that extensive length polymorphisms in the V3 region are a consistent and major cause of differential amplification leading to taxonomic bias in ancient microbiome reconstructions based on amplicon sequencing. We conclude that systematic amplification bias confounds attempts to accurately reconstruct microbiome taxonomic profiles from 16S rRNA V3 amplicon data generated using universal primers. Because in silico analysis indicates that alternative 16S rRNA hypervariable regions will present similar challenges, we advocate for the use of a shotgun metagenomics approach in ancient microbiome reconstructions.