Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Impressive pan-genomic diversity of E. coli from a wild animal community near urban development reflects human impacts

Human and domesticated animal waste infiltrates global freshwater, terrestrial, and marine environments, widely disseminating fecal microbes, antibiotics, and other chemical pollutants. Emerging evidence suggests that guts of wild animals are being invaded by our microbes, including Escherichia coli, which face anthropogenic selective pressures to gain antimicrobial resistance (AMR) and increase virulence. However, wild animal sources remain starkly under-represented among genomic sequence repositories. We sequenced whole genomes of 145 E. coli isolates from 55 wild and 13 domestic animal fecal samples, averaging 2 (ranging 1-7) isolates per sample, on a preserve imbedded in a human-dominated landscape in California Bay Area, USA, to assess AMR, virulence, and pan-genomic diversity. With single nucleotide polymorphism analyses we predict potential transmission routes. We illustrate the usefulness of E. coli to aid our understanding of and ability to surveil the emergence of zoonotic pathogens created by the mixing of human and wild bacteria in the environment.

Detection and characterization of potentially hybrid enteroaggregative Escherichia coli (EAEC) strains isolated from urinary tract infection

Uropathogenic Escherichia coli (UPEC) have the potential to receive the virulence markers of intestinal pathotypes and transform into various important hybrid pathotypes. This study aimed to investigate the frequency and characteristics of hybrid enteroaggregative E. coli (EAEC)/UPEC strains. Out of 202 UPEC strains, nine (4.5%) were detected as hybrid EAEC/UPEC. These strains carried one to four iron uptake systems. Among nine investigated pathogenicity islands (PAIs), PAI IV536, PAI II536, and PAI ICFT073 were found in 9 (100%), 3 (33.3%), and 1 (11.1%) strains, respectively. The chuA and sitA genes were detected in 5 (55.5%) and 3 (33.3%) hybrid strains, respectively. Six hybrid strains were found to be typical extraintestinal pathogenic E. coli (ExPEC) according to their virulence traits. Most of the hybrid strains belonged to the phylogenetic group E (6/9). Among the hybrid strains, seven (7/9) were able to form biofilm and adhere to cells; however, only two strains penetrated into the HeLa cells. Our findings reveal some of the virulence characteristics of hybrid strains that lead to fitness and infection in the urinary tract. These strains, with virulence factors of intestinal and non-intestinal pathotypes, may become emerging pathogens in clinical settings; therefore, further studies are needed to reveal their pathogenicity mechanisms and so that preventive measures can be taken.

Systematic review and meta-analyses of the role of drinking water sources in the environmental dissemination of antibiotic-resistant Escherichia coli in Africa

Escherichia coli are pathogenic and antibiotic-resistant organisms that can spread to humans through water. However, there is sparse synthesised information on the dissemination of antibiotic-resistant E. coli through drinking water in Africa. This review provides an overview of the environmental spread of antimicrobial-resistant E. coli through drinking water in Africa. We performed a systematic review based on PRISMA guidelines, and 40 eligible studies from 12 countries were identified until June 2023. Four electronic databases (PubMed, Elsevier, AJOL, and DOAJ) were searched. Studies that employed phenotypic tests (n = 24/40) in identifying the bacterium outstripped those that utilised genome-based methods (n = 13). Of the 40 studies, nine and five, respectively, assessed the bacterium for antimicrobial resistance (AMR) phenotype and genotype. Multiple antibiotic resistance indices of 0.04-0.1 revealed a low level of antibiotic resistance. The detection of multidrug-resistant E. coli carrying resistance genes in certain water sources suggests that AMR-surveillance expansion should include drinking water.

Diversity of antimicrobial-resistant bacteria isolated from Australian chicken and pork meat

Antimicrobial-resistant bacteria are frequently isolated from retail meat and may infect humans. To determine the diversity of antimicrobial-resistant bacteria in Australian retail meat, bacteria were cultured on selective media from raw chicken (n = 244) and pork (n = 160) meat samples obtained from all four major supermarket chains in the ACT/NSW, Australia, between March and June 2021. Antimicrobial susceptibility testing (AST) was performed for 13 critically and 4 highly important antibiotics as categorised by the World Health Organization (WHO) for a wide range of species detected in the meat samples. A total of 288 isolates underwent whole-genome sequencing (WGS) to identify the presence of antimicrobial resistance (AMR) genes, virulence genes, and plasmids. AST testing revealed that 35/288 (12%) of the isolates were found to be multidrug-resistant (MDR). Using WGS data, 232/288 (81%) of the isolates were found to harbour resistance genes for critically or highly important antibiotics. This study reveals a greater diversity of AMR genes in bacteria isolated from retail meat in Australia than previous studies have shown, emphasising the importance of monitoring AMR in not only foodborne pathogenic bacteria, but other species that are capable of transferring AMR genes to pathogenic bacteria.

Diverse bacteriophages for biocontrol of ESBL- and AmpC-β-lactamase-producing E. coli

Novel solutions are needed to reduce the risk of transmission of extended spectrum β-lactamase (ESBL) and AmpC β-lactamase producing Escherichia coli (ESBL/AmpC E. coli) from livestock to humans. Given that phages are promising biocontrol agents, a collection of 28 phages that infect ESBL/AmpC E. coli were established. Whole genome sequencing showed that all these phages were unique and could be assigned to 15 different genera. Host range analysis showed that 82% of 198 strains, representing the genetic diversity of ESBL/AmpC E. coli, were sensitive to at least one phage. Identifying receptors used for phage binding experimentally as well as in silico predictions, allowed us to combine phages into two different cocktails with broad host range targeting diverse receptors. These phage cocktails efficiently inhibit the growth of ESBL/AmpC E. coli in vitro, thus suggesting the potential of phages as promising biocontrol agents.

Diversification of gene content in the Mycobacterium tuberculosis complex is determined by phylogenetic and ecological signatures

We analyzed the pan-genome and gene content modulation of the most diverse genome data set of the Mycobacterium tuberculosis complex (MTBC) gathered to date. The closed pan-genome of the MTBC was characterized by reduced accessory and strain-specific genomes, compatible with its clonal nature. However, significantly fewer gene families were shared between MTBC genomes as their phylogenetic distance increased. This effect was only observed in inter-species comparisons, not within-species, which suggests that species-specific ecological characteristics are associated with changes in gene content. Gene loss, resulting from genomic deletions and pseudogenization, was found to drive the variation in gene content. This gene erosion differed among MTBC species and lineages, even within M. tuberculosis, where L2 showed more gene loss than L4. We also show that phylogenetic proximity is not always a good proxy for gene content relatedness in the MTBC, as the gene repertoire of Mycobacterium africanum L6 deviated from its expected phylogenetic niche conservatism. Gene disruptions of virulence factors, represented by pseudogene annotations, are mostly not conserved, being poor predictors of MTBC ecotypes. Each MTBC ecotype carries its own accessory genome, likely influenced by distinct selective pressures such as host and geography. It is important to investigate how gene loss confer new adaptive traits to MTBC strains; the detected heterogeneous gene loss poses a significant challenge in elucidating genetic factors responsible for the diverse phenotypes observed in the MTBC. By detailing specific gene losses, our study serves as a resource for researchers studying the MTBC phenotypes and their immune evasion strategies.IMPORTANCEIn this study, we analyzed the gene content of different ecotypes of the Mycobacterium tuberculosis complex (MTBC), the pathogens of tuberculosis. We found that changes in their gene content are associated with their ecological features, such as host preference. Gene loss was identified as the primary driver of these changes, which can vary even among different strains of the same ecotype. Our study also revealed that the gene content relatedness of these bacteria does not always mirror their evolutionary relationships. In addition, some genes of virulence can be variably lost among strains of the same MTBC ecotype, likely helping them to evade the immune system. Overall, our study highlights the importance of understanding how gene loss can lead to new adaptations in these bacteria and how different selective pressures may influence their genetic makeup.

Occurrence of Imipenem-Resistant Uropathogenic Escherichia coli in Pregnant Women: An Insight into Their Virulence Profile and Clonal Structure

