Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae

Spontaneous fermentation mitigates the frequency of genes encoding antimicrobial resistance spreading from the phyllosphere reservoir to the diet

The phyllosphere microbiome of vegetable products constitutes an important reservoir for multidrug resistant bacteria and Antibiotic Resistance Genes (ARG). Vegetable products including fermented products such as Paocai therefore may serve as a shuttle for extrinsic microorganisms with ARGs into the gut of consumers. Here we study the effect of fermentation on Paocai ARG dissemination by metagenomic analysis. Microbial abundance and diversity of the Paocai microbiome were diminished during fermentation, which correlated with the reduction of abundance in ARGs. Specifically, as fermentation progressed, Enterobacterales overtook Pseudomonadales as the predominant ARG carriers, and Lactobacillales and Enterobacteriales became the determinants of Paocai resistome variation. Moreover, the dual effect of microbes and metal resistance genes (MRGs) was the major contributor driving Paocai resistome dynamics. We recovered several metagenome-assembled genomes (MAGs) carrying acquired ARGs in the phyllosphere microbiome. ARGs of potential clinical and epidemiological relevance such as tet M and emrB-qacA, were mainly hosted by non-dominant bacterial genera. Overall, our study provides evidence that changes in microbial community composition by fermentation aid in constraining ARG dispersal from raw ingredients to the human microbiome but does not eliminate them.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Occurrence and characterization of rmtB-harbouring Salmonella and Escherichia coli isolates from a pig farm in the UK

OBJECTIVES

To characterize and elucidate the spread of amikacin-resistant Enterobacteriaceae isolates from environmental samples on a pig farm in the UK, following the previous identification of index Salmonella isolates harbouring the rmtB gene, a 16S rRNA methylase.

METHODS

Environmental samples were collected during two visits to a pig farm in the UK. Isolates were recovered using selective media (amikacin 128 mg/L) followed by real-time PCR and WGS to analyse rmtB-carrying Salmonella and Escherichia coli isolates.

RESULTS

Salmonella and E. coli isolates harbouring the rmtB gene were detected at both farm visits. All Salmonella isolates were found to be monophasic S. enterica serovar Typhimurium variant Copenhagen of ST34. rmtB-harbouring E. coli isolates were found to be one of three STs: ST4089, ST1684 and ST34. Long-read sequencing identified the rmtB gene to be chromosomally located in Salmonella isolates and on IncFII-type plasmids in E. coli isolates. The results showed the rmtB gene to be flanked by IS26 elements and several resistance genes.

CONCLUSIONS

We report on the occurrence of rmtB-harbouring Enterobacteriaceae on a pig farm in the UK. rmtB confers resistance to multiple aminoglycosides and this work highlights the need for surveillance to assess dissemination and risk.

Exploring gene content with pangene graphs

Motivation

The gene content regulates the biology of an organism. It varies between species and between individuals of the same species. Although tools have been developed to identify gene content changes in bacterial genomes, none is applicable to collections of large eukaryotic genomes such as the human pangenome.

Results

We developed pangene, a computational tool to identify gene orientation, gene order and gene copy-number changes in a collection of genomes. Pangene aligns a set of input protein sequences to the genomes, resolves redundancies between protein sequences and constructs a gene graph with each genome represented as a walk in the graph. It additionally finds subgraphs, which we call bibubbles, that capture gene content changes. Applied to the human pangenome, pangene identifies known gene-level variations and reveals complex haplotypes that are not well studied before. Pangene also works with high-quality bacterial pangenome and reports similar numbers of core and accessory genes in comparison to existing tools.

Availability and implementation

Source code at https://github.com/lh3/pangene; pre-built pangene graphs can be downloaded from https://zenodo.org/records/8118576 and visualized at https://pangene.bioinweb.org.

Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

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.

Analysis of the transmission chain of carbapenem-resistant Enterobacter cloacae complex infections in clinical, intestinal and healthcare settings in Zhejiang province, China (2022-2023)

Considerable attention is given to intensive care unit-acquired infections; however, research on the transmission dynamics of multichain carbapenemase-resistant Enterobacter cloacae complex (CRECC) outbreaks remains elusive. A total of 118 non-duplicated CRECC strains were isolated from the clinical, intestinal, and hospital sewage samples collected from Zhejiang province of China during 2022-2023. A total of 64 CRECC strains were isolated from the hospital sewage samples, and their prevalence increased from 10.0 % (95 % confidence interval, CI = 0.52-45.8 %) in 2022 to 63.6 % (95 % CI = 31.6-87.6 %) in 2023. Species-specific identification revealed that Enterobacter hormaechei was the predominant CRECC species isolated in this study (53.4 %, 95 % CI = 44.0-62.6 %). The antimicrobial susceptibility profiles indicated that all 118 CRECC strains conferred high-level resistance to β-lactam antibiotics, ceftacillin/avibactam, and polymyxin. Furthermore, all CRECC strains exhibited resistance to β-lactams, quinolones, and fosfomycin, with a higher colistin resistance rate observed in the hospital sewage samples (67.2 %, 95 % CI = 54.2-78.1 %). Several antibiotic resistance genes were identified in CRECC strains, including Class A carbapenemases (blaKPC-2) and Class B carbapenemases (blaNDM-1/blaIMP), but not Class D carbapenemases. The WGS analysis showed that the majority of the CRECC strains carried carbapenemase-encoding genes, with blaNDM-1 being the most prevalent (86.9 %, 95 % CI = 77.4-92.9 %). Furthermore, sequence typing revealed that the isolated CRECC strains belonged to diverse sequence types (STs), among which ST418 was the most prevalent blaNDM-positive strain. The high risk of carbapenemase-producing ST418 E. hormaechei and the blaNDM-harboring IncFIB-type plasmid (81.4 %, 95 % CI = 72.9-87.7 %) were detected and emphasized in this study. This study provides valuable insights into the prevalence, antimicrobial resistance, genomic characteristics, and plasmid analysis of CRECC strains in diverse populations and environments. The clonal relatedness analysis showed sporadic clonal transmission of ST418 E. hormaechei strains, supporting inter-hospital transmission.

Genomic revisitation and reclassification of the genus Providencia

