Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

Genetic diversity of KPC-2-producing Klebsiella pneumoniae complex from aquatic ecosystems

During the COVID-19 pandemic, the occurrence of carbapenem-resistant Klebsiella pneumoniae increased in human clinical settings worldwide. Impacted by this increase, international high-risk clones harboring carbapenemase-encoding genes have been circulating in different sources, including the environment. The blaKPC gene is the most commonly disseminated carbapenemase-encoding gene worldwide, whose transmission is carried out by different mobile genetic elements. In this study, blaKPC-2-positive Klebsiella pneumoniae complex strains were isolated from different anthropogenically affected aquatic ecosystems and characterized using phenotypic, molecular, and genomic methods. K. pneumoniae complex strains exhibited multidrug-resistant and extensively drug-resistant profiles, spotlighting the resistance to carbapenems, ceftazidime-avibactam, colistin, and tigecycline, which are recognized as last-line antimicrobial treatment options. Molecular analysis showed the presence of several antimicrobial resistance, virulence, and metal tolerance genes. In-depth analysis showed that the blaKPC-2 gene was associated with three different Tn4401 isoforms (i.e., Tn4401a, Tn4401b, and Tn4401i) and NTEKPC elements. Different plasmid replicons were detected and a conjugative IncN-pST15 plasmid harboring the blaKPC-2 gene associated with Tn4401i was highlighted. K. pneumoniae complex strains belonging to international high-risk (e.g., ST11 and ST340) and unusual clones (e.g., ST323, ST526, and ST4216) previously linked to clinical settings. In this context, some clones were reported for the first time in the environmental sector. Therefore, these findings evidence the occurrence of carbapenemase-producing K. pneumoniae complex strains in aquatic ecosystems and contribute to the monitoring of carbapenem resistance worldwide.

Antimicrobial Lasso Peptide Cloacaenodin Utilizes a Unique TonB-Dependent Transporter to Access Susceptible Bacteria

The development of new antimicrobial agents effective against Gram-negative bacteria remains a major challenge in drug discovery. The lasso peptide cloacaenodin has potent antimicrobial activity against multiple strains in the Enterobacter genus, one of the ESKAPE pathogens. Here, we show that cloacaenodin uses a previously uncharacterized TonB-dependent transporter, which we name CloU, to cross the outer membrane (OM) of susceptible bacteria. Inner membrane transport is mediated by the protein SbmA. CloU is distinct from the known OM transporters (FhuA and PupB) utilized by other antimicrobial lasso peptides and thus offers important insight into the spectrum of activity of cloacaenodin. Using knowledge of the transport pathway to predict other cloacaenodin-susceptible strains, we demonstrate the activity of cloacaenodin against clinical isolates of Enterobacter and of a Kluyvera strain. Further, we use molecular dynamics simulations and mutagenesis of CloU to explain the variation in cloacaenodin susceptibility observed across different strains of Enterobacter. This work expands the currently limited understanding of lasso peptide uptake and advances the potential of cloacaenodin as an antibiotic.

Genomic Characterization of a Carbapenem-Resistant Raoultella planticola Strain Co-Harboring blaIMP-4 and blaSHV-12 Genes

Raoultella planticola is an emerging bacterial pathogen responsible for causing infections in both humans and animals. Unfortunately, sporadic reports of carbapenem-resistant R. planticola (CRRP) have been documented worldwide. Here we first reported the complete genome sequence of a CRRP isolate RP_3045 co-carrying blaIMP-4 and blaSHV-12, recovered from a patient in China, and its genetic relatedness to 82 R. planticola strains deposited in the NCBI GenBank database, sourced from humans, animals, and the environment. Whole-genome sequencing was performed using the Illumina NovaSeq 6000 and Oxford Nanopore MinION platforms. Phylogenetic analysis was also performed and visualized using a single nucleotide polymorphism (SNP)-based strategy. The complete genome of R. planticola strain RP_3045 was determined to be 6,312,961 bp in length, comprising five contigs that included one chromosome and four plasmids. RP_3045 was found to be multidrug-resistant and harbored several antimicrobial resistance genes, including both blaIMP-4 and blaSHV-12 genes located on a single plasmid. The most closely related strain was hkcpe63, recovered from humans in Hong Kong, China, in 2014, with 506 SNP differences. R. planticola strains were distributed globally and exhibited strong associations among isolates obtained from different sectors. This study provides evidence for the potential of R. planticola to disseminate carbapenem resistance across different sectors, highlighting the critical need for active and continuous surveillance of CRRP.

