A comprehensive and high-quality collection of Escherichia coli genomes and their genes

Whole-genome sequencing for antimicrobial surveillance: species-specific quality thresholds and data evaluation from the network of the European Union Reference Laboratory for Antimicrobial Resistance genomic proficiency tests of 2021 and 2022

As antimicrobial resistance (AMR) surveillance shifts to genomics, ensuring the quality of whole-genome sequencing (WGS) data produced across laboratories is critical. Participation in genomic proficiency tests (GPTs) not only increases individual laboratories' WGS capacity but also provides a unique opportunity to improve species-specific thresholds for WGS quality control (QC) by repeated resequencing of distinct isolates. Here, we present the results of the EU Reference Laboratory for Antimicrobial Resistance (EURL-AR) network GPTs of 2021 and 2022, which included 25 EU national reference laboratories (NLRs). A total of 392 genomes from 12 AMR-bacteria were evaluated based on WGS QC metrics. Two percent (n = 9) of the data were excluded, due to contamination, and 11% (n = 41) of the remaining genomes were identified as outliers in at least one QC metric and excluded from computation of the adjusted QC thresholds (AQT). Two QC metric correlation groups were identified through linear regression. Eight percent (n = 28) of the submitted genomes, from 11 laboratories, failed one or more of the AQTs. However, only three laboratories (12%) were identified as underperformers, failing across AQTs for uncorrelated QC metrics in at least two genomes. Finally, new species-specific thresholds for "N50" and "number of contigs > 200 bp" are presented for guidance in routine laboratory QC. The continued participation of NRLs in GPTs will reveal WGS workflow flaws and improve AMR surveillance data. GPT data will continue to contribute to the development of reliable species-specific thresholds for routine WGS QC, standardizing sequencing data QC and ensure inter- and intranational laboratory comparability.IMPORTANCEIllumina next-generation sequencing is an integral part of antimicrobial resistance (AMR) surveillance and the most widely used whole-genome sequencing (WGS) platform. The high-throughput, relative low-cost, high discriminatory power, and rapid turnaround time of WGS compared to classical biochemical methods means the technology will likely remain a fundamental tool in AMR surveillance and public health. In this study, we present the current level of WGS capacity among national reference laboratories in the EU Reference Laboratory for AMR network, summarizing applied methodology and statistically evaluating the quality of the obtained sequence data. These findings provide the basis for setting new and revised thresholds for quality metrics used in routine WGS, which have previously been arbitrarily defined. In addition, underperforming participants are identified and encouraged to evaluate their workflows to produce reliable results.

Genotypic and phenotypic characterisation of asymptomatic bacteriuria (ABU) isolates displaying bacterial interference against multi-drug resistant uropathogenic E. Coli

Escherichia coli can colonise the urogenital tract of individuals without causing symptoms of infection, in a condition referred to as asymptomatic bacteriuria (ABU). ABU isolates can protect the host against symptomatic urinary tract infections (UTIs) by bacterial interference against uropathogenic E. coli (UPEC). The aim of this study was to investigate the genotypic and phenotypic characteristics of five ABU isolates from midstream urine samples of adults. Comparative genomic and phenotypic analysis was conducted including an antibiotic resistance profile, pangenome analysis, and a putative virulence profile. Based on the genome analysis, the isolates consisted of one from phylogroup A, three from phylogroup B2, and one from phylogroup D. Two of the isolates, PUTS 58 and SK-106-1, were noted for their lack of antibiotic resistance and virulence genes compared to the prototypic ABU strain E. coli 83,972. This study provides insights into the genotypic and phenotypic profiles of uncharacterised ABU isolates, and how relevant fitness and virulence traits can impact their potential suitability for therapeutic bacterial interference.

Whole genome sequence analysis reveals high genomic diversity and potential host-driven adaptations among multidrug-resistant Escherichia coli from pre-weaned dairy calves