Members of Providencia, although typically opportunistic, can cause severe infections in immunocompromised hosts. Recent advances in genome sequencing provide an opportunity for more precise study of this genus. In this study, we first identified and characterized a novel species named Providencia zhijiangensis sp. nov. It has ≤88.23% average nucleotide identity (ANI) and ≤31.8% in silico DNA-DNA hybridization (dDDH) values with all known Providencia species, which fall significantly below the species-defining thresholds. Interestingly, we found that Providencia stuartii and Providencia thailandensis actually fall under the same species, evidenced by an ANI of 98.59% and a dDDH value of 90.4%. By fusing ANI with phylogeny, we have reclassified 545 genomes within this genus into 20 species, including seven unnamed taxa (provisionally titled Taxon 1-7), which can be further subdivided into 23 lineages. Pangenomic analysis identified 1,550 genus-core genes in Providencia, with coenzymes being the predominant category at 10.56%, suggesting significant intermediate metabolism activity. Resistance analysis revealed that most lineages of the genus (82.61%, 19/23) carry a high number of antibiotic-resistance genes (ARGs) and display diverse resistance profiles. Notably, the majority of ARGs are located on plasmids, underscoring the significant role of plasmids in the resistance evolution within this genus. Three species or lineages (P. stuartii, Taxon 3, and Providencia hangzhouensis L12) that possess the highest number of carbapenem-resistance genes suggest their potential influence on clinical treatment. These findings underscore the need for continued surveillance and study of this genus, particularly due to their role in harboring antibiotic-resistance genes.

IMPORTANCE

The Providencia genus, known to harbor opportunistic pathogens, has been a subject of interest due to its potential to cause severe infections, particularly in vulnerable individuals. Our research offers groundbreaking insights into this genus, unveiling a novel species, Providencia zhijiangensis sp. nov., and highlighting the need for a re-evaluation of existing classifications. Our comprehensive genomic assessment offers a detailed classification of 545 genomes into distinct species and lineages, revealing the rich biodiversity and intricate species diversity within the genus. The substantial presence of antibiotic-resistance genes in the Providencia genus underscores potential challenges for public health and clinical treatments. Our study highlights the pressing need for increased surveillance and research, enriching our understanding of antibiotic resistance in this realm.

The plasmidome associated with Gram-negative bloodstream infections: A large-scale observational study using complete plasmid assemblies

Plasmids carry genes conferring antimicrobial resistance and other clinically important traits, and contribute to the rapid dissemination of such genes. Previous studies using complete plasmid assemblies, which are essential for reliable inference, have been small and/or limited to plasmids carrying antimicrobial resistance genes (ARGs). In this study, we sequenced 1,880 complete plasmids from 738 isolates from bloodstream infections in Oxfordshire, UK. The bacteria had been originally isolated in 2009 (194 isolates) and 2018 (368 isolates), plus a stratified selection from intervening years (176 isolates). We demonstrate that plasmids are largely, but not entirely, constrained to a single host species, although there is substantial overlap between species of plasmid gene-repertoire. Most ARGs are carried by a relatively small number of plasmid groups with biological features that are predictable. Plasmids carrying ARGs (including those encoding carbapenemases) share a putative 'backbone' of core genes with those carrying no such genes. These findings suggest that future surveillance should, in addition to tracking plasmids currently associated with clinically important genes, focus on identifying and monitoring the dissemination of high-risk plasmid groups with the potential to rapidly acquire and disseminate these genes.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Convergence of resistance and evolutionary responses in Escherichia coli and Salmonella enterica co-inhabiting chicken farms in China

Sharing of genetic elements among different pathogens and commensals inhabiting same hosts and environments has significant implications for antimicrobial resistance (AMR), especially in settings with high antimicrobial exposure. We analysed 661 Escherichia coli and Salmonella enterica isolates collected within and across hosts and environments, in 10 Chinese chicken farms over 2.5 years using data-mining methods. Most isolates within same hosts possessed the same clinically relevant AMR-carrying mobile genetic elements (plasmids: 70.6%, transposons: 78%), which also showed recent common evolution. Supervised machine learning classifiers revealed known and novel AMR-associated mutations and genes underlying resistance to 28 antimicrobials, primarily associated with resistance in E. coli and susceptibility in S. enterica. Many were essential and affected same metabolic processes in both species, albeit with varying degrees of phylogenetic penetration. Multi-modal strategies are crucial to investigate the interplay of mobilome, resistance and metabolism in cohabiting bacteria, especially in ecological settings where community-driven resistance selection occurs.

Resistome Diversity in Escherichia coli Isolates of Global Wastewaters

Antimicrobial resistance (AMR) is a global health threat requiring urgent attention and effective strategies for containment. AMR is fueled by wastewater mismanagement and global mobility, disseminating multidrug-resistant (MDR) strains worldwide. While global estimates of AMR burden have been informative, community-level understanding has received little attention despite reports of high AMR prevalence in healthy communities. We assessed the "invasion" of antibiotic resistance genes (ARGs) into the normal human flora by characterizing AMR Escherichia coli in local wastewaters contributed by a healthy youth population. This study estimated 26% (out of 300 isolates) resistant and 59% plasmid-bearing E. coli in local wastewater. Of the 78 AMR isolates, the frequency of mono-resistance was higher against tetracycline (32%), followed by kanamycin (17%) and chloramphenicol (9%). Five isolates were potentially MDR. We further sequenced four MDRs and four sensitive strains to comprehend the genome and resistome diversity in comparison to the global wastewater E. coli (genomes from the PATRIC database). The whole-genome analysis revealed extensive genome similarity among global isolates, suggesting global dissemination and colonization of E. coli. Global wastewater resistome majorly comprised ARGs against aminoglycosides (26%), beta-lactam (17%), sulfonamide (11%), and trimethoprim (8%). Resistance to colistin, a last-resort antibiotic, was prevalent in MDRs of European and South Asian isolates. A systems approach is required to address the AMR crisis on a global scale, reduce antibiotic usage, and increase the efficiency of wastewater management and disinfection.

Rapid Detection of Antimicrobial Resistance Genes in Critically Ill Children Using a Custom TaqMan Array Card

