Pathogen prospecting of museums: Reconstructing malaria epidemiology

Malaria is a disease of global significance. Ongoing changes to the earth's climate, antimalarial resistance, insecticide resistance, and socioeconomic decline test the resilience of malaria prevention programs. Museum insect specimens present an untapped resource for studying vector-borne pathogens, spurring the question: Do historical mosquito collections contain Plasmodium DNA, and, if so, can museum specimens be used to reconstruct the historical epidemiology of malaria? In this Perspective, we explore molecular techniques practical to pathogen prospecting, which, more broadly, we define as the science of screening entomological museum specimens for human, animal, or plant pathogens. Historical DNA and pathogen prospecting provide a means of describing the coevolution of human, vector, and parasite, informing the development of insecticides, diagnostics, therapeutics, and vaccines.

Time to rethink academic publishing: the peer reviewer crisis

There is concern that the time taken to publish academic papers in microbiological science has significantly increased in recent years. While the data do not specifically support this, evidence suggests that editors are having to invite more and more reviewers to identify those willing to perform peer review.

A 14th century CE Brucella melitensis genome and the recent expansion of the Western Mediterranean clade

Brucellosis is a disease caused by the bacterium Brucella and typically transmitted through contact with infected ruminants. It is one of the most common chronic zoonotic diseases and of particular interest to public health agencies. Despite its well-known transmission history and characteristic symptoms, we lack a more complete understanding of the evolutionary history of its best-known species-Brucella melitensis. To address this knowledge gap we fortuitously found, sequenced and assembled a high-quality ancient B. melitensis draft genome from the kidney stone of a 14th-century Italian friar. The ancient strain contained fewer core genes than modern B. melitensis isolates, carried a complete complement of virulence genes, and did not contain any indication of significant antimicrobial resistances. The ancient B. melitensis genome fell as a basal sister lineage to a subgroup of B. melitensis strains within the Western Mediterranean phylogenetic group, with a short branch length indicative of its earlier sampling time, along with a similar gene content. By calibrating the molecular clock we suggest that the speciation event between B. melitensis and B. abortus is contemporaneous with the estimated time frame for the domestication of both sheep and goats. These results confirm the existence of the Western Mediterranean clade as a separate group in the 14th CE and suggest that its divergence was due to human and ruminant co-migration.

Reply to Barton et al: signatures of natural selection during the Black Death

Barton et al.1 raise several statistical concerns regarding our original analyses2 that highlight the challenge of inferring natural selection using ancient genomic data. We show here that these concerns have limited impact on our original conclusions. Specifically, we recover the same signature of enrichment for high FST values at the immune loci relative to putatively neutral sites after switching the allele frequency estimation method to a maximum likelihood approach, filtering to only consider known human variants, and down-sampling our data to the same mean coverage across sites. Furthermore, using permutations, we show that the rs2549794 variant near ERAP2 continues to emerge as the strongest candidate for selection (p = 1.2×10-5), falling below the Bonferroni-corrected significance threshold recommended by Barton et al. Importantly, the evidence for selection on ERAP2 is further supported by functional data demonstrating the impact of the ERAP2 genotype on the immune response to Y. pestis and by epidemiological data from an independent group showing that the putatively selected allele during the Black Death protects against severe respiratory infection in contemporary populations.

Emergence, continuity, and evolution of Yersinia pestis throughout medieval and early modern Denmark

The historical epidemiology of plague is controversial due to the scarcity and ambiguity of available data.^1,2 A common source of debate is the extent and pattern of plague re-emergence and local continuity in Europe during the 14th-18th century CE.³ Despite having a uniquely long history of plague (∼5,000 years), Scandinavia is relatively underrepresented in the historical archives.^4,5 To better understand the historical epidemiology and evolutionary history of plague in this region, we performed in-depth (n = 298) longitudinal screening (800 years) for the plague bacterium Yersinia pestis (Y. pestis) across 13 archaeological sites in Denmark from 1000 to 1800 CE. Our genomic and phylogenetic data captured the emergence, continuity, and evolution of Y. pestis in this region over a period of 300 years (14th-17th century CE), for which the plague-positivity rate was 8.3% (3.3%-14.3% by site). Our phylogenetic analysis revealed that the Danish Y. pestis sequences were interspersed with those from other European countries, rather than forming a single cluster, indicative of the generation, spread, and replacement of bacterial variants through communities rather than their long-term local persistence. These results provide an epidemiological link between Y. pestis and the unknown pestilence that afflicted medieval and early modern Europe. They also demonstrate how population-scale genomic evidence can be used to test hypotheses on disease mortality and epidemiology and help pave the way for the next generation of historical disease research.

Examining pathogen DNA recovery across the remains of a 14th century Italian friar (Blessed Sante) infected with Brucella melitensis

OBJECTIVE

To investigate variation in ancient DNA recovery of Brucella melitensis, the causative agent of brucellosis, from multiple tissues belonging to one individual MATERIALS: 14 samples were analyzed from the mummified remains of the Blessed Sante, a 14 ^th century Franciscan friar from central Italy, with macroscopic diagnosis of probable brucellosis.

METHODS

Shotgun sequencing data from was examined to determine the presence of Brucella DNA.

RESULTS

Three of the 14 samples contained authentic ancient DNA, identified as belonging to B. melitensis. A genome (23.81X depth coverage, 0.98 breadth coverage) was recovered from a kidney stone. Nine of the samples contained reads classified as B. melitensis (7-169), but for many the data quality was insufficient to withstand our identification and authentication criteria.

