MOBscan: Automated Annotation of MOB Relaxases

OXA-48-like carbapenemases in Proteus mirabilis - novel genetic environments and a challenge for detection

OXA-48-like enzymes represent the most frequently detected carbapenemases in Enterobacterales in Western Europe, North Africa and the Middle East. In contrast to other species, the presence of OXA-48-like in Proteus mirabilis leads to an unusually susceptible phenotype with low MICs for carbapenems and piperacillin-tazobactam, which is easily missed in the diagnostic laboratory. So far, there is little data available on the genetic environments of the corresponding genes, blaOXA-48-like, in P. mirabilis. In this study susceptibility phenotypes and genomic data of 13 OXA-48-like-producing P. mirabilis were investigated (OXA-48, n = 9; OXA-181, n = 3; OXA-162, n = 1). Ten isolates were susceptible to meropenem and ertapenem and three isolates were susceptible to piperacillin-tazobactam. The gene blaOXA-48 was chromosomally located in 7/9 isolates. Thereof, in three isolates blaOXA-48 was inserted into a P. mirabilis genomic island. Of the three isolates harbouring blaOXA-181 one was located on an IncX3 plasmid and two were located on a novel MOBF plasmid, pOXA-P12, within the new transposon Tn7713. In 5/6 isolates with plasmidic location of blaOXA-48-like, the plasmids could conjugate to E. coli recipients in vitro. Vice versa, blaOXA-48-carrying plasmids could conjugate from other Enterobacterales into a P. mirabilis recipient. These data show a high diversity of blaOXA-48-like genetic environments compared to other Enterobacterales, where genetic environments are quite homogenous. Given the difficult-to-detect phenotype of OXA-48-like-producing P. mirabilis and the location of blaOXA-48-like on mobile genetic elements, it is likely that OXA-48-like-producing P. mirabilis can disseminate, escape most surveillance systems, and contribute to a hidden spread of OXA-48-like.

Nature-based bioreactors: Tackling antibiotic resistance in urban wastewater treatment

The overuse and misuse of antibiotics have accelerated the selection of antibiotic-resistant bacteria, significantly impacting human, animal, and environmental health. As aquatic environments are vulnerable to antibiotic resistance, suitable management practices should be adopted to tackle this phenomenon. Here we show an effective, nature-based solution for reducing antibiotic resistance from actual wastewater. We utilize a bioreactor that relies on benthic (biofilms) and planktonic microbial communities to treat secondary effluent from a small urban wastewater treatment plant (<10,000 population equivalent). This treated effluent is eventually released into the local aquatic ecosystem. We observe high removal efficiency for genes that provide resistance to commonly used antibiotic families, as well as for mobile genetic elements that could potentially aid in their spread. Importantly, we notice a buildup of sulfonamide (sul1 and sul2) and tetracycline (tet(C), tet(G), and tetR) resistance genes specifically in biofilms. This advancement marks the initial step in considering this bioreactor as a nature-based, cost-effective tertiary treatment option for small UWWTPs facing antibiotic resistance challenges.

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.

Camel Milk Resistome in Kuwait: Genotypic and Phenotypic Characterization

Antimicrobial resistance (AMR) is one of the major global health and economic threats. There is growing concern about the emergence of AMR in food and the possibility of transmission of microorganisms possessing antibiotic resistance genes (ARGs) to the human gut microbiome. Shotgun sequencing and in vitro antimicrobial susceptibility testing were used in this study to provide a detailed characterization of the antibiotic resistance profile of bacteria and their ARGs in dromedary camel milk. Eight pooled camel milk samples, representative of multiple camels distributed in the Kuwait desert, were collected from retail stores and analyzed. The genotypic analysis showed the presence of ARGs that mediate resistance to 18 classes of antibiotics in camel milk, with the highest resistance to fluoroquinolones (12.48%) and disinfecting agents and antiseptics (9%). Furthermore, the results pointed out the possible transmission of the ARGs to other bacteria through mobile genetic elements. The in vitro antimicrobial susceptibility testing indicated that 80% of the isolates were resistant to different classes of antibiotics, with the highest resistance observed against three antibiotic classes: penicillin, tetracyclines, and carbapenems. Multidrug-resistant pathogens including Klebsiella pneumoniae, Escherichia coli, and Enterobacter hormaechei were also revealed. These findings emphasize the human health risks related to the handling and consumption of raw camel milk and highlight the necessity of improving the hygienic practices of farms and retail stores to control the prevalence of ARGs and their transmission.

Complete nucleotide sequence and comparative genomic analysis of microcin B17 plasmid pMccB17

We present a comprehensive sequence and bioinformatic analysis of the prototypical microcin plasmid, pMccb17, which includes a definitive sequence for the microcin operon, mcb. Microcin B17 (MccB17) is a ribosomally synthesized and posttranslationally modified peptide produced by Escherichia coli. It inhibits bacterial DNA gyrase similarly to quinolone antibiotics. The mcb operon, which consists of seven genes encoding biosynthetic and immunity/export functions, was originally located on the low copy number IncFII plasmid pMccB17 in the Escherichia coli strain LP17. It was later transferred to E. coli K-12 through conjugation. In this study, the plasmid was extracted from the E. coli K-12 strain RYC1000 [pMccB17] and sequenced twice using an Illumina short-read method. The first sequencing was conducted with the host bacterial chromosome, and the plasmid DNA was then purified and sequenced separately. After assembly into a single contig, polymerase chain reaction primers were designed to close the single remaining gap via Sanger sequencing. The resulting complete circular DNA sequence is 69,190 bp long and includes 81 predicted genes. These genes were initially identified by Prokka and subsequently manually reannotated using BLAST. The plasmid was assigned to the F2:A-:B- replicon type with a MOBF12 group conjugation system. A comparison with other IncFII plasmids revealed a large proportion of shared genes, particularly in the conjugative plasmid backbone. However, unlike many contemporary IncFII plasmids, pMccB17 lacks transposable elements and antibiotic resistance genes. In addition to the mcb operon, this plasmid carries 25 genes of unknown function.

Transcriptional control of two distinct lactococcal plasmid-encoded conjugation systems

Lactococcal conjugative plasmids are poorly characterized compared to those harbored by numerous other Gram-positive bacteria, despite their significance in dairy fermentations and starter culture development. Furthermore, the transcriptional landscape of these lactococcal conjugation systems and their regulation have not been studied in any detail. Lactococcal plasmids pNP40 and pUC11B possess two genetically distinct and prevalent conjugation systems. Here, we describe the detailed transcriptional analysis of the pNP40 and pUC11B conjugation-associated gene clusters, revealing three and five promoters, respectively, for which the corresponding transcriptional start sites were identified. Regulation of several of these promoters, and therefore conjugation, is shown to involve the individual or concerted activities of the corresponding relaxase and transcriptional repressor(s) encoded by each conjugative plasmid. This work highlights how the conjugative potential of these systems may be unlocked, with significant implications for the starter culture and food fermentation industry.

Unraveling the diversity and dissemination dynamics of antimicrobial resistance genes in Enterobacteriaceae plasmids across diverse ecosystems

AIM

The objective of this study was to investigate the antimicrobial resistance genes (ARGs) in plasmids of Enterobacteriaceae from soil, sewage, and feces of food-producing animals and humans.

METHODS AND RESULTS

