Combined Analysis of Variation in Core, Accessory and Regulatory Genome Regions Provides a Super-Resolution View into the Evolution of Bacterial Populations

Description of blaKPC-carrying Escherichia coli in patients from a Brazilian hospital over a 4-year period

Klebsiella pneumoniae carbapenemase (KPC)-producing Escherichia coli are recognized by the World Health Organization as a critical group of bacterial priority pathogens of public health importance. Thus, increased understanding of the genetic characteristics of KPC-producing E. coli is required. Here, we performed a retrospective study in a Brazilian teaching-hospital to describe the genomic features linked to antimicrobial resistance, virulence, and phylogeny of 40 meropenem-resistant E. coli. All isolates carried the blaKPC-2 gene, but amikacin, tigecycline, colistin, polymyxin B, and fosfomycin showed good activity. Molecular typing by MLST revealed 20 sequence types (STs), with a predominance of ST131. Whole-genome sequencing identified Tn4401 as a mechanism responsible for blaKPC-2 mobilization, a variety of antimicrobial resistance and virulence genes, the predominance of pathogenic phylogroup lineages, and the grouping of genomes belonging to the same ST. KPC-producing E. coli is not a common pathogen, but few treatment alternatives are available against potentially virulent strains.

Comparative Genomic Analysis of Two Vibrio harveyi Strains from Larimichthys crocea with Divergent Virulence Profiles

Vibrio harveyi is a significant pathogen in marine aquaculture, causing vibriosis in various marine species. This study presents a comparative genomic analysis of two V. harveyi strains, N8T11 and 45T2, which exhibit differing virulence profiles. Virulence assays revealed that N8T11 caused 92% mortality in infected fish, while 45T2 resulted in 0% mortality. Whole-genome sequencing revealed that strain N8T11 harbors five plasmids (pN8T11a, pN8T11b, pN8T11c, pN8T11d and pN8T11e) absent in 45T2, encoding genes potentially linked to virulence, such as siderophore-mediated iron acquisition and stress response mechanisms. Pan-genome analysis highlighted substantial genomic plasticity within V. harveyi, with mobile genetic elements, including plasmids and prophages, contributing to horizontal gene transfer. Conjugation experiments demonstrated that all five N8T11 plasmids can transfer to 45T2 with efficiencies up to 87%, with pN8T11b remaining stable across multiple subcultures, enabling the dissemination of virulence-associated genes. These findings suggest that plasmid-mediated gene transfer plays a role in the virulence variability observed between V. harveyi strains. This study contributes to understanding the genomic factors underlying pathogenicity in V. harveyi and provides insights for future research aimed at controlling vibriosis in aquaculture.

Machine learning reveals the dynamic importance of accessory sequences for Salmonella outbreak clustering

Bacterial typing at whole-genome scales is now feasible owing to decreasing costs in high-throughput sequencing and the recent advances in computation. The unprecedented resolution of whole-genome typing is achieved by genotyping the variable segments of bacterial genomes that can fluctuate significantly in gene content. However, due to the transient and hypervariable nature of many accessory elements, the value of the added resolution in outbreak investigations remains disputed. To assess the analytical value of bacterial accessory genomes in clustering epidemiologically related cases, we trained classifiers on a set of genomes collected from 24 Salmonella enterica outbreaks of food, animal, or environmental origin. The models demonstrated high precision and recall on unseen test data with near-perfect accuracy in classifying clonal and short-term outbreaks. Annotating the genomic features important for cluster classification revealed functional enrichment of molecular fingerprints in genes involved in membrane transportation, trafficking, and carbohydrate metabolism. Importantly, we discovered polymorphisms in mobile genetic elements (MGEs) and gain/loss of MGEs to be informative in defining outbreak clusters. To quantify the ability of MGE variations to cluster outbreak clones, we devised a reference-free tree-building algorithm inspired by colored de Bruijn graphs, which enabled topological comparisons between MGE and standard typing methods. Systematic evaluation of clustering MGEs on an unseen dataset of 34 Salmonella outbreaks yielded mixed results that exemplified the power of accessory sequence variations when core genomes of unrelated cases are insufficiently discriminatory, as well as the distortion of outbreak signals by microevolution events or the incomplete assembly of MGEs.

IMPORTANCE

Gene-by-gene typing is widely used to detect clusters of foodborne illnesses that share a common origin. It remains actively debated whether the inclusion of accessory sequences in bacterial typing schema is informative or deleterious for cluster definitions in outbreak investigations due to the potential confounding effects of horizontal gene transfer. By training machine learning models on a curated set of historical Salmonella outbreaks, we revealed an enriched presence of outbreak distinguishing features in a wide range of mobile genetic elements. Systematic comparison of the efficacy of clustering different accessory elements against standard sequence typing methods led to our cataloging of scenarios where accessory sequence variations were beneficial and uninformative to resolving outbreak clusters. The presented work underscores the complexity of the molecular trends in enteric outbreaks and seeks to inspire novel computational ways to exploit whole-genome sequencing data in enteric disease surveillance and management.

High intra-laboratory reproducibility of nanopore sequencing in bacterial species underscores advances in its accuracy