Food-producing animals such as dairy cattle are potential reservoirs of antimicrobial resistance (AMR), with multidrug-resistant (MDR) organisms such as Escherichia coli observed in higher frequency in young calves compared to older cattle. In this study, we characterized the genomes of enteric MDR E. coli from pre-weaned dairy calves with and without diarrhea and evaluated the influence of host-level factors on genomic composition. Whole genome sequence comparative analysis of E. coli (n = 43) revealed substantial genomic diversity that primarily clustered by sequence type and was minimally driven by calf diarrheal disease status (healthy, diarrheic, or recovered), antimicrobial exposure, and dietary zinc supplementation. Diverse AMR genes (ARGs)-including extended-spectrum beta-lactamase genes and quinolone resistance determinants-were identified (n = 40), with unique sets of ARGs co-occurring in gene clusters with large AMR plasmids IncA/C2 and IncFIB(AP001918). Zinc supplementation was not significantly associated with the selection of individual ARGs in E. coli, however analysis of ARG and metal resistance gene pairs identified positive associations between certain aminoglycoside, beta-lactam, sulfonamide, and trimethoprim ARGs with acid, tellurium and mercury resistance genes. Although E. coli in this study lacked the typical virulence factors of diarrheagenic strains, virulence genes overlapping with those in major pathotypes were identified. Among the 103 virulence genes detected, the highest abundance and diversity of genes corresponded to iron acquisition (siderophores and heme uptake). Our findings indicate that the host-level factors evaluated in this study were not key drivers of genomic variability, but that certain accessory genes in enteric MDR E. coli may be enriched. Collectively, this work provides insight into the genomic diversity and host-microbe interface of MDR E. coli from pre-weaned dairy calves.

The inner membrane protein YhiM links copper and CpxAR envelope stress responses in uropathogenic E. coli

Urinary tract infection (UTI) is a ubiquitous infectious condition, and uropathogenic Escherichia coli (UPEC) is the predominant causative agent of UTI. Copper (Cu) is implicated in innate immunity, including against UPEC. Cu is a trace element utilized as a co-factor, but excess Cu is toxic due to mismetalation of non-cognate proteins. E. coli precisely regulates Cu homeostasis via efflux systems. However, Cu import mechanisms into the bacterial cell are not clear. We hypothesized that Cu import defective mutants would exhibit increased resistance to Cu. This hypothesis was tested in a forward genetic screen with transposon (Tn5) insertion mutants in UPEC strain CFT073, and we identified 32 unique Cu-resistant mutants. Transposon and defined mutants lacking yhiM, which encodes a hypothetical inner membrane protein, were more resistant to Cu than parental strain. Loss of YhiM led to decreased cellular Cu content and increased expression of copA, encoding a Cu efflux pump. The CpxAR envelope stress response system was activated in the ΔyhiM mutant as indicated by increased expression of cpxP. Transcription of yhiM was regulated by CueR and CpxR, and the CpxAR system was essential for increased Cu resistance in the ΔyhiM mutant. Importantly, activation of CpxAR system in the ΔyhiM mutant was independent of NlpE, a known activator of this system. YhiM was required for optimal fitness of UPEC in a mouse model of UTI. Our findings demonstrate that YhiM is a critical mediator of Cu homeostasis and links bacterial adaptation to Cu stress with the CpxAR-dependent envelope stress response in UPEC.IMPORTANCEUPEC is a common bacterial infection. Bacterial pathogens are exposed to host-derived Cu during infection, including UTI. Here, we describe detection of genes involved in Cu homeostasis in UPEC. A UPEC mutant lacking YhiM, a membrane protein, exhibited dramatic increase in resistance to Cu. Our study demonstrates YhiM as a nexus between Cu stress and the CpxAR-dependent envelope stress response system. Importantly, our findings establish NlpE-independent activation of CpxAR system during Cu stress in UPEC. Collectively, YhiM emerges as a critical mediator of Cu homeostasis in UPEC and highlights the interlinked nature of bacterial adaptation to survival during Cu and envelope stress.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

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.

Themisto: a scalable colored k-mer index for sensitive pseudoalignment against hundreds of thousands of bacterial genomes

MOTIVATION

Huge datasets containing whole-genome sequences of bacterial strains are now commonplace and represent a rich and important resource for modern genomic epidemiology and metagenomics. In order to efficiently make use of these datasets, efficient indexing data structures-that are both scalable and provide rapid query throughput-are paramount.

RESULTS

Here, we present Themisto, a scalable colored k-mer index designed for large collections of microbial reference genomes, that works for both short and long read data. Themisto indexes 179 thousand Salmonella enterica genomes in 9 h. The resulting index takes 142 gigabytes. In comparison, the best competing tools Metagraph and Bifrost were only able to index 11 000 genomes in the same time. In pseudoalignment, these other tools were either an order of magnitude slower than Themisto, or used an order of magnitude more memory. Themisto also offers superior pseudoalignment quality, achieving a higher recall than previous methods on Nanopore read sets.