The plasmid sequences were obtained from the NCBI database. For the identification of ARG, comprehensive antibiotic resistance database (CARD), and ResFinder were used. Gene conservation and evolution were investigated using DnaSP v.6. The transfer potential of the plasmids was evaluated using oriTfinder and a MOB-based phylogenetic tree was reconstructed using Fastree. We identified a total of 1064 ARGs in all plasmids analyzed, conferring resistance to 15 groups of antibiotics, mostly aminoglycosides, beta-lactams, and sulfonamides. The greatest number of ARGs per plasmid was found in enterobacteria from chicken feces. Plasmids from Escherichia coli carrying multiple ARGs were found in all ecosystems. Some of the most abundant genes were shared among all ecosystems, including aph(6)-Id, aph(3'')-Ib, tet(A), and sul2. A high level of sequence conservation was found among these genes, and tet(A) and sul2 are under positive selective pressure. Approximately 62% of the plasmids carrying at least one ARG were potentially transferable. Phylogenetic analysis indicated a potential co-evolution of Enterobacteriaceae plasmids in nature.

CONCLUSION

The high abundance of Enterobacteriaceae plasmids from diverse ecosystems carrying ARGs reveals their widespread distribution and importance.

A comprehensive list of genes required for the efficient conjugation of plasmid Rts1 was determined by systematic deletion analysis

While conjugation-related genes have been identified in many plasmids by genome sequencing, functional analyses have not yet been performed in most cases, and a full set of conjugation genes has been identified for only a few plasmids. Rts1, a prototype IncT plasmid, is a conjugative plasmid that was originally isolated from Proteus vulgaris. Here, we conducted a systematic deletion analysis of Rts1 to fully understand its conjugation system. Through this analysis along with complementation assays, we identified 32 genes that are required for the efficient conjugation of Rts1 from Escherichia coli to E. coli. In addition, the functions of the 28 genes were determined or predicted; 21 were involved in mating-pair formation, three were involved in DNA transfer and replication, including a relaxase gene belonging to the MOBH12 family, one was involved in coupling, and three were involved in transcriptional regulation. Among the functionally well-analysed conjugation systems, most of the 28 genes showed the highest similarity to those of the SXT element, which is an integrative conjugative element of Vibrio cholerae. The Rts1 conjugation gene set included all 23 genes required for the SXT system. Two groups of plasmids with conjugation systems nearly identical or very similar to that of Rts1 were also identified.

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.

Structure and dispersion of the conjugative mobilome in surface ocean bacterioplankton

Mobile genetic elements (MGEs), collectively referred to as the "mobilome", can have a significant impact on the fitness of microbial communities and therefore on ecological processes. Marine MGEs have mainly been associated with wide geographical and phylogenetic dispersal of adaptative traits. However, whether the structure of this mobilome exhibits deterministic patterns in the natural community is still an open question. The aim of this study was to characterize the structure of the conjugative mobilome in the ocean surface bacterioplankton by searching the publicly available marine metagenomes from the TARA Oceans survey, together with molecular markers, such as relaxases and type IV coupling proteins of the type IV secretion system (T4SS). The T4SS machinery was retrieved in more abundance than relaxases in the surface marine bacterioplankton. Moreover, among the identified MGEs, mobilizable elements were the most abundant, outnumbering self-conjugative sequences. Detection of a high number of incomplete T4SSs provides insight into possible strategies related to trans-acting activity between MGEs, and accessory functions of the T4SS (e.g. protein secretion), allowing the host to maintain a lower metabolic burden in the highly dynamic marine system. Additionally, the results demonstrate a wide geographical dispersion of MGEs throughout oceanic regions, while the Southern Ocean appears segregated from other regions. The marine mobilome also showed a high similarity of functions present in known plasmid databases. Moreover, cargo genes were mostly related to DNA processing, but scarcely associated with antibiotic resistance. Finally, within the MGEs, integrative and conjugative elements showed wider marine geographic dispersion than plasmids.

Plasmid-encoded toxin defence mediates mutualistic microbial interactions

Gut environments harbour dense microbial ecosystems in which plasmids are widely distributed. Plasmids facilitate the exchange of genetic material among microorganisms while enabling the transfer of a diverse array of accessory functions. However, their precise impact on microbial community composition and function remains largely unexplored. Here we identify a prevalent bacterial toxin and a plasmid-encoded resistance mechanism that mediates the interaction between Lactobacilli and Enterococci. This plasmid is widespread across ecosystems, including the rumen and human gut microbiota. Biochemical characterization of the plasmid revealed a defence mechanism against reuterin, a toxin produced by various gut microbes, such as Limosilactobacillus reuteri. Using a targeted metabolomic approach, we find reuterin to be prevalent across rumen ecosystems with impacts on microbial community structure. Enterococcus strains carrying the protective plasmid were isolated and their interactions with L. reuteri, the toxin producer, were studied in vitro. Interestingly, we found that by conferring resistance against reuterin, the plasmid mediates metabolic exchange between the defending and the attacking microbial species, resulting in a beneficial relationship or mutualism. Hence, we reveal here an ecological role for a plasmid-coded defence system in mediating a beneficial interaction.

Plasmids of the incompatibility group FIBK occur in Klebsiella variicola from diverse ecological niches

Plasmids play a fundamental role in the evolution of bacteria by allowing them to adapt to different environments and acquire, through horizontal transfer, genes that confer resistance to different classes of antibiotics. Using the available in vitro and in silico plasmid typing systems, we analyzed a set of isolates and public genomes of K. variicola to study its plasmid diversity. The resistome, the plasmid multilocus sequence typing (pMLST), and molecular epidemiology using the MLST system were also studied. A high frequency of IncF plasmids from human isolates but lower frequency from plant isolates were found in our strain collection. In silico detection revealed 297 incompatibility (Inc) groups, but the IncFIBK (216/297) predominated in plasmids from human and environmental samples, followed by IncFIIK (89/297) and IncFIA/FIA(HI1) (75/297). These Inc groups were associated with clinically important ESBL (CTX-M-15), carbapenemases (KPC-2 and NDM-1), and colistin-resistant genes which were associated with major sequence types (ST): ST60, ST20, and ST10. In silico MOB typing showed 76% (311/404) of the genomes contained one or more of the six relaxase families with MOBF being most abundant. We identified untypeable plasmids carrying blaKPC-2, blaIMP-1, and blaSHV-187 but for which a relaxase was found; this may suggest that novel plasmid structures could be emerging in this bacterial species. The plasmid content in K. variicola has limited diversity, predominantly composed of IncFIBK plasmids dispersed in different STs. Plasmid detection using the replicon and MOB typing scheme provide a broader context of the plasmids in K. variicola. This study showed that whole-sequence-based typing provides current insights of the prevalence of plasmid types and their association with antimicrobial resistant genes in K. variicola obtained from humans and environmental niches.

Characterization of a Tigecycline-, Linezolid- and Vancomycin-Resistant Clinical Enteroccoccus faecium Isolate, Carrying vanA and vanB Genes

INTRODUCTION

Increasing incidence of Enterococcus faecium resistant to key antimicrobials used in therapy of hospitalized patients is a worrisome phenomenon observed worldwide. Our aim was to characterize a tigecycline-, linezolid- and vancomycin-resistant E. faecium isolate with the vanA and vanB genes, originating from a hematoma of a patient hospitalized in an intensive care unit in Poland.

METHODS

Antimicrobial susceptibility (a broad panel) was tested using gradient tests with predefined antibiotic concentrations. The complete genome sequence was obtained from a mixed assembly of Illumina MiSeq and Oxford Nanopore's MinION reads. The genome was analyzed with appropriate tools available at the Center for Genomic Epidemiology, PubMLST and GenBank. Transferability of oxazolidinone, tigecycline and vancomycin resistance genes was investigated by conjugation, followed by PCR screen of transconjugants for antimicrobial resistance genes and plasmid rep genes characteristic for the donor and genomic sequencing of selected transconjugants.

