A curated dataset of complete Enterobacteriaceae plasmids compiled from the NCBI nucleotide database

Genome-wide identification and expression profiling of two-component system (TCS) genes in Brassica oleracea in response to shade stress

The Two-component system (TCS) consists of Histidine kinases (HKs), Phosphotransfers (HPs), and response regulator (RR) proteins. It has an important role in signal transduction to respond to a wide variety of abiotic stresses and hence in plant development. Brassica oleracea (cabbage) is a leafy vegetable, which is used for food and medicinal purposes. Although this system was identified in several plants, it had not been identified in Brassica oleracea yet. This genome-wide study identified 80 BoTCS genes consisting of 21 HKs, 8 HPs, 39 RRs, and 12 PRRs. This classification was done based on conserved domains and motif structure. Phylogenetic relationships of BoTCS genes with Arabidopsis thaliana, Oryza sativa, Glycine max, and Cicer arietinum showed conservation in TCS genes. Gene structure analysis revealed that each subfamily had conserved introns and exons. Both tandem and segmental duplication led to the expansion of this gene family. Almost all of the HPs and RRs were expanded through segmental duplication. Chromosomal analysis showed that BoTCS genes were dispersed across all nine chromosomes. The promoter regions of these genes were found to contain a variety of cis-regulatory elements. The 3D structure prediction of proteins also confirmed the conservation of structure within subfamilies. MicroRNAs (miRNAs) involved in the regulation of BoTCSs were also predicted and their regulatory roles were also evaluated. Moreover, BoTCSs were docked with abscisic acid to evaluate their binding. RNA-seq-based expression analysis and validation by qRT-PCR showed significant variation of expression for BoPHYs, BoERS1.1, BoERS2.1, BoERS2.2, BoRR10.2, and BoRR7.1 suggesting their importance in stress response. These genes showing unique expression can be further used in manipulating the plant's genome to make the plant more resistant the environmental stresses which will ultimately help in the increase of plant's yield. More specifically, these genes have altered expression in shade stress which clearly indicates their importance in biological functions. These findings are important for future functional characterization of TCS genes in generating stress-responsive cultivars.

Horizontal gene transfer enables programmable gene stability in synthetic microbiota

The functions of many microbial communities exhibit remarkable stability despite fluctuations in the compositions of these communities. To date, a mechanistic understanding of this function-composition decoupling is lacking. Statistical mechanisms have been commonly hypothesized to explain such decoupling. Here, we proposed that dynamic mechanisms, mediated by horizontal gene transfer (HGT), also enable the independence of functions from the compositions of microbial communities. We combined theoretical analysis with numerical simulations to illustrate that HGT rates can determine the stability of gene abundance in microbial communities. We further validated these predictions using engineered microbial consortia of different complexities transferring one or more than a dozen clinically isolated plasmids, as well as through the reanalysis of data from the literature. Our results demonstrate a generalizable strategy to program the gene stability of microbial communities.

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

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

Genome-Wide Characterization of Superoxide Dismutase (SOD) Genes in Daucus carota: Novel Insights Into Structure, Expression, and Binding Interaction With Hydrogen Peroxide (H2O2) Under Abiotic Stress Condition

Superoxide dismutase (SOD) proteins are important antioxidant enzymes that help plants to grow, develop, and respond to a variety of abiotic stressors. SOD gene family has been identified in a number of plant species but not yet in Daucus carota. A total of 9 DcSOD genes, comprising 2 FeSODs, 2 MnSODs, and 5 Cu/ZnSODs, are identified in the complete genome of D. carota, which are dispersed in five out of nine chromosomes. Based on phylogenetic analysis, SOD proteins from D. carota were categorized into two main classes (Cu/ZnSODs and MnFeSODs). It was predicted that members of the same subgroups have the same subcellular location. The phylogenetic analysis was further validated by sequence motifs, exon-intron structure, and 3D protein structures, with each subgroup having a similar gene and protein structure. Cis-regulatory elements responsive to abiotic stresses were identified in the promoter region, which may contribute to their differential expression. Based on RNA-seq data, tissue-specific expression revealed that DcCSD2 had higher expression in both xylem and phloem. Moreover, DcCSD2 was differentially expressed in dark stress. All SOD genes were subjected to qPCR analysis after cold, heat, salt, or drought stress imposition. SODs are antioxidants and play a critical role in removing reactive oxygen species (ROS), including hydrogen peroxide (H2O2). DcSODs were docked with H2O2 to evaluate their binding. The findings of this study will serve as a basis for further functional insights into the DcSOD gene family.

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

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

Characterization of the genetic structure of mcr-1 gene among Escherichia coli isolates recovered from surface waters and sediments from Ecuador

Although anthropogenic activities contribute to the selection and spread of antibiotic resistance in aquatic environments, limited information is available from countries with absent or incomplete sewage treatment systems and the impact of their discharges onto water bodies. This study therefore aimed to characterize the genetic structure of colistin resistance (mcr) genes among Escherichia coli isolates recovered from surface waters and sediments in Ecuador. Out of 459 isolates, four Escherichia coli showed multidrug-resistant phenotypes, which harbored the mcr-1 gene and β-lactamases, such as bla_TEM, bla_CTX-M-15, bla_CTX-M-55, or bla_CTX-M-65 genes. Three E. coli isolates (U20, U30 and U144) shared a similar genetic environment surrounding the mcr-1 gene, which was located on plasmids. Only one E. coli isolate (U175) showed that the mcr-1 gene was chromosomally located. Moreover, the core genome multilocus sequence typing (cgMLST) analysis revealed that these isolates belong to different lineages. This study represents the first detection of the mcr-1 gene in multidrug-resistant E. coli isolates from environmental samples in Ecuador.

