A 16th century Escherichia coli draft genome associated with an opportunistic bile infection

Aug 09, 2023

Communications Biology volume 5, Article number: 599 (2022) Cite this article

12k Accesses

557 Altmetric

Metrics details

Escherichia coli – one of the most characterized bacteria and a major public health concern – remains invisible across the temporal landscape. Here, we present the meticulous reconstruction of the first ancient E. coli genome from a 16th century gallstone from an Italian mummy with chronic cholecystitis. We isolated ancient DNA and reconstructed the ancient E. coli genome. It consisted of one chromosome of 4446 genes and two putative plasmids with 52 genes. The E. coli strain belonged to the phylogroup A and an exceptionally rare sequence type 4995. The type VI secretion system component genes appears to be horizontally acquired from Klebsiella aerogenes, however we could not identify any pathovar specific genes nor any acquired antibiotic resistances. A sepsis mouse assay showed that a closely related contemporary E. coli strain was avirulent. Our reconstruction of this ancient E. coli helps paint a more complete picture of the burden of opportunistic infections of the past.

The recovery of ancient pathogen DNA (aDNA) from human victims has almost exclusively focused on historically significant mortality events, such as the Black Death, revealing the evolutionary history of canonical pathogens such as Yersinia pestis, Mycobacterium tuberculosis1, and Variola virus. In contrast, much human morbidity and mortality is the result of opportunistic infections, ones that often remain invisible in the past. Opportunistic pathogens, those without historical records—such as Escherichia coli, Pseudomonas aeruginosa, and Staphylococcus aureus—have been understudied relative to their contemporary burdens on people today2. They are defined by their ability to infect during periods of stress, imbalance, or disturbance while being otherwise commensal or saprophytic3. Opportunistic infections likely played an important role in human mortality in our shared past and thus have had a broader impact on human health than can be or has been measured.

A confounding factor in the identification of historically understudied pathogens is that the resultant infections were likely primarily opportunistic. That is, they colonized their host environments asymptomatically, leaving no identifiable pathologies and as such are not easily identified in human remains4. Ancient DNA studies typically focus on obligatory pathogens such as M. leprae, Salmonella enterica, and Y. pestis that are easily correlated with pathologically distinct, or historically relevant mortality events and as such are easily distinguished as foreign in ancient human metagenomic DNA read data5,6,7. Opportunistic infections have the added burden of authentication due to their ubiquity as environmental contaminants and of modern commensal strains. Historical identification of these pathogens would allow for a careful assessment of their evolutionary history and the commensal-pathogen continuum defined by gene content gain and loss as strains modulate toward or away from host virulence.

E. coli is one such pathogen. It is a common commensal bacteria found in vertebrate gut microbiome8 that can also become an opportunistic pathogen under specific conditions9. E. coli has such a large impact on our health care systems that it is the subject of several vaccine development efforts to mitigate the effects of the most pathogenic strains10. Having been responsible for several food poisoning outbreaks and becoming a leading pathogen for deaths caused by antimicrobial resistance, E. coli is thus a key bacterium used in public health surveillance11.

Global sampling of E. coli strains produces a tree with several unique phylogenetic groups and interspersed pathovars12. While the phylogenetic relationships between strains are constructed based on genetic similarity, pathovars are defined by the virulence traits of their members. In many cases, members of the same pathovar do cluster together in the same clade, however, as virulence genes can be acquired horizontally, some pathovars—like enteroaggregative E. coli—are distributed across multiple phylogenetic groups12. This striking diversity and the evolutionary transitory states among E. coli strains highlight their genomic plasticity and evolutionary versatility along this aforementioned continuum. Ancient E. coli genomes would provide useful insights into the forces that influence the emergence of commensalism and pathogenicity in bacteria. Here, we describe the reconstruction of a 16th-century E. coli genome characterized from the gallstone of an Italian noble—Giovani d’Avalos (1538–1586)—highlighting a genome with commensal characteristics.

In 1983, the mummified remains of several Italian nobles were recovered from the Abbey of Saint Domenico Maggiore in Naples, Italy (Supplemental Fig. 1). A careful paleopathological and histological survey of one of the individuals—Giovani d’Avalos (NASD1), a Neapolitan noble who died in 1586 at the age of 48 (further information can be found in the Supplementary Notes section of the Supplemental Materials)—revealed thickened gallbladder walls, Rokitanski-Aschoff sinuses, and several intact gallstones13. These features suggested that NASD1 may have suffered from chronic cholecystitis14. While not the only cause of cholecystitis, chronic bacterial infections from E. coli, Bacteroides spp., and Clostridium spp. can lead to the formation of gallstones15. These infections are typically indicated by a brown pigment, as can be seen in the gallstone from NASD115 (Fig. 1a, b).

a Liver and gallbladder of Giovani d'Avalos. The gallstones can be seen in the red rectangle. Note its dark brown coloring. Scale bar represents 1 cm b Gallbladder with a thickened wall (a) and Rokitansky-Aschoff sinuses (b) (Hematoxylin-eosin, 3X and 250X). The scale bar for (a) represents 2 cm and 100 μm for (b). c Damage plots of the \({5}^{\prime}\) and \({3}^{\prime}\) ends of mapped reads for E.coli. Colors refer to the digest, not damage type. Mapdamage 2.052 was used to calculate the damage rates. d Fragment length distribution of deduplicated mapped reads from E. coli. A \({\log }_{10}\) scale is used to emphasize the differences between the digests. A minimum length of 30 bp was required for a read to be kept.

Metagenomic profiles generated from the DNA extracts using Kraken 216 provided preliminary evidence for the substantial and increasing presence of Enterobacteriales in digest rounds 2 (outer layer), 3–4 (middle layer), and 5–6 (inner layer) from the gallstone while indicating its absence in the tissue samples of NASD1 (Supplemental Fig. 2). Enterobacteriales reads are also present in digest round 1 however they are comparable with those in the reagent blanks (Supplemental Fig. 2), thus likely to be a mixture of endogenous and contaminant DNA such as Bradyrhizobium. Species-level indication suggested that E. coli made up the largest proportion of reads identified as Enterobacteriales. Fortunately, E. coli-specific reads are virtually absent (≤0.002%) from the blanks and present in the second through sixth digests. In comparison, other taxa such as Alcaligenes and Bradyrhizobium were present in both the experimental samples and the blanks and are likely to be contaminants. Human assigned reads were present in all digests but were also detected in the blanks. Klebsiella aerogenes was also detected, however only in digest rounds 3–4. To assess the authenticity of the Human, E. coli, and K. aerogenes reads we examined the deamination and depurination kinetics (Fig. 1c, d and Supplemental Fig. 3).