RESULTS

The isolate was resistant to most antimicrobials tested; susceptibility to daptomycin, erythromycin and chloramphenicol was significantly reduced, and only oritavancin retained the full activity. The isolate represented sequence type 18 (ST18) and carried vanA, vanB, poxtA, fexB, tet(L), tet(M), aac(6')-aph(2''), ant(6)-Ia and ant(6')-Ii. The vanA, poxtA and tet(M) genes located on ~ 40-kb plasmids were transferable by conjugation yielding transconjugants resistant to vancomycin, linezolid and tigecycline. The substitutions in LiaS, putative histidine kinase, SulP, putative sulfate transporter, RpoB and RpoC were potential determinants of an elevated daptomycin MIC. Comparative analyses of the studied isolate with E. faecium isolates from other countries revealed its similarity to ST18 isolates from Ireland and Uganda from human infections.

CONCLUSIONS

We provide the detailed characteristics of the genomic determinants of antimicrobial resistance of a clinical E. faecium demonstrating the concomitant presence of both vanA and vanB and resistance to vancomycin, linezolid, tigecycline and several other compounds and decreased daptomycin susceptibility. This isolate is a striking example of an accumulation of resistance determinants involving various mechanisms by a single hospital strain.

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

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

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.

Structure and functions of a multireplicon genome of Antarctic Psychrobacter sp. ANT_H3: characterization of the genetic modules suitable for the construction of the plasmid-vectors for cold-active bacteria

Among Psychrobacter spp., there are several multireplicon strains, carrying more than two plasmids. Psychrobacter sp. ANT_H3 carries as many as 11 extrachromosomal replicons, which is the highest number in Psychrobacter spp. Plasmids of this strain were subjected to detailed genomic analysis, which enables an insight into the structure and functioning of this multireplicon genome. The replication and conjugal transfer modules of ANT_H3 plasmids were analyzed functionally to discover their potential for being used as building blocks for the construction of novel plasmid-vectors for cold-active bacteria. It was shown that two plasmids have a narrow host range as they were not able to replicate in species other than Psychrobacter, while remaining plasmids had a wider host range and were functional in various Alpha- and Gammaproteobacteria. Moreover, it was confirmed that mobilization modules of seven plasmids were functional, i.e., could be mobilized for conjugal transfer by the RK2 conjugation system. Auxiliary genes were also distinguished in ANT_H3 plasmids, including these encoding putative DNA-protecting protein DprA, multidrug efflux SMR transporter of EmrE family, glycine cleavage system T protein, MscS small-conductance mechanosensitive channel protein, and two type II restriction-modification systems. Finally, all genome-retrieved plasmids of Psychrobacter spp. were subjected to complex genome- and proteome-based comparative analyses showing that Antarctic replicons are significantly different from plasmids from other locations.

Insights into the physiological and genomic characterization of three bacterial isolates from a highly alkaline, terrestrial serpentinizing system

The terrestrial serpentinite-hosted ecosystem known as "The Cedars" is home to a diverse microbial community persisting under highly alkaline (pH ~ 12) and reducing (Eh < -550 mV) conditions. This extreme environment presents particular difficulties for microbial life, and efforts to isolate microorganisms from The Cedars over the past decade have remained challenging. Herein, we report the initial physiological assessment and/or full genomic characterization of three isolates: Paenibacillus sp. Cedars ('Paeni-Cedars'), Alishewanella sp. BS5-314 ('Ali-BS5-314'), and Anaerobacillus sp. CMMVII ('Anaero-CMMVII'). Paeni-Cedars is a Gram-positive, rod-shaped, mesophilic facultative anaerobe that grows between pH 7-10 (minimum pH tested was 7), temperatures 20-40°C, and 0-3% NaCl concentration. The addition of 10-20 mM CaCl₂ enhanced growth, and iron reduction was observed in the following order, 2-line ferrihydrite > magnetite > serpentinite ~ chromite ~ hematite. Genome analysis identified genes for flavin-mediated iron reduction and synthesis of a bacillibactin-like, catechol-type siderophore. Ali-BS5-314 is a Gram-negative, rod-shaped, mesophilic, facultative anaerobic alkaliphile that grows between pH 10-12 and temperatures 10-40°C, with limited growth observed 1-5% NaCl. Nitrate is used as a terminal electron acceptor under anaerobic conditions, which was corroborated by genome analysis. The Ali-BS5-314 genome also includes genes for benzoate-like compound metabolism. Anaero-CMMVII remained difficult to cultivate for physiological studies; however, growth was observed between pH 9-12, with the addition of 0.01-1% yeast extract. Anaero-CMMVII is a probable oxygen-tolerant anaerobic alkaliphile with hydrogenotrophic respiration coupled with nitrate reduction, as determined by genome analysis. Based on single-copy genes, ANI, AAI and dDDH analyses, Paeni-Cedars and Ali-BS5-314 are related to other species (P. glucanolyticus and A. aestuarii, respectively), and Anaero-CMMVII represents a new species. The characterization of these three isolates demonstrate the range of ecophysiological adaptations and metabolisms present in serpentinite-hosted ecosystems, including mineral reduction, alkaliphily, and siderophore production.

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.

Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches

The Pseudomonas aeruginosa genome can change to adapt to different ecological niches. We compared four genomes from a Mexican hospital and 59 genomes from GenBank from different niches, such as urine, sputum, and environmental. The ST analysis showed that high-risk STs (ST235, ST773, and ST27) were present in the genomes of the three niches from GenBank, and the STs of Mexican genomes (ST167, ST2731, and ST549) differed from the GenBank genomes. Phylogenetic analysis showed that the genomes were clustering according to their ST and not their niche. When analyzing the genomic content, we observed that environmental genomes had genes involved in adapting to the environment not found in the clinics and that their mechanisms of resistance were mutations in antibiotic resistance-related genes. In contrast, clinical genomes from GenBank had resistance genes, in mobile/mobilizable genetic elements in the chromosome, except for the Mexican genomes that carried them mostly in plasmids. This was related to the presence of CRISPR-Cas and anti-CRISPR; however, Mexican strains only had plasmids and CRISPR-Cas. bla_OXA-488 (a variant of bla_OXA50) with higher activity against carbapenems was more prevalent in sputum genomes. The virulome analysis showed that exoS was most prevalent in the genomes of urinary samples and exoU and pldA in sputum samples. This study provides evidence regarding the genetic variability among P. aeruginosa isolated from different niches.

The Plasmidomic Landscape of Clinical Methicillin-Resistant Staphylococcus aureus Isolates from Malaysia

Methicillin-resistant Staphylococcus aureus (MRSA) is a priority nosocomial pathogen with plasmids playing a crucial role in its genetic adaptability, particularly in the acquisition and spread of antimicrobial resistance. In this study, the genome sequences of 79 MSRA clinical isolates from Terengganu, Malaysia, (obtained between 2016 and 2020) along with an additional 15 Malaysian MRSA genomes from GenBank were analyzed for their plasmid content. The majority (90%, 85/94) of the Malaysian MRSA isolates harbored 1-4 plasmids each. In total, 189 plasmid sequences were identified ranging in size from 2.3 kb to ca. 58 kb, spanning all seven distinctive plasmid replication initiator (replicase) types. Resistance genes (either to antimicrobials, heavy metals, and/or biocides) were found in 74% (140/189) of these plasmids. Small plasmids (<5 kb) were predominant (63.5%, 120/189) with a RepL replicase plasmid harboring the ermC gene that confers resistance to macrolides, lincosamides, and streptogramin B (MLS_B) identified in 63 MRSA isolates. A low carriage of conjugative plasmids was observed (n = 2), but the majority (64.5%, 122/189) of the non-conjugative plasmids have mobilizable potential. The results obtained enabled us to gain a rare view of the plasmidomic landscape of Malaysian MRSA isolates and reinforces their importance in the evolution of this pathogen.