Clonal Spread and Intra- and Inter-Species Plasmid Dissemination Associated With Klebsiella pneumoniae Carbapenemase-Producing Enterobacterales During a Hospital Outbreak in Barcelona, Spain

Objectives: The study aimed to characterize the clonal spread of resistant bacteria and dissemination of resistance plasmids among carbapenem-resistant Enterobacterales at a tertiary hospital in Catalonia, Spain. Methods: Isolates were recovered from surveillance rectal swabs and diagnostic samples. Species identification was by matrix-assisted laser desorption ionization-time time of flight mass spectrometry (MALDI-TOF MS). Molecular typing was performed by pulsed-field gel electrophoresis (PFGE) and multi-locus sequence typing (MLST). Antimicrobial susceptibility was assessed by gradient-diffusion and carriage of bla genes was detected by PCR. Plasmid typing, conjugation assays, S1-PFGE studies and long-read sequencing were used to characterize resistance plasmids. Results: From July 2018 to February 2019, 125 Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacterales were recovered from 101 inpatients from surveillance (74.4%) or clinical samples (25.6%), in a tertiary hospital in Barcelona. Clonality studies identified a major clone of Klebsiella pneumoniae belonging to sequence type ST15 and additional isolates of K. pneumoniae, Escherichia coli and Enterobacter sp. from different STs. All isolates but one carried the bla KPC-2 allelic variant. The bla KPC-2 gene was located in an IncFIIk plasmid of circa 106 Kb in a non-classical Tn4401 element designated NTEKPC-pMC-2-1. Whole-genome sequencing revealed different rearrangements of the 106 Kb plasmid while the NTEKPC-pMC-2-1 module was highly conserved. Conclusion: We report a hospital outbreak caused by the clonal dissemination of KPC-producing ST15 K. pneumoniae but also the intra- and inter-species transmission of the bla KPC-2 gene associated with plasmid conjugation and/or transposon dissemination. To our knowledge, this is the first report of an outbreak caused by KPC-producing Enterobacterales isolated from human patients in Catalonia and highlights the relevance of surveillance studies in the early detection and control of antibiotic resistant high-risk clones.

Dissemination of NDM-producing Klebsiella pneumoniae and Escherichia coli high-risk clones in Catalan healthcare institutions

OBJECTIVES

To characterize the clonal spread of carbapenem-resistant Klebsiella pneumoniae and Escherichia coli isolates between different healthcare institutions in Catalonia, Spain.

METHODS

Antimicrobial susceptibility was tested by disc diffusion. MICs were determined by gradient diffusion or broth microdilution. Carbapenemase production was confirmed by lateral flow. PCR and Sanger sequencing were used to identify the allelic variants of resistance genes. Clonality studies were performed by PFGE and MLST. Plasmid typing, conjugation assays, S1-PFGE plus Southern blotting and MinION Oxford Nanopore sequencing were used to characterize resistance plasmids.

RESULTS

Twenty-nine carbapenem-resistant isolates recovered from three healthcare institutions between January and November 2016 were included: 14 K. pneumoniae isolates from a tertiary hospital in the south of Catalonia (hospital A); 2 K. pneumoniae isolates from a nearby healthcare centre; and 12 K. pneumoniae isolates and 1 E. coli isolate from a tertiary hospital in Barcelona (hospital B). The majority of isolates were resistant to all antimicrobial agents, except colistin, and all were NDM producers. PFGE identified a major K. pneumoniae clone (n = 27) belonging to ST147 and co-producing NDM-1 and CTX-M-15, with a few isolates also harbouring blaOXA-48. Two sporadic isolates of K. pneumoniae ST307 and E. coli ST167 producing NDM-7 were also identified. blaNDM-1 was carried in two related IncR plasmid populations and blaNDM-7 in a conjugative 50 kb IncX3 plasmid.

CONCLUSIONS

We report the inter-hospital dissemination of XDR high-risk clones of K. pneumoniae and E. coli associated with the carriage of small, transferable plasmids harbouring blaNDM genes.

Detection of plasmid contigs in draft genome assemblies using customized Kraken databases