The reconstructed ancient E. coli genome had a mean read depth of 28.46[28.31, 28.62] × , mean heterozygosity of 7.66 × 10−4[6.66 × 10−4, 8.66 × 10−4], 4446 genes when mapped to the E. coli pan-genome dataset and a total estimated size of 4050, 429bp when the detected genes are concatenated. Significant differences (two-sided t-test, P = 7.951 × 10−05) were detected when comparing read depths of the core (3122 genes, 28.73[28.61, 28.86] ×) and accessory (1324genes, 27.82[27.39, 28.26] ×) genomes (Fig. 2a and Table 1); the same was true for their mean heterozygosity with the core at 2.74 × 10−4 [2.34 × 10−4, 3.15 × 10−4] and accessory genome with 1.92 × 10−3 [1.61 × 10−3, 2.24 × 10−3](Supplemental Fig. 4). Eighty-three genes were identified as having a substantially (at least two standard deviations) deeper gene coverage. These results were then compared to the closely related A phylogroup E. coli genomes of FSIS11816402, an sequence type (ST) 4995 strain, and K-12 MG1655, an ST10 strain (Fig. 2b and Supplemental Fig. 5). Both genome comparisons returned a lower mean depth with a larger variance of coverage. This is especially true for FSIS11816402 as 25.24% of its contigs had a mean coverage of less than 1 × (Supplemental Fig. 6).

a Distribution of mean gene coverages with a CV ≤1 for the ancient genome. The dashed line indicates the detection threshold at 10×. The black rectangle on the right indicates the area where genes with a high copy number (as defined by \(\bar{x}+2* s\)). b Coverage plot for FSIS11816402. A window of 1% was used for illustration purposes. c Coverage plot for CP019906.1. A window of 0.1% was used for illustration purposes. The first track indicates the coverage with the red line illustrating the overall mean. The second track indicates the number of SNPs over the same window while the third is the GC content. d Gene coverage of the T6SS in K. aerogenes using a 100 bp window. Gene names are included when available.

We identified two potential plasmids within the assembled scaffolds with MOB-suite17: E. coli MDR_56 plasmid unnamed 5 (CP019906.1) and S. flexneri 1a strain 0228 plasmids (CP012732.1). Subsequent mappings of the E. coli reads back to these reference plasmids confirmed their presence in the ancient strain and authenticated their origin (Fig. 2c and Supplemental Figs. 5, 7). One plasmid, CP019906.1 had a read coverage similar to that of K-12 (20.95[20.80, 21.10] ×) and 34 genes, whereas, CP012732.1 had lower coverage (16.46[16.28, 16.64] ×) and only 19 genes. These plasmids contain regions with no read coverage—including a ~6.6 Kbp region in CP012732.1—but we could not confirm if the genes were part of the chromosome.

A subset of the sequenced reads were also classified as belonging to K. aerogenes by Kraken 216. To confirm or refute this result, unmapped reads were aligned against the K. aerogenes (NC_015663.1) reference genome. The overall read depth (0.38[0.37, 0.38] ×) and genome coverage was low (3.43%), suggesting that only a subset of genes from K. aeorgenes was present. However, a 38 Kbp section of the genome (Fig. 2d) had a read mapping depth similar to the ancient E. coli results (20.36[20.23, 20.50] ×).

In order to place the ancient E. coli within the global phylogeny and help identify the strain, we produced a core SNP alignment of 451 E. coli and four Shigella genomes representing the current breadth of E. coli diversity. After removing redundancy, the resulting phylogeny contained 107 genomes. This pruned sample set was then used to create a new alignment of 5007 core SNPs which returned an overall topology resembling previously published results and placed the ancient strain into a strongly supported phylogeny within phylogroup A (Fig. 3a).

a The global phylogeny with bootstrap values and phylogenetic groups with E. coli EC42405 as the out-group. b Phylogeny of the reduced subgroup A0 as defined by the Clermont genotype (+ - - -)61. Tip points represent the sequence type of the strain. IAI1 was the out-group of the phylogeny. c Phylogeny of the ST4995 strains. Red rectangles represent the 95% confidence interval for the topology; labels indicate the median date of divergence. The evolutionary rate for the phylogeny was 2.555 × 10−6[1.567 × 10−6, 3.992 × 10−6] subs/site/year. CFSAN051544 was the out-group. The brown arrows indicate the position of the ancient genome.

To more carefully refine the position of the ancient strain within group A, we generated a reduced ML tree using 94 genomes (Fig. 3b) consisting of 291 SNPs. The ancient genome clustered closely with the ST4995 strains with some statistical support (bootstrap support of 65%), signifying that it is likely part of the same sequence type (Table 2). An even further refined ML tree was then generated using the 22 available ST4995 genomes from Enterobase18 (Fig. 3c). The core SNP alignment contained 16,026 nucleotides, providing a much greater resolution than the global and A0 subgroup alignments. The ancient strain clearly clusters within the ST4995 genomes.

We scanned for temporal signals across all phylogenies using TEMPest19 and found one when the phylogeny was restricted to the ST4995 genomes (Supplemental Fig. 8), a result confirmed a date-randomization test in LSD (P = 0.003). The ST4995 phylogeny (Fig. 3c) suggests an 11th century (1027[787, 1220]CE) tMCRA for ST4995 strains, nearly 500 years before the diversification of the modern strains. The estimated evolutionary rate across the phylogeny was 2.555 × 10−6[1.567 × 10−6, 3.992 × 10−6] subs/site/year, similar to previously published results20.

Multilocus sequence and typing of our ancient E. coli genome confirmed its placement in the ST4995—a rare sequence type with only 22 strains in the Enterobase18 database which encompasses 182,476 entries18. It exhibited an Onovel15:H? serotype, similar to its sister taxa (see the Supplementary Discussion section in the Supplementary Materials for full results). A different fimH variant – fimH86 – was present in our E. coli strain.