CONCLUSIONS

We identified significant variation in the preservation and abundance of B. melitensis DNA present across multiple tissues, with calcified nodules yielding the highest number of authenticated reads. This shows how greatly sample selection can impact pathogen identification.

SIGNIFICANCE

Our results demonstrate variation in the preservation and recovery of pathogen DNA across tissues. This study highlights the importance of sample selection in the reconstruction of infectious disease burden and highlights the importance of a holistic approach to identifying disease.

LIMITATIONS

Study focuses on pathogen recovery in a single individual.

SUGGESTIONS FOR FURTHER RESEARCH

Further analysis of how sampling impacts aDNA recovery will improve pathogen aDNA recovery and advance our understanding of disease in past peoples.

Evolution of immune genes is associated with the Black Death

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.¹ In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.² Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia³-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.^4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,^6-12 in paleodietary reconstructions,^13-19 and facilitated genomic reconstructions from cave contexts.^8,20-22 Murchie et al.^6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.²⁴ Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 30,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

Collapse of the mammoth-steppe in central Yukon as revealed by ancient environmental DNA

The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

Response to Brinkmann et al. "Re-assembly of 19th century smallpox vaccine genomes reveals the contemporaneous use of horsepox and horsepox-related viruses in the United States"

We thank Brinkmann and colleagues for their correspondence and their further investigation into these American Civil War Era vaccination strains. Here, we summarize the difficulties and caveats of work with ancient DNA.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

Pleistocene population dynamics can inform the consequences of current climate change. This phylogeography of 35 complete American mastodon mitochondrial genomes suggests distinct lineages in this species repeatedly expanded northwards and then went locally extinct in response to glacial cycles.

The origins and genomic diversity of American Civil War Era smallpox vaccine strains

Vaccination has transformed public health, most notably including the eradication of smallpox. Despite its profound historical importance, little is known of the origins and diversity of the viruses used in smallpox vaccination. Prior to the twentieth century, the method, source and origin of smallpox vaccinations remained unstandardised and opaque. We reconstruct and analyse viral vaccine genomes associated with smallpox vaccination from historical artefacts. Significantly, we recover viral molecules through non-destructive sampling of historical materials lacking signs of biological residues. We use the authenticated ancient genomes to reveal the evolutionary relationships of smallpox vaccination viruses within the poxviruses as a whole.

Ancient Mitogenomes Reveal the Evolutionary History and Biogeography of Sloths

Living sloths represent two distinct lineages of small-sized mammals that independently evolved arboreality from terrestrial ancestors. The six extant species are the survivors of an evolutionary radiation marked by the extinction of large terrestrial forms at the end of the Quaternary. Until now, sloth evolutionary history has mainly been reconstructed from phylogenetic analyses of morphological characters. Here, we used ancient DNA methods to successfully sequence 10 extinct sloth mitogenomes encompassing all major lineages. This includes the iconic continental ground sloths Megatherium, Megalonyx, Mylodon, and Nothrotheriops and the smaller endemic Caribbean sloths Parocnus and Acratocnus. Phylogenetic analyses identify eight distinct lineages grouped in three well-supported clades, whose interrelationships are markedly incongruent with the currently accepted morphological topology. We show that recently extinct Caribbean sloths have a single origin but comprise two highly divergent lineages that are not directly related to living two-fingered sloths, which instead group with Mylodon. Moreover, living three-fingered sloths do not represent the sister group to all other sloths but are nested within a clade of extinct ground sloths including Megatherium, Megalonyx, and Nothrotheriops. Molecular dating also reveals that the eight newly recognized sloth families all originated between 36 and 28 million years ago (mya). The early divergence of recently extinct Caribbean sloths around 35 mya is consistent with the debated GAARlandia hypothesis postulating the existence at that time of a biogeographic connection between northern South America and the Greater Antilles. This new molecular phylogeny has major implications for reinterpreting sloth morphological evolution, biogeography, and diversification history.

Genetic resiliency and the Black Death: No apparent loss of mitogenomic diversity due to the Black Death in medieval London and Denmark

OBJECTIVES

In the 14th century AD, medieval Europe was severely affected by the Great European Famine as well as repeated bouts of disease, including the Black Death, causing major demographic shifts. This high volatility led to increased mobility and migration due to new labor and economic opportunities, as evidenced by documentary and stable isotope data. This study uses ancient DNA (aDNA) isolated from skeletal remains to examine whether evidence for large-scale population movement can be gleaned from the complete mitochondrial genomes of 264 medieval individuals from England (London) and Denmark.

MATERIALS AND METHODS

Using a novel library-conserving approach to targeted capture, we recovered 264 full mitochondrial genomes from the petrous portion of the temporal bones and teeth and compared genetic diversity across the medieval period within and between English (London) and Danish populations and with contemporary populations through population pairwise Φ_ST analysis.

RESULTS

We find no evidence of significant differences in genetic diversity spatially or temporally in our dataset, yet there is a high degree of haplotype diversity in our medieval samples with little exact sequence sharing.

DISCUSSION

The mitochondrial genomes of both medieval Londoners and medieval Danes suggest high mitochondrial diversity before, during and after the Black Death. While our mitochondrial genomic data lack geographically correlated signals, these data could be the result of high, continual female migration before and after the Black Death or may simply indicate a large female effective population size unaffected by the upheaval of the medieval period. Either scenario suggests a genetic resiliency in areas of northwestern medieval Europe.