Plasmids play an important role in bacterial evolution and mediate horizontal transfer of genes including virulence and antimicrobial resistance genes. Although short-read sequencing technologies have enabled large-scale bacterial genomics, the resulting draft genome assemblies are often fragmented into hundreds of discrete contigs. Several tools and approaches have been developed to identify plasmid sequences in such assemblies, but require trade-off between sensitivity and specificity. Here we propose using the Kraken classifier, together with a custom Kraken database comprising known chromosomal and plasmid sequences of Klebsiella pneumoniae species complex (KpSC), to identify plasmid-derived contigs in draft assemblies. We assessed performance using Illumina-based draft genome assemblies for 82 KpSC isolates, for which complete genomes were available to supply ground truth. When benchmarked against five other classifiers (Centrifuge, RFPlasmid, mlplasmids, PlaScope and Platon), Kraken showed balanced performance in terms of overall sensitivity and specificity (90.8 and 99.4 %, respectively, for contig count; 96.5 and >99.9 %, respectively, for cumulative contig length), and the highest accuracy (96.8% vs 91.8-96.6% for contig count; 99.8% vs 99.0-99.7 % for cumulative contig length), and F1-score (94.5 % vs 84.5-94.1 %, for contig count; 98.0 % vs 88.9-96.7 % for cumulative contig length). Kraken also achieved consistent performance across our genome collection. Furthermore, we demonstrate that expanding the Kraken database with additional known chromosomal and plasmid sequences can further improve classification performance. Although we have focused here on the KpSC, this methodology could easily be applied to other species with a sufficient number of completed genomes.

Vertical and Horizontal Transmission of ESBL Plasmid from Escherichia coli O104:H4

Multidrug resistance (MDR) often results from the acquisition of mobile genetic elements (MGEs) that encode MDR gene(s), such as conjugative plasmids. The spread of MDR plasmids is founded on their ability of horizontal transference, as well as their faithful inheritance in progeny cells. Here, we investigated the genetic factors involved in the prevalence of the IncI conjugative plasmid pESBL, which was isolated from the Escherichia coli O104:H4 outbreak strain in Germany in 2011. Using transposon-insertion sequencing, we identified the pESBL partitioning locus (par). Genetic, biochemical and microscopic approaches allowed pESBL to be characterized as a new member of the Type Ib partitioning system. Inactivation of par caused mis-segregation of pESBL followed by post-segregational killing (PSK), resulting in a great fitness disadvantage but apparent plasmid stability in the population of viable cells. We constructed a variety of pESBL derivatives with different combinations of mutations in par, conjugational transfer (oriT) and pnd toxin-antitoxin (TA) genes. Only the triple mutant exhibited plasmid-free cells in viable cell populations. Time-lapse tracking of plasmid dynamics in microfluidics indicated that inactivation of pnd improved the survival of plasmid-free cells and allowed oriT-dependent re-acquisition of the plasmid. Altogether, the three factors—active partitioning, toxin-antitoxin and conjugational transfer—are all involved in the prevalence of pESBL in the E. coli population.

Sequencing of five poultry strains elucidates phylogenetic relationships and divergence in virulence genes in Morganella morganii

Background

M. morganii is a bacterium frequently associated with urinary infections in humans. While many human strains are sequenced, only the genomes of few poultry strains are available. Here, we performed a detailed characterization of five highly resistant Morganella morganii strains isolated in association with Escherichia coli from diseased domestic Austrian poultry flocks, namely geese, turkeys and chicken layers. Additionally, we sequenced the genomes of these strains by NGS and analyzed phylogenetic clustering, resistance and virulence genes in the context of host-specificity.

Results

Two strains were identified to be Extended Spectrum Beta Lactamase (ESBL) and one as AmpC beta-lactamases (AMP-C) phenotype, while two were ESBL negative. By integrating the genome sequences of these five poultry strains with all the available M. morganii genomes, we constructed a phylogenetic tree that clearly separates the Morganella genus into two clusters (M1 and M2), which approximately reflect the proposed subspecies classification (morganii and sibonii). Additionally, we found no association between phylogenetic structure and host, suggesting interspecies transmission. All five poultry strains contained genes for resistance to aminocoumarins, beta-lactams, colistin, elfamycins, fluoroquinolones, phenicol, rifampin and tetracycline. A comparative genomics analysis of virulence genes showed acquisition of novel virulence genes involved in secretion system and adherence in cluster M2. We showed that some of these genes were acquired by horizontal gene transfer from closely related Morganellaceae species and propose that novel virulence genes could be responsible for expansion of tissue tropism in M. morganii. Finally, we detected variability in copy number and high sequence divergence in toxin genes and provided evidence for positive selection in insecticidal toxins genes, likely reflecting host-related adaptations.

Conclusions

In summary, this study describes i) the first isolation and characterization of M. morganii from goose and turkey, ii) a large-scale genetic analysis of M. morganii and an attempt to generate a global picture of the M. morganii intraspecific phylogenetic structure.

Microbiome Diagnostics

BACKGROUND

During the past decade, breakthroughs in sequencing technology and computational biology have provided the basis for studies of the myriad ways in which microbial communities ("microbiota") in and on the human body influence human health and disease. In almost every medical specialty, there is now a growing interest in accurate and replicable profiling of the microbiota for use in diagnostic and therapeutic application.

CONTENT

This review provides an overview of approaches, challenges, and considerations for diagnostic applications borrowing from other areas of molecular diagnostics, including clinical metagenomics. Methodological considerations and evolving approaches for microbiota profiling from mitochondrially encoded 16S rRNA-based amplicon sequencing to metagenomics and metatranscriptomics are discussed. To improve replicability, at least the most vulnerable steps in testing workflows will need to be standardized and continuous efforts needed to define QC standards. Challenges such as purity of reagents and consumables, improvement of reference databases, and availability of diagnostic-grade data analysis solutions will require joint efforts across disciplines and with manufacturers.

SUMMARY