Bacteria are identified in only 22% of critically ill children with respiratory infections treated with antimicrobial therapy. Once an organism is isolated, antimicrobial susceptibility results (phenotypic testing) can take another day. A rapid diagnostic test identifying antimicrobial resistance (AMR) genes could help clinicians make earlier, informed antimicrobial decisions. Here we aimed to validate a custom AMR gene TaqMan Array Card (AMR-TAC) for the first time and assess its feasibility as a screening tool in critically ill children. An AMR-TAC was developed using a combination of commercial and bespoke targets capable of detecting 23 AMR genes. This was validated using isolates with known phenotypic resistance. The card was then tested on lower respiratory tract and faecal samples obtained from mechanically ventilated children in a single-centre observational study of respiratory infection. There were 82 children with samples available, with a median age of 1.2 years. Major comorbidity was present in 29 (35%) children. A bacterial respiratory pathogen was identified in 13/82 (16%) of children, of which 4/13 (31%) had phenotypic AMR. One AMR gene was detected in 49/82 (60%), and multiple AMR genes were detected in 14/82 (17%) children. Most AMR gene detections were not associated with the identification of phenotypic AMR. AMR genes are commonly detected in samples collected from mechanically ventilated children with suspected respiratory infections. AMR-TAC may have a role as an adjunct test in selected children in whom there is a high suspicion of antimicrobial treatment failure.

Distinct increase in antimicrobial resistance genes among Vibrio parahaemolyticus in recent decades worldwide

Vibrio parahaemolyticus is a common pathogen, and has emerged with multiple antimicrobial resistance (AMR). However, few studies have conducted large-scale investigations of AMR and virulence trends of V. parahaemolyticus worldwide. This study longitudinally monitored antibiotic resistance genes (ARGs) and virulence factor genes (VFGs) trends of 1540 V. parahaemolyticus isolates isolated from 1951 to 2021. The number of ARGs in V. parahaemolyticus isolates distinctly increased over the years (P = 5.9e-10), while the number of VFGs decreased significantly (P < 2.2e-16). However, the number of VFGs of isolates isolated from humans has not changed significantly over the years (R = 0.013, P = 0.74), suggesting that the pathogenic risk to humans has not been reduced. Besides, mobile genetic elements are important contributors to ARGs in V. parahaemolyticus (R = 0.34, P < 2.2e-16), but have no promoting effect on VFGs (P = 0.50). The structural equation model illustrated that the human development index promoted the consumption of antibiotics, thereby indirectly promoting an increase in the AMR of the V. parahaemolyticus isolates. Finally, the random forest was performed to predict the ARG and VFG risks of global terrestrial V. parahaemolyticus isolates, and successfully map these threats with over 80% accuracy. This study aimed to evaluate the global risks posed by AMR and virulence, which helps to develop methods specifically targeting V. parahaemolyticus to mitigate these threats.

Systematic analysis of plasmids of the Serratia marcescens complex using 142 closed genomes

Plasmids play important roles in bacterial genome diversification. In the Serratia marcescens complex (SMC), a notable contribution of plasmids to genome diversification was also suggested by our recent analysis of >600 draft genomes. As accurate analyses of plasmids in draft genomes are difficult, in this study we analysed 142 closed genomes covering the entire complex, 67 of which were obtained in this study, and identified 132 plasmids (1.9-244.4 kb in length) in 77 strains. While the average numbers of plasmids in clinical and non-clinical strains showed no significant difference, strains belonging to clade 2 (one of the two hospital-adapted lineages) contained more plasmids than the others. Pangenome analysis revealed that of the 28 954 genes identified, 12.8 % were plasmid-specific, and 1.4 % were present in plasmids or chromosomes depending on the strain. In the latter group, while transposon-related genes were most prevalent (31.4 % of the function-predicted genes), genes related to antimicrobial resistance and heavy metal resistance accounted for a notable proportion (22.7 %). Mash distance-based clustering separated the 132 plasmids into 23 clusters and 50 singletons. Most clusters/singletons showed notably different GC contents compared to those of host chromosomes, suggesting their recent or relatively recent appearance in the SMC. Among the 23 clusters, 17 were found in only clinical or only non-clinical strains, suggesting the possible preference of their distribution on the environmental niches of host strains. Regarding the host strain phylogeny, 16 clusters were distributed in two or more clades, suggesting their interclade transmission. Moreover, for many plasmids, highly homologous plasmids were found in other species, indicating the broadness of their potential host ranges, beyond the genus, family, order, class or even phylum level. Importantly, highly homologous plasmids were most frequently found in Klebsiella pneumoniae and other species in the family Enterobacteriaceae, suggesting that this family, particularly K. pneumoniae, is the main source for plasmid exchanges with the SMC. These results highlight the power of closed genome-based analysis in the investigation of plasmids and provide important insights into the nature of plasmids distributed in the SMC.

plASgraph2: using graph neural networks to detect plasmid contigs from an assembly graph

Identification of plasmids from sequencing data is an important and challenging problem related to antimicrobial resistance spread and other One-Health issues. We provide a new architecture for identifying plasmid contigs in fragmented genome assemblies built from short-read data. We employ graph neural networks (GNNs) and the assembly graph to propagate the information from nearby nodes, which leads to more accurate classification, especially for short contigs that are difficult to classify based on sequence features or database searches alone. We trained plASgraph2 on a data set of samples from the ESKAPEE group of pathogens. plASgraph2 either outperforms or performs on par with a wide range of state-of-the-art methods on testing sets of independent ESKAPEE samples and samples from related pathogens. On one hand, our study provides a new accurate and easy to use tool for contig classification in bacterial isolates; on the other hand, it serves as a proof-of-concept for the use of GNNs in genomics. Our software is available at https://github.com/cchauve/plasgraph2 and the training and testing data sets are available at https://github.com/fmfi-compbio/plasgraph2-datasets.

Uncharacterized and lineage-specific accessory genes within the Proteus mirabilis pan-genome landscape