Resolving the phylogenetic position of Darwin's extinct ground sloth (Mylodon darwinii) using mitogenomic and nuclear exon data

Mylodon darwinii is the extinct giant ground sloth named after Charles Darwin, who first collected its remains in South America. We have successfully obtained a high-quality mitochondrial genome at 99-fold coverage using an Illumina shotgun sequencing of a 12 880-year-old bone fragment from Mylodon Cave in Chile. Low level of DNA damage showed that this sample was exceptionally well preserved for an ancient subfossil, probably the result of the dry and cold conditions prevailing within the cave. Accordingly, taxonomic assessment of our shotgun metagenomic data showed a very high percentage of endogenous DNA with 22% of the assembled metagenomic contigs assigned to Xenarthra. Additionally, we enriched over 15 kb of sequence data from seven nuclear exons, using target sequence capture designed against a wide xenarthran dataset. Phylogenetic and dating analyses of the mitogenomic dataset including all extant species of xenarthrans and the assembled nuclear supermatrix unambiguously place Mylodon darwinii as the sister-group of modern two-fingered sloths, from which it diverged around 22 million years ago. These congruent results from both the mitochondrial and nuclear data support the diphyly of the two modern sloth lineages, implying the convergent evolution of their unique suspensory behaviour as an adaption to arboreality. Our results offer promising perspectives for whole-genome sequencing of this emblematic extinct taxon.

Dorset Pre-Inuit and Beothuk foodways in Newfoundland, ca. AD 500-1829

Archaeological research on the Canadian island of Newfoundland increasingly demonstrates that the island's subarctic climate and paucity of terrestrial food resources did not restrict past Pre-Inuit (Dorset) and Native American (Beothuk) hunter-gatherer populations to a single subsistence pattern. This study first sought to characterize hunter-gatherer diets over the past 1500 years; and second, to assess the impact of European colonization on Beothuk lifeways by comparing the bone chemistry of Beothuk skeletal remains before and after the intensification of European settlement in the early 18th century. We employed radiocarbon dating and stable carbon and nitrogen isotope ratio analysis of bulk bone collagen from both Dorset (n = 9) and Beothuk (n = 13) cultures, including a naturally mummified 17th century Beothuk individual. Carbon and nitrogen isotope analysis of 108 faunal samples from Dorset and Beothuk archaeological sites around the island were used as a dietary baseline for the humans. We combined our results with previously published isotope data and radiocarbon dates from Dorset (n = 12) and Beothuk (n = 18) individuals and conducted a palaeodietary analysis using Bayesian modelling, cluster analysis and comparative statistical tests. Dorset diets featured more marine protein than those of the Beothuk, and the diets of Beothuk after the 18th century featured less high trophic level marine protein than those of individuals predating the 18th century. Despite inhabiting the same island, Dorset and Beothuk cultures employed markedly different dietary strategies, consistent with interpretations of other archaeological data. Significantly, European colonization had a profound effect on Beothuk lifeways, as in response to the increasing European presence on the coast, the Beothuk relied more extensively on the limited resources of the island's boreal forests and rivers.

The paradox of HBV evolution as revealed from a 16th century mummy

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.

Hepatitis 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.

Adaptation in a Fibronectin Binding Autolysin of Staphylococcus saprophyticus

Human-pathogenic bacteria are found in a variety of niches, including free-living, zoonotic, and microbiome environments. Identifying bacterial adaptations that enable invasive disease is an important means of gaining insight into the molecular basis of pathogenesis and understanding pathogen emergence. Staphylococcus saprophyticus, a leading cause of urinary tract infections, can be found in the environment, food, animals, and the human microbiome. We identified a selective sweep in the gene encoding the Aas adhesin, a key virulence factor that binds host fibronectin. We hypothesize that the mutation under selection (aas_2206A>C) facilitates colonization of the urinary tract, an environment where bacteria are subject to strong shearing forces. The mutation appears to have enabled emergence and expansion of a human-pathogenic lineage of S. saprophyticus. These results demonstrate the power of evolutionary genomic approaches in discovering the genetic basis of virulence and emphasize the pleiotropy and adaptability of bacteria occupying diverse niches. IMPORTANCEStaphylococcus saprophyticus is an important cause of urinary tract infections (UTI) in women; such UTI are common, can be severe, and are associated with significant impacts to public health. In addition to being a cause of human UTI, S. saprophyticus can be found in the environment, in food, and associated with animals. After discovering that UTI strains of S. saprophyticus are for the most part closely related to each other, we sought to determine whether these strains are specially adapted to cause disease in humans. We found evidence suggesting that a mutation in the gene aas is advantageous in the context of human infection. We hypothesize that the mutation allows S. saprophyticus to survive better in the human urinary tract. These results show how bacteria found in the environment can evolve to cause disease.

Comment: Characterization of Two Historic Smallpox Specimens from a Czech Museum

The complete genome sequences of two strains of variola virus (VARV) sampled from human smallpox specimens present in the Czech National Museum, Prague, were recently determined, with one of the sequences estimated to date to the mid-19th century. Using molecular clock methods, the authors of this study go on to infer that the currently available strains of VARV share an older common ancestor, at around 1350 AD, than some recent estimates based on other archival human samples. Herein, we show that the two Czech strains exhibit anomalous branch lengths given their proposed age, and by assuming a constant rate of evolutionary change across the rest of the VARV phylogeny estimate that their true age in fact lies between 1918 and 1937. We therefore suggest that the age of the common ancestor of currently available VARV genomes most likely dates to late 16th and early 17th centuries and not ~1350 AD.