AVAILABILITY AND IMPLEMENTATION

Themisto is available and documented as a C++ package at https://github.com/algbio/themisto available under the GPLv2 license.

Evolutionary and functional history of the Escherichia coli K1 capsule

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

Genealogical inference and more flexible sequence clustering using iterative-PopPUNK

Bacterial genome data are accumulating at an unprecedented speed due to the routine use of sequencing in clinical diagnoses, public health surveillance, and population genetics studies. Genealogical reconstruction is fundamental to many of these uses; however, inferring genealogy from large-scale genome data sets quickly, accurately, and flexibly is still a challenge. Here, we extend an alignment- and annotation-free method, PopPUNK, to increase its flexibility and interpretability across data sets. Our method, iterative-PopPUNK, rapidly produces multiple consistent cluster assignments across a range of sequence identities. By constructing a partially resolved genealogical tree with respect to these clusters, users can select a resolution most appropriate for their needs. We showed the accuracy of clusters at all levels of similarity and genealogical inference of iterative-PopPUNK based on simulated data and obtained phylogenetically concordant results in real data sets from seven bacterial species. Using two example sets of Escherichia/Shigella and Vibrio parahaemolyticus genomes, we show that iterative-PopPUNK can achieve cluster resolutions ranging from phylogroup down to sequence typing (ST). The iterative-PopPUNK algorithm is implemented in the "PopPUNK_iterate" program, available as part of the PopPUNK package.

Genomic analysis of extended-spectrum beta-lactamase (ESBL) producing Escherichia coli colonising adults in Blantyre, Malawi reveals previously undescribed diversity

Escherichia coli is one of the most prevalent Gram-negative species associated with drug resistant infections. Strains that produce extended-spectrum beta-lactamases (ESBLs) or carbapenemases are both particularly problematic and disproportionately impact resource limited healthcare settings where last-line antimicrobials may not be available. A large number of E. coli genomes are now available and have allowed insights into pathogenesis and epidemiology of ESBL E. coli but genomes from sub-Saharan Africa (sSA) are significantly underrepresented. To reduce this gap, we investigated ESBL-producing E. coli colonising adults in Blantyre, Malawi to assess bacterial diversity and AMR determinants and to place these isolates in the context of the wider population structure. We performed short-read whole-genome sequencing of 473 colonising ESBL E. coli isolated from human stool and contextualised the genomes with a previously curated multi-country collection of 10 146 E. coli genomes and sequence type (ST)-specific collections for our three most commonly identified STs. These were the globally successful ST131, ST410 and ST167, and the dominant ESBL genes were bla CTX-M, mirroring global trends. However, 37 % of Malawian isolates did not cluster with any isolates in the curated multicountry collection and phylogenies were consistent with locally spreading monophyletic clades, including within the globally distributed, carbapenemase-associated B4/H24RxC ST410 lineage. A single ST2083 isolate in this collection harboured a carbapenemase gene. Long read sequencing demonstrated the presence of a globally distributed ST410-associated carbapenemase carrying plasmid in this isolate, which was absent from the ST410 strains in our collection. We conclude there is a risk that carbapenem resistance in E. coli could proliferate rapidly in Malawi under increasing selection pressure, and that both ongoing antimicrobial stewardship and genomic surveillance are critical as local carbapenem use increases.

E. coli ST11 (O157:H7) does not encode a functional AcrF efflux pump

Escherichia coli is a facultative anaerobe found in a wide range of environments. Commonly described as the laboratory workhorse, E. coli is one of the best characterized bacterial species to date, however much of our understanding comes from studies involving the laboratory strain E. coli K-12. Resistance-nodulation-division efflux pumps are found in Gram-negative bacteria and can export a diverse range of substrates, including antibiotics. E. coli K-12 has six RND pumps; AcrB, AcrD, AcrF, CusA, MdtBC and MdtF, and it is frequently reported that all E. coli strains possess these six pumps. However, this is not true of E. coli ST11, a lineage of E. coli, which is primarily composed of the highly virulent important human pathogen, E. coli O157:H7. Here we show that acrF is absent from the pangenome of ST11 and that this lineage of E. coli has a highly conserved insertion within the acrF gene, which when translated encodes 13 amino acids and two stop codons. This insertion was found to be present in 97.59 % of 1787 ST11 genome assemblies. Non-function of AcrF in ST11 was confirmed in the laboratory as complementation with acrF from ST11 was unable to restore AcrF function in E. coli K-12 substr. MG1655 ΔacrB ΔacrF. This shows that the complement of RND efflux pumps present in laboratory bacterial strains may not reflect the situation in virulent strains of bacterial pathogens.

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.