Proteus mirabilis is a Gram-negative bacterium recognized for its unique swarming motility and urease activity. A previous proteomic report on four strains hypothesized that, unlike other Gram-negative bacteria, P. mirabilis may not exhibit significant intraspecies variation in gene content. However, there has not been a comprehensive analysis of large numbers of P. mirabilis genomes from various sources to support or refute this hypothesis. We performed comparative genomic analysis on 2,060 Proteus genomes. We sequenced the genomes of 893 isolates recovered from clinical specimens from three large US academic medical centers, combined with 1,006 genomes from NCBI Assembly and 161 genomes assembled from Illumina reads in the public domain. We used average nucleotide identity (ANI) to delineate species and subspecies, core genome phylogenetic analysis to identify clusters of highly related P. mirabilis genomes, and pan-genome annotation to identify genes of interest not present in the model P. mirabilis strain HI4320. Within our cohort, Proteus is composed of 10 named species and 5 uncharacterized genomospecies. P. mirabilis can be subdivided into three subspecies; subspecies 1 represented 96.7% (1,822/1,883) of all genomes. The P. mirabilis pan-genome includes 15,399 genes outside of HI4320, and 34.3% (5,282/15,399) of these genes have no putative assigned function. Subspecies 1 is composed of several highly related clonal groups. Prophages and gene clusters encoding putatively extracellular-facing proteins are associated with clonal groups. Uncharacterized genes not present in the model strain P. mirabilis HI4320 but with homology to known virulence-associated operons can be identified within the pan-genome. IMPORTANCE Gram-negative bacteria use a variety of extracellular facing factors to interact with eukaryotic hosts. Due to intraspecies genetic variability, these factors may not be present in the model strain for a given organism, potentially providing incomplete understanding of host-microbial interactions. In contrast to previous reports on P. mirabilis, but similar to other Gram-negative bacteria, P. mirabilis has a mosaic genome with a linkage between phylogenetic position and accessory genome content. P. mirabilis encodes a variety of genes that may impact host-microbe dynamics beyond what is represented in the model strain HI4320. The diverse, whole-genome characterized strain bank from this work can be used in conjunction with reverse genetic and infection models to better understand the impact of accessory genome content on bacterial physiology and pathogenesis of infection.

Antimicrobial Peptides in the Global Microbiome: Biosynthetic Genes and Resistance Determinants

Antimicrobial peptides are a promising new class of antimicrobials that could address the antibiotic resistance crisis, which poses a major threat to human health. These peptides are present in all kingdoms of life, but especially in microorganisms, having multiple origins in diverse taxa. To date, there has been no global study on the diversity of antimicrobial peptides, the hosts in which these occur, and the potential for resistance to these agents. Here, we investigated the diversity and number of antimicrobial peptides in four main habitats (aquatic, terrestrial, human, and engineered) by analyzing 52,515 metagenome-assembled genomes. The number of antimicrobial peptides was higher in the human gut microbiome than in other habitats, and most hosts of antimicrobial peptides were habitat-specific. The relative abundance of genes that confer resistance to antimicrobial peptides varied between habitats and was generally low, except for the built environment and on human skin. The horizontal transfer of potential resistance genes among these habitats was probably constrained by ecological barriers. We systematically quantified the risk of each resistance determinant to human health and found that nearly half of them pose a threat, especially those that confer resistance to multiple AMPs and polymyxin B. Our results help identify the biosynthetic potential of antimicrobial peptides in the global microbiome, further identifying peptides with a low risk of developing resistance.

A longitudinal study reveals persistence of antimicrobial resistance on livestock farms is not due to antimicrobial usage alone

Introduction

There are concerns that antimicrobial usage (AMU) is driving an increase in multi-drug resistant (MDR) bacteria so treatment of microbial infections is becoming harder in humans and animals. The aim of this study was to evaluate factors, including usage, that affect antimicrobial resistance (AMR) on farm over time.

Methods

A population of 14 cattle, sheep and pig farms within a defined area of England were sampled three times over a year to collect data on AMR in faecal Enterobacterales flora; AMU; and husbandry or management practices. Ten pooled samples were collected at each visit, with each comprising of 10 pinches of fresh faeces. Up to 14 isolates per visit were whole genome sequenced to determine presence of AMR genes.

Results

Sheep farms had very low AMU in comparison to the other species and very few sheep isolates were genotypically resistant at any time point. AMR genes were detected persistently across pig farms at all visits, even on farms with low AMU, whereas AMR bacteria was consistently lower on cattle farms than pigs, even for those with comparably high AMU. MDR bacteria was also more commonly detected on pig farms than any other livestock species.

Discussion

The results may be explained by a complex combination of factors on pig farms including historic AMU; co-selection of AMR bacteria; variation in amounts of antimicrobials used between visits; potential persistence in environmental reservoirs of AMR bacteria; or importation of pigs with AMR microbiota from supplying farms. Pig farms may also be at increased risk of AMR due to the greater use of oral routes of group antimicrobial treatment, which were less targeted than cattle treatments; the latter mostly administered to individual animals. Also, farms which exhibited either increasing or decreasing trends of AMR across the study did not have corresponding trends in their AMU. Therefore, our results suggest that factors other than AMU on individual farms are important for persistence of AMR bacteria on farms, which may be operating at the farm and livestock species level.

Aeromonas spp. from hospital sewage act as a reservoir of genes resistant to last-line antibiotics

BACKGROUND

Aeromonas species are opportunistic pathogens distributed widely in the ecosystem. They are known to be capable of acquiring antibiotic resistance genes, including those encoding proteins against last-line antibiotics, such as the tmexCD-toprJ, mcr and carbapenemase genes. We investigated the genomic and phenotypic characteristics of tmexCD-toprJ-positive Aeromonas strains collected from human, animals, and water samples, particularly those from hospital wastewater in China.

METHODS

Samples were collected from living animals, meat, water and human. Aeromonas strains in these samples were isolated in selective media. Antimicrobial resistance profiles of all Aeromonas strains were tested by the broth microdilution method. The presence of tmexCD-toprJ was verified by polymerase chain reaction (PCR). All tmexCD-toprJ-positive (n = 36) and selected tmexCD-toprJ-negative (n = 18) Aeromonas strains were subjected to whole genome sequencing. Carriage of antimicrobial resistance genes, the genetic environment of tmexCD-toprJ and genetic diversity of tmexCD-toprJ-positive Aeromonas strains were determined by bioinformatics analysis. Phylogenetic tree of the Aeromonas strains was built by using the Harvest Suite.

