COPLA, a taxonomic classifier of plasmids

Genomic insights into epidemic plasmids carrying bla CTX-M and mcr-1 genes in Escherichia coli from Lebanese broiler production

Background

In a previous nationwide survey in the Lebanese broiler production, multidrug-resistant CTX-M-producing E. coli were found to carry the mobile colistin resistance gene mcr-1.

Objectives

To investigate the mobile genetic supports responsible for the spread of these resistance genes among E. coli in healthy broilers in Lebanon.

Methods

Thirty-three bla CTX-M and mcr-1 positive E. coli of various sequence types from 17 broilers farms were subjected to conjugation assays. Long-read sequencing (Oxford Nanopore Technologies) and hybrid assembly were performed to determine complete plasmid sequences and their phylogenetic diversity.

Results

Twenty-nine conjugative IncFII plasmids harboured the extended-spectrum β-lactamase genes bla CTX-M-3 (n = 25) or bla CTX-M-55 (n = 4). Highly related IncF2:A-:B-/bla CTX-M-3 plasmids differing only through IS-mediated genetic rearrangements in antibiotic resistance gene clusters were found in genetically diverse E. coli strains isolated from distant farms. The mobile colistin resistance genes mcr-1.1 and mcr-1.26 were carried by IncX4 and IncI2 plasmids. Worryingly, in one isolate, the ISEcp1-bla CTX-M-55 transposable unit was found integrated in a mcr-1.26-carrying IncX4 plasmid. Beside expanded cephalosporins and colistin resistances, all E. coli isolates were multidrug-resistant with different additional resistances against aminoglycosides, (fluoro)quinolones, fosfomycin, phenicols, sulphonamides, tetracycline and trimethoprim.

Conclusions

Closely related blaCTX-M-3/55-borne IncF2:A-:B- plasmids harbouring variable MDR regions and mcr-1 carrying IncX4 plasmids are widely disseminated in the E. coli population of healthy broilers in Lebanon. Further surveillance programmes of antimicrobial resistance and interventions to reduce the abusive use of medically important antibiotics are necessary to limit the spread of resistances in food-producing animals in Lebanon.

Towards quantifying plasmid similarity

Plasmids are extrachromosomal replicons which can quickly spread resistance and virulence genes between clinical pathogens. From the tens of thousands of currently available plasmid sequences we know that overall plasmid diversity is structured, with related plasmids sharing a largely conserved 'backbone' of genes while being able to carry very different genetic cargo. Moreover, plasmid genomes can be structurally plastic and undergo frequent rearrangements. So, how can we quantify plasmid similarity? Answering this question requires practical efforts to sample natural variation as well as theoretical considerations of what defines a group of related plasmids. Here we consider the challenges of analysing and rationalising the current plasmid data deluge to define appropriate similarity thresholds.

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

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

Azithromycin-resistant mph(A)-positive Salmonella enterica serovar Typhi in the United States

OBJECTIVES

The United States Centers for Disease Control and Prevention (CDC) conducts active surveillance for typhoid fever cases caused by Salmonella enterica serovar Typhi (Typhi). Here we describe the characteristics of the first two cases of mph(A)-positive azithromycin-resistant Typhi identified through US surveillance.

METHODS

Isolates were submitted to public health laboratories, sequenced, and screened for antimicrobial resistance determinants and plasmids, as part of CDC PulseNet's routine genomic surveillance. Antimicrobial susceptibility testing and long-read sequencing were also performed. Basic case information (age, sex, travel, outcome) was collected through routine questionnaires; additional epidemiological data was requested through follow-up patient interviews.

RESULTS

The patients are related and both reported travel to India (overlapping travel dates) before illness onset. Both Typhi genomes belong to the GenoTyphi lineage 4.3.1.1 and carry the azithromycin-resistance gene mph(A) on a PTU-FE (IncFIA/FIB/FII) plasmid. These strains differ genetically from mph(A)-positive Typhi genomes recently reported from Pakistan, suggesting independent emergence of azithromycin resistance in India.

CONCLUSIONS

Cases of typhoid fever caused by Typhi strains resistant to all available oral treatment options are cause for concern and support the need for vaccination of travellers to Typhi endemic regions. US genomic surveillance serves as an important global sentinel for detection of strains with known and emerging antimicrobial resistance profiles, including strains from areas where routine surveillance is not conducted.

Genomic and metagenomic analysis reveals shared resistance genes and mobile genetic elements in E. coli and Klebsiella spp. isolated from hospital patients and hospital wastewater at intra- and inter-genus level

Antimicrobial resistance (AMR) is a global problem that gives serious cause for concern. Hospital wastewater (HWW) is an important link between the clinical setting and the natural environment, and an escape route for pathogens that cause hospital infections, including urinary tract infections (UTI). Bacteria of the genera Escherichia and Klebsiella are common etiological factors of UTI, especially in children, and they can cause short-term infections, as well as chronic conditions. ESBL-producing Escherichia and Klebsiella have also emerged as potential indicators for estimating the burden of antimicrobial resistance under environmental conditions and the spread of AMR between clinical settings and the natural environment. In this study, whole-genome sequencing and the nanopore technology were used to analyze the complete genomes of ESBL-producing E.coli and Klebsiella spp. and the HWW metagenome, and to characterize the mechanisms of AMR. The similarities and differences in the encoded mechanisms of AMR in clinical isolates (causing UTI) and environmental strains (isolated from HWW and the HWW metagenome) were analyzed. Special attention was paid to the genetic context and the mobility of antibiotic resistance genes (ARGs) to determine the common sources and potential transmission of these genes. The results of this study suggest that the spread of drug resistance from healthcare facilities via HWW is not limited to the direct transmission of resistant clonal lines that are typically found in the clinical setting, but it also involves the indirect transfer of mobile elements carrying ARGs between bacteria colonizing various environments. Hospital wastewater could offer a supportive environment for plasmid evolution through the insertion of new ARGs, including typical chromosomal regions. These results indicate that interlined environments (hospital patients - HWW) should be closely monitored to evaluate the potential transmission routes of drug resistance in bacteria.

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Plasmid genomic epidemiology of carbapenem-hydrolysing class D β-lactamase (CDHL)-producing Enterobacterales in Canada, 2010-2021

Carbapenems are last-resort antibiotics for treatment of infections caused by multidrug-resistant Enterobacterales, but carbapenem resistance is a rising global threat due to the acquisition of carbapenemase genes. Oxacillinase-48 (bla OXA-48)-type carbapenemases are increasing in abundance in Canada and elsewhere; these genes are frequently found on mobile genetic elements and are associated with specific transposons. This means that alongside clonal dissemination, bla OXA-48-type genes can spread through plasmid-mediated horizontal gene transfer. We applied whole genome sequencing to characterize 249 bla OXA-48-type-producing Enterobacterales isolates collected by the Canadian Nosocomial Infection Surveillance Program from 2010 to 2021. Using a combination of short- and long-read sequencing, we obtained 70 complete and circular bla OXA-48-type-encoding plasmids. Using MOB-suite, four major plasmids clustered were identified, and we further estimated a plasmid cluster for 91.9 % (147/160) of incomplete bla OXA-48-type-encoding contigs. We identified different patterns of carbapenemase mobilization across Canada, including horizontal transmission of bla OXA-181/IncX3 plasmids (75/249, 30.1 %) and bla OXA-48/IncL/M plasmids (47/249, 18.9 %), and both horizontal transmission and clonal transmission of bla OXA-232 for Klebsiella pneumoniae ST231 on ColE2-type/ColKP3 plasmids (25/249, 10.0 %). Our findings highlight the diversity of OXA-48-type plasmids and indicate that multiple plasmid clusters and clonal transmission have contributed to bla OXA-48-type spread and persistence in Canada.