A molecular portrait of maternal sepsis from Byzantine Troy

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

DOI:http://dx.doi.org/10.7554/eLife.20983.001

Why and how have some bacteria evolved to cause illness in humans? One way to study bacterial evolution is to search for ancient samples of bacteria and use DNA sequencing technology to investigate how modern bacteria have changed from their ancestors. Understanding the evolution process may help researchers to understand how some bacteria become resistant to the antibiotics designed to kill them.

Complications that occur during pregnancy, including bacterial infections, have long been a major cause of death for women. Now, Devault, Mortimer et al. have been able to sequence the DNA of bacteria found in tissue collected from a woman buried 800 years ago in a cemetery in Troy. Some of the woman’s tissues had been well preserved because they had calcified (probably as the result of infection), which preserved their structure in a mineralized layer. Two mineralized “nodules” in the body appear to be the remains of abscesses. Some of the human DNA in the nodules came from a male, suggesting that the woman was pregnant with a boy and that the abscesses formed in placental tissue.

Sequencing the DNA of the bacteria in the abscess allowed Devault, Mortimer et al. to diagnose the woman’s infection, which was caused by two types of bacteria. One species, called Gardnerella vaginalis, is found in modern pregnancy-related infections. The DNA of the ancient samples was similar to that of modern bacteria. The other bacteria species was an ancient form of Staphylococcus saprophyticus, a type of bacteria that causes urinary tract infections. However, the DNA of the ancient S. saprophyticus bacteria is quite different to that of the bacteria found in modern humans. Instead, their DNA sequence appears more similar to forms of the bacteria that infect currently livestock. As humans lived closely with their livestock at the time the woman lived, her infection may be due to a type of bacteria that passed easily between humans and animals.

Overall, the results suggest that the disease-causing properties of bacteria can arise from a wide range of sources. In addition, Devault, Mortimer et al. have demonstrated that certain types of tissue found in archeological remains are a potential gold mine of information about the evolution of bacteria and other microbes found in the human body.

DOI:http://dx.doi.org/10.7554/eLife.20983.002

Eighteenth century Yersinia pestis genomes reveal the long-term persistence of an historical plague focus

The 14th–18th century pandemic of Yersinia pestis caused devastating disease outbreaks in Europe for almost 400 years. The reasons for plague’s persistence and abrupt disappearance in Europe are poorly understood, but could have been due to either the presence of now-extinct plague foci in Europe itself, or successive disease introductions from other locations. Here we present five Y. pestis genomes from one of the last European outbreaks of plague, from 1722 in Marseille, France. The lineage identified has not been found in any extant Y. pestis foci sampled to date, and has its ancestry in strains obtained from victims of the 14th century Black Death. These data suggest the existence of a previously uncharacterized historical plague focus that persisted for at least three centuries. We propose that this disease source may have been responsible for the many resurgences of plague in Europe following the Black Death.

DOI:http://dx.doi.org/10.7554/eLife.12994.001

A bacterium called Yersina pestis is responsible for numerous human outbreaks of plague throughout history. It is carried by rats and other rodents and can spread to humans causing what we conventionally refer to as plague. The most notorious of these plague outbreaks – the Black Death – claimed millions of lives in Europe in the mid-14^th century. Several other plague outbreaks emerged in Europe over the next 400 years. Then, there was a large gap before the plague re-emerged as threat in the 19^th century and it continues to infect humans today, though on a smaller scale.

Scientists have extensively studied Y. pestis to understand its origin and how it evolved to become such a deadly threat. These studies led to the assumption that the plague outbreaks of the 14–18^th centuries likely originated in rodents in Asia and spread along trade routes to other parts of the world. However, it is not clear why the plague persisted in Europe for 400 years after the Black Death. Could the bacteria have gained a foothold in local rodents instead of being reintroduced from Asia each time? If it did, why did it then disappear for such a long period from the end of the 18^th century?

To help answer these questions, Bos, Herbig et al. sequenced the DNA of Y. pestis samples collected from the teeth of five individuals who died of plague during the last major European outbreak of plague in 1722 in Marseille, France. The DNA sequences of these bacterial samples were then compared with the DNA sequences of modern day Y. pestis and other historical samples of the bacteria. The results showed the bacteria in the Marseille outbreak likely evolved from the strain that caused the Black Death back in the 14^th century.

The comparisons showed that the strain isolated from the teeth is not found today, and may be extinct. This suggests that a historical reservoir for plague existed somewhere, perhaps in Asia, or perhaps in Europe itself, and was able to cause outbreaks up until the 18^th century.Bos, Herbig et al.’s findings may help researchers trying to control the current outbreaks of the plague in Madagascar and other places.

DOI:http://dx.doi.org/10.7554/eLife.12994.002

Shotgun Mitogenomics Provides a Reference Phylogenetic Framework and Timescale for Living Xenarthrans