Clinical and Genomic Characterization of a Cohort of Patients With Klebsiella pneumoniae Bloodstream Infection

BACKGROUND

The clinical and microbial factors associated with Klebsiella pneumoniae bloodstream infections (BSIs) are not well characterized. Prior studies have focused on highly resistant or hypervirulent isolates, limiting our understanding of K. pneumoniae strains that commonly cause BSI. We performed a record review and whole-genome sequencing to investigate the clinical characteristics, bacterial diversity, determinants of antimicrobial resistance, and risk factors for in-hospital death in a cohort of patients with K. pneumoniae BSI.

METHODS

We identified 562 patients at Massachusetts General Hospital with K. pneumoniae BSIs between 2016 and 2022. We collected data on comorbid conditions, infection source, clinical outcomes, and antibiotic resistance and performed whole-genome sequencing on 108 sequential BSI isolates from 2021 to 2022.

RESULTS

Intra-abdominal infection was the most common source of infection accounting for 34% of all BSIs. A respiratory tract source accounted for 6% of BSIs but was associated with a higher in-hospital mortality rate (adjusted odds ratio, 5.4 [95% confidence interval, 2.2-12.8]; P < .001 for comparison with other sources). Resistance to the first antibiotic prescribed was also associated with a higher risk of death (adjusted odds ratio, 5.2 [95% confidence interval, 2.2-12.4]; P < .001). BSI isolates were genetically diverse, and no clusters of epidemiologically and genetically linked cases were observed. Virulence factors associated with invasiveness were observed at a low prevalence, although an unexpected association between O-antigen type and the source of infection was found.

CONCLUSIONS

These observations demonstrate the versatility of K. pneumoniae as an opportunistic pathogen and highlight the need for new approaches for surveillance and the rapid identification of patients with invasive antimicrobial-resistant K. pneumoniae infection.

The carbapenem inoculum effect provides insight into the molecular mechanisms underlying carbapenem resistance in Enterobacterales

Carbapenem-resistant Enterobacterales (CRE) are important pathogens that can develop resistance via multiple molecular mechanisms, including hydrolysis or reduced antibiotic influx. Identifying these mechanisms can improve pathogen surveillance, infection control, and patient care. We investigated how resistance mechanisms influence the carbapenem inoculum effect (IE), a phenomenon where inoculum size affects antimicrobial susceptibility testing (AST). We demonstrated that seven different carbapenemases impart a meropenem IE in Escherichia coli. Across 110 clinical CRE isolates, the carbapenem IE strictly depended on resistance mechanism: all carbapenemase-producing CRE (CP-CRE) exhibited a strong IE, whereas porin-deficient CRE displayed none. Concerningly, 50% and 24% of CP-CRE isolates changed susceptibility classification to meropenem and ertapenem, respectively, across the allowable inoculum range in clinical guidelines. The meropenem IE, and the ratio of ertapenem to meropenem minimal inhibitory concentration (MIC) at standard inoculum, reliably identified CP-CRE. Understanding how resistance mechanisms affect AST could improve diagnosis and guide therapies for CRE infections.

SHIP: identifying antimicrobial resistance gene transfer between plasmids

MOTIVATION

Plasmids are carriers for antimicrobial resistance (AMR) genes and can exchange genetic material with other structures, contributing to the spread of AMR. There is no reliable approach to identify the transfer of AMR genes across plasmids. This is mainly due to the absence of a method to assess the phylogenetic distance of plasmids, as they show large DNA sequence variability. Identifying and quantifying such transfer can provide novel insight into the role of small mobile elements and resistant plasmid regions in the spread of AMR.

RESULTS

