Nested Russian Doll-Like Genetic Mobility Drives Rapid Dissemination of the Carbapenem Resistance Gene blaKPC

Development of two quantitative PCR assays for the detection of emerging opportunistic human pathogens belonging to the genus Phytobacter in routine diagnostics

The genus Phytobacter (fam. Enterobacteriaceae) includes species like Phytobacter diazotrophicus and Phytobacter ursingii, which have emerged as opportunistic human pathogens, particularly affecting vulnerable populations such as pre-term infants and immunocompromised patients. Traditional biochemical and molecular methods have struggled to accurately identify Phytobacter species in clinical diagnostics. This study addresses the issue by developing and validating two quantitative PCR (qPCR) assays using SYBR® Green I and TaqMan® technologies, targeting the nitrogen fixation regulatory gene (nifL) of Phytobacter spp. The SYBR® Green I assay showed a detection limit of a single cell per reaction, while the TaqMan® assay was easier to interpret due to the absence of background noise. These assays, validated with clinical isolates from Brazil, identified multiple new Phytobacter isolates, including a potentially novel species, providing improved diagnostic tools for detecting Phytobacter spp. and aiding in better clinical management.

Antibiotic resistance plasmids in Enterobacteriaceae isolated from fresh produce in northern Germany

In this study, the genomes of 22 Enterobacteriaceae isolates from fresh produce and herbs obtained from retail markets in northern Germany were completely sequenced with MiSeq short-read and MinION long-read sequencing and assembled using a Unicycler hybrid assembly. The data showed that 17 of the strains harbored between one and five plasmids, whereas in five strains, only the circular chromosomal DNA was detected. In total, 38 plasmids were identified. The size of the plasmids detected varied between ca. 2,000 and 326,000 bp, and heavy metal resistance genes were found on seven (18.4%) of the plasmids. Eleven plasmids (28.9%) showed the presence of antibiotic resistance genes. Among large plasmids (>32,000 bp), IncF plasmids (specifically, IncFIB and IncFII) were the most abundant replicon types, while all small plasmids were Col-replicons. Six plasmids harbored unit and composite transposons carrying antibiotic resistance genes, with IS26 identified as the primary insertion sequence. Class 1 integrons carrying antibiotic resistance genes were also detected on chromosomes of two Citrobacter isolates and on four plasmids. Mob-suite analysis revealed that 36.8% of plasmids in this study were found to be conjugative, while 28.9% were identified as mobilizable. Overall, our study showed that Enterobacteriaceae from fresh produce possess antibiotic resistance genes on both chromosome and plasmid, some of which are considered to be transferable. This indicates the potential for Enterobacteriaceae from fresh produce that is usually eaten in the raw state to contribute to the transfer of resistance genes to bacteria of the human gastrointestinal system.

IMPORTANCE

This study showed that Enterobacteriaceae from raw vegetables carried plasmids ranging in size from 2,715 to 326,286 bp, of which about less than one-third carried antibiotic resistance genes encoding resistance toward antibiotics such as tetracyclines, aminoglycosides, fosfomycins, sulfonamides, quinolones, and β-lactam antibiotics. Some strains encoded multiple resistances, and some encoded extended-spectrum β-lactamases. The study highlights the potential of produce, which may be eaten raw, as a potential vehicle for the transfer of antibiotic-resistant bacteria.

Real-time plasmid transmission detection pipeline

The spread of antimicrobial resistance among bacteria by horizontal plasmid transmissions poses a major challenge for clinical microbiology. Here, we evaluate a new real-time plasmid transmission detection pipeline implemented in the SeqSphere+ (Ridom GmbH, Münster, Germany) software. Within the pipeline, a local Mash plasmid database is created, and Mash searches with a distance threshold of 0.001 are used to trigger plasmid transmission early warning alerts (EWAs). Clonal transmissions are detected using core-genome multi-locus sequence typing allelic differences. The tools MOB-suite, NCBI AMRFinderPlus, CGE MobileElementFinder, pyGenomeViz, and MUMmer, integrated in SeqSphere+, are used to characterize plasmids and for visual pairwise plasmid comparisons, respectively. We evaluated the pipeline using published hybrid assemblies (Oxford Nanopore Technology/Illumina) of a surveillance and outbreak data set with plasmid transmissions. To emulate prospective usage, samples were imported in chronological order of sampling date. Different combinations of the user-adjustable parameters sketch size (1,000 vs 10,000) and plasmid size correction were tested, and discrepancies between resulting clusters were analyzed with Quast. When using a sketch size of 1,000 with size correction turned on, the SeqSphere+ pipeline agreed with the published data and produced the same clonal and carbapenemase-carrying plasmid clusters. EWAs were in the correct chronological order. In summary, the developed pipeline presented here is suitable for integration into clinical microbiology settings with limited bioinformatics knowledge due to its automated analyses and alert system, which are combined with the GUI-based SeqSphere+ platform. Thus, with its integrated sample database, (near) real-time plasmid transmission detection is within reach in bacterial routine-diagnostic settings when long-read sequencing is employed.

IMPORTANCE

Plasmid-mediated spread of antimicrobial resistance is a major challenge for clinical microbiology, and monitoring of potential plasmid transmissions is essential to combat further dissemination. Whole-genome sequencing is often used to surveil nosocomial transmissions but usually limited to the detection of clonal transmissions (based on chromosomal markers). Recent advances in long-read sequencing technologies enable full reconstruction of plasmids and the detection of very similar plasmids, but so far, easy-to-use bioinformatic tools for this purpose have been missing. Here, we present an evaluation of an innovative real-time plasmid transmission detection pipeline. It is integrated into the GUI-based SeqSphere+ software, which already offers core-genome multi-locus sequence typing-based pathogen outbreak detection. It requires very limited bioinformatics knowledge, and its database, automated analyses, and alert system make it well suited for prospective clinical application.

Various plasmid strategies limit the effect of bacterial restriction-modification systems against conjugation

In bacteria, genes conferring antibiotic resistance are mostly carried on conjugative plasmids, mobile genetic elements that spread horizontally between bacterial hosts. Bacteria carry defence systems that defend them against genetic parasites, but how effective these are against plasmid conjugation is poorly understood. Here, we study to what extent restriction-modification (RM) systems-by far the most prevalent bacterial defence systems-act as a barrier against plasmids. Using 10 different RM systems and 13 natural plasmids conferring antibiotic resistance in Escherichia coli, we uncovered variation in defence efficiency ranging from none to 105-fold protection. Further analysis revealed genetic features of plasmids that explain the observed variation in defence levels. First, the number of RM recognition sites present on the plasmids generally correlates with defence levels, with higher numbers of sites being associated with stronger defence. Second, some plasmids encode methylases that protect against restriction activity. Finally, we show that a high number of plasmids in our collection encode anti-restriction genes that provide protection against several types of RM systems. Overall, our results show that it is common for plasmids to encode anti-RM strategies, and that, as a consequence, RM systems form only a weak barrier for plasmid transfer by conjugation.

Applying rearrangement distances to enable plasmid epidemiology with pling