Plasmid-Encoded Traits Vary across Environments

Plasmids are key mobile genetic elements in bacterial evolution and ecology as they allow the rapid adaptation of bacteria under selective environmental changes. However, the genetic information associated with plasmids is usually considered separately from information about their environmental origin. To broadly understand what kinds of traits may become mobilized by plasmids in different environments, we analyzed the properties and accessory traits of 9,725 unique plasmid sequences from a publicly available database with known bacterial hosts and isolation sources. Although most plasmid research focuses on resistance traits, such genes made up <1% of the total genetic information carried by plasmids. Similar to traits encoded on the bacterial chromosome, plasmid accessory trait compositions (including general Clusters of Orthologous Genes [COG] functions, resistance genes, and carbon and nitrogen genes) varied across seven broadly defined environment types (human, animal, wastewater, plant, soil, marine, and freshwater). Despite their potential for horizontal gene transfer, plasmid traits strongly varied with their host's taxonomic assignment. However, the trait differences across environments of broad COG categories could not be entirely explained by plasmid host taxonomy, suggesting that environmental selection acts on the plasmid traits themselves. Finally, some plasmid traits and environments (e.g., resistance genes in human-related environments) were more often associated with mobilizable plasmids (those having at least one detected relaxase) than others. Overall, these findings underscore the high level of diversity of traits encoded by plasmids and provide a baseline to investigate the potential of plasmids to serve as reservoirs of adaptive traits for microbial communities. IMPORTANCE Plasmids are well known for their role in the transmission of antibiotic resistance-conferring genes. Beyond human and clinical settings, however, they disseminate many other types of genes, including those that contribute to microbially driven ecosystem processes. In this study, we identified the distribution of traits genetically encoded by plasmids isolated from seven broadly categorized environments. We find that plasmid trait content varied with both bacterial host taxonomy and environment and that, on average, half of the plasmids were potentially mobilizable. As anthropogenic activities impact ecosystems and the climate, investigating and identifying the mechanisms of how microbial communities can adapt will be imperative for predicting the impacts on ecosystem functioning.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Multiple relaxases contribute to the horizontal transfer of the virulence plasmids from the tumorigenic bacterium Pseudomonas syringae pv. savastanoi NCPPB 3335

Pseudomonas syringae pv. savastanoi NCPPB 3335 is the causal agent of olive knot disease and contains three virulence plasmids: pPsv48A (pA), 80 kb; pPsv48B (pB), 45 kb, and pPsv48C (pC), 42 kb. Here we show that pB contains a complete MPF_T (previously type IVA secretion system) and a functional origin of conjugational transfer adjacent to a relaxase of the MOB_P family; pC also contains a functional oriT-MOB_P array, whereas pA contains an incomplete MPF_I (previously type IVB secretion system), but not a recognizable oriT. Plasmid transfer occurred on solid and in liquid media, and on leaf surfaces of a non-host plant (Phaseolus vulgaris) with high (pB) or moderate frequency (pC); pA was transferred only occasionally after cointegration with pB. We found three plasmid-borne and three chromosomal relaxase genes, although the chromosomal relaxases did not contribute to plasmid dissemination. The MOB_P relaxase genes of pB and pC were functionally interchangeable, although with differing efficiencies. We also identified a functional MOB_Q mobilization region in pC, which could only mobilize this plasmid. Plasmid pB could be efficiently transferred to strains of six phylogroups of P. syringae sensu lato, whereas pC could only be mobilized to two strains of phylogroup 3 (genomospecies 2). In two of the recipient strains, pB was stably maintained after 21 subcultures in liquid medium. The carriage of several relaxases by the native plasmids of P. syringae impacts their transfer frequency and, by providing functional diversity and redundancy, adds robustness to the conjugation system.

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.

Epidemiology of Plasmid Lineages Mediating the Spread of Extended-Spectrum Beta-Lactamases among Clinical Escherichia coli

The prevalence of extended-spectrum beta-lactamases (ESBLs) among clinical isolates of Escherichia coli has been increasing, with this spread driven by ESBL-encoding plasmids. However, the epidemiology of ESBL-disseminating plasmids remains understudied, obscuring the roles of individual plasmid lineages in ESBL spread. To address this, we performed an in-depth genomic investigation of 149 clinical ESBL-like E. coli isolates from a tertiary care hospital. We obtained high-quality assemblies for 446 plasmids, revealing an extensive map of plasmid sharing that crosses time, space, and bacterial sequence type boundaries. Through a sequence-based network, we identified specific plasmid lineages that are responsible for the dissemination of major ESBLs. Notably, we demonstrate that IncF plasmids separate into 2 distinct lineages that are enriched for different ESBLs and occupy distinct host ranges. Our work provides a detailed picture of plasmid-mediated spread of ESBLs, demonstrating the extensive sequence diversity within identified lineages, while highlighting the genetic elements that underlie the persistence of these plasmids within the clinical E. coli population. IMPORTANCE The increasing incidence of nosocomial infections with extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli represents a significant threat to public health, given the limited treatment options available for such infections. The rapid ESBL spread is suggested to be driven by localization of the resistance genes on conjugative plasmids. Here, we identify the contributions of different plasmid lineages in the nosocomial spread of ESBLs. We provide further support for plasmid-mediated spread of ESBLs but demonstrate that some ESBL genes rely on dissemination through plasmids more than the others. We identify key plasmid lineages that are enriched in major ESBL genes and highlight the encoded genetic elements that facilitate the transmission and stable maintenance of these plasmid groups within the clinical E. coli population. Overall, our work provides valuable insight into the dissemination of ESBLs through plasmids, furthering our understating of factors underlying the increased prevalence of these genes in nosocomial settings.

ICEscreen: a tool to detect Firmicute ICEs and IMEs, isolated or enclosed in composite structures

Mobile Genetic Elements (MGEs) are integrated in bacterial genomes and key elements that drive prokaryote genome evolution. Among them are Integrative and Conjugative Elements (ICEs) and Integrative Mobilizable Elements (IMEs) which are important for bacterial fitness since they frequently carry genes participating in important bacterial adaptation phenotypes such as antibiotic resistance, virulence or specialized metabolic pathways. Although ICEs and IMEs are widespread, they are as yet almost never annotated in public bacterial genomes. To address the need of dedicated strategies for the annotation of these elements, we developed ICEscreen, a tool that introduces two new features to detect ICEs and IMEs in Firmicute genomes. First, ICEscreen uses an efficient strategy to detect Signature Proteins of ICEs and IMEs based on a database dedicated to Firmicutes and composed of manually curated proteins and Hidden Markov Models (HMM) profiles. Second, ICEscreen includes a new original algorithm that detects composite structures of ICEs and IMEs that are frequent in genomes of Firmicutes but are currently not resolved by any other tool. We benchmarked ICEscreen on experimentally supported elements and on a public dataset of 246 manually annotated elements including the genomes of 40 Firmicutes and demonstrate its efficiency to detect ICEs and IMEs.