The body of literature supporting important links between the microbiota at different anatomic sites with human health and disease is expanding rapidly and therapeutic manipulation of the intestinal microbiota is becoming routine. The next decade will likely see implementation of microbiome diagnostics in diagnostic laboratories to fully capitalize on technological and scientific advances and apply them in routine medical practice.

Large-scale network analysis captures biological features of bacterial plasmids

Many bacteria can exchange genetic material through horizontal gene transfer (HGT) mediated by plasmids and plasmid-borne transposable elements. Here, we study the population structure and dynamics of over 10,000 bacterial plasmids, by quantifying their genetic similarities and reconstructing a network based on their shared k-mer content. We use a community detection algorithm to assign plasmids into cliques, which correlate with plasmid gene content, bacterial host range, GC content, and existing classifications based on replicon and mobility (MOB) types. Further analysis of plasmid population structure allows us to uncover candidates for yet undescribed replicon genes, and to identify transposable elements as the main drivers of HGT at broad phylogenetic scales. Our work illustrates the potential of network-based analyses of the bacterial 'mobilome' and opens up the prospect of a natural, exhaustive classification framework for bacterial plasmids.

PLSDB: a resource of complete bacterial plasmids

The study of bacterial isolates or communities requires the analysis of the therein included plasmids in order to provide an extensive characterization of the organisms. Plasmids harboring resistance and virulence factors are of especial interest as they contribute to the dissemination of antibiotic resistance. As the number of newly sequenced bacterial genomes is growing a comprehensive resource is required which will allow to browse and filter the available plasmids, and to perform sequence analyses. Here, we present PLSDB, a resource containing 13 789 plasmid records collected from the NCBI nucleotide database. The web server provides an interactive view of all obtained plasmids with additional meta information such as sequence characteristics, sample-related information and taxonomy. Moreover, nucleotide sequence data can be uploaded to search for short nucleotide sequences (e.g. specific genes) in the plasmids, to compare a given plasmid to the records in the collection or to determine whether a sample contains one or multiple of the known plasmids (containment analysis). The resource is freely accessible under https://ccb-microbe.cs.uni-saarland.de/plsdb/.

Comparative Whole-Genome Phylogeny of Animal, Environmental, and Human Strains Confirms the Genogroup Organization and Diversity of the Stenotrophomonas maltophilia Complex

The Stenotrophomonas maltophilia complex (Smc) comprises opportunistic environmental Gram-negative bacilli responsible for a variety of infections in both humans and animals. Beyond its large genetic diversity, its genetic organization in genogroups was recently confirmed through the whole-genome sequencing of human and environmental strains. As they are poorly represented in these analyses, we sequenced the whole genomes of 93 animal strains to determine their genetic background and characteristics. Combining these data with 81 newly sequenced human strains and the genomes available from RefSeq, we performed a genomic analysis that included 375 nonduplicated genomes with various origins (animal, 104; human, 226; environment, 30; unknown, 15). Phylogenetic analysis and clustering based on genome-wide average nucleotide identity confirmed and specified the genetic organization of Smc in at least 20 genogroups. Two new genogroups were identified, and two previously described groups were further divided into two subgroups each. Comparing the strains isolated from different host types and their genogroup affiliation, we observed a clear disequilibrium in certain groups. Surprisingly, some antimicrobial resistance genes, integrons, and/or clusters of attC sites lacking integron-integrase (CALIN) sequences targeting antimicrobial compounds extensively used in animals were mainly identified in animal strains. We also identified genes commonly found in animal strains coding for efflux systems. The result of a large whole-genome analysis performed by us supports the hypothesis of the putative contribution of animals as a reservoir of Stenotrophomonas maltophilia complex strains and/or resistance genes for strains in humans.IMPORTANCE Given its naturally large antimicrobial resistance profile, the Stenotrophomonas maltophilia complex (Smc) is a set of emerging pathogens of immunosuppressed and cystic fibrosis patients. As it is group of environmental microorganisms, this adaptation to humans is an opportunity to understand the genetic and metabolic selective mechanisms involved in this process. The previously reported genomic organization was incomplete, as data from animal strains were underrepresented. We added the missing piece of the puzzle with whole-genome sequencing of 93 strains of animal origin. Beyond describing the phylogenetic organization, we confirmed the genetic diversity of the Smc, which could not be estimated through routine phenotype- or matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF)-based laboratory tests. Animals strains seem to play a key role in the diversity of Smc and could act as a reservoir for mobile resistance genes. Some genogroups seem to be associated with particular hosts; the genetic support of this association and the role of the determinants/corresponding genes need to be explored.

Genomic Epidemiology of Complex, Multispecies, Plasmid-Borne bla KPC Carbapenemase in Enterobacterales in the United Kingdom from 2009 to 2014