Distinct ecological fitness factors coordinated by a conserved Escherichia coli regulator during systemic bloodstream infection

The ability of bacterial pathogens to adapt to host niches is driven by the carriage and regulation of genes that benefit pathogenic lifestyles. Genes that encode virulence or fitness-enhancing factors must be regulated in response to changing host environments to allow rapid response to challenges presented by the host. Furthermore, this process can be controlled by preexisting transcription factors (TFs) that acquire new roles in tailoring regulatory networks, specifically in pathogens. However, the mechanisms underlying this process are poorly understood. The highly conserved Escherichia coli TF YhaJ exhibits distinct genome-binding dynamics and transcriptome control in pathotypes that occupy different host niches, such as uropathogenic E. coli (UPEC). Here, we report that this important regulator is required for UPEC systemic survival during murine bloodstream infection (BSI). This advantage is gained through the coordinated regulation of a small regulon comprised of both virulence and metabolic genes. YhaJ coordinates activation of both Type 1 and F1C fimbriae, as well as biosynthesis of the amino acid tryptophan, by both direct and indirect mechanisms. Deletion of yhaJ or the individual genes under its control leads to attenuated survival during BSI. Furthermore, all three systems are up-regulated in response to signals derived from serum or systemic host tissue, but not urine, suggesting a niche-specific regulatory trigger that enhances UPEC fitness via pleiotropic mechanisms. Collectively, our results identify YhaJ as a pathotype-specific regulatory aide, enhancing the expression of key genes that are collectively required for UPEC bloodstream pathogenesis.

Strong pathogen competition in neonatal gut colonisation

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

Escherichia coli Strains from Patients with Inflammatory Bowel Diseases have Disease-specific Genomic Adaptations

BACKGROUND AND AIMS

Escherichia coli is over-abundant in the gut microbiome of patients with inflammatory bowel disease [IBD]. Here, we aimed to identify IBD-specific genomic functions of diverse E. coli lineages.

METHODS

We investigated E. coli genomes from patients with ulcerative colitis [UC], Crohn's disease [CD] or a pouch, and healthy subjects. The majority of genomes were reconstructed from metagenomic samples, including newly sequenced faecal metagenomes. Clinical metadata were collected. Functional analysis at the gene and mutation level were performed and integrated with IBD phenotypes and biomarkers.

RESULTS

Overall, 530 E. coli genomes were analysed. The E. coli B2 lineage was more prevalent in UC compared with other IBD phenotypes. Genomic metabolic capacities varied across E. coli lineages and IBD phenotypes. Host mucin utilisation enzymes were present in a single lineage and depleted in patients with a pouch, whereas those involved in inulin hydrolysis were enriched in patients with a pouch. E. coli strains from patients with UC were twice as likely to encode the genotoxic molecule colibactin than strains from patients with CD or a pouch. Strikingly, patients with a pouch showed the highest inferred E. coli growth rates, even in the presence of antibiotics. Faecal calprotectin did not correlate with the relative abundance of E. coli. Finally, we identified multiple IBD-specific non-synonymous mutations in E. coli genes encoding for bacterial cell envelope components.

CONCLUSIONS

Comparative genomics indicates that E. coli is a commensal species adapted to the overactive mucosal immune milieu in IBD, rather than causing it. Our results reveal mutations that may lead to attenuated antigenicity in some E. coli strains.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

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.

Carbapenemase-Producing Extraintestinal Pathogenic Escherichia coli From Argentina: Clonal Diversity and Predominance of Hyperepidemic Clones CC10 and CC131