Prevalence and Population Diversity of Listeria monocytogenes Isolated from Dairy Cattle Farms in the Cantabria Region of Spain

Simple Summary

The origin and prevalence of Listeria monocytogenes was studied in dairy cattle farms in order to examine its diversity and determine its possible persistence in manure. The utilization of manure for agricultural purposes is common in many countries. While properly treated and managed manure is an effective and safe fertilizer, foodborne illness outbreaks can occur, as many of the most prominent foodborne pathogens are carried by healthy livestock. It is, therefore, necessary to study the origin and persistence of zoonotic agents in general and of L. monocytogenes in particular, in order to avoid recirculation in farms and reduce risk for human populations.

Abstract

Listeria monocytogenes is an opportunistic pathogen that is widely distributed in the environment. Here we show the prevalence and transmission of L. monocytogenes in dairy farms in the Cantabria region, on the northern coast of Spain. A total of 424 samples was collected from 14 dairy farms (5 organic and 9 conventional) and 211 L. monocytogenes isolates were recovered following conventional microbiological methods. There were no statistically significant differences in antimicrobial resistance ratios between organic and conventional farms. A clonal relationship among the isolates was assessed by pulsed field gel electrophoresis (PFGE) analysis and 64 different pulsotypes were obtained. Most isolates (89%, n = 187) were classified as PCR serogroup IVb by using a multiplex PCR assay. In this case, 45 isolates of PCR serogroup IVb were whole genome-sequenced to perform a further analysis at genomic level. In silico MLST analysis showed the presence of 12 sequence types (ST), of which ST1, ST54 and ST666 were the most common. Our data indicate that the environment of cattle farms retains a high incidence of L. monocytogenes, including subtypes involved in human listeriosis reports and outbreaks. This pathogen is shed in the feces and could easily colonize dairy products, as a result of fecal contamination. Effective herd and manure management are needed in order to prevent possible outbreaks.

A Plasmid Carrying blaIMP-56 in Pseudomonas aeruginosa Belonging to a Novel Resistance Plasmid Family

bla_IMP and bla_VIM are the most detected plasmid-encoded carbapenemase genes in Pseudomonas aeruginosa. Previous studies have reported plasmid sequences carrying bla_IMP variants, except bla_IMP-56. In this study, we aimed to characterize a plasmid carrying bla_IMP-56 in a P. aeruginosa strain isolated from a Mexican hospital. The whole genome of P. aeruginosa strain PE52 was sequenced using Illumina Miseq 2 × 150 bp, with 5 million paired-end reads. We characterized a 27 kb plasmid (pPE52IMP) that carried bla_IMP-56. The phylogenetic analysis of RepA in pPE52IMP and 33 P. aeruginosa plasmids carrying resistance genes reported in the GenBank revealed that pPE52IMP and four plasmids (pMATVIM-7, unnamed (FDAARGOS_570), pD5170990, and pMRVIM0713) were in the same clade. These closely related plasmids belonged to the MOB_P11 subfamily and had similar backbones. Another plasmid (p4130-KPC) had a similar backbone to pPE52IMP; however, its RepA was truncated. In these plasmids, the resistance genes bla_KPC-2, bla_VIM variants, aac(6′)-Ib4, bla_OXA variants, and bla_IMP-56 were inserted between phd and resolvase genes. This study describes a new family of plasmids carrying resistance genes, with a similar backbone, the same RepA, and belonging to the MOB_P11 subfamily in P. aeruginosa. In addition, our characterized plasmid harboring bla_IMP-56 (pPE52IMP) belongs to this family.

Assessing the genomic composition, putative ecological relevance and biotechnological potential of plasmids from sponge bacterial symbionts

Plasmid-mediated transfer of genes can have direct consequences in several biological processes within sponge microbial communities. However, very few studies have attempted genomic and functional characterization of plasmids from marine host-associated microbial communities in general and those of sponges in particular. In the present study, we used an endogenous plasmid isolation method to obtain plasmids from bacterial symbionts of the marine sponges Stylissa carteri and Paratetilla sp. and investigated the genomic composition, putative ecological relevance and biotechnological potential of these plasmids. In total, we isolated and characterized three complete plasmids, three plasmid prophages and one incomplete plasmid. Our results highlight the importance of plasmids to transfer relevant genetic traits putatively involved in microbial symbiont adaptation and host-microbe and microbe-microbe interactions. For example, putative genes involved in bacterial response to chemical stress, competition, metabolic versatility and mediation of bacterial colonization and pathogenicity were detected. Genes coding for enzymes and toxins of biotechnological potential were also detected. Most plasmid prophage coding sequences were, however, hypothetical proteins with unknown functions. Overall, this study highlights the ecological relevance of plasmids in the marine sponge microbiome and provides evidence that plasmids of sponge bacterial symbionts may represent an untapped resource of genes of biotechnological interest.

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.

Genomic Analysis of Two MDR Isolates of Salmonella enterica Serovar Infantis from a Spanish Hospital Bearing the blaCTX-M-65 Gene with or without fosA3 in pESI-like Plasmids

Salmonella enterica serovar Infantis (S. Infantis) is a broiler-associated pathogen which ranks in the fourth position as a cause of human salmonellosis in the European Union. Here, we report a comparative genomic analysis of two clinical S. Infantis isolates recovered in Spain from children who just returned from Peru. The isolates were selected on the basis of resistance to cefotaxime, one of the antibiotics of choice for treatment of S. enterica infections. Antimicrobial susceptibility testing demonstrated that they were resistant to eight classes of antimicrobial agents: penicillins, cephalosporins, phenicols, aminoglycosides, tetracyclines, inhibitors of folate synthesis, (fluoro)quinolones and nitrofurans, and one of them was also resistant to fosfomycin. As shown by whole-genome sequence analysis, each isolate carried a pESI-like megaplasmid of ca. 300 kb harboring multiple resistance genes [bla_CTX-M-65, aph(4)-Ia, aac(3)-IVa, aph(3′)-Ia, floR, dfrA14, sul1, tet(A), aadA1 ± fosA3], as well as genes for resistance to heavy metals and disinfectants (mer, ars and qacEΔ1). These genes were distributed in two complex regions, separated by DNA belonging to the plasmid backbone, and associated with a wealth of transposable elements. The two isolates had a D87Y amino acid substitution in the GyrA protein, and truncated variants of the nitroreductase genes nfsA and nsfB, accounting for chromosomally encoded resistances to nalidixic acid and nitrofurantoin, respectively. The two S. Infantis isolates were assigned to sequence type ST32 by in silico multilocus sequence typing (MLST). Phylogenetic analysis revealed that they were closely related, differing only by 12 SNPs, although they were recovered from different children two years apart. They were also genetically similar to bla_CTX-M-65-positive ± fosA3 isolates obtained from humans and along the poultry production chain in the USA, South America, as well as from humans in several European countries, usually associated with a travel history to America. However, this is the first time that the S. Infantis bla_CTX-M-65 ± fosA3 MDR clone has been reported in Spain.

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.

Genomic Surveillance of Vancomycin-Resistant Enterococcus faecium Reveals Spread of a Linear Plasmid Conferring a Nutrient Utilization Advantage

