key: cord-0962076-3h59uw8r
authors: Duchêne, Sebastián; Ho, Simon Y.W.; Carmichael, Ann G.; Holmes, Edward C.; Poinar, Hendrik
title: The Recovery, Interpretation and Use of Ancient Pathogen Genomes
date: 2020-10-05
journal: Curr Biol
DOI: 10.1016/j.cub.2020.08.081
sha: f82da78af5fe13515641614909aebc163d67af10
doc_id: 962076
cord_uid: 3h59uw8r

The ability to sequence genomes from ancient biological material has provided a rich source of information for evolutionary biology and engaged considerable public interest. Although most studies of ancient genomes have focused on vertebrates, particularly archaic humans, newer technologies allow the capture of microbial pathogens and microbiomes from ancient and historical human and non-human remains. This coming of age has been made possible by techniques that allow the preferential capture and amplification of discrete genomes from a background of predominantly host and environmental DNA. There are now near-complete ancient genome sequences for three pathogens of considerable historical interest — pre-modern bubonic plague (Yersinia pestis), smallpox (Variola virus) and cholera (Vibrio cholerae) — and for three equally important endemic human disease agents — Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy) and Treponema pallidum pallidum (syphilis). Genomic data from these pathogens have extended earlier work by paleopathologists. There have been efforts to sequence the genomes of additional ancient pathogens, with the potential to broaden our understanding of the infectious disease burden common to past populations from the Bronze Age to the early 20(th) century. In this review we describe the state-of-the-art of this rapidly developing field, highlight the contributions of ancient pathogen genomics to multidisciplinary endeavors and describe some of the limitations in resolving questions about the emergence and long-term evolution of pathogens.

The evolutionary and epidemiological history of human pathogens reflects a combination of short-and long-term interactions among biological, ecological and cultural factors. Although these processes have been studied using a variety of approaches, work over the past 20 years has shown how studies of ancient DNA and ancient RNA can provide important insights into the origins and history of infectious disease. For example, ancient DNA has led to a new understanding of endemic human diseases that span millennia, such as malaria, leprosy and tuberculosis [1] [2] [3] , as well as sporadic and devastating pandemic diseases, such as the recurrences of plague across Eurasia [4] . Reconstructing the evolutionary history of these pathogens with the help of ancient DNA has illuminated novel aspects of their evolution and epidemiology and has provided unexpected insights into the nature of disease-associated microbes [5] [6] [7] [8] [9] [10] [11] [12] [13] .