FINDINGS

Among the 636 Aeromonas strains isolated from different sources, 36 were positive for tmexCD-toprJ, with the highest prevalence of tmexCD-toprJ being found in fishes (8.8%, 95 CI% 3.6-17.2%), followed by hospital wastewater (6.5%, 95 CI% 4.3-9.3%), river water (2.0%, 0.1-10.9) and duck (1.2%, 95 CI% 3.6-17.2%). All tmexCD-toprJ-positive Aeromonas strains carried multiple antimicrobial resistance genes and exhibited resistance to different classes of antibiotics. Co-existence of tmexCD-toprJ, mcr and bla_KPC-2 were identified in 21 strains. The tmexCD-toprJ-positive Aeromonas strains were genetically diverse and found to belong to four different species that could be clustered into three major lineages. The tmexCD-toprJ gene clusters were predominantly located in the chromosome (35/36) of Aeromonas spp., with only one strain carrying the plasmid-borne tmexCD-toprJ cluster. The tmexCD-toprJ genes were associated with seven different types of genetic environments, each of which carried distinct types of mobile elements that may be responsible for mediating transmission of this gene cluster.

How do interactions between mobile genetic elements affect horizontal gene transfer?

Horizontal gene transfer is central to bacterial adaptation and is facilitated by mobile genetic elements (MGEs). Increasingly, MGEs are being studied as agents with their own interests and adaptations, and the interactions MGEs have with one another are recognised as having a powerful effect on the flow of traits between microbes. Collaborations and conflicts between MGEs are nuanced and can both promote and inhibit the acquisition of new genetic material, shaping the maintenance of newly acquired genes and the dissemination of important adaptive traits through microbiomes. We review recent studies that shed light on this dynamic and oftentimes interlaced interplay, highlighting the importance of genome defence systems in mediating MGE-MGE conflicts, and outlining the consequences for evolutionary change, that resonate from the molecular to microbiome and ecosystem levels.

Genetic Determinants of Escherichia coli Survival in Beach Sand

Escherichia coli contain a high level of genetic diversity and are generally associated with the guts of warm-blooded animals but have also been isolated from secondary habitats outside hosts. We used E. coli isolates from previous in situ microcosm experiments conducted under actual beach conditions and performed population-level genomic analysis to identify accessory genes associated with survival within the beach sand environment. E. coli strains capable of surviving had been selected for by seeding isolates originating from sand, sewage, and gull waste (n = 528; 176 from each source) into sand, which was sealed in microcosm chambers and buried for 45 days in the backshore beach of Lake Michigan. In the current work, survival-associated genes were identified by comparing the pangenome of viable E. coli populations at the end of the microcosm experiment with the original isolate collection and identifying loci enriched in the out put samples. We found that environmental survival was associated with a wide variety of genetic factors, with the majority corresponding to metabolism enzymes and transport proteins. Of the 414 unique functions identified, most were present across E. coli phylogroups, except B2 which is often associated with human pathogens. Gene modules that were enriched in surviving populations included a betaine biosynthesis pathway, which produces an osmoprotectant, and the GABA (gamma-aminobutyrate) biosynthesis pathway, which aids in pH homeostasis and nutrient use versatility. Overall, these results demonstrate that the genetic flexibility within this species allows for survival in the environment for extended periods. IMPORTANCE Escherichia coli is commonly used as an indicator of recent fecal pollution in recreational water despite its known ability to survive in secondary environments, such as beach sand. These long-term survivors from sand reservoirs can be introduced into the water column through wave action or runoff during precipitation events, thereby impacting the perception of local water quality. Current beach monitoring methods cannot differentiate long-term environmental survivors from E. coli derived from recent fecal input, resulting in inaccurate monitoring results and unnecessary beach closures. This work identified the genetic factors that are associated with long-term survivors, providing insight into the mechanistic basis for E. coli accumulation in beach sand. A greater understanding of the intrinsic ability of E. coli to survive long-term and conditions that promote such survival will provide evidence of the limitations of beach water quality assessments using this indicator.

Comparison of R9.4.1/Kit10 and R10/Kit12 Oxford Nanopore flowcells and chemistries in bacterial genome reconstruction

Complete, accurate, cost-effective, and high-throughput reconstruction of bacterial genomes for large-scale genomic epidemiological studies is currently only possible with hybrid assembly, combining long- (typically using nanopore sequencing) and short-read (Illumina) datasets. Being able to use nanopore-only data would be a significant advance. Oxford Nanopore Technologies (ONT) have recently released a new flowcell (R10.4) and chemistry (Kit12), which reportedly generate per-read accuracies rivalling those of Illumina data. To evaluate this, we sequenced DNA extracts from four commonly studied bacterial pathogens, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, using Illumina and ONT's R9.4.1/Kit10, R10.3/Kit12, R10.4/Kit12 flowcells/chemistries. We compared raw read accuracy and assembly accuracy for each modality, considering the impact of different nanopore basecalling models, commonly used assemblers, sequencing depth, and the use of duplex versus simplex reads. 'Super accuracy' (sup) basecalled R10.4 reads - in particular duplex reads - have high per-read accuracies and could be used to robustly reconstruct bacterial genomes without the use of Illumina data. However, the per-run yield of duplex reads generated in our hands with standard sequencing protocols was low (typically <10 %), with substantial implications for cost and throughput if relying on nanopore data only to enable bacterial genome reconstruction. In addition, recovery of small plasmids with the best-performing long-read assembler (Flye) was inconsistent. R10.4/Kit12 combined with sup basecalling holds promise as a singular sequencing technology in the reconstruction of commonly studied bacterial genomes, but hybrid assembly (Illumina+R9.4.1 hac) currently remains the highest throughput, most robust, and cost-effective approach to fully reconstruct these bacterial genomes.

Genomic Insights into the Mobilome and Resistome of Sentinel Microorganisms Originating from Farms of Two Different Swine Production Systems