We developed SHIP, a novel method to quantify plasmid similarity based on the dynamics of plasmid evolution. This allowed us to find conserved fragments containing AMR genes in structurally different and phylogenetically distant plasmids, which is evidence for lateral transfer. Our results show that regions carrying AMR genes are highly mobilizable between plasmids through transposons, integrons, and recombination events, and contribute to the spread of AMR. Identified transferred fragments include a multi-resistant complex class 1 integron in Escherichia coli and Klebsiella pneumoniae, and a region encoding tetracycline resistance transferred through recombination in Enterococcus faecalis.

AVAILABILITY AND IMPLEMENTATION

The code developed in this work is available at https://github.com/AbeelLab/plasmidHGT.

Phylogenetic lineages and antimicrobial resistance determinants of clinical Klebsiella oxytoca spanning local to global scales

Klebsiella oxytoca is an opportunistic pathogen causing serious nosocomial infections. Knowledge about the population structure and diversity of healthcare-associated K. oxytoca from a genomic standpoint remains limited. Here, we characterized the phylogenetic relationships and genomic characteristics of 20 K. oxytoca sensu stricto isolates recovered from bloodstream infections at the Dartmouth-Hitchcock Medical Center, New Hampshire, USA from 2017 to 2021. Results revealed a diverse population consisting of 15 sequence types (STs) that together harbored 10 variants of the intrinsic beta-lactamase gene bla OXY-2, conferring resistance to penicillins. Similar sets of antimicrobial resistance (AMR) determinants reside in multiple distinct lineages, with no one lineage dominating the local population. To place the New Hampshire K. oxytoca in a broader context, we compared them to 304 publicly available genomes of clinical isolates from 18 countries. This global clinical K. oxytoca sensu stricto population is represented by over 65 STs that together harbored resistance genes against 14 antimicrobial classes, including eight bla OXY-2 variants. Three dominant STs in the global population (ST2, ST176, ST199) circulate across multiple countries and were also present in the New Hampshire population. The global K. oxytoca population is genetically diverse, but there is evidence for broad dissemination of a few lineages carrying distinct set of AMR determinants. Our findings reveal the clinical diversity of K. oxytoca sensu stricto and its importance in surveillance efforts aimed at monitoring the evolution of this drug-resistant nosocomial pathogen. IMPORTANCE The opportunistic pathogen Klebsiella oxytoca has been increasingly implicated in patient morbidity and mortality worldwide, including several outbreaks in healthcare settings. The emergence and spread of antimicrobial resistant strains exacerbate the disease burden caused by this species. Our study showed that clinical K. oxytoca sensu stricto is phylogenetically diverse, harboring various antimicrobial resistance determinants and bla OXY-2 variants. Understanding the genomic and population structure of K. oxytoca is important for international initiatives and local epidemiological efforts for surveillance and control of drug-resistant K. oxytoca.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Empirically derived sequence similarity thresholds to study the genomic epidemiology of plasmids shared among healthcare-associated bacterial pathogens

BACKGROUND

Healthcare-associated bacterial pathogens frequently carry plasmids that contribute to antibiotic resistance and virulence. The horizontal transfer of plasmids in healthcare settings has been previously documented, but genomic and epidemiologic methods to study this phenomenon remain underdeveloped. The objectives of this study were to apply whole-genome sequencing to systematically resolve and track plasmids carried by nosocomial pathogens in a single hospital, and to identify epidemiologic links that indicated likely horizontal plasmid transfer.

METHODS

We performed an observational study of plasmids circulating among bacterial isolates infecting patients at a large hospital. We first examined plasmids carried by isolates sampled from the same patient over time and isolates that caused clonal outbreaks in the same hospital to develop thresholds with which horizontal plasmid transfer within a tertiary hospital could be inferred. We then applied those sequence similarity thresholds to perform a systematic screen of 3074 genomes of nosocomial bacterial isolates from a single hospital for the presence of 89 plasmids. We also collected and reviewed data from electronic health records for evidence of geotemporal links between patients infected with bacteria encoding plasmids of interest.

FINDINGS