Xenarthra (armadillos, sloths, and anteaters) constitutes one of the four major clades of placental mammals. Despite their phylogenetic distinctiveness in mammals, a reference phylogeny is still lacking for the 31 described species. Here we used Illumina shotgun sequencing to assemble 33 new complete mitochondrial genomes, establishing Xenarthra as the first major placental clade to be fully sequenced at the species level for mitogenomes. The resulting data set allowed the reconstruction of a robust phylogenetic framework and timescale that are consistent with previous studies conducted at the genus level using nuclear genes. Incorporating the full species diversity of extant xenarthrans points to a number of inconsistencies in xenarthran systematics and species definition. We propose to split armadillos into two distinct families Dasypodidae (dasypodines) and Chlamyphoridae (euphractines, chlamyphorines, and tolypeutines) to better reflect their ancient divergence, estimated around 42 Ma. Species delimitation within long-nosed armadillos (genus Dasypus) appeared more complex than anticipated, with the discovery of a divergent lineage in French Guiana. Diversification analyses showed Xenarthra to be an ancient clade with a constant diversification rate through time with a species turnover driven by high but constant extinction. We also detected a significant negative correlation between speciation rate and past temperature fluctuations with an increase in speciation rate corresponding to the general cooling observed during the last 15 My. Biogeographic reconstructions identified the tropical rainforest biome of Amazonia and the Guiana Shield as the cradle of xenarthran evolutionary history with subsequent dispersions into more open and dry habitats.

Correction to 'Ancient DNA and the tropics: a rodent's tale'

The present erratum is in regards to our article entitled ‘Ancient DNA and the tropics: a rodent's tale’. We were made aware of problems with some of the ancient sequences submitted to GenBank and conducted a systematic review of all the files used in our study. We discovered that, unfortunately, an incorrect file was sent to GenBank and was also used in some of our downstream analyses. We immediately contacted GenBank, explained the situation and corrected the file. We have redone some analyses with the correct file and describe these changes below.

Ancient DNA and the tropics: a rodent's tale

Most genetic studies of Holocene fauna have been performed with ancient samples from dry and cold regions, in which preservation of fossils is facilitated and molecular damage is reduced. Ancient DNA work from tropical regions has been precluded owing to factors that limit DNA preservation (e.g. temperature, hydrolytic damage). We analysed ancient DNA from rodent jawbones identified as Ototylomys phyllotis, found in Holocene and Late Pleistocene stratigraphic layers from Loltún, a humid tropical cave located in the Yucatan peninsula. We extracted DNA and amplified six short overlapping fragments of the cytochrome b gene, totalling 666 bp, which represents an unprecedented success considering tropical ancient DNA samples. We performed genetic, phylogenetic and divergence time analyses, combining sequences from ancient and modern O. phyllotis, in order to assess the ancestry of the Loltún samples. Results show that all ancient samples fall into a unique clade that diverged prior to the divergence of the modern O. phyllotis, supporting it as a distinct Pleistocene form of the Ototylomys genus. Hence, this rodent's tale suggests that the sister group to modern O. phyllotis arose during the Miocene-Pliocene, diversified during the Pleistocene and went extinct in the Holocene.

Ancient human genomics: the methodology behind reconstructing evolutionary pathways

High-throughput sequencing (HTS) has radically altered approaches to human evolutionary research. Recent contributions highlight that HTS is able to reach depths of the human lineage previously thought to be impossible. In this paper, we outline the methodological advances afforded by recent developments in DNA recovery, data output, scalability, speed, and resolution of the current sequencing technology. We review and critically evaluate the 'DNA pipeline' for ancient samples: from DNA extraction, to constructing immortalized sequence libraries, to enrichment strategies (e.g., polymerase chain reaction [PCR] and hybridization capture), and finally, to bioinformatic analyses of sequence data. We argue that continued evaluations and improvements to this process are essential to ensure sequence data validity. Also, we highlight the role of contamination and authentication in ancient DNA-HTS, which is particularly relevant to ancient human genomics, since sequencing the genomes of hominins such as Homo erectus and Homo heidelbergensis may soon be within the realm of possibility.

Ancient pathogen DNA in archaeological samples detected with a Microbial Detection Array

Ancient human remains of paleopathological interest typically contain highly degraded DNA in which pathogenic taxa are often minority components, making sequence-based metagenomic characterization costly. Microarrays may hold a potential solution to these challenges, offering a rapid, affordable, and highly informative snapshot of microbial diversity in complex samples without the lengthy analysis and/or high cost associated with high-throughput sequencing. Their versatility is well established for modern clinical specimens, but they have yet to be applied to ancient remains. Here we report bacterial profiles of archaeological and historical human remains using the Lawrence Livermore Microbial Detection Array (LLMDA). The array successfully identified previously-verified bacterial human pathogens, including Vibrio cholerae (cholera) in a 19th century intestinal specimen and Yersinia pestis (“Black Death” plague) in a medieval tooth, which represented only minute fractions (0.03% and 0.08% alignable high-throughput shotgun sequencing reads) of their respective DNA content. This demonstrates that the LLMDA can identify primary and/or co-infecting bacterial pathogens in ancient samples, thereby serving as a rapid and inexpensive paleopathological screening tool to study health across both space and time.

Ancient whole genome enrichment using baits built from modern DNA

We report metrics from complete genome capture of nuclear DNA from extinct mammoths using biotinylated RNAs transcribed from an Asian elephant DNA extract. Enrichment of the nuclear genome ranged from 1.06- to 18.65-fold, to an apparent maximum threshold of ∼80% on-target. This projects an order of magnitude less costly complete genome sequencing from long-dead organisms, even when a reference genome is unavailable for bait design.