Nanopore sequencing is a third-generation technology known for its portability, real-time analysis and ability to generate long reads. It has great potential for use in clinical diagnostics, but thorough validation is required to address accuracy concerns and ensure reliable and reproducible results. In this study, we automated an open-source workflow (freely available at https://gitlab.com/FLI_Bioinfo/nanobacta) for the assembly of Oxford Nanopore sequencing data and used it to investigate the reproducibility of assembly results under consistent conditions. We used a benchmark dataset of five bacterial reference strains and generated eight technical sequencing replicates of the same DNA using the Ligation and Rapid Barcoding kits together with the Flongle and MinION flow cells. We assessed reproducibility by measuring discrepancies such as substitution and insertion/deletion errors, analysing plasmid recovery results and examining genetic markers and clustering information. We compared the results of genome assemblies with and without short-read polishing. Our results show an average reproducibility accuracy of 99.999955% for nanopore-only assemblies and 99.999996% when the short reads were used for polishing. The genomic analysis results were highly reproducible for the nanopore-only assemblies without short read in the following areas: identification of genetic markers for antimicrobial resistance and virulence, classical MLST, taxonomic classification, genome completeness and contamination analysis. Interestingly, the clustering information results from the core genome SNP and core genome MLST analyses were also highly reproducible for the nanopore-only assemblies, with pairwise differences of up to two allele differences in core genome MLST and two SNPs in core genome SNP across replicates. After polishing the assemblies with short reads, the pairwise differences for cgMLST were 0 and for cgSNP were 0-1 SNP across replicates. These results highlight the advances in sequencing accuracy of nanopore data without the use of short reads.

Characterization of antibiotic determinants and heavy metal resistance genes in Escherichia coli from pigs in Catalonia

More antibiotics are administered to livestock animals than to treat human infections. Industrialization, large animal densities and early weaning mean pigs are exposed to more antibiotics than any other livestock animal. Consequently, antimicrobial resistance (AMR) is common among commensal and pathogenic bacteria. Heavy metals (HMs) are also often used as feed additives for growth promotion and infection prevention alongside antimicrobials, and increased exposure to copper, zinc and cadmium can further encourage AMR through co-selection. In this study, we sequenced an archived collection of 112 Escherichia coli isolates from pigs in Catalonia using short- and long-read sequencing methods to detect AMR and HM tolerance genes. The most common AMR genes were mdfA (84.8%), aph(3″)-Ib (52.7%), bla TEM-1B (45.6%) and aph(6)-Id (45.6%). Genes relevant to public health, such as the extended-spectrum β-lactamases (15.4%), bla CTX-M type or bla SHV, or mobile colistin resistance (mcr) genes (13.4%), such as mcr-1, were also found. HM tolerance genes were present in almost every genome but were rarely located in plasmids, and, in most cases, AMR and HM tolerance genes were not located on the same plasmids. Of the genes predicted to increase tolerance to HMs, only those with activity to mercury were co-located on plasmids alongside other AMR determinants. However, mercury is rarely used in pig farming and does not support a scenario where AMR and HM genes are co-selected. Finally, we identified the exclusive association between mcr-4 and ColE10 plasmid, which may help target interventions to curtail its spread among pig Escherichia coli.

Comparative Pan- and Phylo-Genomic Analysis of Ideonella and Thermobifida Strains: Dissemination of Biodegradation Potential and Genomic Divergence

Ideonella and Thermobifida were the most promising bacterial candidates for degrading plastic polymers. A comparative pan- and phylogenomic analysis of 33 Ideonella and Thermobifida strains was done to determine their plastic degradation potential, niche adaptation and speciation. Our study disclosed that more accessory genes in the strains showed phenotypic plasticity, according to the BPGA data. Pan and core genes were employed for the phylogenetic reconstruction. Pathway enrichment analyses scrutinized the functional roles of the core and adaptive-associated genes. KEGG annotation revealed that most genes were associated with the metabolism of amino acids and carbohydrates. The detailed COG analysis disclosed that approximately 40% of the pan genes performed metabolic functions. The unique gene pool consisted of genes chiefly involved in "general function prediction" and "amino acid transport and metabolism". Our in silico study revealed that these strains could assist in agronomic applications in the future since they devour nitrogen compounds and their central metabolic pathways are involved in amino acid metabolism. The rational selection of strains of Ideonella is far more effective at depolymerising plastics than Thermobifida. A greater number of unique genes, 1701 and 692, were identified for Ideonella sakaiensis 201-F6 and Thermobifida alba DSM-43795, respectively. Furthermore, we examined the singletons involved in xenobiotic catabolism. The unique singleton data were used to construct a supertree. To characterize the conserved patterns, we used SMART and MEME to identify domain and transmembrane regions in the unique protein sequences. Therefore, our study unraveled the genomic insights into the ecology-driven speciation of Ideonella and Thermobifida.

Linkage-based ortholog refinement in bacterial pangenomes with CLARC

Bacterial genomes exhibit significant variation in gene content and sequence identity. Pangenome analyses explore this diversity by classifying genes into core and accessory clusters of orthologous groups (COGs). However, strict sequence identity cutoffs can misclassify divergent alleles as different genes, inflating accessory gene counts. CLARC (Connected Linkage and Alignment Redefinition of COGs) [https://github.com/IndraGonz/CLARC] improves pangenome analyses by condensing accessory COGs using functional annotation and linkage information. Through this approach, orthologous groups are consolidated into more practical units of selection. Analyzing 8,000+ Streptococcus pneumoniae genomes, CLARC reduced accessory gene estimates by more than 30% and improved evolutionary predictions based on accessory gene frequencies. By refining COG definitions, CLARC offers critical insights into bacterial evolution, aiding genetic studies across diverse populations.

Predicting sepsis mortality into an era of pandrug-resistant E. coli through modeling

BACKGROUND

Infections caused by antibiotic-resistant bacteria are increasingly frequent, burdening healthcare systems worldwide. As pathogens acquire resistance to all known antibiotics - i.e., become pan-resistant - treatment of the associated infections will become exceedingly difficult. We hypothesized that the emergence of pan-resistant bacterial pathogens will result in a sharp increase in human mortality.

METHODS

We tested this hypothesis by modeling the impact of a single hypothetical pan-resistant Escherichia coli strain on sepsis deaths in the United States. We used long-term data on sepsis incidence, mortality rates, strain dynamics, and treatment outcomes to parameterize a set of models encompassing a range of plausible future scenarios. All models accounted for historical and projected temporal changes in population size and age distribution.

RESULTS

The models suggest that sepsis deaths could increase 18- to 46-fold within 5 years of the emergence of a single pan-resistant E. coli strain. This large and rapid change contrasts sharply with the current expectation of gradual change under continuing multidrug-resistance.

CONCLUSIONS

Failure to prevent the emergence of pan-resistance would have dire consequences for public health.

Genomic and growth fitness study of extended-spectrum β-lactamase-producing Escherichia coli from bloodstream infections after introduction of a national 4C antimicrobial stewardship policy in Scotland

BACKGROUND

Extended-spectrum β-lactamase-producing Escherichia coli remains a major cause of hospital-acquired bloodstream infections in countries with high antimicrobial stewardship compliance.

METHODS

Isolates from bloodstream infections that occurred between 2010 and 2020 in a tertiary-level hospital in North-East Scotland soon after introduction of the '4C' antimicrobial stewardship policy were analysed for phylogenetic structure, antimicrobial resistance, plasmid, and virulence gene carriage. Growth fitness was measured in kinetic assays. Non-metric-multidimensional-scaling was used to evaluate clonal relationships, antimicrobial resistance, and virulence profiles in early and later years after the 4C policy introduction. Clonal and fitness trends over the study period were determined by generalised additive modelling. The relationship between clonal type, antimicrobial resistance, and fitness was evaluated by linear regression.

RESULTS

Three hierarchical phylogenetic clusters were identified, with the most dominant cluster (O25:H4/fimH30) including all, and nearly exclusively, Clade C ST131 isolates as well as minor non-ST131 sequence types. The prevalence of ST131 was largely stable over the study period. Resistance to aminoglycosides and aztreonam in ST131 was lower (P = 0.019 and P = 0.004, respectively) during later years (2016-2020) by 28% on average compared to early years soon after 4C policy implementation (2010-2014). Carriage of virulence factors involved in bacterial adaptation was higher in ST131 compared to non-ST131 but mostly stable in early vs. later years. Growth fitness of ST131 was lower than non-ST131 and declined steadily in later years (P < 0.0001).

CONCLUSIONS

Despite stable virulence factor carriage, population structure, and resistance to cephalosporins, we show increased susceptibility of ST131 to aminoglycosides and aztreonam and concurrent fitness decline years after the introduction of the 4C policy.

Mobile genetic elements define the non-random structure of the Salmonella enterica serovar Typhi pangenome

Bacterial relatedness measured using select chromosomal loci forms the basis of public health genomic surveillance. While approximating vertical evolution through this approach has proven exceptionally valuable for understanding pathogen dynamics, it excludes a fundamental dimension of bacterial evolution-horizontal gene transfer. Incorporating the accessory genome is the logical remediation and has recently shown promise in expanding epidemiological resolution for enteric pathogens. Employing k-mer-based Jaccard index analysis, and a novel genome length distance metric, we computed pangenome (i.e., core and accessory) relatedness for the globally important pathogen Salmonella enterica serotype Typhi (Typhi), and graphically express both vertical (homology-by-descent) and horizontal (homology-by-admixture) evolutionary relationships in a reticulate network of over 2,200 U.S. Typhi genomes. This analysis revealed non-random structure in the Typhi pangenome that is driven predominantly by the gain and loss of mobile genetic elements, confirming and expanding upon known epidemiological patterns, revealing novel plasmid dynamics, and identifying avenues for further genomic epidemiological exploration. With an eye to public health application, this work adds important biological context to the rapidly improving ways of analyzing bacterial genetic data and demonstrates the value of the accessory genome to infer pathogen epidemiology and evolution.IMPORTANCEGiven bacterial evolution occurs in both vertical and horizontal dimensions, inclusion of both core and accessory genetic material (i.e., the pangenome) is a logical step toward a more thorough understanding of pathogen dynamics. With an eye to public, and indeed, global health relevance, we couple contemporary tools for genomic analysis with decades of research on mobile genetic elements to demonstrate the value of the pangenome, known and unknown, annotated, and hypothetical, for stratification of Salmonella enterica serovar Typhi (Typhi) populations. We confirm and expand upon what is known about Typhi epidemiology, plasmids, and antimicrobial resistance dynamics, and offer new avenues of exploration to further deduce Typhi ecology and evolution, and ultimately to reduce the incidence of human disease.

Longitudinal molecular analysis of clinical and fecal Escherichia coli isolates at a Veterans Affairs Medical Center in Minnesota, USA, 2012-2019

Introduction

Extraintestinal Escherichia coli infections represent a growing public health threat, However, current studies often overlook important factors such as temporal patterns of infection, phylogenetic and clonal background, or the host gut E. coli population, despite their likely significance.

Methods

In this study, we analyzed >7000 clinical E. coli isolates from patients at the Minneapolis Veterans Affairs Health Care System (2012-2019), and concurrent fecal E. coli from uninfected veterans. We assessed phylogenetic group distribution, membership in selected sequence types (STs), and subsets thereof-including the pandemic, resistance-associated ST131-H30R, and ST1193 lineages-and strain type, as defined by pulsed-field gel electrophoresis. We then analyzed these features alongside the temporal patterns of infection in individual hosts.

Results

The H30R lineage emerged as the leading lineage, both overall and among fluoroquinolone-resistant isolates, with ST1193 following among fluoroquinolone-resistant isolates. Recurrences were common, occurring in 31% of subjects and 41% of episodes, and often multiple and delayed/prolonged (up to 23 episodes per subject; up to 2655d post-index). Remarkably, these recurrences typically involved the subject's index strain (63% of recurrences), even when affecting extra-urinary sites. ST131, H30R, ST1193, and fluoroquinolone-resistant strains generally caused significantly more recurrences than did other strains, despite similar recurrence intervals. ST131 strain types shifted significantly over the study period. Infection-causing strains were commonly detectable in host feces at times other than during an infection episode; the likelihood of detection varied with surveillance intensity and proximity to the infection. H30R and ST1193 were prominent causes of fecal-clinical clonal overlap.

Discussion

These findings provide novel insights into the temporal and clonal characteristics of E. coli infections in veterans and support efforts to develop anti-colonization interventions.

Evolutionary dynamics of the accessory genomes of Staphylococcus aureus

Staphylococcus aureus is a ubiquitous commensal and opportunistic bacterial pathogen that can cause a wide gamut of infections, which are exacerbated by the presence of multidrug-resistant and methicillin-resistant S. aureus. S. aureus is genetically heterogeneous and consists of numerous distinct lineages. Using 558 complete genomes of S. aureus, we aim to determine how the accessory genome content among phylogenetic lineages of S. aureus is structured and has evolved. Bayesian hierarchical clustering identified 10 sequence clusters, of which seven contained major sequence types (ST 1, 5, 8, 30, 59, 239, and 398). The seven sequence clusters differed in their accessory gene content, including genes associated with antimicrobial resistance and virulence. Focusing on the two largest clusters, BAPS8 and BAPS10, and each consisting mostly of ST5 and ST8, respectively, we found that the structure and connected components in the co-occurrence networks of accessory genomes varied between them. These differences are explained, in part, by the variation in the rates at which the two sequence clusters gained and lost accessory genes, with the highest rate of gene accumulation occurring recently in their evolutionary histories. We also identified a divergent group within BAPS10 that has experienced high gene gain and loss early in its history. Together, our results show highly variable and dynamic accessory genomes in S. aureus that are structured by the history of the specific lineages that carry them.IMPORTANCEStaphylococcus aureus is an opportunistic, multi-host pathogen that can cause a variety of benign and life-threatening infections. Our results revealed considerable differences in the structure and evolution of the accessory genomes of major lineages within S. aureus. Such genomic variation within a species can have important implications on disease epidemiology, pathogenesis of infection, and interactions with the vertebrate host. Our findings provide important insights into the underlying genetic basis for the success of S. aureus as a highly adaptable and resistant pathogen, which will inform current efforts to control and treat staphylococcal diseases.

Sequence-subtype association of multi-drug-resistant diarrheagenic Escherichia coli

Background and Objectives

Multi-drug-resistant pathogens pose a significant threat as they can rapidly spread, leading to severe healthcare-associated invasive infections. In developing countries, diarrheagenic Escherichia coli (DEC) is a major bacterial pathogen responsible for causing diarrhea. However, the outbreak of resistant strains has made the treatment of DEC infections much more challenging. This study aimed to investigate the relationship between antibiotic resistance genes and other virulence categories in E. coli strains that cause diarrhea, particularly DEC.

Materials and Methods

The phylogenetic grouping was defined using PCR and multi-locus sequence type (MLST) methods.

Results

Among the isolates analyzed, 14 were identified as resistant and were classified into eight distinct sequence types: ST3, ST53, ST77, ST483, ST512, ST636, ST833, and ST774, indicating genetic diversity among the resistant strains. Certain sequence types, notably ST512 and ST636, were found to be associated with multiple antibiotic resistance in DEC. Regarding antibiotic susceptibility, strains showed the highest resistance to amoxicillin, suggesting that this antibiotic may not be effective in treating DEC infections. On the other hand, the isolates demonstrated susceptibility to amikacin and chloramphenicol, implying that these antibiotics could be more suitable treatment options for DEC infections.

Conclusion

The findings underscore the importance of promptly identifying antibiotic resistance patterns and their correlation with specific pathogenic virulence categories, as this knowledge can aid in selecting the most appropriate antibiotics for treating DEC infections. Considering the antibiotic resistance profiles and associated resistance genes is crucial in managing and containing diarrheal outbreaks and in selecting effective antibiotic therapies for DEC infections.

De novo mutations mediate phenotypic switching in an opportunistic human lung pathogen

Bacteria evolving within human hosts encounter selective tradeoffs that render mutations adaptive in one context and deleterious in another. Here, we report that the cystic fibrosis-associated pathogen Burkholderia dolosa overcomes in-human selective tradeoffs by acquiring successive point mutations that alternate phenotypes. We sequenced the whole genomes of 931 respiratory isolates from two recently infected patients and an epidemiologically-linked, chronically-infected patient. These isolates are contextualized using 112 historical genomes from the same outbreak strain. Within both newly infected patients, diverse parallel mutations that disrupt O-antigen expression quickly arose, comprising 29% and 63% of their B. dolosa communities by 3 years. The selection for loss of O-antigen starkly contrasts with our previous observation of parallel O-antigen-restoring mutations after many years of chronic infection in the historical outbreak. Experimental characterization revealed that O-antigen loss increases uptake in immune cells while decreasing competitiveness in the mouse lung. We propose that the balance of these pressures, and thus whether O-antigen expression is advantageous, depends on tissue localization and infection duration. These results suggest that mutation-driven alternation during infection may be more frequent than appreciated and is underestimated without dense temporal sampling.

High-risk lineages of extended spectrum cephalosporinase producing Escherichia coli from Eurasian griffon vultures (Gyps fulvus) foraging in landfills in north-eastern Spain

Extended-spectrum cephalosporinase producing (ESC) E. coli are regarded as key indicator microorganisms of antimicrobial resistance (AMR), calling for a One Health integrated global surveillance strategy. Wildlife is exposed to antibiotic contaminants and/or resistant bacteria that have been released into the environment, potentially acting as reservoirs and spreaders of resistance genes as well as sentinels of anthropogenic pressure. Monitoring AMR in wildlife has become crucial in determining anthropogenic environmental impacts as well as transmission routes. In this study, we determined the occurrence and potential sources of ESC E. coli in 218 Eurasian griffon vultures (Gyps fulvus) foraging regularly on human waste disposed at a dumpsite in north-eastern Spain. Minimal inhibitory concentration for 14 different antimicrobials was performed to evaluate the phenotype of the isolates, and whole genome sequencing was carried out to investigate lineages and plasmids harbouring ESC genes. Our sequences were compared to previously published Spanish sequences of human, animal, and wildlife origin. We report a high prevalence of CTX-M-15, as well as the presence of other resistance genes such as OXA-10, CTX-M-27, and CTX-M-65 which are rarely described in European livestock, suggesting a human origin. The isolates also carried a diverse range of additional AMR genes for a broad spectrum of drug families, with the majority being multi-drug resistant. The phylogenomic analyses suggests the transmission of high-risk lineages from humans to vultures, with 49 % of our isolates matching the most common extraintestinal pathogenic E. coli (ExPEC) lineages described in humans worldwide, including ST131, ST10 and ST58. We conclude that anthropogenically altered habitats, such as landfills, are hotspots for the acquisition and spread of high-risk ESC E. coli lineages associated with hospital infections. Measures must be implemented to limit their spread into natural environments.

Laboratory Evolution of Antimicrobial Resistance in Bacteria to Develop Rational Treatment Strategies

Laboratory evolution studies, particularly with Escherichia coli, have yielded invaluable insights into the mechanisms of antimicrobial resistance (AMR). Recent investigations have illuminated that, with repetitive antibiotic exposures, bacterial populations will adapt and eventually become tolerant and resistant to the drugs. Through intensive analyses, these inquiries have unveiled instances of convergent evolution across diverse antibiotics, the pleiotropic effects of resistance mutations, and the role played by loss-of-function mutations in the evolutionary landscape. Moreover, a quantitative analysis of multidrug combinations has shed light on collateral sensitivity, revealing specific drug combinations capable of suppressing the acquisition of resistance. This review article introduces the methodologies employed in the laboratory evolution of AMR in bacteria and presents recent discoveries concerning AMR mechanisms derived from laboratory evolution. Additionally, the review outlines the application of laboratory evolution in endeavors to formulate rational treatment strategies.

Commercial farmed swine harbour a variety of pathogenic bacteria and antimicrobial resistance genes

Introduction. The northern region of Thailand serves as a crucial area for swine production, contributing to the Thai community food supply. Previous studies have highlighted the presence of foodborne bacterial pathogens originating from swine farms in this region, posing a threat to both human and animal health.Gap statement. Multiple swine bacterial pathogens have been studied at a species level, but the distribution and co-occurrence of bacterial pathogens in agricultural swine has not been well established.Aim. Our study employed the intestinal scraping technique to directly examine the bacterial micro-organisms interacting with the swine host.Methodology. We used shotgun metagenomic sequencing to analyse the bacterial pathogens inhabiting the caecal microbiome of swine from five commercial farms in northern Thailand.Results. A variety of pathogenic and opportunistic bacteria were identified, including Escherichia coli, Clostridium botulinum, Staphylococcus aureus and the Corynebacterium genus. From a One Health perspective, these species are important foodborne and opportunistic pathogens in both humans and agricultural animals, making swine a critical pathogen reservoir that can cause illness in humans, especially farm workers. Additionally, the swine caecal microbiome contains commensal bacteria such as Bifidobacterium, Lactobacillus and Faecalibacterium, which are associated with normal physiology and feed utilization in healthy swine. Antimicrobial resistance genes were also detected in all samples, specifically conferring resistance to tetracycline and aminoglycosides, which have historically been used extensively in swine farming.Conclusion. The findings further support the need for improved sanitation standards in swine farms, and additional monitoring of agricultural animals and farm workers to reduce contamination and improved produce safety for human consumption.

Functional diversification despite structural congruence in the HipBST toxin-antitoxin system of Legionella pneumophila

IMPORTANCE

Toxin-antitoxin (TA) systems are parasitic genetic elements found in almost all bacterial genomes. They are exchanged horizontally between cells and are typically poorly conserved across closely related strains and species. Here, we report the characterization of a tripartite TA system in the bacterial pathogen Legionella pneumophila that is highly conserved across Legionella species genomes. This system (denoted HipBSTLp) is a distant homolog of the recently discovered split-HipA system in Escherichia coli (HipBSTEc). We present bioinformatic, molecular, and structural analyses of the divergence between these two systems and the functionality of this newly described TA system family. Furthermore, we provide evidence to refute previous claims that the toxin in this system (HipTLp) possesses bifunctionality as an L. pneumophila virulence protein. Overall, this work expands our understanding of the split-HipA system architecture and illustrates the potential for undiscovered biology in these abundant genetic elements.

Multidrug-resistant E. coli encoding high genetic diversity in carbohydrate metabolism genes displace commensal E. coli from the intestinal tract

Extra-intestinal pathogenic Escherichia coli (ExPEC) can cause a variety of infections outside of the intestine and are a major causative agent of urinary tract infections. Treatment of these infections is increasingly frustrated by antimicrobial resistance (AMR) diminishing the number of effective therapies available to clinicians. Incidence of multidrug resistance (MDR) is not uniform across the phylogenetic spectrum of E. coli. Instead, AMR is concentrated in select lineages, such as ST131, which are MDR pandemic clones that have spread AMR globally. Using a gnotobiotic mouse model, we demonstrate that an MDR E. coli ST131 is capable of out-competing and displacing non-MDR E. coli from the gut in vivo. This is achieved in the absence of antibiotic treatment mediating a selective advantage. In mice colonised with non-MDR E. coli strains, challenge with MDR E. coli either by oral gavage or co-housing with MDR E. coli colonised mice results in displacement and dominant intestinal colonisation by MDR E. coli ST131. To investigate the genetic basis of this superior gut colonisation ability by MDR E. coli, we assayed the metabolic capabilities of our strains using a Biolog phenotypic microarray revealing altered carbon metabolism. Functional pangenomic analysis of 19,571 E. coli genomes revealed that carriage of AMR genes is associated with increased diversity in carbohydrate metabolism genes. The data presented here demonstrate that independent of antibiotic selective pressures, MDR E. coli display a competitive advantage to colonise the mammalian gut and points to a vital role of metabolism in the evolution and success of MDR lineages of E. coli via carriage and spread.

Genomic and functional characterization of carbapenem-resistant Klebsiella pneumoniae from hospital wastewater

BACKGROUND

The emergence of carbapenem-resistant Klebsiella pneumoniae (CRKP) attracted extensive attention. Information on CRKP from hospital wastewater (HWW) is limited. The aims of this study were to investigate the genomic characteristics and to evaluate the survivability characteristics of 11 CRKP from HWW in a Chinese teaching hospital in Fujian province.

RESULTS

A total of 11 CRKP from HWW were recovered in this study. All CRKP from HWW were resistant to most antibiotics. Comparative genetic analysis demonstrated that all CRKP isolates were clustered into the three distinct phylogenetic clades and clade 2 and clade 3 were mixtures of samples collected from both HWW and clinical settings. Varieties of resistance genes, virulence genes and plasmid replicon types were detected in CRKP from HWW. In vitro transfer of bla_KPC-2 was successful for 3 bla_KPC-2-positive CRKP from HWW with high conjugation frequency. Our study demonstrated that the genetic environments of bla_KPC-2 shared core structure with ISKpn27-bla_KPC-2-ISKpn6. Group analysis showed that CRKP from HWW had a lower survivability in serum compared to clinical CRKP (p < 005); and CRKP from HWW had no significant difference in survivability in HWW compared to clinical CRKP (p > 005).

CONCLUSIONS

We analyzed the genomic and survivability characteristics of CRKP from HWW in a Chinese teaching hospital. These genomes represent a significant addition of genomic data from the genus and could serve as a valuable resource for future genomic studies about CRKP from HWW.

Transiently silent acquired antimicrobial resistance: an emerging challenge in susceptibility testing

Acquisition and expression of antimicrobial resistance (AMR) mechanisms in bacteria are often associated with a fitness cost. Thus, evolutionary adaptation and fitness cost compensation may support the advance of subpopulations with a silent resistance phenotype when the antibiotic selection pressure is absent. However, reports are emerging on the transient nature of silent acquired AMR, describing genetic alterations that can change the expression of these determinants to a clinically relevant level of resistance, and the association with breakthrough infections causing treatment failures. This phenomenon of transiently silent acquired AMR (tsaAMR) is likely to increase, considering the overall expansion of acquired AMR in bacterial pathogens. Moreover, the augmented use of genotypic methods in combination with conventional phenotypic antimicrobial susceptibility testing (AST) will increasingly enable the detection of genotype and phenotype discrepancy. This review defines tsaAMR as acquired antimicrobial resistance genes with a corresponding phenotype within the wild-type distribution or below the clinical breakpoint for susceptibility for which genetic alterations can mediate expression to a clinically relevant level of resistance. References to in vivo resistance development and therapeutic failures caused by selected resistant subpopulations of tsaAMR in Gram-positive and Gram-negative pathogens are given. We also describe the underlying molecular mechanisms, including alterations in the expression, reading frame or copy number of AMR determinants, and discuss the clinical relevance concerning challenges for conventional AST.

Genomic epidemiology of clinical ESBL-producing Enterobacteriaceae in a German hospital suggests infections are primarily community- and regionally-acquired

Clinical Enterobacteriaceae isolates that produce extended-spectrum β-lactamases (ESBLs) have been increasingly reported at a global scale. However, comprehensive data on the molecular epidemiology of ESBL-producing strains are limited and few studies have been conducted in non-outbreak situations.We used whole-genome sequencing to describe the population structure of 294 ESBL-producing Escherichia coli and Klebsiella pneumoniae isolates that were recovered from a German community hospital throughout a 1 year sampling period in a non-outbreak situation.We found a high proportion of E. coli isolates (61.5 %) belonged to the globally disseminated extraintestinal pathogenic ST131, whereas a wider diversity of STs was observed among K. pneumoniae isolates. The E. coli ST131 population in this study was shaped by multiple introductions of strains as demonstrated by contextual genomic analysis including ST131 strains from other geographical sources. While no recent common ancestor of the isolates of the current study and other international isolates was found, our clinical isolates clustered with those previously recovered in the region. Furthermore, we found that the isolation of ESBL-producing clinical strains in hospitalized patients could only rarely be associated with likely patient-to-patient transmission, indicating primarily a community and regional acquisition of strains.Further genomic analyses of clinical, carriage and environmental isolates is needed to uncover hidden transmissions and thus discover the most common sources of ESBL-producing pathogen infections in our hospitals.

Genomics Analysis to Identify Multiple Genetic Determinants That Drive the Global Transmission of the Pandemic ST95 Lineage of Extraintestinal Pathogenic Escherichia coli (ExPEC)

Extraintestinal pathogenic Escherichia coli (ExPEC) is a pathogen that causes host extraintestinal diseases. The ST95 E. coli lineage is one of the dominant ExPEC lineages in humans and poultry. In this study, we took advantage of extensive E. coli genomes available through public open-access databases to construct a detailed understanding of the phylogeny and evolution of ST95. We used a high variability of accessory genomes to highlight the diversity and dynamic traits of ST95. Isolates from diverse hosts and geographic sources were randomly located on the phylogenetic tree, which suggested that there is no host specificity for ST95. The time-scaled phylogeny showed that ST95 is an ancient and long-lasting lineage. The virulence genes, resistance genes, and pathogenicity islands (PAIs) were characterized in ST95 pan-genomes to provide novel insights into the pathogenicity and multidrug resistance (MDR) genotypes. We found that a pool of large plasmids drives virulence and MDR. Based on the unique genes in the ST95 pan-genome, we designed a novel multiplex PCR reaction to rapidly detect ST95. Overall, our study addressed a gap in the current understanding of ST95 ExPEC genomes, with significant implications for recognizing the success and spread of ST95.

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.

Informing plasmid compatibility with bacterial hosts using protein-protein interaction data

The compatibility of plasmids with new host cells is significant given their role in spreading antimicrobial resistance (AMR) and virulence factor genes. Evaluating this using in vitro screening is laborious and can be informed by computational analyses of plasmid-host compatibility through rates of protein-protein interactions (PPIs) between plasmid and host cell proteins. We identified large excesses of such PPIs in eight important plasmids, including pOXA-48, using most known bacteria (n = 4363). 23 species had high rates of interactions with four blaOXA-48-positive plasmids. We also identified 48 species with high interaction rates with plasmids common in Escherichia coli. We found a strong association between one plasmid and the fimbrial adhesin operon pil, which could enhance host cell adhesion in aqueous environments. An excess rate of PPIs could be a sign of host-plasmid compatibility, which is important for AMR control given that plasmids like pOXA-48 move between species with ease.

Complete genomic analysis of ST117 lineage extraintestinal pathogenic Escherichia coli (ExPEC) to reveal multiple genetic determinants to drive its global transmission: ST117 E. coli as an emerging multidrug-resistant foodborne ExPEC with zoonotic potential

Avian pathogenic Escherichia coli (APEC) is recognized as a primary source of foodborne extraintestinal pathogenic E. coli (ExPEC), which poses a significant risk of extraintestinal infections in humans. The potential of human infection with ST117 lineage APEC/ExPEC from poultry is particularly concerning. However, relatively few whole-genome studies have focused on ST117 as an emerging ExPEC lineage. In this study, the complete genomes of 11 avian ST117 isolates and the draft genomes of 20 ST117 isolates in China were sequenced to reveal the genomic islands and large plasmid composition of ST117 APEC. With reference to the extensive E. coli genomes available in public databases, large-scale comprehensive genomic analysis of the ST117 lineage APEC/ExPEC was performed to reveal the features of the ST117 pan-genome and population. The high variability of the accessory genome emphasized the diversity and dynamic traits of the ST117 pan-genome. ST117 isolates recovered from different hosts and geographic sources were randomly located on a phylogeny tree, suggesting that ST117 E. coli lacked host specificity. A time-scaled phylogeny tree showed that ST117 was a recent E. coli lineage with a relatively short evolutionary period. Further characterization of a wide diversity of ExPEC-related virulence genes, pathogenicity islands (PAIs), and resistance genes of the ST117 pan-genome provided insights into the virulence and resistance of ST117 APEC/ExPEC. The results suggested zoonotic potential of ST117 APEC/ExPEC between birds and humans. Moreover, genomic analysis showed that a pool of diverse plasmids drove the virulence and multidrug resistance of ST117 APEC/ExPEC. Several types of large plasmids were scattered across the ST117 isolates, but there was no strong plasmid-clade adaptation. Combined with the pan-genome analysis, a double polymerase chain reaction (PCR) method was designed for rapid and cost-effective detection of ST117 isolates from various avian and human APEC/ExPEC isolates. Overall, this study addressed a gap in current knowledge about the ST117 APEC/ExPEC genome, with significant implications to understand the success and spread of ST117 APEC/ExPEC.

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.

Whole-Genome Sequence Comparisons of Listeria monocytogenes Isolated from Meat and Fish Reveal High Inter- and Intra-Sample Diversity

Interpretation of whole-genome sequencing (WGS) data for foodborne outbreak investigations is complex, as the genetic diversity within processing plants and transmission events need to be considered. In this study, we analyzed 92 food-associated Listeria monocytogenes isolates by WGS-based methods. We aimed to examine the genetic diversity within meat and fish production chains and to assess the applicability of suggested thresholds for clustering of potentially related isolates. Therefore, meat-associated isolates originating from the same samples or processing plants as well as fish-associated isolates were analyzed as distinct sets. In silico serogrouping, multilocus sequence typing (MLST), core genome MLST (cgMLST), and pangenome analysis were combined with screenings for prophages and genetic traits. Isolates of the same subtypes (cgMLST types (CTs) or MLST sequence types (STs)) were additionally compared by SNP calling. This revealed the occurrence of more than one CT within all three investigated plants and within two samples. Analysis of the fish set resulted in predominant assignment of isolates from pangasius catfish and salmon to ST2 and ST121, respectively, potentially indicating persistence within the respective production chains. The approach not only allowed the detection of distinct subtypes but also the determination of differences between closely related isolates, which need to be considered when interpreting WGS data for surveillance.

Assessing the Load, Virulence and Antibiotic-Resistant Traits of ESBL/Ampc E. coli from Broilers Raised on Conventional, Antibiotic-Free, and Organic Farms

Poultry is the most likely source of livestock-associated Extended Spectrum Beta-Lactamase (ESBL) and plasmid-mediated AmpC (pAmpC)-producing E. coli (EC) for humans. We tested the hypothesis that farming methods have an impact on the load of ESBL/pAmpC-EC in the gut of broilers at slaughter. Isolates (n = 156) of antibiotic-free (AF), organic (O), and conventional (C) animals were characterized for antibiotic susceptibility and antibiotic resistance genes. Thirteen isolates were whole-genome sequenced. The average loads of ESBL/pAmpC-EC in cecal contents were 4.17 Log CFU/g for AF; 2.85 Log CFU/g for O; and 3.88 Log CFU/g for C type (p < 0.001). ESBL/pAmpC-EC isolates showed resistance to antibiotic classes historically used in poultry, including penicillins, tetracyclines, quinolones, and sulfonamides. Isolates from O and AF farms harbored a lower proportion of resistance to antibiotics than isolates from C farms. Among the determinants for ESBL/pAmpC, CTX-M-1 prevailed (42.7%), followed by TEM-type (29%) and SHV (19.8%). Avian pathogenic E. coli (APEC), belonging to ST117 and ST349, were identified in the collection. These data confirm the possible role of a broiler as an ESBL/AmpC EC and APEC reservoir for humans. Overall, our study suggests that antibiotic-free and organic production may contribute to a reduced exposure to ESBL/AmpC EC for the consumer.

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.

Whole Genome Multi-Locus Sequence Typing and Genomic Single Nucleotide Polymorphism Analysis for Epidemiological Typing of Pseudomonas aeruginosa From Indonesian Intensive Care Units

We have previously studied carbapenem non-susceptible Pseudomonas aeruginosa (CNPA) strains from intensive care units (ICUs) in a referral hospital in Jakarta, Indonesia (Pelegrin et al., 2019). We documented that CNPA transmissions and acquisitions among patients were variable over time and that these were not significantly reduced by a set of infection control measures. Three high risk international CNPA clones (sequence type (ST)235, ST823, ST357) dominated, and carbapenem resistance was due to carbapenemase-encoding genes and mutations in the porin OprD. Pelegrin et al. (2019) reported core genome analysis of these strains. We present a more refined and detailed whole genome-based analysis of major clones represented in the same dataset. As per our knowledge, this is the first study reporting Single Nucleotide Polymorphisms (wgSNP) analysis of Pseudomonas strains. With whole genome-based Multi Locus Sequence Typing (wgMLST) of the 3 CNPA clones (ST235, ST357 and ST823), three to eleven subgroups with up to 200 allelic variants were observed for each of the CNPA clones. Furthermore, we analyzed these CNPA clone clusters for the presence of wgSNP to redefine CNPA transmission events during hospitalization. A maximum number 35350 SNPs (including non-informative wgSNPs) and 398 SNPs (ST-specific_informative-wgSNPs) were found in ST235, 34,570 SNPs (including non-informative wgSNPs) and 111 SNPs (ST-specific_informative-wgSNPs) in ST357 and 26,443 SNPs (including non-informative SNPs) and 61 SNPs (ST-specific_informative-wgSNPs) in ST823. ST-specific_Informative-wgSNPs were commonly noticed in sensor-response regulator genes. However, the majority of non-informative wgSNPs was found in conserved hypothetical proteins or in uncharacterized proteins. Of note, antibiotic resistance and virulence genes segregated according to the wgSNP analyses. A total of 8 transmission chains for ST235 strains followed by 9 and 4 possible transmission chains for ST357 and ST823 were traceable on the basis of pairwise distances of informative-wgSNPs (0 to 4 SNPs) among the strains. The present study demonstrates the value of detailed whole genome sequence analysis for highly refined epidemiological analysis of P. aeruginosa.

The Neisseria gonorrhoeae Accessory Genome and Its Association with the Core Genome and Antimicrobial Resistance

The bacterial accessory genome provides the genetic flexibility needed to facilitate environment and host adaptation. In Neisseria gonorrhoeae, known accessory elements include plasmids which can transfer and mediate antimicrobial resistance (AMR); however, chromosomal accessory genes could also play a role in AMR. Here, the gonococcal accessory genome was characterized using gene-by-gene approaches and its association with the core genome and AMR were assessed. The gonococcal accessory gene pool consisted of 247 genes, which were mainly genes located on large mobile genetic elements, phage associated genes, or genes encoding putative secretion systems. Accessory elements showed similar synteny across genomes, indicating either a predisposition for particular genomic locations or ancestral inheritance that are conserved during strain expansion. Significant associations were found between the prevalence of accessory elements and core genome multi-locus sequence types (cgMLST), consistent with a structured gonococcal population despite frequent horizontal gene transfer (HGT). Increased prevalence of putative DNA exchange regulators was significantly associated with AMR, which included a putative secretion system, methyltransferases and a toxin-antitoxin system. Although frequent HGT results in high genetic diversity in the gonococcus, we found that this is mediated by a small gene pool. In fact, a highly organized genome composition was identified with a strong association between the accessory and core genome. Increased prevalence of DNA exchange regulators in antimicrobial resistant isolates suggests that genetic material exchange plays a role in the development or maintenance of AMR. These findings enhance our understanding of gonococcal genome architecture and have important implications for gonococcal population biology. IMPORTANCE The emergence of antimicrobial resistance (AMR) against third generation cephalosporins in Neisseria gonorrhoeae is a major public health concern, as these are antibiotics of last resort for the effective treatment of gonorrhea. Although the resistance mechanisms against this class of antibiotics have not been entirely resolved, resistance against other classes of antibiotics, such as tetracyclines, is known to be mediated through plasmids, which are known gonococcal extra-chromosomal accessory elements. A complete assessment of the chromosomal accessory genome content and its role in AMR has not yet been undertaken. Here, we comprehensively characterize the gonococcal accessory genome to better understand genome architecture as well as the evolution and mechanisms of AMR in this species.

Novel Multiplex PCR Method and Genome Sequence-Based Analog for High-Resolution Subclonal Assignment and Characterization of Escherichia coli Sequence Type 131 Isolates

Escherichia coli sequence type 131 (ST131) is a pandemic, multidrug-resistant extraintestinal pathogen. The multiple distinctive ST131 subclones differ for rfb and fliC alleles (O and H antigens), fimH allele (type-1 fimbriae adhesin), resistance phenotype and genotype, clinical correlates, and host predilection. Current PCR assays for detecting ST131 and its main subclones offer limited sub-ST characterization. Here we combined 22 novel and 14 published primers for a multiplex PCR assay to detect and extensively characterize ST131 isolates. The primers target mdh36, gyrB47, trpA72, sbmA, plsB, nupC, rmuC, kefC, ybbW, the O16 and O25b rfb variants, five fimH alleles (fimH22, fimH27, fimH30, fimH35, and fimH41), two fliC alleles (H4 and H5), a (subclone-specific) fluoroquinolone resistance-associated parC allele, and a (subclone-specific) prophage marker. The resulting amplicons resolve 15 molecular subsets within ST131, including 3 within clade A (H41 subclone), 5 within clade B (H22 subclone), and 7 within clade C (H30 subclone), which includes subclones C0 (H30S: 2 subsets), C1 and C1-M27 (H30R1: 2 subsets), and C2 (H30Rx: 3 subsets). Validation in three laboratories showed that this assay provides a rapid, accurate, and portable method for rapidly detecting and characterizing E. coli ST131 and its key subsets. Additionally, for users with whole genome sequencing (WGS) capability, we developed a command-line executable called ST131Typer, an in silico version of the extended multiplex PCR assay. Its accuracy was 87.8%, with most issues due to incomplete or fragmented input genome assemblies. These two novel assays should facilitate detailed ST131 subtyping using either endpoint PCR or WGS. IMPORTANCE These novel assays provide greater subclonal resolution and characterization of E. coli ST131 isolates than do the available comparable PCR assays, plus offer a novel sequence-based alternative to PCR. They may prove useful for molecular epidemiological studies, surveillance, and, potentially, clinical management.

Genomic Diversity of Hospital-Acquired Infections Revealed through Prospective Whole-Genome Sequencing-Based Surveillance

Healthcare-associated infections (HAIs) cause mortality, morbidity, and waste of health care resources. HAIs are also an important driver of antimicrobial resistance, which is increasing around the world. Beginning in November 2016, we instituted an initiative to detect outbreaks of HAIs using prospective whole-genome sequencing-based surveillance of bacterial pathogens collected from hospitalized patients. Here, we describe the diversity of bacteria sampled from hospitalized patients at a single center, as revealed through systematic analysis of bacterial isolate genomes. We sequenced the genomes of 3,004 bacterial isolates from hospitalized patients collected over a 25-month period. We identified bacteria belonging to 97 distinct species, which were distributed among 14 groups of related species. Within these groups, isolates could be distinguished from one another by both average nucleotide identity (ANI) and principal-component analysis of accessory genes (PCA-A). Core genome genetic distances and rates of evolution varied among species, which has practical implications for defining shared ancestry during outbreaks and for our broader understanding of the origins of bacterial strains and species. Finally, antimicrobial resistance genes and putative mobile genetic elements were frequently observed, and our systematic analysis revealed patterns of occurrence across the different species sampled from our hospital. Overall, this study shows how understanding the population structure of diverse pathogens circulating in a single health care setting can improve the discriminatory power of genomic epidemiology studies and can help define the processes leading to strain and species differentiation. IMPORTANCE Hospitalized patients are at increased risk of becoming infected with antibiotic-resistant organisms. We used whole-genome sequencing to survey and compare over 3,000 clinical bacterial isolates collected from hospitalized patients at a large medical center over a 2-year period. We identified nearly 100 different bacterial species, which we divided into 14 different groups of related species. When we examined how genetic relatedness differed between species, we found that different species were likely evolving at different rates within our hospital. This is significant because the identification of bacterial outbreaks in the hospital currently relies on genetic similarity cutoffs, which are often applied uniformly across organisms. Finally, we found that antibiotic resistance genes and mobile genetic elements were abundant and were shared among the bacterial isolates we sampled. Overall, this study provides an in-depth view of the genomic diversity and evolutionary processes of bacteria sampled from hospitalized patients, as well as genetic similarity estimates that can inform hospital outbreak detection and prevention efforts.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Population genomics of Escherichia coli in livestock-keeping households across a rapidly developing urban landscape

Quantitative evidence for the risk of zoonoses and the spread of antimicrobial resistance remains lacking. Here, as part of the UrbanZoo project, we sampled Escherichia coli from humans, livestock and peri-domestic wildlife in 99 households across Nairobi, Kenya, to investigate its distribution among host species in this rapidly developing urban landscape. We performed whole-genome sequencing of 1,338 E. coli isolates and found that the diversity and sharing patterns of E. coli were heavily structured by household and strongly shaped by host type. We also found evidence for inter-household and inter-host sharing and, importantly, between humans and animals, although this occurs much less frequently. Resistome similarity was differently distributed across host and household, consistent with being driven by shared exposure to antimicrobials. Our results indicate that a large, epidemiologically structured sampling framework combined with WGS is needed to uncover strain-sharing events among different host populations in complex environments and the major contributing pathways that could ultimately drive the emergence of zoonoses and the spread of antimicrobial resistance.

Phylogeography and phylogenomic analyses of E. coli isolates collected from humans and domesticated and wild animals across 99 households in Nairobi reveal strong intra-household, and lower but detectable inter-household and inter-host, strain-sharing patterns.

A role for ColV plasmids in the evolution of pathogenic Escherichia coli ST58

Escherichia coli ST58 has recently emerged as a globally disseminated uropathogen that often progresses to sepsis. Unlike most pandemic extra-intestinal pathogenic E. coli (ExPEC), which belong to pathogenic phylogroup B2, ST58 belongs to the environmental/commensal phylogroup B1. Here, we present a pan-genomic analysis of a global collection of 752 ST58 isolates from diverse sources. We identify a large ST58 sub-lineage characterized by near ubiquitous carriage of ColV plasmids, which carry genes encoding virulence factors, and by a distinct accessory genome including genes typical of the Yersiniabactin High Pathogenicity Island. This sub-lineage includes three-quarters of all ExPEC sequences in our study and has a broad host range, although poultry and porcine sources predominate. By contrast, strains isolated from cattle often lack ColV plasmids. Our data indicate that ColV plasmid acquisition contributed to the divergence of the major ST58 sub-lineage, and different sub-lineages inhabit poultry, swine and cattle.

Genomic epidemiology links Burkholderia pseudomallei from individual human cases to B. pseudomallei from targeted environmental sampling in Northern Australia

Each case of melioidosis results from a single event when a human is infected by the environmental bacterium Burkholderia pseudomallei. Darwin in tropical northern Australia has the highest incidences of melioidosis globally and the Darwin Prospective Melioidosis Study (DPMS) commenced in 1989, documenting all culture confirmed melioidosis cases. From 2000-2019 we sampled DPMS patient's environments for B. pseudomallei when a specific location was considered to have been where infection occurred. With the aim to use genomic epidemiology to understand B. pseudomallei transmission and infecting scenarios. Environmental sampling was performed at 98 DPMS patient sites, where we collected 975 environmental samples (742 soil; 233 water). Genotyping matched the clinical and epidemiologically linked environmental B. pseudomallei for 19 patients (19%), with the environmental isolates cultured from soil (n=11) or water (n=8) sources. B. pseudomallei isolates from patients and their local environments that matched on genotyping were whole genome sequenced (WGS). Of the 19 patients with a clinical-environmental genotype match, 17 pairs clustered on a Darwin core genome single-nucleotide polymorphism (SNP) phylogeny, later confirmed by single ST phylogenies and pairwise comparative genomics. When related back to patient clinical scenarios, the matched clinical and environmental B. pseudomallei pairs informed likely modes of infection: percutaneous inoculation, inhalation, and ingestion. Targeted environmental sampling for B. pseudomallei can inform infecting scenarios for melioidosis and dangerous occupational and recreational activities and identify hot spots of B. pseudomallei presence. However, WGS and careful genomics are required to avoid overcalling the relatedness between clinical and environmental isolates of B. pseudomallei.

The cnf1 gene is associated with an expanding Escherichia coli ST131 H30Rx/C2 subclade and confers a competitive advantage for gut colonization

Epidemiological projections point to acquisition of ever-expanding multidrug resistance (MDR) by Escherichia coli, a commensal of the digestive tract and a source of urinary tract pathogens. Bioinformatics analyses of a large collection of E. coli genomes from EnteroBase, enriched in clinical isolates of worldwide origins, suggest the Cytotoxic Necrotizing Factor 1 (CNF1)-toxin encoding gene, cnf1, is preferentially distributed in four common sequence types (ST) encompassing the pandemic E. coli MDR lineage ST131. This lineage is responsible for a majority of extraintestinal infections that escape first-line antibiotic treatment, with known enhanced capacities to colonize the gastrointestinal tract. Statistical projections based on this dataset point to a global expansion of cnf1-positive multidrug-resistant ST131 strains from subclade H30Rx/C2, accounting for a rising prevalence of cnf1-positive strains in ST131. Despite the absence of phylogeographical signals, cnf1-positive isolates segregated into clusters in the ST131-H30Rx/C2 phylogeny, sharing a similar profile of virulence factors and the same cnf1 allele. The suggested dominant expansion of cnf1-positive strains in ST131-H30Rx/C2 led us to uncover the competitive advantage conferred by cnf1 for gut colonization to the clinical strain EC131GY ST131-H30Rx/C2 versus cnf1-deleted isogenic strain. Complementation experiments showed that colon tissue invasion was compromised in the absence of deamidase activity on Rho GTPases by CNF1. Hence, gut colonization factor function of cnf1 was confirmed for another clinical strain ST131-H30Rx/C2. In addition, functional analysis of the cnf1-positive clinical strain EC131GY ST131-H30Rx/C2 and a cnf1-deleted isogenic strain showed no detectable impact of the CNF1 gene on bacterial fitness and inflammation during the acute phase of bladder monoinfection. Together these data argue for an absence of role of CNF1 in virulence during UTI, while enhancing gut colonization capacities of ST131-H30Rx/C2 and suggested expansion of cnf1-positive MDR isolates in subclade ST131-H30Rx/C2.

Genomic comparisons of Escherichia coli ST131 from Australia

Escherichia coli ST131 is a globally dispersed extraintestinal pathogenic E. coli lineage contributing significantly to hospital and community acquired urinary tract and bloodstream infections. Here we describe a detailed phylogenetic analysis of the whole genome sequences of 284 Australian ST131 E. coli isolates from diverse sources, including clinical, food and companion animals, wildlife and the environment. Our phylogeny and the results of single nucleotide polymorphism (SNP) analysis show the typical ST131 clade distribution with clades A, B and C clearly displayed, but no niche associations were observed. Indeed, interspecies relatedness was a feature of this study. Thirty-five isolates (29 of human and six of wild bird origin) from clade A (32 fimH41, 2 fimH89, 1 fimH141) were observed to differ by an average of 76 SNPs. Forty-five isolates from clade C1 from four sources formed a cluster with an average of 46 SNPs. Within this cluster, human sourced isolates differed by approximately 37 SNPs from isolates sourced from canines, approximately 50 SNPs from isolates from wild birds, and approximately 52 SNPs from isolates from wastewater. Many ST131 carried resistance genes to multiple antibiotic classes and while 41 (14 %) contained the complete class one integron-integrase intI1, 128 (45 %) isolates harboured a truncated intI1 (462-1014 bp), highlighting the ongoing evolution of this element. The module intI1-dfrA17-aadA5-qacEΔ1-sul1-ORF-chrA-padR-IS1600-mphR-mrx-mphA, conferring resistance to trimethoprim, aminoglycosides, quaternary ammonium compounds, sulphonamides, chromate and macrolides, was the most common structure. Most (73 %) Australian ST131 isolates carry at least one extended spectrum β-lactamase gene, typically bla_CTX-M-15 and bla_CTX-M-27. Notably, dual parC-1aAB and gyrA-1AB fluoroquinolone resistant mutations, a unique feature of clade C ST131 isolates, were identified in some clade A isolates. The results of this study indicate that the the ST131 population in Australia carries diverse antimicrobial resistance genes and plasmid replicons and indicate cross-species movement of ST131 strains across diverse reservoirs.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Combining pangenome analysis to identify potential cross-protective antigens against avian pathogenic Escherichia coli

AbstractAvian pathogenic Escherichia coli (APEC) is the bacterial pathogen of poultry colibacillosis, which causes significant economic losses to the poultry industry. The lack of an effective vaccine against multiple serotypes and the emergence of multi-resistant isolates have made the control of avian colibacillosis troublesome. To identify conserved potential vaccine candidates, 58 genomes of APEC were obtained (54 sequenced by our laboratory and 4 downloaded from NCBI). A reverse vaccinology (RV) method based on the pangenome - called Pan-RV analysis - was performed in APEC-protective protein mining for the first time. Finally, four proteins were selected, and their immunoreactivity with anti-O1, O2, and O78 serum was verified by western blotting. Our in silico method of analysis will pave the way for rapid screening of vaccine candidates and will lay the foundation for the development of a highly effective subunit vaccine controlling APEC infection.

Ecological and evolutionary solutions to the plasmid paradox

The 'plasmid paradox' arises because, although plasmids are common features of bacterial genomes, theoretically they should not exist: rates of conjugation were believed insufficient to allow plasmids to persist by infectious transmission, whereas the costs of plasmid maintenance meant that plasmids should be purged by negative selection regardless of whether they encoded beneficial accessory traits because these traits should eventually be captured by the chromosome, enabling the loss of the redundant plasmid. In the decade since the plasmid paradox was described, new data and theory show that a range of ecological and evolutionary mechanisms operate in bacterial populations and communities to explain the widespread distribution and stable maintenance of plasmids. We conclude, therefore, that multiple solutions to the plasmid paradox are now well understood. The current challenge for the field, however, is to better understand how these solutions operate in natural bacterial communities to explain and predict the distribution of plasmids and the dynamics of the horizontal gene transfer that they mediate in bacterial (pan)genomes.

Phylogenetic Analysis Reveals Distinct Evolutionary Trajectories of the Fluoroquinolones-Resistant Escherichia coli ST1193 From Fuzhou, China

Escherichia coli (E. coli) ST1193 is an emerging fluoroquinolones-resistant and virulent lineage. Large gaps remain in our understanding of the evolutionary processes and differences of this lineage. Therefore, we used 76 E. coli ST1193 genomes to detect strain-level genetic diversity and phylogeny of this lineage globally. All E. coli ST1193 possessed fimH64, filCH5, and fumC14. There was 94.7% of isolates classified as O-type O75. There was 9.33% of E. coli ST1193 that possessed K5 capsular, while 90.67% of isolates possessed K1 capsular. The core genome analysis revealed that all isolates were divided into two phylogenetic clades (clade A and B). Clade A included 25 non-Chinese E. coli ST1193, and clade B contained all isolates collected from Fuzhou, China, respectively. The results of comparative genomics indicated Indels were identified in 150 clade-specific genes, which were enriched into the biological process and molecular function. Accessory genome phylogenetic tree showed a high degree of correlation between accessory genome clusters and core genome clades. There was significant difference in antibiotic resistance genes (ARGs) [bla _CTX-M-55 , bla _TEM-1 , sul2, tet(B), tet(R), APH(6)-Id, and AAC(3)-IId], virulence factors (cia, neuC, gad, and traT), and plasmid replicon types (IncQ1, Col156, and IncB/O/K/Z) between clade A (non-Chinese isolates) and clade B (Chinese isolates) (p < 0.05). Further analysis of the genetic environments of bla _CTX-M-55 demonstrated that the flanking contexts of bla _CTX-M-55 were diverse. In conclusion, our results reveal the distinct evolutionary trajectories of the spread of E. coli ST1193 in Fuzhou, China and non-China regions. This supports both global transmission and localized lineage expansion of this lineage following specific introductions into a geographic locality.

The role of potentiating mutations in the evolution of pandemic Escherichia coli clones

The Escherichia coli species exhibits a vast array of variable lifestyles, including environmental, commensal, and pathogenic organisms. Many of these E. coli contribute significantly to the global threat of antimicrobial resistance (AMR). Multidrug-resistant (MDR) clones of E. coli have arisen multiple times over varying timescales. The repeated emergence of successful pandemic clones, including the notorious ST131 lineage, highlights a desperate need to further study the evolutionary processes underlying their emergence and success. Here, we review the evolutionary emergence of E. coli ST131 pandemic clones and draw parallels between their evolutionary trajectories and those of other lineages. From colonization and expansion to the acquisition of multidrug resistance plasmids, potentiating mutations are present at each stage, leading to a proposed sequence of events that may result in the formation of an antimicrobial-resistant pandemic clone.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified (Butyrivibrio fibrisolvens, Butyrivibrio hungatei, Butyrivibrio proteoclasticus, Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio. Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations (P. ruminis, P. xylanivorans, B. fibrisolvens, Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens, B. hungatei and P. ruminis, based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Comparative activity of plazomicin against extended-spectrum cephalosporin-resistant Escherichia coli clinical isolates (2012-2017) in relation to phylogenetic background, sequence type 131 subclones, blaCTX-M genotype, and resistance to comparator agents

Extended-spectrum cephalosporin-resistant Escherichia coli (ESCREC) are a growing threat. Leading ESCREC lineages include sequence type ST131, especially its (blaCTX-M-15-associated) H30Rx subclone and (blaCTX-M-27-associated) C1-M27 subset within the H30R1 subclone. The comparative activity against such strains of alternative antimicrobial agents, including the recently developed aminoglycoside plazomicin, is undefined, so was investigated here. We assessed plazomicin and 11 comparators for activity against 216 well-characterized ESCREC isolates (Minnesota, 2012-2017) and then compared broth microdilution MICs with phylogenetic and clonal background, beta-lactamase genotype (blaCTX-M; group 1 and 9 variants), and co-resistance. Percent susceptible was > 99% for plazomicin, meropenem, imipenem, and tigecycline; 96-98% for amikacin and ertapenem; and ≤ 75% for the remaining comparators. For most comparators, MICs varied significantly in relation to multiple bacterial characteristics, in agent-specific patterns. By contrast, for plazomicin, the only bacterial characteristic significantly associated with MICs was ST131 subclone: plazomicin MICs were lowest among O16 ST131 isolates and highest among ST131-H30R1 C1-M27 subclone isolates. Additionally, plazomicin MICs varied significantly in relation to resistance vs. susceptibility to comparator agents only for amikacin and levofloxacin. For most study agents, antimicrobial activity against ESCREC varied extensively in relation to multiple bacterial characteristics, including clonal background, whereas for plazomicin, it varied only by ST131 subclone (C1-M27 isolates least susceptible, O16 isolates most susceptible). These findings support plazomicin as a reliable alternative for treating ESCREC infections and urge continued attention to the C1-M27 ST131 subclone.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.