Uropathogenic Escherichia coli (UPEC) is the predominant pathogen in Urinary Tract Infection (UTI) in pregnant and non-pregnant women. Limited studies were initiated to explore UPEC from pregnant women with respect to imipenem resistance, pathogenicity, and their clonal lineage. In this study, imipenem resistance, phylogenetic background, virulence-associated genes, and clonal characteristics in UPECs isolated from pregnant and non-pregnant cohorts were investigated. E. coli was identified biochemically from urine culture-positive samples from pregnant and non-pregnant women. Carbapenem (meropenem, ertapenem, imipenem) susceptibility was determined by Kirby-Bauer disk diffusion test. The pathogenic determinants were identified by PCR. MEGA 11 was used to interpret clonal lineages from MLST. GraphPad Prism 8.0 and SPSS 26.0 were used for statistical interpretation. Results indicated highest resistance against imipenem compared to meropenem and ertapenem in UPECs isolated from pregnant (UPECp; 63.89%) and non-pregnant (UPECnp; 87.88%) women. Although phylogroup E was predominant in both imipenem-resistant isolates, acquisition of virulence factors was higher among UPECnp than UPECp. Akin to this observation, the presence of PAI III536 and PAI IV536 was statistically significant (p < 0.05) in the former. MLST analysis revealed similar clonal lineages between UPECnp and UPECp, which showed an overall occurrence of ST405 followed by ST101, ST410, ST131, and ST1195 in UPECnp and ST167 in UPECp, respectively, with frequent occurrence of CC131, CC405. Therefore, imipenem-resistant UPECp although discrete with respect to their virulence determinants when compared to UPECnp shared similar STs and CCs, which implied common evolutionary history. Thus, empiric treatment must be restricted in UTIs to especially protect maternal and fetal health.

Genome-wide screen of genetic determinants that govern Escherichia coli growth and persistence in lake water

Although enteric bacteria normally reside within the animal intestine, the ability to persist extraintestinally is an essential part of their overall lifestyle, and it might contribute to transmission between hosts. Despite this potential importance, few genetic determinants of extraintestinal growth and survival have been identified, even for the best-studied model, Escherichia coli. In this work, we thus used a genome-wide library of barcoded transposon insertions to systematically identify functional clusters of genes that are crucial for E. coli fitness in lake water. Our results revealed that inactivation of pathways involved in maintaining outer membrane integrity, nucleotide biosynthesis, and chemotaxis negatively affected E. coli growth or survival in this extraintestinal environment. In contrast, inactivation of another group of genes apparently benefited E. coli growth or persistence in filtered lake water, resulting in higher abundance of these mutants. This group included rpoS, which encodes the general stress response sigma factor, as well as genes encoding several other global transcriptional regulators and RNA chaperones, along with several poorly annotated genes. Based on this co-enrichment, we identified these gene products as novel positive regulators of RpoS activity. We further observed that, despite their enhanced growth, E. coli mutants with inactive RpoS had reduced viability in lake water, and they were not enriched in the presence of the autochthonous microbiota. This highlights the duality of the general stress response pathway for E. coli growth outside the host.

Vaccines against extraintestinal pathogenic Escherichia coli (ExPEC): progress and challenges

The emergence of antimicrobial resistance (AMR) is a principal global health crisis projected to cause 10 million deaths annually worldwide by 2050. While the Gram-negative bacteria Escherichia coli is commonly found as a commensal microbe in the human gut, some strains are dangerously pathogenic, contributing to the highest AMR-associated mortality. Strains of E. coli that can translocate from the gastrointestinal tract to distal sites, called extraintestinal E. coli (ExPEC), are particularly problematic and predominantly afflict women, the elderly, and immunocompromised populations. Despite nearly 40 years of clinical trials, there is still no vaccine against ExPEC. One reason for this is the remarkable diversity in the ExPEC pangenome across pathotypes, clades, and strains, with hundreds of genes associated with pathogenesis including toxins, adhesins, and nutrient acquisition systems. Further, ExPEC is intimately associated with human mucosal surfaces and has evolved creative strategies to avoid the immune system. This review summarizes previous and ongoing preclinical and clinical ExPEC vaccine research efforts to help identify key gaps in knowledge and remaining challenges.

Scaling of Protein Function across the Tree of Life

Scaling laws are a powerful way to compare genomes because they put all organisms onto a single curve and reveal nontrivial generalities as genomes change in size. The abundance of functional categories across genomes has previously been found to show power law scaling with respect to the total number of functional categories, suggesting that universal constraints shape genomic category abundance. Here, we look across the tree of life to understand how genome evolution may be related to functional scaling. We revisit previous observations of functional genome scaling with an expanded taxonomy by analyzing 3,726 bacterial, 220 archaeal, and 79 unicellular eukaryotic genomes. We find that for some functional classes, scaling is best described by multiple exponents, revealing previously unobserved shifts in scaling as genome-encoded protein annotations increase or decrease. Furthermore, we find that scaling varies between phyletic groups at both the domain and phyla levels and is less universal than previously thought. This variability in functional scaling is not related to taxonomic phylogeny resolved at the phyla level, suggesting that differences in cell plan or physiology outweigh broad patterns of taxonomic evolution. Since genomes are maintained and replicated by the functional proteins encoded by them, these results point to functional degeneracy between taxonomic groups and unique evolutionary trajectories toward these. We also find that individual phyla frequently span scaling exponents of functional classes, revealing that individual clades can move across scaling exponents. Together, our results reveal unique shifts in functions across the tree of life and highlight that as genomes grow or shrink, proteins of various functions may be added or lost.

Promoter capture drives the emergence of proto-genes in Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli Long-Term Evolution Experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time-span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, thereby serving as raw material for new gene emergence. Most proto-genes result either from insertion element activity or chromosomal translocations that fused pre-existing regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, persist stably, and can serve as potential substrates for new gene formation.

Semantic search using protein large language models detects class II microcins in bacterial genomes

Class II microcins are antimicrobial peptides that have shown some potential as novel antibiotics. However, to date only ten class II microcins have been described, and discovery of novel microcins has been hampered by their short length and high sequence divergence. Here, we ask if we can use numerical embeddings generated by protein large language models to detect microcins in bacterial genome assemblies and whether this method can outperform sequence-based methods such as BLAST. We find that embeddings detect known class II microcins much more reliably than does BLAST and that any two microcins tend to have a small distance in embedding space even though they typically are highly diverged at the sequence level. In datasets of Escherichia coli , Klebsiella spp., and Enterobacter spp. genomes, we further find novel putative microcins that were previously missed by sequence-based search methods.

Bacterial genome-wide association study substantiates papGII of Escherichia coli as a major risk factor for urosepsis

BACKGROUND

Urinary tract infections (UTIs) are among the most common bacterial infections worldwide, often caused by uropathogenic Escherichia coli. Multiple bacterial virulence factors or patient characteristics have been linked separately to progressive, more invasive infections. In this study, we aim to identify pathogen- and patient-specific factors that drive the progression to urosepsis by jointly analysing bacterial and host characteristics.

METHODS

We analysed 1076 E. coli strains isolated from 825 clinical cases with UTI and/or bacteraemia by whole-genome sequencing (Illumina). Sequence types (STs) were determined via srst2 and capsule loci via fastKaptive. We compared the isolates from urine and blood to confirm clonality. Furthermore, we performed a bacterial genome-wide association study (bGWAS) (pyseer) using bacteraemia as the primary clinical outcome. Clinical data were collected by an electronic patient chart review. We concurrently analysed the association of the most significant bGWAS hit and important patient characteristics with the clinical endpoint bacteraemia using a generalised linear model (GLM). Finally, we designed qPCR primers and probes to detect papGII-positive E. coli strains and prospectively screened E. coli from urine samples (n = 1657) at two healthcare centres.

RESULTS