The majority (95%) of the E. coli strains contained 3190 core genes, 3122 of which were detected in our ancient genome (Supplemental Fig. 9). Enrichment analysis of the missing core genes revealed a significant lack of flagellar assembly genes (P = 0.0445), however, the protein-protein network contained the expected number of interactions (P = 0.11). An analysis of regions with mapped ancient reads in CP019906.1 revealed the presence of genes involved in biofilm formation (NP_311382.1) and a multidrug efflux pump (acrD). For CP012732.1, in contrast, reads mapped to pseudogenes are likely related to osmoprotection and transcription regulations.

A presence/absence (P/A) analysis of the accessory genome recapitulates the results from the ML phylogeny results (Fig. 4a) and previously published results21. The accessory PCoA indicates that the genome is a member of phylogroup A and an ancestral member of ST4995. The latter is confirmed by a P/A analysis of an ST4995-only pan-genome (Supplemental Fig. 10).

The accessory genome (a, c) and identified virulence genes (b, d) were clustered using a binary distance. Phylogroups were identified using Clermont Typing59 while pathovars were determined using available metadata from the Patric49 database. The ancient genome is indicated by the black arrow. In both binary distances were calculated prior to create the PCoA.

A total of 91 virulence factors were identified in the ancient genome, 37 of which were not found in E. coli K-12 MG1655 (see Table 3 for gene families with more than one hit). Type VI secretion system (T6SS) components consisted of the majority of these genes and contained high copy numbers. The T6SS— formed by the tss gene family as well as hcp, vasK, and vgrG—mediates antagonistic interactions between competing bacteria. In addition to its role in bacterial killing, T6SS is involved in interbacterial signaling, biofilm formation, and phage defense22. Parts of a type III secretion system were also detected at much lower levels.

We also observed incomplete virulence gene complements involved in several mechanisms. Several genes belonging to the ecp and fim fimbriae families were found in the ancient genome. While they are involved in virulence9,23,24, they can also be present in commensal strains. Fimbriae known to be associated with virulence such as cfaB and csg were also identified23 yet they do not provide enough information to pinpoint the specific pathovar of the ancient strain9,23,24. The key virulence factors for the E. coli pathovars such as eae and stx for Shiga toxin-producing E. coli (STEC), heat-labile and heat-stable genes for Enterotoxigenic E. coli (ETEC), and the type III secretion system in EPEC and EHEC were absent9.

In contrast astA, a key virulence factor in Enteroaggregative E. coli (EAEC)25, was detected in the ancient genome. In combination with evidence that T6SS is commonly found in EAEC strains9, this may indicate that the ancient strain was been a member of this pathovar. The remaining EAEC virulence factors in the pan-genome—aggA, aggR, and aap—however, were not detected25. The ancient genome also contains a complete enterobactin system which is responsible for iron uptake23. The enterobactin system and astA gene were also found in E. coli K-12 MG1655, once again indicating that they do not confer enough virulence on their own.

To evaluate the extraintestinal virulence potential of ST4995 we tested a close relative (507 core SNPs separating the two strains) of the ancient genome, the reference strain ATCC11229 (see Fig. 3b and Supplemental Fig. 11 for similarities to the modern strain) in a well-calibrated mouse sepsis assay where ten mice were inoculated with the strain and their death recorded26. None of the mice were killed by ATCC11229, whereas all of the mice infected by the positive control B2 phylogroup CFT073 strain were killed. The phenotype of the strain which killed the mice is linked to the presence of specific extraintestinal virulence genes with a major role in the iron capture system encoding genes as those bared by the high-pathogenicity island (HPI)27,28. This result is in agreement with the absence of typical extraintestinal virulence genes present in the ATCC11229 and ancient strain.

We searched our ancient genome for the presence of antimicrobial resistance (AMR) genes using the Resistance Gene Identifier29 and found 47 distinct genes in the scaffold assembly of our ancient E. coli genome. Of these, 35 genes were also identified in the global pan-genome P/A analysis. Eight genes contained duplicates with mdtB being the most commonly detected. The remaining subset spanned thirteen antimicrobial drug resistance classes with nine classes represented more than once (Table 4). Resistance cassettes targeting the five most common drug classes were detected in the ancient E. coli strain. The majority of these genes are multidrug efflux pumps, which is typical of E. coli30. There were no unexpected AMR genes present in our ancient genome with all 35 being found in E. coli K-12 MG1655.

The DNA isolated from the stone showed clear evidence of DNA damage. Specifically, the deamination plots for the ancient reads mapping to the E. coli pan-genome reference contained characteristic aDNA profiles which indicated increased deamination rates at or near the terminal bases (Fig. 1c). In addition depurination kinetics plots of the E. coli mapped reads from libraries with sufficient E. coli read depth (digests 2, 3–4, 5–6) showed that fragments from internal sections of the gallstone (better protected from hydrolysis) were on average longer (34 bp median in digest 2 to 37 bp median in digests 5–6) (Fig. 1d). We obtained similar results for the K. aerogenes data across digests two through six. Interestingly, the mapped human reads lack these signature features (Supplemental Fig. 3). More importantly, however, is that they confirm that the E. coli and K. aerogenes sequences are indeed ancient and were not the result of modern contamination.

Two potential explanations exist for the lack of deamination signal in the mapped human reads. The first is that the sample likely contains a mixture of both contaminating and to a lesser degree authentic ancient human DNA. This is almost certainly the case for the first digest, as it represents the outer layer of the gallstone, which would have been subjected to contamination from handling, storage, sampling, etc. However, as we moved to deeper layers of the gallstone, presumably better protected from contamination, the lack of endogenous human DNA was surprising. This is especially striking with the edit distances for both human and E. coli reads from across the different digests (Supplemental Fig. 12). The more likely explanation is that there are simply not enough human DNA reads for meaningful damage analysis. Gallstones are typically formed from a combination of cholesterol, bile salts, and phosphatidylcholine15. Furthermore, in the highly acidic environment of the gallbladder, only bacteria that are specialized in surviving these conditions can thrive and potentially form a stone15. Since gallstones are not formed from direct human cellular components its interior is likely virtually devoid of human DNA, which is very different from the results we obtained from the DNA of ancient abscesses31.