The evolution of antibiotic resistance islands occurs within the framework of plasmid lineages

Bacterial pathogens carrying multidrug resistance (MDR) plasmids are a major threat to human health. The acquisition of antibiotic resistance genes (ARGs) in plasmids is often facilitated by mobile genetic elements that copy or translocate ARGs between DNA molecules. The agglomeration of mobile elements in plasmids generates resistance islands comprising multiple ARGs. However, whether the emergence of resistance islands is restricted to specific MDR plasmid lineages remains understudied. Here we show that the agglomeration of ARGs in resistance islands is biased towards specific large plasmid lineages. Analyzing 6784 plasmids in 2441 Escherichia, Salmonella, and Klebsiella isolates, we quantify that 84% of the ARGs in MDR plasmids are found in resistance islands. We furthermore observe rapid evolution of ARG combinations in resistance islands. Most regions identified as resistance islands are shared among closely related plasmids but rarely among distantly related plasmids. Our results suggest the presence of barriers for the dissemination of ARGs between plasmid lineages, which are related to plasmid genetic properties, host range and the plasmid evolutionary history. The agglomeration of ARGs in plasmids is attributed to the workings of mobile genetic elements that operate within the framework of existing plasmid lineages.

Comparison of CTX-M encoding plasmids present during the early phase of the ESBL pandemic in western Sweden

Plasmids encoding blaCTX-M genes have greatly shaped the evolution of E. coli producing extended-spectrum beta-lactamases (ESBL-E. coli) and adds to the global threat of multiresistant bacteria by promoting horizontal gene transfer (HGT). Here we screened the similarity of 47 blaCTX-M -encoding plasmids, from 45 epidemiologically unrelated and disperse ESBL-E. coli strains, isolated during the early phase (2009-2014) of the ESBL pandemic in western Sweden. Using optical DNA mapping (ODM), both similar and rare plasmids were identified. As many as 57% of the plasmids formed five ODM-plasmid groups of at least three similar plasmids per group. The most prevalent type (28%, IncIl, pMLST37) encoded blaCTX-M-15 (n = 10), blaCTX-M-3 (n = 2) or blaCTX-M-55 (n = 1). It was found in isolates of various sequence types (STs), including ST131. This could indicate ongoing local HGT as whole-genome sequencing only revealed similarities with a rarely reported, IncIl plasmid. The second most prevalent type (IncFII/FIA/FIB, F1:A2:B20) harboring blaCTX-M-27, was detected in ST131-C1-M27 isolates, and was similar to plasmids previously reported for this subclade. The results also highlight the need for local surveillance of plasmids and the importance of temporospatial epidemiological links so that detection of a prevalent plasmid is not overestimated as a potential plasmid transmission event in outbreak investigations.

Antibiotic resistance begets more resistance: chromosomal resistance mutations mitigate fitness costs conferred by multi-resistant clinical plasmids

Plasmids are the primary vectors of horizontal transfer of antibiotic resistance genes among bacteria. Previous studies have shown that the spread and maintenance of plasmids among bacterial populations depend on the genetic makeup of both the plasmid and the host bacterium. Antibiotic resistance can also be acquired through mutations in the bacterial chromosome, which not only confer resistance but also result in changes in bacterial physiology and typically a reduction in fitness. However, it is unclear whether chromosomal resistance mutations affect the interaction between plasmids and the host bacteria. To address this question, we introduced 13 clinical plasmids into a susceptible Escherichia coli strain and three different congenic mutants that were resistant to nitrofurantoin (ΔnfsAB), ciprofloxacin (gyrA, S83L), and streptomycin (rpsL, K42N) and determined how the plasmids affected the exponential growth rates of the host in glucose minimal media. We find that though plasmids confer costs on the susceptible strains, those costs are fully mitigated in the three resistant mutants. In several cases, this results in a competitive advantage of the resistant strains over the susceptible strain when both carry the same plasmid and are grown in the absence of antibiotics. Our results suggest that bacteria carrying chromosomal mutations for antibiotic resistance could be a better reservoir for resistance plasmids, thereby driving the evolution of multi-drug resistance.IMPORTANCEPlasmids have led to the rampant spread of antibiotic resistance genes globally. Plasmids often carry antibiotic resistance genes and other genes needed for its maintenance and spread, which typically confer a fitness cost on the host cell observed as a reduced growth rate. Resistance is also acquired via chromosomal mutations, and similar to plasmids they also reduce bacterial fitness. However, we do not know whether resistance mutations affect the bacterial ability to carry plasmids. Here, we introduced 13 multi-resistant clinical plasmids into a susceptible and three different resistant E. coli strains and found that most of these plasmids do confer fitness cost on susceptible cells, but these costs disappear in the resistant strains which often lead to fitness advantage for the resistant strains in the absence of antibiotic selection. Our results imply that already resistant bacteria are a more favorable reservoir for multi-resistant plasmids, promoting the ascendance of multi-resistant bacteria.

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

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

Phage-plasmids promote recombination and emergence of phages and plasmids

Phages and plasmids are regarded as distinct types of mobile genetic elements that drive bacterial evolution by horizontal gene transfer. However, the distinction between both types is blurred by the existence of elements known as prophage-plasmids or phage-plasmids, which transfer horizontally between cells as viruses and vertically within cellular lineages as plasmids. Here, we study gene flow between the three types of elements. We show that the gene repertoire of phage-plasmids overlaps with those of phages and plasmids. By tracking recent recombination events, we find that phage-plasmids exchange genes more frequently with plasmids than with phages, and that direct gene exchange between plasmids and phages is less frequent in comparison. The results suggest that phage-plasmids can mediate gene flow between plasmids and phages, including exchange of mobile element core functions, defense systems, and antibiotic resistance. Moreover, a combination of gene transfer and gene inactivation may result in the conversion of elements. For example, gene loss turns P1-like phage-plasmids into integrative prophages or into plasmids (that are no longer phages). Remarkably, some of the latter have acquired conjugation-related functions to became mobilisable by conjugation. Thus, our work indicates that phage-plasmids can play a key role in the transfer of genes across mobile elements within their hosts, and can act as intermediates in the conversion of one type of element into another.

Plasmid content of carbapenem resistant Acinetobacter baumannii isolates belonging to five International Clones collected from hospitals of Alexandria, Egypt