Healthcare-associated outbreaks of vancomycin-resistant Enterococcus faecium (VREfm) are a worldwide problem with increasing prevalence. The genomic plasticity of this hospital-adapted pathogen contributes to its efficient spread despite infection control measures. Here, we aimed to identify the genomic and phenotypic determinants of health care-associated transmission of VREfm. We assessed the VREfm transmission networks at the tertiary-care University Hospital of Zurich (USZ) between October 2014 and February 2018 and investigated microevolutionary dynamics of this pathogen. We performed whole-genome sequencing for the 69 VREfm isolates collected during this time frame and assessed the population structure and variability of the vancomycin resistance transposon. Phylogenomic analysis allowed us to reconstruct transmission networks and to unveil external or wider transmission networks undetectable by routine surveillance. Notably, it unveiled a persistent clone, sampled 31 times over a 29-month period. Exploring the evolutionary dynamics of this clone and characterizing the phenotypic consequences revealed the spread of a variant with decreased daptomycin susceptibility and the acquired ability to utilize N-acetyl-galactosamine (GalNAc), one of the primary constituents of the human gut mucins. This nutrient utilization advantage was conferred by a novel plasmid, termed pELF_USZ, which exhibited a linear topology. This plasmid, which was harbored by two distinct clones, was transferable by conjugation. Overall, this work highlights the potential of combining epidemiological, functional genomic, and evolutionary perspectives to unveil adaptation strategies of VREfm. IMPORTANCE Sequencing microbial pathogens causing outbreaks has become a common practice to characterize transmission networks. In addition to the signal provided by vertical evolution, bacterial genomes harbor mobile genetic elements shared horizontally between clones. While macroevolutionary studies have revealed an important role of plasmids and genes encoding carbohydrate utilization systems in the adaptation of Enterococcus faecium to the hospital environment, mechanisms of dissemination and the specific function of many of these genetic determinants remain to be elucidated. Here, we characterize a plasmid providing a nutrient utilization advantage and show evidence for its clonal and horizontal spread at a local scale. Further studies integrating epidemiological, functional genomics, and evolutionary perspectives will be critical to identify changes shaping the success of this pathogen.

Beyond the ABCs—Discovery of Three New Plasmid Types in Rhodobacterales (RepQ, RepY, RepW)

Copiotrophic marine bacteria of the Roseobacter group (Rhodobacterales, Alphaproteobacteria) are characterized by a multipartite genome organization. We sequenced the genomes of Sulfitobacter indolifex DSM 14862^T and four related plasmid-rich isolates in order to investigate the composition, distribution, and evolution of their extrachromosomal replicons (ECRs). A combination of long-read PacBio and short-read Illumina sequencing was required to establish complete closed genomes that comprised up to twelve ECRs. The ECRs were differentiated in stably evolving chromids and genuine plasmids. Among the chromids, a diagnostic RepABC-8 replicon was detected in four Sulfitobacter species that likely reflects an evolutionary innovation that originated in their common ancestor. Classification of the ECRs showed that the most abundant plasmid system is RepABC, followed by RepA, DnaA-like, and RepB. However, the strains also contained three novel plasmid types that were designated RepQ, RepY, and RepW. We confirmed the functionality of their replicases, investigated the genetic inventory of the mostly cryptic plasmids, and retraced their evolutionary origin. Remarkably, the RepY plasmid of S. pontiacus DSM 110277 is the first high copy-number plasmid discovered in Rhodobacterales.

Mobile Genetic Elements Drive Antimicrobial Resistance Gene Spread in Pasteurellaceae Species

Mobile genetic elements (MGEs) and antimicrobial resistance (AMR) drive important ecological relationships in microbial communities and pathogen-host interaction. In this study, we investigated the resistome-associated mobilome in 345 publicly available Pasteurellaceae genomes, a large family of Gram-negative bacteria including major human and animal pathogens. We generated a comprehensive dataset of the mobilome integrated into genomes, including 10,820 insertion sequences, 2,939 prophages, and 43 integrative and conjugative elements. Also, we assessed plasmid sequences of Pasteurellaceae. Our findings greatly expand the diversity of MGEs for the family, including a description of novel elements. We discovered that MGEs are comparable and dispersed across species and that they also co-occur in genomes, contributing to the family's ecology via gene transfer. In addition, we investigated the impact of these elements in the dissemination and shaping of AMR genes. A total of 55 different AMR genes were mapped to 721 locations in the dataset. MGEs are linked with 77.6% of AMR genes discovered, indicating their important involvement in the acquisition and transmission of such genes. This study provides an uncharted view of the Pasteurellaceae by demonstrating the global distribution of resistance genes linked with MGEs.

Development of a Genomics-Based Approach To Identify Putative Hypervirulent Nontyphoidal Salmonella Isolates: Salmonella enterica Serovar Saintpaul as a Model

While differences in human virulence have been reported across nontyphoidal Salmonella (NTS) serovars and associated subtypes, a rational and scalable approach to identify Salmonella subtypes with differential ability to cause human diseases is not available. Here, we used NTS serovar Saintpaul (S. Saintpaul) as a model to determine if metadata and associated whole-genome sequence (WGS) data in the NCBI Pathogen Detection (PD) database can be used to identify (i) subtypes with differential likelihoods of causing human diseases and (ii) genes and single nucleotide polymorphisms (SNPs) potentially responsible for such differences. S. Saintpaul SNP clusters (n = 211) were assigned different epidemiology types (epi-types) based on statistically significant over- or underrepresentation of human clinical isolates, including human associated (HA; n = 29), non-human associated (NHA; n = 23), and other (n = 159). Comparative genomic analyses identified 384 and 619 genes overrepresented among isolates in 5 HA and 4 NHA SNP clusters most significantly associated with the respective isolation source. These genes included 5 HA-associated virulence genes previously reported to be present on Gifsy-1/Gifsy-2 prophages. Additionally, premature stop codons in 3 and 7 genes were overrepresented among the selected HA and NHA SNP clusters, respectively. Tissue culture experiments with strains representing 4 HA and 3 NHA SNP clusters did not reveal evidence for enhanced invasion or intracellular survival for HA strains. However, the presence of sodCI (encoding a superoxide dismutase), found in 4 HA and 1 NHA SNP clusters, was positively correlated with intracellular survival in macrophage-like cells. Post hoc analyses also suggested a possible difference in intracellular survival among S. Saintpaul lineages. IMPORTANCE Not all Salmonella isolates are equally likely to cause human disease, and Salmonella control strategies may unintentionally focus on serovars and subtypes with high prevalence in source populations but are rarely associated with human clinical illness. We describe a framework leveraging WGS data in the NCBI PD database to identify Salmonella subtypes over- and underrepresented among human clinical cases. While we identified genomic signatures associated with HA/NHA SNP clusters, tissue culture experiments failed to identify consistent phenotypic characteristics indicative of enhanced human virulence of HA strains. Our findings illustrate the challenges of defining hypo- and hypervirulent S. Saintpaul and potential limitations of phenotypic assays when evaluating human virulence, for which in vivo experiments are essential. Identification of sodCI, an HA-associated virulence gene associated with enhanced intracellular survival, however, illustrates the potential of the framework and is consistent with prior work identifying specific genomic features responsible for enhanced or reduced virulence of nontyphoidal Salmonella.

Structural and biochemical characterization of the relaxosome auxiliary proteins encoded on the Bacillus subtilis plasmid pLS20