Our patient cohort had a median age of 75.3 years (range: 18.00-103.1) and was predominantly female (574/825, 69.6%). The bacterial phylogroups B2 (60.6%; 500/825) and D (16.6%; 137/825), which are associated with extraintestinal infections, represent the majority of the strains in our collection, many of which encode a polysaccharide capsule (63.4%; 525/825). The most frequently observed STs were ST131 (12.7%; 105/825), ST69 (11.0%; 91/825), and ST73 (10.2%; 84/825). Of interest, in 12.3% (13/106) of cases, the E. coli pairs in urine and blood were only distantly related. In line with previous bGWAS studies, we identified the gene papGII (p-value < 0.001), which encodes the adhesin subunit of the E. coli P-pilus, to be associated with 'bacteraemia' in our bGWAS. In our GLM, correcting for patient characteristics, papGII remained highly significant (odds ratio = 5.27, 95% confidence interval = [3.48, 7.97], p-value < 0.001). An independent cohort of cases which we screened for papGII-carrying E. coli at two healthcare centres further confirmed the increased relative frequency of papGII-positive strains causing invasive infection, compared to papGII-negative strains (p-value = 0.033, chi-squared test).

CONCLUSIONS

This study builds on previous work linking papGII with invasive infection by showing that it is a major risk factor for progression from UTI to bacteraemia that has diagnostic potential.

Circulation of enterotoxigenic Escherichia coli (ETEC) isolates expressing CS23 from the environment to clinical settings

IMPORTANCE

The importance of clean water cannot be overstated. It is a vital resource for maintaining health and well-being. Unfortunately, water sources contaminated with fecal discharges from animal and human origin due to a lack of wastewater management pose a significant risk to communities, as they can become a means of transmission of pathogenic bacteria like enterotoxigenic E. coli (ETEC). ETEC is frequently found in polluted water in countries with a high prevalence of diarrheal diseases, such as Bolivia. This study provides novel insights into the circulation of ETEC between diarrheal cases and polluted water sources in areas with high rates of diarrheal disease. These findings highlight the Choqueyapu River as a potential reservoir for emerging pathogens carrying antibiotic-resistance genes, making it a crucial area for monitoring and intervention. Furthermore, the results demonstrate the feasibility of a low-cost, high-throughput method for tracking bacterial pathogens in low- and middle-income countries, making it a valuable tool for One Health monitoring efforts.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

Stress Resistance and Virulence Gene Profiles Associated with Phylogeny and Phenotypes of Escherichia coli from Cattle

Seven serogroups of E. coli (Top seven E. coli) are frequently implicated in foodborne outbreaks in North America, largely due to their carriage of Shiga toxin genes (stx). This study aimed to profile resistance genes and virulence factors (VF), and their potential association with phylogeny and phenotypes of Top seven E. coli originating from cattle in Canada. 155 Top seven E. coli isolates previously characterized for heat and acid resistance and biofilm-forming ability were whole-genome sequenced and analyzed for phylogeny, VF, and stress resistance genes. The 155 E. coli strains belonged to six phylogroups: A (n = 32), B1 (n = 93), C (n = 3), D (n = 11), E (n = 15), and G (n = 1). Different phylogroups were clearly separated on the core genome tree, with strains of the same serotype closely clustered. The carriage of stx and the transmissible locus of stress tolerance (tLST), the extreme heat resistance marker, was mutually exclusive, in 33 and 15 genomes, respectively. A novel O84:H2 strain carrying stx1a was also identified. In total, 70, 41, and 32 VF, stress resistance genes and antibiotic resistance genes were identified. The stress resistance genes included those for metal (n = 29), biocides/acid (n = 4), and heat (n = 8) resistance. All heat resistance genes and most metal-resistance genes that were differentially distributed among the phylogroups were exclusively in phylogroup A. VF were least and most present in phylogroups A and D, respectively. No specific genes associated with acid resistance or biofilm formation phenotypes were identified. VF were more abundant (P < 0.05) in the non-biofilm-forming population and acid-resistant population.

Comparative Genomic Analyses of Escherichia coli from a Meat Processing Environment in Relation to Their Biofilm Formation and Persistence

We investigated the phylogeny of biofilm forming (BF) and nonbiofilm forming (NBF) Escherichia coli (n = 114) from a beef processing environment as well as genetic elements in their BF and persistence via a comparative genomic analysis. Phylogroup B1 made up the largest proportion of both the BF (73.8%) and NBF (50.9%) groups. E. coli from all of the sources that were examined had mixed phylogroups, except for those that were recovered from equipment after cleaning, which were exclusively from phylogroup B1. Both the core genome and gene content trees showed a tree-wide spread of BF strains, with clusters, including both BF and NBF strains. Genome-wide association studies (GWAS) via Scoary or Pyseer did not find any genes or mutations that were overrepresented in the BF group. A retrospective analysis of phenotypes found a significant correlation (P < 0.05) between BF ability and curli production, cellulose synthesis, and/or mobility. However, the BF group also included strains that were negative for curli and cellulose and/or missing encoding genes for the two traits. All curli and cellulose encoding genes were present in most genomes, regardless of their BF status. The degree of motility was correlated with both curli and cellulose production, and 80 common genes were overrepresented in all three of the trait-positive groups. A PTS enzyme II, a subsidiary gluconate catabolism pathway, and an iron-dicitrate transport system were more abundant in the persisting E. coli group. These findings suggest gene function redundancy in E. coli for biofilm formation as well as additional substrate utilization and iron acquisition in its persistence. IMPORTANCE The persistence of potentially hazardous bacteria is a major challenge for meat processing environments, which are conducive for biofilm formation. Marker genes/phenotypes are commonly used to differentiate biofilm forming E. coli strains from their nonbiofilm forming counterparts. We took a comparative genomic analysis approach to analyze E. coli strains that were from the same environment but were differentiated by their biofilm forming ability. A diversification of the genes involved in the biofilm formation of E. coli was observed. Even though there is a correlation on the population level between biofilm formation and the expression of curli and cellulose, uncertainties exist on the individual strain level. Novel substrate utilization and iron acquisition could contribute to the persistence of E. coli. These findings not only advance our understanding of the ecology of E. coli with respect to its persistence but also show that a marker gene/phenotype driven approach for the biofilm control of E. coli may not be prudent.

PanGraph: scalable bacterial pan-genome graph construction

The genomic diversity of microbes is commonly parameterized as SNPs relative to a reference genome of a well-characterized, but arbitrary, isolate. However, any reference genome contains only a fraction of the microbial pangenome, the total set of genes observed in a given species. Reference-based approaches are thus blind to the dynamics of the accessory genome, as well as variation within gene order and copy number. With the widespread usage of long-read sequencing, the number of high-quality, complete genome assemblies has increased dramatically. In addition to pangenomic approaches that focus on the variation in the sets of genes present in different genomes, complete assemblies allow investigations of the evolution of genome structure and gene order. This latter problem, however, is computationally demanding with few tools available that shed light on these dynamics. Here, we present PanGraph, a Julia-based library and command line interface for aligning whole genomes into a graph. Each genome is represented as a path along vertices, which in turn encapsulate homologous multiple sequence alignments. The resultant data structure succinctly summarizes population-level nucleotide and structural polymorphisms and can be exported into several common formats for either downstream analysis or immediate visualization.

Under-Appreciated Phylogroup Diversity of Escherichia coli within and between Animals at the Urban-Wildland Interface