Carbapenem resistance in Enterobacterales is a public health threat. Klebsiella pneumoniae carbapenemase (encoded by alleles of the bla _KPC family) is one of the most common transmissible carbapenem resistance mechanisms worldwide. The dissemination of bla _KPC historically has been associated with distinct K. pneumoniae lineages (clonal group 258 [CG258]), a particular plasmid family (pKpQIL), and a composite transposon (Tn4401). In the United Kingdom, bla _KPC has represented a large-scale, persistent management challenge for some hospitals, particularly in North West England. The dissemination of bla _KPC has evolved to be polyclonal and polyspecies, but the genetic mechanisms underpinning this evolution have not been elucidated in detail; this study used short-read whole-genome sequencing of 604 bla _KPC-positive isolates (Illumina) and long-read assembly (PacBio)/polishing (Illumina) of 21 isolates for characterization. We observed the dissemination of bla _KPC (predominantly bla _KPC-2; 573/604 [95%] isolates) across eight species and more than 100 known sequence types. Although there was some variation at the transposon level (mostly Tn4401a, 584/604 [97%] isolates; predominantly with ATTGA-ATTGA target site duplications, 465/604 [77%] isolates), bla _KPC spread appears to have been supported by highly fluid, modular exchange of larger genetic segments among plasmid populations dominated by IncFIB (580/604 isolates), IncFII (545/604 isolates), and IncR (252/604 isolates) replicons. The subset of reconstructed plasmid sequences (21 isolates, 77 plasmids) also highlighted modular exchange among non-bla _KPC and bla _KPC plasmids and the common presence of multiple replicons within bla _KPC plasmid structures (>60%). The substantial genomic plasticity observed has important implications for our understanding of the epidemiology of transmissible carbapenem resistance in Enterobacterales for the implementation of adequate surveillance approaches and for control.

Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multireplicon and Extensive Diversity of IncF Plasmids

Plasmids are genetic elements that enable rapid adaptation and evolution by transferring genes conferring selective advantages to their hosts. Conjugative plasmids are predominantly responsible for the global dissemination of antimicrobial resistance, representing an important threat to global health. As the number of plasmid sequences grows exponentially, it becomes critical to depict the global diversity and decipher the distribution of circulating plasmids in the bacterial community. To this end, we created COMPASS, a novel and comprehensive database compiling 12,084 complete plasmids with associated metadata from 1571 distinct species isolated worldwide over more than 100 years. The curation of the database allowed us to identify identical plasmids across different bacteria revealing mainly intraspecies dissemination and rare cases of horizontal transmission. We outlined and analyzed all relevant features, plasmid properties, host range and characterized their replication and mobilization systems. After an exhaustive comparison of PlasmidFinder and MOB-typer, the MOB-typer-based analysis revealed that the current knowledge embedded in the current typing schemes fails to classify all the plasmid sequences collected in COMPASS. We were able to categorize 6828 and 5229 plasmids by replicon and MOB typing, respectively, mostly associated with Proteobacteria and Firmicutes. We then searched for the presence of multiple core genes involved in replication and propagation. Our results showed that 2403 plasmids carried multiple replicons that were distributed in 206 bacterial species. The co-integration of replicon types from different incompatibility (Inc) groups is an adaptive mechanism, which plays an important role in plasmid survival and dissemination by extending their host range. Our results highlight the crucial role of IncF alleles (present in 56% of all multireplicons) and revealed that IncH, IncR, and IncU replicons were also frequently carried in multireplicons. Here, we provided a comprehensive picture of the different IncF subtypes by identifying 20 different profiles in 849 IncF multireplicons, which were mostly associated with Enterobacteriaceae. These results could provide the basis for a novel IncF plasmid nomenclature based on different allelic profiles.

Environmental and genetic determinants of plasmid mobility in pathogenic Escherichia coli

Plasmids are key vehicles of horizontal gene transfer (HGT), mobilizing antibiotic resistance, virulence, and other traits among bacterial populations. The environmental and genetic forces that drive plasmid transfer are poorly understood, however, due to the lack of definitive quantification coupled with genomic analysis. Here, we integrate conjugative phenotype with plasmid genotype to provide quantitative analysis of HGT in clinical Escherichia coli pathogens. We find a substantial proportion of these pathogens (>25%) able to readily spread resistance to the most common classes of antibiotics. Antibiotics of varied modes of action had less than a 5-fold effect on conjugation efficiency in general, with one exception displaying 31-fold promotion upon exposure to macrolides and chloramphenicol. In contrast, genome sequencing reveals plasmid incompatibility group strongly correlates with transfer efficiency. Our findings offer new insights into the determinants of plasmid mobility and have implications for the development of treatments that target HGT.

The Transformation of Reference Microbiology Methods and Surveillance for Salmonella With the Use of Whole Genome Sequencing in England and Wales

The use of whole genome sequencing (WGS) as a method for supporting outbreak investigations, studying Salmonella microbial populations and improving understanding of pathogenicity has been well-described (1–3). However, performing WGS on a discrete dataset does not pose the same challenges as implementing WGS as a routine, reference microbiology service for public health surveillance. Challenges include translating WGS data into a useable format for laboratory reporting, clinical case management, Salmonella surveillance, and outbreak investigation as well as meeting the requirement to communicate that information in an understandable and universal language for clinical and public health action. Public Health England have been routinely sequencing all referred presumptive Salmonella isolates since 2014 which has transformed our approach to reference microbiology and surveillance. Here we describe an overview of the integrated methods for cross-disciplinary working, describe the challenges and provide a perspective on how WGS has impacted the laboratory and surveillance processes in England and Wales.

Antibiotic-Induced, Increased Conjugative Transfer Is Common to Diverse Naturally Occurring ESBL Plasmids in Escherichia coli

