Research articlesAncient Yersinia pestis genomes lack the virulence-associated YpfΦ prophage present in modern pandemic strains
Abstract
Yersinia pestis is the causative agent of at least three major plague pandemics (Justinianic, Medieval and Modern). Previous studies on ancient Y. pestis genomes revealed that several genomic alterations had occurred approximately 5000–3000 years ago and contributed to the remarkable virulence of this pathogen. How a subset of strains evolved to cause the Modern pandemic is less well-understood. Here, we examined the virulence-associated prophage (YpfΦ), which had been postulated to be exclusively present in the genomes of strains associated with the Modern pandemic. The analysis of two new Y. pestis genomes from medieval/early modern Denmark confirmed that the phage is absent from the genome of strains dating to this time period. An extended comparative genome analysis of over 300 strains spanning more than 5000 years showed that the prophage is found in the genomes of modern strains only and suggests an integration into the genome during recent Y. pestis evolution. The phage-encoded Zot protein showed structural homology to a virulence factor of Vibrio cholerae. Similar to modern Y. pestis, we observed phages with a common origin to YpfΦ in individual strains of other bacterial species. Our findings present an updated view on the prevalence of YpfΦ, which might contribute to our understanding of the host spectrum, geographical spread and virulence of Y. pestis responsible for the Modern pandemic.
1. Introduction
Yersinia pestis is the pathogenic agent of plague—a zoonotic disease that can be transmitted from rodents to humans via a bite of an infected flea. This route of infection leads to bubonic plague. Yersinia pestis can also spread between humans resulting in the pneumonic plague. When the disease is untreated, the bacterium enters the bloodstream causing sepsis, which is also referred to as septicaemic plague. Although mortality rates of plague vary depending on the clinical form (30%–100%), Y. pestis infection is almost always fatal when no antibiotics are administered in a timely manner. The bacterium is responsible for at least three major pandemics in human history: the Justinianic plague (6th–8th c. AD), the Medieval plague (started in the 14th c. AD, subsequent sporadic outbreaks occurred until the 18th c. AD) and the Modern plague (end of 19th–mid-20th c. AD) [1–3]. Nowadays, Y. pestis still persists in environmental reservoirs around the globe [4] remaining a major threat to public health, as evidenced by the recent outbreaks of plague in Africa [5–7] and Asia [8].
Y. pestis evolved from Yersinia pseudotuberculosis approximately 7000 years ago [9] and gained high virulence by several genomic alterations. These changes include the inactivation of genes and the acquisition of additional virulence genes [10–12]. The impact of those genomic alterations on disease is most drastically demonstrated by their effect on the transmission routes of Y. pestis and the affected organs in humans. Acquiring Yersinia murine toxin (ymt), together with inactivating several other genes (pde2, pde3, ureD, rcsA, flhD), allowed efficient transmission of the pathogen via fleas from rodents to humans and resulted in infection of the lymphatic system [12,13]. The acquisition of the plasminogen activator pla enabled the bacterium to be transmitted by droplets from humans to humans resulting in airway infections and pneumonic plague [14,15]. Taken together, these findings describe the early stages of Y. pestis evolution – dating all the way back to the Neolithic period (approx. 3300–1200 BCE) [9,10,12,13,16]. As much as genetic alterations can aid the rise of pandemics, they might contribute to their fall and the extinction of strains. We recently reported the depletion of pla in medieval post-Black Death strains that have no modern descendants and seem to have gone extinct [17].
Unlike the genetic changes affecting Y. pestis virulence in its early phase of evolution [18], the more recent genome alterations that occurred between the Medieval and Modern pandemics are less well-understood. Among multiple diverse lineages of modern Y. pestis, only the strains of the 1.ORI phylogenetic group are thought to have caused the Modern plague pandemic [4]. Interestingly, the chromosome of 1.ORI strains encodes YpfΦ or CUS-2 filamentous prophage [19–22]. By contrast to other types of bacteriophages, the filamentous phages do not kill the host and pose a relatively limited burden on the bacterium in a form of phage DNA replication and production of its proteins. Assembled virions are then egressed into the outside environment using pore-like channels (secretins) in the bacterium's wall [23]. In turn, the prophage-encoded genes can contribute to the host's virulence [19–21,24]. For instance, the CTX prophage of Vibrio cholerae encodes zonula occludens toxin (Zot) that increases permeability of epithelial barriers in the gut, leading to gastroenteritis [25]. YpfΦ was also shown to enhance Y. pestis virulence, although the exact mechanism and its function remain unknown [21].
In this study, we reconstructed two new Y. pestis genomes from medieval/early modern Denmark, which were subjected to a comparative analysis together with over 300 previously published ancient, medieval and modern Y. pestis genomes. The examined bacterial genomes represent a timespan of over 5000 years. We focused on the virulence-associated filamentous phage (YpfΦ) [19–21], which was absent from the genomes of all ancient and medieval strains.
2. Results
(a) Ancient, medieval and early modern Y. pestis strains lack the virulence-associated prophage YpfΦ
To characterize medieval and early modern Y. pestis genomes, we screened skeletal remains of 42 individuals (see Material) for molecular evidence of Y. pestis infection that were sampled from two medieval cemeteries in Denmark: Sct Trinitatis/Drotten in Viborg and Ødekirkegård in Sejet. Two individuals were excluded due to contamination, leaving 40 individuals in the analysis (17 females and 23 males) (electronic supplementary material, table S1). Yersinia pestis reads were noted in 25% (10/40) of individuals (4 females and 6 males) (electronic supplementary material, table S1). We successfully reconstructed two new pathogen genomes (table 1).
| skeleton no. (Dating [cal AD]) | site name (location no.) | genomic region | no. of aligned reads | mean coverage | coverage ≥1× [%] | coverage ≥2× [%] | coverage ≥3× [%] |
|---|---|---|---|---|---|---|---|
| X52 | Sct. Trinitatis/Drotten, Viborg (VSM F902) | chromosome | 516 891 | 7.0 | 94.8 | 93.1 | 89.6 |
| [1432–1469] | pCD1 | 20 202 | 18.6 | 94.0 | 93.3 | 92.9 | |
| pPCP1 | 107 07 | 72.7 | 80.9 | 80.8 | 80.8 | ||
| pMT1 | 17 892 | 11.6 | 94.0 | 94.0 | 93.3 | ||
| X3003 | Sejet Ødekirkegård (HOM 1046) | chromosome | 1 215 553 | 17.42 | 95.2 | 94.9 | 94.7 |
| [1490–1646] | pCD1 | 53 758 | 52.04 | 95.3 | 94.9 | 94.4 | |
| pPCP1 | 14 843 | 103.82 | 81.3 | 80.9 | 80.9 | ||
| pMT1 | 24 813 | 16.79 | 94.2 | 94.1 | 94.1 |
To explore genomic differences between modern and medieval/early modern Y. pestis, the sequence alignments between the Danish isolates and the modern reference CO92 were inspected for gaps in coverage. An 8734bp-long genomic region was found in the modern CO92 chromosome that was not covered by any reads of the early strains (figure 1a). This region encompasses 13 genes (electronic supplementary material, table S2) and was previously identified as Y. pestis phage YpfΦ or CUS-2 [19–22]. To test if the lack of YpfΦ is specific to the Danish genomes or common among early modern, medieval and ancient Y. pestis, available sequencing reads of 45 medieval/early modern (6th–18th c. AD, electronic supplementary material, table S3) and 9 ancient (Neolithic and Bronze Age, electronic supplementary material, table S4) strains were additionally analysed. None of the ancient or medieval/early modern genomes carried YpfΦ (figure 1b). This finding indicates that the absence of the phage seems to be a generalizable pattern among isolates spanning over five thousand years, including the first two plague pandemics.
Figure 1. YpfΦ of the modern strain CO92 is absent from ancient and medieval Y. pestis strains. (a) YpfΦ of modern Y. pestis strain CO92. Sequencing reads of the medieval strain Sejet Ødekirkegård X3003 (as a representative strain) map to adjacent genes (coloured in grey) but not to the phage (coloured in orange) of the modern reference strain. Individual genes are labelled with their names or locus tags. The positions of the visualized genomic region are indicated. The G + C ratio was calculated using the online GC Content Calculator [26]. (b) Ancient and medieval strains lack YpfΦ (empty circles). Individual strains are indicated along a timeline based on their estimated dating. The medieval Danish strains generated in this study (Sejet Ødekirkegård X3003 and Viborg Sct. Trinitatis/Drotten X52) are marked in bold. BP, before present; c. century; AD, anno domini.
(b) A subset of modern Y. pestis strains stably integrated YpfΦ in their genome
To examine the prevalence of YpfΦ among modern Y. pestis strains, 255 published genomes were screened for the phage (electronic supplementary material, table S3). The screening revealed YpfΦ in the chromosomes of a subset of strains (n = 45; electronic supplementary material, table S5), which clustered together in one subbranch of branch 1 in the phylogenetic tree (figure 2; electronic supplementary material, figure S1). All isolates with YpfΦ belonged to the phylogenetic group 1.ORI that is thought to be responsible for the Modern plague pandemic. All remaining branch 1 strains (1.ANT and 1.IN groups) did not carry the phage in their genomes and hereby will be referred to as the 1.ANT + IN strains. Although two previous studies suggested an extrachromosomal presence of YpfΦ in several modern strains from different phylogenetic groups based on a PCR-analysis [21,27], our genome-based approach did not detect such extrachromosomal elements. This discrepancy could be explained by differences in the applied methods and the strains examined. Unfortunately, genome sequences are not available of the non-1.ORI strains that had been analysed by PCR.
Figure 2. Maximum-likelihoodphylogeny of 255 modern and 47 medieval/early modern Y. pestis strains. The two Danish strains from Sejet Ødekirkegård and Viborg (Sct. Trinitatis/Drotten) are marked in bold. Medieval Black Death strains are shown in violet and the post-Black Death strains in purple. All strains carrying YpfΦ in their genome belong to the modern phylogenetic group 1.ORI (orange). Remaining modern branch 1 strains form groups 1.ANT and 1.IN (1.ANT + IN). The tree is based on the SNP alignment of 10 315 positions with a bootstrap of 500 replicates and includes Y. pseudotuberculosis as an outgroup. The full tree with uncollapsed branches and bootstrap values can be found in the supplementary material (electronic supplementary material, figure S1).
Next, we analysed the intraspecific diversity of YpfΦ among the different strains. Our results revealed the existence of multiple YpfΦ variants in the existing assembled sequences (electronic supplementary material, figure S2, table S6). Some strains encoded two copies of the phage in a row (e.g. Java9) whereas other strains encoded one copy of the entire phage followed by an incomplete second copy (e.g. IP275). Moreover, pseudogenes in YpfΦ of several strains were observed. In strains with two copies of YpfΦ, one copy often contained a pseudogene whereas the other copy was marked as functional. This observation might indicate selection for one fully functional YpfΦ.
In the phylogeny, the two new Danish genomes (Sejet Ødekirkegård X3003 and Sct. Trinitatis/Drotten X52, Viborg (VSM F902)) without YpfΦ clustered within the known diversity of medieval Y. pestis strains (figure 2). Interestingly, although the C14 dating (electronic supplementary material, Excel sheet S1) places both strains in the post-Black Death period of the pandemic (1354 AD to 18th c. AD), Sejet X3003 clustered together with Y. pestis isolates responsible for the Black Death (1346–1353 AD). In addition, unlike Sct. Drotten X52, Sejet X3003 did not exhibit depletion of the pla region in the pPCP1 plasmid that was shown to be characteristic for the strains of the post-Black Death period [17]. These findings suggest that the Black Death Y. pestis lineage persisted in Denmark for at least over 150 years.
(c) Common ancestry of YpfΦ in Y. pestis and other Enterobacteria
To determine a possible origin of YpfΦ in Y. pestis, the NCBI database was screened for bacterial strains carrying analogous sequences. We found a highly similar genomic region in the Escherichia coli strain MOD1-EC6770 (figure 3a). Unlike in previous observations [20], sequencing reads and contigs of the identified E. coli isolate aligned to full-length YpfΦ of Y. pestis CO92 with high nucleotide identities (figure 3a, S3). This finding suggests a common ancestry of both phages in E. coli and Y. pestis.
Figure 3. YpfΦ is encoded by other Enterobacteria. (a) Alignment of YpfΦ from E.coli MOD1-EC6770 (three contigs) and Y. pestis CO92 (NC_003143.1). SNPs are marked in colours along the contigs. A graphical depiction of the phage is shown above the alignment in orange. The numbers on the gene symbols have the prefix ‘YPO_RS12’ that was removed for clarity. In the NCBI the annotation for zot and gspd2 are YPO_RS12400 and YPO_RS12405, respectively. (b) S. enterica PNUSAS053175 carries a putative phage with a similar genetic architecture to YpfΦ. A graphical depiction of the respective genes and annotation is shown. The numbers on the gene symbols have the prefix ‘D3500_09’ that was removed for clarity. Genes in the phage of S. enterica that encode proteins exhibiting at least 30% amino acid identity in pairwise alignments with the respective proteins of Y. pestis are marked in orange. Amino acid sequence identities (AAI [%]) are indicated. HP stands for hypothetical protein.
Slightly truncated forms of YpfΦ are also found in other Enterobacteria such as Enterobacter ludwigii and hormarcheri, Morganella morganii, Klebsiella aerogenes, Salmonella enterica, Citrobacter amalonaticus, koseri and portucalensis as well as E. coli and Cedecea davisiae. Those truncated forms encode 11 genes and lack equivalents to the two open reading frames at the 5′ end of Y. pestis YpfΦ, namely YPO_RS12355 and YPO_RS12360. In Citrobacter freundi, Enterobacter cloacae and Shigella sonnei strains, we also found multiple copies of the phage. An alignment of the YpfΦs of various species reveals over 98.7% nucleotide identity over the 11 conserved genes (electronic supplementary material, table S7) and suggests high conservation across a wide range of enterobacterial isolates of mostly human origin.
The YpfΦ of Y. pestis might belong to a large family of phages with a modular genomic architecture. The YpfΦ-like phage in S. enterica PNUSAS053175 encodes common and distinct features to YpfΦ of Y. pestis (figure 3b). Whereas the first five genes (supposedly involved in phage regulation and replication) share high homology between the two phages, the proteins encoded in the remaining eight predicted open reading frames (mediating phage morphogenesis and secretion) differ strongly and share less than 30% amino acid identity. Interestingly, the S. enterica strain does not encode a homologue to zot and future work will need to resolve whether the proteins encoded by the YpfΦ-like phages also act as toxins. In the following, we examine the genes encoded by YpfΦ to better understand its potential role in Y. pestis.
(d) Ypfφ encodes proteins with structural homology to zonula occludens toxin and T2SS secretin
Interestingly, one of the YpfΦ genes (zot) encodes a protein homologous to zonula occludens toxin (Zot) (electronic supplementary material, figure S4), which had recently been detected in the Y. pestis CO92 strain in an independent study [28]. 3D structure prediction of the Y. pestis and V. cholerae Zot proteins shows a partial structure similarity (figure 4a). Both proteins have the C-terminal (CT) and N-terminal (NT) domains that for V. cholerae were shown to reside in the bacterial periplasm and cytoplasm, respectively [29]. In V. cholerae, the CT domain contains the biologically active region (FCIGRL sequence at the 288–293 amino acid position), which modulates tight junctions of the epithelial cells [30]. Like other human pathogens, such as Campylobacter concisus and Vibrio parahaemolyticus, Y. pestis's Zot lacks the FCIGRL (electronic supplementary material, figure S4). The presence of the FCIGRL sequence is, however, unnecessary for the Zot-mediated disruption of tight junctions for those bacterial species [31]. Thus, the activity of Y. pestis's Zot might also be mediated via a different active site.
Figure 4. YpfΦ-encoded molecules exhibit homology to Zot—a virulence factor of V. cholerae (a) and secretin of the T2SS (b). Predicted structures are coloured based on the pLDDT confidence measure. (a) Predicted three-dimensional protein structure of Y. pestis's Zot reveals a partial structural similarity to Zot of V. cholerae. Both proteins have a transmembrane region (TR) as well as C-terminal (ZOT-CT) and N-terminal (ZOT-NT) domains. The Y. pestis’ Zot lacks the active region (FCIGRL) present in V. cholerae. (b) Predicted three-dimensional protein structures reveal differences in the protein core as well as in the periplasmic domain that is relevant for substrate specificity. The main domains (Nx and S) are labelled. The nomenclature is based on the resolved structure of the GspD of Klebsiella pneumoniae [31].
In addition to Zot, YpfΦ encodes a gene for a protein that is annotated as a type II secretion system (T2SS) secretin and that we call GspD2 (figure 4b). The chromosome of Y. pestis already encodes a protein called GspD in a gene outside of YpfΦ (electronic supplementary material, figure S5A). To assess the difference between the two proteins, their size was compared and an alignment of the amino acid sequences of GspD and GspD2 was performed (electronic supplementary material, figure S5B). While GspD is built of 640 amino acids (aa), GspD2 is smaller (414 aa). Only 26% amino acid identity was detected between the two proteins in the alignment.
To further test if YpfΦ-encoded GspD2 and GspD share similar three-dimensional structures despite sequence dissimilarities, the structures of the two Y. pestis proteins were modelled (figure 4b). The three-dimensional models of the two proteins look surprisingly similar given the low amino acid sequence identity. Both proteins form a core with an inner and outer barrel as well as the N-terminal extensions of alpha helices. In comparison to GspD, GspD2 lacks two periplasmic domains (N1 and N2) and exhibits an outer barrel of reduced complexity. However, the pLDDT score, that is a measure of prediction accuracy, shows a relatively low confidence for the model of the GspD2 outer barrel. Apart from the GspD2 outer barrel, however, the pLDDT indicates a good prediction for both proteins. Both GspD and GspD2 consist of domains characteristic for an outer membrane secretin channel of T2SS, including the periplasmic domain and the domains forming the pore structure (the barrels and the S domain) [32]. In sum, the high similarity in overall structure, reflected in the presence of domains typical for T2SS secretin, suggests that GspD and GspD2 may both function as such, although with different substrate specificities as determined by the N-terminal domains in the periplasm. Based on a possible role of Zot in phage assembly and release [33] (in reference 33 Zot is annotated as YPO2279), GspD2 might form a complex with Zot to facilitate the secretion of YpfΦ.
(e) Yersinia pestis strains with and without YpfΦ differ in their host spectrum
To better understand the potential impact of the prophage YpfΦ on bacteria beyond its known association with virulence [21], we determined features that are associated with modern YpfΦ-positive strains (1.ORI) and are absent from modern YpfΦ-negative branch 1 strains (1.ANT + IN). Of particular interest is the source of the analysed isolates, as Y. pestis is found in multiple animal hosts besides humans. These animal hosts serve as natural reservoirs of this pathogenic bacterium and play a key role during transmission. When comparing the host spectrum between strains with and without the phage, we found YpfΦ-positive strains (1.ORI) predominantly among animals associated with human habitats (no. of hosts = 9), like rats, mice, cats and dogs, whereas YpfΦ-negative isolates are found among more diverse animals (no. of hosts = 15) (figure 5a; electronic supplementary material, table S8). Furthermore, the proportion of the human host was higher for the YpfΦ-positive strains (20/44, 45.45%) relative to the YpfΦ-negative strains (7/31, 22.6%) (figure 5a). Although a strong trend was noted, it was not statistically significant (p = 0.0527). However, the proportion of human isolates among YpfΦ-positive bacteria that is seen here is likely an underestimate, considering that the analysed isolates were chosen to represent maximal host diversity. The influence of sampling bias cannot be excluded, as in many cases the isolate source is unknown or the host's natural habitat is undefined, i.e. the animal can be found in both human-associated and wild environment (electronic supplementary material, table S8). Nevertheless, we observe an altered host range among isolates with and without the phage in the analysed sample (figure 5b). Another feature of the YpfΦ-positive strains is that they all belong to the phylogenetic group 1.ORI which has been linked with the Modern plague pandemic and subsequent disease outbreaks (figure 5b).
Figure 5. Host spectrum of YpfΦ-negative and YpfΦ-positive strains of branch 1. (a) Proportion of different hosts of 1.ANT + IN (YpfΦ-negative, n = 31) and 1.ORI (YpfΦ-positive, n = 45) strains analysed in this study. Isolation sources are indicated with respective colours. More details about the sources can be found in electronic supplementary material, table S8. (b) Differences between the YpfΦ-positive and YpfΦ-negative Y. pestis strains with regards to the host spectrum and pandemic potential. The number of different hosts for 1.ORI was 9, while for 1.ANT + IN it was 15. Illustrations were created with BioRender.com.
3. Discussion
In this study, we confirm that YpfΦ is associated with Y. pestis responsible for the Modern plague pandemic (1.ORI) (figure 2). Furthermore, the presence of YpfΦ in modern Y. pestis strains of the 1.ORI group correlates with their seemingly altered host spectrum and pandemic potential relative to other branch 1 modern Y. pestis strains (YpfΦ-negative) (figure 5). YpfΦ was also absent in the genomes of strains which were responsible for the Medieval and Justinianic plague pandemics as well as those in the Neolithic or Bronze Age (figure 1b). The lack of coverage of YpfΦ genes in medieval strains can also be seen in the data published by Seguin-Orlando and colleagues (fig. 4A in Seguin-Orlando et al., [34]). As the phage is only present in one modern phylogenetic sub-branch, it is possible to estimate the approximate time of its incorporation based on the molecular dating of the splits in the phylogeny. Molecular dating depends on various factors, such as the included strains and the methods used for the analysis, and thus there seems to be no consensus in the literature concerning the exact date of the 1.ORI branch split [35]. Although the 1.ORI branch could have emerged anytime between approximately 550 and 150 years ago [9,35–37], it is likely that both the split and the genomic acquisition of YpfΦ occurred before the Modern plague pandemic.
As previous functional studies in mice showed that YpfΦ enables better colonization of the host and confers increased fitness during infection of mammals [21], the chromosomal acquisition of the phage likely influenced the pathogenicity of the bacterium towards humans. Changes in pathogenicity could possibly affect the clinical manifestation of plague, suggesting possible differences between the ancient (Justinianic and Medieval) and Modern pandemics. Differences in pathogenicity do not necessarily have to be reflected in the estimated mortality rates of the pandemics, as various factors, such as proximity to zoonotic reservoirs, climate, settlement type/crowding, health status and behaviour-related exposure, can influence the outbreak outcomes (e.g. [38–42]). This finding highlights the introduction of potential bias when interpreting various historical sources on plague epidemiology in the past based on the inference from knowledge about modern plague.
Y. pestis has evolved its remarkable virulence and transmission abilities via horizontal acquisition of large pieces of foreign DNA, such as plasmids (pPCP1, pPMT1) and pathogenicity islands. These alien genetic elements, while non-essential for survival, were probably crucial for the emergence of pathogenic Yersinia in general [43,44]. For instance, the virulence of pathogenic Yersinia (Y. pestis, Y. pseudotuberculosis and Y. enterocolitica) strictly depends on the presence of the pCD1 plasmid [44]. Furthermore, within these three species, certain subgroups have acquired a genomic region which allows a systemic dissemination of the bacterium and confers the high-virulence phenotype. Therefore, the region is called a high pathogenicity island (HPI) [45]. Both HPI and, to a lesser extent, YpfΦ are linked to an increased fitness during the infection process and their genomic acquisition represents rapid modification in bacterial pathogenicity [21,44–46]. By contrast to the relatively well described role of the HPI and the Y. pestis plasmids in the bacterium's virulence and transmission, the mechanisms responsible for the increased dissemination of the pathogen that is conferred by YpfΦ have so far remained unknown. Here, we show that YpfΦ encodes a protein homologous to zonula occludens toxin (Zot) that is a virulence factor in other pathogenic bacteria, such as Neisseria meningitidis, Acinetobacter baumanii, Salmonella enterica, Burkholderia cenocepacia and Vibrio cholerae, in which it had first been identified [25,28]. In these pathogens, Zot increases permeability of the gut epithelia [25], the mucosa [47] and endothelial cells in the brain [48]. Based on our findings, we propose that proteins encoded in the phage (i.e. Zot and GspD2) represent good candidates for further experimental studies on the molecular mechanism of phage-mediated virulence.
Similar to HPI [44,46], YpfΦ is found across various bacterial genera (electronic supplementary material, table S6) [20,28]. Presence of complete YpfΦ, as seen in Y. pestis, in the genome of E. coli MOD1-EC6770 suggests this species as a possible source of the YpfΦ in Y. pestis. Such horizontal spread of large genetic elements, like prophages, plasmids and pathogenicity islands, often cause ‘quantum leaps’ in evolution [49], affecting the bacterium's metabolism, transmission and virulence. For Y. pestis, the acquisition of pPCP1 and pMT1 plasmids as well as HPI allowed for a flea- and airborne transmission and rapid systemic dissemination of the pathogen—a deadly combination that kills up to 100% of infected individuals [50]. In recent years, several Y. pestis strains acquired additional plasmids (pIP1202, pIP203 and pIP2180H) that confer the ability to resist antibiotic treatment [51]. YpfΦ likely represents another ‘leap’ in the evolution of the pathogen. A better knowledge on the recent evolution of Y. pestis is key to the understanding of the bacterium's pathogenicity. Acquisition of foreign genetic material, which can include virulence factors, might lead to unusual clinical forms of plague that can be challenging to diagnose. For instance, during a plague-related outbreak of gastroenteritis in Afghanistan, 20.5% of the infected individuals (17/83) died before the appropriate treatment was applied [52]. If the YpfΦ-encoded Zot indeed disrupts the tight junctions of epithelial cells, like in several other bacteria [25,31], the gastroenteritis can be a result of an infection with the phage-positive Y. pestis.
4. Material
The skeletal material used for Y. pestis screening comprised 74 bones and teeth belonging to 42 individuals from two medieval/early modern parish cemeteries (1000–1575 AD) [53–56] in Denmark (electronic supplementary material, table S1). The skeletal material analysed in this study is stored at the ADBOU skeletal collection (Department of Forensic Medicine, University of Southern Denmark).
The YpfΦ distribution across the Y. pestis strains was analysed among previously published 255 modern, 45 medieval/early modern (6th–18th century AD; electronic supplementary material, table S3) and 9 ancient (Neolithic and Bronze Age, electronic supplementary material, table S4) strains. Two new genomes (Sejet Ødekirkegård X3003 and Sct. Trinitatis/Drotten X52, Viborg (VSM F902)) were also analysed. Sequencing data was downloaded from online repositories such as the European Nucleotide Archive (ENA) and the National Center for Biotechnology Information (NCBI).
5. Methods
(a) Processing of the metagenomic medieval/early modern samples
All DNA samples were extracted and processed in a dedicated ancient DNA facility at the University of Kiel following the guidelines on contamination control in ancient DNA [57–59], according to a previously published protocol for the non-UGD treated samples [60]. Shotgun sequencing was performed on the Illumina HiSeq 6000 (2 × 100) platform of the Institute of Clinical Molecular Biology in Kiel.
Adapter sequences were removed and paired-end reads were merged with ClipAndMerge v1.7.7. [61]. Shotgun sequence data was mapped to the human genome (build hg19) using BWA v0.7.12 [62] with a reduced mapping stringency parameter ‘-n 0.01’ to account for mismatches in aDNA. Duplicated reads were removed with DeDup v0.12.2 [61].
To confirm the ancient origin of the sequences, terminal damage of the reads (C to T substitutions) was assessed with DamageProfiler [63]. After the validation, the first two positions from the 5′ and 3′-ends of the reads were trimmed. Furthermore, X-chromosome and mitochondrial DNA contamination were assessed with ANGSD and Schmutzi, respectively [64,65].
(b) Identification of Y. pestis-positive samples in medieval and early modern individuals
Initial screening of all samples for the presence of Y. pestis DNA was performed with Megan Alignment Tool 0.3.0 (MALT) [66] (SemiGlobal alignment mode, identity threshold = 90%), using a custom database containing bacterial genomes available at the NCBI platform (24.01.2019). Output alignments were inspected visually in MEGAN 6 [67]. The samples that contained reads aligning to Y. pestis were further evaluated.
The Y. pestis-positive status was based on detection of reads unique for the pathogen in the sample. Y. pestis-specific reads were obtained in competitive mapping against Y. pestis (NC_003143.1, NC_003131.1, NC_003134.1, NC_003132.1) and Yersinia pseudotuberculosis (NC_006155.1) reference genomes, using Burrows-Wheeler Aligner (BWA) v0.7.12 (n = 0.01, l = 300) [62]. Output BAM files were then filtered for quality 30 with SAMtools [68]. Subsequently, the number of Y. pestis-specific reads was noted with samtools idxstats.
(c) Alignment and detection of YpfΦ
Sequencing data from all Y. pestis strains (electronic supplementary material, table S3 and S4) were mapped against the CO92 reference genome (NC_003143.1, NC003131.1, NC_003132.1, NC_003134.1) with BWA v0.7.12 [62]. For ancient data, a reduced stringency parameter was used (-n 0.01) to account for mismatches in ancient DNA and two positions from the 5′ and 3′-ends of the reads were trimmed. Duplicated reads were removed with DeDup v.0.12.2 [61]. Regions with zero coverage were identified and the samples with no reads mapping to the NC_003143.1 : 2554178–2562912 region were classified as strains without YpfΦ. Moreover, sequencing data of an example strain CMCC10012 (not carrying the phage) was mapped in a competitive mapping against the CO92 reference genome and the CO92 chromosome without the 2554178–2562912 region. This way, unique gap-bridging reads were identified for the CMCC10012, confirming the lack of YpfΦ.
(d) Analysis of YpfΦ
Reads mapping to the CO92 YpfΦ region were extracted from a randomly chosen example strain (EV76). The phage reads were then mapped against previously published genome assemblies and contigs of the YpfΦ-positive Y. pestis strains (electronic supplementary material, table S5) with BWA v0.7.12 (n = 0.01, l = 300) [62] to locate the position of YpfΦ within the chromosome of each strain. Knowing the phage locus for the strains, genome annotation was inspected in the graphic view panel of the sequences in NCBI to identify the YpfΦ variants. Association with a particular phylogenetic group was made based on available data [37] or the position in the phylogeny. To identify similar sequences in other bacterial species, the extracted phage reads were blasted with the BLASTn online tool (default parameters) against the nucleotide collection (nt).
Differences in host diversity between modern 1.ORI and 1.ANT + IN strains were assessed with Fisher's exact test using IBM SPSS Statistics (v. 26).
(e) Comparative analysis of Zot and secretins
Three-dimensional (3D) structures of CO92 Y. pestis secretins (GspD and GspD2) as well as Y. pestis and V. cholerae Zot molecules were predicted using Alpha Fold with default parameters [69]. MUSCLE online tool [70] (ClustalW) was used for multiple sequence alignment of amino acid sequences.
(f) Phylogenetic analysis
Phylogenetic analysis was performed with RAxML [71] using the GTRGAMMA model with 500 bootstrap replicates. MultiVCFAnalyzer [72] was used to generate a SNP-based multiple alignment of 255 previously published modern and 47 medieval/early modern Y. pestis strains (including the two Danish genomes reconstructed in this study: Sejet Ødekirkegård X3003 and Sct. Trinitatis/Drotten X52, Viborg (VSM F902) with Y. pseudotuberculosis (NZ_CP008943.1) as an outgroup (electronic supplementary material, table S3). The input VCF files were generated with the UnifiedGenotyper module from the Genome Analysis Toolkit (GATK) v3.6 [73]. A SNP was called if the position was covered by at least three reads, the genotype quality was at least 30 and the fraction of mapped reads containing the SNP was at least 90%.
Ethics
Samples extracted from human skeletal remains were analysed in this study. Due to the archaeological nature of the material, no ethical approval was required. Approval for the study was given by the respective museums and the curators of the skeletal collections who are co-authors of this study. The letters of approval can be found in the electronic supplementary material, Information.
Data accessibility
Sequences used in the reconstruction and subsequent analysis of Y. pestis genomes from Sejet X3003 and Viborg X52 are available through the European Nucleotide Archive (http://www.ebi.ac.uk/ena/browser/view/) under Accession Number PRJEB60595.
The data are provided in electronic supplementary material [74].
Authors' contributions
J.H.B.: conceptualization, formal analysis, investigation, methodology, visualization, writing—original draft; J.S.: methodology; B.K.-K.: resources, writing—review and editing; D.D.P.: data curation, writing—review and editing; J.B.: data curation; L.A.L.: data curation, writing—review and editing; L.S.: data curation, writing—review and editing; A.N.: supervision, writing—review and editing; D.U.: conceptualization, methodology, supervision, visualization, writing—original draft.
All authors gave final approval for publication and agreed to be held accountable for the work performed therein.
Conflict of interest declaration
We declare we have no competing interests.
Funding
Open access funding provided by the Max Planck Society.
J.H.B. was funded by the International Max Planck Research School for Evolutionary Biology and the Excellence Cluster Precision Medicine in Chronic Inflammation (PMI) (excellence strategy – EXC 2167–390884018). Work in the Unterweger Lab is supported by the German Federal Ministry for Education and Research (grant 01KI2020).
Acknowledgements
We thank the Kiel Evolution Center for facilitating this collaborative project and the Museums of Horsens and Viborg for permission to sample the skeletal material. We are indebted to Prof. Dr Holger Sondermann and Dr Steffi Jimmy from the Centre for Structural Systems Biology (CSSB, Hamburg, Germany) as well as Ekaterina Ovchinnikova from the Institute for Experimental Medicine at Kiel University for their help during the revision process.
References
- 1.
Haensch S 2010 Distinct clones of Yersinia pestis caused the black death. PLoS Pathog. 6, e1001134. (doi:10.1371/journal.ppat.1001134) Crossref, PubMed, Web of Science, Google Scholar - 2.
Harbeck M 2013 Yersinia pestis DNA from skeletal remains from the 6th century AD reveals insights into Justinianic plague. PLoS Pathog. 9, e1003349. (doi:10.1371/journal.ppat.1003349) Crossref, PubMed, Web of Science, Google Scholar - 3.
Pryor E . 1975 The great plague of Hong Kong. J. Hong Kong Branch R. Asiat. Soc. 15, 61-70. PubMed, Google Scholar - 4.
Vogler AJ, Chan F, Nottingham R, Andersen G, Drees K, Beckstrom-Sternberg SM, Wagner DM, Chanteau S, Keim P . 2013 A decade of plague in Mahajanga, Madagascar: insights into the global maritime spread of pandemic plague. mBio 4, e00623-12. (doi:10.1128/mBio.00623-12) Crossref, PubMed, Web of Science, Google Scholar - 5.
Abedi AA 2018 Ecologic features of plague outbreak areas, democratic republic of the Congo, 2004–2014. Emerg. Infect. Dis. 24, 210-220. (doi:10.3201/eid2402.160122) Crossref, PubMed, Web of Science, Google Scholar - 6.
Randremanana R . 2019 Epidemiological characteristics of an urban plague epidemic in Madagascar, August–November, 2017: an outbreak report. Lancet Infect. Dis. 19, 537-545. (doi:10.1016/S1473-3099(18)30730-8) Crossref, PubMed, Web of Science, Google Scholar - 7.
Respicio-Kingry LB 2016 Two distinct Yersinia pestis populations causing plague among humans in the West Nile region of Uganda. PLoS Negl. Trop. Dis. 10, e0004360. (doi:10.1371/journal.pntd.0004360) Crossref, PubMed, Web of Science, Google Scholar - 8.
Shi L 2018 Reemergence of human plague in Yunnan, China in 2016. PLoS ONE 13, e0198067. (doi:10.1371/journal.pone.0198067) Crossref, PubMed, Web of Science, Google Scholar - 9.
Susat J 2021 A 5,000-year-old hunter-gatherer already plagued by Yersinia pestis. Cell Rep. 35, 109278. (doi:10.1016/j.celrep.2021.109278) Crossref, PubMed, Web of Science, Google Scholar - 10.
Rasmussen S 2015 Early divergent strains of Yersinia pestis in Eurasia 5,000 years ago. Cell 163, 571-582. (doi:10.1016/j.cell.2015.10.009) Crossref, PubMed, Web of Science, Google Scholar - 11.
Spyrou MA 2019 Phylogeography of the second plague pandemic revealed through analysis of historical Yersinia pestis genomes. Nat. Commun. 10, 4470. (doi:10.1038/s41467-019-12154-0) Crossref, PubMed, Web of Science, Google Scholar - 12.
Sun Y-C, Jarrett CO, Bosio CF, Hinnebusch BJ . 2014 Retracing the evolutionary path that led to flea-borne transmission of Yersinia pestis. Cell Host Microbe 15, 578-586. (doi:10.1016/j.chom.2014.04.003) Crossref, PubMed, Web of Science, Google Scholar - 13.
Chouikha I, Hinnebusch BJ . 2014 Silencing urease: a key evolutionary step that facilitated the adaptation of Yersinia pestis to the flea-borne transmission route. Proc. Natl. Acad. Sci. USA 111, 18 709-18 714. (doi:10.1073/pnas.1413209111) Crossref, Web of Science, Google Scholar - 14.
Lathem WW, Price PA, Miller VL, Goldman WE . 2007 A plasminogen-activating protease specifically controls the development of primary pneumonic plague. Science 315, 509-513. (doi:10.1126/science.1137195) Crossref, PubMed, Web of Science, Google Scholar - 15.
Perry RD, Fetherston JD . 1997 Yersinia pestis: etiologic agent of plague. Clin. Microbiol. Rev. 10, 35-66. (doi:10.1128/CMR.10.1.35) Crossref, PubMed, Web of Science, Google Scholar - 16.
Zimbler DL, Schroeder JA, Eddy JL, Lathem WW . 2015 Early emergence of Yersinia pestis as a severe respiratory pathogen. Nat. Commun. 6, 7487. (doi:10.1038/ncomms8487) Crossref, PubMed, Web of Science, Google Scholar - 17.
Susat J, Bonczarowska JH, Pētersone-Gordina E, Immel A, Nebel A, Gerhards G, Krause-Kyora B . 2020 Yersinia pestis strains from Latvia show depletion of the pla virulence gene at the end of the second plague pandemic. Sci. Rep. 10, 14628. (doi:10.1038/s41598-020-71530-9) Crossref, PubMed, Web of Science, Google Scholar - 18.
Demeure CE, Dussurget O, Mas Fiol G, le Guern A-S, Savin C, Pizarro-Cerdá J . 2019 Yersinia pestis and plague: an updated view on evolution, virulence determinants, immune subversion, vaccination, and diagnostics. Genes Immun. 20, 357-370. (doi:10.1038/s41435-019-0065-0) Crossref, PubMed, Web of Science, Google Scholar - 19.
Gonzalez MD, Lichtensteiger CA, Vimr ER . 2001 Adaptation of signature-tagged mutagenesis to Escherichia coli K1 and the infant-rat model of invasive disease. FEMS Microbiol. Lett. 198, 125-128. (doi:10.1111/j.1574-6968.2001.tb10630.x) Crossref, PubMed, Web of Science, Google Scholar - 20.
Gonzalez MD, Lichtensteiger CA, Caughlan R, Vimr ER . 2002 Conserved filamentous prophage in Escherichia coli O18:K1:H7 and Yersinia pestis Biovar orientalis. J. Bacteriol. 184, 6050-6055. (doi:10.1128/JB.184.21.6050-6055.2002) Crossref, PubMed, Web of Science, Google Scholar - 21.
Derbise A, Chenal-Francisque V, Pouillot F, Fayolle C, Prévost M, Médigue C, Hinnebusch BJ, Carniel E . 2007 A horizontally acquired filamentous phage contributes to the pathogenicity of the plague bacillus. Mol. Microbiol. 63, 1145-1157. (doi:10.1111/j.1365-2958.2006.05570.x) Crossref, PubMed, Web of Science, Google Scholar - 22.
Hay ID, Lithgow T . 2019 Filamentous phages: masters of a microbial sharing economy. EMBO Rep. 20, e47427. (doi:10.15252/embr.201847427) Crossref, PubMed, Web of Science, Google Scholar - 23.
Conners R 2021 CryoEM structure of the outer membrane secretin channel pIV from the f1 filamentous bacteriophage. Nat. Commun. 12, 6316. (doi:10.1038/s41467-021-26610-3) Crossref, PubMed, Web of Science, Google Scholar - 24.
Rice SA 2009 The biofilm life cycle and virulence of Pseudomonas aeruginosa are dependent on a filamentous prophage. ISME J. 3, 271-282. (doi:10.1038/ismej.2008.109) Crossref, PubMed, Web of Science, Google Scholar - 25.
Fasano A, Baudry B, Pumplin DW, Wasserman SS, Tall BD, Ketley JM, Kaper JB . 1991 Vibrio cholerae produces a second enterotoxin, which affects intestinal tight junctions. Proc. Natl. Acad. Sci. USA. 88, 5242-5246. (doi:10.1073/pnas.88.12.5242) Crossref, PubMed, Web of Science, Google Scholar - 26.
Guerra ÁP, Calvo EP, Wasserman M, Chaparro-Olaya J . 2016 Production of recombinant proteins from Plasmodium falciparum in Escherichia coli. Biomédica. 36, 97-108. (doi:10.7705/biomedica.v36i3.3011) PubMed, Web of Science, Google Scholar - 27.
Li Y 2008 Different region analysis for genotyping Yersinia pestis isolates from China. PLoS ONE 3, e2166. (doi:10.1371/journal.pone.0002166) Crossref, PubMed, Web of Science, Google Scholar - 28.
Yap KP, Gan HM, Teh CSJ, Chai LC, Thong KL . 2014 Comparative genomics of closely related serovar Typhi strains reveals genome dynamics and the acquisition of novel pathogenic elements. BMC Genomics 15, 1007. (doi:10.1186/1471-2164-15-1007) Crossref, PubMed, Web of Science, Google Scholar - 29.
Schmidt E, Kelly SM, van der Walle CF . 2007 Tight junction modulation and biochemical characterisation of the zonula occludens toxin C- and N-termini. FEBS Lett. 581, 2974-2980. (doi:10.1016/j.febslet.2007.05.051) Crossref, PubMed, Web of Science, Google Scholar - 30.
Goldblum SE 2011 The active Zot domain (aa 288–293) increases ZO-1 and myosin 1C serine/threonine phosphorylation, alters interaction between ZO-1 and its binding partners, and induces tight junction disassembly through proteinase activated receptor 2 activation. FASEB J. 25, 144-158. (doi:10.1096/fj.10-158972) Crossref, PubMed, Web of Science, Google Scholar - 31.
Pérez-Reytor D 2020 Analysis of the Zonula occludens toxin found in the genome of the Chilean non-toxigenic Vibrio parahaemolyticus strain PMC53.7. Front. Cell Infect. Microbiol. 10, 482. (doi:10.3389/fcimb.2020.00482) Crossref, PubMed, Google Scholar - 32.
Chernyatina AA, Low HH . 2019 Core architecture of a bacterial type II secretion system. Nat. Commun. 10, 5437. (doi:10.1038/s41467-019-13301-3) Crossref, PubMed, Web of Science, Google Scholar - 33.
Chouikha I, Charrier L, Filali S, Derbise A, Carniel E . 2010 Insights into the infective properties of YpfΦ, the Yersinia pestis filamentous phage. Virology 407, 43-52. (doi:10.1016/j.virol.2010.07.048) Crossref, PubMed, Web of Science, Google Scholar - 34.
Seguin-Orlando A 2021 No particular genomic features underpin the dramatic economic consequences of 17th century plague epidemics in Italy. iScience 24, 102383. (doi:10.1016/j.isci.2021.102383) Crossref, PubMed, Web of Science, Google Scholar - 35.
Eaton K, Featherstone L, Duchene S, Carmichael AG, Varlık N, Golding GB, Holmes EC, Poinar HN . 2023 Plagued by a cryptic clock: insight and issues from the global phylogeny of Yersinia pestis. Commun. Biol. 6, 23. (doi:10.1038/s42003-022-04394-6) Crossref, PubMed, Web of Science, Google Scholar - 36.
Pisarenko SV, Evchenko A, Kovalev DA, Evchenko YМ, Bobrysheva OV, Shapakov NA, Volynkina AS, Kulichenko AN . 2021 Yersinia pestis strains isolated in natural plague foci of Caucasus and Transcaucasia in the context of the global evolution of species. Genomics 113, 1952-1961. (doi:10.1016/j.ygeno.2021.04.021) Crossref, PubMed, Web of Science, Google Scholar - 37.
Cui Y 2013 Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis. Proc. Natl. Acad. Sci. USA. 110, 577-582. (doi:10.1073/pnas.1205750110) Crossref, PubMed, Web of Science, Google Scholar - 38.
Slavin P . 2021 Out of the West: formation of a permanent plague reservoir in South-Central Germany (1349–1356) and its implications. Past Present 252, 3-51. (doi:10.1093/pastj/gtaa028) Crossref, Google Scholar - 39.
Stenseth NC 2006 Plague dynamics are driven by climate variation. Proc. Natl. Acad. Sci. USA. 103, 13 110-13 115. (doi:10.1073/pnas.0602447103) Crossref, Web of Science, Google Scholar - 40.
Izdebski A 2022 Palaeoecological data indicates land-use changes across Europe linked to spatial heterogeneity in mortality during the Black Death pandemic. Nat. Ecol. Evol. 6, 297-306. (doi:10.1038/s41559-021-01652-4) Crossref, PubMed, Web of Science, Google Scholar - 41.
DeWitte SN . 2010 Sex differentials in frailty in medieval England. Am. J. Phys. Anthropol. 143, 285-297. (doi:10.1002/ajpa.21316) Crossref, PubMed, Web of Science, Google Scholar - 42.
Cleri DJ, Vernaleo JR, Lombardi LJ, Rabbat MS, Mathew A, Marton R, Reyelt MC . 1997 Plague pneumonia disease caused by Yersinia pestis. Semin. Respir. Infect. 12, 12-23. PubMed, Google Scholar - 43.
Zhao X, Skurnik M . 2016 Bacteriophages of Yersinia pestis. In Yersinia pestis: retrospective and perspective advances in experimental medicine and biology (edsYang R, Anisimov A ), pp. 361-375. Dordrecht, The Netherlands: Springer. Google Scholar - 44.
Carniel E . 2002 Plasmids and Pathogenicity Islands of Yersinia. In Pathogenicity islands and the evolution of pathogenic microbes. Current topics in microbiology and immunology, vol 264/2 (edsHacker J, Kaper JB ), pp. 89-108. Berlin, Germany: Springer. Google Scholar - 45.
Carniel E, Guilvout I, Prentice M . 1996 Characterization of a large chromosomal ‘high-pathogenicity island’ in biotype 1B Yersinia enterocolitica. J. Bacteriol. 178, 6743-6751. (doi:10.1128/jb.178.23.6743-6751.1996) Crossref, PubMed, Web of Science, Google Scholar - 46.
Schubert S, Rakin A, Heesemann J . 2004 The Yersinia high-pathogenicity island (HPI): evolutionary and functional aspects. Intl. J. Med. Microbiol. 294, 83-94. (doi:10.1016/j.ijmm.2004.06.026) Crossref, PubMed, Web of Science, Google Scholar - 47.
Marinaro M, di Tommaso A, Uzzau S, Fasano A, de Magistris MT . 1999 Zonula Occludens toxin is a powerful mucosal adjuvant for intranasally delivered antigens. Infect. Immun. 67, 1287-1291. (doi:10.1128/IAI.67.3.1287-1291.1999) Crossref, PubMed, Web of Science, Google Scholar - 48.
Karyekar CS, Fasano A, Raje S, Lu R, Dowling TC, Eddington ND . 2003 Zonula occludens toxin increases the permeability of molecular weight markers and chemotherapeutic agents across the bovine brain microvessel endothelial cells. J. Pharm. Sci. 92, 414-423. (doi:10.1002/jps.10310) Crossref, PubMed, Web of Science, Google Scholar - 49.
Groisman EA, Ochman H . 1996 Pathogenicity islands: bacterial evolution in quantum leaps. Cell. 87, 791-794. (doi:10.1016/S0092-8674(00)81985-6) Crossref, PubMed, Web of Science, Google Scholar - 50.
Prentice MB, Rahalison L . 2007 Plague. Lancet 369, 1196-1207. (doi:10.1016/S0140-6736(07)60566-2) Crossref, PubMed, Web of Science, Google Scholar - 51.
Lei C, Kumar S . 2022 Yersinia pestis antibiotic resistance: a systematic review. Osong Public Health Res. Perspect. 13, 24-36. (doi:10.24171/j.phrp.2021.0288) Crossref, PubMed, Google Scholar - 52.
Leslie T 2011 Outbreak of gastroenteritis caused by Yersinia pestis in Afghanistan. Epidemiol. Infect. 139, 728-735. (doi:10.1017/S0950268810001792) Crossref, PubMed, Web of Science, Google Scholar - 53.
Kjærgård A . 2006 Sejet Ødekirkegård, Gl. Sognevej 14b, Sejet. Horsens, Denmark: Horsens Museum. Google Scholar - 54.
Pedersen V . 1999 VSM F902: Sct. Trinitatis/Drotten, Sct. Leonis Gade 10, Viborg. Viborg, Denmark: Viborg Museum. Google Scholar - 55.
Kristensen HK . 1987 Middelalderbyen Viborg, pp. 47-56. Copenhagen, Denmark: SAGA Egmont. Google Scholar - 56.
Keyes G . 2007 VSM 09264: Sct. Trinitatis/Drotten, Sct. Leonis Gade 12, Viborg. Viborg, Denmark: Viborg Museum. Google Scholar - 57.
Knapp M, Clarke AC, Horsburgh KA, Matisoo-Smith EA . 2012 Setting the stage: building and working in an ancient DNA laboratory. Ann. Anat. 194, 3-6. (doi:10.1016/j.aanat.2011.03.008) Crossref, PubMed, Web of Science, Google Scholar - 58.
Pilli E 2013 Monitoring DNA contamination in handled vs. directly excavated ancient human skeletal remains. PLoS ONE 8, e52524. (doi:10.1371/journal.pone.0052524) Crossref, PubMed, Web of Science, Google Scholar - 59.
Yang DY, Watt K . 2005 Contamination controls when preparing archaeological remains for ancient DNA analysis. J. Archaeol. Sci. 32, 331-336. (doi:10.1016/j.jas.2004.09.008) Crossref, Web of Science, Google Scholar - 60.
Krause-Kyora B 2018 Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans. Nat. Commun. 9, 1569. (doi:10.1038/s41467-018-03857-x) Crossref, PubMed, Web of Science, Google Scholar - 61.
Peltzer A, Jäger G, Herbig A, Seitz A, Kniep C, Krause J, Nieselt K . 2016 EAGER: efficient ancient genome reconstruction. Genome Biol. 17, 60. (doi:10.1186/s13059-016-0918-z) Crossref, PubMed, Web of Science, Google Scholar - 62.
Li H, Durbin R . 2009 Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760. (doi:10.1093/bioinformatics/btp324) Crossref, PubMed, Web of Science, Google Scholar - 63.
Neukamm J, Peltzer A, Nieselt K . 2021 DamageProfiler: fast damage pattern calculation for ancient DNA. Bioinformatics 37, 3652-3653. (doi:10.1093/bioinformatics/btab190) Crossref, PubMed, Web of Science, Google Scholar - 64.
Korneliussen TS, Albrechtsen A, Nielsen R . 2014 ANGSD: analysis of next generation sequencing data. BMC Bioinf. 15, 356. (doi:10.1186/s12859-014-0356-4) Crossref, PubMed, Web of Science, Google Scholar - 65.
Renaud G, Slon V, Duggan AT, Kelso J . 2015 Schmutzi: estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA. Genome Biol. 16, 224. (doi:10.1186/s13059-015-0776-0) Crossref, PubMed, Web of Science, Google Scholar - 66.
Herbig A, Maixner F, Bos K, Zink A, Krause J, Huson D . 2016 MALT: Fast alignment and analysis of metagenomic DNA sequence data applied to the Tyrolean Iceman. bioRxiv. (doi:10.1101/050559) Google Scholar - 67.
Huson DH, Auch AF, Qi J, Schuster SC . 2007 MEGAN analysis of metagenomic data. Genome Res. 17, 377-386. (doi:10.1101/gr.5969107) Crossref, PubMed, Web of Science, Google Scholar - 68.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R . 2009 The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25, 2079. (doi:10.1093/bioinformatics/btp352) Crossref, Web of Science, Google Scholar - 69.
Jumper J 2021 Highly accurate protein structure prediction with AlphaFold. Nature 596, 583-589. (doi:10.1038/s41586-021-03819-2) Crossref, PubMed, Web of Science, Google Scholar - 70.
Edgar Robert C . 2004 MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5, 113. (doi:10.1186/1471-2105-5-113) Crossref, PubMed, Web of Science, Google Scholar - 71.
Stamatakis A . 2014 RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313. (doi:10.1093/bioinformatics/btu033) Crossref, PubMed, Web of Science, Google Scholar - 72.
Bos KI 2014 Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis. Nature 514, 494-497. (doi:10.1038/nature13591) Crossref, PubMed, Web of Science, Google Scholar - 73.
McKenna A 2010 The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303. (doi:10.1101/gr.107524.110) Crossref, PubMed, Web of Science, Google Scholar - 74.
Bonczarowska JH, Susat J, Krause-Kyora B, Dangvard Pedersen D, Boldsen J, Agersnap Larsen L, Seeberg L, Nebel A, Unterweger D . 2023Ancient Yersinia pestis genomes lack the virulence-associated ypfΦ prophage present in modern pandemic strains . Figshare. (doi:10.6084/m9.figshare.c.6723841) Google Scholar
