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

Plagued by a cryptic clock: insight and issues from the global phylogeny of Yersinia pestis.

Plague has an enigmatic history as a zoonotic pathogen. This infectious disease will unexpectedly appear in human populations and disappear just as suddenly. As a result, a long-standing line of inquiry has been to estimate when and where plague appeared in the past. However, there have been significant disparities between phylogenetic studies of the causative bacterium, Yersinia pestis, regarding the timing and geographic origins of its reemergence. Here, we curate and contextualize an updated phylogeny of Y. pestis using 601 genome sequences sampled globally. Through a detailed Bayesian evaluation of temporal signal in subsets of these data we demonstrate that a Y. pestis-wide molecular clock is unstable. To resolve this, we developed a new approach in which each Y. pestis population was assessed independently, enabling us to recover substantial temporal signal in five populations, including the ancient pandemic lineages which we now estimate may have emerged decades, or even centuries, before a pandemic was historically documented from European sources. Despite this methodological advancement, we only obtain robust divergence dates from populations sampled over a period of at least 90 years, indicating that genetic evidence alone is insufficient for accurately reconstructing the timing and spread of short-term plague epidemics.

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.

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.

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.1 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.2 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 Beringia3-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.24 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 40,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.

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.

Capturing the Resistome: a Targeted Capture Method To Reveal Antibiotic Resistance Determinants in Metagenomes.

Identification of the nucleotide sequences encoding antibiotic resistance elements and determination of their association with antibiotic resistance are critical to improve surveillance and monitor trends in antibiotic resistance. Current methods to study antibiotic resistance in various environments rely on extensive deep sequencing or laborious culturing of fastidious organisms, both of which are heavily time-consuming operations. An accurate and sensitive method to identify both rare and common resistance elements in complex metagenomic samples is needed. Referencing the sequences in the Comprehensive Antibiotic Resistance Database, we designed a set of 37,826 probes to specifically target over 2,000 nucleotide sequences associated with antibiotic resistance in clinically relevant bacteria. Testing of this probe set on DNA libraries generated from multidrug-resistant bacteria to selectively capture resistance genes reproducibly produced higher numbers of reads on target at a greater length of coverage than shotgun sequencing. We also identified additional resistance gene sequences from human gut microbiome samples that sequencing alone was not able to detect. Our method to capture the resistome enables a sensitive means of gene detection in diverse environments where genes encoding antibiotic resistance represent less than 0.1% of the metagenome.

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.

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.

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.

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

[This corrects the article DOI: 10.1371/journal.ppat.1006750.].

A Single Amino Acid Change in the Response Regulator PhoP, Acquired during Yersinia pestis Evolution, Affects PhoP Target Gene Transcription and Polymyxin B Susceptibility.

Yersinia pestis, the causative agent of plague, evolved from the closely related pathogen Yersinia pseudotuberculosis During its emergence, Y. pestis is believed to have acquired its unique pathogenic characteristics through numerous gene gains/losses, genomic rearrangements, and single nucleotide polymorphism (SNP) changes. One such SNP creates a single amino acid variation in the DNA binding domain of PhoP, the response regulator in the PhoP/PhoQ two-component system. Y. pseudotuberculosis and the basal human-avirulent strains of Y. pestis harbor glycines at position 215 of PhoP, whereas the modern human-virulent strains (e.g., KIM and CO92) harbor serines at this residue. Since PhoP plays multiple roles in the adaptation of Y. pestis to stressful host conditions, we tested whether this amino acid substitution affects PhoP activity or the ability of Y. pestis to survive in host environments. Compared to the parental KIM6+ strain carrying the modern allele of phoP (phoP-S215), a derivative carrying the basal allele (phoP-G215) exhibited slightly defective growth under a low-Mg2+ condition and decreased transcription of a PhoP target gene, ugd, as well as an ∼8-fold increase in the susceptibility to the antimicrobial peptide polymyxin B. The phoP-G215 strain showed no apparent defect in flea colonization, although a phoP-null mutant showed decreased flea infectivity in competition experiments. Our results suggest that the amino acid variation at position 215 of PhoP causes subtle changes in the PhoP activity and raise the possibility that the change in this residue have contributed to the evolution of increased virulence in Y. pestisIMPORTANCEY. pestis acquired a single nucleotide polymorphism (SNP) in phoP when the highly human-virulent strains diverged from less virulent basal strains, resulting in an amino acid substitution in the DNA binding domain of the PhoP response regulator. We show that Y. pestis carrying the modern phoP allele has an increased ability to induce the PhoP-regulated ugd gene and resist antimicrobial peptides compared to an isogenic strain carrying the basal allele. Given the important roles PhoP plays in host adaptation, the results raise an intriguing possibility that this amino acid substitution contributed to the evolution of increased virulence in Y. pestis Additionally, we present the first evidence that phoP confers a survival fitness advantage to Y. pestis inside the flea midgut.