Previously, we showed that cefotaxime (CTX) exposure increases conjugative transfer of a bla CTX-M- 1 encoding IncI1 plasmid (IncI1/pST49/CTX-M-1) in Escherichia coli in a SOS-independent manner. This study aimed at investigating whether the observation was unique for that plasmid/strain/antibiotic combination or whether antibiotic-induced plasmid transfer (PT) is a more general phenomenon among plasmids in E. coli. Whole genome sequences of 25 E. coli strains were analyzed to identify different extended spectrum beta-lactamases (ESBL) plasmids enabling selection of a diverse collection of plasmids. Experiments were performed following exposure of these strains to 1/2 minimal inhibitory concentration (MIC) of CTX, ampicillin (AMP), or ciprofloxacin (CIP) before conjugation experiments. The frequency of PT was measured and compared to that of donors not exposed to antibiotics. Reverse-transcribed-quantitative real time polymerase chain reaction (RT-qPCR) was used to measure mRNA levels of five PT genes and two SOS response genes in donors exposed to antibiotics. The PT of eight strains (30.8% of strains tested) with IncI1/pST7/CTX-M-1, IncI1/pST49/CTX-M-1, IncI1/pST3/CTX-M-1, IncI1/pST293/CTX-M-1, IncI1/pST295/CTX-M-1, IncI1/pST16/CTX-M-55, and IncFII/CTX-M-14 (n = 2) plasmids was significantly increased following antibiotic exposure. CTX increased PT in all of these eight strain/plasmid combinations, AMP and CIP increased the PT in six and three strains, respectively. RT-qPCR showed that PT genes were up-regulated in the presence of the three antibiotics, whereas SOS-response genes were up-regulated only following CIP exposure. Our findings reveal that antibiotics can increase PT in E. coli strains with various ESBL plasmids. Thus, antibiotic-induced conjugative transfer of ESBL plasmids appears to be a common phenomenon in E. coli, having important implications for assessing the risks of antibiotic use.

PlaScope: a targeted approach to assess the plasmidome from genome assemblies at the species level

Plasmid prediction may be of great interest when studying bacteria of medical importance such as Enterobacteriaceae as well as Staphylococcus aureus or Enterococcus. Indeed, many resistance and virulence genes are located on such replicons with major impact in terms of pathogenicity and spreading capacities. Beyond strain outbreak, plasmid outbreaks have been reported in particular for some extended-spectrum beta-lactamase- or carbapenemase-producing Enterobacteriaceae. Several tools are now available to explore the ‘plasmidome’ from whole-genome sequences with various approaches, but none of them are able to combine high sensitivity and specificity. With this in mind, we developed PlaScope, a targeted approach to recover plasmidic sequences in genome assemblies at the species or genus level. Based on Centrifuge, a metagenomic classifier, and a custom database containing complete sequences of chromosomes and plasmids from various curated databases, PlaScope classifies contigs from an assembly according to their predicted location. Compared to other plasmid classifiers, PlasFlow and cBar, it achieves better recall (0.87), specificity (0.99), precision (0.96) and accuracy (0.98) on a dataset of 70 genomes of Escherichia coli containing plasmids. In a second part, we identified 20 of the 21 chromosomal integrations of the extended-spectrum beta-lactamase coding gene in a clinical dataset of E. coli strains. In addition, we predicted virulence gene and operon locations in agreement with the literature. We also built a database for Klebsiella and correctly assigned the location for the majority of resistance genes from a collection of 12 Klebsiella pneumoniae strains. Similar approaches could also be developed for other well-characterized bacteria.

Identification of a Novel Plasmid Lineage Associated With the Dissemination of Metallo-β-Lactamase Genes Among Pseudomonads

Acquisition of metallo-β-lactamases (MBLs) represents one of most relevant resistance mechanisms to all β-lactams, including carbapenems, ceftolozane and available β-lactamase inhibitors, in Pseudomonas spp. VIM-type enzymes are the most common acquired MBLs in Pseudomonas aeruginosa and, to a lesser extent, in other Pseudomonas species. Little is known about the acquisition dynamics of these determinants, that are usually carried on integrons embedded into chromosomal mobile genetic elements. To date, few MBL-encoding plasmids have been described in Pseudomonas spp., and their diversity and role in the dissemination of these MBLs remains largely unknown. Here we report on the genetic features of the VIM-1-encoding plasmid pMOS94 from P. mosselii AM/94, the earliest known VIM-1-producing strain, and of related elements involved in dissemination of MBL. Results of plasmid DNA sequencing showed that pMOS94 had a modular organization, consisting of backbone modules associated with replication, transfer and antibiotic resistance. Plasmid pMOS94, although not typable according to the PBRT scheme, was classifiable either in MOB_F11 or MPF_T plasmid families. The resistance region included the class I integron In70, carrying bla _{V IM-1}, in turn embedded in a defective Tn402-like transposon. Comparison with pMOS94-like elements led to the identification of a defined plasmid lineage circulating in different Pseudomonas spp. of clinical and environmental origin and spreading different MBL-encoding genes, including bla _IMP-63, bla _BIM, and bla _{V IM}-type determinants. Genetic analysis revealed that this plasmid lineage likely shared a common ancestor and had evolved through the acquisition and recombination of different mobile elements, including the MBL-encoding transposons. Our findings provide new insights about the genetic diversity of MBL-encoding plasmids circulating among Pseudomonas spp., potentially useful for molecular epidemiology purposes, and revealed the existence and persistence of a successful plasmid lineage over a wide spatio-temporal interval, spanning over five different countries among two continents and over 20-years.