Extraintestinal pathogenic Escherichia coli (ExPEC) causes infections outside the intestine. Particular ExPEC clones, such as clonal complex (CC)/sequence type (ST)131, have been known to sequentially accumulate antimicrobial resistance that starts with chromosomal mutations against fluoroquinolones, followed with the acquisition of bla_CTX–M–15 and, more recently, carbapenemases. Here we aimed to investigate the distribution of global epidemic clones of carbapenemase-producing ExPEC from Argentina in representative clinical isolates recovered between July 2008 and March 2017. Carbapenemase-producing ExPEC (n = 160) were referred to the Argentinean reference laboratory. Of these, 71 were selected for genome sequencing. Phenotypic and microbiological studies confirmed the presence of carbapenemases confirmed as KPC-2 (n = 52), NDM-1 (n = 16), IMP-8 (n = 2), and VIM-1 (n = 1) producers. The isolates had been recovered mainly from urine, blood, and abdominal fluids among others, and some were from screening samples. After analyzing the virulence gene content, 76% of the isolates were considered ExPEC, although non-ExPEC isolates were also obtained from extraintestinal sites. Pan-genome phylogeny and clonal analysis showed great clonal diversity, although the first phylogroup in abundance was phylogroup A, harboring CC10 isolates, followed by phylogroup B2 with CC/ST131, mostly H30Rx, the subclone co-producing CTX-M-15. Phylogroups D, B1, C, F, and E were also detected with fewer strains. CC10 and CC/ST131 were found throughout the country. In addition, CC10 nucleated most metalloenzymes, such as NDM-1. Other relevant international clones were identified, such as CC/ST38, CC155, CC14/ST1193, and CC23. Two isolates co-produced KPC-2 and OXA-163 or OXA-439, a point mutation variant of OXA-163, and three isolates co-produced MCR-1 among other resistance genes. To conclude, in this work, we described the molecular epidemiology of carbapenemase-producing ExPEC in Argentina. Further studies are necessary to determine the plasmid families disseminating carbapenemases in ExPEC in this region.

Genomic and Phenotypic Characterization of a Colistin-Resistant Escherichia coli Isolate Co-Harboring blaNDM-5, blaOXA-1, and blaCTX-M-55 Isolated from Urine

Background

Materials and Methods

Results

Conclusion

Considerations for Constructing a Protein Sequence Database for Metaproteomics

Mass spectrometry-based metaproteomics has emerged as a prominent technique for interrogating the functions of specific organisms in microbial communities, in addition to total community function. Identifying proteins by mass spectrometry requires matching mass spectra of fragmented peptide ions to a database of protein sequences corresponding to the proteins in the sample. This sequence database determines which protein sequences can be identified from the measurement, and as such the taxonomic and functional information that can be inferred from a metaproteomics measurement. Thus, the construction of the protein sequence database directly impacts the outcome of any metaproteomics study. Several factors, such as source of sequence information and database curation, need to be considered during database construction to maximize accurate protein identifications traceable to the species of origin. In this review, we provide an overview of existing strategies for database construction and the relevant studies that have sought to test and validate these strategies. Based on this review of the literature and our experience we provide a decision tree and best practices for choosing and implementing database construction strategies.

Genomic Epidemiology and Antimicrobial Susceptibility Profile of Enterotoxigenic Escherichia coli From Outpatients With Diarrhea in Shenzhen, China, 2015-2020

Enterotoxigenic Escherichia coli (ETEC) is the leading cause of severe diarrhea in children and the most common cause of diarrhea in travelers. However, most ETEC infections in Shenzhen, China were from indigenous adults. In this study, we characterized 106 ETEC isolates from indigenous outpatients with diarrhea (77% were adults aged >20 years) in Shenzhen between 2015 and 2020 by whole-genome sequencing and antimicrobial susceptibility testing. Shenzhen ETEC isolates showed a remarkable high diversity, which belonged to four E. coli phylogroups (A: 71%, B1: 13%, E: 10%, and D: 6%) and 15 ETEC lineages, with L11 (25%, O159:H34/O159:H43, ST218/ST3153), novel L2/4 (21%, O6:H16, ST48), and L4 (15%, O25:H16, ST1491) being major lineages. Heat-stable toxin (ST) was most prevalent (76%, STh: 60% STp: 16%), followed by heat-labile toxin (LT, 17%) and ST + LT (7%). One or multiple colonization factors (CFs) were identified in 68 (64%) isolates, with the common CFs being CS21 (48%) and CS6 (34%). Antimicrobial resistance mutation/gene profiles of genomes were concordant with the phenotype testing results of 52 representative isolates, which revealed high resistance rate to nalidixic acid (71%), ampicillin (69%), and ampicillin/sulbactam (46%), and demonstrated that the novel L2/4 was a multidrug-resistant lineage. This study provides novel insight into the genomic epidemiology and antimicrobial susceptibility profile of ETEC infections in indigenous adults for the first time, which further improves our understanding on ETEC epidemiology and has implications for the development of vaccine and future surveillance and prevention of ETEC infections.