Wild animals have been implicated as reservoirs and even "melting pots" of pathogenic and antimicrobial-resistant bacteria of concern to human health. Though Escherichia coli is common among vertebrate guts and plays a role in the propagation of such genetic information, few studies have explored its diversity beyond humans nor the ecological factors that influence its diversity and distribution in wild animals. We characterized an average of 20 E. coli isolates per scat sample (n = 84) from a community of 14 wild and 3 domestic species. The phylogeny of E. coli comprises 8 phylogroups that are differentially associated with pathogenicity and antibiotic resistance, and we uncovered all of them in one small biological preserve surrounded by intense human activity. Challenging previous assumptions that a single isolate is representative of within-host phylogroup diversity, 57% of individual animals sampled carried multiple phylogroups simultaneously. Host species' phylogroup richness saturated at different levels across species and encapsulated vast within-sample and within-species variation, indicating that distribution patterns are influenced both by isolation source and laboratory sampling depth. Using ecological methods that ensure statistical relevance, we identify trends in phylogroup prevalence associated with host and environmental factors. The vast genetic diversity and broad distribution of E. coli in wildlife populations has implications for biodiversity conservation, agriculture, and public health, as well as for gauging unknown risks at the urban-wildland interface. We propose critical directions for future studies of the "wild side" of E. coli that will expand our understanding of its ecology and evolution beyond the human environment. IMPORTANCE To our knowledge, neither the phylogroup diversity of E. coli within individual wild animals nor that within an interacting multispecies community have previously been assessed. In doing so, we uncovered the globally known phylogroup diversity from an animal community on a preserve imbedded in a human-dominated landscape. We revealed that the phylogroup composition in domestic animals differed greatly from that in their wild counterparts, implying potential human impacts on the domestic animal gut. Significantly, many wild individuals hosted multiple phylogroups simultaneously, indicating the potential for strain-mixing and zoonotic spillback, especially as human encroachment into wildlands increases in the Anthropocene. We reason that due to extensive anthropogenic environmental contamination, wildlife is increasingly exposed to our waste, including E. coli and antibiotics. The gaps in the ecological and evolutionary understanding of E. coli thus necessitate a significant uptick in research to better understand human impacts on wildlife and the risk for zoonotic pathogen emergence.

High prevalence of β-lactam and fluoroquinolone resistance in various phylotypes of Escherichia coli isolates from urinary tract infections in Jiroft city, Iran

BACKGROUND

Urinary tract infection (UTI) is one of the most prevalent infectious diseases with worldwide health threatening. Antimicrobial resistant strains of Escherichia coli (E. coli) are a common cause of UTI which were identified as a treatment challenge. This study aimed to assay the prevalence of common β-lactam resistance genes including bla_TEM, bla_SHV, bla_CTX-M and bla_CMY and phenotypic resistance to commonly used β-lactam and fluoroquinolone antibiotics in UTIs. These factors were evaluated in various phylogenetic groups (phylotypes) of E. coli isolates. Real-time PCR was applied to detect β-lactam resistance genes and conventional PCR was used to determine the phylotypes. Phenotypic resistance against β-lactams (ceftazidime, cefotaxime, aztreonam and ceftriaxone) and fluoroquinolones (ciprofloxacin) were identified by the disc diffusion technique. The ability of extended spectrum β-lactamases (ESBLs) production in E. coli isolates was detected using the combined disc diffusion method.

RESULTS

The prevalence of resistance genes were 89.6% for bla_TEM, 44.3% for bla_CTX-M, 6.6% for bla_SHV and 0.9% for bla_CMY. The two high prevalent phylotypes were B2 (29.2%) and D (17.9%) followed by E (14.1%), F (9.4%), C (6.6%) and 10.3% of isolates were unknown in phylotyping. Disc diffusion results showed high prevalence of antibiotic resistance to cefotaxime (88.6%), aztreonam (83%), ceftireaxon (77.3%), ceftazidime (76.4%) and ciprofloxacin (55.6%). Totally, 52.8% of isolates were found as phenotypical ESBL-producers.

CONCLUSIONS

This study's results confirmed an explosion of antibiotic resistance amongst E. coli isolates from UTI against β-lactams and fluoroquinolones. Findings explain the necessity of deep changes in quantity and quality of drug resistance diagnosis and antibiotic therapy strategies. More studies are suggested to better and confident evaluations.

Horizontal gene transfer among host-associated microbes

Horizontal gene transfer is an important evolutionary force, facilitating bacterial diversity. It is thought to be pervasive in host-associated microbiomes, where bacterial densities are high and mobile elements are frequent. These genetic exchanges are also key for the rapid dissemination of antibiotic resistance. Here, we review recent studies that have greatly extended our knowledge of the mechanisms underlying horizontal gene transfer, the ecological complexities of a network of interactions involving bacteria and their mobile elements, and the effect of host physiology on the rates of genetic exchanges. Furthermore, we discuss other, fundamental challenges in detecting and quantifying genetic exchanges in vivo, and how studies have contributed to start overcoming these challenges. We highlight the importance of integrating novel computational approaches and theoretical models with experimental methods where multiple strains and transfer elements are studied, both in vivo and in controlled conditions that mimic the intricacies of host-associated environments.

Genomic diversity of non-diarrheagenic fecal Escherichia coli from children in sub-Saharan Africa and south Asia and their relatedness to diarrheagenic E. coli

Escherichia coli is a frequent member of the healthy human gastrointestinal microbiota, as well as an important human pathogen. Previous studies have focused on the genomic diversity of the pathogenic E. coli and much remains unknown about the non-diarrheagenic E. coli residing in the human gut, particularly among young children in low and middle income countries. Also, gaining additional insight into non-diarrheagenic E. coli is important for understanding gut health as non-diarrheagenic E. coli can prevent infection by diarrheagenic bacteria. In this study we examine the genomic diversity of non-diarrheagenic fecal E. coli from male and female children with or without diarrhea from countries in sub-Saharan Africa and south Asia as part of the Global Enteric Multicenter Study (GEMS). We find that these E. coli exhibit considerable genetic diversity as they were identified in all E. coli phylogroups and an Escherichia cryptic clade. Although these fecal E. coli lack the characteristic virulence factors of diarrheagenic E. coli pathotypes, many exhibit remarkable genomic similarity to previously described diarrheagenic isolates with differences attributed to mobile elements. This raises an important question of whether these non-diarrheagenic fecal E. coli may have at one time possessed the mobile element-encoded virulence factors of diarrheagenic pathotypes or may have the potential to acquire these virulence factors.

A Representative Collection of Commensal Extended-Spectrum- and AmpC-β-Lactamase-Producing Escherichia coli of Animal Origin for Phage Sensitivity Studies

Introduction

Extended-spectrum β-lactamase (ESBL)- and AmpC β-lactamase (AmpC)-producing Escherichia coli from livestock and meat represent a zoonotic risk and biocontrol solutions are needed to prevent transmission to humans.

Methods

In this study, we established a representative collection of animal-origin ESBL/AmpC E. coli as target to test the antimicrobial potential of bacteriophages.

Results

Bioinformatic analysis of whole-genome sequence data of 198 ESBL/AmpC E. coli from pigs, broilers, and broiler meat identified strains belonging to all known E. coli phylogroups and 65 multilocus sequence types. Various ESBL/AmpC genes and plasmid types were detected with expected source-specific patterns. Plaque assay using 15 phages previously isolated using the E. coli reference collection demonstrated that Warwickvirus phages showed the broadest host range, killing up to 26 strains.

Conclusions

154/198 strains were resistant to infection by all phages tested, suggesting a need for isolating phages specific for ESBL/AmpC E. coli. The strain collection described in this study is a useful resource fulfilling such need.

Migration Rates on Swim Plates Vary between Escherichia coli Soil Isolates: Differences Are Associated with Variants in Metabolic Genes