Illumina short-read and MinION long-read WGS to characterize the molecular epidemiology of an NDM-1 Serratia marcescens outbreak in Romania

Background and Objectives

Serratia marcescens is an emerging nosocomial pathogen, and the carbapenemase bla_NDM has been reported in several surveys in Romania. We aimed to investigate the molecular epidemiology of S. marcescens in two Romanian hospitals over 2010–15, including a neonatal NDM-1 S. marcescens outbreak.

Methods

Isolates were sequenced using Illumina technology together with carbapenem-non-susceptible NDM-1-positive and NDM-1-negative Klebsiella pneumoniae and Enterobacter cloacae to provide genomic context. A subset was sequenced with MinION to fully resolve NDM-1 plasmid structures. Resistance genes, plasmid replicons and ISs were identified in silico for all isolates; an annotated phylogeny was reconstructed for S. marcescens. Fully resolved study NDM-1 plasmid sequences were compared with the most closely related publicly available NDM-1 plasmid reference.

Results

44/45 isolates were successfully sequenced (S. marcescens, n = 33; K. pneumoniae, n = 7; E. cloacae, n = 4); 10 with MinION. The S. marcescens phylogeny demonstrated several discrete clusters of NDM-1-positive and -negative isolates. All NDM-1-positive isolates across species harboured a pKOX_NDM1-like plasmid; more detailed comparisons of the plasmid structures demonstrated a number of differences, but highlighted the largely conserved plasmid backbones across species and hospital sites.

Conclusions

The molecular epidemiology is most consistent with the importation of a pKOX_NDM1-like plasmid into Romania and its dissemination amongst K. pneumoniae/E. cloacae and subsequently S. marcescens across hospitals. The data suggested multiple acquisitions of this plasmid by S. marcescens in the two hospitals studied; transmission events within centres, including a large outbreak on the Targu Mures neonatal unit; and sharing of the pKOX_NDM1-like plasmid between species within outbreaks.

Phylogenetic barriers to horizontal transfer of antimicrobial peptide resistance genes in the human gut microbiota

The human gut microbiota has adapted to the presence of antimicrobial peptides (AMPs), which are ancient components of immune defence. Despite its medical importance, it has remained unclear whether AMP resistance genes in the gut microbiome are available for genetic exchange between bacterial species. Here, we show that AMP resistance and antibiotic resistance genes differ in their mobilization patterns and functional compatibilities with new bacterial hosts. First, whereas AMP resistance genes are widespread in the gut microbiome, their rate of horizontal transfer is lower than that of antibiotic resistance genes. Second, gut microbiota culturing and functional metagenomics have revealed that AMP resistance genes originating from phylogenetically distant bacteria have only a limited potential to confer resistance in Escherichia coli, an intrinsically susceptible species. Taken together, functional compatibility with the new bacterial host emerges as a key factor limiting the genetic exchange of AMP resistance genes. Finally, our results suggest that AMPs induce highly specific changes in the composition of the human microbiota, with implications for disease risks.

One Health Genomic Surveillance of Escherichia coli Demonstrates Distinct Lineages and Mobile Genetic Elements in Isolates from Humans versus Livestock

Livestock have been proposed as a reservoir for drug-resistant Escherichia coli that infect humans. We isolated and sequenced 431 E. coli isolates (including 155 extended-spectrum β-lactamase [ESBL]-producing isolates) from cross-sectional surveys of livestock farms and retail meat in the East of England. These were compared with the genomes of 1,517 E. coli bacteria associated with bloodstream infection in the United Kingdom. Phylogenetic core genome comparisons demonstrated that livestock and patient isolates were genetically distinct, suggesting that E. coli causing serious human infection had not directly originated from livestock. In contrast, we observed highly related isolates from the same animal species on different farms. Screening all 1,948 isolates for accessory genes encoding antibiotic resistance revealed 41 different genes present in variable proportions in human and livestock isolates. Overall, we identified a low prevalence of shared antimicrobial resistance genes between livestock and humans based on analysis of mobile genetic elements and long-read sequencing. We conclude that within the confines of our sampling framework, there was limited evidence that antimicrobial-resistant pathogens associated with serious human infection had originated from livestock in our region.IMPORTANCE The increasing prevalence of E. coli bloodstream infections is a serious public health problem. We used genomic epidemiology in a One Health study conducted in the East of England to examine putative sources of E. coli associated with serious human disease. E. coli from 1,517 patients with bloodstream infections were compared with 431 isolates from livestock farms and meat. Livestock-associated and bloodstream isolates were genetically distinct populations based on core genome and accessory genome analyses. Identical antimicrobial resistance genes were found in livestock and human isolates, but there was limited overlap in the mobile elements carrying these genes. Within the limitations of sampling, our findings do not support the idea that E. coli causing invasive disease or their resistance genes are commonly acquired from livestock in our region.

Reconsidering plasmid maintenance factors for computational plasmid design