Antimicrobial resistance (AMR) is a threat to public health due to long-term antimicrobial use (AMU), which promotes the bacterial acquisition of antimicrobial resistance determinants (ARDs). Within food-producing animals, organic and extensive Iberian swine production is based on sustainable and eco-friendly management systems, providing an excellent opportunity to evaluate how sustained differences in AMU impact the development and spread of AMR. Here, through a whole-genome sequencing approach, we provide an in-depth characterization of the resistome and mobilome and their interaction in 466 sentinel bacteria, namely, Escherichia coli, Enterococcus spp., Campylobacter coli, and Staphylococcus spp., recovered from 37 intensive and organic-extensive pig farms. Both ARDs and mobile genetic elements (MGEs) were primarily taxon-associated, with higher similarities among bacteria which were closely phylogenetically related. E. coli exhibited the most diverse resistome and mobilome, with 85.4% mobilizable ARDs, 50.3% of which were plasmid-associated. Staphylococcus spp. exhibited a broad repertoire of ARDs and MGEs, with 52.3% of its resistome being mobilizable. Although Enterococcus spp. carried the highest number of ARDs per isolate and its plasmidome was similar in size to that of E. coli, 43.7% of its resistome was mobilizable. A narrow spectrum of ARDs constituted the C. coli resistome, with point mutations as its main AMR driver. A constrained AMU, as observed in organic-extensive herds, determined a reduction in the quantitative composition of the resistome and the complexity of the resistome-mobilome interaction. These results demonstrate taxon-associated AMR-MGE interactions and evidence that responsible AMU can contribute to reducing AMR pressure in the food chain. IMPORTANCE This study provides the first integral genomic characterization of the resistome and mobilome of sentinel microorganisms for antimicrobial resistance (AMR) surveillance from two different swine production systems. Relevant differences were observed among taxa in the resistomes and mobilomes they harbored, revealing their distinctive risk in AMR dissemination and spread. Thus, Escherichia coli and, to a lesser extent, Staphylococcus spp. constituted the main reservoirs of mobilizable antimicrobial resistance genes, which were predominantly plasmid-associated; in contrast to Campylobacter coli, whose resistome was mainly determined by point mutations. The reduced complexity of mobilome-resistome interaction in Enterococcus spp. suggested its limited role in AMR dissemination from swine farms. The significant differences in antimicrobial use among the studied farms allowed us to assess the suitability of whole-genome sequencing as a rapid and efficient technique for the assessment of mid- to long-term on-farm interventions for the reduction of antimicrobial use and the evaluation of AMR status.

Global Distribution and Diversity of Prevalent Sewage Water Plasmidomes

Sewage water from around the world contains an abundance of short plasmids, several of which harbor antimicrobial resistance genes (ARGs). The global dynamics of plasmid-derived antimicrobial resistance and functions are only starting to be unveiled. Here, we utilized a previously created data set of 159,332 assumed small plasmids from 24 different global sewage samples. The detailed phylogeny, as well as the interplay between their protein domains, ARGs, and predicted bacterial host genera, were investigated to understand sewage plasmidome dynamics globally. A total of 58,429 circular elements carried genes encoding plasmid-related features, and MASH distance analyses showed a high degree of diversity. A single (yet diverse) cluster of 520 predicted Acinetobacter plasmids was predominant among the European sewage water. Our results suggested a prevalence of plasmid-backbone gene combinations over others. This could be related to selected bacterial genera that act as bacterial hosts. These combinations also mirrored the geographical locations of the sewage samples. Our functional domain network analysis identified three groups of plasmids. However, these backbone domains were not exclusive to any given group, and Acinetobacter was the dominant host genus among the theta-replicating plasmids, which contained a reservoir of the macrolide resistance gene pair msr(E) and mph(E). Macrolide resistance genes were the most common in the sewage plasmidomes and were found in the largest number of unique plasmids. While msr(E) and mph(E) were limited to Acinetobacter, erm(B) was disseminated among a range of Firmicutes plasmids, including Staphylococcus and Streptococcus, highlighting a potential reservoir of antibiotic resistance for these pathogens from around the globe. IMPORTANCE Antimicrobial resistance is a global threat to human health, as it inhibits our ability to treat infectious diseases. This study utilizes sewage water plasmidomes to identify plasmid-derived features and highlights antimicrobial resistance genes, particularly macrolide resistance genes, as abundant in sewage water plasmidomes in Firmicutes and Acinetobacter hosts. The emergence of macrolide resistance in these bacteria suggests that macrolide selective pressure exists in sewage water and that the resident bacteria can readily acquire macrolide resistance via small plasmids.

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.

Systematic review of wastewater surveillance of antimicrobial resistance in human populations

OBJECTIVES

We systematically reviewed studies using wastewater for AMR surveillance in human populations, to determine: (i) evidence of concordance between wastewater-human AMR prevalence estimates, and (ii) methodological approaches which optimised identifying such an association, and which could be recommended as standard. We used Lin's concordance correlation coefficient (CCC) to quantify concordance between AMR prevalence estimates in wastewater and human compartments (where CCC = 1 reflects perfect concordance), and logistic regression to identify study features (e.g. sampling methods) associated with high agreement studies (defined as >70% of within-study wastewater-human AMR prevalence comparisons within ±10%).

RESULTS

Of 8,867 records and 441 full-text methods reviewed, 33 studies were included. AMR prevalence data was extractable from 24 studies conducting phenotypic-only (n = 7), genotypic-only (n = 1) or combined (n = 16) AMR detection. Overall concordance of wastewater-human AMR prevalence estimates was reasonably high for both phenotypic (CCC = 0.85 [95% CI 0.8-0.89]) and genotypic approaches (CCC = 0.88 (95% CI 0.84-0.9)) despite diverse study designs, bacterial species investigated and phenotypic/genotypic targets. No significant relationships between methodological approaches and high agreement studies were identified using logistic regression; however, this was limited by inconsistent reporting of study features, significant heterogeneity in approaches and limited sample size. Based on a secondary, descriptive synthesis, studies conducting composite sampling of wastewater influent, longitudinal sampling >12 months, and time-/location-matched sampling of wastewater and human compartments generally had higher agreement.

CONCLUSION