By the early 1990s, ancient DNA work on tuberculosis and leprosy had stoked optimism that genomic techniques could provide informative new data on the origin and evolution of infections still common today. Nevertheless, most human pathogens remained undetected in ancient samples until novel highthroughput sequencing and targeted enrichment offered enhanced resolution. For example, this genomic-scale work has led to revised narratives on the temporal-spatial emergence and transmission of M. tuberculosis and M. leprae, generated indisputable evidence of the bacterium responsible for yaws (Treponema palladium pertenue) being present in Western Europe before the journeys of Columbus and revealed the global dispersal of ubiquitous human pathogens, including hepatitis B virus (Box 1) and Helicobacter pylori. In combination, this work will continue to revise our understanding of the history of infectious disease [4, [13] [14] [15] .

The causes of differential mortality within populations had previously received little more than rudimentary attention in most histories of pre-modern disease. Ancient DNA provided a new perspective, by showing endemic infectious diseases increased the risk of dying during times of major epidemics such as the Black Death [16, 17] . For instance, Salmonella enterica contributed to catastrophic population collapse in Mesoamerica during the 16 th century European conquest, but also was part of the background disease landscape in 13 th century Norway [11, 18, 19] . Similarly, pathogens fortuitously retrieved from individual skeletal remains provide important evidence related to intervals of known, but poorly understood, mortality crises. For example, historical evidence of the toll that typically survivable infections could cause in crisis circumstances is illustrated by the Great Irish Famine (Box 2). Incidental and unexpected findings are also becoming increasingly routine with ancient pathogen genomics, enriching our biological and cultural understanding of the past. As a case in point, some of the best and more surprising sources of pathogen ancient DNA are calcified abscesses and kidney, urinary and placental stones [12, 20] . The analysis of vaginal flora isolated from a calcified abscess in the skeletal remains of a female from 13 th century CE Troy reminds us that maternal sepsis was, and still is, a common cause of death among adult females.

Despite the evident advances stemming from ancient DNA, it is also clear that increased efforts to include additional disciplines (e.g., historians and ethicists) and stakeholders (e.g., descendants, relatives or community relations and museums) are needed for the true value of ancient pathogen genomes to be properly contextualized [21] . In particular, an immediate challenge that we will discuss below is how to place genetic evidence into a more inclusive and expansive historical context, augmented by climatological, evolutionary and archaeological frameworks [22] .

A new understanding of past human disease experience is built on the preservation and extraction of unstable and damageprone nucleic acids. Ancient DNA and ancient RNA are relative terms used to characterize the damaged nucleic acids isolated from variously aged samples. DNA damage begins immediately post mortem during necrosis, via host enzymes and microbiomes, as well as via chemical pathways such as hydrolysis, alkylation, condensation and oxidation [23, 24] . Once the body is buried, stored or treated (the rate of which matters in each case), further damage to DNA occurs via exogenous microbial and fungal contaminants which leach into the remains over time ( Figure 1 ). More broadly, the conditions under which DNA is degraded or preserved can vary dramatically by the regional geology, climate and weather, by proximate environmental conditions (cave versus open-air sites, burial depth and humidity) and at the microscopic level (skeletal elements, mineral content and exposure) [23] [24] [25] . Given the tremendous variation in preservation potential, there is no simple correlation between the level of DNA damage and the age of a sample or its origin. Abiotic environmental constraints limit access to ancient DNA from historically relevant contexts (Africa, the Middle East and Central and South America), but exceptions, such as cave sites, may prove successful.

Ancient DNA molecules are short, usually between 30 to 60 base pairs long, the result of hydrolytic damage to both the phosphodiester bond (the weakest link in DNA) and the glycosidic bond [26, 27] , as well as oxidative damage to the sugar rings and pyrimidines [27, 28] . Genome stretches break into increasingly shorter DNA molecules, in a manner that is largely dependent on the availability of free water and the movement of reactants [23] . Double-stranded fragments often contain short single-stranded overhangs that are more susceptible to modification. Ancient DNA termini contain modified phosphates that are likely to be inaccessible to current library-preparation methods [28] . In addition, amine-containing nucleobases, such as cytosines (and adenine and guanine to a lesser extent), found in these single-stranded overhangs are subject to hydrolytic In the 1950s, Baruch S. Blumberg, who was interested in human variation in susceptibility and resistance to disease, began a global-scale collection of human serum. He identified a specific antigen present in the blood of hepatitis victims, called the ''Australia antigen'' because it was also present in an indigenous Australian. In the 1970s, David S. Dane identified the presence of virus-like particles with antigen-positive serum, leading to the isolation of hepatitis B virus (HBV) and now the identification of 10 genotypes (A through J), with geographic distributions differing among genomes (for example, some genotypes are restricted to the Americas). Until recently the historical evidence for HBV infection in humans spanned less than 100 years, although the global distribution and high prevalence of the virus in many populations suggested that it might have been associated with humans since their dispersal out of Africa [125] . Molecular clock analyses of contemporary genomes were unable to paint a clear picture of HBV origins, nor the timescale over which it has occurred. However, the integration of ancient DNA shows that the virus evolves slowly in the long term: for example, there was no evidence for temporal signal in the data even with the addition of ancient genomes recovered from 16 th century child mummies from Italy and Korea [9, 126] . Consequently, some researchers have attempted to calibrate the evolutionary timescale of HBV by assuming codivergence of the viral genotypes with their human host populations [111, 127] . Analyses of HBV genomes from the two child mummies showed that the genomes were nested within the contemporary diversity of the virus in these locations, indicating that the present genotype distribution is at least 400 years old. However, some ancient HBV genomes reconstructed from Bronze Age (5000 years ago) and Iron Age (1500 years ago) human remains from Eastern and Western Europe are actually closer to lineages circulating in chimpanzees and gorillas today [128, 129] (Figure 3 ). Although the extent of temporal signal is still uncertain even with the inclusion of the Bronze Age genomes ( Figure 3C ), the phylogenetic position of the ancient samples indicates that HBV has had a long association with human populations, spanning at least seven thousand years, and that its present distribution can largely be explained by normal and/or forced human migration [130] . Nevertheless, the exact patterns of migration, including when and how the virus arrived in the Americas, remain unclear. deamination, which results in uracil. Uracils that are subsequently copied as thymines in amplification reactions lead to C-to-T and G-to-A miscoding lesions, the tell-tale sign of ancient DNA authenticity [24, 29, 30] .

Ancient DNA isolated from skeletal remains is a complex mix of endogenous (target) with exogenous (non-target) DNAs that have worked their way into the remains over their burial history. As endogenous DNA of the sample degrades and leaches into the sedimentary environment, exogenous DNAs from the environment do the opposite. The equilibrium between these DNAs is guided mostly by the availability and movement of free water, bacteria, biofilms and fungal hyphae. Finally, exogenous DNAs can infiltrate the sample or extracts during or after excavation, handling, storage and laboratory-based methodologies or activities [31] . Between 60% and 90% of ancient DNA is lost during extraction, library preparation and the numerous purification steps that render the extracts free of their often extensive, enzyme-inhibiting coeluates [24, 32] . The target constituent of ancient DNA extracts typically ranges between 1% and 10%, with a detectable pathogen fraction usually far less than 0.5% of the total metagenome [33] [34] [35] , although exceptions have been found in calcified nodules, ear ossicles and permafrost samples [12, 20, 36] .

The even faster degradation of RNA molecules places limitations on evolutionary studies of pathogens with RNA genomes, such as influenza virus, measles virus, HIV and yellow fever virus. The 2' hydroxyl group of RNA increases the strength of the glycosidic bond, slowing depurination but weakening the phosphodiester bond. The nearly ubiquitous presence of RNases explains the further rapid degradation of RNA. Little work has been done on the preservation and degradation of RNA from ancient samples, although recent studies suggest that RNA can survive under favorable conditions for millennia [37, 38] . On the other hand, formalin fixation, particularly problematic for DNA, appears to preserve RNA and this has led to exciting recoveries of RNA viruses from stored archival tissues [39] [40] [41] [42] .

Pathogen Fortitude Differential preservation of pathogens plays a key role in their survival and thus their direct association with some historically important epidemics and infections, but not others. Bacterial endospores fall on one end of the spectrum, with species such as Bacillus anthracis and Clostridium difficile the most suited for long-term DNA preservation [43] . Endospores have been revived after decades of dormancy, yet they have not been detected nor are they being sought from ancient remains.

Most bacteria have a single wall and membrane composed entirely of peptidoglycan. Gram-positive bacteria have a twoto ten-fold thicker cell wall than Gram-negative bacteria, which may play a role in the longer preservation of their DNA. This is unfortunate, as Gram-negative bacteria account for the bulk of pathogenic bacteria. Mycobacterium tuberculosis and M. leprae (both Gram-positive) produce thick lipid exudates or nodules that are likely to help repel water and enhance preservation, which may account for the reduced signal of deamination that was seen in an ancient M. leprae genome [44] .

Treponema pallidum pallidum and pertenue (causative agents of syphilis and yaws, respectively) sit at the other end of the spectrum because they lack an outer cell membrane altogether. This feature, coupled with the low pathogenic load of Treponema in late stages of infection [45] , largely precluded its identification in ancient remains [45] [46] [47] until the recent application of highthroughput sequencing methods [15, 17, 48] .

Some viruses have been shown to survive and still be infectious after freeze-drying and being at room temperature for decades, but how this correlates with long-term preservation of their nucleic acids in the geological record remains unknown [49] . Both Gram-negative bacteria and RNA viruses are The plant pathogen Phytophthora infestans, a eukaryotic oomycete that causes potato late blight, arrived in Europe by 1845 and destroyed potato crops across much of the continent [131] . The resultant Great Irish Famine became a demographic catastrophe because a quarter of Ireland's people subsisted on potatoes and a great many others relied on the potato as a primary staple. Over a million persons emigrated, leaving the net demographic decline near 30% of Ireland's population (worse than what Ireland suffered during the Black Death), though relatively few people died from outright starvation. Demographic and epidemiological analysis of the causes of death in this subsistence crisis [132] shows patterns common in historical mortality crises. Cause-of-death reports show that many different infectious diseases circulated and were made more lethal by progressive undernutrition and malnutrition. Still-able hunger victims migrated to towns and cities in search of food or work. Tickborne typhus (Rickettsia prowaszecki) and relapsing fever (Borrelia recurrentis) in particular added to the poorly defined 'famine fevers' that were the leading cause of death and illness. Deaths from diarrhea and dysentery during these years were amplified by the second wave of pandemic cholera, making it difficult to distinguish the burden of each disease. Transatlantic outbreaks of virulent streptococcal infections might have contributed to deaths from 'rheumatism', but would not have been distinctive among rashes from scurvy or the various epidemic fevers. Botanists at the time described the pathogen and its effects on potato plants. Ancient genomic DNA recovered from preserved leaves in several different European herbaria illustrates some of the uses and limitations of recent research [133] [134] [135] [136] [137] . Phylogenetic analyses of complete nuclear and mitochondrial genomes from P. infestans have painted a complex picture of host switching and hybridization with P. andina, an emerging pathogen of solanaceous crops [135] . The genomic evidence supports an origin for P. infestans in Mexico or the South American Andes [134] [135] [136] [137] . The HERB-1 lineage of the pathogen, which was responsible for the Irish potato famine, had either multiple introductions into Europe or a single introduction of multiple haplotypes in the 19th century [135] . Although the HERB-1 lineage was believed to have gone extinct [137] , analyses of strains from the Americas have revealed that it persists today in Central and South America [134] .

particularly prone to rapid degradation, yet are the causative agents of many infectious diseases; they remain underrepresented in current ancient DNA studies. Further clarification on how the physical characteristics of some bacteria and viruses act to enhance or detract from nucleic acid preservation, along with improvements to ancient DNA methodologies, may help in their further targeting, detection and reconstruction.

As the individual response to disease is complex [14, 50, 51] , the presence of pathological lesions can be an unreliable way to sample ancient pathogen DNA. It may be better to sample from the population at large, or from a locality where an infectious agent is presumed present, such as a hospital site or mass grave. Only a handful of ancient pathogen DNA studies have relied on skeletal indicators as a sampling strategy and almost all revolve around where skeletal evidence strongly, but not solely, suggests a diagnosis of tuberculosis or leprosy [1, 2, 52] .

In general, ancient pathogen DNA detection strongly depends on many important factors: host susceptibility to a particular infection; level of pathogen in the blood or tissue (i.e., bacteremia, viremia or parasitemia); differential immune responses; presence and level of infection at time of death; nutritional state (result of famine or social status); and prior or ongoing co-infection or co-morbidity. Although a presumed pathogen may be associated with a well-dated burial assemblage (e.g., Y. pestis at East Smithfield cemetery, London, UK), the level of infection at death and the burial conditions will impact the relative abundance and detectability of ancient pathogen DNA within each skeletal element. Because yields of pathogen ancient DNA can be low even in optimal sites for recovery, there are ethical ramifications that need to be carefully considered prior to beginning the work [21] . Relatively little attention has thus far been paid to the sheer bulk of finite, irreplaceable resources required in obtaining the genome sequences of ancient pathogens.

The archaeological record comprises, in decreasing order of frequency, skeletal elements (bones and teeth), mummified soft tissue, hair and human-associated trace fossils (e.g., faeces, abscesses and associated sediment). The historical period includes archival samples stored in medical museums [42, 53, 54] . Microbial DNA has now been isolated from almost all of these sources [55] . The successful recovery of pathogen ancient DNA depends on taking into account disease pathogenesis (i.e., localization or dissemination of microbes throughout infection and immunity responses) to guide selection of the most appropriate sampling strategy.

Pathogens that cause infections confined to soft tissue usually cannot be recovered from skeletal material, as they cause varied Degradation and preservation of host DNA (including the microbiome and circulating pathogens) occurs rapidly after death. DNA degrades in an environmentally informed manner, where water, temperature, pH, microbial soil content and other factors shorten the remaining endogenous DNA but most is lost. At the same time, environmental (exogenous) DNA swamps out the original signal, with fragment-length distributions that both overlap and extend beyond the ancient DNA in the sample. Medieval peoples were well acquainted with death as the 15 th century manuscript illustration on the left suggests (image: ª British Library Board: Add MS 37049). Here, a recently buried noblewoman (endogenous) sits atop the worms and other vermin (exogenous), who will soon set upon her corpse. The image accompanies a Middle English poem relating a discussion between the decaying woman and eager worms [138] .

clinical manifestations that primarily affect the intestinal system (e.g., inflammation, vomiting, fever or diarrhea) [56] . For these pathogens, acquiring ancient DNA requires soft-tissue samples, faecal remains, gut-associated sediment samples or archival medical or histologically preserved samples [57] . In exceptional cases, enteric pathogens can become blood borne, which might explain the recovery of Salmonella enterica serovar Paratyphi C (causing bacterial enteric fever) from teeth of victims from a 16 th century Mesoamerican epidemic [18, 58] as well as from 13 th century Norway [11, 19] . Hematogenous infections, such as bubonic plague, are more likely to be detectable in the blood-fed root or pulp cavity of teeth than in bones. Commensal microorganisms, which form part of the human microbiome, inhabit diverse niches within the human body and can become opportunistic pathogens [59] . Reconstructing commensal pathogen genomes and our ancient microbiomes will clearly provide a more nuanced comprehension of health status and disease susceptibility in the past. Recent work on reconstructing these metagenomes has largely focused on using the constituents of dental calculus (calcified dental plaque) from human remains. These studies have shown the shifting composition of the oral microbiota during major transitionssuch as from hunting and gathering to farming 10,000 years ago [60] -as well as the antiquity of putative antibiotic resistance genes associated with the oral microbiota [61] . The presence, preservation and ultimate detection of ancient pathogen DNA or RNA depends on many complicating factors that include the epidemiology of the disease, the biology of pathogen and its load at time of death. These factors coupled with the fickle nature of preservation make their isolation from ancient remains both surprising and fortuitous.

Until 2011, pathogen detection from archaeological remains was almost entirely based on polymerase chain reaction (PCR) [24] . While PCR benefits from high sensitivity (theoretically able to detect single molecules) and specificity, most amplicons are long (100-200 base pairs) and target a tiny fraction of the available DNA (Figure 1 ). When adding the complications of preservation and pathogenesis, it is unsurprising that many attempts at detecting pathogens have met with negative results. Authenticating PCR-based studies of pathogen DNA requires a number of experimental controls and replication [24] . Discriminating true human pathogens from closely related environmental relatives can be complicated; it is particularly problematic with the use of repetitive elements for the identification of Mycobacterium species, because these mobile elements are similar to those in other mycobacterial agents or extraneous soil microorganisms. Today PCR assays are used as a quick and easy but resourcecostly way to prioritize the best libraries from a large pool of potential samples for a more targeted or deep shotgun sequencing.

With the advent of high-throughput sequencing came the move from the simple identification of pathogens at a certain place and time to genome-level reconstruction and characterization, enabling a shift to hypothesis-driven research carried out within evolutionary frameworks [62] [63] [64] [65] . The methodological tweaks to DNA extractions and novel library preparations designed to target short, damaged and sometimes single-stranded DNA have enhanced access to the target component of ancient extracts (Box 3) [66, 67] . High-throughput sequencing methods require conversion (and therefore loss) to platform-specific DNA libraries. This can involve end-repair and polymerase fillin reactions, which will occur on only a small proportion of actual surviving and less damaged, repairable DNA molecules (those containing appropriate 5' and 3' termini) [68] . Single-stranded library preparations make the most theoretical sense because they mitigate some of the many problems associated with ancient damaged DNAs [67, 69] .

Pathogen-screening workflows include a shallow shotgunsequencing approach (i.e., 1-2 million DNA reads), followed by various bioinformatic pipelines to ascertain, at first, the likelihood of pathogen-specific DNA markers. There is a growing need for new algorithms to confirm the presence of a pathogen in ancient remains, given that individual strains of microorganism contain few genomic regions that are specific at taxonomic levels and that DNA is damaged and at low quantities. Various methods exist for taxonomic binning, sequence alignment and visualization, such as MEGAN [70] , HOPS [71] , BLAST [72] and KRAKEN [73] . For a good overview of the applicability of these and other methods to ancient DNA, see [74] .

As pathogen-specific endogenous molecules will usually be at extremely low abundance or undetectable at shallow sequencing depths, costs can quickly become prohibitive to most researchers. Shotgun sequencing as a means for the recovery of full genomes is only possible with exceptionally preserved specimens, where endogenous pathogen DNA content is high [12, 20] . An alternative way to identify the few pathogen reads in a DNA extract is through the use of pathogen-enrichment arrays designed to enrich short genome-specific regions, thereby increasing the proportion of the target molecules [35, 75, 76] . Probe sets are scalable (hundreds of kilobases to megabases), flexible and cost-effective; however, custom microarrays can obtain higher probe density (hundreds of thousands to millions of probes) to enable the capture of larger target regions [77] . Anecdotal evidence suggests that these methods differ considerably in efficiency, yet there remains tremendous room for optimization [78, 79] .

Current limitations of enrichment include biased database representation, leading to an inability to capture distantly related, extinct, divergent species/strains, novel genomic material, indels, plasmids, or exceptionally divergent regions. Ancestrally designed genomes (based on a phylogenetic tree) can be used to capture more distantly related taxa, those that are perceived to have diverged over longer time periods, or those that evolve more rapidly (RNA viruses).

In the future, advances in the recovery of nucleic acids from ancient remains and in the detection and capture of ancient pathogen DNA will involve reduction in costs and in the amount of material required. Such advances will also lead to improvements on success rates, which are currently below 0.01%.

Establishing the Presence of Specific Pathogens in Historical Populations Disease symptoms, especially when inferred from historic records or skeletal remains are rarely diagnostic and have long generated doubts and debates about the nature of past disease and epidemic experience ( Table 1 ). Demonstrating that a given ll Current Biology 30, R1215-R1231, October 5, 2020 R1219 Review microbe is present at a particular time and place can help to include that pathogen as a possible driver of a past mortality event. Characterizing pathogens genetically can also identify epidemic-associated, but clinically unconfirmed, opportunistic infections [12] . A dramatic illustration of the new ability to identify ancient DNA in a mass burial site came with the isolation of Y. pestis in a cemetery precisely dated to the Black Death (between the 1340s and 1350s) [5] , generating new questions about the origins, persistence and transmission of recurrent plagues over subsequent centuries [80] . Similarly, isolates from mass burial sites related to the Roman Justinianic plague cycle (540-750) assisted the identification of the pathogen associated with these important historical outbreaks [6, [80] [81] [82] . Finally, plague ancient DNA recovered from western European sites of the 16 th through 18 th centuries showed the persistence of the initial Black Death strain of Y. pestis over several centuries [7, 83, 84] . Hence, the recovery and firm identification of ancient Y. pestis helped focus debates about the causative agent of these outbreaks that had developed due to the inevitable ambiguities of historical accounts.

Other notable ancient DNA findings linked to major mortality events include isolation of Salmonella enterica in association with cocoliztli, an indigenous Aztec term used to distinguish catastrophic die-offs from 1545 to the 1570s from the initial epidemics attributed to smallpox or measles [25] . A reconstructed strain of Vibrio cholerae, classical subtype, provided a more expected ancient DNA recovery from the second global pandemic of cholera in 1849 [25] . Important insights into the evolution of recent viral pandemics have been provided by influenza A/H1N1 virus recovered from paraffin-imbedded lung tissue of soldiers who died in the 1918/1919 influenza pandemic [25] and from strains of HIV subtype D from African victims during the early stages (1950s-1960s) of what would become the global HIV/AIDS epidemic [39, 41] .

Evidence that links a given pathogen to a specific historical context is rare, particularly as multiple pathogens can co-circulate in populations. Indeed, attention to co-circulating pathogens in a given population or environment may provide a more complete understanding of past epidemics than does focussing on individual pathogenic causes [85] . Such an approach, known as 'pathocoenosis', has been applied to the study of ancient Plasmodium parasites [86, 87] .

Outbreaks of a specific disease can vary enormously in the extent of associated disease morbidity and mortality. Comparison of the genomes of microbes associated with very different case fatality rates provides a unique opportunity to explore whether these differences are due to genetically encoded virulence traits or environmental or host factors (or both), the presence or absence of key genetic features (e.g., antimicrobial resistance) and how virulence determinants evolve [88] . At a genomic level, virulence can be driven by the presence of plasmids, mobile genetic elements or single-nucleotide polymorphisms [88, 89] .

Increased access to ancient microbial genomes will help place constraints on the origin and timing of certain alleles in the population and enable, in some cases, the direct testing of their effect on pathogen adaptation [89] . Work on Y. pestis has shown A range of methods can be employed to authenticate ancient pathogen sequences. These fall into eight categories: 1. DNA damage. Ancient DNA molecules are expected to be shorter (medians of 40-60 base pairs) than those from recent contaminating sources, although the distributions overlap. Current ancient DNA extractions and library preparations select for shorter and less damaged DNA molecules and thus the pool includes contaminating DNAs. Deamination patterns marked as C-to-T and G-to-A changes appear on the ends of authentic DNA molecules. 2. Random genome coverage. Mapped reads should be randomly and equally distributed across the genome (enriched at GC regions). Stacked, mapped reads are a tell-tale sign of multiple contributing sources of DNA. 3. Human DNA as benchmark. As most pathogen DNA is being reconstructed from the remains of human victims, using the presence and damage patterns of human DNA as a surrogate for their preservation is not without merit. 4. Mismatch distributions. Simple mismatch distributions, percent identities, or edit distances can be used to assess the proximity of the ancient reads to closely related contemporary genomes. This metric is database reliant. The expectation is that the majority of reads (barring damage) should be zero to a few mismatches from contemporary strains. 5. Haploidy of mapped reads. If sequencing depth is extensive, one can assess the overall contribution of singular or multiple circulating strains or even related pathogens through a comparison of the small number of variant sites using a site-frequency spectrum. 6. Comparison to contamination. There are ways to estimate the amount of contaminating modern human mitochondrial and nuclear DNA in ancient libraries [139] . Species identifications based on single-nucleotide polymorphisms or k-mers can help to distinguish soil-dwelling contaminants from real signals that can be further mined for genomic information. 7. Surrogates of microbial presence. Recovery of pathogen-specific protein or lipid markers from skeletal remains can support the identification of non-DNA products of pathogen infections (e.g., mycolic acids, proteins), immunological detection of antigens and/or the identification of organic biomolecules such as lipids [140] . 8. Phylogenetics. Phylogenetic analysis can determine the relationship between ancient genomes and the modern diversity of the pathogen and a molecular-clock analysis can potentially be used to validate the age of the ancient sample itself [9, 42] . Ancient genomes from lineages that are presumed to be now extinct will tend to sit outside the diversity of modern samples ( Figure 2 ) [10, 19] . In contrast, ancient genomes that are nested within the current diversity should be carefully assessed as possible contaminants. Paleopathological diagnosis presumed smallpox via pox-like rash on skin and electron microscopic images of viral particles as well as immunoassays from a 16 th century Italian mummy [141] . Contemporary with urban smallpox epidemics and catastrophic presumed smallpox in the Americas.

Discussion of early medical descriptions of smallpox epidemics in historical context, and the appearance of pocks on mummified human skin, were convincing evidence of VARV in the past. HBV, not VARV, found in a 16th century Italian mummy. The strain of HBV recovered belongs to genotype D, common in the Mediterranean region today [9] .

Partial mummified remains of a child from Lithuanian crypt, late 17 th century (image from [10] Review that it differs from Yersinia pseudotuberculosis, its closest modern relative, by deletions to the genome and the acquisition of two plasmids (pMT1 and pPCP1), which enabled its transmission from arthropods to mammalian hosts as bubonic plague.

The finding that a 2,900-year-old Y. pestis genome already had these adaptations allows us some resolution on the origin and timing of the virulence mechanisms that led to high human fatality rates [11, 90] . In contrast, ancient virulence factors found in a 

Calcified nodule, 12 th century Troy (image: Caitlin Pepperell).

Much perinatal maternal mortality in the past assumed to be due to acute infections with common skin pathogens, Streptococcus spp. and Staphylococcus spp.

Both S. saprophyticus and G. vaginalis found in an abcess resulting from ectopic calcification, 13th century Troy. Opportunistic pathogens in the past enhanced overall mortality and contributed to differential mortality experience documented in many historical societies [12] .

12 th century genome of S. enterica Paratyphi C contained the same pathogenicity islands in mobile genetic elements as do modern strains, suggesting that the bacterium has caused similar disease over its 3,000-year history [19] . With both S. enterica and Y. pestis, the ancient presence of virulence factors indicates that subsequent decline in their incidence and disease severity was most likely driven by changes in human behaviour and pathogen ecology, rather than a decrease in their ability to spread and cause disease.

One of the main goals of genomic and historical analyses of pathogens is to reconstruct the timescale of their emergence and spread. Human pathogens commonly emerge by jumping across a species barrier from animals (zoonoses) [91] , although most of these jumps result in dead-end emergences and extinctions as was the case with the strain of Y. pestis responsible for the Plague of Justinian. A different strain of Y. pestis affected all areas of western Eurasia and north Africa during the second pandemic of plague, between the 14 th and 19 th centuries from regionally diverse rodent host foci. Occasionally, however, newly emerged pathogens establish sustained transmission cycles. Estimating the timing of the initial host switch, while interesting from a historical perspective, is challenging because this event could have happened at any point in the window of time between the divergence of the pathogen from its sister lineage and the most recent common ancestor of the sampled genomes of the pathogen (Figure 2 ) [92] . The bounds of this window can be estimated using molecular-clock analyses of genomic sequence data but can be very uncertain, which limits specific historical contextualization.

Reconstructing the evolutionary timescale of a pathogen relies on accurate estimation of the rate of its evolution. Pathogen genomes can undergo measurable amounts of evolutionary change between sampling times, which can span a period of weeks, years or even millennia [93, 94] . Adding ancient genomes to data sets widens the sampling period, allowing a greater amount of evolutionary change to be evaluated ( Figure 3 ) [95, 96] . The assumption behind this approach is that evolutionary change occurs at a predictable rate over time, known as the 'molecular clock'. Unfortunately, pathogens with complicated epidemiology might evolve at highly variable rates and thus need careful evaluation before their genomes are used in evolutionary analyses, if molecular dating is to be attempted at all. An important step in the phylogenomic analysis of pathogens is testing whether the genomic data carry a temporal signal. The simplest way to do this is to estimate a phylogenetic tree with branch lengths that are proportional to the observed number of substitutions along that branch and then to conduct a regression of the distance from the root to each of the tips as a function of the sampling times [97] . A positive slope suggests that the data set has a temporal signal (Figure 3 ), given that we expect samples with ancient DNA to have a smaller distance from the root than those with modern DNA (variola virus genomes show this nicely, for example). More sophisticated methods for evaluating the temporal signal and for estimating evolutionary rates have been developed over the past two decades [98] . The flexibility of these phylogenetic methods also allows the dating of ancient DNA samples for which the true age might be unknown [99] , or to validate estimates of sample ages obtained using independent techniques [9, 27] . Schematic phylogenetic tree showing how ancient genomes can provide information on key various aspects of the evolutionary history of pathogens. An ancient pathogen genome can potentially be placed as (A) an extinct sister lineage to the modern diversity (e.g., ancient variola virus from 7 th -10 th century [103] ; image: teeth from an East Smithfield individual used to indicate source from which ancient variola-like sequences have been isolated, courtesy of Sharon DeWitte, Museum of London); (B) an extinct sister lineage to a modern haplogroup or genotype, but still falling within the modern clade of the pathogen (e.g., variola virus from the 17 th century; image: VD21 child mummy from Vilnius, Lithuania, from which smallpox was detected, from [10] ); (C) belonging to a present-day haplogroup or genotype (e.g., hepatitis B virus from the 16 th century [9, 125] ; image: child mummy from Naples, Italy, with HBV detected, ª 2018 Patterson Ross et al. CC BY 4.0); or (D) at the base of the modern clade, possibly as a direct ancestor (e.g., Y. pestis from 1348 [5] ; image: skull of an individual from East Smithfield that yielded a Black Death genome, courtesy of Jelena Bekvalac, Museum of London (MIN86)). The square symbol denotes the common ancestor of modern pathogen samples. Importantly, the emergence of a pathogen in its present host could have occurred at any point along the branch between the divergence from its closest extant or extinct relative and the most recent common ancestor of the sampled isolates (branch in blue). Ancient pathogen genomes can help narrow this window of emergence [90] while also shedding light on any gain or loss of function along the branches leading up to the modern clade (branch in orange).

The main determinants of evolutionary rates are the type of nucleic acid (i.e., RNA or DNA), the genome size and the host. Viruses with single-stranded RNA genomes usually evolve orders of magnitude more rapidly than those with double-stranded DNA genomes [100] . The upshot of such evolutionary rate variation in microbial genomes is that the sampling period should be carefully considered according to the pathogen being studied. Sequence data collected from rapidly evolving RNA viruses over a few years are usually sufficiently informative to allow their evolutionary rate to be estimated. In contrast, for more slowly evolving bacteria such as M. tuberculosis, an ancient genome is sometimes essential to obtain sufficient temporal signal to estimate an evolutionary rate and timescale [101, 102] .

The divergence between a pathogen and its closest known relative places a maximum bound on the timing of its emergence. For example, sequences of Variola virus recovered from human remains from the Early Middle Ages (793-1066) fall outside the diversity of the lineages associated with modern smallpox and are separated from them by a long branch [103] . This phylogenetic pattern indicates that the initial cross-species transmission to humans must have occurred between the earliest detection of this virus -dating back to at least the Viking era -and its closest known animal relatives, camelpox from camels and taterapox from African gerbils [10] (Figure 2) . Importantly, the 'Viking' age and modern smallpox Variola strains differ in the inactivated genes they possess. This feature suggests that smallpox, as we understand it epidemiologically today, is likely to have evolved at some point after the divergence of the Viking age Variola strains and the common ancestor of modern smallpox. Unfortunately, estimating this date can be seriously misled by the effects of purifying selection and misspecification of the evolutionary model. In particular, the removal of deleterious genetic variants by purifying selection can produce a timedependent bias in evolutionary rate estimates, leading to underestimation of deep divergence times [102, 104] . Adding ancient genomes to the data set widens the sampling period, potentially allowing for better inference of the evolutionary rate of the pathogen. However, if the ancient genomes represent divergent lineages that differ in their biological properties, the inclusion of these genomes can also introduce considerable heterogeneity that may require careful consideration of the molecular-clock model and abandoning hope of obtaining meaningful date estimates [105, 106] .

Ancient genomes are often the only way to uncover the existence of now-extinct pathogen lineages. For some pathogens, the most recent common ancestor of modern lineages is surprisingly young compared with perceptions of the age of the disease. This is especially true for pathogens that have experienced substantial reductions in genetic diversity, such as viruses that have been targeted by extensive vaccination programs [10, 42] Review analyses can still be useful in placing bounds on when pathogens emerged and diversified. In this way, molecular dating provides a valuable means of testing hypotheses about pathogen origins and evolution.

The ancestral movements and present-day distributions of pathogens can be inferred by phylogeographic analysis, which combines geographic and evolutionary information. Simple phylogeographic methods overlay geographic information on the tips of phylogenetic trees to interpret how migration events could have given rise to the observed distribution. Formal statistical approaches are used to model past movements between geographic locations, which can be specified as a discrete (e.g., country or city) or a continuous trait (e.g., GPS coordinates) [107] . Importantly, phylogeographic models are sensitive to sampling bias; accurate inference of ancestral locations typically requires the inclusion of samples from the ancestral location itself [108] , which is certainly unknown in the ancient context described herein. One of the benefits of ancient DNA is to provide evidence for or against the presumed ancestral origin of a particular strain or pathogen (Box 1). In some cases, ancient pathogen phylogeography has been used as a proxy for ancient human dispersal [109, 110] , sometimes providing results contrary to expectations defined by archaeological evidence alone [111] .

The spatial distribution and emergence of pathogens that spread rapidly in host populations tend to have complex geographic histories that would be impossible to untangle without ancient DNA. How Y. pestis spread through rodent populations and caused repeated human spill-over epidemics in Europe and the Middle East between the 1340s and the 1840s has received substantial attention. As more than 100 ancient genomes of Y. pestis genomes have now been sequenced, it may be possible to inform historical debates through analyses of ancient DNA. For example, did European plague epidemics originate from repeated re-emergences from local rodent reservoirs [112] or did successive transmission waves reach Europe, initiated by great gerbil epizootics in the central Eurasian grasslands [113] ? The discovery of Y. pestis genomes in mass burial sites in 18 th century France revealed a lineage that is not present in the modern diversity of this bacterium. Instead, the lineage is linked to the 14 th century one from England and thus supports the hypothesis of a presently extinct Mediterranean reservoir rather than the traditional narrative, that plague always 'came from the east' [84, 85] .

Many pathogens evolve at timescales similar to those over which they accumulate genetic changes, and their phylogenetic patterns provide meaningful insights into their epidemiological processes. These approaches have been combined in the field of phylodynamics [114] , where various branching models have been proposed to infer epidemiological parameters. Most of their models are based on the coalescent, which models the ancestry of genes sampled from a population, or the birthdeath process, which models speciation and extinction through time [115] . The inclusion of ancient DNA data brings an opportunity to infer the tempo and mode of previous pandemics, but it requires careful consideration of the model in question and its assumptions.

Coalescent and birth-death models with constant dynamics perform well in analyses of recently emerging outbreaks, but they are not realistic for the complex epidemiological dynamics of infectious diseases with long histories. For example, the persistent re-emergence of Y. pestis would not conform to such constant dynamics. A more flexible alternative is offered by skyline methods that can allow the epidemic growth rate, population size and other parameters to vary across different time intervals [116] . Recent developments in skyline techniques that can explicitly model the sampling process [117] and those that allow the inclusion of occurrence data with low-quality or no sequence information [118] offer promising avenues for coherent phylodynamic analyses of ancient DNA. Further research into the performance of these methods is needed to assess their accuracy in reconstructing past epidemics.

Ethics and the Hunt for Ancient Pathogen DNA Despite the insights provided by ancient pathogen DNA, it is important to ask what the inherent value of an ancient pathogen genome is when it comes at the cost of partial destruction of a skeletal feature. As of this writing, 236 pathogen genomes have been reconstructed from 12,733 individual remains sampled, a success rate of 0.01%. Does this represent a reasonable rate of success, thereby justifying the extensive sampling? If not, what does and how is this value to be determined? While this ethical quandary is nothing new for the field of ancient DNA, nor archaeology more generally, it remains an important and unresolved issue. Curators rightly ask to see the results of a pilot project to determine the likelihood of overall success but measuring success can be elusive. Is it the isolation of human DNA from a certain percentage of human remains and, if so, what percentage is deemed a success in a pilot study and opens the way to broader sampling? Or do researchers need to find the presence of a notable pathogen prior to additional sampling? Neither of these guarantees the presence or retrievability of actual ancient pathogen genomes. Most ancient pathogen DNAs detected have come from scanning the 'left-over', non-human DNA reads of deeply sequenced libraries, which avoids additional destructive sampling but comes at a high financial cost.

With an increasing number of researchers in the field, competing for access to the same collections, the question of who gets to sample has increasing political ramifications. Should sampling permission be given to researchers who work (or propose to do so) on ancient human DNA, providing the raw sequencing data for others interested in pathogen research and vice versa, or only to research groups that are doing both? One thing that is clear is that before sampling can begin, stakeholders in any particular project or present-day descendants of the human remains being studied must agree on best practices, what and how much to sample and how best to integrate the DNA findings with existing written or oral histories. We also need to shed outdated historical narratives and the patterns of resource collection and storage that once drove Western colonial medicine and science. Fortunately, there is a strong precedent of how this can work, as illustrated by the SING consortium that trains indigenous scientists in genomics while bringing the indigenous perspective to the interpretation of the data [119] . However, unlike work for ancient human DNA, most pathogen DNA work remains mysterious as pathogen discoveries have come almost entirely as surprise findings, making it difficult but presumably necessary to justify destructive sampling for hypothesis-driven projects.

The recovery of ancient DNA has proven valuable for our understanding of how and when pathogens emerge and evolve. Early on, ancient DNA research relied on PCR and was used to query the dispersals of humans and other vertebrate species, adding some investigation of presumed tuberculosis or leprosy seen in skeletal remains. By the end of the 1990s, ancient DNA research had shifted to the phylogenetic relationships of particular pathogens prominent since the Neolithic [11, 120] . The recent sequencing of ancient RNA viral genomes [42] is an exciting development that may revolutionize our understanding of the emergence and origin of many other RNA-based pathogens. Moreover, newer analytical tools to query huge data sets from the modern era can now provide a more comprehensive understanding of ancient pathogen diversity. These advances were unimaginable at the beginning of the current millennium but now challenge the field to consider how best to conserve some part of our collective past resources for future generations. In the future, to better understand and contextualize the dynamics that drove pathogen evolution, it will be essential to improve collaboration with scholars who study the remote to the recent past, including palaeoclimatologists, archaeologists and historians.

The past can never be fully retrieved, no matter how advanced our technologies. For example, despite the generation of unprecedented amounts of genome sequence data during the COVID-19 pandemic (75,000 genomes in the GISAID database [121] ) and a detailed understanding of the spread of its causative agent, the SARS-CoV-2 virus [122] , the evolutionary history of the virus remains mysterious because of sparse sampling, multiple possible reservoirs and a complex history of recombination [123, 124] . Rewind a few centuries and pathogen identification, origins and evolution become clouded even further. Fortunately, by accessing the genomes of ancient pathogens that were responsible for historical morbidity and mortality events, ancient DNA can provide a glimpse into the evolutionary and epidemiological past, amplifying our understanding of infectious disease emergence and thereby placing pandemics like COVID-19 into their broader context. 

Ancient genomes reveal a high diversity of Mycobacterium leprae in medieval Europe

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis

Plasmodium falciparum malaria in 1st-2nd century CE southern Italy

The globalisations of disease

A draft genome of Yersinia pestis from victims of the Black Death

Yersinia pestis and the Plague of Justinian 541-543 AD: a genomic analysis

Historical Y. pestis genomes reveal the European Black Death as the source of ancient and modern plague pandemics

Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans

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

17th century variola virus reveals the recent history of smallpox

Emergence of human-adapted Salmonella enterica is linked to the Neolithization process

A molecular portrait of maternal sepsis from Byzantine Troy

The 5300-year-old Helicobacter pylori genome of the Iceman

Palaeopathology and genes: investigating the genetics of infectious diseases in excavated human skeletal remains and mummies from past populations

Ancient bacterial genomes reveal a formerly unknown diversity of Treponema pallidum strains in early modern Europe

Stature and frailty during the Black Death: the effect of stature on risks of epidemic mortality in London, A.D. 1348-1350

A treponemal genome from an historic plague victim supports a recent emergence of yaws and its presence in 15th century

Salmonella enterica genomes from victims of a major sixteenth-century epidemic in Mexico

Pan-genome analysis of ancient and modern Salmonella enterica demonstrates genomic stability of the invasive Para C lineage for millennia

Recovery of a medieval Brucella melitensis genome using shotgun metagenomics

To curate the molecular past, museums need a carefully considered set of best practices

History meets palaeoscience: consilience and collaboration in studying past societal responses to environmental change

Molecular preservation

Genetic analyses from ancient DNA

Comparing ancient DNA preservation in petrous bone and tooth cementum

Rate of depurination of native deoxyribonucleic acid

Instability and decay of the primary structure of DNA

Ancient DNA: extraction, characterization, molecular cloning and enzymatic amplification

DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA

mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters

Reagent and laboratory contamination can critically impact sequencebased microbiome analyses

On the circumstances surrounding the preservation and analysis of very old DNA

Timing of human protein evolution as revealed by massively parallel capture of Neandertal nuclear DNA sequences

Pulling out the 1%: whole-genome capture for the targeted enrichment of ancient DNA sequencing libraries

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

Human auditory ossicles as an alternative optimal source of ancient DNA

Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch

Ancient RNA from Late Pleistocene permafrost and historical canids shows tissue-specific transcriptome survival

An African HIV-1 sequence from 1959 and implications for the origin of the epidemic

Initial genetic characterization of the 1918 ''Spanish'' influenza virus

A near full-length HIV-1 genome from 1966 recovered from formalin-fixed paraffin-embedded tissue

Measles virus and rinderpest virus divergence dated to the sixth century

Spore resistance properties

Genome-wide comparison of medieval and modern Mycobacterium leprae

Secondary syphilis in Cali, Colombia: new concepts in disease pathogenesis

Identification of Treponema pallidum subspecies pallidum in a 200-year-old skeletal specimen

Digging deeper into the limits of ancient DNA research on syphilis

Historic Treponema pallidum genomes from Colonial Mexico retrieved from archaeological remains

Virion structure and composition

The osteological paradox 20 years later: past perspectives, future directions

The bioarchaeology of health crisis: infectious disease in the past

The bioarchaeology of leprosy and tuberculosis

Plasmodium vivax malaria viewed through the lens of an eradicated European Strain

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

Ancient pathogen genomics as an emerging tool for infectious disease research

Enteric pathogens through life stages

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

Salmonella: a review on pathogenesis, epidemiology and antibiotic resistance

Colonizing opportunistic pathogens (COPs): the beasts in all of us

Sequencing ancient calcified dental plaque shows changes in oral microbiota with dietary shifts of the Neolithic and Industrial revolutions

Pathogens and host immunity in the ancient human oral cavity

Ancient pathogens through human history: a paleogenomic perspective

Metagenomics to paleogenomics: large-scale sequencing of mammoth DNA

Next generation sequencing of ancient DNA: requirements, strategies and perspectives

New developments in ancient genomics

Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments

Single-stranded DNA library preparation from highly degraded DNA using T4 DNA ligase

Massive influence of DNA isolation and library preparation approaches on palaeogenomic sequencing data

A ligation-based singlestranded library preparation method to analyze cell-free DNA and synthetic oligos

MEGAN analysis of metagenomic data

HOPS: automated detection and authentication of pathogen DNA in archaeological remains

Basic local alignment search tool

Kraken: ultrafast metagenomic sequence classification using exact alignments

A robust framework for microbial archaeology

Enhanced virome sequencing using targeted sequence capture

Parallel detection of ancient pathogens via array-based DNA capture

Reconstructing ancient genomes and epigenomes

Comparative performance of two wholegenome capture methodologies on ancient DNA Illumina libraries

Experimental conditions improving in-solution target enrichment for ancient DNA

A high-coverage Yersinia pestis genome from a sixth-century Justinianic plague victim

Ancient Yersinia pestis genomes from across Western Europe reveal early diversification during the First Pandemic (541-750)

Yersinia pestis DNA from skeletal remains from the 6th century AD reveals insights into Justinianic plague

Phylogeography of the second plague pandemic revealed through analysis of historical Yersinia pestis genomes

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

Paleoparasitology and paleopathology. Synergies for reconstructing the past of human infectious diseases and their pathocenosis

Pr eliminaire d'une etude historique des maladies

Malarial pathocoenosis: beneficial and deleterious interactions between malaria and other human diseases

Microbial genomics of ancient plagues and outbreaks

A single amino acid change in the response regulator PhoP, acquired during Yersinia pestis evolution, affects PhoP target gene transcription and Polymyxin B Susceptibility

Early divergent strains of Yersinia pestis in Eurasia 5,000 years ago

Risk factors for human disease emergence

Dating the emergence of human pathogens

Measurably evolving populations

Measurably evolving pathogens in the genomic era

Bayesian evaluation of temporal signal in measurably evolving populations

The performance of the date-randomization test in phylogenetic analyses of time-structured virus data

Evolution of human influenza A viruses over 50 years: rapid, uniform rate of change in NS gene

Molecular-clock methods for estimating evolutionary rates and timescales

A Bayesian phylogenetic method to estimate unknown sequence ages

Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates

The molecular clock of Mycobacterium tuberculosis

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age

Genome-scale rates of evolutionary change in bacteria

Time-dependent rate phenomenon in viruses

Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis

Efficiently analysing large viral data sets in computational phylogenomics

Sampling bias and model choice in continuous phylogeography: getting lost on a random walk

The peopling of the Pacific from a bacterial perspective

How old are bacterial pathogens?

Tracing ancient human migrations into Sahul using Hepatitis B virus genomes

Plague persistence in Western Europe: a hypothesis

Climate-driven introduction of the Black Death and successive plague reintroductions into

Unifying the epidemiological and evolutionary dynamics of pathogens

Getting to the root of epidemic spread with phylodynamic analysis of genomic data

Evolutionary analysis of the dynamics of viral infectious disease

Jointly inferring the dynamics of population size and sampling intensity from molecular sequences

Inferring infectious disease phylodynamics with notification data

Advancing the ethics of paleogenomics

Yersinia pestis, the cause of plague, is a recently emerged clone of Yersinia pseudotuberculosis

Data, disease and diplomacy: GISAID's innovative contribution to global health

Temporal signal and the phylodynamic threshold of SARS-CoV-2. bioRxiv

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

The proximal origin of SARS-CoV-2

Hepatitis B: The Hunt for a Killer Virus

Tracing hepatitis B virus to the 16th century in a Korean mummy

Dating the origin and dispersal of hepatitis B virus infection in humans and primates

Ancient hepatitis B viruses from the Bronze Age to the Medieval period

Neolithic and medieval virus genomes reveal complex evolution of hepatitis B

Origin and health status of firstgeneration Africans from early colonial Mexico

The potato murrain on the European continent and the revolutions of 1848

What do people die of during famines: the Great Irish Famine in comparative perspective

Reconstructing genome evolution in historic samples of the Irish potato famine pathogen

Persistence of the mitochondrial lineage responsible for the Irish Potato Famine in extant New World Phytophthora infestans

Genomic characterization of a South American Phytophthora hybrid mandates reassessment of the geographic origins of Phytophthora infestans

The rise and fall of the Phytophthora infestans lineage that triggered the Irish potato famine

The Irish potato famine pathogen Phytophthora infestans originated in central Mexico rather than the Andes

A Disputacioun Betwyx 5e Body and Wormes: a translation

Schmutzi: estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA

Beyond ancient microbial DNA: nonnucleotidic biomolecules for paleomicrobiology

Intact smallpox virus particles in an Italian mummy of sixteenth century