Genomic and Phenotypic Analysis of Heat and Sanitizer Resistance in Escherichia coli from Beef in Relation to the Locus of Heat Resistance

The locus of heat resistance (LHR) can confer heat resistance to Escherichia coli to various extents. This study investigated the phylogenetic relationships and the genomic and phenotypic characteristics of E. coli with or without LHR recovered from beef by direct plating or from enrichment broth at 42°C. LHR-positive E. coli isolates (n = 24) were subjected to whole-genome sequencing by short and long reads. LHR-negative isolates (n = 18) from equivalent sources as LHR-positive isolates were short-read sequenced. All isolates were assessed for decimal reduction time at 60°C (D_60°C) and susceptibility to the sanitizers E-SAN and Perox-E. Selected isolates were evaluated for growth at 42°C. The LHR-positive and -negative isolates were well separated on the core genome tree, with 22/24 positive isolates clustering into three clades. Isolates within clade 1 and 2, despite their different D_60°C values, were clonal, as determined by subtyping (multilocus sequence typing [MLST], core genome MLST, and serotyping). Isolates within each clade are of one serotype. The LHR-negative isolates were genetically diverse. The LHR-positive isolates had a larger (P < 0.001) median genome size by 0.3 Mbp (5.0 versus 4.7 Mbp) and overrepresentation of genes related to plasmid maintenance, stress response, and cryptic prophages but underrepresentation of genes involved in epithelial attachment and virulence. All LHR-positive isolates harbored a chromosomal copy of LHR, and all clade 2 isolates had an additional partial copy of LHR on conjugative plasmids. The growth rates at 42°C were 0.71 ± 0.02 and 0.65 ± 0.02 log(OD) h^-1 for LHR-positive and -negative isolates, respectively. No meaningful difference in sanitizer susceptibility was noted between LHR-positive and -negative isolates. IMPORTANCE Resistant bacteria are serious food safety and public health concerns. Heat resistance conferred by the LHR varies largely among different strains of E. coli. The findings in this study show that genomic background and composition of LHR, in addition to the presence of LHR, play an important role in the degree of heat resistance in E. coli and that strains with certain genetic backgrounds are more likely to acquire and maintain the LHR. Also, caution should be exercised when recovering E. coli at elevated temperatures, as the presence of LHR may confer growth advantages to some strains. Interestingly, the LHR-harboring strains seem to have evolved further from their primary animal host to adapt to their secondary habitat, as reflected by fewer genes involved in virulence and epithelial attachment. The phylogenetic relationships among the isolates point toward multiple mechanisms for acquisition of LHR by E. coli, likely prior to its being deposited on meat.

Different evolutionary trends form the twilight zone of the bacterial pan-genome

The pan-genome is defined as the combined set of all genes in the gene pool of a species. Pan-genome analyses have been very useful in helping to understand different evolutionary dynamics of bacterial species: an open pan-genome often indicates a free-living lifestyle with metabolic versatility, while closed pan-genomes are linked to host-restricted, ecologically specialized bacteria. A detailed understanding of the species pan-genome has also been instrumental in tracking the phylodynamics of emerging drug resistance mechanisms and drug-resistant pathogens. However, current approaches to analyse a species’ pan-genome do not take the species population structure into account, nor do they account for the uneven sampling of different lineages, as is commonplace due to over-sampling of clinically relevant representatives. Here we present the application of a population structure-aware approach for classifying genes in a pan-genome based on within-species distribution. We demonstrate our approach on a collection of 7500 Escherichia coli genomes, one of the most-studied bacterial species and used as a model for an open pan-genome. We reveal clearly distinct groups of genes, clustered by different underlying evolutionary dynamics, and provide a more biologically informed and accurate description of the species’ pan-genome.

Machine Learning Prediction of Resistance to Subinhibitory Antimicrobial Concentrations from Escherichia coli Genomes