Wastewater-based surveillance of AMR appears promising, with high overall concordance between wastewater and human AMR prevalence estimates in studies irrespective of heterogenous approaches. However, our review suggests future work would benefit from: time-/location-matched sampling of wastewater and human populations, composite sampling of influent, and sampling >12 months for longitudinal studies. Further research and clear and consistent reporting of study methods is required to identify optimal practice.

Assessment of global health risk of antibiotic resistance genes

Antibiotic resistance genes (ARGs) have accelerated microbial threats to human health in the last decade. Many genes can confer resistance, but evaluating the relative health risks of ARGs is complex. Factors such as the abundance, propensity for lateral transmission and ability of ARGs to be expressed in pathogens are all important. Here, an analysis at the metagenomic level from various habitats (6 types of habitats, 4572 samples) detects 2561 ARGs that collectively conferred resistance to 24 classes of antibiotics. We quantitatively evaluate the health risk to humans, defined as the risk that ARGs will confound the clinical treatment for pathogens, of these 2561 ARGs by integrating human accessibility, mobility, pathogenicity and clinical availability. Our results demonstrate that 23.78% of the ARGs pose a health risk, especially those which confer multidrug resistance. We also calculate the antibiotic resistance risks of all samples in four main habitats, and with machine learning, successfully map the antibiotic resistance threats in global marine habitats with over 75% accuracy. Our novel method for quantitatively surveilling the health risk of ARGs will help to manage one of the most important threats to human and animal health.

Antibiotic resistance genes (ARGs) have accelerated microbial threats to human health. Here, Zhang et al. analyze 4572 metagenomic samples to illustrate the global patterns of ARG distribution in diverse habitats. They quantitatively evaluate the health risk to humans of 2561 ARGs by integrating human accessibility, mobility, pathogenicity and clinical availability. With the machine learning, they map the antibiotic resistance threats in global marine habitats.

Role of mobile genetic elements in the global dissemination of the carbapenem resistance gene blaNDM

The mobile resistance gene blaNDM encodes the NDM enzyme which hydrolyses carbapenems, a class of antibiotics used to treat some of the most severe bacterial infections. The blaNDM gene is globally distributed across a variety of Gram-negative bacteria on multiple plasmids, typically located within highly recombining and transposon-rich genomic regions, which leads to the dynamics underlying the global dissemination of blaNDM to remain poorly resolved. Here, we compile a dataset of over 6000 bacterial genomes harbouring the blaNDM gene, including 104 newly generated PacBio hybrid assemblies from clinical and livestock-associated isolates across China. We develop a computational approach to track structural variants surrounding blaNDM, which allows us to identify prevalent genomic contexts, mobile genetic elements, and likely events in the gene's global spread. We estimate that blaNDM emerged on a Tn125 transposon before 1985, but only reached global prevalence around a decade after its first recorded observation in 2005. The Tn125 transposon seems to have played an important role in early plasmid-mediated jumps of blaNDM, but was overtaken in recent years by other elements including IS26-flanked pseudo-composite transposons and Tn3000. We found a strong association between blaNDM-carrying plasmid backbones and the sampling location of isolates. This observation suggests that the global dissemination of the blaNDM gene was primarily driven by successive between-plasmid transposon jumps, with far more restricted subsequent plasmid exchange, possibly due to adaptation of plasmids to their specific bacterial hosts.

Remarkable genomic diversity among Escherichia isolates recovered from healthy chickens

The genus Escherichia has been extensively studied and it is known to encompass a range of commensal and pathogenic bacteria that primarily inhabit the gastrointestinal tracts of warm-blooded vertebrates. However, the presence of E. coli as a model organism and potential pathogen has diverted attention away from commensal strains and other species in the genus. To investigate the diversity of Escherichia in healthy chickens, we collected fecal samples from antibiotic-free Lohmann Brown layer hens and determined the genome sequences of 100 isolates, 81 of which were indistinguishable at the HC0 level of the Hierarchical Clustering of Core Genome Multi-Locus Sequence Typing scheme. Despite initial selection on CHROMagar Orientation medium, which is considered selective for E. coli, in silico phylotyping and core genome single nucleotide polymorphism analysis revealed the presence of at least one representative of all major clades of Escherichia, except for E. albertii, Shigella, and E. coli phylogroup B2 and cryptic clade I. The most frequent phylogenomic groups were E. coli phylogroups A and B1 and E. ruysiae (clades III and IV). We compiled a collection of reference strains isolated from avian sources (predominantly chicken), representing every Escherichia phylogroup and species, and used it to confirm the phylogeny and diversity of our isolates. Overall, the isolates carried low numbers of the virulence and antibiotic resistance genes typically seen in avian pathogenic E. coli. Notably, the clades not recovered are ones that have been most strongly associated with virulence by other studies.

Beyond to the Stable: Role of the Insertion Sequences as Epidemiological Descriptors in Corynebacterium striatum

In recent years, epidemiological studies of infectious agents have focused mainly on the pathogen and stable components of its genome. The use of these stable components makes it possible to know the evolutionary or epidemiological relationships of the isolates of a particular pathogen. Under this approach, focused on the pathogen, the identification of resistance genes is a complementary stage of a bacterial characterization process or an appendix of its epidemiological characterization, neglecting its genetic components’ acquisition or dispersal mechanisms. Today we know that a large part of antibiotic resistance is associated with mobile elements. Corynebacterium striatum, a bacterium from the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. Despite the different existing mobile genomic elements, there is evidence that acquired resistance genes are coupled to insertion sequences in C. striatum. This perspective article reviews the insertion sequences linked to resistance genes, their relationship to evolutionary lineages, epidemiological characteristics, and the niches the strains inhabit. Finally, we evaluate the potential of the insertion sequences for their application as a descriptor of epidemiological scenarios, allowing us to anticipate the emergence of multidrug-resistant lineages.

The potential of using E. coli as an indicator for the surveillance of antimicrobial resistance (AMR) in the environment