The gene depth results of the ancient E. coli genome confirm the lack of modern contaminants as they are part of a unimodal normal distribution. The sole exception to these coverage stats are 83 genes with substantially larger mean depths (Supplemental Data 2). These genes consist of primarily hypothetical proteins (24.10%), intergrases and transposons (19.28%), and the T6SS (13.25%). They also contain a significantly (P = 5.03 × 10−11) greater mean heterozygosity (1.23 × 10−2[9.25 × 10−3, 1.53 × 10−2]) than the rest of the ancient genome. These high copy number genes, however, are mainly accessory genes (96.34%) which suggests—in combination with the greater heterozygosity and presence of transposons—that they were located in mobile elements shared between multiple bacterial species.

The plasmids we detected in the scaffold assembly, CP019906.1 and CP012732.1, had lower gene depth coverages than that for our reconstructed chromosome (Table 1). Most plasmids are typically present in greater copy numbers than chromosomes32, however, this was not the case as the ancient chromosome had a significant (P = 0 using Tukey's method to compare coverage estimates between the plasmids and ancient chromosome) deeper mean coverage than CP019906.1 and CP012732.1. One possibility is that these two plasmids are actually not present, but rather the genes to which the ancient reads mapped were actually located in the ancient chromosome. This is supported by the fact that the mean gene depths for the plasmids fall firmly within the coverage distribution of the chromosome (Supplemental Table 2). This could explain large regions of the plasmids that lack any mapped reads, as well as the lack of reads mapping to intergenic regions (see Fig. 2c and Supplementary Fig. 5B).

The reads mapping to the K. aeorgenes genome could represent a similar situation to the plasmid data. Interestingly, the 38 Kbp region in the K. aerogenes chromosome which had a read depth similar yet significantly different (P ≈ 0 using Tukey's method to compare coverage estimates between the K. aerogenes T6SS and the ancient E. coli chromosome) from the plasmids and E. coli chromosome (20.36 ± 13.56 × ) contains the Type VI Secretion System (T6SS), which is ubiquitous across Gram-negative bacteria and plays an important role in antagonistic interactions22. Given that T6SS is widespread among bacteria, its presence is not an exclusive marker for K. aerogenes. However, competitive read mapping between E. coli and K. aerogenes does suggest that this T6SS does indeed belong to K. aerogenes rather than the homologous system in E. coli suggesting a recent horizontal transfer, which is supported by a comparison of their relative GC contents33 (Supplemental Fig. 13).

The structure of our global phylogeny resembles those previously published34 and places the ancient genome in phylogroup A with meaningful support (100% bootstrap for a clade of three strains). Interestingly, other historical (~1800s) E. coli strains have also been placed in this phylogroup35. Strains that have been isolated from humans living in less industrialized and more rural communities typically harbor strains falling within phylogroup A8. The late Medieval date of our sample, along with its phylogenetic placement helps confirm its authenticity. In addition, a study of E. coli strains isolated from bile infections showed that they belong mainly to the phylogenetic A group36.

We identified a temporal signal in the reduced phylogeny (Fig. 3c) which provides additional context for the evolutionary history of ST4995 strains. The majority of the sequence type can trace their lineage to a strain from approximately 1592[1479, 1690] CE, however, the divergence of our ancient genome and two others – ESC_AA9618AA_AS and ESC_VA4573AA_AS – has a deeper tMRCA ca. 956[689, 1172] CE. Given the disparity in ages between the two clades, it is likely that ST4995 is actually composed of at least two subgroups. Evidence supporting this can be seen in an MDS of the P/A analysis of all existing ST4995 genomes (Supplementary Fig. 9).

The P/A analysis of the accessory genes recapitulates the results from the ML phylogenies (Fig. 4a). The accessory PCoA confirms that the ancient strain is a member of phylogenetic group A while a subsequent P/A analysis with an ST4995-only pan-genome supports the proposition that the genome is a member of ST4995 (Supplemental Fig. 9). Phylogenetic group A—the group from which ST4995 stems—is a well-described clade and is the most commonly isolated type from human commensals8. E. coli strain diversity is driven primarily by socioeconomic factors as individuals living in modern industrialized countries are more likely to carry B2 strains while those in less industrialized and rural communities harbor mostly A8. Given the time period NASD1 lived, our ancient E. coli is a member of phylogenetic group A.

No acquired AMR genes were detected in the ancient genome confirming the age of the strain. Post removal of the multidrug efflux pumps from our genome, no other AMR genes were found when compared to E. coli K-12 MG165530. The latter is extremely pertinent as strong resistance to a particular drug class requires several mechanisms37. Thus, it is likely that only mild resistances to these antimicrobial compounds were needed30,38.

The identified ancient genome likely exhibited commensal characteristics as it shares most of its virulence factors with E. coli K-12 MG1655. The ancient strain lacks the canonical virulence factors for STEC and ETEC while also missing the components of the type III secretion system used in EHEC. Components of EAEC—astA and the T6SS9—were detected, yet other key virulence factors were missing and astA was also detected in K-12 MG1655. Results from the P/A analysis of the accessory genes with pathovar labels (Fig. 4c) and the virulence genes (Fig. 4d) present a similar picture. The ancient genome is positioned in a relatively general cluster containing a variety of different pathovars which include EHEC, ETEC, ExPEC, and STEC. The mouse sepsis model further reinforces that the ancient genome was an opportunistic pathogen as the proxy strain—ATCC 11229—did not infect any mice. This information, in combination with the identified virulence genes and P/A clustering of the ancient E. coli suggests that the ancient genome is an opportunistic pathogen that acquired a K. aerogenes T6SS gene cassette during expansion in the gallbladder, consistent with previous work showing Klebsiella sp often associated with E. coli in contemporary bile infections39.

We performed DNA extractions on a single gallstone—that had been isolated from the mummified remains of NASD1—in the clean-room facilities of the McMaster Ancient DNA Centre (McMaster University, Hamilton, Ontario, Canada) using methodologies specific to ancient DNA. We hypothesized that the successive mineralized layers of the stone would contain increasingly better-preserved DNA and thus subjected the gallstone to six successive rounds of demineralization and digestion40. Subsequently, DNA was extracted from 250 μL of the supernatant stemming from each round, which was then purified and prepared into dual-indexed libraries, size-selected to remove artifacts, and sequenced on the Illumina HiSeq 1500 platform with 90 bp paired-end reads41,42. Three additional tissue types from the same individual (bladder, small intestine, and lung) were processed and sequenced as above for comparative purposes. Full details of the methods can be found in the Supplementary Methods section in the Supplementary Materials.

Adapters were identified through AdapterRemoval43 with trimming and merging performed by leeHom44 using its aDNA setting. Sequencing lanes were pooled when applicable and Kraken 216 determined the overall metagenomic composition of the samples. This was done to characterize the metagenome and identify taxa of interest. A standard Kraken 216 database with some modifications was used—a kmer size of 25 bp and no minimizer—to account for the smaller read lengths of aDNA fragments. Reads underwent string deduplication using prinseq45 prior to the metagenomic analysis.

Samples were then mapped against the human genome GRCh3846 with BWA47 with a maximum edit distance of 0.01 (-n 0.01), at most two gap openings (-o 2), and seeding effectively disabled (-l 16500)48. A minimum length of 30 bp and mapping quality of 30 was required for a read to be successfully matched. Fragments that did not successfully map were then compared with an E. coli pan-genome using the same settings. The pan-genome was created with E. coli genomes from PATRIC49 and four Shigella genomes from NCBI. Specifically, genomes which were found in humans, chickens, cows, dogs, mice, or pigs and were determined to have a "good" genome quality (as defined by PATRIC) were selected resulting in 451 strains. These sequences were annotated with Prokka50 using the proteins from E. coli K-12 MG1655 (NC_000913.3) obtained from NCBI as trusted sequences. The resulting annotations were compiled by Roary51 to create the pan-genome. Paralogs were not split and a minimum blastP identity of 90% were used as additional settings. Unsuccessfully mapped human reads were also compared with K. aerogenes (NC_015663.1) as it was substantially— proportional abundance ≥1%—present in one of the digests while also being excluded from the blanks.

Successfully mapped human, E. coli, and K. aerogenes reads were deduplicated based on their \({5}^{\prime}\) and \({3}^{\prime}\) coordinates with bam-rmdup (https://bitbucket.org/ustenzel/biohazard-tools/src/master/). The deduplicated samples were pooled together based on their digest of origin (1, 2, 3–4, or 5–6). MapDamage 2.052 was used to estimate the amount of deamination across the mapped reads. Fragment length distributions (FLD) and mapping mismatches were also extracted. Heterozygosity was tested by calling single nucleotide polymorphisms (SNPs) using bcftools with ploidy set to two53 and a quality threshold of 100.

Following authentication, an assembly was created with the complete set of mapped E. coli and K. aerogenes reads (i.e., not deduplicated) using SPAdes 3.14.1 in isolate mode54 with customized kmer lengths—9, 19, and 29. The resulting assembly was an incomplete and highly fragmented scaffold. This was mostly due to the short median fragment lengths (35 bp) of the input reads and the required kmer sizes. Nevertheless, the N50 of the assembly was 939 bp with a total length of 3, 917, 339 bp which is substantially shorter than the complete E. coli chromosomes in PATRIC (5.00[4.96, 5.05]Mbp). MOB-suite was used to identify potential plasmids in the assembly, however, the final identification was done by mapping to the reference plasmids17.

We created an SNP alignment using the E. coli genomes from PATRIC and the Shigella genomes using Snippy (https://github.com/tseemann/snippy) with E. coli EC42405 (CP043414.1) as the reference strain. Non-canonical base pairs were replaced with an N and recombinant regions were removed with gubbins55. As the ancient genome contained several gaps the "filter_percentage" flag was increased to 32% to ensure that the assembly was not excluded. Once recombination sites were identified and removed a maximum likelihood phylogeny was created using IQ-TREE 256. ModelFinder was used to select the appropriate phylogenetic model with ascertainment bias57 and the tree was bootstrapped 1000 times. Redundant genomes were removed from the phylogeny using Treemmer making sure to retain 95% of the diversity—measured by Root-to-Tip length –, which returned 107 genomes58. A new SNP alignment and phylogeny was then created using these pruned genomes. The modern strains were labeled based on known E. coli phylogenetic groups using ClermonTyping59. Four E. coli strains which were identified as part of A were relabeled as potential Shigella strains. This is due to the presence of an ipaH gene in two of the genomes and their positioning in the phylogeny60.

Upon determining the phylogroup of the ancient genome, the subgroup A0 E. coli genomes were used to carefully position the ancient strain within the same clade. This subgroup is defined by the Clermont genotype (+ - - -)61. An additional seven enterotoxigenic E. coli strains obtained from von Mentzer et al. 2014 were also included in the phylogeny62. This phylogeny was created using the same process noted previously and was rooted using IAI1 B1 (NC_011741.1). The phylogeny of the ST4995 strains was created with the 22 ST4995 genomes available in Enterobase18 and rooted using the ST325 strain CFSAN051544 (SAMN05414627). To determine if the data contained any temporal signature, sampling dates or sequencing dates were assigned to the pruned phylogeny and a Root-to-tip (RTT) regression was performed for the global, A0 subgroup, and ST4995 phylogenies. The presence of a temporal signal was then assessed via a date-randomization test and with TEMPest19. A molecular clock was fit to the data using a least-squares dating approach using LSD263,64.

The dated ST4995 phylogeny was created by dating the bootstrap trees with LSD2 (using -r l n-s 1500000 as settings)63. These trees were then summarized using treeannotator65 with the height of the nodes set by the mean of the bootstrap distribution. The resulting phylogeny was then used to plot the dated phylogeny with the node height confidence intervals included.

To determine what genes were present in the ancient strain, we converted the read depths of the ancient E. coli pan-genome into a gene presence/absence (P/A) metric. We then used a conservative approach, identifying a gene as present if it had an average sequence depth of at least 10 × with a coefficient of variation (CV) ≤1. A CV filter was used to remove genes, which contained substantial regions of stacked reads. In comparison to a percent coverage threshold, a CV is more permissive with unmapped areas if the rest of the gene coverage is relatively consistent. An additional threshold of \(\bar{x}+2* s\) was used to identify genes with a high copy number. While still part of the ancient genome these genes could potentially belong to plasmids32 or represent gene amplification66.

The P/A matrix for modern strains was generated by Roary51. A 95% core genome of the modern strains was identified and genes missing from our ancient strain were submitted to STRING67 to be functionally annotated using a network consisting of "high confidence" interactions. Virulence factors were identified using a curated set of known virulence genes68. Ambiguously named genes (i.e., group_1000) were run through blastX69 against the RefSeq non-redundant protein database from 2019-11-29 with a maximum E-value of 10−5 to determine candidates. If multiple hits were found for a gene the smallest E-value was selected. Ties were resolved by choosing the gene with the highest bitscore followed by the highest percent identity. A minimum percent identity of 90% was required for a match to be kept and the ambiguous gene renamed. The descriptions of the previously ambiguous genes were searched for the keywords "secretion, invasion, and enterotoxin". Genes with hits from this search were included with the virulence genes. To determine the pathovar of the ancient genome, relevant information was extracted from the PATRIC metadata.

We performed a P/A analysis on the core, accessory genomes, and the virulence genes. The strains were transformed into a binary distance matrix and clustered using a PCoA. To ensure that redundant genes were not included, two exclusion criteria were created: genes that were ubiquitously present were excluded, and genes not present in at least five genomes were removed from the accessory genome analysis. Of the 451 E. coli genomes, four were removed from the pan-genome analysis due to low gene counts (see Supplemental Fig. 14 for gene counts and Supplemental Data 1 for removed genomes). The Resistance Gene Identifier (RGI)29 with default settings identified potential antimicrobial resistances in the ancient strain. Only RGI entries returned due to gene homologies and present in the pan-genome were accepted. Both the identified virulence and resistance genes were compared to the archetypal E. coli strain K-12 MG1655. It is avirulent and susceptible to antibiotics70.

Ten 14–16 g (4-week-old) female mice OF1 from Charles River® (L’Arbresle, France) per strain received a subcutaneous injection of 0.2 ml of bacterial suspension in the neck (2 × 108 colony-forming unit). Time to death was recorded during the following 7 days. Mice surviving more than 7 days were considered cured and sacrificed. The E. coli CFT073 strain was used as a positive control killing all the inoculated mice whereas the E. coli K-12 MG1655 strain was used as a negative control for which all the inoculated mice survive26. The protocol (n∘APAFIS#4948) was approved by the French Ministry of Research and by the ethical committee for animal experiments. The ATCC11229 strain (AMC 198) was obtained from the Institut Pasteur collection (CIP 103795) and whole-genome sequenced using Illumina technology to verify the identity of the strain. This strain is of commensal origin and is commonly used in bacterial resistance testing.

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

The raw sequencing data has been uploaded to NCBI as PRJNA810725 with supporting metadata located in Supplementary Table 1. Access to the sample itself can be arranged by contacting G.S.L.

The scripts and additional data used to analyze the data after pan-genome mapping are available at https://github.com/longg2/AncientEcoli71. The scripts used for the initial mapping and pan-genome creation is available at https://github.com/longg2/LongBioinformatics/tree/vAncEcoli72.

Duchêne, S., Ho, S. Y., Carmichael, A. G., Holmes, E. C. & Poinar, H. The recovery, interpretation and use of ancient pathogen genomes. Curr. Biol. 30, R1215–R1231 (2020).

Article PubMed PubMed Central CAS Google Scholar

Parkins, M. D., Somayaji, R. & Waters, V. J. Epidemiology, biology, and impact of clonal Pseudomonas aeruginosa infections in cystic fibrosis. Clin. Microbiol. Rev. 31, e00019-18 (2018).

Achermann, Y., Goldstein, E. J., Coenye, T. & Shirtliff, M. E. Propionibacterium acnes: from commensal to opportunistic biofilm-associated implant pathogen. Clin. Microbiol. Rev. 27, 419–440 (2014).

Article PubMed PubMed Central Google Scholar

Wood, J. W. et al. The osteological paradox: problems of inferring prehistoric health from skeletal samples. Curr. Anthropol. 33, 343–370 (1992).

Article Google Scholar

Krause-Kyora, B. et al. Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans. Nat. Commun. 9, 1–11 (2018).

Article CAS Google Scholar

Wu, X. et al. A 3,000-year-old, basal S. enterica lineage from Bronze Age Xinjiang suggests spread along the Proto-Silk Road. PLoS Pathog. 17, e1009886 (2021).

Article CAS PubMed PubMed Central Google Scholar

Wagner, D. M. et al. Yersinia pestis and the Plague of Justinian 541-543 AD: a genomic analysis. Lancet Infect. Dis. 14, 319–326 (2014).

Article PubMed Google Scholar

Tenaillon, O., Skurnik, D., Picard, B. & Denamur, E. The population genetics of commensal Escherichia coli. Nat. Rev. Microbiol. 8, 207–217 (2010).

Article CAS PubMed Google Scholar

Croxen, M. A. & Finlay, B. B. Molecular mechanisms of Escherichia coli pathogenicity. Nat. Rev. Microbiol. 8, 26–38 (2010).

Article CAS PubMed Google Scholar

Rojas-Lopez, M., Monterio, R., Pizza, M., Desvaux, M. & Rosini, R. Intestinal pathogenic Escherichia coli: insights for vaccine development. Front. Microbiol. 9, 440 (2018).

Article PubMed PubMed Central Google Scholar

Murray, C. J. et al. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet 399, 629–655 (2022).

Denamur, E., Clermont, O., Bonacorsi, S. & Gordon, D. The population genetics of pathogenic Escherichia coli. Nat. Rev. Microbiol. 19, 37–54 (2021).

Article CAS PubMed Google Scholar

Fornaciari, G., Pollina, L., Tornaboni, D. & Tognetti, A. Pulmonary and hepatic pathologies in the series of mummies of s. domenico maggiore at naples (xvi century). In Proc. VII European Meeting of the Paleopathology Association (Lyon, September 1988) 89–92. Marino Solfanelli Editore (Journal of Paleopathology, Monographic Publications, 1989).

Knab, L. M., Boller, A.-M. & Mahvi, D. M. Cholecystitis. Surg. Clin. North Am. 94, 455–470 (2014).

Article PubMed Google Scholar

Venneman, N. G. & van Erpecum, K. J. Pathogenesis of gallstones. Gastroenterolog. Clin. North Am. 39, 171–183 (2010).

Article Google Scholar

Wood, D. E., Lu, J. & Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 20, 257 (2019).

Article CAS PubMed PubMed Central Google Scholar

Robertson, J. & Nash, J. H. MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb. Genom. 4, e000206 (2018).

Zhou, Z. et al. The EnteroBase user's guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity. Genome Res. 30, 138–152 (2020).

Article CAS PubMed PubMed Central Google Scholar

Rambaut, A., Lam, T. T., Carvalho, L. M. & Pybus, O. G. Pybus Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2, vew007 (2016).

Lee, H., Popodi, E., Tang, H. & Foster, P. L. Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing. Proc. Natl Acad. Sci. 109, E2774–E2783 (2012).

Article CAS PubMed PubMed Central Google Scholar

Touchon, M. et al. Phylogenetic background and habitat drive the genetic diversification of Escherichia coli. PLoS Genet. 16, e1008866 (2020).

Article CAS PubMed PubMed Central Google Scholar

Russell, A. B., Peterson, S. B. & Mougous, J. D. Type VI secretion system effectors: poisons with a purpose. Nat. Rev. Microbiol. 12, 137–148 (2014).

Article CAS PubMed PubMed Central Google Scholar

Liu, B., Zheng, D., Jin, Q., Chen, L. & Yang, J. VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 47, D687–D692 (2019).

Article CAS PubMed Google Scholar

Morales-Espinosa, R. et al. UPEC strain characterization isolated from Mexican patients with recurrent urinary infections. J. Infect. Dev. Ctries. 10, 317–328 (2016).

Article CAS PubMed Google Scholar

Aslani, M. M., Alikhani, M. Y., Zavari, A., Yousefi, R. & Zamani, A. R. Characterization of enteroaggregative Escherichia coli (EAEC) clinical isolates and their antibiotic resistance pattern. Int. J. Infect. Dis. 15, e136–e139 (2011).

Article CAS PubMed Google Scholar

Johnson, J. R. et al. Experimental mouse lethality of Escherichia coli isolates, in relation to accessory traits, phylogenetic group, and ecological source. J. Infect. Dis. 194, 1141–1150 (2006).

Article CAS PubMed Google Scholar

Johnson, J. R. et al. Accessory traits and phylogenetic background predict Escherichia coli extraintestinal virulence better than does ecological source. J. Infect. Dis. 219, 121–132 (2019).

CAS PubMed Google Scholar

Galardini, M. et al. Major role of iron uptake systems in the intrinsic extra-intestinal virulence of the genus Escherichia revealed by a genome-wide association study. PLoS Genet. 16, e1009065 (2020).

Article CAS PubMed PubMed Central Google Scholar

Alcock, B. P. et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 48, D517–D525 (2020).

CAS PubMed Google Scholar

Goldstone, R. J. & Smith, D. G. A population genomics approach to exploiting the accessory ’resistome’ of Escherichia coli. Microb. Genom. 3, e000108 (2017).

Devault, A. M. et al. A molecular portrait of maternal sepsis from byzantine troy. eLife 6, e20983 (2017).

Article PubMed PubMed Central Google Scholar

Friehs, K. Plasmid copy number and plasmid stability. Adv. Bioch. Eng. Biotechnol. 86, 47–82 (2004).

Daubin, V., Lerat, E. & Perrière, G. The source of laterally transferred genes in bacterial genomes. Genome Biol. 4, 1–12 (2003).

Article Google Scholar

Chaudhuri, R. R. & Henderson, I. R. The evolution of the Escherichia coli phylogeny. Infect. Genet. Evol. 12, 214–226 (2012).

Article PubMed Google Scholar

Desroches, M. et al. The odyssey of the ancestral Escherich strain through culture collections: an example of allopatric diversification. mSphere 3, e00553-17 (2018).

Wang, M.-C. et al. Different roles of host and bacterial factors in Escherichia coli extra-intestinal infections. Clin. Microbiol. Infect. 15, 372–379 (2009).

Article PubMed Google Scholar

Benveniste, R. & Davies, J. Mechanisms of antibiotic resistance in bacteria. Ann. Rev. Biochem. 42, 471–506 (1973).

Article CAS PubMed Google Scholar

Piddock, L. J. Multidrug-resistance efflux pumps? not just for resistance. Nat. Rev. Microbiol. 4, 629–636 (2006).

Article CAS PubMed Google Scholar

Brook, I. Aerobic and anaerobic microbiology of biliary tract disease. J. Clin. Microbiol. 27, 2373–2375 (1989).

Article CAS PubMed PubMed Central Google Scholar

Schwarz, C. et al. New insights from old bones: DNA preservation and degradation in permafrost preserved mammoth remains. Nucleic Acids Res. 37, 3215–3229 (2009).

Article CAS PubMed PubMed Central Google Scholar

Dabney, J. et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort dna fragments. Proc. Natl Acad Sci. USA 110, 15758–15763 (2013).

Glocke, I. & Meyer, M. Extending the spectrum of dna sequences retrieved from ancient bones and teeth. Genome Res. 27, 1230–1237 (2017).

Schubert, M., Lindgreen, S. & Orlando, L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res. Notes 9, 88 (2016).

Article PubMed PubMed Central Google Scholar

Renaud, G., Stenzel, U. & Kelso, J. leeHom: adaptor trimming and merging for Illumina sequencing reads. Nucleic Acids Res. 42, e141–e141 (2014).

Article PubMed PubMed Central CAS Google Scholar

Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864 (2011).

Article CAS PubMed PubMed Central Google Scholar

Schneider, V. A. et al. Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly. Genome Res. 27, 849–864 (2017).

Article CAS PubMed PubMed Central Google Scholar

Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).

Article CAS PubMed PubMed Central Google Scholar

Duggan, A. T. et al. 17th century variola virus reveals the recent history of smallpox. Curr. Biol. 26, 3407–3412 (2016).

Article CAS PubMed PubMed Central Google Scholar

Wattam, A. R. et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 45, D535–D542 (2016).

Article PubMed PubMed Central CAS Google Scholar

Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).

Article CAS PubMed Google Scholar

Page, A. J. et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693 (2015).

Article CAS PubMed PubMed Central Google Scholar

Jónsson, H., Ginolhac, A., Schubert, M., Johnson, P. L. & Orlando, L. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 29, 1682–1684 (2013).

Article PubMed PubMed Central CAS Google Scholar

Seguin-Orlando, A. et al. No particular genomic features underpin the dramatic economic consequences of 17th century plague epidemics in italy. iScience 24, 102383 (2021).

Article CAS PubMed PubMed Central Google Scholar

Nurk, S. et al. Assembling genomes and mini-metagenomes from highly chimeric reads. In Annual International Conference on Research in Computational Molecular Biology (Springer, 2013).

Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using gubbins. Nucleic Acids Res. 43, e15–e15 (2015).

Article PubMed CAS Google Scholar

Minh, B. Q. et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534 (2020).

Article CAS PubMed PubMed Central Google Scholar

Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K., Von Haeseler, A. & Jermiin, L. S. Modelfinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589 (2017).

Article CAS PubMed PubMed Central Google Scholar

Menardo, F. et al. Treemmer: a tool to reduce large phylogenetic datasets with minimal loss of diversity. BMC Bioinformatics 19, 1–8 (2018).

Article Google Scholar

Beghain, J., Bridier-Nahmias, A., Le Nagard, H., Denamur, E. & Clermont, O. Clermontyping: an easy-to-use and accurate in silico method for escherichia genus strain phylotyping. Microb. Genom. 4, e000192 (2018).

Hazen, T. H. et al. Investigating the relatedness of enteroinvasive Escherichia coli to other E. coli and Shigella isolates by using comparative genomics. Infect. Immun. 84, 2362–2371 (2016).

Clermont, O., Christenson, J. K., Denamur, E. & Gordon, D. M. The Clermont Escherichia coli phylo-typing method revisited: improvement of specificity and detection of new phylo-groups. Environ. Microbiol. Rep. 5, 58–65 (2013).

Article CAS PubMed Google Scholar

von Mentzer, A. et al. Identification of enterotoxigenic Escherichia coli (ETEC) clades with long-term global distribution. Nat. Genet. 46, 1321 (2014).

Article CAS Google Scholar

To, T.-H., Jung, M., Lycett, S. & Gascuel, O. Fast dating using least-squares criteria and algorithms. Syst. Biol. 65, 82–97 (2016).

Article CAS PubMed Google Scholar

Duchene, S. et al. Inferring demographic parameters in bacterial genomic data using bayesian and hybrid phylogenetic methods. BMC Evol. Biol. 18, 1–11 (2018).

Article Google Scholar

Suchard, M. A. et al. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4, vey016 (2018).

Article PubMed PubMed Central Google Scholar

Riehle, M. M., Bennett, A. F. & Long, A. D. Genetic architecture of thermal adaptation in Escherichia coli. Proc. Natl Acad. Sci. USA 98, 525–530 (2001).

Article CAS PubMed PubMed Central Google Scholar

Szklarczyk, D. et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613 (2019).

Article CAS PubMed Google Scholar

de Lastours, V. et al. Mortality in Escherichia coli bloodstream infections: antibiotic resistance still does not make it. J. Antimicrob. Chemother. 75, 2334–2343 (2020).

Article PubMed CAS Google Scholar

Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinformatics 10, 421 (2009).

Article PubMed PubMed Central CAS Google Scholar

Touchon, M. et al. Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths. PLoS Genet. 5, e1000344 (2009).

Article PubMed PubMed Central CAS Google Scholar

longg2. longg2/AncientEcoli: published version. Zenodo https://doi.org/10.5281/zenodo.6403348 (2022).

longg2. longg2/LongBioinformatics: ancient Ecoli version. Zenodo https://doi.org/10.5281/zenodo.6403351 (2022).

Download references

We would like to thank Wael Elhenawy and Brian Coombes for their correspondence as well as Sara Dion for her technical assistance with the mouse model. Funding for this work was supported by an SSHRC Insight Grant, research funding from CIFAR's human and the microbiome program and the Boris Family Fund to H.P. as well as the NSERC grant RGPIN-2020-05733 awarded to G.B.G. This work was also partially supported by the "Fondation pour la Recherche Médicale" Equipe FRM 2016, grant number DEQ20161136698 to E.D.

These authors contributed equally: George S. Long, Jennifer Klunk.

Department of Biology, McMaster University, Hamilton, Canada

George S. Long, Jennifer Klunk & G. Brian Golding

McMaster Ancient DNA Centre, Departments of Anthropology and Biochemistry, McMaster University, Hamilton, Canada

George S. Long, Jennifer Klunk, Ana T. Duggan, Madeline Tapson & Hendrik Poinar

Daicel Arbor Biosciences, 5840 Interface Drive, Suite 101, Ann Arbor, MI, 48103, USA

Jennifer Klunk & Madeline Tapson

Department of Anthropology, McMaster University, Hamilton, ON, Canada

Ana T. Duggan

Division of Paleopathology, Department of Translational Research and New Technologies in Medicine and Surgery, University of Pisa, Via Roma 57, 56126, Pisa, Italy

Valentina Giuffra

Department of Human Science (DISUM), University of Catania, Piazza Dante 32, 95124, Catania, Italy

Lavinia Gazzè

Department of Civilisations and Forms of Knowledge, University of Pisa, Via Trieste 40, 56126, Pisa, Italy

Antonio Fornaciari & Gino Fornaciari

Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC, Australia

Sebastian Duchene

Université Paris Cité, IAME, UMR 1137, INSERM, 75018, Paris, France

Olivier Clermont & Erick Denamur

Laboratoire de Génétique Moléculaire, Hôpital Bichat, APHP, 75018, Paris, France

Erick Denamur

Michael G. DeGroote Institute for Infectious Disease Research and CIFAR Humans and the Microbiome Program, Toronto, ON, Canada

Hendrik Poinar

CIFAR Humans and the Microbiome Program, Toronto, ON, Canada

Hendrik Poinar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

G.S.L. and J.K. were responsible for the investigation, methodology, formal analysis, and original draft of this study. G.S.L. also performed the data curation. M.T., A.T.D., V.G., L.G., and A.F. aided in the investigation. A.T.D. also contributed to the methodology. G.F. provided the sample. E.D. and O.C. validated the E. coli function and typing results. E.D. also provided funding for the mouse sepsis assay. S.D. validated the temporal signal and aided in its investigation. G.B.G. supervised and provided computational resources and funding. H.P. conceived and administered the project, provided funding, and also supervised. All authors approved the final manuscript.

Correspondence to George S. Long, Erick Denamur or Hendrik Poinar.

The authors declare no competing interests.

Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editor: Luke R. Grinham. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Long, G.S., Klunk, J., Duggan, A.T. et al. A 16th century Escherichia coli draft genome associated with an opportunistic bile infection. Commun Biol 5, 599 (2022). https://doi.org/10.1038/s42003-022-03527-1

Download citation

Received: 01 March 2022

Accepted: 23 May 2022

Published: 16 June 2022

DOI: https://doi.org/10.1038/s42003-022-03527-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.