Plasmids are genetic parasites of microorganisms. The genomes of naturally occurring plasmids are expected to be polished via natural selection to achieve long-term persistence in the microbial cell population. However, plasmid genomes are extremely diverse, and the rules governing plasmid genomes are not fully understood. Therefore, computationally designing plasmid genomes optimized for model and nonmodel organisms remains challenging. Here, we summarize current knowledge of the plasmid genome organization and the factors that can affect plasmid persistence, with the aim of constructing synthetic plasmids for use in gram-negative bacteria. Then, we introduce publicly available resources, plasmid data, and bioinformatics tools that are useful for computational plasmid design.

Contemporary IncI1 plasmids involved in the transmission and spread of antimicrobial resistance in Enterobacteriaceae

IncI1 has become one of the most common plasmid families in contemporary Enterobacteriaceae from both human and animal sources. In clinical epidemiology, this plasmid type ranks first as the confirmed vehicle of transmission of extended spectrum beta-lactamase and plasmid AmpC genes in isolates from food-producing animals. In this review, we describe the epidemiology and evolution of IncI1 plasmids and closely related IncIγ plasmids. We highlight the emergence of epidemic plasmids circulating among different bacterial hosts in geographically distant countries, and we address the phylogeny of the IncI1 and IncIγ family based on plasmid Multilocus Sequence Typing.

Hotspot mutations and ColE1 plasmids contribute to the fitness of Salmonella Heidelberg in poultry litter

Salmonella enterica subsp. enterica serovar Heidelberg (S. Heidelberg) is a clinically-important serovar linked to food-borne illness, and commonly isolated from poultry. Investigations of a large, multistate outbreak in the USA in 2013 identified poultry litter (PL) as an important extra-intestinal environment that may have selected for specific S. Heidelberg strains. Poultry litter is a mixture of bedding materials and chicken excreta that contains chicken gastrointestinal (GI) bacteria, undigested feed, feathers, and other materials of chicken origin. In this study, we performed a series of controlled laboratory experiments which assessed the microevolution of two S. Heidelberg strains (SH-2813 and SH-116) in PL previously used to raise 3 flocks of broiler chickens. The strains are closely related at the chromosome level, differing from the reference genome by 109 and 89 single nucleotide polymorphisms/InDels, respectively. Whole genome sequencing was performed on 86 isolates recovered after 0, 1, 7 and 14 days of microevolution in PL. Only strains carrying an IncX1 (37kb), 2 ColE1 (4 and 6kb) and 1 ColpVC (2kb) plasmids survived more than 7 days in PL. Competition experiments showed that carriage of these plasmids was associated with increased fitness. This increased fitness was associated with an increased copy number of IncX1 and ColE1 plasmids. Further, all Col plasmid-bearing strains had hotspot mutations in 37 loci on the chromosome and in 3 loci on the IncX1 plasmid. Additionally, we observed a decrease in susceptibility to tobramycin, kanamycin, gentamicin, neomycin and fosfomycin for Col plasmid-bearing strains. Our study demonstrates how positive selection from poultry litter can change the evolutionary path of S. Heidelberg.

Ordering the mob: Insights into replicon and MOB typing schemes from analysis of a curated dataset of publicly available plasmids

Plasmid typing can provide insights into the epidemiology and transmission of plasmid-mediated antibiotic resistance. The principal plasmid typing schemes are replicon typing and MOB typing, which utilize variation in replication loci and relaxase proteins respectively. Previous studies investigating the proportion of plasmids assigned a type by these schemes (‘typeability’) have yielded conflicting results; moreover, thousands of plasmid sequences have been added to NCBI in recent years, without consistent annotation to indicate which sequences represent complete plasmids. Here, a curated dataset of complete Enterobacteriaceae plasmids from NCBI was compiled, and used to assess the typeability and concordance of in silico replicon and MOB typing schemes. Concordance was assessed at hierarchical replicon type resolutions, from replicon family-level to plasmid multilocus sequence type (pMLST)-level, where available. We found that 85% and 65% of the curated plasmids could be replicon and MOB typed, respectively. Overall, plasmid size and the number of resistance genes were significant independent predictors of replicon and MOB typing success. We found some degree of non-concordance between replicon families and MOB types, which was only partly resolved when partitioning plasmids into finer-resolution groups (replicon and pMLST types). In some cases, non-concordance was attributed to ambiguous boundaries between MOBP and MOBQ types; in other cases, backbone mosaicism was considered a more plausible explanation. β-lactamase resistance genes tended not to show fidelity to a particular plasmid type, though some previously reported associations were supported. Overall, replicon and MOB typing schemes are likely to continue playing an important role in plasmid analysis, but their performance is constrained by the diverse and dynamic nature of plasmid genomes.

•

92% of clinically-relevant plasmids could be replicon typed, compared with reports of 100% typeability in 2014.

•

Replicon and MOB typing were partly concordant; partitioning plasmids into finer-resolution groups increased concordance.

•

Overlap between MOBP and MOBQ relaxase families complicates assignment of MOBP and MOBQ types.

•

Generally, resistance genes showed low fidelity towards particular plasmid backbones.

•

PacBio sequencing has driven increased availability of complete plasmid sequences, but retrieved datasets require curation.