Multidrug resistant Acinetobacter baumannii is one of the most important nosocomial pathogens worldwide. During the last decades it has become a major threat for healthcare settings due to the high antibiotic resistance rates among these isolates. Many resistance determinants are coded by conjugative or mobilizable plasmids, facilitating their dissemination. The majority of plasmids harbored by Acinetobacter species are less than 20 Kb, however, high molecular weight elements are the most clinically relevant since they usually contain antibiotic resistance genes. The aim of this work was to describe, classify and determine the genetic content of plasmids harbored by carbapem resistant A. baumannii isolates belonging to predominant clonal lineages circulating in hospitals from Alexandria, Egypt. The isolates were subjected to S1-Pulsed Field Gel Electrophoresis experiments to identify high molecular weight plasmids. To further analyze the plasmid content and the genetic localization of the antibiotic resistance genes, isolates were sequenced by Illumina Miseq and MinION Mk1C and a hybrid assembly was performed using Unicycler v0.5.0. Plasmids were detected with MOBsuite 3.0.3 and Copla.py v.1.0 and mapped using the online software Proksee.ca. Replicase genes were further analyzed through a BLAST against the Acinetobacter Plasmid Typing database. Eleven plasmids ranging in size from 4.9 to 205.6 Kb were characterized and mapped. All isolates contained plasmids, and, in many cases, more than two elements were identified. Antimicrobial resistance genes such as bla OXA-23, bla GES-like, aph(3')-VI and qacEΔ1 were found in likely conjugative large plasmids; while virulence determinants such as septicolysin or TonB-dependent receptors were identified in plasmids of small size. Some of these resistance determinants were, in turn, located within transposons and class 1 integrons. Among the identified plasmids, the majority encoded proteins belonging to the Rep_3 family, but replicases of the RepPriCT_1 (Aci6) family were also identified. Plasmids are of high interest as antibiotic resistance control tools, since they may be used as genetic markers for antibiotic resistance and virulence, and also as targets for the development of compounds that can inhibit transfer processes.

Inter-plasmid transfer of antibiotic resistance genes accelerates antibiotic resistance in bacterial pathogens

Antimicrobial resistance is a major threat for public health. Plasmids play a critical role in the spread of antimicrobial resistance via horizontal gene transfer between bacterial species. However, it remains unclear how plasmids originally recruit and assemble various antibiotic resistance genes (ARGs). Here, we track ARG recruitment and assembly in clinically relevant plasmids by combining a systematic analysis of 2420 complete plasmid genomes and experimental validation. Results showed that ARG transfer across plasmids is prevalent, and 87% ARGs were observed to potentially transfer among various plasmids among 8229 plasmid-borne ARGs. Interestingly, recruitment and assembly of ARGs occur mostly among compatible plasmids within the same bacterial cell, with over 88% of ARG transfers occurring between compatible plasmids. Integron and insertion sequences drive the ongoing ARG acquisition by plasmids, especially in which IS26 facilitates 63.1% of ARG transfer events among plasmids. In vitro experiment validated the important role of IS26 involved in transferring gentamicin resistance gene aacC1 between compatible plasmids. Network analysis showed four beta-lactam genes (blaTEM-1, blaNDM-4, blaKPC-2, and blaSHV-1) shuffling among 1029 plasmids and 45 clinical pathogens, suggesting that clinically alarming ARGs transferred accelerate the propagation of antibiotic resistance in clinical pathogens. ARGs in plasmids are also able to transmit across clinical and environmental boundaries, in terms of the high-sequence similarities of plasmid-borne ARGs between clinical and environmental plasmids. This study demonstrated that inter-plasmid ARG transfer is a universal mechanism for plasmid to recruit various ARGs, thus advancing our understanding of the emergence of multidrug-resistant plasmids.

IMG/PR: a database of plasmids from genomes and metagenomes with rich annotations and metadata

Plasmids are mobile genetic elements found in many clades of Archaea and Bacteria. They drive horizontal gene transfer, impacting ecological and evolutionary processes within microbial communities, and hold substantial importance in human health and biotechnology. To support plasmid research and provide scientists with data of an unprecedented diversity of plasmid sequences, we introduce the IMG/PR database, a new resource encompassing 699 973 plasmid sequences derived from genomes, metagenomes and metatranscriptomes. IMG/PR is the first database to provide data of plasmid that were systematically identified from diverse microbiome samples. IMG/PR plasmids are associated with rich metadata that includes geographical and ecosystem information, host taxonomy, similarity to other plasmids, functional annotation, presence of genes involved in conjugation and antibiotic resistance. The database offers diverse methods for exploring its extensive plasmid collection, enabling users to navigate plasmids through metadata-centric queries, plasmid comparisons and BLAST searches. The web interface for IMG/PR is accessible at https://img.jgi.doe.gov/pr. Plasmid metadata and sequences can be downloaded from https://genome.jgi.doe.gov/portal/IMG_PR.

MOBFinder: a tool for mobilization typing of plasmid metagenomic fragments based on a language model

BACKGROUND

Mobilization typing (MOB) is a classification scheme for plasmid genomes based on their relaxase gene. The host ranges of plasmids of different MOB categories are diverse, and MOB is crucial for investigating plasmid mobilization, especially the transmission of resistance genes and virulence factors. However, MOB typing of plasmid metagenomic data is challenging due to the highly fragmented characteristics of metagenomic contigs.

RESULTS

We developed MOBFinder, an 11-class classifier, for categorizing plasmid fragments into 10 MOB types and a nonmobilizable category. We first performed MOB typing to classify complete plasmid genomes according to relaxase information and then constructed an artificial benchmark dataset of plasmid metagenomic fragments (PMFs) from those complete plasmid genomes whose MOB types are well annotated. Next, based on natural language models, we used word vectors to characterize the PMFs. Several random forest classification models were trained and integrated to predict fragments of different lengths. Evaluating the tool using the benchmark dataset, we found that MOBFinder outperforms previous tools such as MOBscan and MOB-suite, with an overall accuracy approximately 59% higher than that of MOB-suite. Moreover, the balanced accuracy, harmonic mean, and F1-score reached up to 99% for some MOB types. When applied to a cohort of patients with type 2 diabetes (T2D), MOBFinder offered insights suggesting that the MOBF type plasmid, which is widely present in Escherichia and Klebsiella, and the MOBQ type plasmid might accelerate antibiotic resistance transmission in patients with T2D.

CONCLUSIONS

To the best of our knowledge, MOBFinder is the first tool for MOB typing of PMFs. The tool is freely available at https://github.com/FengTaoSMU/MOBFinder.

The contribution of plasmids to trait diversity in a soil bacterium

Plasmids are so closely associated with pathogens and antibiotic resistance that their potential for conferring other traits is often overlooked. Few studies consider how the full suite of traits encoded by plasmids is related to a host's environmental adaptation, particularly for Gram-positive bacteria. To investigate the role that plasmid traits might play in microbial communities from natural ecosystems, we identified plasmids carried by isolates of Curtobacterium (phylum Actinomycetota) from a variety of soil environments. We found that plasmids were common, but not ubiquitous, in the genus and varied greatly in their size and genetic diversity. There was little evidence of phylogenetic conservation among Curtobacterium plasmids even for closely related bacterial strains within the same ecotype, indicating that horizontal transmission of plasmids is common. The plasmids carried a wide diversity of traits that were not a random subset of the host chromosome. Furthermore, the composition of these plasmid traits was associated with the environmental context of the host bacterium. Together, the results indicate that plasmids contribute substantially to the microdiversity of a soil bacterium and that this diversity may play a role in niche differentiation and a bacterium's adaptation to its local environment.

Microbial evolution through horizontal gene transfer by mobile genetic elements

Mobile genetic elements (MGEs) are crucial for horizontal gene transfer (HGT) in bacteria and facilitate their rapid evolution and adaptation. MGEs include plasmids, integrative and conjugative elements, transposons, insertion sequences and bacteriophages. Notably, the spread of antimicrobial resistance genes (ARGs), which poses a serious threat to public health, is primarily attributable to HGT through MGEs. This mini-review aims to provide an overview of the mechanisms by which MGEs mediate HGT in microbes. Specifically, the behaviour of conjugative plasmids in different environments and conditions was discussed, and recent methodologies for tracing the dynamics of MGEs were summarised. A comprehensive understanding of the mechanisms underlying HGT and the role of MGEs in bacterial evolution and adaptation is important to develop strategies to combat the spread of ARGs.

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.

Tracing the Evolutionary Pathways of Serogroup O78 Avian Pathogenic Escherichia coli

Avian pathogenic E. coli (APEC) causes severe economic losses in the poultry industry, and O78 serogroup APEC strains are prevalent in chickens. In this study, we aimed to understand the evolutionary pathways and relationships between O78 APEC and other E. coli strains. To trace these evolutionary pathways, we classified 3101 E. coli strains into 306 subgenotypes according to the numbers and types of single nucleotide polymorphisms (RST0 to RST63-1) relative to the consensus sequence (RST0) of the RNA polymerase beta subunit gene and performed network analysis. The E. coli strains showed four apparently different evolutionary pathways (I-1, I-2, I-3, and II). The thirty-two Korean O78 APEC strains tested in this study were classified into RST4-4 (45.2%), RST3-1 (32.3%), RST21-1 (12.9%), RST4-5 (3.2%), RST5-1 (3.2%), and RST12-6 (3.2%), and all RSTs except RST21-1 (I-2) may have evolved through the same evolutionary pathway (I-1). A comparative genomic study revealed the highest relatedness between O78 strains of the same RST in terms of genome sequence coverage/identity and the spacer sequences of CRISPRs. The early-appearing RST3-1 and RST4-4 prevalence among O78 APEC strains may reflect the early settlement of O78 E. coli in chickens, after which these bacteria accumulated virulence and antibiotic resistance genes to become APEC strains. The zoonotic risk of the conventional O78 APEC strains is low at present, but the appearance of genetically distinct and multiple virulence gene-bearing RST21-1 O78 APEC strains may alert us to a need to evaluate their virulence in chickens as well as their zoonotic risk.

Genotypic characterization of Streptococcus didelphis causative of fatal infection in white-eared opossums

Streptococcus didelphis was once reported as related to severe infections in opossums. Thus, we present the first comprehensive whole-genome characterization of clinical S. didelphis strains isolated from white-eared opossums (Didelphis albiventris). Long-read whole-genome sequencing was performed using the MinION platform, which allowed the prediction of several genomic features. We observed that S. didelphis genomes harbor a cluster for streptolysin biosynthesis and a conserved genomic island with genes involved in transcriptional regulation (arlR) and transmembrane transport (bcrA). Antimicrobial resistance genes for several drug classes were found, including beta-lactam, which is the main antimicrobial class used in Streptococcus spp. infections; however, no phenotypical resistance was observed. In addition, we predicted the presence of 33 virulence factors in the analyzed genomes. High phylogenetic similarity was observed between clinical and reference strains, yet no clonality was suggested. We also proposed dnaN, gki, pros, and xpt as housekeeping candidates to be used in S. didelphis sequence typing. This is the first whole-genome characterization of S. didelphis, whose data provide important insights into its pathogenicity.

Development of an antimicrobial resistance plasmid transfer gene database for enteric bacteria

Introduction: Type IV secretion systems (T4SSs) are integral parts of the conjugation process in enteric bacteria. These secretion systems are encoded within the transfer (tra) regions of plasmids, including those that harbor antimicrobial resistance (AMR) genes. The conjugal transfer of resistance plasmids can lead to the dissemination of AMR among bacterial populations. Methods: To facilitate the analyses of the conjugation-associated genes, transfer related genes associated with key groups of AMR plasmids were identified, extracted from GenBank and used to generate a plasmid transfer gene dataset that is part of the Virulence and Plasmid Transfer Factor Database at FDA, serving as the foundation for computational tools for the comparison of the conjugal transfer genes. To assess the genetic feature of the transfer gene database, genes/proteins of the same name (e.g., traI/TraI) or predicted function (VirD4 ATPase homologs) were compared across the different plasmid types to assess sequence diversity. Two analyses tools, the Plasmid Transfer Factor Profile Assessment and Plasmid Transfer Factor Comparison tools, were developed to evaluate the transfer genes located on plasmids and to facilitate the comparison of plasmids from multiple sequence files. To assess the database and associated tools, plasmid, and whole genome sequencing (WGS) data were extracted from GenBank and previous WGS experiments in our lab and assessed using the analysis tools. Results: Overall, the plasmid transfer database and associated tools proved to be very useful for evaluating the different plasmid types, their association with T4SSs, and increased our understanding how conjugative plasmids contribute to the dissemination of AMR genes.

Machine Learning Suggests That Small Size Helps Broaden Plasmid Host Range

Plasmids mediate gene exchange across taxonomic barriers through conjugation, shaping bacterial evolution for billions of years. While plasmid mobility can be harnessed for genetic engineering and drug-delivery applications, rapid plasmid-mediated spread of resistance genes has rendered most clinical antibiotics useless. To solve this urgent and growing problem, we must understand how plasmids spread across bacterial communities. Here, we applied machine-learning models to identify features that are important for extending the plasmid host range. We assembled an up-to-date dataset of more than thirty thousand bacterial plasmids, separated them into 1125 clusters, and assigned each cluster a distribution possibility score, taking into account the host distribution of each taxonomic rank and the sampling bias of the existing sequencing data. Using this score and an optimized plasmid feature pool, we built a model stack consisting of DecisionTreeRegressor, EvoTreeRegressor, and LGBMRegressor as base models and LinearRegressor as a meta-learner. Our mathematical modeling revealed that sequence brevity is the most important determinant for plasmid spread, followed by P-loop NTPases, mobility factors, and β-lactamases. Ours and other recent results suggest that small plasmids may broaden their range by evading host defenses and using alternative modes of transfer instead of autonomous conjugation.

Strain belonging to an emerging, virulent sublineage of ST131 Escherichia coli isolated in fresh spinach, suggesting that ST131 may be transmissible through agricultural products

Food contamination with pathogenic Escherichia coli can cause severe disease. Here, we report the isolation of a multidrug resistant strain (A23EC) from fresh spinach. A23EC belongs to subclade C2 of ST131, a virulent clone of Extraintestinal Pathogenic E. coli (ExPEC). Most A23EC virulence factors are concentrated in three pathogenicity islands. These include PapGII, a fimbrial tip adhesin linked to increased virulence, and CsgA and CsgB, two adhesins known to facilitate spinach leaf colonization. A23EC also bears TnMB1860, a chromosomally-integrated transposon with the demonstrated potential to facilitate the evolution of carbapenem resistance among non-carbapenemase-producing enterobacterales. This transposon consists of two IS26-bound modular translocatable units (TUs). The first TU carries aac(6')-lb-cr, bla OXA-1, ΔcatB3, aac(3)-lle, and tmrB, and the second one harbors bla CXT-M-15. A23EC also bears a self-transmissible plasmid that can mediate conjugation at 20°C and that has a mosaic IncF [F(31,36):A(4,20):B1] and Col156 origin of replication. Comparing A23EC to 86 additional complete ST131 sequences, A23EC forms a monophyletic cluster with 17 other strains that share the following four genomic traits: (1) virotype E (papGII+); (2) presence of a PAI II536-like pathogenicity island with an additional cnf1 gene; (3) presence of chromosomal TnMB1860; and (4) frequent presence of an F(31,36):A(4,20):B1 plasmid. Sequences belonging to this cluster (which we named "C2b sublineage") are highly enriched in septicemia samples and their associated genetic markers align with recent reports of an emerging, virulent sublineage of the C2 subclade, suggesting significant pathogenic potential. This is the first report of a ST131 strain belonging to subclade C2 contaminating green leafy vegetables. The detection of this uropathogenic clone in fresh food is alarming. This work suggests that ST131 continues to evolve, gaining selective advantages and new routes of transmission. This highlights the pressing need for rigorous epidemiological surveillance of ExPEC in vegetables with One Health perspective.

Phenotypic and genomic characterization of Pseudomonas aeruginosa isolates recovered from catheter-associated urinary tract infections in an Egyptian hospital

Catheter-associated urinary tract infections (CAUTIs) represent one of the major healthcare-associated infections, and Pseudomonas aeruginosa is a common Gram-negative bacterium associated with catheter infections in Egyptian clinical settings. The present study describes the phenotypic and genotypic characteristics of 31 P. aeruginosa isolates recovered from CAUTIs in an Egyptian hospital over a 3 month period. Genomes of isolates were of good quality and were confirmed to be P. aeruginosa by comparison to the type strain (average nucleotide identity, phylogenetic analysis). Clonal diversity among the isolates was determined; eight different sequence types were found (STs 244, 357, 381, 621, 773, 1430, 1667 and 3765), of which ST357 and ST773 are considered to be high-risk clones. Antimicrobial resistance (AMR) testing according to European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines showed that the isolates were highly resistant to quinolones [ciprofloxacin (12/31, 38.7 %) and levofloxacin (9/31, 29 %) followed by tobramycin (10/31, 32.5 %)] and cephalosporins (7/31, 22.5 %). Genotypic analysis of resistance determinants predicted all isolates to encode a range of AMR genes, including those conferring resistance to aminoglycosides, β-lactamases, fluoroquinolones, fosfomycin, sulfonamides, tetracyclines and chloramphenicol. One isolate was found to carry a 422 938 bp pBT2436-like megaplasmid encoding OXA-520, the first report from Egypt of this emerging family of clinically important mobile genetic elements. All isolates were able to form biofilms and were predicted to encode virulence genes associated with adherence, antimicrobial activity, anti-phagocytosis, phospholipase enzymes, iron uptake, proteases, secretion systems and toxins. The present study shows how phenotypic analysis alongside genomic analysis may help us understand the AMR and virulence profiles of P. aeruginosa contributing to CAUTIs in Egypt.

Type II bacterial toxin-antitoxins: hypotheses, facts, and the newfound plethora of the PezAT system

Toxin-antitoxin (TA) systems are entities found in the prokaryotic genomes, with eight reported types. Type II, the best characterized, is comprised of two genes organized as an operon. Whereas toxins impair growth, the cognate antitoxin neutralizes its activity. TAs appeared to be involved in plasmid maintenance, persistence, virulence, and defence against bacteriophages. Most Type II toxins target the bacterial translational machinery. They seem to be antecessors of Higher Eukaryotes and Prokaryotes Nucleotide-binding (HEPN) RNases, minimal nucleotidyltransferase domains, or CRISPR-Cas systems. A total of four TAs encoded by Streptococcus pneumoniae, RelBE, YefMYoeB, Phd-Doc, and HicAB, belong to HEPN-RNases. The fifth is represented by PezAT/Epsilon-Zeta. PezT/Zeta toxins phosphorylate the peptidoglycan precursors, thereby blocking cell wall synthesis. We explore the body of knowledge (facts) and hypotheses procured for Type II TAs and analyse the data accumulated on the PezAT family. Bioinformatics analyses showed that homologues of PezT/Zeta toxin are abundantly distributed among 14 bacterial phyla mostly in Proteobacteria (48%), Firmicutes (27%), and Actinobacteria (18%), showing the widespread distribution of this TA. The pezAT locus was found to be mainly chromosomally encoded whereas its homologue, the tripartite omega-epsilon-zeta locus, was found mostly on plasmids. We found several orphan pezT/zeta toxins, unaccompanied by a cognate antitoxin.

Mge-cluster: a reference-free approach for typing bacterial plasmids

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

Reoccurring Escherichia coli O157:H7 Strain Linked to Leafy Greens-Associated Outbreaks, 2016-2019

Genomic characterization of an Escherichia coli O157:H7 strain linked to leafy greens-associated outbreaks dates its emergence to late 2015. One clade has notable accessory genomic content and a previously described mutation putatively associated with increased arsenic tolerance. This strain is a reoccurring, emerging, or persistent strain causing illness over an extended period.

Genetic rearrangements in Pseudomonas amygdali pathovar aesculi shape coronatine plasmids

Plant pathogenic Pseudomonas species use multiple classes of toxins and virulence factors during host infection. The genes encoding these pathogenicity factors are often located on plasmids and other mobile genetic elements, suggesting that they are acquired through horizontal gene transfer to confer an evolutionary advantage for successful adaptation to host infection. However, the genetic rearrangements that have led to mobilization of the pathogenicity genes are not fully understood. In this study, we have sequenced and analyzed the complete genome sequences of four Pseudomonas amygdali pv. aesculi (Pae), which infect European horse chestnut trees (Aesculus hippocastanum) and belong to phylogroup 3 of the P. syringae species complex. The four investigated genomes contain six groups of plasmids that all encode pathogenicity factors. Effector genes were found to be mostly associated with insertion sequence elements, suggesting that virulence genes are generally mobilized and potentially undergo horizontal gene transfer after transfer to a conjugative plasmid. We show that the biosynthetic gene cluster encoding the phytotoxin coronatine was recently transferred from a chromosomal location to a mobilizable plasmid that subsequently formed a co-integrate with a conjugative plasmid.

Characterization of the Environmental Plasmidome of the Red Sea

Plasmids contribute to microbial diversity and adaptation, providing microorganisms with the ability to thrive in a wide range of conditions in extreme environments. However, while the number of marine microbiome studies is constantly increasing, very little is known about marine plasmids, and they are very poorly represented in public databases. To extend the repertoire of environmental marine plasmids, we established a pipeline for the de novo assembly of plasmids in the marine environment by analyzing available microbiome metagenomic sequencing data. By applying the pipeline to data from the Red Sea, we identified 362 plasmid candidates. We showed that the distribution of plasmids corresponds to environmental conditions, particularly, depth, temperature, and physical location. At least 7 of the 362 candidates are most probably real plasmids, based on a functional analysis of their open reading frames (ORFs). Only one of the seven has been described previously. Three plasmids were identified in other public marine metagenomic data from different locations all over the world; these plasmids contained different cassettes of functional genes at each location. Analysis of antibiotic and metal resistance genes revealed that the same positions that were enriched with genes encoding resistance to antibiotics were also enriched with resistance to metals, suggesting that plasmids contribute site-dependent phenotypic modules to their ecological niches. Finally, half of the ORFs (50.8%) could not be assigned to a function, emphasizing the untapped potential of the unique marine plasmids to provide proteins with multiple novel functions. IMPORTANCE Marine plasmids are understudied and hence underrepresented in databases. Plasmid functional annotation and characterization is complicated but, if successful, may provide a pool of novel genes and unknown functions. Newly discovered plasmids and their functional repertoire are potentially valuable tools for predicting the dissemination of antimicrobial resistance, providing vectors for molecular cloning and an understanding of plasmid-bacterial interactions in various environments.

Genomic dissection of endemic carbapenem resistance reveals metallo-beta-lactamase dissemination through clonal, plasmid and integron transfer

Infections caused by metallo-beta-lactamase-producing organisms (MBLs) are a global health threat. Our understanding of transmission dynamics and how MBLs establish endemicity remains limited. We analysed two decades of blaIMP-4 evolution in a hospital using sequence data from 270 clinical and environmental isolates (including 169 completed genomes) and identified the blaIMP-4 gene across 7 Gram-negative genera, 68 bacterial strains and 7 distinct plasmid types. We showed how an initial multi-species outbreak of conserved IncC plasmids (95 genomes across 37 strains) allowed endemicity to be established through the ability of blaIMP-4 to disseminate in successful strain-genetic setting pairs we termed propagators, in particular Serratia marcescens and Enterobacter hormaechei. From this reservoir, blaIMP-4 persisted through diversification of genetic settings that resulted from transfer of blaIMP-4 plasmids between bacterial hosts and of the integron carrying blaIMP-4 between plasmids. Our findings provide a framework for understanding endemicity and spread of MBLs and may have broader applicability to other carbapenemase-producing organisms.

Co-Existence of blaNDM-1, blaOXA-23, blaOXA-64, blaPER-7 and blaADC-57 in a Clinical Isolate of Acinetobacter baumannii from Alexandria, Egypt

The increasing rates of antimicrobial resistance among carbapenem-resistant Acinetobacter baumannii in the Middle East and North Africa are one of the major concerns for healthcare settings. We characterised the first A. baumannii isolate harbouring five β-lactamases identified in Egypt. The isolate Ale25 was obtained from an ICU patient of a hospital from Alexandria. The isolate was phenotypically and genotypically screened for carbapenemase genes. The isolate was resistant to carbapenems, aminoglycosides, fluoroquinolones and cefiderocol. Whole-Genome Sequencing identified five β-lactamase genes, blaNDM-1, blaOXA-23, blaOXA-64, blaPER-7 and blaADC-57, together with other antibiotic resistance genes, conferring resistance to sulfonamides, macrolides, tetracyclines, rifamycin and chloramphenicol. Virulome analysis showed the presence of genes involved in adhesion and biofilm production, type II and VI secretion systems, exotoxins, etc. Multi-Locus Sequence Typing analysis identified the isolate as Sequence Types 113Pas and 2246Oxf, belonging to International Clone 7. Sequencing experiments revealed the presence of four plasmids of 2.7, 22.3, 70.4 and 240.8 Kb. All the β-lactamase genes were located in the chromosome, except the blaPER-7, gene which was found within the plasmid of 240.8 Kb. This study highlights the threat of the emergence and dissemination of these types of isolates.

A mathematician's guide to plasmids: an introduction to plasmid biology for modellers

Plasmids, extrachromosomal DNA molecules commonly found in bacterial and archaeal cells, play an important role in bacterial genetics and evolution. Our understanding of plasmid biology has been furthered greatly by the development of mathematical models, and there are many questions about plasmids that models would be useful in answering. In this review, we present an introductory, yet comprehensive, overview of the biology of plasmids suitable for modellers unfamiliar with plasmids who want to get up to speed and to begin working on plasmid-related models. In addition to reviewing the diversity of plasmids and the genes they carry, their key physiological functions, and interactions between plasmid and host, we also highlight selected plasmid topics that may be of particular interest to modellers and areas where there is a particular need for theoretical development. The world of plasmids holds a great variety of subjects that will interest mathematical biologists, and introducing new modellers to the subject will help to expand the existing body of plasmid theory.

Phylogenomics of Globally Spread Clonal Groups 14 and 15 of Klebsiella pneumoniae

Klebsiella pneumoniae sequence type 14 (ST14) and ST15 caused outbreaks of CTX-M-15 and/or carbapenemase producers worldwide, but their phylogeny and global dynamics remain unclear. We clarified the evolution of K. pneumoniae clonal group 14 (CG14) and CG15 by analyzing the capsular locus (KL), resistome, virulome, and plasmidome of public genomes (n = 481) and de novo sequences (n = 9) representing main sublineages circulating in Portugal. CG14 and CG15 evolved independently within 6 main subclades defined according to the KL and the accessory genome. The CG14 (n = 65) clade was structured in two large monophyletic subclades, CG14-I (KL2, 86%) and CG14-II (KL16, 14%), whose emergences were dated to 1932 and 1911, respectively. Genes encoding extended-spectrum β-lactamase (ESBL), AmpC, and/or carbapenemases were mostly observed in CG14-I (71% versus 22%). CG15 clade (n = 170) was segregated into subclades CG15-IA (KL19/KL106, 9%), CG15-IB (variable KL types, 6%), CG15-IIA (KL24, 43%) and CG15-IIB (KL112, 37%). Most CG15 genomes carried specific GyrA and ParC mutations and emerged from a common ancestor in 1989. CTX-M-15 was especially prevalent in CG15 (68% CG15 versus 38% CG14) and in CG15-IIB (92%). Plasmidome analysis revealed 27 predominant plasmid groups (PG), including particularly pervasive and recombinant F-type (n = 10), Col (n = 10), and new plasmid types. While blaCTX-M-15 was acquired multiple times by a high diversity of F-type mosaic plasmids, other antibiotic resistance genes (ARGs) were dispersed by IncL (blaOXA-48) or IncC (blaCMY/TEM-24) plasmids. We first demonstrate an independent evolutionary trajectory for CG15 and CG14 and how the acquisition of specific KL, quinolone-resistance determining region (QRDR) mutations (CG15), and ARGs in highly recombinant plasmids could have shaped the expansion and diversification of particular subclades (CG14-I and CG15-IIA/IIB). IMPORTANCE Klebsiella pneumoniae represents a major threat in the burden of antibiotic resistance (ABR). Available studies to explain the origin, the diversity, and the evolution of certain ABR K. pneumoniae populations have mainly been focused on a few clonal groups (CGs) using phylogenetic analysis of the core genome, the accessory genome being overlooked. Here, we provide unique insights into the phylogenetic evolution of CG14 and CG15, two poorly characterized CGs which have contributed to the global dissemination of genes responsible for resistance to first-line antibiotics such as β-lactams. Our results point out an independent evolution of these two CGs and highlight the existence of different subclades structured by the capsular type and the accessory genome. Moreover, the contribution of a turbulent flux of plasmids (especially multireplicon F type and Col) and adaptive traits (antibiotic resistance and metal tolerance genes) to the pangenome reflect the exposure and adaptation of K. pneumoniae under different selective pressures.

Lociq provides a loci-seeking approach for enhanced plasmid subtyping and structural characterization

Antimicrobial resistance (AMR) monitoring for public health is relying more on whole genome sequencing to characterize and compare resistant strains. This requires new approaches to describe and track AMR that take full advantage of the detailed data provided by genomic technologies. The plasmid-mediated transfer of AMR genes is a primary concern for AMR monitoring because plasmid rearrangement events can integrate new AMR genes into the plasmid backbone or promote hybridization of multiple plasmids. To better monitor plasmid evolution and dissemination, we developed the Lociq subtyping method to classify plasmids by variations in the sequence and arrangement of core plasmid genetic elements. Subtyping with Lociq provides an alpha-numeric nomenclature that can be used to denominate plasmid population diversity and characterize the relevant features of individual plasmids. Here we demonstrate how Lociq generates typing schema to track and characterize the origin, evolution and epidemiology of multidrug resistant plasmids.

Plasmid classifications

Plasmids are universally present in bacteria and play key roles in the dissemination of genes such as antibiotic resistance determinants. Major concepts in Plasmid Biology derive from the efforts to classify plasmids. Here, we review the main plasmid classification systems, starting by phenotype-based methods, such as fertility inhibition and incompatibility, followed by schemes based on a single gene (replicon type and MOB class), and finishing with recently developed approaches that use genetic distances between whole plasmid sequences. A comparison of the latter highlights significant differences between them. We further discuss the need for an operational definition of plasmid species that reveals their biological features, akin to plasmid taxonomic units (PTUs).

Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation

BACKGROUND

The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer.

RESULTS

SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation.

CONCLUSION

Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

Lytic bacteriophage vB_KmiS-Kmi2C disrupts biofilms formed by members of the Klebsiella oxytoca complex, and represents a novel virus family and genus

AIMS

This study aimed to characterize the lytic phage vB_KmiS-Kmi2C, isolated from sewage water on a GES-positive strain of Klebsiella michiganensis.

METHODS AND RESULTS

Comparative phylogenetic and network-based analyses were used to characterize the genome of phage vB_KmiS-Kmi2C (circular genome of 42 234 bp predicted to encode 55 genes), demonstrating it shared little similarity with other known phages. The phage was lytic on clinical strains of K. oxytoca (n = 2) and K. michiganensis (n = 4), and was found to both prevent biofilm formation and disrupt established biofilms produced by these strains.

CONCLUSIONS

We have identified a phage capable of killing clinically relevant members of the K. oxytoca complex (KoC). The phage represents a novel virus family (proposed name Dilsviridae) and genus (proposed name Dilsvirus).

Rapid emergence of extensively drug-resistant Shigella sonnei in France

Shigella sonnei, the main cause of bacillary dysentery in high-income countries, has become increasingly resistant to antibiotics. We monitored the antimicrobial susceptibility of 7121 S. sonnei isolates collected in France between 2005 and 2021. We detected a dramatic increase in the proportion of isolates simultaneously resistant to ciprofloxacin (CIP), third-generation cephalosporins (3GCs) and azithromycin (AZM) from 2015. Our genomic analysis of 164 such extensively drug-resistant (XDR) isolates identified 13 different clusters within CIP-resistant sublineage 3.6.1, which was selected in South Asia ∼15 years ago. AZM resistance was subsequently acquired, principally through IncFII (pKSR100-like) plasmids. The last step in the development of the XDR phenotype involved various extended-spectrum beta-lactamase genes (bla_CTX-M-3, bla_CTX-M-15, bla_CTX-M-27, bla_CTX-M-55, and bla_CTX-M-134) carried by different plasmids (IncFII, IncI1, IncB/O/K/Z) or even integrated into the chromosome, and encoding resistance to 3GCs. This rapid emergence of XDR S. sonnei, including an international epidemic strain, is alarming, and good laboratory-based surveillance of shigellosis will be crucial for informed decision-making and appropriate public health action.

High prevalence of GR2 and GR4 plasmids in Acinetobacter baumannii strains from Brazil

Acinetobacter baumannii is Gram-negative pathogen with extensive role in healthcare-associated infections (HAIs). Plasmids in this species are important carriers of antimicrobial resistance genes. In this work, we investigated the plasmids of 227 Brazilian A. baumannii genomes. A total of 389 plasmid sequences with 424 Rep proteins typed to 22 different homology groups (GRs) were identified. The GR2 plasmid group was the most predominant (40.6%), followed by the GR4 group (16.7%), representing ∼57% of all plasmids. There is a wide distribution of plasmids among the isolates and most strains carry more than one plasmid. Our analyses revealed a significant prevalence of GR4 plasmids in Brazilian A. baumannii genomes carrying several antimicrobial resistance genes, notably to carbapenem (39.43%). These plasmids harbor a MOBQ relaxase that might confer increased spreading potential in the environment. Most plasmids of the predominant groups belong to the same plasmid taxonomic unit (PTU-Pse7) and have a AbkA/AbkB toxin-antitoxin system that has a role in plasmid stability and dissemination of carbapenem resistance genes. The results of this work should contribute to our understanding of the molecular content of plasmids in a large and populous country, highlighting the importance of genomics for enhanced epidemiological surveillance.

Origins of transfer establish networks of functional dependencies for plasmid transfer by conjugation

Plasmids can be transferred between cells by conjugation, thereby driving bacterial evolution by horizontal gene transfer. Yet, we ignore the molecular mechanisms of transfer for many plasmids because they lack all protein-coding genes required for conjugation. We solved this conundrum by identifying hundreds of plasmids and chromosomes with conjugative origins of transfer in Escherichia coli and Staphylococcus aureus. These plasmids (pOriT) hijack the relaxases of conjugative or mobilizable elements, but not both. The functional dependencies between pOriT and other plasmids explain their co-occurrence: pOriT are abundant in cells with many plasmids, whereas conjugative plasmids are the most common in the others. We systematically characterized plasmid mobility in relation to conjugation and alternative mechanisms of transfer and can now propose a putative mechanism of transfer for ∼90% of them. In most cases, plasmid mobility seems to involve conjugation. Interestingly, the mechanisms of mobility are important determinants of plasmid-encoded accessory traits, since pOriTs have the highest densities of antimicrobial resistance genes, whereas plasmids lacking putative mechanisms of transfer have the lowest. We illuminate the evolutionary relationships between plasmids and suggest that many pOriT may have arisen by gene deletions in other types of plasmids. These results suggest that most plasmids can be transferred by conjugation.

Globally Disseminated Multidrug Resistance Plasmids Revealed by Complete Assembly of Multidrug Resistant Escherichia coli and Klebsiella pneumoniae Genomes from Diarrheal Disease in Botswana

Antimicrobial resistance is a disseminated global health challenge because many of the genes that cause resistance can transfer horizontally between bacteria. Despite the central role of extrachromosomal DNA elements called plasmids in driving the spread of resistance, the detection and surveillance of plasmids remains a significant barrier in molecular epidemiology. We assessed two DNA sequencing platforms alone and in combination for laboratory diagnostics in Botswana by annotating antibiotic resistance genes and plasmids in extensively drug resistant bacteria from diarrhea in Botswana. Long-read Nanopore DNA sequencing and high accuracy basecalling effectively estimated the architecture and gene content of three plasmids in Escherichia coli HUM3355 and two plasmids in Klebsiella pneumoniae HUM7199. Polishing the assemblies with Illumina reads increased base calling precision with small improvements to gene prediction. All five plasmids encoded one or more antibiotic resistance genes, usually within gene islands containing multiple antibiotic and metal resistance genes, and four plasmids encoded genes associated with conjugative transfer. Two plasmids were almost identical to antibiotic resistance plasmids sequenced in Europe and North America from human infection and a pig farm. These One Health connections demonstrate how low-, middle-, and high-income countries collectively benefit from increased whole genome sequencing capacity for surveillance and tracking of infectious diseases and antibiotic resistance genes that can transfer between animal hosts and move across continents.

Dissemination of metaldehyde catabolic pathways is driven by mobile genetic elements in Proteobacteria

Bioremediation of metaldehyde from drinking water using metaldehyde-degrading strains has recently emerged as a promising alternative. Whole-genome sequencing was used to obtain full genomes for metaldehyde degraders Acinetobacter calcoaceticus E1 and Sphingobium CMET-H. For the former, the genetic context of the metaldehyde-degrading genes had not been explored, while for the latter, none of the degrading genes themselves had been identified. In A. calcoaceticus E1, IS91 and IS6-family insertion sequences (ISs) were found surrounding the metaldehyde-degrading gene cluster located in plasmid pAME76. This cluster was located in closely-related plasmids and associated to identical ISs in most metaldehyde-degrading β- and γ-Proteobacteria, indicating horizontal gene transfer (HGT). For Sphingobium CMET-H, sequence analysis suggested a phytanoyl-CoA family oxygenase as a metaldehyde-degrading gene candidate due to its close homology to a previously identified metaldehyde-degrading gene known as mahX. Heterologous gene expression in Escherichia coli alongside degradation tests verified its functional significance and the degrading gene homolog was henceforth called mahS. It was found that mahS is hosted within the conjugative plasmid pSM1 and its genetic context suggested a crossover between the metaldehyde and acetoin degradation pathways. Here, specific replicons and ISs responsible for maintaining and dispersing metaldehyde-degrading genes in α, β and γ-Proteobacteria through HGT were identified and described. In addition, a homologous gene implicated in the first step of metaldehyde utilisation in an α-Proteobacteria was uncovered. Insights into specific steps of this possible degradation pathway are provided.

Emergence of a Novel Lineage and Wide Spread of a blaCTX-M-15/IncHI2/ST1 Plasmid among Nosocomial Enterobacter in Guadeloupe

Between April 2018 and August 2019, a total of 135 strains of Enterobacter cloacae complex (ECC) were randomly collected at the University Hospital Center of Guadeloupe to investigate the structure and diversity of the local bacterial population. These nosocomial isolates were initially identified genetically by the hsp60 typing method, which revealed the clinical relevance of E. xiangfangensis (n = 69). Overall, 57/94 of the third cephalosporin-resistant strains were characterized as extended-spectrum-β-lactamase (ESBL) producers, and their whole-genome was sequenced using Illumina technology to determine the clonal relatedness and diffusion of resistance genes. We found limited genetic diversity among sequence types (STs). ST114 (n = 13), ST1503 (n = 9), ST53 (n = 5) and ST113 (n = 4), which belong to three different Enterobacter species, were the most prevalent among the 57 ESBL producers. The blaCTXM-15 gene was the most prevalent ESBL determinant (56/57) and was in most cases associated with IncHI2/ST1 plasmid replicon carriage (36/57). To fully characterize this predominant blaCTXM-15/IncHI2/ST1 plasmid, four isolates from different lineages were also sequenced using Oxford Nanopore sequencing technology to generate long-reads. Hybrid sequence analyses confirmed the circulation of a well-conserved plasmid among ECC members. In addition, the novel ST1503 and its associated species (ECC taxon 4) were analyzed, in view of its high prevalence in nosocomial infections. These genetic observations confirmed the overall incidence of nosocomial ESBL Enterobacteriaceae infections acquired in this hospital during the study period, which was clearly higher in Guadeloupe (1.59/1000 hospitalization days) than in mainland France (0.52/1,000 hospitalization days). This project revealed issues and future challenges for the management and surveillance of nosocomial and multidrug-resistant Enterobacter in the Caribbean.

Three Distinct Annotation Platforms Differ in Detection of Antimicrobial Resistance Genes in Long-Read, Short-Read, and Hybrid Sequences Derived from Total Genomic DNA or from Purified Plasmid DNA

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

Genomic analysis of plasmid content in food isolates of E. coli strongly supports its role as a reservoir for the horizontal transfer of virulence and antibiotic resistance genes

The link between E. coli strains contaminating foods and human disease is unclear, with some reports supporting a direct transmission of pathogenic strains via food and others highlighting their role as reservoirs for resistance and virulence genes. Here we take a genomics approach, analyzing a large set of fully-assembled genomic sequences from E. coli available in GenBank. Most of the strains isolated in food are more closely related to each other than to clinical strains, arguing against a frequent direct transmission of pathogenic strains from food to the clinic. We also provide strong evidence of genetic exchanges between food and clinical strains that are facilitated by plasmids. This is based on an overlapped representation of virulence and resistance genes in plasmids isolated from these two sources. We identify clusters of phylogenetically-related plasmids that are largely responsible for the observed overlap and see evidence of specialization, with some food plasmid clusters preferentially transferring virulence factors over resistance genes. Consistent with these observations, food plasmids have a high mobilization potential based on their plasmid taxonomic unit classification and on an analysis of mobilization gene content. We report antibiotic resistance genes of high clinical relevance and their specific incompatibility group associations. Finally, we also report a striking enrichment for adhesins in food plasmids and their association with specific IncF replicon subtypes. The identification of food plasmids with specific markers (Inc and PTU combinations) as mediators of horizontal transfer between food and clinical strains opens new research avenues and should assist with the design of surveillance strategies.

Genome Sequences of 18 Salmonella enterica Serotype Hadar Strains Collected from Patients in the United States

Despite being linked to a number of recent poultry-associated outbreaks in the United States, few reference genomes are available for Salmonella enterica serotype Hadar. Here, we address this need by reporting 18 Salmonella Hadar genomes from samples collected from patients in the United States between 2014 and 2020.

Five Complete Salmonella enterica Serotype Reading Genomes Recovered from Patients in the United States

Between 2018 and 2019, Salmonella enterica serotype Reading caused a large, multistate outbreak linked to contact with raw turkey products in the United States. Here, we provide five Salmonella Reading reference genomes collected from US patients between 2016 and 2018.

Evolution of plasmid mobility: origin and fate of conjugative and non-conjugative plasmids

Conjugation drives the horizontal transfer of adaptive traits across prokaryotes. One-fourth of the plasmids encode the functions necessary to conjugate autonomously, the others being eventually mobilizable by conjugation. To understand the evolution of plasmid mobility, we studied plasmid size, gene repertoires, and conjugation-related genes. Plasmid gene repertoires were found to vary rapidly in relation to the evolutionary rate of relaxases, for example, most pairs of plasmids with 95% identical relaxases have fewer than 50% of homologs. Among 249 recent transitions of mobility type, we observed a clear excess of plasmids losing the capacity to conjugate. These transitions are associated with even greater changes in gene repertoires, possibly mediated by transposable elements, including pseudogenization of the conjugation locus, exchange of replicases reducing the problem of incompatibility, and extensive loss of other genes. At the microevolutionary scale of plasmid taxonomy, transitions of mobility type sometimes result in the creation of novel taxonomic units. Interestingly, most transitions from conjugative to mobilizable plasmids seem to be lost in the long term. This suggests a source-sink dynamic, where conjugative plasmids generate nonconjugative plasmids that tend to be poorly adapted and are frequently lost. Still, in some cases, these relaxases seem to have evolved to become efficient at plasmid mobilization in trans, possibly by hijacking multiple conjugative systems. This resulted in specialized relaxases of mobilizable plasmids. In conclusion, the evolution of plasmid mobility is frequent, shapes the patterns of gene flow in bacteria, the dynamics of gene repertoires, and the ecology of plasmids.

PlasmidHostFinder: Prediction of Plasmid Hosts Using Random Forest

Plasmids play a major role facilitating the spread of antimicrobial resistance between bacteria. Understanding the host range and dissemination trajectories of plasmids is critical for surveillance and prevention of antimicrobial resistance. Identification of plasmid host ranges could be improved using automated pattern detection methods compared to homology-based methods due to the diversity and genetic plasticity of plasmids. In this study, we developed a method for predicting the host range of plasmids using machine learning-specifically, random forests. We trained the models with 8,519 plasmids from 359 different bacterial species per taxonomic level; the models achieved Matthews correlation coefficients of 0.662 and 0.867 at the species and order levels, respectively. Our results suggest that despite the diverse nature and genetic plasticity of plasmids, our random forest model can accurately distinguish between plasmid hosts. This tool is available online through the Center for Genomic Epidemiology (https://cge.cbs.dtu.dk/services/PlasmidHostFinder/). IMPORTANCE Antimicrobial resistance is a global health threat to humans and animals, causing high mortality and morbidity while effectively ending decades of success in fighting against bacterial infections. Plasmids confer extra genetic capabilities to the host organisms through accessory genes that can encode antimicrobial resistance and virulence. In addition to lateral inheritance, plasmids can be transferred horizontally between bacterial taxa. Therefore, detection of the host range of plasmids is crucial for understanding and predicting the dissemination trajectories of extrachromosomal genes and bacterial evolution as well as taking effective countermeasures against antimicrobial resistance.