Second-pandemic strain of Vibrio cholerae from the Philadelphia cholera outbreak of 1849

In the 19th century, there were several major cholera pandemics in the Indian subcontinent, Europe, and North America. The causes of these outbreaks and the genomic strain identities remain a mystery. We used targeted high-throughput sequencing to reconstruct the Vibrio cholerae genome from the preserved intestine of a victim of the 1849 cholera outbreak in Philadelphia, part of the second cholera pandemic. This O1 biotype strain has 95 to 97% similarity with the classical O395 genome, differing by 203 single-nucleotide polymorphisms (SNPs), lacking three genomic islands, and probably having one or more tandem cholera toxin prophage (CTX) arrays, which potentially affected its virulence. This result highlights archived medical remains as a potential resource for investigations into the genomic origins of past pandemics.

Amplicon pyrosequencing late Pleistocene permafrost: the removal of putative contaminant sequences and small-scale reproducibility

DNA sequencing of ancient permafrost samples can be used to reconstruct past plant, animal and bacterial communities. In this study, we assess the small-scale reproducibility of taxonomic composition obtained from sequencing four molecular markers (mitochondrial 12S ribosomal DNA (rDNA), prokaryote 16S rDNA, mitochondrial cox1 and chloroplast trnL intron) from two soil cores sampled 10 cm apart. In addition, sequenced control reactions were used to produce a contaminant library that was used to filter similar sequences from sample libraries. Contaminant filtering resulted in the removal of 1% of reads or 0.3% of operational taxonomic units. We found similar richness, overlap, abundance and taxonomic diversity from the 12S, 16S and trnL markers from each soil core. Jaccard dissimilarity across the two soil cores was highest for metazoan taxa detected by the 12S and cox1 markers. Taxonomic community distances were similar for each marker across the two soil cores when the chi-squared metric was used; however, the 12S and cox1 markers did not cluster well when the Goodall similarity metric was used. A comparison of plant macrofossil vs. read abundance corroborates previous work that suggests eastern Beringia was dominated by grasses and forbs during cold stages of the Pleistocene, a habitat that is restricted to isolated sites in the present-day Yukon.

Yersinia pestis: new evidence for an old infection

The successful reconstruction of an ancient bacterial genome from archaeological material presents an important methodological advancement for infectious disease research. The reliability of evolutionary histories inferred by the incorporation of ancient data, however, are highly contingent upon the level of genetic diversity represented in modern genomic sequences that are publicly accessible, and the paucity of available complete genomes restricts the level of phylogenetic resolution that can be obtained. Here we add to our original analysis of the Yersinia pestis strain implicated in the Black Death by consolidating our dataset for 18 modern genomes with single nucleotide polymorphism (SNP) data for an additional 289 strains at over 600 positions. The inclusion of this additional data reveals a cluster of Y. pestis strains that diverge at a time significantly in advance of the Black Death, with divergence dates roughly coincident with the Plague of Justinian (6^th to 8^th century AD). In addition, the analysis reveals further clues regarding potential radiation events that occurred immediately preceding the Black Death, and the legacy it may have left in modern Y. pestis populations. This work reiterates the need for more publicly available complete genomes, both modern and ancient, to achieve an accurate understanding of the history of this bacterium.

Isolation of nucleic acids and cultures from fossil ice and permafrost

Owing to their constant low temperatures, glacial ice and permafrost might contain the oldest nucleic acids and microbial cells on Earth, which could prove key to reconstructing past ecosystems and for the planning of missions to other planets. However, recent claims concerning viable cells and microbial nucleic acids obtained from ice- and permafrost cores from hundreds of thousands to millions of years old are not properly authenticated and the findings could be the result of contamination. Here, we discuss the processes that restrict the long-term survival of DNA and/or RNA molecules in ice and permafrost, and highlight sources of contamination that could result in false claims. Additionally, we present a set of precautions, controls and criteria to help ensure that future cultures and sequences are authentic.

Herald waves of cholera in nineteenth century London

Deaths from cholera in London, UK, were recorded weekly from 1824 to 1901. Three features of the time series stand out: (i) cholera deaths were strongly seasonal, with peak mortality almost always in the summer, (ii) the only non-summer outbreaks occurred in the spring of 1832, the autumn of 1848 and the winter of 1853, and (iii) extraordinarily severe summer outbreaks occurred in 1832, 1849, 1854 and 1866 (the four 'great' cholera years). The non-summer outbreaks of 1832, 1848 and 1853 appear to have been herald waves of newly invading cholera strains. In addition, a simple mathematical model confirms that a non-summer introduction of a new cholera strain can result in an initial herald wave, followed by a severe outbreak the following summer. Through the analysis of the genomes of nineteenth-century specimens, it may be possible to identify the strains that caused these herald waves and the well-known cholera epidemics that followed.

Ancient DNA from lake sediments: bridging the gap between paleoecology and genetics

Background

Quaternary plant ecology in much of the world has historically relied on morphological identification of macro- and microfossils from sediments of small freshwater lakes. Here, we report new protocols that reliably yield DNA sequence data from Holocene plant macrofossils and bulk lake sediment used to infer ecological change. This will allow changes in census populations, estimated from fossils and associated sediment, to be directly associated with population genetic changes.

Results