To understand the dynamics of antimicrobial resistance (AMR), in a One-Health perspective, surveillance play an important role. Monitoring systems already exist in the human health and livestock sectors, but there are no environmental monitoring programs. Therefore there is an urgent need to initiate environmental AMR monitoring programs nationally and globally, which will complement existing systems in different sectors. However, environmental programs should not only identify anthropogenic influences and levels of AMR, but they should also allow for identification of transmissions to and from human and animal populations. In the current review we therefore propose using antimicrobial resistant Escherichia coli as indicators for monitoring occurrence and levels of AMR in the environment, including wildlife.

A genomic epidemiological study shows that prevalence of antimicrobial resistance in Enterobacterales is associated with the livestock host, as well as antimicrobial usage

Enterobacterales from livestock are potentially important reservoirs for antimicrobial resistance (AMR) to pass through the food chain to humans, thereby increasing the AMR burden and affecting our ability to tackle infections. In this study 168 isolates from four genera of the order Enterobacterales, primarily Escherichia coli, were purified from livestock (cattle, pigs and sheep) faeces from 14 farms in the United Kingdom. Their genomes were resolved using long- and short-read sequencing to analyse AMR genes and their genetic context, as well as to explore the relationship between AMR burden and on-farm antimicrobial usage (AMU), in the three months prior to sampling. Although E. coli isolates were genomically diverse, phylogenetic analysis using a core-genome SNP tree indicated pig isolates to generally be distinct from sheep isolates, with cattle isolates being intermediates. Approximately 28 % of isolates harboured AMR genes, with the greatest proportion detected in pigs, followed by cattle then sheep; pig isolates also harboured the highest number of AMR genes per isolate. Although 90 % of sequenced isolates harboured diverse plasmids, only 11 % of plasmids (n=58 out of 522) identified contained AMR genes, with 91 % of AMR plasmids being from pig, 9 % from cattle and none from sheep isolates; these results indicated that pigs were a principle reservoir of AMR genes harboured by plasmids and likely to be involved in their horizontal transfer. Significant associations were observed between AMU (mg kg⁻¹) and AMR. As both the total and the numbers of different antimicrobial classes used on-farm increased, the risk of multi-drug resistance (MDR) in isolates rose. However, even when AMU on pig farms was comparatively low, pig isolates had increased likelihood of being MDR; harbouring relatively more resistances than those from other livestock species. Therefore, our results indicate that AMR prevalence in livestock is not only influenced by recent AMU on-farm but also livestock-related factors, which can influence the AMR burden in these reservoirs and its plasmid mediated transmission.

Isolation of Ciprofloxacin and Ceftazidime-Resistant Enterobacterales From Vegetables and River Water Is Strongly Associated With the Season and the Sample Type

The dissemination of antibiotic-resistant bacteria (ARB) from water used for crop irrigation to vegetables is poorly studied. During a year, five farmer markets in a city in Central Chile were visited, and 478 vegetable samples (parsleys, corianders, celeries, lettuces, chards, and beets) were collected. Simultaneously, 32 water samples were collected from two rivers which are used to irrigate the vegetables produced in the area. Resistant Enterobacterales were isolated and identified. Colistin resistance gene mcr-1 and extended spectrum β-lactamases (ESBL) were molecularly detected. The association of environmental factors was evaluated, with the outcomes being the presence of Enterobacterales resistant to four antibiotic families and the presence of multidrug resistance (MDR) phenotypes. Parsley, coriander, and celery showed the highest prevalence of resistant Enterobacterales (41.9% for ciprofloxacin and 18.5% for ceftazidime). A total of 155 isolates were obtained, including Escherichia coli (n=109), Citrobacter sp. (n=20), Enterobacter cloacae complex (n=8), Klebsiella pneumoniae (n=8), and Klebsiella aerogenes (n=1). Resistance to ampicillin (63.2%) and ciprofloxacin (74.2%) was most frequently found; 34.5% of the isolates showed resistance to third-generation cephalosporins, and the MDR phenotype represented 51.6% of the isolates. In two E. coli isolates (1.29%), the gene mcr-1 was found and ESBL genes were found in 23/62 isolates (37%), with bla_CTX-M being the most frequently found in 20 isolates (32%). Resistant Enterobacterales isolated during the rainy season were less likely to be MDR as compared to the dry season. Understanding environmental associations represent the first step toward an improved understanding of the public health impact of ARB in vegetables and water.

Ten-year longitudinal molecular epidemiology study of Escherichia coli and Klebsiella species bloodstream infections in Oxfordshire, UK

Background

The incidence of Gram-negative bloodstream infections (BSIs), predominantly caused by Escherichia coli and Klebsiella species, continues to increase; however, the causes of this are unclear and effective interventions are therefore hard to design.

Methods

In this study, we sequenced 3468 unselected isolates over a decade in Oxfordshire (UK) and linked this data to routinely collected electronic healthcare records and mandatory surveillance reports. We annotated genomes for clinically relevant genes, contrasting the distribution of these within and between species, and compared incidence trends over time using stacked negative binomial regression.

Results

We demonstrate that the observed increases in E. coli incidence were not driven by the success of one or more sequence types (STs); instead, four STs continue to dominate a stable population structure, with no evidence of adaptation to hospital/community settings. Conversely in Klebsiella spp., most infections are caused by sporadic STs with the exception of a local drug-resistant outbreak strain (ST490). Virulence elements are highly structured by ST in E. coli but not Klebsiella spp. where they occur in a diverse spectrum of STs and equally across healthcare and community settings. Most clinically hypervirulent (i.e. community-onset) Klebsiella BSIs have no known acquired virulence loci. Finally, we demonstrate a diverse but largely genus-restricted mobilome with close associations between antimicrobial resistance (AMR) genes and insertion sequences but not typically specific plasmid replicon types, consistent with the dissemination of AMR genes being highly contingent on smaller mobile genetic elements (MGEs).

Conclusions

Our large genomic study highlights distinct differences in the molecular epidemiology of E. coli and Klebsiella BSIs and suggests that no single specific pathogen genetic factors (e.g. AMR/virulence genes/sequence type) are likely contributing to the increasing incidence of BSI overall, that association with AMR genes in E. coli is a contributor to the increasing number of E. coli BSIs, and that more attention should be given to AMR gene associations with non-plasmid MGEs to try and understand horizontal gene transfer networks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00947-2.

Genomic network analysis of environmental and livestock F-type plasmid populations

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum β-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.