Plasmids are a key vector of antibiotic resistance, but the current bioinformatics toolkit is not well suited to tracking them. The rapid structural changes seen in plasmid genomes present considerable challenges to evolutionary and epidemiological analysis. Typical approaches are either low resolution (replicon typing) or use shared k-mer content to define a genetic distance. However, this distance can both overestimate plasmid relatedness by ignoring rearrangements, and underestimate by over-penalizing gene gain/loss. Therefore a model is needed which captures the key components of how plasmid genomes evolve structurally - through gene/block gain or loss, and rearrangement. A secondary requirement is to prevent promiscuous transposable elements (TEs) leading to over-clustering of unrelated plasmids. We choose the 'Double Cut and Join Indel' (DCJ-Indel) model, in which plasmids are studied at a coarse level, as a sequence of signed integers (representing genes or aligned blocks), and the distance between two plasmids is the minimum number of rearrangement events or indels needed to transform one into the other. We show how this gives much more meaningful distances between plasmids. We introduce a software workflow pling (https://github.com/iqbal-lab-org/pling), which uses the DCJ-Indel model, to calculate distances between plasmids and then cluster them. In our approach, we combine containment distances and DCJ-Indel distances to build a TE-aware plasmid network. We demonstrate superior performance and interpretability to other plasmid clustering tools on the 'Russian Doll' dataset and a hospital transmission dataset.

Genomic epidemiology and longitudinal sampling of ward wastewater environments and patients reveals complexity of the transmission dynamics of bla KPC-carbapenemase-producing Enterobacterales in a hospital setting

Background

Healthcare-associated wastewater and asymptomatic patient reservoirs colonized by carbapenemase-producing Enterobacterales (CPE) contribute to nosocomial CPE dissemination, but the characteristics and dynamics of this remain unclear.

Methods

We systematically sampled wastewater sites (n = 4488 samples; 349 sites) and patients (n = 1247) across six wards over 6-12 months to understand blaKPC-associated CPE (KPC-E) diversity within these reservoirs and transmission in a healthcare setting. Up to five KPC-E-positive isolates per sample were sequenced (Illumina). Recombination-adjusted phylogenies were used to define genetically related strains; assembly and mapping-based approaches were used to characterize antimicrobial resistance genes, insertion sequences (ISs) and Tn4401 types/target site sequences. The accessory genome was evaluated in some of the largest clusters, and those crossing reservoirs.

Results

Wastewater site KPC-E-positivity was substantial [101/349 sites (28.9%); 228/5601 (4.1%) patients cultured]. Thirteen KPC-E species and 109 strains were identified using genomics, and 24% of wastewater and 26% of patient KPC-E-positive samples harboured one or more strains. Most diversity was explained by the individual niche, suggesting localized factors are important in selection and spread. Tn4401 + flanking target site sequence diversity was greater in wastewater sites (P < 0.001), which might favour Tn4401-associated transposition/evolution. Shower/bath- and sluice/mop-associated sites were more likely to be KPC-E-positive (adjusted OR = 2.69; 95% CI: 1.44-5.01; P = 0.0019; and adjusted OR = 2.60; 95% CI: 1.04-6.52; P = 0.0410, respectively). Different strains had different blaKPC dissemination dynamics.

Conclusions

We identified substantial and diverse KPC-E colonization of wastewater sites and patients in this hospital setting. Reservoir and niche-specific factors (e.g. microbial interactions, selection pressures), and different strains and mobile genetic elements likely affect transmission dynamics. This should be considered in surveillance and control strategies.

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 instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of ortholog groups for individual protein-encoding genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of evolutionary statistics for each inferred ortholog group. Importantly, in comparison to tools for visual exploration of homologous relationships between gene clusters, zol can scale to thousands of gene cluster instances and produce detailed reports that are easy to digest. To showcase fai and zol, we apply them for: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-level genetic insights of two common BGCs in the fungal species Aspergillus flavus, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Phytobacter spp: the emergence of a new genus of healthcare-associated Enterobacterales encoding carbapenemases in Argentina: a case series

Members of the genus Phytobacter (order Enterobacterales) are isolated from the natural environment and clinical settings. Identification of Phytobacter strains based on biochemical characteristics is complicated due to taxonomic confusion, and they are often misidentified by automated identification systems in laboratories. In this study we describe the first three clinical cases associated with Phytobacter spp. reported in Argentina. We describe the identification, the molecular analysis using whole genome sequencing and the potential clinical relevance.

Bacterial lifestyle shapes pangenomes

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

Hybracter: enabling scalable, automated, complete and accurate bacterial genome assemblies

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants. They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.

Hybracter: Enabling Scalable, Automated, Complete and Accurate Bacterial Genome Assemblies

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants (SNVs). They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance (AMR) genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic, and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.

Resistome, mobilome, and virulome explored in clinical isolates derived from acne patients in Egypt: unveiling unique traits of an emerging coagulase-negative Staphylococcus pathogen

Coagulase-negative staphylococci (CoNS) are a group of gram-positive staphylococcal species that naturally inhabit the healthy human skin and mucosa. The clinical impact of CoNS-associated infections has recently been regarded as a challenge for diagnosis and therapeutic options. CoNS-associated infections are primarily caused by bacterial resistance to antibiotics and biofilm formation. As antibiotics are still the most used treatment, this problem will likely persist in the future. The present study aimed to investigate the resistance and virulence of CoNS recovered from various acne lesions and explore their genetic basis. Skin swab samples were collected from participants with acne and healthy skin. All samples underwent conventional culture for the isolation of CoNS, MALDI-TOF confirmation, antibiotic susceptibility, and biofilm formation testing. A total of 85 CoNS isolates were recovered from the samples and preliminarily identified as Staphylococcus epidermidis. Isolates from the acne group (n = 60) showed the highest rates of resistance to penicillin (73%), cefoxitin (63%), clindamycin (53.3%), and erythromycin (48%), followed by levofloxacin (36.7%) and gentamycin (31.7%). The lowest rates of resistance were observed against tetracycline (28.3%), doxycycline (11.7%), and minocycline (8.3%). CoNS isolated from mild, moderate acne and healthy isolates did not show strong biofilm formation, whereas the isolates from the severe cases of the acne group showed strong biofilm formation (76.6%). Four extensively drug-resistant and strong biofilm-forming staphylococcal isolates recovered from patients with severe acne were selected for whole-genome sequencing (WGS), and their genomes were investigated using bioinformatics tools. Three of the sequenced genomes were identified as S. epidermidis; however, isolate 29AM was identified as Staphylococcus warneri, which is a newly emerging pathogen that is not commonly associated with acne and was not detected by MALDI-TOF. All the sequenced strains were multidrug-resistant and carried multiple resistance genes, including blaZ, mecA, tet(K), erm(C), lnuA, vgaA, dfrC, fusB, fosBx1, norA, and vanT, which were found to be located on plasmids and chromosomes. Virulence features were detected in all genomes in the presence of genes involved in adherence and biofilm formation (icaA, icaB, icaC, sdrG, sdrH, atl, ebh, and ebp). Only the S. warneri isolate 29AM contained immune evasion genes (capB, capC, acpXL, and manA), an anti-phagocytosis gene (cdsA), and other unique features. As a result of their potential pathogenicity and antibiotic resistance, CoNS must be monitored as an emerging pathogen associated with acne infections. To the best of our knowledge, this is the first report to isolate, identify, and correlate S. warneri with severe acne infections among Egyptian patients using WGS and bioinformatic analysis.

Synergistic horizontal transfer of antibiotic resistance genes and transposons in the infant gut microbial genome

Transposons, plasmids, bacteriophages, and other mobile genetic elements facilitate horizontal gene transfer in the gut microbiota, allowing some pathogenic bacteria to acquire antibiotic resistance genes (ARGs). Currently, the relationship between specific ARGs and specific transposons in the comprehensive infant gut microbiome has not been elucidated. In this study, ARGs and transposons were annotated from the Unified Human Gastrointestinal Genome (UHGG) and the Early-Life Gut Genomes (ELGG). Association rules mining was used to explore the association between specific ARGs and specific transposons in UHGG, and the robustness of the association rules was validated using the external database in ELGG. Our results suggested that ARGs and transposons were more likely to be relevant in infant gut microbiota compared to adult gut microbiota, and nine robust association rules were identified, among which Klebsiella pneumoniae, Enterobacter hormaechei_A, and Escherichia coli_D played important roles in this association phenomenon. The emphasis of this study is to investigate the synergistic transfer of specific ARGs and specific transposons in the infant gut microbiota, which can contribute to the study of microbial pathogenesis and the ARG dissemination dynamics.IMPORTANCEThe transfer of transposons carrying antibiotic resistance genes (ARGs) among microorganisms accelerates antibiotic resistance dissemination among infant gut microbiota. Nonetheless, it is unclear what the relationship between specific ARGs and specific transposons within the infant gut microbiota. K. pneumoniae, E. hormaechei_A, and E. coli_D were identified as key players in the nine robust association rules we discovered. Meanwhile, we found that infant gut microorganisms were more susceptible to horizontal gene transfer events about specific ARGs and specific transposons than adult gut microorganisms. These discoveries could enhance the understanding of microbial pathogenesis and the ARG dissemination dynamics within the infant gut microbiota.

Emergence and genomic characteristics of multi-drug-resistant Salmonella in pet turtles and children with diarrhoea

Pet turtles are a well-recognized source of human salmonellosis, posing a threat to human health, particularly children who commonly keep pet turtles. To date, the genomic characteristics of Salmonella among pet turtles and children has not been well described. We investigated the prevalence, antimicrobial resistance (AMR) and genomic characteristics of Salmonella from pet turtles in Beijing, China. In total, 9.6 % (46/480) of pet turtles were positive for Salmonella with S. Thompson being the dominant serovar (19/46) in 2019. Moreover, 80.4 % of Salmonella were multi-drug resistant (MDR) and 60.7 % were resistant to ampicillin, streptomycin, sulfonamides and tetracycline (ASSuT). We further compared the genomes of S. Thompson isolates from pet turtles (n=19) with those from children with diarrhoea (n=28) in the same region and year, most of which were sequence type (ST)26, with one novel ST7937 identified from a child-associated isolate. S. Thompson isolates from children with diarrhoea exhibited less resistance than isolates from pet turtles. Most MDR isolates possessed multiple AMR genes, including the AmpC β-lactamase-encoding genes bla DHA-15 and bla DHA-1 which co-occurred with the IncA/C and IncHI plasmid replicon types. To the best of our knowledge, this is the first time that the bla DHA-15 gene has been detected from Salmonella. Several pet turtle-associated S. Thompson isolates comprised phylogenetically close clusters with those from children with diarrhoea (<20 SNP differences). Bayesian analysis demonstrated that the Chinese ST26 S. Thompson strains had a recent evolutionary history and evolved into two major clades, with one clade acquiring various resistant plasmids. Our findings revealed the emergence of MDR Salmonella among pet turtles in China and provided evidence for the interspecies transmission of S. Thompson.

Carbapenem-resistant Klebsiella pneumoniae: the role of plasmids in emergence, dissemination, and evolution of a major clinical challenge

INTRODUCTION

Klebsiella pneumoniae is a major agent of healthcare-associated infections and a cause of some community-acquired infections, including severe bacteremic infections associated with metastatic abscesses in liver and other organs. Clinical relevance is compounded by its outstanding propensity to evolve antibiotic resistance. In particular, the emergence and dissemination of carbapenem resistance in K. pneumoniae has posed a major challenge due to the few residual treatment options, which have only recently been expanded by some new agents. The epidemiological success of carbapenem-resistant K. pneumoniae (CR-Kp) is mainly linked with clonal lineages that produce carbapenem-hydrolyzing enzymes (carbapenemases) encoded by plasmids.

AREAS COVERED

Here, we provide an updated overview on the mechanisms underlying the emergence and dissemination of CR-Kp, focusing on the role that plasmids have played in this phenomenon and in the co-evolution of resistance and virulence in K. pneumoniae.

EXPERT OPINION

CR-Kp have disseminated on a global scale, representing one of the most important contemporary public health issues. These strains are almost invariably associated with complex multi-drug resistance (MDR) phenotypes, which can also include recently approved antibiotics. The heterogeneity of the molecular bases responsible for these phenotypes poses significant hurdles for therapeutic and diagnostic purposes.

Plasmids, a molecular cornerstone of antimicrobial resistance in the One Health era

Antimicrobial resistance (AMR) poses a substantial threat to human health. The widespread prevalence of AMR is, in part, due to the horizontal transfer of antibiotic resistance genes (ARGs), typically mediated by plasmids. Many of the plasmid-mediated resistance genes in pathogens originate from environmental, animal or human habitats. Despite evidence that plasmids mobilize ARGs between these habitats, we have a limited understanding of the ecological and evolutionary trajectories that facilitate the emergence of multidrug resistance (MDR) plasmids in clinical pathogens. One Health, a holistic framework, enables exploration of these knowledge gaps. In this Review, we provide an overview of how plasmids drive local and global AMR spread and link different habitats. We explore some of the emerging studies integrating an eco-evolutionary perspective, opening up a discussion about the factors that affect the ecology and evolution of plasmids in complex microbial communities. Specifically, we discuss how the emergence and persistence of MDR plasmids can be affected by varying selective conditions, spatial structure, environmental heterogeneity, temporal variation and coexistence with other members of the microbiome. These factors, along with others yet to be investigated, collectively determine the emergence and transfer of plasmid-mediated AMR within and between habitats at the local and global scale.

Mobilome-driven partitions of the resistome in Salmonella

IMPORTANCE

Antimicrobial resistance (AMR) has become a significant global challenge, with an estimated 10 million deaths annually by 2050. The emergence of AMR is mainly attributed to mobile genetic elements (MGEs or mobilomes), which accelerate wide dissemination among pathogens. The interaction between mobilomes and AMR genes (or resistomes) in Salmonella, a primary cause of diarrheal diseases that results in over 90 million cases annually, remains poorly understood. The available fragmented or incomplete genomes remain a significant limitation in investigating the relationship between AMR and MGEs. Here, we collected the most extensive closed Salmonella genomes (n = 1,817) from various sources across 58 countries. Notably, our results demonstrate that resistome transmission between Salmonella lineages follows a specific pattern of MGEs and is influenced by external drivers, including certain socioeconomic factors. Therefore, targeted interventions are urgently needed to mitigate the catastrophic consequences of Salmonella AMR.

Plasmid genomic epidemiology of blaKPC carbapenemase-producing Enterobacterales in Canada, 2010-2021

Carbapenems are considered last-resort antibiotics for the treatment of infections caused by multidrug-resistant Enterobacterales, but carbapenem resistance due to acquisition of carbapenemase genes is a growing threat that has been reported worldwide. Klebsiella pneumoniae carbapenemase (blaKPC) is the most common type of carbapenemase in Canada and elsewhere; it can hydrolyze penicillins, cephalosporins, aztreonam, and carbapenems and is frequently found on mobile plasmids in the Tn4401 transposon. This means that alongside clonal expansion, blaKPC can disseminate through plasmid- and transposon-mediated horizontal gene transfer. We applied whole genome sequencing to characterize the molecular epidemiology of 829 blaKPC carbapenemase-producing isolates collected by the Canadian Nosocomial Infection Surveillance Program from 2010 to 2021. Using a combination of short-read and long-read sequencing, we obtained 202 complete and circular blaKPC-encoding plasmids. Using MOB-suite, 10 major plasmid clusters were identified from this data set which represented 87% (175/202) of the Canadian blaKPC-encoding plasmids. We further estimated the genomic location of incomplete blaKPC-encoding contigs and predicted a plasmid cluster for 95% (603/635) of these. We identified different patterns of carbapenemase mobilization across Canada related to different plasmid clusters, including clonal transmission of IncF-type plasmids (108/829, 13%) in K. pneumoniae clonal complex 258 and novel repE(pEh60-7) plasmids (44/829, 5%) in Enterobacter hormaechei ST316, and horizontal transmission of IncL/M (142/829, 17%) and IncN-type plasmids (149/829, 18%) across multiple genera. Our findings highlight the diversity of blaKPC genomic loci and indicate that multiple, distinct plasmid clusters have contributed to blaKPC spread and persistence in Canada.

Distribution of blaCTX-M-gene variants in E. coli from different origins in Ecuador

The increasing abundance of extended spectrum (β-lactamase (ESBL) genes in E. coli, and other commensal and pathogenic bacteria, endangers the utility of third or more recent generation cephalosporins, which are major tools for fighting deadly infections. The role of domestic animals in the transmission of ESBL carrying bacteria has been recognized, especially in low- and middle-income countries, however the horizontal gene transfer of these genes is difficult to assess. Here we investigate blaCTX-M gene diversity (and flanking nucleotide sequences) in E. coli from chicken and humans, in an Ecuadorian rural community and from chickens in another location in Ecuador. The blaCTX-M associated sequences in isolates from humans and chickens in the same remote community showed greater similarity than those found in E. coli in a chicken industrial operation 200 km away. Our study may provide evidence of blaCTX-M transfer between chickens and humans in the community.

Diverse mobile genetic elements shaped the evolution of Streptomyces virulence

Mobile genetic elements can innovate bacteria with new traits. In plant pathogenic Streptomyces, frequent and recent acquisition of integrative and conjugative or mobilizable genetic elements is predicted to lead to the emergence of new lineages that gained the capacity to synthesize Thaxtomin, a phytotoxin neccesary for induction of common scab disease on tuber and root crops. Here, we identified components of the Streptomyces-potato pathosystem implicated in virulence and investigated them as a nested and interacting system to reevaluate evolutionary models. We sequenced and analysed genomes of 166 strains isolated from over six decades of sampling primarily from field-grown potatoes. Virulence genes were associated to multiple subtypes of genetic elements differing in mechanisms of transmission and evolutionary histories. Evidence is consistent with few ancient acquisition events followed by recurrent loss or swaps of elements carrying Thaxtomin A-associated genes. Subtypes of another genetic element implicated in virulence are more distributed across Streptomyces. However, neither the subtype classification of genetic elements containing virulence genes nor taxonomic identity was predictive of pathogenicity on potato. Last, findings suggested that phytopathogenic strains are generally endemic to potato fields and some lineages were established by historical spread and further dispersed by few recent transmission events. Results from a hierarchical and system-wide characterization refine our understanding by revealing multiple mechanisms that gene and bacterial dispersion have had on shaping the evolution of a Gram-positive pathogen in agricultural settings.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Intensive care unit sinks are persistently colonized with multidrug resistant bacteria and mobilizable, resistance-conferring plasmids

Contamination of hospital sinks with microbial pathogens presents a serious potential threat to patients, but our understanding of sink colonization dynamics is largely based on infection outbreaks. Here, we investigate the colonization patterns of multidrug-resistant organisms (MDROs) in intensive care unit sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. Using culture-based methods, we recovered 822 bacterial isolates representing 104 unique species and genomospecies. Genomic analyses revealed long-term colonization by Pseudomonas spp. and Serratia marcescens strains across multiple rooms. Nanopore sequencing uncovered examples of long-term persistence of resistance-conferring plasmids in unrelated hosts. These data indicate that antibiotic resistance (AR) in Pseudomonas spp. is maintained both by strain colonization and horizontal gene transfer (HGT), while HGT maintains AR within Acinetobacter spp. and Enterobacterales, independent of colonization. These results emphasize the importance of proactive, genomic-focused surveillance of built environments to mitigate MDRO spread. IMPORTANCE Hospital sinks are frequently linked to outbreaks of antibiotic-resistant bacteria. Here, we used whole-genome sequencing to track the long-term colonization patterns in intensive care unit (ICU) sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. We analyzed 822 bacterial genomes, representing over 100 different species. We identified long-term contamination by opportunistic pathogens, as well as transient appearance of other common pathogens. We found that bacteria recovered from the ICU had more antibiotic resistance genes (ARGs) in their genomes compared to matched community spaces. We also found that many of these ARGs are harbored on mobilizable plasmids, which were found shared in the genomes of unrelated bacteria. Overall, this study provides an in-depth view of contamination patterns for common nosocomial pathogens and identifies specific targets for surveillance.

Long-read sequencing reveals genomic diversity and associated plasmid movement of carbapenemase-producing bacteria in a UK hospital over 6 years

Healthcare-associated infections (HCAIs) affect the most vulnerable people in society and are increasingly difficult to treat in the face of mounting antimicrobial resistance (AMR). Routine surveillance represents an effective way of understanding the circulation and burden of bacterial resistance and transmission in hospital settings. Here, we used whole-genome sequencing (WGS) to retrospectively analyse carbapenemase-producing Gram-negative bacteria from a single hospital in the UK over 6 years (n=165). We found that the vast majority of isolates were either hospital-onset (HAI) or HCAI. Most carbapenemase-producing organisms were carriage isolates, with 71 % isolated from screening (rectal) swabs. Using WGS, we identified 15 species, the most common being Escherichia coli and Klebsiella pneumoniae. Only one significant clonal outbreak occurred during the study period and involved a sequence type (ST)78 K. pneumoniae carrying bla NDM-1 on an IncFIB/IncHI1B plasmid. Contextualization with public data revealed little evidence of this ST outside of the study hospital, warranting ongoing surveillance. Carbapenemase genes were found on plasmids in 86 % of isolates, the most common types being bla NDM- and bla OXA-type alleles. Using long-read sequencing, we determined that approximately 30 % of isolates with carbapenemase genes on plasmids had acquired them via horizontal transmission. Overall, a national framework to collate more contextual genomic data, particularly for plasmids and resistant bacteria in the community, is needed to better understand how carbapenemase genes are transmitted in the UK.

PanGraph: scalable bacterial pan-genome graph construction

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

Whole Genome Sequencing-Based Tracing of a 2022 Introduction and Outbreak of Xanthomonas hortorum pv. pelargonii

Globalization has made agricultural commodities more accessible, available, and affordable. However, their global movement increases the potential for invasion by pathogens and necessitates development and implementation of sensitive, rapid, and scalable surveillance methods. Here, we used 35 strains, isolated by multiple diagnostic laboratories, as a case study for using whole genome sequence data in a plant disease diagnostic setting. Twenty-seven of the strains were isolated in 2022 and identified as Xanthomonas hortorum pv. pelargonii. Eighteen of these strains originated from material sold by a plant breeding company that had notified clients following a release of infected geranium cuttings. Analyses of whole genome sequences revealed epidemiological links among the 27 strains from different growers that confirmed a common source of the outbreak and uncovered likely secondary spread events within facilities that housed plants originating from different plant breeding companies. Whole genome sequencing data were also analyzed to reveal how preparatory and analytical methods can impact conclusions on outbreaks of clonal pathogenic strains. The results demonstrate the potential power of using whole genome sequencing among a network of diagnostic labs and highlight how sharing such data can help shorten response times to mitigate outbreaks more expediently and precisely than standard methods.

Phenotypic and genotypic discrepancies for carbapenemase-producing Citrobacter freundii in multiple isolates from a single patient

BACKGROUND

Carbapenemase-producing gram-negative organisms continue to be a significant healthcare concern and a therapeutic challenge. Members of the genus Citrobacter have emerged as increasingly multidrug resistant and versatile healthcare-associated pathogens. In this study we investigated five KPC-producing Citrobacter freundii isolates, from the same patient, that presented unusual phenotypic characteristics including false susceptibility to carbapenems detection by culture-based methods.

METHODS

The isolates were tested for antimicrobial susceptibility using broth microdilution and disk diffusion. Production of serine carbapenemase was confirmed with the mCIM (modified carbapenem inactivation method) test. Genotypes were determined by PCR and whole genome sequencing analysis.

RESULTS

The five isolates were susceptible to meropenem by broth microdilution and presented varying colonial morphologies and levels of susceptibility to carbapenems by multiple phenotypic methods, despite being positive for carbapenemase production by mCIM and positive for bla_KPC by PCR. Whole genome sequence analysis showed that three of the five highly related isolates harbor an additional gene cassette, including bla_CARB-2, ant(2''), aadA2, dfrA19, catB3, cmlA1, mph(E), msr(E), and qnrA1. The presence of these genes explains the difference in phenotypes observed.

CONCLUSION

Failure to detect and completely eradicate the carbapenemase-producing C. freundii in the urine with ertapenem therapy, likely due to the presence of a heterogeneous population, resulted in the phenotypic and genotypic adaptations of the organism as it disseminated to the bloodstream and kidneys. The fact that carbapenemase-producing C. freundii can elude detection by phenotypic methods and can so easily acquire and transfer resistance gene cassettes is of concern.

Worldwide Dissemination of blaKPC Gene by Novel Mobilization Platforms in Pseudomonas aeruginosa: A Systematic Review

The dissemination of blaKPC-harboring Pseudomonas aeruginosa (KPC-Pa) is considered a serious public health problem. This study provides an overview of the epidemiology of these isolates to try to elucidate novel mobilization platforms that could contribute to their worldwide spread. A systematic review in PubMed and EMBASE was performed to find articles published up to June 2022. In addition, a search algorithm using NCBI databases was developed to identify sequences that contain possible mobilization platforms. After that, the sequences were filtered and pair-aligned to describe the blaKPC genetic environment. We found 691 KPC-Pa isolates belonging to 41 different sequence types and recovered from 14 countries. Although the blaKPC gene is still mobilized by the transposon Tn4401, the non-Tn4401 elements (NTEKPC) were the most frequent. Our analysis allowed us to identify 25 different NTEKPC, mainly belonging to the NTEKPC-I, and a new type (proposed as IVa) was also observed. This is the first systematic review that consolidates information about the behavior of the blaKPC acquisition in P. aeruginosa and the genetic platforms implied in its successful worldwide spread. Our results show high NTEKPC prevalence in P. aeruginosa and an accelerated dynamic of unrelated clones. All information collected in this review was used to build an interactive online map.

Dissemination of carbapenemase-producing Enterobacterales in Ireland from 2012 to 2017: a retrospective genomic surveillance study

The spread of carbapenemase-producing Enterobacterales (CPE) is of major public health concern. The transmission dynamics of CPE in hospitals, particularly at the national level, are not well understood. Here, we describe a retrospective nationwide genomic surveillance study of CPE in Ireland between 2012 and 2017. We sequenced 746 national surveillance CPE samples obtained between 2012 and 2017. After clustering the sequences, we used thresholds based on pairwise SNPs, and reported within-host diversity along with epidemiological data to infer recent putative transmissions. All clusters in circulating clones, derived from high-resolution phylogenies, of a species (Klebsiella pneumoniae, Escherichia coli, Klebsiella oxytoca, Enterobacter cloacae, Enterobacter hormaechei and Citrobacter freundii) were individually examined for evidence of transmission. Antimicrobial resistance trends over time were also assessed. We identified 352 putative transmission events in six species including widespread and frequent transmissions in three species. We detected putative outbreaks in 4/6 species with three hospitals experiencing prolonged outbreaks. The bla OXA-48 gene was the main cause of carbapenem resistance in Ireland in almost all species. An expansion in the number of sequence types carrying bla OXA-48 was an additional cause of the increasing prevalence of carbapenemase-producing K. pneumoniae and E. coli.

Distribution and Transmission of Colistin Resistance Genes mcr-1 and mcr-3 among Nontyphoidal Salmonella Isolates in China from 2011 to 2020

Mobile colistin resistance (mcr) genes are present mainly in plasmids and can disseminate clonally or horizontally via either plasmids or insertion sequences in different genomic locations among the Enterobacteriaceae. A nationwide large-scale study on mcr prevalence and transmission in nontyphoidal Salmonella isolates is still lacking. Here, we identified 140 mcr-positive Salmonella isolates out of 7,106 isolates from 29 provinces in China from 2011 to 2020. We aligned short reads to putative plasmids from long-read hybrid assemblies and predicted the plasmid backbones of non-long-read sequencing isolates to elucidate mcr transmission patterns. The mcr-1 and mcr-3 genes are transmitted on similar high-risk clones (sequence type 34 [ST34]) but through plasmids of various replicon types. Furthermore, the ban on colistin use in food animals can lead to a decrease in the mcr-positive Salmonella prevalence among diarrheal patients, related mainly to IncHI2A_IncHI2 plasmids. We provide a framework for plasmid data incorporation into genomic surveillance systems, contributing to a better understanding of mcr spread and transmission. IMPORTANCE Nontyphoidal Salmonella is one of four major causative agents of diarrheal diseases globally, with most cases of salmonellosis being mild. Antimicrobial treatments are required for cases of life-threatening infections, and colistin is one of the last-line antibiotics for the treatment of multidrug-resistant Salmonella infections. However, the efficacy of colistin has been compromised by the emergence of various mcr genes. To elucidate the transmission of mcr genes in Salmonella isolates, our study analyzed 7,106 Salmonella strains from 29 provinces in China from 2011 to 2020. The results showed that mcr genes are transmitted on similar high-risk clones (ST34) but through plasmids of various replicon types. In addition, our data illustrated that the ban on the use of colistin in food animals led to a significant decrease in mcr-positive isolates. Our findings offer an essential step toward a more comprehensive understanding of the spread and transmission of mcr genes.

The DNA Phosphorothioation Restriction-Modification System Influences the Antimicrobial Resistance of Pathogenic Bacteria

Bacterial defense barriers, such as DNA methylation-associated restriction-modification (R-M) and the CRISPR-Cas system, play an important role in bacterial antimicrobial resistance (AMR). Recently, a novel R-M system based on DNA phosphorothioate (PT) modification has been shown to be widespread in the kingdom of Bacteria as well as Archaea. However, the potential role of the PT R-M system in bacterial AMR remains unclear. In this study, we explored the role of PT R-Ms in AMR with a series of common clinical pathogenic bacteria. By analyzing the distribution of AMR genes related to mobile genetic elements (MGEs), it was shown that the presence of PT R-M effectively reduced the distribution of horizontal gene transfer (HGT)-derived AMR genes in the genome, even in the bacteria that did not tend to acquire AMR genes by HGT. In addition, unique gene variation analysis based on pangenome analysis and MGE prediction revealed that the presence of PT R-M could suppress HGT frequency. Thus, this is the first report showing that the PT R-M system has the potential to repress HGT-derived AMR gene acquisition by reducing the HGT frequency. IMPORTANCE In this study, we demonstrated the effect of DNA PT modification-based R-M systems on horizontal gene transfer of AMR genes in pathogenic bacteria. We show that there is no apparent association between the genetic background of the strains harboring PT R-Ms and the number of AMR genes or the kinds of gene families. The strains equipped with PT R-M harbor fewer plasmid-derived, prophage-derived, or integrating mobile genetic element (iMGE)-related AMR genes and have a lower HGT frequency, but the degree of inhibition varies among different bacteria. In addition, compared with Salmonella enterica and Escherichia coli, Klebsiella pneumoniae prefers to acquire MGE-derived AMR genes, and there is no coevolution between PT R-M clusters and bacterial core genes.

Association of antibiotics and heavy metal arsenic to horizontal gene transfer from multidrug-resistant clinical strains to antibiotic-sensitive environmental strains

The emergence of multidrug-resistant bacteria is currently posing a significant threat to global public health. By testing for resistance to different antibiotic classes, we discovered that the majority of clinical bacteria are multidrug-resistant. These clinical multidrug-resistant species have antibiotic resistance genes on their plasmids that can be horizontally transferred to various antibiotic susceptible environmental bacterial species, resulting in antibiotic-resistant transconjugates. Furthermore, we discovered that the presence of an optimal concentration of antibiotics or heavy metal (arsenic) facilitates horizontal gene transfer through the formation of transconjugants. Notably, the addition of a conjugation inhibitor (2-hexadecynoic acid, a synthetic fatty acid) completely blocked the formation of antibiotic or arsenic-induced transconjugants. We discovered a high level of arsenic in water from the Shukratal region, Uttarakhand, India, which corresponded to a high serum level of arsenic in clinically infected individuals from the Shukratal region compared to other locations in Uttarakhand. Importantly, bacterial strains isolated from infected people who drink water from the Shukratal region with high arsenic levels were found to be more antibiotic-resistant than strains isolated from other sites. We discovered that bacterial strains isolated from individuals with high serum arsenic levels are significantly more resistant to antibiotics than individuals with low serum arsenic levels within the Shurkratal. This research sheds light on imminent threats to global health in which improper clinical, industrial, and other waste disposal, increased antibiotic concentrations in the environment, and increased human interference can easily transform commensal and pathogenic bacteria found in environmental niches into life-threatening multidrug-resistant superbugs.

Transmission of Carbapenem-Resistant Klebsiella pneumoniae in US Hospitals

BACKGROUND

Carbapenem-resistant Klebsiella pneumoniae (CRKp) is the most prevalent carbapenem-resistant Enterobacterales in the United States. We evaluated CRKp clustering in patients in US hospitals.

METHODS

From April 2016 to August 2017, 350 patients with clonal group 258 CRKp were enrolled in the Consortium on Resistance Against Carbapenems in Klebsiella and other Enterobacteriaceae, a prospective, multicenter, cohort study. A maximum likelihood tree was constructed using RAxML. Static clusters shared ≤21 single-nucleotide polymorphisms (SNP) and a most recent common ancestor. Dynamic clusters incorporated SNP distance, culture timing, and rates of SNP accumulation and transmission using the R program TransCluster.

RESULTS

Most patients were admitted from home (n = 150, 43%) or long-term care facilities (n = 115, 33%). Urine (n = 149, 43%) was the most common isolation site. Overall, 55 static and 47 dynamics clusters were identified involving 210 of 350 (60%) and 194 of 350 (55%) patients, respectively. Approximately half of static clusters were identical to dynamic clusters. Static clusters consisted of 33 (60%) intrasystem and 22 (40%) intersystem clusters. Dynamic clusters consisted of 32 (68%) intrasystem and 15 (32%) intersystem clusters and had fewer SNP differences than static clusters (8 vs 9; P = .045; 95% confidence interval [CI]: -4 to 0). Dynamic intersystem clusters contained more patients than dynamic intrasystem clusters (median [interquartile range], 4 [2, 7] vs 2 [2, 2]; P = .007; 95% CI: -3 to 0).

CONCLUSIONS

Widespread intrasystem and intersystem transmission of CRKp was identified in hospitalized US patients. Use of different methods for assessing genetic similarity resulted in only minor differences in interpretation.

Conjugation of plasmid harboring bla NDM-1 in a clinical Providencia rettgeri strain through the formation of a fusion plasmid

Providencia rettgeri has recently gained increased importance owing to the New Delhi metallo-β-lactamase (NDM) and other β-lactamases produced by its clinical isolates. These enzymes reduce the efficiency of antimicrobial therapy. Herein, we reported the findings of whole-genome sequence analysis and a comprehensive pan-genome analysis performed on a multidrug-resistant P. rettgeri 18004577 clinical strain recovered from the urine of a hospitalized patient in Shandong, China, in 2018. Providencia rettgeri 18004577 was found to have a genome assembly size of 4.6 Mb with a G + C content of 41%; a circular plasmid p18004577_NDM of 273.3 Kb, harboring an accessory multidrug-resistant region; and a circular, stable IncT plasmid p18004577_Rts of 146.2 Kb. Additionally, various resistance genes were identified in its genome, including bla _NDM-1, bla _OXA-10, bla _PER-4, aph(3')-VI, ant(2'')-Ia, ant(3')-Ia, sul1, catB8, catA1, mph(E), and tet. Conjugation experiments and whole-genome sequencing revealed that the bla _NDM-1 gene could be transferred to the transconjugant via the formation of pJ18004577_NDM, a novel hybrid plasmid. Based on the genetic comparison, the main possible formation process for pJ18004577_NDM was the insertion of the [ΔISKox2-IS26-ΔISKox2]-aph(3')-VI-bla _NDM-1 translocatable unit module from p18004577_NDM into plasmid p18004577_Rts in the Russian doll insertion structure (ΔISKox2-IS26-ΔISKox2), which played a role similar to that of IS26 using the "copy-in" route in the mobilization of [aph(3')-VI]-bla _NDM-1. The array, multiplicity, and diversity of the resistance and virulence genes in this strain necessitate stringent infection control, antibiotic stewardship, and periodic resistance surveillance/monitoring policies to preempt further horizontal and vertical spread of the resistance genes. Roary analysis based on 30 P. rettgeri strains pan genome identified 415 core, 756 soft core, 5,744 shell, and 12,967 cloud genes, highlighting the "close" nature of P. rettgeri pan-genome. After a comprehensive pan-genome analysis, representative biological information was revealed that included phylogenetic distances, presence or absence of genes across the P. rettgeri bacteria clade, and functional distribution of proteins. Moreover, pan-genome analysis has been shown to be an effective approach to better understand P. rettgeri bacteria because it helps develop various tailored therapeutic strategies based on their biological similarities and differences.

Implication of different replicons in the spread of the VIM-1-encoding integron, In110, in Enterobacterales from Czech hospitals

Background

VIM metallo-β-lactamases are enzymes characterized by the ability to hydrolyze all β-lactams. Usually, bla _VIM-like genes are carried by class 1 integrons. In the Czech Republic, only sporadic cases of VIM-producing Enterobacterales have been reported in which those isolates carried the VIM-1 carbapenemase-encoding integron In110. However, during 2019-2020, an increased number was reported. Therefore, the aim of the current study was to characterize the genetic elements involved in the increased spread of bla _VIM genes.

Materials and methods

32 VIM-producing Enterobacterales collected between 2019 and 2020 were subjected to: antimicrobial susceptibility testing, integron analysis, and short reads sequencing. Based on the results, 19 isolates were selected as representative and sequenced using Sequel I platform.

Results

The 32 VIM-producing isolates exhibited variations in the MICs of carbapenems. Based on short-read data, 26 of the 32 sequenced isolates harbored the bla _VIM-1 allele while six isolates carried the bla _VIM-4 gene. The most prevalent was the In110 integron (n = 24) and two isolates carried the In4873 class 1 integron. The bla _VIM-4 allele was identified in class 1 integrons In1174 (n = 3), In416 (n = 1), In2143 (n = 1) and In2150. Long reads sequencing revealed that the bla _VIM was carried by: pKPC-CAV1193-like (n = 6), HI1 (pNDM-CIT; n = 4), HI2 (n = 3), FIB (pECLA; n = 2) and N (n = 1) incompatibility groups. Two bla _VIM-carrying plasmids could not be typed by the database, while another one was integrated into the chromosome.

Conclusion

We observed the spread of VIM-encoding integrons, mainly of In110, among Enterobacterales isolated from Czech hospitals, but also an increased number of novel elements underlining the ongoing evolution.

Evaluation of five commercial DNA extraction kits using Salmonella as a model for implementation of rapid Nanopore sequencing in routine diagnostic laboratories

Oxford Nanopore long-read sequencing offers advantages over Illumina short reads for the identification and characterization of bacterial pathogens for outbreak detection and surveillance activities within a diagnostic public health laboratory context. Compared to Illumina, Nanopore is more cost-effective for small batches, has a lower capital cost and has a faster turnaround time, in addition to the ability to assemble complete bacterial genomes. The quantity and quality of DNA required for Nanopore sequencing are greater than for Illumina, and the DNA extraction methods recommended for obtaining high-molecular-weight DNA are different from those typically used in diagnostic laboratories. Using a Salmonella isolate with a previously closed PacBio genome as a model Enterobacteriaceae organism, we evaluated the quantity, quality and fragmentation of five commercial DNA extraction kits. Nanopore sequencing performance was evaluated for the top three methods: Qiagen EZ1 DNA Tissue, Qiagen DNeasy Blood and Tissue, and a modified, in-house version of the MasterPure Complete DNA and RNA purification. To evaluate the effect of post-extraction DNA purification methods, we subjected extracted DNA from the three selected extraction methods to purification by AMPure beads or ethanol precipitation and compared these outputs with untreated DNA as a control. All methods are suitable for routine whole-genome sequencing (WGS), since all 60 replicates had very high genome recovery rates, with ≥98 % of the reference genome covered by mapped Nanopore reads. For 85 % of the replicates, assembly was able to produce a complete, circular chromosome using either Flye or Canu. In most cases, it is recommended to move directly from extraction to sequencing, as untreated DNA had the highest rates of genome closure regardless of extraction method. Using our evaluation criteria, the Qiagen DNeasy Blood and Tissue kit was found to be the best overall method due to its low cost, ability to scale from single tubes to 96-well plates, and high consistency in yield and sequencing performance.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor^® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

The effect of sequencing and assembly on the inference of horizontal gene transfer on chromosomal and plasmid phylogenies

The spread of antibiotic resistance genes on plasmids is a threat to human and animal health. Phylogenies of bacteria and their plasmids contain clues regarding the frequency of plasmid transfer events, as well as the co-evolution of plasmids and their hosts. However, whole genome sequencing data from diverse ecological or clinical bacterial samples are rarely used to study plasmid phylogenies and resistance gene transfer. This is partially due to the difficulty of extracting plasmids from short-read sequencing data. Here, we use both short- and long-read sequencing data of 24 clinical extended-spectrum [Formula: see text]-lactamase (ESBL)-producing Escherichia coli to estimate chromosomal and plasmid phylogenies. We compare the impact of different sequencing and assembly methodologies on these phylogenies and on the inference of horizontal gene transfer. We find that chromosomal phylogenies can be estimated robustly with all methods, whereas plasmid phylogenies have more variable topology and branch lengths across the methods used. Specifically, hybrid methods that use long reads to resolve short-read assemblies (HybridSPAdes and Unicycler) perform better than those that started from long reads during assembly graph generation (Canu). By contrast, the inference of plasmid and antibiotic resistance gene transfer using a parsimony-based criterion is mostly robust to the choice of sequencing and assembly method. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

ESBL plasmids in Klebsiella pneumoniae: diversity, transmission and contribution to infection burden in the hospital setting

Background

Resistance to third-generation cephalosporins, often mediated by extended-spectrum beta-lactamases (ESBLs), is a considerable issue in hospital-associated infections as few drugs remain for treatment. ESBL genes are often located on large plasmids that transfer horizontally between strains and species of Enterobacteriaceae and frequently confer resistance to additional drug classes. Whilst plasmid transmission is recognised to occur in the hospital setting, the frequency and impact of plasmid transmission on infection burden, compared to ESBL + strain transmission, is not well understood.

Methods

We sequenced the genomes of clinical and carriage isolates of Klebsiella pneumoniae species complex from a year-long hospital surveillance study to investigate ESBL burden and plasmid transmission in an Australian hospital. Long-term persistence of a key transmitted ESBL + plasmid was investigated via sequencing of ceftriaxone-resistant isolates during 4 years of follow-up, beginning 3 years after the initial study.

Results

We found 25 distinct ESBL plasmids. We identified one plasmid, which we called Plasmid A, that carried bla_CTX-M-15 in an IncF backbone similar to pKPN-307. Plasmid A was transmitted at least four times into different Klebsiella species/lineages and was responsible for half of all ESBL episodes during the initial 1-year study period. Three of the Plasmid A-positive strains persisted locally 3–6 years later, and Plasmid A was detected in two additional strain backgrounds. Overall Plasmid A accounted for 21% of ESBL + infections in the follow-up period.

Conclusions

Here, we systematically surveyed ESBL strain and plasmid transmission over 1 year in a single hospital network. Whilst ESBL plasmid transmission events were rare in this setting, they had a significant and sustained impact on the burden of ceftriaxone-resistant and multidrug-resistant infections. If onward transmission of Plasmid A-carrying strains could have been prevented, this may have reduced the number of opportunities for Plasmid A to transmit and create novel ESBL + strains, as well as reducing overall ESBL infection burden.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01103-0.

The Sixth Element: a 102-kb RepABC Plasmid of Xenologous Origin Modulates Chromosomal Gene Expression in Dinoroseobacter shibae

The model organism Dinoroseobacter shibae and many other marine Rhodobacterales (Roseobacteraceae, Alphaproteobacteria) are characterized by a multipartite genome organization. Here, we show that the original isolate (Dshi-6) contained six extrachromosomal replicons (ECRs), whereas the strain deposited at the DSMZ (Dshi-5) lacked a 102-kb plasmid. To determine the role of the sixth plasmid, we investigated the genomic and physiological differences between the two strains. Therefore, both genomes were (re)sequenced, and gene expression, growth, and substrate utilization were examined. For comparison, we included additional plasmid-cured strains in the analysis. In the Dshi-6 population, the conjugative 102-kb RepABC-9 plasmid was present in only about 50% of the cells, irrespective of its experimentally validated stability. In the presence of the sixth plasmid, copy number changes of other ECRs, in particular, a decrease of the 86-kb plasmid, were observed. The most conspicuous finding was the strong influence of plasmids on chromosomal gene expression, especially the repression of the CtrA regulon and the activation of the denitrification gene cluster. Expression is inversely controlled by either the presence of the 102-kb plasmid or the absence of the 86-kb plasmid. We identified regulatory genes on both plasmids, i.e., a sigma 70 factor and a quorum sensing synthase, that might be responsible for these major changes. The tremendous effects that were probably even underestimated challenge the current understanding of the relevance of volatile plasmids not only for the original host but also for new recipients after conjugation.

IMPORTANCE Plasmids are small DNA molecules that replicate independently of the bacterial chromosome. The common view of the role of plasmids is dominated by the accumulation of resistance genes, which is responsible for the antibiotic crisis in health care and livestock breeding. Beyond rapid adaptations to a changing environment, no general relevance for the host cell’s regulome was attributed to these volatile ECRs. The current study shows for the model organism D. shibae that its chromosomal gene expression is strongly influenced by two plasmids. We provide evidence that the gain or loss of plasmids not only results in minor alterations of the genetic repertoire but also can have tremendous effects on bacterial physiology. The central role of some plasmids in the regulatory network of the host could also explain their persistence despite fitness costs, which has been described as the “plasmid paradox.”

Contrasting long-term dynamics of antimicrobial resistance and virulence plasmids in Salmonella Typhimurium from animals

Plasmids are mobile elements that can carry genes encoding traits of clinical concern, including antimicrobial resistance (AMR) and virulence. Population-level studies of Enterobacterales, including Escherichia coli, Shigella and Klebsiella, indicate that plasmids are important drivers of lineage expansions and dissemination of AMR genes. Salmonella Typhimurium is the second most common cause of salmonellosis in humans and livestock in the UK and Europe. The long-term dynamics of plasmids between S. Typhimurium were investigated using isolates collected through national surveillance of animals in England and Wales over a 25-year period. The population structure of S. Typhimurium and its virulence plasmid (where present) were inferred through phylogenetic analyses using whole-genome sequence data for 496 isolates. Antimicrobial resistance genes and plasmid markers were detected in silico. Phenotypic plasmid characterization, using the Kado and Liu method, was used to confirm the number and size of plasmids. The differences in AMR and plasmids between clades were striking, with livestock clades more likely to carry one or more AMR plasmid and be multi-drug-resistant compared to clades associated with wildlife and companion animals. Multiple small non-AMR plasmids were distributed across clades. However, all hybrid AMR–virulence plasmids and most AMR plasmids were highly clade-associated and persisted over decades, with minimal evidence of horizontal transfer between clades. This contrasts with the role of plasmids in the short-term dissemination of AMR between diverse strains in other Enterobacterales in high-antimicrobial-use settings, with implications for predicting plasmid dissemination amongst S. Typhimurium.

Genome Analysis of Enterobacter asburiae and Lelliottia spp. Proliferating in Oligotrophic Drinking Water Reservoirs and Lakes

Surface waters are one of the main sources for drinking water production, and thus microbial contamination should be as minimal as possible. However, high concentrations of coliform bacteria were detected in reservoirs and lakes used for drinking water production during summer months due to autochthonous proliferation processes. Here, we present the genomic analyses of 17 strains of Enterobacter asburiae and Lelliottia spp. proliferating in reservoirs and lakes with special focus on the hygienic relevance, antibiotic resistance, and adaptations to the oligotrophic environments. The genomes contain neither genes for the type III secretion system nor cytotoxins or hemolysins, which are considered typical virulence factors. Examination of antibiotic resistance genes revealed mainly efflux pumps and β-lactamase class C (ampC) genes. Phenotypically, single isolates of Enterobacter asburiae showed resistance to fosfomycin and ceftazidime. The genome analyses further suggest adaptations to oligotrophic and changing environmental conditions in reservoirs and lakes, e.g., genes to cope with low nitrate and phosphate levels and the ability to utilize substances released by algae, like amino acids, chitin, alginate, rhamnose, and fucose. This leads to the hypothesis that the proliferation of the coliform bacteria could occur at the end of summer due to algae die-off. IMPORTANCE Certain strains of coliform bacteria have been shown to proliferate in the oligotrophic water of drinking water reservoirs and lakes, reaching values above 10⁴ per 100 mL. Such high concentrations challenge drinking water treatment, and occasionally the respective coliform bacteria have been detected in the treated drinking water. Thus, the question of their hygienic relevance is of high importance for water suppliers and authorities. Our genomic analyses suggest that the strains are not hygienically relevant, as typical virulence factors are absent and antibiotic resistance genes in the genomes most likely are of natural origin. Furthermore, their presence in the water is not related to fecal contamination. The proliferation in reservoirs and lakes during stable summer stratification is an autochthonic process of certain E. asburiae and Lelliottia strains that are well adapted to the surrounding oligotrophic environment.

Dynamics of Antimicrobial Resistance and Genomic Epidemiology of Multidrug-Resistant Salmonella enterica Serovar Indiana ST17 from 2006 to 2017 in China

The genetic features of foodborne Salmonella have changed in recent years as multidrug-resistant (MDR) strains have become prevalent among various serovars. The recent expansion of MDR Salmonella enterica serovar Indiana sequence type 17 (ST17) poses an increasing threat to global public health, as 24.3% (61/251) of S. Indiana isolates in this study exhibited resistance to three clinically important antimicrobial agents: fluoroquinolones (ciprofloxacin), extended-spectrum β-lactams (cephalosporin), and macrolides (azithromycin). Both the evolutionary histories and antimicrobial resistance (AMR) profiles of this serovar remain to be described. Bioinformatic analysis revealed multiple lineages have coexisted and spread throughout China. Specifically, emergence of a predominant lineage appears to be associated with accumulated various substitutions in the chromosomal quinolone resistance-determining regions (GyrA S83F D87N and ParC T57S S80R) (141 [56.2%]), as well as acquisition of an extended-spectrum β-lactamase (ESBL)-producing IncHI2 plasmid that has subsequently undergone extensive rearrangement and an IncX1 plasmid that contains mph(A), conferring resistance to azithromycin. Several other evolutionary events influencing the trajectory of this drug-resistant serovar were also identified, including sporadic acquisitions of bla_CTX-M-carrying plasmids, along with chromosomal integration of bla_CTX-M within subclusters. Most human isolates reside in clusters containing isolates from animals, mainly from chickens, indicating the close relationship of human isolates with those from food animals. These data demonstrate that MDR S. Indiana ST17 is already widespread and capable of acquiring resistance traits against the clinical important antimicrobial agents, suggesting it should be considered a high-risk global MDR pathogen. The complexity of its evolutionary history has implications for AMR surveillance, epidemiological analysis, and control of emerging clinical lineages.

IMPORTANCE The emergence and worldwide spread of AMR Salmonella constitute great public health concerns. S. enterica serovar Indiana is a typical MDR serovar characterized by sporadic reports. However, comprehensive population genomics studies have not been performed on this serovar. This study provides a detailed and comprehensive insight into the rapid evolution of AMR in this important Salmonella serovar in the past 15 years in eight provinces of China. We documented diverse contributory genetic processes, including stable chromosomal integrations of resistance genes, the persistence and evolution of mobile resistance elements within sublineages, and sporadic acquisition of different resistance determinants that occur at all genetic levels (genes, genetic contexts, plasmids, and host strains). There are different mechanisms of antimicrobial resistance in S. enterica serovar Indiana from those of other serovars. This study sheds light on the formation of MDR S. enterica serovar Indiana with chickens as its potential reservoirs and paves the way to curb its further expansion among food animals.

Evolution of horizontal transmission in antimicrobial resistance plasmids

Mobile genetic elements (MGEs) are one of the main vectors for the spread of antimicrobial resistance (AMR) across bacteria, due to their ability to move horizontally between bacterial lineages. Horizontal transmission of AMR can increase AMR prevalence at multiple scales, from increasing the prevalence of infections by resistant bacteria to pathogen epidemics and worldwide spread of AMR across species. Among MGEs, conjugative plasmids are the main contributors to the spread of AMR. This review discusses the selective pressures acting on MGEs and their hosts to promote or limit the horizontal transmission of MGEs, the mechanisms by which transmission rates can evolve, and their implications for limiting the spread of AMR, with a focus on AMR plasmids.

Carbapenemase producing Enterobacteriaceae: Environmental reservoirs as primary targets for control and prevention strategies

Carbapenemase-producing Enterobacteriaceae (CPE) have become one of the greatest public health challenges globally. In the past decade, antimicrobial resistance (AMR) was viewed as a clinical problem in many parts of the world; hence, the role and magnitude of the contribution of the environment were not well appreciated. This review article was done with online published articles extracted from different databases using search terms related to the work. Evidence has shown that there exists the presence of carbapenemase genes in the environment, consequently fuelling the dissemination with alarming consequences. CPE when acquired causes life-threatening infections in humans. The health and economic impact of these infections are numerous, including treatment failure due to limited therapeutic options which hamper the containment of infectious diseases, further contaminating the environment and worsening the public health challenge. It is a well-known fact that the rate of emergence of resistant genes has outpaced the production of new antimicrobial agents, so it is pertinent to institute effective environmental measures to combat the spread of AMR organisms before it will completely gain a foothold and take us back to 'the pre-antibiotic era'. Environmental sources and reservoirs of resistant genes should therefore be amongst the primary targets for the control and prevention of the spread of resistant genes in the environment. This calls for the effective implementation of the 'one health' strategy with stakeholders committed to the design and enforcement of environmental mitigation policies and guidelines.

Population structure of blaKPC-harbouring IncN plasmids at a New York City medical centre and evidence for multi-species horizontal transmission

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are highly concerning MDR pathogens. Horizontal transfer of broad-host-range IncN plasmids may contribute to the dissemination of the Klebsiella pneumoniae carbapenemase (KPC), spreading carbapenem resistance among unrelated bacteria. However, the population structure and genetic diversity of IncN plasmids has not been fully elucidated.

OBJECTIVES

We reconstructed blaKPC-harbouring IncN plasmid genomes to characterize shared gene content, structural variability, and putative horizontal transfer within and across patients and diverse bacterial clones.

METHODS

We performed short- and long-read sequencing and hybrid assembly on 45 CRE isolates with blaKPC-harbouring IncN plasmids. Eight serial isolates from two patients were included to assess intra-patient plasmid dynamics. Comparative genomic analysis was performed to assess structural and sequence similarity across plasmids. Within IncN sublineages defined by plasmid MLST and kmer-based clustering, phylogenetic analysis was used to identify closely related plasmids.

RESULTS

Comparative analysis of IncN plasmid genomes revealed substantial heterogeneity including large rearrangements in serial patient plasmids and differences in structure and content across plasmid clusters. Within plasmid sublineages, core genome content and resistance gene regions were largely conserved. Closely related plasmids (≤1 SNP) were found in highly diverse isolates, including ten pST6 plasmids found in eight bacterial clones from three different species.

CONCLUSIONS

Genomic analysis of blaKPC-harbouring IncN plasmids revealed the presence of several distinct sublineages as well as substantial host diversity within plasmid clusters suggestive of frequent mobilization. This study reveals complex plasmid dynamics within a single plasmid family, highlighting the challenge of tracking plasmid-mediated transmission of blaKPC in clinical settings.

Whole genome sequencing reveals hidden transmission of carbapenemase-producing Enterobacterales

Carbapenemase-producing Enterobacterales (CPE) infection control practices are based on the paradigm that detected carriers in the hospital transmit to other patients who stay in the same ward. The role of plasmid-mediated transmission at population level remains largely unknown. In this retrospective cohort study over 4.7 years involving all multi-disciplinary public hospitals in Singapore, we analysed 779 patients who acquired CPE (1215 CPE isolates) detected by clinical or surveillance cultures. 42.0% met putative clonal transmission criteria, 44.8% met putative plasmid-mediated transmission criteria and 13.2% were unlinked. Only putative clonal transmissions associated with direct ward contact decreased in the second half of the study. Both putative clonal and plasmid-mediated transmission associated with indirect (no temporal overlap in patients’ admission period) ward and hospital contact did not decrease during the study period. Indirect ward and hospital contact were identified as independent risk factors associated with clonal transmission. In conclusion, undetected CPE reservoirs continue to evade hospital infection prevention measures. New measures are needed to address plasmid-mediated transmission, which accounted for 50% of CPE dissemination.

Carbapenemase-producing Enterobacterales cause healthcare-associated infections but modes of transmission are not well understood. Here, the authors find evidence of transmission without direct patient contact, indicating presence of undetected environmental reservoirs, whilst half of the transmission events are likely due to plasmid-mediated transmission.

KPC-3-, GES-5-, and VIM-1-Producing Enterobacterales Isolated from Urban Ponds

Carbapenems are antibiotics of pivotal importance in human medicine, the efficacy of which is threatened by the increasing prevalence of carbapenem-resistant Enterobacterales (CRE). Urban ponds may be reservoirs of CRE, although this hypothesis has been poorly explored. We assessed the proportion of CRE in urban ponds over a one-year period and retrieved 23 isolates. These were submitted to BOX-PCR, PFGE, 16S rDNA sequencing, antibiotic susceptibility tests, detection of carbapenemase-encoding genes, and conjugation assays. Isolates were affiliated with Klebsiella (n = 1), Raoultella (n = 11), Citrobacter (n = 8), and Enterobacter (n = 3). Carbapenemase-encoding genes were detected in 21 isolates: bla_KPC (n = 20), bla_GES-5 (n = 6), and bla_VIM (n = 1), with 7 isolates carrying two carbapenemase genes. Clonal isolates were collected from different ponds and in different campaigns. Citrobacter F6, Raoultella N9, and Enterobacter N10 were predicted as pathogens from whole-genome sequence analysis, which also revealed the presence of several resistance genes and mobile genetic elements. We found that bla_KPC-3 was located on Tn4401b (Citrobacter F6 and Enterobacter N10) or Tn4401d (Raoultella N9). The former was part of an IncFIA-FII pBK30683-like plasmid. In addition, bla_GES-5 was in a class 3 integron, either chromosomal (Raoultella N9) or plasmidic (Enterobacter N10). Our findings confirmed the role of urban ponds as reservoirs and dispersal sites for CRE.

High plasmidome diversity of extended-spectrum beta-lactam-resistant Escherichia coli isolates collected during one year in one community hospital

Plasmid-encoded antibiotic resistance encompasses many classes of currently used antibiotics. In globally distributed Escherichia coli lineages plasmids, which spread via horizontal gene transfer, are responsible for the dissemination of genes encoding extended-spectrum β-lactamases (ESBL). In this study, we combined 2nd and 3rd generation sequencing techniques to reconstruct the plasmidome of overall 97 clinical ESBL-E. coli isolates. Our results highlight the enormous plasmid diversity in respect to size, replicon-type and genetic content. Furthermore, we emphasize the diverse plasmid distribution patterns among the clinical isolates and the high intra- and extracellular mobility potential of resistance conferring genes. While the majority of resistance conferring genes were located on large plasmids of known replicon type, small cryptic plasmids seem to be underestimated resistance gene vectors. Our results contribute to a better understanding of the dissemination of resistance-conferring genes through horizontal gene transfer as well as clonal spread.