Escherichia coli is an important cause of bacterial infections worldwide, with multidrug-resistant strains incurring substantial costs on human lives. Besides therapeutic concentrations of antimicrobials in health care settings, the presence of subinhibitory antimicrobial residues in the environment and in clinics selects for antimicrobial resistance (AMR), but the underlying genetic repertoire is less well understood. Here, we used machine learning to predict the population doubling time and cell growth yield of 1,407 genetically diverse E. coli strains expanding under exposure to three subinhibitory concentrations of six classes of antimicrobials from single-nucleotide genetic variants, accessory gene variation, and the presence of known AMR genes. We predicted cell growth yields in the held-out test data with an average correlation (Spearman's ρ) of 0.63 (0.36 to 0.81 across concentrations) and cell doubling times with an average correlation of 0.59 (0.32 to 0.92 across concentrations), with moderate increases in sample size unlikely to improve predictions further. This finding points to the remaining missing heritability of growth under antimicrobial exposure being explained by effects that are too rare or weak to be captured unless sample size is dramatically increased, or by effects other than those conferred by the presence of individual single-nucleotide polymorphisms (SNPs) and genes. Predictions based on whole-genome information were generally superior to those based only on known AMR genes and were accurate for AMR resistance at therapeutic concentrations. We pinpointed genes and SNPs determining the predicted growth and thereby recapitulated many known AMR determinants. Finally, we estimated the effect sizes of resistance genes across the entire collection of strains, disclosing the growth effects for known resistance genes in each individual strain. Our results underscore the potential of predictive modeling of growth patterns from genomic data under subinhibitory concentrations of antimicrobials, although the remaining missing heritability poses a challenge for achieving the accuracy and precision required for clinical use. IMPORTANCE Predicting bacterial growth from genome sequences is important for a rapid characterization of strains in clinical diagnostics and to disclose candidate novel targets for anti-infective drugs. Previous studies have dissected the relationship between bacterial growth and genotype in mutant libraries for laboratory strains, yet no study so far has examined the predictive power of genome sequence in natural strains. In this study, we used a high-throughput phenotypic assay to measure the growth of a systematic collection of natural Escherichia coli strains and then employed machine learning models to predict bacterial growth from genomic data under nontherapeutic subinhibitory concentrations of antimicrobials that are common in nonclinical settings. We found a moderate to strong correlation between predicted and actual values for the different collected data sets. Moreover, we observed that the known resistance genes are still effective at sublethal concentrations, pointing to clinical implications of these concentrations.

Phylum barrier and Escherichia coli intra-species phylogeny drive the acquisition of antibiotic-resistance genes

Escherichia coli is a ubiquitous bacterium that has been widely exposed to antibiotics over the last 70 years. It has adapted by acquiring different antibiotic-resistance genes (ARGs), the census of which we aim to characterize here. To do so, we analysed 70 301 E. coli genomes obtained from the EnteroBase database and detected 1 027 651 ARGs using the AMRFinder, Mustard and ResfinderFG ARG databases. We observed a strong phylogroup and clonal lineage specific distribution of some ARGs, supporting the argument for epistasis between ARGs and the strain genetic background. However, each phylogroup had ARGs conferring a similar antibiotic class resistance pattern, indicating phenotypic adaptive convergence. The G+C content or the type of ARG was not associated with the frequency of the ARG in the database. In addition, we identified ARGs from anaerobic, non-Proteobacteria bacteria in four genomes of E. coli, supporting the hypothesis that the transfer between anaerobic bacteria and E. coli can spontaneously occur but remains exceptional. In conclusion, we showed that phylum barrier and intra-species phylogenetic history are major drivers of the acquisition of a resistome in E. coli.

Long-read-sequenced reference genomes of the seven major lineages of enterotoxigenic Escherichia coli (ETEC) circulating in modern time

Enterotoxigenic Escherichia coli (ETEC) is an enteric pathogen responsible for the majority of diarrheal cases worldwide. ETEC infections are estimated to cause 80,000 deaths annually, with the highest rates of burden, ca 75 million cases per year, amongst children under 5 years of age in resource-poor countries. It is also the leading cause of diarrhoea in travellers. Previous large-scale sequencing studies have found seven major ETEC lineages currently in circulation worldwide. We used PacBio long-read sequencing combined with Illumina sequencing to create high-quality complete reference genomes for each of the major lineages with manually curated chromosomes and plasmids. We confirm that the major ETEC lineages all harbour conserved plasmids that have been associated with their respective background genomes for decades, suggesting that the plasmids and chromosomes of ETEC are both crucial for ETEC virulence and success as pathogens. The in-depth analysis of gene content, synteny and correct annotations of plasmids will elucidate other plasmids with and without virulence factors in related bacterial species. These reference genomes allow for fast and accurate comparison between different ETEC strains, and these data will form the foundation of ETEC genomics research for years to come.