Our analyses determined that 95% of analyzed genomes maintained roughly 95% of their plasmid genetic content and accumulated fewer than 15 SNPs per 100 kb of plasmid sequence. Applying these similarity thresholds to identify horizontal plasmid transfer identified 45 plasmids that potentially circulated among clinical isolates. Ten highly preserved plasmids met criteria for geotemporal links associated with horizontal transfer. Several plasmids with shared backbones also encoded different additional mobile genetic element content, and these elements were variably present among the sampled clinical isolate genomes.

INTERPRETATION

Evidence suggests that the horizontal transfer of plasmids among nosocomial bacterial pathogens appears to be frequent within hospitals and can be monitored with whole genome sequencing and comparative genomics approaches. These approaches should incorporate both nucleotide identity and reference sequence coverage to study the dynamics of plasmid transfer in the hospital.

FUNDING

This research was supported by the US National Institute of Allergy and Infectious Disease (NIAID) and the University of Pittsburgh School of Medicine.

Comparison of long- and short-read metagenomic assembly for low-abundance species and resistance genes

Recent technological and computational advances have made metagenomic assembly a viable approach to achieving high-resolution views of complex microbial communities. In previous benchmarking, short-read (SR) metagenomic assemblers had the highest accuracy, long-read (LR) assemblers generated the most contiguous sequences and hybrid (HY) assemblers balanced length and accuracy. However, no assessments have specifically compared the performance of these assemblers on low-abundance species, which include clinically relevant organisms in the gut. We generated semi-synthetic LR and SR datasets by spiking small and increasing amounts of Escherichia coli isolate reads into fecal metagenomes and, using different assemblers, examined E. coli contigs and the presence of antibiotic resistance genes (ARGs). For ARG assembly, although SR assemblers recovered more ARGs with high accuracy, even at low coverages, LR assemblies allowed for the placement of ARGs within longer, E. coli-specific contigs, thus pinpointing their taxonomic origin. HY assemblies identified resistance genes with high accuracy and had lower contiguity than LR assemblies. Each assembler type's strengths were maintained even when our isolate was spiked in with a competing strain, which fragmented and reduced the accuracy of all assemblies. For strain characterization and determining gene context, LR assembly is optimal, while for base-accurate gene identification, SR assemblers outperform other options. HY assembly offers contiguity and base accuracy, but requires generating data on multiple platforms, and may suffer high misassembly rates when strain diversity exists. Our results highlight the trade-offs associated with each approach for recovering low-abundance taxa, and that the optimal approach is goal-dependent.

Domestication and microbiome succession may drive pathogen spillover

Emerging infectious diseases have posed growing medical, social and economic threats to humanity. The biological background of pathogen spillover or host switch, however, still has to be clarified. Disease ecology finds pathogen spillovers frequently but struggles to explain at the molecular level. Contrarily, molecular biological traits of host-pathogen relationships with specific molecular binding mechanisms predict few spillovers. Here we aim to provide a synthetic explanation by arguing that domestication, horizontal gene transfer even between superkingdoms as well as gradual exchange of microbiome (microbiome succession) are essential in the whole scenario. We present a new perspective at the molecular level which can explain the observations of frequent pathogen spillover events at the ecological level. This proposed rationale is described in detail, along with supporting evidence from the peer-reviewed literature and suggestions for testing hypothesis validity. We also highlight the importance of systematic monitoring of virulence genes across taxonomical categories and in the whole biosphere as it helps prevent future epidemics and pandemics. We conclude that that the processes of domestication, horizontal gene transfer and microbial succession might be important mechanisms behind the many spillover events driven and accelerated by climate change, biodiversity loss and globalization.

Cloacaenodin, an Antimicrobial Lasso Peptide with Activity against Enterobacter

Using genome mining and heterologous expression, we report the discovery and production of a new antimicrobial lasso peptide from species related to the Enterobacter cloacae complex. Using NMR and mass spectrometric analysis, we show that this lasso peptide, named cloacaenodin, employs a threaded lasso fold which imparts proteolytic resistance that its unthreaded counterpart lacks. Cloacaenodin has selective, low micromolar, antimicrobial activity against species related to the E. cloacae complex, including species implicated in nosocomial infections and against clinical isolates of carbapenem-resistant Enterobacterales. We further used site-directed mutagenesis to probe the importance of specific residues to the peptide's biosynthesis, stability, and bioactivity.