Bacterial conjugation is an important route for horizontal gene transfer. The initial step in this process involves a macromolecular protein-DNA complex called the relaxosome, which in plasmids consists of the origin of transfer (oriT) and several proteins that prepare the transfer. The relaxosome protein named relaxase introduces a nick in one of the strands of the oriT to initiate the process. Additional relaxosome proteins can exist. Recently, several relaxosome proteins encoded on the Bacillus subtilis plasmid pLS20 were identified, including the relaxase, named Rel_pLS20, and two auxiliary DNA-binding factors, named Aux1_pLS20 and Aux2_pLS20. Here, we extend this characterization in order to define their function. We present the low-resolution SAXS envelope of the Aux1_pLS20 and the atomic X-ray structure of the C-terminal domain of Aux2_pLS20. We also study the interactions between the auxiliary proteins and the full-length Rel_pLS20, as well as its separate domains. The results show that the quaternary structure of the auxiliary protein Aux1_pLS20 involves a tetramer, as previously determined. The crystal structure of the C-terminal domain of Aux2_pLS20 shows that it forms a tetramer and suggests that it is an analog of TraM_pF of plasmid F. This is the first evidence of the existence of a TraM_pF analog in gram positive conjugative systems, although, unlike other TraM_pF analogs, Aux2_pLS20 does not interact with the relaxase. Aux1_pLS20 interacts with the C-terminal domain, but not the N-terminal domain, of the relaxase Rel_pLS20. Thus, the pLS20 relaxosome exhibits some unique features despite the apparent similarity to some well-studied G- conjugation systems.

The Facts and Family Secrets of Plasmids That Replicate via the Rolling-Circle Mechanism

Plasmids are self-replicative DNA elements that are transferred between bacteria. Plasmids encode not only antibiotic resistance genes but also adaptive genes that allow their hosts to colonize new niches. Plasmid transfer is achieved by conjugation (or mobilization), phage-mediated transduction, and natural transformation. Thousands of plasmids use the rolling-circle mechanism for their propagation (RCR plasmids). They are ubiquitous, have a high copy number, exhibit a broad host range, and often can be mobilized among bacterial species. Based upon the replicon, RCR plasmids have been grouped into several families, the best known of them being pC194 and pUB110 (Rep_1 family), pMV158 and pE194 (Rep_2 family), and pT181 and pC221 (Rep_trans family). Genetic traits of RCR plasmids are analyzed concerning (i) replication mediated by a DNA-relaxing initiator protein and its interactions with the cognate DNA origin, (ii) lagging-strand origins of replication, (iii) antibiotic resistance genes, (iv) mobilization functions, (v) replication control, performed by proteins and/or antisense RNAs, and (vi) the participating host-encoded functions. The mobilization functions include a relaxase initiator of transfer (Mob), an origin of transfer, and one or two small auxiliary proteins. There is a family of relaxases, the MOB_V family represented by plasmid pMV158, which has been revisited and updated. Family secrets, like a putative open reading frame of unknown function, are reported. We conclude that basic research on RCR plasmids is of importance, and our perspectives contemplate the concept of One Earth because we should incorporate bacteria into our daily life by diminishing their virulence and, at the same time, respecting their genetic diversity.

Evolution of ColE1-like plasmids across γ-Proteobacteria: From bacteriocin production to antimicrobial resistance

Antimicrobial resistance is one of the major threats to Public Health worldwide. Understanding the transfer and maintenance of antimicrobial resistance genes mediated by mobile genetic elements is thus urgent. In this work, we focus on the ColE1-like plasmid family, whose distinctive replication and multicopy nature has given rise to key discoveries and tools in molecular biology. Despite being massively used, the hosts, functions, and evolutionary history of these plasmids remain poorly known. Here, we built specific Hidden Markov Model (HMM) profiles to search ColE1 replicons within genomes. We identified 1,035 ColE1 plasmids in five Orders of γ-Proteobacteria, several of which are described here for the first time. The phylogenetic analysis of these replicons and their characteristic MOB_P5/HEN relaxases suggest that ColE1 plasmids have diverged apart, with little transfer across orders, but frequent transfer across families. Additionally, ColE1 plasmids show a functional shift over the last decades, losing their characteristic bacteriocin production while gaining several antimicrobial resistance genes, mainly enzymatic determinants and including several extended-spectrum betalactamases and carbapenemases. Furthermore, ColE1 plasmids facilitate the intragenomic mobilization of these determinants, as various replicons were identified co-integrated with large non-ColE1 plasmids, mostly via transposases. These results illustrate how families of plasmids evolve and adapt their gene repertoires to bacterial adaptive requirements.

The extraordinary adaptability of bacteria and the massive prevalence of mobile genetic elements within populations has turned antimicrobial resistance into a growing threat to Public Health. Among all the mobile genetic elements, plasmids have been the focus of attention as these extrachromosomal molecules of DNA are able to mobilize several antimicrobial resistance genes at once through conjugation. However, although small mobilizable and non-conjugative replicons have been traditionally overlooked when analyzing plasmid-mediated antimicrobial resistance, they have recently been described as important carriers of AMR genes. In this work, we have analyzed the ColE1-like plasmid family, whose study has been neglected even if they are one of the main groups of small plasmids in natural populations of Proteobacteria. We observed that these plasmids have evolved for a long time within γ-Proteobacteria acquiring different genetic features in specific hosts, being major players in the spread of antimicrobial resistance determinants.

PLASmid TAXonomic PCR (PlasTax-PCR), a Multiplex Relaxase MOB Typing to Assort Plasmids into Taxonomic Units

Plasmids transmissible by conjugation are responsible for disseminating antibiotic-resistance genes, making plasmid detection relevant for pathogen tracking. We describe the use of a multiplex PCR method for the experimental identification of specific plasmid taxonomic units (PTUs) of transmissible plasmids. The PCR primers were designed to target conserved segments of the relaxase MOB gene of PTUs encoding adaptive traits for enterobacteria (antimicrobial resistance, virulence, and metabolism). In this way, PlasTax-PCR detects the presence of these plasmids and allows their direct assignation to a PTU.

Assessment of the plasmidome of an extremophilic microbial community from the Diamante Lake, Argentina

Diamante Lake located at 4589 m.a.s.l. in the Andean Puna constitutes an extreme environment. It is exposed to multiple extreme conditions such as an unusually high concentration of arsenic (over 300 mg L^-1) and low oxygen pressure. Microorganisms thriving in the lake display specific genotypes that facilitate survival, which include at least a multitude of plasmid-encoded resistance traits. Hence, the genetic information provided by the plasmids essentially contributes to understand adaptation to different stressors. Though plasmids from cultivable organisms have already been analyzed to the sequence level, the impact of the entire plasmid-borne genetic information on such microbial ecosystem is not known. This study aims at assessing the plasmidome from Diamante Lake, which facilitates the identification of potential hosts and prediction of gene functions as well as the ecological impact of mobile genetic elements. The deep-sequencing analysis revealed a large fraction of previously unknown DNA sequences of which the majority encoded putative proteins of unknown function. Remarkably, functions related to the oxidative stress response, DNA repair, as well as arsenic- and antibiotic resistances were annotated. Additionally, all necessary capacities related to plasmid replication, mobilization and maintenance were detected. Sequences characteristic for megaplasmids and other already known plasmid-associated genes were identified as well. The study highlights the potential of the deep-sequencing approach specifically targeting plasmid populations as it allows to evaluate the ecological impact of plasmids from (cultivable and non-cultivable) microorganisms, thereby contributing to the understanding of the distribution of resistance factors within an extremophilic microbial community.

BioF is a novel B2 metallo-β-lactamase from Pseudomonas sp. isolated from an on-farm biopurification system