We successfully sequenced DNA from 64 samples (out of 126) comprised of bulk sediment and seeds, leaf fragments, budscales, and samaras extracted from Holocene lake sediments in the western Great Lakes region of North America. Overall, DNA yields were low. However, we were able to reliably amplify samples with as few as 10 copies of a short cpDNA fragment with little detectable PCR inhibition. Our success rate was highest for sediments < 2000 years old, but we were able to successfully amplify DNA from samples up to 4600 years old. DNA sequences matched the taxonomic identity of the macrofossil from which they were extracted 79% of the time. Exceptions suggest that DNA molecules from surrounding nearby sediments may permeate or adhere to macrofossils in sediments.

Conclusions

An ability to extract ancient DNA from Holocene sediments potentially allows exciting new insights into the genetic consequences of long-term environmental change. The low DNA copy numbers we found in fossil material and the discovery of multiple sequence variants from single macrofossil extractions highlight the need for careful experimental and laboratory protocols. Further application of these protocols should lead to better understanding of the ecological and evolutionary consequences of environmental change.

Comparison of methods in the recovery of nucleic acids from archival formalin-fixed paraffin-embedded autopsy tissues

Archival formalin-fixed paraffin-embedded (FFPE) human tissue collections are typically in poor states of storage across the developing world. With advances in biomolecular techniques, these extraordinary and virtually untapped resources have become an essential part of retrospective epidemiological studies. To successfully use such tissues in genomic studies, scientists require high nucleic acid yields and purity. In spite of the increasing number of FFPE tissue kits available, few studies have analyzed their applicability in recovering high-quality nucleic acids from archived human autopsy samples. Here we provide a study involving 10 major extraction methods used to isolate total nucleic acid from FFPE tissues ranging in age from 3 to 13years. Although all 10 methods recovered quantifiable amounts of DNA, only 6 recovered quantifiable RNA, varying considerably and generally yielding lower DNA concentrations. Overall, we show quantitatively that TrimGen's WaxFree method and our in-house phenol-chloroform extraction method recovered the highest yields of amplifiable DNA, with considerable polymerase chain reaction (PCR) inhibition, whereas Ambion's RecoverAll method recovered the most amplifiable RNA.

Out of America: ancient DNA evidence for a new world origin of late quaternary woolly mammoths

Although the iconic mammoth of the Late Pleistocene, the woolly mammoth (Mammuthus primigenius), has traditionally been regarded as the end point of a single anagenetically evolving lineage, recent paleontological and molecular studies have shown that successive allopatric speciation events must have occurred within Pleistocene Mammuthus in Asia, with subsequent expansion and hybridization between nominal taxa [1, 2]. However, the role of North American mammoth populations in these events has not been adequately explored from an ancient-DNA standpoint. To undertake this task, we analyzed mtDNA from a large data set consisting of mammoth samples from across Holarctica (n = 160) and representing most of radiocarbon time. Our evidence shows that, during the terminal Pleistocene, haplotypes originating in and characteristic of New World populations replaced or succeeded those endemic to Asia and western Beringia. Also, during the Last Glacial Maximum, mammoth populations do not appear to have suffered an overall decline in diversity, despite differing responses on either side of the Bering land bridge. In summary, the "Out-of-America" hypothesis holds that the dispersal of North American woolly mammoths into other parts of Holarctica created major phylogeographic structuring within Mammuthus primigenius populations, shaping the last phase of their evolutionary history before their demise.

Ancient biomolecules from deep ice cores reveal a forested southern Greenland

It is difficult to obtain fossil data from the 10% of Earth's terrestrial surface that is covered by thick glaciers and ice sheets, and hence, knowledge of the paleoenvironments of these regions has remained limited. We show that DNA and amino acids from buried organisms can be recovered from the basal sections of deep ice cores, enabling reconstructions of past flora and fauna. We show that high-altitude southern Greenland, currently lying below more than 2 kilometers of ice, was inhabited by a diverse array of conifer trees and insects within the past million years. The results provide direct evidence in support of a forested southern Greenland and suggest that many deep ice cores may contain genetic records of paleoenvironments in their basal sections.

Metagenomics to Paleogenomics: Large-Scale Sequencing of Mammoth DNA

We sequenced 28 million base pairs of DNA in a metagenomics approach, using a woolly mammoth (Mammuthus primigenius) sample from Siberia. As a result of exceptional sample preservation and the use of a recently developed emulsion polymerase chain reaction and pyrosequencing technique, 13 million base pairs (45.4%) of the sequencing reads were identified as mammoth DNA. Sequence identity between our data and African elephant (Loxodonta africana) was 98.55%, consistent with a paleontologically based divergence date of 5 to 6 million years. The sample includes a surprisingly small diversity of environmental DNAs. The high percentage of endogenous DNA recoverable from this single mammoth would allow for completion of its genome, unleashing the field of paleogenomics.

Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes

Analysis of mitochondrial DNA sequence variation has been used extensively to study the evolutionary relationships of individuals and populations, both within and across species. So ubiquitous and easily acquired are mtDNA data that it has been suggested that such data could serve as a taxonomic 'barcode' for an objective species classification scheme. However, there are technical pitfalls associated with the acquisition of mtDNA data. One problem is the presence of translocated pieces of mtDNA in the nuclear genome of many taxa that may be mistaken for authentic organellar mtDNA. We assessed the extent to which such 'numt' sequences may pose an overlooked problem in analyses of mtDNA from humans and apes. Using long-range polymerase chain reaction (PCR), we generated necessarily authentic mtDNA sequences for comparison with sequences obtained using typical methods for a segment of the mtDNA control region in humans, chimpanzees, bonobos, gorillas and orangutans. Results revealed that gorillas are notable for having such a variety of numt sequences bearing high similarity to authentic mtDNA that any analysis of mtDNA using standard approaches is rendered impossible. Studies on humans, chimpanzees, bonobos or orangutans are apparently less problematic. One implication is that explicit measures need to be taken to authenticate mtDNA sequences in newly studied taxa or when any irregularities arise. Furthermore, some taxa may not be amenable to analysis of mtDNA variation at all.

The genetic secrets some fossils hold

Most animals that once lived have gone extinct. The remains of a few of these can be found in museum collections worldwide. As modern evolutionary biology is limited to the use of extant taxa, retrieving DNA from extinct or subfossil organisms can add significant insight into past population history and resolve phylogenies that can be tentative by morphology alone. DNA is a relatively weak molecule, comparatively speaking, yet under certain conditions it persists in the fossil record, despite what in vitro chemistry predicts. While most fossil remains do not contain DNA, museum specimens can be screened for the presence of conditions that would be conducive for nucleic acid preservation by measuring the extent of amino acid racemization and by looking at the extent of protein hydrolysis by pyrolysis gas chromatography/mass spectrometry. Results from these types of analyses suggest that the preservation of DNA is linked to the temperature and its constancy at a site rather than its age. Chemical analyses of coprolites from extinct herbivores from the late Pleistocene, as well as Archaic Native Americans, show the presence of compounds from the Maillard reaction. Upon the cleaving of these products, the defecator can be identified and his diet analyzed.

Molecular analysis of an 11,700-year-old rodent midden from the Atacama Desert, Chile

DNA was extracted from an 11,700-year-old rodent midden from the Atacama Desert, Chile and the chloroplast and animal mitochondrial DNA (mtDNA) gene sequences were analysed to investigate the floral environment surrounding the midden, and the identity of the midden agent. The plant sequences, together with the macroscopic identifications, suggest the presence of 13 plant families and three orders that no longer exist today at the midden locality, and thus point to a much more diverse and humid climate 11,700 years ago. The mtDNA sequences suggest the presence of at least four different vertebrates, which have been putatively identified as a camelid (vicuna), two rodents (Phyllotis and Abrocoma), and a cardinal bird (Passeriformes). To identify the midden agent, DNA was extracted from pooled faecal pellets, three small overlapping fragments of the mitochondrial cytochrome b gene were amplified and multiple clones were sequenced. These results were analysed along with complete cytochrome b sequences for several modern Phyllotis species to place the midden sequence phylogenetically. The results identified the midden agent as belonging to an ancestral P. limatus. Today, P. limatus is not found at the midden locality but it can be found 100 km to the north, indicating at least a small range shift. The more extensive sampling of modern Phyllotis reinforces the suggestion that P. limatus is recently derived from a peripheral isolate.

Ancient DNA

DNA that has been recovered from archaeological and palaeontological remains makes it possible to go back in time and study the genetic relationships of extinct organisms to their contemporary relatives. This provides a new perspective on the evolution of organisms and DNA sequences. However, the field is fraught with technical pitfalls and needs stringent criteria to ensure the reliability of results, particularly when human remains are studied.

A molecular analysis of dietary diversity for three archaic Native Americans

DNA was extracted from three fecal samples, more than 2,000 years old, from Hinds Cave, Texas. Amplification of human mtDNA sequences showed their affiliation with contemporary Native Americans, while sequences from pronghorn antelope, bighorn sheep, and cottontail rabbit allowed these animals to be identified as part of the diet of these individuals. Furthermore, amplification of chloroplast DNA sequences identified eight different plants as dietary elements. These archaic humans consumed 2-4 different animal species and 4-8 different plant species during a short time period. The success rate for retrieval of DNA from paleofeces is in strong contrast to that from skeletal remains where the success rate is generally low. Thus, human paleofecal remains represent a source of ancient DNA that significantly complements and may in some cases be superior to that from skeletal tissue.

A molecular analysis of ground sloth diet through the last glaciation

DNA was extracted from five coprolites, excavated in Gypsum Cave, Nevada and radiocarbon dated to approximately 11 000, 20 000 and 28 500 years BP. All coprolites contained mitochondrial DNA sequences identical to a DNA sequence determined from a bone of the extinct ground sloth Nothrotheriops shastensis. A 157-bp fragment of the chloroplast gene for the large subunit of the ribulosebisphosphate carboxylase (rbcL) was amplified from the boluses and several hundred clones were sequenced. In addition, the same DNA fragment was sequenced from 99 plant species that occur in the vicinity of Gypsum Cave today. When these were compared to the DNA sequences in GenBank, 69 were correctly (two incorrectly) assigned to taxonomic orders. The plant sequences from the five coprolites as well as from one previously studied coprolite were compared to rbcL sequences in GenBank and the contemporary plant species. Thirteen families or orders of plants that formed part of the diet of the Shasta ground sloth could be identified, showing that the ground sloth was feeding on trees as well as herbs and grasses. The plants in the boluses further indicate that the climate 11 000 years BP was dryer than 20 000 and 28 500 years BP. However, the sloths seem to have visited water sources more frequently at 11 000 BP than at earlier times.