This study investigates migration phenotypes of 265 Escherichia coli soil isolates from the Buffalo River basin in Minnesota, USA. Migration rates on semisolid tryptone swim plates ranged from nonmotile to 190% of the migration rate of a highly motile E. coli K-12 strain. The nonmotile isolate, LGE0550, had mutations in flagellar and chemotaxis genes, including two IS3 elements in the flagellin-encoding gene fliC. A genome-wide association study (GWAS), associating the migration rates with genetic variants in specific genes, yielded two metabolic variants (rygD-serA and metR-metE) with previous implications in chemotaxis. As a novel way of confirming GWAS results, we used minimal medium swim plates to confirm the associations. Other variants in metabolic genes and genes that are associated with biofilm were positively or negatively associated with migration rates. A determination of growth phenotypes on Biolog EcoPlates yielded differential growth for the 10 tested isolates on d-malic acid, putrescine, and d-xylose, all of which are important in the soil environment. IMPORTANCE E. coli is a Gram-negative, facultative anaerobic bacterium whose life cycle includes extra host environments in addition to human, animal, and plant hosts. The bacterium has the genomic capability of being motile. In this context, the significance of this study is severalfold: (i) the great diversity of migration phenotypes that we observed within our isolate collection supports previous (G. NandaKafle, A. A. Christie, S. Vilain, and V. S. Brözel, Front Microbiol 9:762, 2018, https://doi.org/10.3389/fmicb.2018.00762; Y. Somorin, F. Abram, F. Brennan, and C. O'Byrne, Appl Environ Microbiol 82:4628-4640, 2016, https://doi.org/10.1128/AEM.01175-16) ideas of soil promoting phenotypic heterogeneity, (ii) such heterogeneity may facilitate bacterial growth in the many different soil niches, and (iii) such heterogeneity may enable the bacteria to interact with human, animal, and plant hosts.

About the dark corners in the gene function space of Escherichia coli remaining without illumination by scientific literature

BACKGROUND

Although Escherichia coli (E. coli) is the most studied prokaryote organism in the history of life sciences, many molecular mechanisms and gene functions encoded in its genome remain to be discovered. This work aims at quantifying the illumination of the E. coli gene function space by the scientific literature and how close we are towards the goal of a complete list of E. coli gene functions.

RESULTS

The scientific literature about E. coli protein-coding genes has been mapped onto the genome via the mentioning of names for genomic regions in scientific articles both for the case of the strain K-12 MG1655 as well as for the 95%-threshold softcore genome of 1324 E. coli strains with known complete genome. The article match was quantified with the ratio of a given gene name's occurrence to the mentioning of any gene names in the paper. The various genome regions have an extremely uneven literature coverage. A group of elite genes with ≥ 100 full publication equivalents (FPEs, FPE = 1 is an idealized publication devoted to just a single gene) attracts the lion share of the papers. For K-12, ~ 65% of the literature covers just 342 elite genes; for the softcore genome, ~ 68% of the FPEs is about only 342 elite gene families (GFs). We also find that most genes/GFs have at least one mentioning in a dedicated scientific article (with the exception of at least 137 protein-coding transcripts for K-12 and 26 GFs from the softcore genome). Whereas the literature growth rates were highest for uncharacterized or understudied genes until 2005-2010 compared with other groups of genes, they became negative thereafter. At the same time, literature for anyhow well-studied genes started to grow explosively with threshold T10 (≥ 10 FPEs). Typically, a body of ~ 20 actual articles generated over ~ 15 years of research effort was necessary to reach T10. Lineage-specific co-occurrence analysis of genes belonging to the accessory genome of E. coli together with genomic co-localization and sequence-analytic exploration hints previously completely uncharacterized genes yahV and yddL being associated with osmotic stress response/motility mechanisms.

CONCLUSION

If the numbers of scientific articles about uncharacterized and understudied genes remain at least at present levels, full gene function lists for the strain K-12 MG1655 and the E. coli softcore genome are in reach within the next 25-30 years. Once the literature body for a gene crosses 10 FPEs, most of the critical fundamental research risk appears overcome and steady incremental research becomes possible.

Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Pseudomonas aeruginosa is an opportunistic pathogen that causes difficult-to-treat infections. Two well-studied divergent P. aeruginosa strain types, PAO1 and PA14, have significant genomic heterogeneity, including diverse accessory genes present in only some strains. Genome content comparisons find core genes that are conserved across both PAO1 and PA14 strains and accessory genes that are present in only a subset of PAO1 and PA14 strains. Here, we use recently assembled transcriptome compendia of publicly available P. aeruginosa RNA sequencing (RNA-seq) samples to create two smaller compendia consisting of only strain PAO1 or strain PA14 samples with each aligned to their cognate reference genome. We confirmed strain annotations and identified other samples for inclusion by assessing each sample's median expression of PAO1-only or PA14-only accessory genes. We then compared the patterns of core gene expression in each strain. To do so, we developed a method by which we analyzed genes in terms of which genes showed similar expression patterns across strain types. We found that some core genes had consistent correlated expression patterns across both compendia, while others were less stable in an interstrain comparison. For each accessory gene, we also determined core genes with correlated expression patterns. We found that stable core genes had fewer coexpressed neighbors that were accessory genes. Overall, this approach for analyzing expression patterns across strain types can be extended to other groups of genes, like phage genes, or applied for analyzing patterns beyond groups of strains, such as samples with different traits, to reveal a deeper understanding of regulation. IMPORTANCE Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. Understanding how these different strain-level traits manifest is important for identifying targets that regulate different traits of interest. With the availability of thousands of PAO1 and PA14 samples, we created two strain-specific RNA-seq compendia where each one contains hundreds of samples from PAO1 or PA14 strains and used them to compare the expression patterns of core genes that are conserved in both strain types and to determine which core genes have expression patterns that are similar to those of accessory genes that are unique to one strain or the other using an approach that we developed. We found a subset of core genes with different transcriptional patterns across PAO1 and PA14 strains and identified those core genes with expression patterns similar to those of strain-specific accessory genes.

Comparative Genome Analysis of Enterococcus cecorum Reveals Intercontinental Spread of a Lineage of Clinical Poultry Isolates

Enterococcus cecorum is an emerging pathogen responsible for osteomyelitis, spondylitis, and femoral head necrosis causing animal suffering and mortality and requiring antimicrobial use in poultry. Paradoxically, E. cecorum is a common inhabitant of the intestinal microbiota of adult chickens. Despite evidence suggesting the existence of clones with pathogenic potential, the genetic and phenotypic relatedness of disease-associated isolates remains little investigated. Here, we sequenced and analyzed the genomes and characterized the phenotypes of more than 100 isolates, the majority of which were collected over the last 10 years from 16 French broiler farms. Comparative genomics, genome-wide association studies, and the measured susceptibility to serum, biofilm-forming capacity, and adhesion to chicken type II collagen were used to identify features associated with clinical isolates. We found that none of the tested phenotypes could discriminate the origin of the isolates or the phylogenetic group. Instead, we found that most clinical isolates are grouped phylogenetically, and our analyses selected six genes that discriminate 94% of isolates associated with disease from those that are not. Analysis of the resistome and the mobilome revealed that multidrug-resistant clones of E. cecorum cluster into a few clades and that integrative conjugative elements and genomic islands are the main carriers of antimicrobial resistance. This comprehensive genomic analysis shows that disease-associated clones of E. cecorum belong mainly to one phylogenetic clade. IMPORTANCE Enterococcus cecorum is an important pathogen of poultry worldwide. It causes a number of locomotor disorders and septicemia, particularly in fast-growing broilers. Animal suffering, antimicrobial use, and associated economic losses require a better understanding of disease-associated E. cecorum isolates. To address this need, we performed whole-genome sequencing and analysis of a large collection of isolates responsible for outbreaks in France. By providing the first data set on the genetic diversity and resistome of E. cecorum strains circulating in France, we pinpoint an epidemic lineage that is probably also circulating elsewhere that should be targeted preferentially by preventive strategies in order to reduce the burden of E. cecorum-related diseases.

Identification of Subunits for Novel Universal Vaccines against Three Predominant Serogroups and the Emerging O145 among Avian Pathogenic Escherichia coli by Pan-RV Pipeline

Avian pathogenic Escherichia coli, a causative agent of avian colibacillosis, has been causing serious economic losses in the poultry industry. The increase in multidrug-resistant isolates and the complexity of the serotypes of this pathogen, especially the recently reported emergence of a newly predominant serogroup of O145, make the control of this disease difficult. To address this challenge, a high-throughput screening approach, called Pan-RV (Reverse vaccinology based on pangenome analysis), is proposed to search for universal protective antigens against the three traditional serogroups and the newly emerged O145. Using this approach, a total of 61 proteins regarded as probable antigens against the four important serogroups were screened from the core genome of 127 Avian pathogenic Escherichia coli (APEC) genomes, and six were verified by Western blots using antisera. Overall, our research will provide a foundation for the development of an APEC subunit vaccine against avian colibacillosis. Given the exponential growth of whole-genome sequencing (WGS) data, our Pan-RV pipeline will make screening of bacterial vaccine candidates inexpensive, rapid, and efficient. IMPORTANCE With the emergence of drug resistance and the newly predominant serogroup O145, the control of Avian pathogenic Escherichia coli is facing a serious challenge; an efficient immunological method is urgently needed. Here, for the first time, we propose a high-throughput screening approach to search for universal protective antigens against the three traditional serogroups and the newly emerged O145. Importantly, using this approach, a total of 61 proteins regarded as probable antigens against the four important serogroups were screened, and three were shown to be immunoreactive with all antisera (covering the four serogroups), thereby providing a foundation for the development of APEC subunit vaccines against avian colibacillosis. Further, our Pan-RV pipeline will provide immunological control strategies for pathogens with complex and variable genetic backgrounds such as Escherichia coli and will make screening of bacterial vaccine candidates more inexpensive, rapid, and efficient.

A high-throughput sequencing approach identifies immunotherapeutic targets for bacterial meningitis in neonates

BACKGROUND

Worldwide, Escherichia coli is the leading cause of neonatal Gram-negative bacterial meningitis, but full understanding of the pathogenesis of this disease is not yet achieved. Moreover, to date, no vaccine is available against bacterial neonatal meningitis.

METHODS

Here, we used Transposon Sequencing of saturated banks of mutants (TnSeq) to evaluate E. coli K1 genetic fitness in murine neonatal meningitis. We identified E. coli K1 genes encoding for factors important for systemic dissemination and brain infection, and focused on products with a likely outer-membrane or extra-cellular localization, as these are potential vaccine candidates. We used in vitro and in vivo models to study the efficacy of active and passive immunization.

RESULTS

We selected for further study the conserved surface polysaccharide Poly-β-(1-6)-N-Acetyl Glucosamine (PNAG), as a strong candidate for vaccine development. We found that PNAG was a virulence factor in our animal model. We showed that both passive and active immunization successfully prevented and/or treated meningitis caused by E. coli K1 in neonatal mice. We found an excellent opsonophagocytic killing activity of the antibodies to PNAG and in vitro these antibodies were also able to decrease binding, invasion and crossing of E. coli K1 through two blood brain barrier cell lines. Finally, to reinforce the potential of PNAG as a vaccine candidate in bacterial neonatal meningitis, we demonstrated that Group B Streptococcus, the main cause of neonatal meningitis in developed countries, also produced PNAG and that antibodies to PNAG could protect in vitro and in vivo against this major neonatal pathogen.

INTERPRETATION

Altogether, these results indicate the utility of a high-throughput DNA sequencing method to identify potential immunotherapy targets for a pathogen, including in this study a potential broad-spectrum target for prevention of neonatal bacterial infections.

FUNDINGS

ANR Seq-N-Vaq, Charles Hood Foundation, Hearst Foundation, and Groupe Pasteur Mutualité.

Identification of highly variable sequence fragments in unmapped reads for rapid bacterial genotyping

BACKGROUND

Bacterial genotyping is a crucial process in outbreak investigation and epidemiological studies. Several typing methods such as pulsed-field gel electrophoresis, multilocus sequence typing (MLST) and whole genome sequencing are currently used in routine clinical practice. However, these methods are costly, time-consuming and have high computational demands. An alternative to these methods is mini-MLST, a quick, cost-effective and robust method based on high-resolution melting analysis. Nevertheless, no standardized approach to identify markers suitable for mini-MLST exists. Here, we present a pipeline for variable fragment detection in unmapped reads based on a modified hybrid assembly approach using data from one sequencing platform.

RESULTS

In routine assembly against the reference sequence, high variable reads are not aligned and remain unmapped. If de novo assembly of them is performed, variable genomic regions can be located in created scaffolds. Based on the variability rates calculation, it is possible to find a highly variable region with the same discriminatory power as seven housekeeping gene fragments used in MLST. In the work presented here, we show the capability of identifying one variable fragment in de novo assembled scaffolds of 21 Escherichia coli genomes and three variable regions in scaffolds of 31 Klebsiella pneumoniae genomes. For each identified fragment, the melting temperatures are calculated based on the nearest neighbor method to verify the mini-MLST's discriminatory power.

CONCLUSIONS

A pipeline for a modified hybrid assembly approach consisting of reference-based mapping and de novo assembly of unmapped reads is presented. This approach can be employed for the identification of highly variable genomic fragments in unmapped reads. The identified variable regions can then be used in efficient laboratory methods for bacterial typing such as mini-MLST with high discriminatory power, fully replacing expensive methods such as MLST. The results can and will be delivered in a shorter time, which allows immediate and fast infection monitoring in clinical practice.

Evidence for Widespread Class II Microcins in Enterobacterales Genomes

Microcins are a class of antimicrobial peptides produced by certain Gram-negative bacterial species to kill or inhibit the growth of competing bacteria. Only 10 unique, experimentally validated class II microcins have been identified, and the majority of these come from Escherichia coli. Although the current representation of microcins is sparse, they exhibit a diverse array of molecular functionalities, uptake mechanisms, and target specificities. This broad diversity from such a small representation suggests that microcins may have untapped potential for bioprospecting peptide antibiotics from genomic data sets. We used a systematic bioinformatics approach to search for verified and novel class II microcins in E. coli and other species within its family, Enterobacteriaceae. Nearly one-quarter of the E. coli genome assemblies contained one or more microcins, where the prevalence of hits to specific microcins varied by isolate phylogroup. E. coli isolates from human extraintestinal and poultry meat sources were enriched for microcins, while those from freshwater were depleted. Putative microcins were found in various abundances across all five distinct phylogenetic lineages of Enterobacteriaceae, with a particularly high prevalence in the "Klebsiella" clade. Representative genome assemblies from species across the Enterobacterales order, as well as a few outgroup species, also contained putative microcin sequences. This study suggests that microcins have a complicated evolutionary history, spanning far beyond our limited knowledge of the currently validated microcins. Efforts to functionally characterize these newly identified microcins have great potential to open a new field of peptide antibiotics and microbiome modulators and elucidate the ways in which bacteria compete with each other. IMPORTANCE Class II microcins are small bacteriocins produced by strains of Gram-negative bacteria in the Enterobacteriaceae. They are generally understood to play a role in interbacterial competition, although direct evidence of this is limited, and they could prove informative in developing new peptide antibiotics. However, few examples of verified class II microcins exist, and novel microcins are difficult to identify due to their sequence diversity, making it complicated to study them as a group. Here, we overcome this limitation by developing a bioinformatics pipeline to detect microcins in silico. Using this pipeline, we demonstrate that both verified and novel class II microcins are widespread within and outside the Enterobacteriaceae, which has not been systematically shown previously. The observed prevalence of class II microcins suggests that they are ecologically important, and the elucidation of novel microcins provides a resource that can be used to expand our knowledge of the structure and function of microcins as antibacterials.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

When bacteria are phage playgrounds: interactions between viruses, cells, and mobile genetic elements

Studies of viral adaptation have focused on the selective pressures imposed by hosts. However, there is increasing evidence that interactions between viruses, cells, and other mobile genetic elements are determinant to the success of infections. These interactions are often associated with antagonism and competition, but sometimes involve cooperation or parasitism. We describe two key types of interactions - defense systems and genetic regulation - that allow the partners of the interaction to destroy or control the others. These interactions evolve rapidly by genetic exchanges, including among competing partners. They are sometimes followed by functional diversification. Gene exchanges also facilitate the emergence of cross-talk between elements in the same bacterium. In the end, these processes produce multilayered networks of interactions that shape the outcome of viral infections.

Origins of transfer establish networks of functional dependencies for plasmid transfer by conjugation

Plasmids can be transferred between cells by conjugation, thereby driving bacterial evolution by horizontal gene transfer. Yet, we ignore the molecular mechanisms of transfer for many plasmids because they lack all protein-coding genes required for conjugation. We solved this conundrum by identifying hundreds of plasmids and chromosomes with conjugative origins of transfer in Escherichia coli and Staphylococcus aureus. These plasmids (pOriT) hijack the relaxases of conjugative or mobilizable elements, but not both. The functional dependencies between pOriT and other plasmids explain their co-occurrence: pOriT are abundant in cells with many plasmids, whereas conjugative plasmids are the most common in the others. We systematically characterized plasmid mobility in relation to conjugation and alternative mechanisms of transfer and can now propose a putative mechanism of transfer for ∼90% of them. In most cases, plasmid mobility seems to involve conjugation. Interestingly, the mechanisms of mobility are important determinants of plasmid-encoded accessory traits, since pOriTs have the highest densities of antimicrobial resistance genes, whereas plasmids lacking putative mechanisms of transfer have the lowest. We illuminate the evolutionary relationships between plasmids and suggest that many pOriT may have arisen by gene deletions in other types of plasmids. These results suggest that most plasmids can be transferred by conjugation.

Role of adherent invasive Escherichia coli in pathogenesis of inflammatory bowel disease

Gut microbiota imbalances play an important role in inflammatory bowel disease (IBD), but no single pathogenic microorganism critical to IBD that is specific to the IBD terminal ileum mucosa or can invade intestinal epithelial cells has been found. Invasive Escherichia coli (E. coli) adhesion to macrophages is considered to be closely related to the pathogenesis of inflammatory bowel disease. Further study of the specific biological characteristics of adherent invasive E. coli (AIEC) may contribute to a further understanding of IBD pathogenesis. This review explores the relationship between AIEC and the intestinal immune system, discusses the prevalence and relevance of AIEC in Crohn's disease and ulcerative colitis patients, and describes the relationship between AIEC and the disease site, activity, and postoperative recurrence. Finally, we highlight potential therapeutic strategies to attenuate AIEC colonization in the intestinal mucosa, including the use of phage therapy, antibiotics, and anti-adhesion molecules. These strategies may open up new avenues for the prevention and treatment of IBD in the future.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

Lactic Acid Resistance and Population Structure of Escherichia coli from Meat Processing Environment

To explore the effect of beef processing on Escherichia coli populations in relation to lactic acid resistance, this study investigated the links among acid response, phylogenetic structure, genome diversity, and genotypes associated with acid resistance of meat plant E. coli. Generic E. coli isolates (n = 700) were from carcasses, fabrication equipment, and beef products. Acid treatment was carried out in Luria-Bertani broth containing 5.5% lactic acid (pH 2.9). Log reductions of E. coli ranged from <0.5 to >5 log CFU/mL (median: 1.37 log). No difference in lactic acid resistance was observed between E. coli populations recovered before and after a processing step or antimicrobial interventions. E. coli from the preintervention carcasses were slightly more resistant than E. coli isolated from equipment, differing by <0.5 log unit. Acid-resistant E. coli (log reduction <1, n = 45) had a higher prevalence of genes related to energy metabolism (ydj, xap, ato) and oxidative stress (fec, ymjC) than the less resistant E. coli (log reduction >1, n = 133). The ydj and ato operons were abundant in E. coli from preintervention carcasses. In contrast, fec genes were abundant in E. coli from equipment surfaces. The preintervention E. coli contained phylogroups A and B1 in relatively equal proportions. Phylogroup B1 predominated (95%) in the population from equipment. Of note, E. coli collected after sanitation shared either the antigens of O8 or H21. Additionally, genome diversity decreased after chilling and equipment sanitation. Overall, beef processing did not select for E. coli resistant to lactic acid but shaped the population structure. IMPORTANCE Antimicrobial interventions have significantly reduced the microbial loads on carcasses/meat products; however, the wide use of chemical and physical biocides has raised concerns over their potential for selecting resistant populations in the beef processing environment. Phenotyping of acid resistance and whole-genome analysis described in this study demonstrated beef processing practices led to differences in acid resistance, genotype, and population structure between carcass- and equipment-associated E. coli but did not select for the acid-resistant population. Results indicate that genes coding for the metabolism of long-chain sugar acids (ydj) and short-chain fatty acids (ato) were more prevalent in carcass-associated than equipment-associated E. coli. These results suggest E. coli from carcasses and equipment surfaces have been exposed to different selective pressures. The findings improve our understanding of the microbial ecology of E. coli in food processing environments and in general.

Genome-scale metabolic network reconstructions of diverse Escherichia strains reveal strain-specific adaptations

Bottom-up approaches to systems biology rely on constructing a mechanistic basis for the biochemical and genetic processes that underlie cellular functions. Genome-scale network reconstructions of metabolism are built from all known metabolic reactions and metabolic genes in a target organism. A network reconstruction can be converted into a mathematical format and thus lend itself to mathematical analysis. Genome-scale models (GEMs) of metabolism enable a systems approach to characterize the pan and core metabolic capabilities of the Escherichia genus. In this work, GEMs were constructed for 222 representative strains of Escherichia across HC1100 levels spanning the known Escherichia phylogeny. The models were used to study Escherichia metabolic diversity and speciation on a large scale. The results show that unique strain-specific metabolic capabilities correspond to different species and nutrient niches. This work is a first step towards a curated reconstruction of pan-Escherichia metabolism.

This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

Selfish, promiscuous and sometimes useful: how mobile genetic elements drive horizontal gene transfer in microbial populations

Horizontal gene transfer (HGT) drives microbial adaptation but is often under the control of mobile genetic elements (MGEs) whose interests are not necessarily aligned with those of their hosts. In general, transfer is costly to the donor cell while potentially beneficial to the recipients. The diversity and plasticity of cell–MGEs interactions, and those among MGEs, result in complex evolutionary processes where the source, or even the existence of selection for maintaining a function in the genome, is often unclear. For example, MGE-driven HGT depends on cell envelope structures and defense systems, but many of these are transferred by MGEs themselves. MGEs can spur periods of intense gene transfer by increasing their own rates of horizontal transmission upon communicating, eavesdropping, or sensing the environment and the host physiology. This may result in high-frequency transfer of host genes unrelated to the MGE. Here, we review how MGEs drive HGT and how their transfer mechanisms, selective pressures and genomic traits affect gene flow, and therefore adaptation, in microbial populations. The encoding of many adaptive niche-defining microbial traits in MGEs means that intragenomic conflicts and alliances between cells and their MGEs are key to microbial functional diversification.

This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Implications of Escherichia coli community diversity in free-ranging Australian pinniped pups

Escherichia coli is a widely studied bacterium, commonly used as an indicator of faecal contamination. Investigations into the structure and diversity of E. coli in free-ranging wildlife species has been limited. The objective of this study was to characterise intra-individual and inter-species E. coli phylotype and B2 sub-type diversity in free-ranging Australian pinniped pups, to determine whether a single E. coli colony is representative of the phylotype and B2 sub-type diversity in these hosts. Faecal samples were collected from free-ranging Australian fur seal (Arctocephalus pusillus doriferus), Australian sea lion (Neophoca cinerea) and long-nosed fur seal (Arctocephalus forsteri) pups from three breeding colonies between 2018 and 2021. Faecal swabs from thirty randomly selected pups (n = 10 from each species) were cultured and ten E. coli colonies were selected from each culture based on morphology and separation between colonies on agar plates. Molecular screening techniques were utilised to assign isolates to phylotypes and B2 sub-types. There was no significant difference (p > 0.05) in either intra-individual or inter-species E. coli phylotype and B2 sub-type diversity. The B2 phylotype was the most dominant, with 78% of isolates (n = 234) assigned to this phylotype. Host factors (species, weight [kg] and standard length [cm]) did not significantly affect phylotype diversity. The absence of intra-individual and inter-species differences in E. coli diversity at a phylotype level suggests that a single E. coli colony could be used as an indicator of overall diversity of E. coli at a phylotype level in A. p. doriferus, N. cinerea and A. forsteri pups. These findings can be used to simplify and improve the efficiency of sampling protocols for ongoing monitoring of human-associated E. coli phylotypes in free-ranging pinniped populations.

Spatiotemporal dynamics of Escherichia coli presence and magnitude across a national groundwater monitoring network, Republic of Ireland, 2011-2020

Groundwater is a vital drinking water resource and its protection from microbiological contamination is paramount to safeguard public health. The Republic of Ireland (RoI) is characterised by the highest incidence of verocytotoxigenic Escherichia coli (VTEC) enteritis in the European Union (EU), linked to high reliance on unregulated groundwater sources (~16% of the population). Yet, the spatio-temporal factors influencing the frequency and magnitude of microbial contamination remain largely unknown, with past studies typically constrained to spatio-temporally 'limited' sampling campaigns. Accordingly, the current investigation sought to analyse an extensive spatially distributed time-series (2011-2020) of groundwater monitoring data in the RoI. The dataset, compiled by the Environmental Protection Agency (EPA), showed 'high' contamination rates, with 66.7% (88/132) of supplies testing positive for E. coli, and 29.5% (39/132) exceeding concentrations of 10MPN/100 ml (i.e. gross contamination) at least once during the 10-year monitoring period. Seasonal decomposition analyses indicate that E. coli detection rates peak during late autumn/early winter, coinciding with increases in annual rainfall, while gross contamination peaks in spring (May) and late-summer (August), likely reflecting seasonal shifts in agricultural practices. Mixed effects logistic regression modelling indicates that monitoring sources located in karst limestone are statistically associated with E. coli presence (OR = 2.76, p = 0.03) and gross contamination (OR = 2.54, p = 0.037) when compared to poorly productive aquifers (i.e., transmissivity below 10m²/d). Moreover, 5-day and 30-day antecedent rainfall increased the likelihood of E. coli contamination (OR = 1.027, p < 0.001 and OR = 1.005, p = 0.016, respectively), with the former also being associated with gross contamination (OR = 1.042, p < 0.001). As such, it is inferred that preferential flow and direct ingress of surface runoff are the most likely ingress mechanisms associated with E. coli groundwater supply contamination. The results presented are expected to inform policy change around groundwater source protection and provide insight for the development of groundwater monitoring programmes in geologically heterogeneous regions.

Comparative Genomic Analysis of Antimicrobial-Resistant Escherichia coli from South American Camelids in Central Germany

South American camelids (SAC) are increasingly kept in Europe in close contact with humans and other livestock species and can potentially contribute to transmission chains of epizootic, zoonotic and antimicrobial-resistant (AMR) agents from and to livestock and humans. Consequently, SAC were included as livestock species in the new European Animal Health Law. However, the knowledge on bacteria exhibiting AMR in SAC is too scarce to draft appropriate monitoring and preventive programs. During a survey of SAC holdings in central Germany, 39 Escherichia coli strains were isolated from composite fecal samples by selecting for cephalosporin or fluoroquinolone resistance and were here subjected to whole-genome sequencing. The data were bioinformatically analyzed for strain phylogeny, detection of pathovars, AMR genes and plasmids. Most (33/39) strains belonged to phylogroups A and B1. Still, the isolates were highly diverse, as evidenced by 28 multi-locus sequence types. More than half of the isolates (23/39) were genotypically classified as multidrug resistant. Genes mediating resistance to trimethoprim/sulfonamides (22/39), aminoglycosides (20/39) and tetracyclines (18/39) were frequent. The most common extended-spectrum-β-lactamase gene was bla_CTX-M-1 (16/39). One strain was classified as enteropathogenic E. coli. The positive results indicate the need to include AMR bacteria in yet-to-be-established animal disease surveillance protocols for SAC.

The Population Genomics of Increased Virulence and Antibiotic Resistance in Human Commensal Escherichia coli over 30 Years in France

Escherichia coli is a commensal species of the lower intestine but is also a major pathogen causing intestinal and extraintestinal infections that is increasingly prevalent and resistant to antibiotics. Most studies on genomic evolution of E. coli used isolates from infections. Here, instead, we whole-genome sequenced a collection of 403 commensal E. coli isolates from fecal samples of healthy adult volunteers in France (1980 to 2010). These isolates were distributed mainly in phylogroups A and B2 (30% each) and belonged to 152 sequence types (STs), the five most frequent being ST10 (phylogroup A; 16.3%), ST73 and ST95 (phylogroup B2; 6.3 and 5.0%, respectively), ST69 (phylogroup D; 4.2%), and ST59 (phylogroup F; 3.9%), and 224 O:H serotypes. ST and serotype diversity increased over time. The O1, O2, O6, and O25 groups used in bioconjugate O-antigen vaccine against extraintestinal infections were found in 23% of the strains of our collection. The increase in frequency of virulence-associated genes and antibiotic resistance was driven by two evolutionary mechanisms. Evolution of virulence gene frequency was driven by both clonal expansion of STs with more virulence genes ("ST-driven") and increases in gene frequency within STs independent of changes in ST frequencies ("gene-driven"). In contrast, the evolution of resistance was dominated by increases in frequency within STs ("gene-driven"). This study provides a unique picture of the phylogenomic evolution of E. coli in its human commensal habitat over 30 years and will have implications for the development of preventive strategies. IMPORTANCE Escherichia coli is an opportunistic pathogen with the greatest burden of antibiotic resistance, one of the main causes of bacterial infections and an increasing concern in an aging population. Deciphering the evolutionary dynamics of virulence and antibiotic resistance in commensal E. coli is important to understand adaptation and anticipate future changes. The gut of vertebrates is the primary habitat of E. coli and probably where selection for virulence and resistance takes place. Unfortunately, most whole-genome-sequenced strains are isolated from pathogenic conditions. Here, we whole-genome sequenced 403 E. coli commensals isolated from healthy French subjects over a 30-year period. Virulence genes increased in frequency by both clonal expansion of clones carrying them and increases in frequency within clones, whereas resistance genes increased by within-clone increased frequency. Prospective studies of E. coli commensals should be performed worldwide to have a broader picture of evolution and adaptation of this species.

Pre-existing chromosomal polymorphisms in pathogenic E. coli potentiate the evolution of resistance to a last-resort antibiotic

Bacterial pathogens show high levels of chromosomal genetic diversity, but the influence of this diversity on the evolution of antibiotic resistance by plasmid acquisition remains unclear. Here, we address this problem in the context of colistin, a 'last line of defence' antibiotic. Using experimental evolution, we show that a plasmid carrying the MCR-1 colistin resistance gene dramatically increases the ability of Escherichia coli to evolve high-level colistin resistance by acquiring mutations in lpxC, an essential chromosomal gene involved in lipopolysaccharide biosynthesis. Crucially, lpxC mutations increase colistin resistance in the presence of the MCR-1 gene, but decrease the resistance of wild-type cells, revealing positive sign epistasis for antibiotic resistance between the chromosomal mutations and a mobile resistance gene. Analysis of public genomic datasets shows that lpxC polymorphisms are common in pathogenic E. coli, including those carrying MCR-1, highlighting the clinical relevance of this interaction. Importantly, lpxC diversity is high in pathogenic E. coli from regions with no history of MCR-1 acquisition, suggesting that pre-existing lpxC polymorphisms potentiated the evolution of high-level colistin resistance by MCR-1 acquisition. More broadly, these findings highlight the importance of standing genetic variation and plasmid/chromosomal interactions in the evolutionary dynamics of antibiotic resistance.

Molecular and metabolic characteristics of wastewater associated Escherichia coli strains

We previously characterized the genetic diversity of Escherichia coli strains isolated from septic tanks in the Canberra region, Australia. In this study, we used repetitive element palindromic (REP) PCR fingerprinting to identify dominant REP-types belonging to phylogroups A and B1 strains across septic tanks. Subsequently, 76 E. coli strains were selected for whole-genome sequencing and phenotype microarrays. Comparative genome analysis was performed to compare septic tank E. coli genomes with a collection of 433 E. coli isolates from different hosts and freshwater. Clonal complexes (CCs) 10 (n = 15) and 399 (n = 10) along with sequence type (ST) 401 (n = 9) were the common lineages in septic tanks. CC10 strains have been detected from animal hosts and freshwater, whereas CC399 and ST401 strains appeared to be associated with septic tanks as they were uncommon in isolates from other sources. Comparative genome analysis revealed that CC399 and ST401 were genetically distinct from other isolates and carried an abundance of niche-specific traits involved in environmental adaptation. These strains also showed distinct metabolic characteristics, such as the ability to utilize pectin, which may provide a fitness advantage under nutrient-limited conditions. The results of this study characterized the adaptive mechanisms allowing E. coli to persist in wastewater.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

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

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 16^th 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.

Ancient DNA from an Italian mummy’s gallstone provides insight into opportunistic E. coli infection.