Antimicrobial resistance represents a major global health concern and environmental bacteria are considered a source of resistance genes. Carbapenems are often used as the last antibiotic option to treat multidrug-resistant bacteria. Metallo-β-lactamases (MBLs) are able to render resistance to almost all β-lactam antibiotics, including carbapenems. Unfortunately, there are no inhibitors against MBLs for clinical use. Subclass B2 MBLs are the only enzymes working as strict carbapenemases, under-represented, encoded in chromosome genes and only functional as mono-zinc enzymes. Despite current efforts in MBLs inhibitor development, B2 carbapenemase activity is especially difficult to suppress, even in vitro. In this study we characterized BioF, a novel subclass B2 MBL identified in a new environmental Pseudomonas sp. strain isolated from an on-farm biopurification system (BPS). Although bla_BioF is most likely a chromosomal gene, it is found in a genomic island and may represent a step previous to the horizontal transmission of B2 genes. The new B2 MBL is active as a mono-zinc enzyme and is a potent carbapenemase with incipient activity against some cephalosporins. BioF activity is not affected by excess zinc and is only inhibited at high metal chelator concentrations. The discovery and characterization of B2 MBL BioF as a potent carbapenemase in a BPS bacterial isolate emphasizes the importance of exploring antibiotic resistances existing in the environmental microbiota under the influence of human activities before they could emerge clinically.

Genetic Analysis of mcr-1-Carrying Plasmids From Gram-Negative Bacteria in a Dutch Tertiary Care Hospital: Evidence for Intrapatient and Interspecies Transmission Events

The role of plasmids in the complex pandemic of antimicrobial resistance is increasingly being recognized. In this respect, multiple mobile colistin resistance (mcr) gene-carrying plasmids have been described. However, the characteristics and epidemiology of these plasmids within local healthcare settings are largely unknown. We retrospectively characterized the genetic composition and epidemiology of plasmids from mcr-1-positive bacterial isolates identified from patients from a large academic hospital in the Netherlands. Clinical Gram-negative bacteria with an MIC > 2 μg/mL for colistin, obtained from patients hospitalized at the Erasmus MC University Medical Center Rotterdam during the years 2010-2018, were screened for presence of the mcr-1 gene. Extracted plasmids from mcr-1-positive isolates were sequenced using a combination of short- and long-read sequencing platforms, characterized by incompatibility type and genetic composition and compared to publicly available mcr-1-carrying plasmid sequences. In 21 isolates from 14 patients, mcr-1 was located on a plasmid. These plasmids were of diverse genetic background involving Inc types IncX4, IncI2(delta), IncHI2, as well as double Inc types IncHI2/IncN and IncHI2/IncQ. mcr-1-carrying plasmids were found in Escherichia coli, Klebsiella pneumoniae, and Kluyvera georgiana, and within the chromosome of an ST147 K. pneumoniae isolate. In depth analysis indicated intrapatient, interpatient, and interspecies transmission events of mcr-1-carrying plasmids. In addition, our results show that the mcr-1 gene resides in a rich environment full of other (mcr-1 negative) plasmids and of many different Inc types, enabling interplasmidal transfer events and facilitating widespread dissemination of the mcr-1 gene. Multiple mcr-1-carrying plasmid transmission events had likely occurred among isolates from hospitalized patients. Recognition and identification of plasmid transmission events within hospitals is necessary in order to design and implement effective infection control measures.

COPLA, a taxonomic classifier of plasmids

Background

Plasmids are mobile genetic elements, key in the dissemination of antibiotic resistance, virulence determinants and other adaptive traits in bacteria. Obtaining a robust method for plasmid classification is necessary to better understand the genetics and epidemiology of many pathogens. Until now, plasmid classification systems focused on specific traits, which limited their precision and universality. The definition of plasmid taxonomic units (PTUs), based on average nucleotide identity metrics, allows the generation of a universal plasmid classification scheme, applicable to all bacterial taxa. Here we present COPLA, a software able to assign plasmids to known and novel PTUs, based on their genomic sequence.

Results

We implemented an automated pipeline able to assign a given plasmid DNA sequence to its cognate PTU, and assessed its performance using a sample of 1000 unclassified plasmids. Overall, 41% of the samples could be assigned to a previously defined PTU, a number that reached 63% in well-known taxa such as the Enterobacterales order. The remaining plasmids represent novel PTUs, indicating that a large fraction of plasmid backbones is still uncharacterized.

Conclusions

COPLA is a bioinformatic tool for universal, species-independent, plasmid classification. Offered both as an automatable pipeline and an open web service, COPLA will help bacterial geneticists and clinical microbiologists to quickly classify plasmids.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04299-x.

Comparative Genomics Suggests Mechanisms of Genetic Adaptation toward the Catabolism of the Phenylurea Herbicide Linuron in Variovorax

Biodegradation of the phenylurea herbicide linuron appears a specialization within a specific clade of the Variovorax genus. The linuron catabolic ability is likely acquired by horizontal gene transfer but the mechanisms involved are not known. The full-genome sequences of six linuron-degrading Variovorax strains isolated from geographically distant locations were analyzed to acquire insight into the mechanisms of genetic adaptation toward linuron metabolism. Whole-genome sequence analysis confirmed the phylogenetic position of the linuron degraders in a separate clade within Variovorax and indicated that they unlikely originate from a common ancestral linuron degrader. The linuron degraders differentiated from Variovorax strains that do not degrade linuron by the presence of multiple plasmids of 20–839 kb, including plasmids of unknown plasmid groups. The linuron catabolic gene clusters showed 1) high conservation and synteny and 2) strain-dependent distribution among the different plasmids. Most of them were bordered by IS1071 elements forming composite transposon structures, often in a multimeric array configuration, appointing IS1071 as a key element in the recruitment of linuron catabolic genes in Variovorax. Most of the strains carried at least one (catabolic) broad host range plasmid that might have been a second instrument for catabolic gene acquisition. We conclude that clade 1 Variovorax strains, despite their different geographical origin, made use of a limited genetic repertoire regarding both catabolic functions and vehicles to acquire linuron biodegradation.

Co-harboring of Novel bla KPC-2 Plasmid and Integrative and Conjugative Element Carrying Tn6203 in Multidrug-Resistant Pseudomonas aeruginosa

Many strains of the opportunistic pathogen Pseudomonas aeruginosa have acquired resistance to multiple antibiotics. Carbapenem-resistant P. aeruginosa poses a global healthcare problem due to limited therapeutic options for the treatment of infections. Plasmids and integrative and conjugative elements (ICEs) are the major vectors of antibiotic-resistance gene transfer. In our study, four carbapenem-resistant strains of P. aeruginosa were isolated from the same patient in a tertiary referral hospital in China, one of these was resistant to gentamicin and tobramycin. In this strain P33, we observed a non-transferable plasmid, pP33-2, carrying a novel bla_KPC−2 gene segment (ISKpn27-bla_KPC−2-ISKpn6-korC-ORF-klcA-IS26), which we concluded to have been formed by IS26-mediated gene cluster translocation. In addition, by comparing the chromosomes of the P. aeruginosa strains that belong to the same sequence type, we identified an ICE, ICEP33, adjacent to a prophage. The attL site of ICEP33 is identical to the terminal part of the attR site of the prophage. The ICEP33 element contains the transposon Tn6203, which encodes antibiotic and metal resistance genes. The insertion of ICEP33 in the chromosome mediates resistance to multiple antibiotics. Our study contributes to the understanding of the acquisition of antibiotic resistance in P. aeruginosa facilitated by mobile genetic elements.