MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies

Mobility of β-lactam resistance under ampicillin treatment in gut microbiota suffering from pre-disturbance

Ingestion of food- or waterborne antibiotic-resistant bacteria may lead to dissemination of antibiotic resistance genes (ARGs) in the gut microbiota. The gut microbiota often suffers from various disturbances. It is not clear whether and how disturbed microbiota may affect ARG mobility under antibiotic treatments. For proof of concept, in the presence or absence of streptomycin pre-treatment, mice were inoculated orally with a β-lactam-susceptible Salmonella enterica serovar Heidelberg clinical isolate (recipient) and a β-lactam resistant Escherichia coli O80:H26 isolate (donor) carrying a blaCMY-2 gene on an IncI2 plasmid. Immediately following inoculation, mice were treated with or without ampicillin in drinking water for 7 days. Faeces were sampled, donor, recipient and transconjugant were enumerated, blaCMY-2 abundance was determined by quantitative PCR, faecal microbial community composition was determined by 16S rRNA amplicon sequencing and cecal samples were observed histologically for evidence of inflammation. In faeces of mice that received streptomycin pre-treatment, the donor abundance remained high, and the abundance of S. Heidelberg transconjugant and the relative abundance of Enterobacteriaceae increased significantly during the ampicillin treatment. Co-blooming of the donor, transconjugant and commensal Enterobacteriaceae in the inflamed intestine promoted significantly (P<0.05) higher and possibly wider dissemination of the blaCMY-2 gene in the gut microbiota of mice that received the combination of streptomycin pre-treatment and ampicillin treatment (Str-Amp) compared to the other mice. Following cessation of the ampicillin treatment, faecal shedding of S. Heidelberg transconjugant persisted much longer from mice in the Str-Amp group compared to the other mice. In addition, only mice in the Str-Amp group shed a commensal E. coli O2:H6 transconjugant, which carries three copies of the blaCMY-2 gene, one on the IncI2 plasmid and two on the chromosome. The findings highlight the significance of pre-existing gut microbiota for ARG dissemination and persistence during and following antibiotic treatments of infectious diseases.

RFPlasmid: predicting plasmid sequences from short-read assembly data using machine learning

Antimicrobial-resistance (AMR) genes in bacteria are often carried on plasmids and these plasmids can transfer AMR genes between bacteria. For molecular epidemiology purposes and risk assessment, it is important to know whether the genes are located on highly transferable plasmids or in the more stable chromosomes. However, draft whole-genome sequences are fragmented, making it difficult to discriminate plasmid and chromosomal contigs. Current methods that predict plasmid sequences from draft genome sequences rely on single features, like k-mer composition, circularity of the DNA molecule, copy number or sequence identity to plasmid replication genes, all of which have their drawbacks, especially when faced with large single-copy plasmids, which often carry resistance genes. With our newly developed prediction tool RFPlasmid, we use a combination of multiple features, including k-mer composition and databases with plasmid and chromosomal marker proteins, to predict whether the likely source of a contig is plasmid or chromosomal. The tool RFPlasmid supports models for 17 different bacterial taxa, including Campylobacter, Escherichia coli and Salmonella, and has a taxon agnostic model for metagenomic assemblies or unsupported organisms. RFPlasmid is available both as a standalone tool and via a web interface.

Evaluation of whole-genome sequencing-based subtyping methods for the surveillance of Shigella spp. and the confounding effect of mobile genetic elements in long-term outbreaks

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232  Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST<SNVPhyl<cgMLST<wgMLST (range: 0.970–1.000). The genetic differences observed for non-MSM-related Shigella spp. outbreaks respect the standard 0–10 allele/SNV guideline; however, mobile genetic element (MGE)-encoded loci caused inflated genetic variation and discrepant phylogenies for prolonged MSM-related S. sonnei outbreaks via wgMLST. The S. sonnei correlation coefficients of wgMLST were also the lowest at 0.680, 0.703 and 0.712 for SNVPhyl, coreMLST and cgMLST, respectively. Plasmid maintenance, mobilization and conjugation-associated genes were found to be the main source of genetic distance inflation in addition to prophage-related genes. Duplicated alleles arising from the repeated nature of IS elements were also responsible for many false cg/wgMLST differences. The coreMLST approach was shown to be the most robust, followed by SNVPhyl and wgMLST for inter-laboratory comparability. Our results highlight the need for validating species-specific subtyping methods based on microbial genome plasticity and outbreak dynamics in addition to the importance of filtering confounding MGEs for cluster detection.

Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification

Command-line annotation software tools have continuously gained popularity compared to centralized online services due to the worldwide increase of sequenced bacterial genomes. However, results of existing command-line software pipelines heavily depend on taxon-specific databases or sufficiently well annotated reference genomes. Here, we introduce Bakta, a new command-line software tool for the robust, taxon-independent, thorough and, nonetheless, fast annotation of bacterial genomes. Bakta conducts a comprehensive annotation workflow including the detection of small proteins taking into account replicon metadata. The annotation of coding sequences is accelerated via an alignment-free sequence identification approach that in addition facilitates the precise assignment of public database cross-references. Annotation results are exported in GFF3 and International Nucleotide Sequence Database Collaboration (INSDC)-compliant flat files, as well as comprehensive JSON files, facilitating automated downstream analysis. We compared Bakta to other rapid contemporary command-line annotation software tools in both targeted and taxonomically broad benchmarks including isolates and metagenomic-assembled genomes. We demonstrated that Bakta outperforms other tools in terms of functional annotations, the assignment of functional categories and database cross-references, whilst providing comparable wall-clock runtimes. Bakta is implemented in Python 3 and runs on MacOS and Linux systems. It is freely available under a GPLv3 license at https://github.com/oschwengers/bakta. An accompanying web version is available at https://bakta.computational.bio.

A genomic epidemiological study shows that prevalence of antimicrobial resistance in Enterobacterales is associated with the livestock host, as well as antimicrobial usage

Enterobacterales from livestock are potentially important reservoirs for antimicrobial resistance (AMR) to pass through the food chain to humans, thereby increasing the AMR burden and affecting our ability to tackle infections. In this study 168 isolates from four genera of the order Enterobacterales , primarily Escherichia coli , were purified from livestock (cattle, pigs and sheep) faeces from 14 farms in the United Kingdom. Their genomes were resolved using long- and short-read sequencing to analyse AMR genes and their genetic context, as well as to explore the relationship between AMR burden and on-farm antimicrobial usage (AMU), in the three months prior to sampling. Although E. coli isolates were genomically diverse, phylogenetic analysis using a core-genome SNP tree indicated pig isolates to generally be distinct from sheep isolates, with cattle isolates being intermediates. Approximately 28 % of isolates harboured AMR genes, with the greatest proportion detected in pigs, followed by cattle then sheep; pig isolates also harboured the highest number of AMR genes per isolate. Although 90 % of sequenced isolates harboured diverse plasmids, only 11 % of plasmids (n =58 out of 522) identified contained AMR genes, with 91 % of AMR plasmids being from pig, 9 % from cattle and none from sheep isolates; these results indicated that pigs were a principle reservoir of AMR genes harboured by plasmids and likely to be involved in their horizontal transfer. Significant associations were observed between AMU (mg kg−1) and AMR. As both the total and the numbers of different antimicrobial classes used on-farm increased, the risk of multi-drug resistance (MDR) in isolates rose. However, even when AMU on pig farms was comparatively low, pig isolates had increased likelihood of being MDR; harbouring relatively more resistances than those from other livestock species. Therefore, our results indicate that AMR prevalence in livestock is not only influenced by recent AMU on-farm but also livestock-related factors, which can influence the AMR burden in these reservoirs and its plasmid mediated transmission.

Infective endocarditis caused by Enterobacteriaceae: phenotypic and molecular characterization of Escherichia coli and Klebsiella pneumoniae in Rio de Janeiro, Brazil

The etiological agent for infective endocarditis (IE), a life-threatening disease, is usually gram-positive bacteria. However, gram-negative bacteria can rarely cause IE and 4% of cases are associated with morbidity and mortality. This study aimed to characterize Escherichia coli and Klebsiella pneumoniae isolates from the blood of patients with IE. The characteristics of blood isolates were compared with those of urinary isolates from patients with urinary tract infections (UTIs). The results of this study revealed that K. pneumoniae isolates from patients with IE were phylogenetically related to those from patients with UTI. Additionally, the resistance phenotype, resistance gene, virulence gene, and plasmid profiles were similar between the blood and urinary isolates. The isolates belonging to the sequence types (STs) 76, 36, 101 (K. pneumoniae), and 69 (E. coli) are reported to be associated with drug resistance. The Enterobacteriaceae isolates from patients with IE did not produce extended-spectrum β-lactamase or carbapenemase. Additionally, this study investigated the virulence phenotype, biofilm formation ability, and the ability to adhere to the epithelial cells in vitro of the isolates. The isolates from patients with IE exhibited weaker biofilm formation ability than the urinary isolates. All isolates from patients with IE could adhere to the renal epithelial cells. However, three isolates from patients with UTIs could not adhere to the epithelial cells. The closely related K. pneumoniae isolates (648, KP1, KP2, KP3, and KP4) could not form biofilms or adhere to the epithelial cells. In summary, the molecular analysis revealed that the genetic characteristics of IE-causing K. pneumoniae and E. coli were similar to those of UTI-causing isolates. These isolates belonged to the STs that are considered treatable. Genetically similar isolates did not exhibit the same virulence phenotype. Thus, these non-hypervirulent clones must be monitored as they can cause complex infections in susceptible hosts.

The complete genome of 2,6-dichlorobenzamide (BAM) degrader Aminobacter sp. MSH1 suggests a polyploid chromosome, phylogenetic reassignment, and functions of plasmids

Aminobacter sp. MSH1 (CIP 110285) can use the pesticide dichlobenil and its recalcitrant transformation product, 2,6-dichlorobenzamide (BAM), as sole source of carbon, nitrogen, and energy. The concentration of BAM in groundwater often exceeds the threshold limit for drinking water, requiring additional treatment in drinking water treatment plants or closure of the affected abstraction wells. Biological treatment with MSH1 is considered a potential sustainable alternative to remediate BAM-contamination in drinking water production. We present the complete genome of MSH1, which was determined independently in two institutes at Aarhus University and KU Leuven. Divergences were observed between the two genomes, i.e. one of them lacked four plasmids compared to the other. Besides the circular chromosome and the two previously described plasmids involved in BAM catabolism, pBAM1 and pBAM2, the genome of MSH1 contained two megaplasmids and three smaller plasmids. The MSH1 substrain from KU Leuven showed a reduced genome lacking a megaplasmid and three smaller plasmids and was designated substrain MK1, whereas the Aarhus variant with all plasmids was designated substrain DK1. A plasmid stability experiment indicate that substrain DK1 may have a polyploid chromosome when growing in R2B medium with more chromosomes than plasmids per cell. Finally, strain MSH1 is reassigned as Aminobacter niigataensis MSH1.

Occurrence of Colibacillosis in Broilers and Its Relationship With Avian Pathogenic Escherichia coli (APEC) Population Structure and Molecular Characteristics

Avian pathogenic Escherichia coli (APEC) causes colibacillosis, the disease with the highest economic loss for the broiler industry. However, studies focusing on the prevalence and population structure of APEC in the broiler production pyramid are scarce. Here, we used genotyping and serotyping data to elucidate the APEC population structure and its changes in different broiler production stages along with whole-genome sequencing (WGS) in a subset of APEC isolates to determine transmission patterns amongst dominant APEC sequence types (STs) and characterize them in detail. Comparison of genotypes encountered in both APEC and avian fecal E. coli (AFEC) provided further insights. Overall, APEC-related mortality, as the proportion of the total sampled mortality in the broiler production, was high (35%), while phylogroup C and serogroup O78 were predominant amongst APEC isolates. We found a low (34.0%) and high (53.3%) incidence of colibacillosis in chicks and end-cycle broilers, respectively, which may be related to a shift in APEC genotypes, suggesting a trend from commensalism to pathogenicity across different broiler production stages. Despite considerable APEC genotypic diversity, there was substantial genotype overlap (40.9%, overall) over the production stages and convergence of STs to the four clusters. Within these clusters, WGS data provided evidence of clonal transmission events and revealed an enriched virulence and resistance APEC repertoire. More specifically, sequenced APEC were assigned to defined pathotypes based on their virulence gene content while the majority (86%) was genotypically multi-drug resistant. Interestingly, WGS-based phylogeny showed that a subset of APEC, which are cephalosporin-resistant, may originate directly from cephalosporin-resistant AFEC. Finally, exploration of the APEC plasmidome indicated that the small fraction of the APEC virulome carried by IncF plasmids is pivotal for the manifestation of the APEC pathotype; thus, plasmid exchange can promote pathogenicity in strains that are at the edge of the commensal and pathogenic states.

Flanker: a tool for comparative genomics of gene flanking regions

Analysing the flanking sequences surrounding genes of interest is often highly relevant to understanding the role of mobile genetic elements (MGEs) in horizontal gene transfer, particular for antimicrobial-resistance genes. Here, we present Flanker, a Python package that performs alignment-free clustering of gene flanking sequences in a consistent format, allowing investigation of MGEs without prior knowledge of their structure. These clusters, known as 'flank patterns' (FPs), are based on Mash distances, allowing for easy comparison of similarity across sequences. Additionally, Flanker can be flexibly parameterized to fine-tune outputs by characterizing upstream and downstream regions separately, and investigating variable lengths of flanking sequence. We apply Flanker to two recent datasets describing plasmid-associated carriage of important carbapenemase genes (blaOXA-48 and blaKPC-2/3) and show that it successfully identifies distinct clusters of FPs, including both known and previously uncharacterized structural variants. For example, Flanker identified four Tn4401 profiles that could not be sufficiently characterized using TETyper or MobileElementFinder, demonstrating the utility of Flanker for flanking-gene characterization. Similarly, using a large (n=226) European isolate dataset, we confirm findings from a previous smaller study demonstrating association between Tn1999.2 and blaOXA-48 upregulation and demonstrate 17 FPs (compared to the 5 previously identified). More generally, the demonstration in this study that FPs are associated with geographical regions and antibiotic-susceptibility phenotypes suggests that they may be useful as epidemiological markers. Flanker is freely available under an MIT license at https://github.com/wtmatlock/flanker.

Whole genome analysis unveils genetic diversity and potential virulence determinants in Vibrio parahaemolyticus associated with disease outbreak among cultured Litopenaeus vannamei (Pacific white shrimp) in India

Vibrio parahaemolyticus has caused widespread mortality in Indian shrimp aquaculture in recent years. However, there are insufficient genome data for the isolates from Indian shrimp vibriosis to analyze genetic diversity and track the acquisition of genetic features that could be involved in virulence and fitness. In this study, we have performed genome analysis of V. parahaemolyticus isolated from moribund shrimps collected from shrimp farms along coastal Karnataka, India, for better understanding of their diversity and virulence. Five newly sequenced genomes of V. parahaemolyticus along with 40 genomes retrieved from NCBI were subjected to comparative genome analysis. The sequenced genomes had an overall genome size of 5.2 Mb. MLST analysis and core genome phylogenomic analysis revealed considerable genetic diversity among the isolates obtained from the moribund shrimps. Interestingly, none of the V. parahaemolyticus isolates possessed the classical features (PirAB) of the strains associated with Acute Hepatopancreatic Necrosis Disease (AHPND). This study also revealed the presence of multiple virulence attributes, including ZOT, ACE and RTX toxins, secretion systems, and mobile genetic elements. The findings of this study provide insights into the possible transition of an environmental V. parahaemolyticus to emerge as pathogens of aquaculture species by increasing its virulence and host adaptation. Future studies focusing on continuous genomic surveillance of V. parahaemolyticus are required to study the evolution and transmission of new variants in shrimp aquaculture, as well as to design and implement biosecurity programs to prevent disease outbreaks.

Universal whole-sequence-based plasmid typing and its utility to prediction of host range and epidemiological surveillance

Bacterial plasmids play a large role in allowing bacteria to adapt to changing environments and can pose a significant risk to human health if they confer virulence and antimicrobial resistance (AMR). Plasmids differ significantly in the taxonomic breadth of host bacteria in which they can successfully replicate, this is commonly referred to as 'host range' and is usually described in qualitative terms of 'narrow' or 'broad'. Understanding the host range potential of plasmids is of great interest due to their ability to disseminate traits such as AMR through bacterial populations and into human pathogens. We developed the MOB-suite to facilitate characterization of plasmids and introduced a whole-sequence-based classification system based on clustering complete plasmid sequences using Mash distances (https://github.com/phac-nml/mob-suite). We updated the MOB-suite database from 12 091 to 23 671 complete sequences, representing 17 779 unique plasmids. With advances in new algorithms for rapidly calculating average nucleotide identity (ANI), we compared clustering characteristics using two different distance measures - Mash and ANI - and three clustering algorithms on the unique set of plasmids. The plasmid nomenclature is designed to group highly similar plasmids together that are unlikely to have multiple representatives within a single cell. Based on our results, we determined that clusters generated using Mash and complete-linkage clustering at a Mash distance of 0.06 resulted in highly homogeneous clusters while maintaining cluster size. The taxonomic distribution of plasmid biomarker sequences for replication and relaxase typing, in combination with MOB-suite whole-sequence-based clusters have been examined in detail for all high-quality publicly available plasmid sequences. We have incorporated prediction of plasmid replication host range into the MOB-suite based on observed distributions of these sequence features in combination with known plasmid hosts from the literature. Host range is reported as the highest taxonomic rank that covers all of the plasmids which share replicon or relaxase biomarkers or belong to the same MOB-suite cluster code. Reporting host range based on these criteria allows for comparisons of host range between studies and provides information for plasmid surveillance.

Characteristics of fosA-carrying plasmids in E. coli and Klebsiella spp. isolates originating from food and environmental samples

OBJECTIVES

Fosfomycin is an important antibiotic for the treatment of MDR Enterobacteriaceae infections. High susceptibility rates are, however, threatened by the spread of plasmids encoding fosfomycin-modifying enzymes. In this study, we sought to characterize the genetic context of fosA in plasmids from Escherichia coli and Klebsiella spp. isolates recovered from food, wastewater and surface water in Switzerland.

METHODS

E. coli and Klebsiella spp. isolates collected between 2012 and 2019 in Switzerland were screened for fosfomycin resistance. Presence of fosA was verified by PCR and sodium phosphonoformate (PPF) disc potentiation testing, and transferability was tested using conjugation assays. Whole-genome sequences including complete fosA-containing plasmids were determined using long- and short-read sequencing.

RESULTS

In 11 E. coli and two Klebsiella spp. isolates, high-level fosfomycin resistance was mediated by plasmids containing fosA3 (n = 12) or fosA8 (n = 1). Four isolates harboured a near-identical 45 kb IncN plasmid with fosA3, while replicon types varied in the remaining plasmids. The fosA genes were typically embedded in IS26-bounded transposition units and frequently located in the proximity of blaCTX-M transposition units.

CONCLUSIONS

Although fosfomycin resistance rates are currently low, the presence of fosA-encoding plasmids circulating in the Enterobacteriaceae population suggests that fosfomycin resistance may rapidly spread upon increased selection pressure. Transposition mobility of fosA and co-location on plasmids with other resistance genes may further promote its dissemination.

Genomic network analysis of environmental and livestock F-type plasmid populations

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum β-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.

Whole Genome-Based Characterization of Listeria monocytogenes Isolates Recovered From the Food Chain in South Africa

Listeria monocytogenes is an important foodborne pathogen which has the ability to adapt and survive in food and food processing facilities where it can persist for years. In this study, a total of 143 L. monocytogenes isolates in South Africa (SA) were characterized for their strain’s genetic relatedness, virulence profiles, stress tolerance and resistance genes associated with L. monocytogenes. The Core Genome Multilocus Sequence Typing (cgMLST) analysis revealed that the most frequent serogroups were IVb and IIa; Sequence Types (ST) were ST204, ST2, and ST1; and Clonal Complexes (CC) were CC204, CC1, and CC2. Examination of genes involved in adaptation and survival of L. monocytogenes in SA showed that ST1, ST2, ST121, ST204, and ST321 are well adapted in food processing environments due to the significant over-representation of Benzalkonium chloride (BC) resistance genes (bcrABC cassette, ermC, mdrL and Ide), stress tolerance genes (SSI-1 and SSI-2), Prophage (φ) profiles (LP_101, vB LmoS 188, vB_LmoS_293, and B054 phage), plasmids profiles (N1-011A, J1776, and pLM5578) and biofilm formation associated genes. Furthermore, the L. monocytogenes strains that showed hyper-virulent potential were ST1, ST2 and ST204, and hypo-virulent were ST121 and ST321 because of the presence and absence of major virulence factors such as LIPI-1, LIPI-3, LIPI-4 and the internalin gene family members including inlABCEFJ. The information provided in this study revealed that hyper-virulent strains ST1, ST2, and ST204 could present a major public health risk due to their association with meat products and food processing environments in SA.

PlasForest: a homology-based random forest classifier for plasmid detection in genomic datasets

Background

Plasmids are mobile genetic elements that often carry accessory genes, and are vectors for horizontal transfer between bacterial genomes. Plasmid detection in large genomic datasets is crucial to analyze their spread and quantify their role in bacteria adaptation and particularly in antibiotic resistance propagation. Bioinformatics methods have been developed to detect plasmids. However, they suffer from low sensitivity (i.e., most plasmids remain undetected) or low precision (i.e., these methods identify chromosomes as plasmids), and are overall not adapted to identify plasmids in whole genomes that are not fully assembled (contigs and scaffolds).

Results

We developed PlasForest, a homology-based random forest classifier identifying bacterial plasmid sequences in partially assembled genomes. Without knowing the taxonomical origin of the samples, PlasForest identifies contigs as plasmids or chromosomes with a F1 score of 0.950. Notably, it can detect 77.4% of plasmid contigs below 1 kb with 2.8% of false positives and 99.9% of plasmid contigs over 50 kb with 2.2% of false positives.

Conclusions

PlasForest outperforms other currently available tools on genomic datasets by being both sensitive and precise. The performance of PlasForest on metagenomic assemblies are currently well below those of other k-mer-based methods, and we discuss how homology-based approaches could improve plasmid detection in such datasets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04270-w.

On-Farm Anaerobic Digestion of Dairy Manure Reduces the Abundance of Antibiotic Resistance-Associated Gene Targets and the Potential for Plasmid Transfer

The present study investigated the impact of on-farm anaerobic digestion on the abundance of enteric bacteria, antibiotic resistance-associated gene targets, and the horizontal transfer potential of extended-spectrum β-lactamase (ESBL) genes. Samples of raw and digested manure were obtained from six commercial dairy farms in Ontario, Canada. Digestion significantly abated populations of viable coliforms in all six farms. Conjugative transfer of plasmids carrying β-lactamase genes from manure bacteria enriched overnight with buffered peptone containing 4 mg/liter cefotaxime into a β-lactam-sensitive green fluorescent protein (GFP)-labeled Escherichia coli recipient strain was evaluated in patch matings. Digestion significantly decreased the frequency of the horizontal transfer of ESBL genes. Twenty-five transconjugants were sequenced, revealing six distinct plasmids, ranging in size from 40 to 180 kb. A variety of ESBL genes were identified: bla_CTX-M-1, bla_CTX-M-14, bla_CTX-M-15, bla_CTX-M-27, bla_CTX-M-55, and bla_PER-1. bla_CTX-M-15 was the most prevalent ESBL gene detected on plasmids harbored by transconjugants. Various mobile genetic elements were found located proximal to resistance genes. Ten gene targets, including sul1, str(A), str(B), erm(B), erm(F), intI1, aadA, incW, bla_PSE, and bla_OXA-20, were quantified by quantitative PCR on a subset of 18 raw and 18 digested samples. Most targets were significantly more abundant in raw manure; however, erm(B) and erm(F) targets were more abundant in digested samples. Overall, on-farm digestion of dairy manure abated coliform bacteria, a number of antibiotic resistance-associated gene targets, and the potential for in vitro conjugation of plasmids conferring resistance to extended-spectrum β-lactams and other classes of antibiotics into E. coli CV601. IMPORTANCE Using livestock manure for fertilization can entrain antibiotic-resistant bacteria into soil. Manure on some dairy farms is anaerobically digested before being land applied. Recommending the widespread implementation of the practice should be founded on understanding the impact of this treatment on various endpoints of human health concern. Although lab-scale anaerobic treatments have shown potential for reducing the abundance of antibiotic resistance genes, there are very few data from commercial farms. Anaerobic digestion of manure on six dairy farms efficiently abated coliform bacteria, E. coli, and a majority of antibiotic resistance-associated gene targets. In addition, the conjugation potential of plasmids carrying ESBL genes into introduced E. coli strain CV601 was reduced. Overall, anaerobic digestion abated coliform bacteria, the genes that they carry, and the potential for ESBL-carrying plasmid transfer.

Conjugative plasmids interact with insertion sequences to shape the horizontal transfer of antimicrobial resistance genes

It is well established that plasmids play an important role in the dissemination of antimicrobial resistance (AMR) genes; however, little is known about the role of the underlying interactions between different plasmid categories and other mobile genetic elements (MGEs) in shaping the promiscuous spread of AMR genes. Here, we developed a tool designed for plasmid classification, AMR gene annotation, and plasmid visualization and found that most plasmid-borne AMR genes, including those localized on class 1 integrons, are enriched in conjugative plasmids. Notably, we report the discovery and characterization of a massive insertion sequence (IS)-associated AMR gene transfer network (245 combinations covering 59 AMR gene subtypes and 53 ISs) linking conjugative plasmids and phylogenetically distant pathogens, suggesting a general evolutionary mechanism for the horizontal transfer of AMR genes mediated by the interaction between conjugative plasmids and ISs. Moreover, our experimental results confirmed the importance of the observed interactions in aiding the horizontal transfer and expanding the genetic range of AMR genes within complex microbial communities.

The Spatiotemporal Dynamics and Microevolution Events That Favored the Success of the Highly Clonal Multidrug-Resistant Monophasic Salmonella Typhimurium Circulating in Europe

The European epidemic monophasic variant of Salmonella enterica serovar Typhimurium (S. 1,4,[5],12:i:-) characterized by the multi locus sequence type ST34 and the antimicrobial resistance ASSuT profile has become one of the most common serovars in Europe (EU) and the United States (US). In this study, we reconstructed the time-scaled phylogeny and evolution of this Salmonella in Europe. The epidemic S. 1,4,[5],12:i:- ST34 emerged in the 1980s by an acquisition of the Salmonella Genomic Island (SGI)-4 at the 3' end of the phenylalanine phe tRNA locus conferring resistance to copper and arsenic toxicity. Subsequent integration of the Tn21 transposon into the fljAB locus gave resistance to mercury toxicity and several classes of antibiotics used in food-producing animals (ASSuT profile). The second step of the evolution occurred in the 1990s, with the integration of mTmV and mTmV-like prophages carrying the perC and/or sopE genes involved in the ability to reduce nitrates in intestinal contents and facilitate the disruption of the junctions of the host intestinal epithelial cells. Heavy metals are largely used as food supplements or pesticide for cultivation of seeds intended for animal feed so the expansion of the epidemic S. 1,4,[5],12:i:- ST34 was strongly related to the multiple-heavy metal resistance acquired by transposons, integrative and conjugative elements and facilitated by the escape until 2011 from the regulatory actions applied in the control of S. Typhimurium in Europe. The genomic plasticity of the epidemic S. 1,4,[5],12:i:- was demonstrated in our study by the analysis of the plasmidome. We were able to identify plasmids harboring genes mediating resistance to phenicols, colistin, and fluoroquinolone and also describe for the first time in six of the analyzed genomes the presence of two plasmids (pERR1744967-1 and pERR2174855-2) previously described only in strains of enterotoxigenic Escherichia coli and E. fergusonii.

Fatal exudative dermatitis in island populations of red squirrels (Sciurus vulgaris): spillover of a virulent Staphylococcus aureus clone (ST49) from reservoir hosts

Fatal exudative dermatitis (FED) is a significant cause of death of red squirrels (Sciurus vulgaris) on the island of Jersey in the Channel Islands where it is associated with a virulent clone of Staphylococcus aureus, ST49. S. aureus ST49 has been found in other hosts such as small mammals, pigs and humans, but the dynamics of carriage and disease of this clone, or any other lineage in red squirrels, is currently unknown. We used whole-genome sequencing to characterize 228 isolates from healthy red squirrels on Jersey, the Isle of Arran (Scotland) and Brownsea Island (England), from red squirrels showing signs of FED on Jersey and the Isle of Wight (England) and a small number of isolates from other hosts. S. aureus was frequently carried by red squirrels on the Isle of Arran with strains typically associated with small ruminants predominating. For the Brownsea carriage, S. aureus was less frequent and involved strains associated with birds, small ruminants and humans, while for the Jersey carriage S. aureus was rare but ST49 predominated in diseased squirrels. By combining our data with publicly available sequences, we show that the S. aureus carriage in red squirrels largely reflects frequent but facile acquisitions of strains carried by other hosts sharing their habitat ('spillover'), possibly including, in the case of ST188, humans. Genome-wide association analysis of the ruminant lineage ST133 revealed variants in a small number of mostly bacterial-cell-membrane-associated genes that were statistically associated with squirrel isolates from the Isle of Arran, raising the possibility of specific adaptation to red squirrels in this lineage. In contrast there is little evidence that ST49 is a common carriage isolate of red squirrels and infection from reservoir hosts such as bank voles or rats, is likely to be driving the emergence of FED in red squirrels.

Genomics of Environmental Salmonella: Engaging Students in the Microbiology and Bioinformatics of Foodborne Pathogens

We have developed and implemented an undergraduate microbiology course in which students isolate, characterize, and perform whole genome assembly and analysis of Salmonella enterica from stream sediments and poultry litter. In the development of the course and over three semesters, successive teams of undergraduate students collected field samples and performed enrichment and isolation techniques specific for the detection of S. enterica. Eighty-eight strains were confirmed using standard microbiological methods and PCR of the invA gene. The isolates' genomes were Illumina-sequenced by the Center for Food Safety and Applied Nutrition at the FDA and the Virginia state Division of Consolidated Laboratory Services as part of the GenomeTrakr program. Students used GalaxyTrakr and other web- and non-web-based platforms and tools to perform quality control on raw and assembled sequence data, assemble, and annotate genomes, identify antimicrobial resistance and virulence genes, putative plasmids, and other mobile genetic elements. Strains with putative plasmid-borne antimicrobial resistance genes were further sequenced by students in our research lab using the Oxford Nanopore MinIONTM platform. Strains of Salmonella that were isolated include human infectious serotypes such as Typhimurium and Infantis. Over 31 of the isolates possessed antibiotic resistance genes, some of which were located on large, multidrug resistance plasmids. Plasmid pHJ-38, identified in a Typhimurium isolate, is an apparently self-transmissible 183 kb IncA/C2 plasmid that possesses multiple antimicrobial resistance and heavy-metal resistance genes. Plasmid pFHS-02, identified in an Infantis isolate, is an apparently self-transmissible 303 kb IncF1B plasmid that also possesses numerous heavy-metal and antimicrobial resistance genes. Using direct and indirect measures to assess student outcomes, results indicate that course participation contributed to cognitive gains in relevant content knowledge and research skills such as field sampling, molecular techniques, and computational analysis. Furthermore, participants self-reported a deeper interest in scientific research and careers as well as psychosocial outcomes (e.g., sense of belonging and self-efficacy) commonly associated with student success and persistence in STEM. Overall, this course provided a powerful combination of field, wet lab, and computational biology experiences for students, while also providing data potentially useful in pathogen surveillance, epidemiological tracking, and for the further study of environmental reservoirs of S. enterica.

Importance of anthropogenic sources at shaping the antimicrobial resistance profile of a peri-urban mesocarnivore

Anthropogenically derived antimicrobial-resistant bacteria (ARB) and antimicrobial resistance genes (ARG) have been detected in wildlife. The likelihood of detecting ARB and ARG in wildlife increases with wildlife exposure to anthropogenic sources of antimicrobial resistance (AMR). Whether anthropogenic sources also increase the risk for AMR to spread in bacteria of wildlife is not well understood. The spread of AMR in bacteria of wildlife can be estimated by examining the richness of ARB and ARG, and the prevalence of ARB that have mobilizable ARG (i.e., ARG that can be transferred across bacteria via plasmids). Here, we investigated whether raccoons (Procyon lotor), with different exposures to anthropogenic sources, differed in prevalence and richness of extended-spectrum cephalosporin-resistant (ESC-R) Escherichia coli, richness of ARG present in ESC-R E. coli, and prevalence of ESC-R E. coli with plasmid-associated ARG. Sampling took place over the course of 10 months at seven sites in Chicago, USA. ESC-R E. coli were isolated from over half of the 211 raccoons sampled and were more likely to be isolated from urban than suburban raccoons. When examining the whole-genome sequences of ESC-R E. coli, 56 sequence types were identified, most of which were associated with the ARG bla_CMY and bla_CTX-M. A greater richness of ESC-R E. coli sequence types was found at sites with a wastewater treatment plant (WWTP) than without, but no difference was detected based on urban context. ARG richness in ESC-R E. coli did not significantly vary by urban context nor with presence of a WWTP. Importantly, ESC-R E. coli carrying plasmid-associated bla_CTX-M and bla_CMY ARG were more likely to be isolated from raccoons sampled at sites with a WWTP than without. Our findings indicate that anthropogenic sources may shape the AMR profile of wildlife, reinforcing the need to prevent dissemination of AMR into the environment.

Genomic Characterization of Salmonella enterica Isolates From Retail Meat in Beijing, China

Salmonella enterica remains one of the leading causes of foodborne bacterial disease. Retail meat is a major source of human salmonellosis. However, comparative genomic analyses of S. enterica isolates from retail meat from different sources in China are lacking. A total of 341 S. enterica strains were isolated from retail meat in sixteen districts of Beijing, China, at three different time points (January 1st, May 1st, and October 1st) in 2017. Comparative genomics was performed to investigate the genetic diversity, virulence and antimicrobial resistance gene (ARG) profiles of these isolates. The most common serotype was S. Enteritidis (203/341, 59.5%), which dominated among isolates from three different time points during the year. Laboratory retesting confirmed the accuracy of the serotyping results predicted by the Salmonella In Silico Typing Resource (SISTR) (96.5%). The pangenome of the 341 S. enterica isolates contained 13,931 genes, and the core genome contained 3,635 genes. Higher Salmonella phage 118970 sal3 (219/341, 64.2%) and Gifsy-2 (206/341, 60.4%) prevalence contributed to the diversity of the accessory genes, especially those with unknown functions. IncFII(S), IncX1, and IncFIB(S) plasmid replicons were more common in these isolates and were major sources of horizontally acquired foreign genes. The virulence gene profile showed fewer virulence genes associated with type III secretion systems in certain isolates from chicken. A total of 88 different ARGs were found in the 341 isolates. Three beta-lactamases, namely, bla_CTX–M–55 (n = 15), bla_CTX–M–14 (n = 11), and bla_CTX–M–65 (n = 11), were more prevalent in retail meats. The emergence of qnrE1 and bla_{CTX–M–123} indicated a potential increase in the prevalence of retail meats. After the prohibition of colistin in China, three and four isolates were positive for the colistin resistance genes mcr-1.1 and mcr-9, respectively. Thus, we explored the evolution and genomic features of S. enterica isolates from retail meats in Beijing, China. The diverse ARGs of these isolates compromise food security and are a clinical threat.

Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae

Escherichia coli and other Enterobacteriaceae are diverse species with "open" pangenomes, where genes move intra- and interspecies via horizontal gene transfer. However, most analyses focus on clinical isolates. The pangenome dynamics of natural populations remain understudied, despite their suggested role as reservoirs for antimicrobial resistance (AMR) genes. Here, we analyze near-complete genomes for 827 Enterobacteriaceae (553 Escherichia and 274 non-Escherichia spp.) with 2292 circularized plasmids in total, collected from 19 locations (livestock farms and wastewater treatment works in the United Kingdom) within a 30-km radius at three time points over a year. We find different dynamics for chromosomal and plasmid-borne genes. Plasmids have a higher burden of AMR genes and insertion sequences, and AMR-gene-carrying plasmids show evidence of being under stronger selective pressure. Environmental niche and local geography both play a role in shaping plasmid dynamics. Our results highlight the importance of local strategies for controlling the spread of AMR.

Resistance determinants and their genetic context in enterobacteria from a longitudinal study of pigs reared under various husbandry conditions

Pigs are major reservoirs of resistant Enterobacteriaceae that can reach humans through consumption of contaminated meat or vegetables grown in manure-fertilized soil. Samples were collected from sows during lactation and their piglets at five time points spanning the production cycle. Cefotaxime-resistant bacteria were quantified and isolated from feed, feces, manures and carcasses of pigs reared with penicillin-using or antibiotic-free husbandries. The isolates were characterized by antibiotic susceptibility testing, whole genome sequencing and conjugation assays. The extended spectrum β-lactamase (ESBL) phenotype was more frequent in isolates originating from antibiotic-free animals, while the bacteria isolated from penicillin-using animals were on average resistant to a greater number of antibiotics. The ESBL-encoding genes identified were bla _CTX-M-1, bla _CTX-M-15 and bla _CMY-2 and they co-localised on plasmids with various genes encoding resistance to ß-lactams, co-trimoxazole, phenicols and tetracycline, all antibiotics used in pig production. Groups of genes conferring the observed resistance and the mobile elements disseminating multidrug resistance were determined. The observed resistance to ß-lactams was mainly due to the complementary actions of penicillin-binding proteins, an efflux pump and ß-lactamases. Most resistance determinants were shared by animals raised with or without antimicrobials. This suggests a key contribution of indigenous enterobacteria maternally transmitted along the sow lineage, regardless of antimicrobial use. It is unclear if the antimicrobial resistance observed in the enterobacteria populations of the commercial pig herds studied were present before the use of antibiotics, or the extent to which historical antimicrobial use exerted a selective pressure defining the resistant bacterial populations in farms using penicillin prophylaxis.Importance: Antimicrobial resistance is a global threat that needs to be fought on numerous fronts along the One Health continuum. Vast quantities of antimicrobials are used in agriculture to ensure animal welfare and productivity, and are arguably a driving force for the persistence of environmental and food-borne resistant bacteria. This study evaluated the impact of conventional, organic and other antibiotic-free husbandry practices on the frequency and nature of antimicrobial resistance genes and multidrug resistant enterobacteria. It provides knowledge about the relative contribution of specific resistance determinants to observed antibiotic resistance. It also showed the clear co-selection of genes coding for extended-spectrum beta-lactamases and genes coding for the resistance to antibiotics commonly used for prophylaxis or in curative treatments in pig operations.

The extent of carbapenemase-encoding genes in public genome sequences

Genome sequences provide information on the genetic elements present in an organism, and currently there are databases containing hundreds of thousands of bacterial genome sequences. These repositories allow for mining patterns concerning antibiotic resistance gene occurrence in both pathogenic and non-pathogenic bacteria in e.g. natural or animal environments, and link these to relevant metadata such as bacterial host species, country and year of isolation, and co-occurrence with other resistance genes. In addition, the advances in the prediction of mobile genetic elements, and discerning chromosomal from plasmid DNA, broadens our view on the mechanism mediating dissemination. In this study we utilize the vast amount of data in the public database PATRIC to investigate the dissemination of carbapenemase-encoding genes (CEGs), the emergence and spread of which is considered a grave public health concern. Based on publicly available genome sequences from PATRIC and manually curated CEG sequences from the beta lactam database, we found 7,964 bacterial genomes, belonging to at least 70 distinct species, that carry in total 9,892 CEGs, amongst which bla _NDM, bla _OXA, bla _VIM, bla _IMP and bla _KPC. We were able to distinguish between chromosomally located resistance genes (4,137; 42%) and plasmid-located resistance genes (5,753; 58%). We found that a large proportion of the identified CEGs were identical, i.e. displayed 100% nucleotide similarity in multiple bacterial species (8,361 out of 9,892 genes; 85%). For example, the New Delhi metallo-beta-lactamase NDM-1 was found in 42 distinct bacterial species, and present in seven different environments. Our data show the extent of carbapenem-resistance far beyond the canonical species Acetinobacter baumannii, Klebsiella pneumoniae or Pseudomonas aeruginosa. These types of data complement previous systematic reviews, in which carbapenem-resistant Enterobacteriaceae were found in wildlife, livestock and companion animals. Considering the widespread distribution of CEGs, we see a need for comprehensive surveillance and transmission studies covering more host species and environments, akin to previous extensive surveys that focused on extended spectrum beta-lactamases. This may help to fully appreciate the spread of CEGs and improve the understanding of mechanisms underlying transmission, which could lead to interventions minimizing transmission to humans.

Outcome of Different Sequencing and Assembly Approaches on the Detection of Plasmids and Localization of Antimicrobial Resistance Genes in Commensal Escherichia coli

Antimicrobial resistance (AMR) is a major threat to public health worldwide. Currently, AMR typing changes from phenotypic testing to whole-genome sequence (WGS)-based detection of resistance determinants for a better understanding of the isolate diversity and elements involved in gene transmission (e.g., plasmids, bacteriophages, transposons). However, the use of WGS data in monitoring purposes requires suitable techniques, standardized parameters and approved guidelines for reliable AMR gene detection and prediction of their association with mobile genetic elements (plasmids). In this study, different sequencing and assembly strategies were tested for their suitability in AMR monitoring in Escherichia coli in the routines of the German National Reference Laboratory for Antimicrobial Resistances. To assess the outcomes of the different approaches, results from in silico predictions were compared with conventional phenotypic- and genotypic-typing data. With the focus on (fluoro)quinolone-resistant E.coli, five qnrS-positive isolates with multiple extrachromosomal elements were subjected to WGS with NextSeq (Illumina), PacBio (Pacific BioSciences) and ONT (Oxford Nanopore) for in depth characterization of the qnrS1-carrying plasmids. Raw reads from short- and long-read sequencing were assembled individually by Unicycler or Flye or a combination of both (hybrid assembly). The generated contigs were subjected to bioinformatics analysis. Based on the generated data, assembly of long-read sequences are error prone and can yield in a loss of small plasmid genomes. In contrast, short-read sequencing was shown to be insufficient for the prediction of a linkage of AMR genes (e.g., qnrS1) to specific plasmid sequences. Furthermore, short-read sequencing failed to detect certain duplications and was unsuitable for genome finishing. Overall, the hybrid assembly led to the most comprehensive typing results, especially in predicting associations of AMR genes and mobile genetic elements. Thus, the use of different sequencing technologies and hybrid assemblies currently represents the best approach for reliable AMR typing and risk assessment.

A high prevalence of blaOXA-48 in Klebsiella (Raoultella) ornithinolytica and related species in hospital wastewater in South West England

Klebsiella species occupy a wide range of environmental and animal niches, and occasionally cause opportunistic infections that are resistant to multiple antibiotics. In particular, Klebsiella pneumoniae (Kpne) has gained notoriety as a major nosocomial pathogen, due principally to the rise in non-susceptibility to carbapenems and other beta-lactam antibiotics. Whilst it has been proposed that the urban water cycle facilitates transmission of pathogens between clinical settings and the environment, the level of risk posed by resistant Klebsiella strains in hospital wastewater remains unclear. We used whole genome sequencing (WGS) to compare Klebsiella species in contemporaneous samples of wastewater from an English hospital and influent to the associated wastewater treatment plant (WWTP). As we aimed to characterize representative samples of Klebsiella communities, we did not actively select for antibiotic resistance (other than for ampicillin), nor for specific Klebsiella species. Two species, Kpne and K. (Raoultella) ornithinolytica (Korn), were of equal dominance in the hospital wastewater, and four other Klebsiella species were present in low abundance in this sample. In contrast, despite being the species most closely associated with healthcare settings, Kpne was the dominant species within the WWTP influent. In total, 29 % of all isolates harboured the blaOXA-48 gene on a pOXA-48-like plasmid, and these isolates were almost exclusively recovered from the hospital wastewater. This gene was far more common in Korn (68 % of isolates) than in Kpne (3.4 % of isolates). In general plasmid-borne, but not chromosomal, resistance genes were significantly enriched in the hospital wastewater sample. These data implicate hospital wastewater as an important reservoir for antibiotic-resistant Klebsiella, and point to an unsuspected role of species within the Raoultella group in the maintenance and dissemination of plasmid-borne blaOXA-48. This article contains data hosted by Microreact.

Genomic Markers for Quaternary Ammonium Compound Resistance as a Persistence Indicator for Listeria monocytogenes Contamination in Food Manufacturing Environments

ABSTRACT

Persistent contamination of food manufacturing environments by Listeria monocytogenes is an important public health risk, because such contamination events defy standard sanitization protocols, for example, the application of quaternary ammonium compounds such as benzalkonium chloride (BC), providing a source for prolonged dissemination of the bacteria in food products. We performed whole genome sequencing analyses of 1,279 well-characterized L. monocytogenes isolates from various foods and food manufacturing environments and identified the bcrABC gene cassette associated with BC resistance in 531 (41.5%) isolates. The bcrABC cassette was significantly associated with L. monocytogenes isolates belonging to clonal complex (CC) 321, CC155, CC204, and CC199, which are among the 10 most prevalent genotypes recovered from foods and food production environments. All but 1 of the 177 CC321 isolates harbored the bcrABC cassette. In addition, 384 (38.6%) of the 994 isolates recovered from foods representing 67 different CCs and 119 (59.2%) of isolates from food manufacturing environmental samples representing 26 different CCs were found to harbor the intact bcrABC cassette. A representative set of 69 isolates with and without bcrABC was assayed for the ability to grow in the presence of BC, and 34 of 35 isolates harboring the bcrABC cassette exhibited MICs of ≥10 μg/mL BC. Determination of bcrABC in isolates could be achieved using both PCR and whole genome sequencing techniques, providing food testing laboratories with options for the characterization of isolates. The ability to determine markers of quaternary ammonium compound resistance such as bcrABC and epidemiologic lineage may provide risk managers with a tool to assess the potential for persistent contamination of the food manufacturing environment and the need for more targeted surveillance to ensure the efficacy of mitigation actions.

HIGHLIGHTS

Whole Genome Sequencing Applied to Pathogen Source Tracking in Food Industry: Key Considerations for Robust Bioinformatics Data Analysis and Reliable Results Interpretation

Whole genome sequencing (WGS) has arisen as a powerful tool to perform pathogen source tracking in the food industry thanks to several developments in recent years. However, the cost associated to this technology and the degree of expertise required to accurately process and understand the data has limited its adoption at a wider scale. Additionally, the time needed to obtain actionable information is often seen as an impairment for the application and use of the information generated via WGS. Ongoing work towards standardization of wet lab including sequencing protocols, following guidelines from the regulatory authorities and international standardization efforts make the technology more and more accessible. However, data analysis and results interpretation guidelines are still subject to initiatives coming from distinct groups and institutions. There are multiple bioinformatics software and pipelines developed to handle such information. Nevertheless, little consensus exists on a standard way to process the data and interpret the results. Here, we want to present the constraints we face in an industrial setting and the steps we consider necessary to obtain high quality data, reproducible results and a robust interpretation of the obtained information. All of this, in a time frame allowing for data-driven actions supporting factories and their needs.

Sequence characterisation and novel insights into bovine mastitis-associated Streptococcus uberis in dairy herds

Streptococcus uberis is one of the most frequent mastitis-causing pathogens isolated from dairy cows. Further understanding of S. uberis genetics may help elucidate the disease pathogenesis. We compared the genomes of S. uberis isolates cultured from dairy cows located in distinctly different geographic regions of Australia. All isolates had novel multi locus sequence types (MLST) indicating a highly diverse population of S. uberis. Global clonal complexes (GCC) were more conserved. GCC ST86 and GCC ST143 represented 30% of the total isolates (n = 27) and were clustered within different geographic regions. Core genome phylogeny revealed low phylogenetic clustering by region, isolation source, and MLST. Identification of putative sortase (srtA) substrates and generation of a custom putative virulence factor database revealed genes which may explain the affinity of S. uberis for mammary tissue, evasion of antimicrobial efforts and disease pathogenesis. Of 27 isolates, four contained antibiotic resistance genes including an antimicrobial resistance cluster containing mel/mef(A), mrsE, vatD, lnuD, and transposon-mediated lnuC was also identified. These are novel genes for S. uberis, which suggests interspecies lateral gene transfer. The presence of resistance genes across the two geographic regions tested within one country supports the need for a careful, tailored, implementation and monitoring of antimicrobial stewardship.

Complete Genome Sequences of 12 Quinolone-Resistant Escherichia coli Strains Containing qnrS1 Based on Hybrid Assemblies

In total, 12 quinolone-resistant Escherichia coli (QREC) strains containing qnrS1 were submitted to long-read sequencing using a FLO-MIN106 flow cell on a MinION device. The long reads were assembled with short reads (Illumina) and analyzed using the MOB-suite pipeline. Six of these QREC genome sequences were closed after hybrid assembly.

In total, 12 quinolone-resistant Escherichia coli (QREC) strains containing qnrS1 were submitted to long-read sequencing using a FLO-MIN106 flow cell on a MinION device. The long reads were assembled with short reads (Illumina) and analyzed using the MOB-suite pipeline. Six of these QREC genome sequences were closed after hybrid assembly.

Collateral sensitivity associated with antibiotic resistance plasmids

Collateral sensitivity (CS) is a promising alternative approach to counteract the rising problem of antibiotic resistance (ABR). CS occurs when the acquisition of resistance to one antibiotic produces increased susceptibility to a second antibiotic. Recent studies have focused on CS strategies designed against ABR mediated by chromosomal mutations. However, one of the main drivers of ABR in clinically relevant bacteria is the horizontal transfer of ABR genes mediated by plasmids. Here, we report the first analysis of CS associated with the acquisition of complete ABR plasmids, including the clinically important carbapenem-resistance conjugative plasmid pOXA-48. In addition, we describe the conservation of CS in clinical E. coli isolates and its application to selectively kill plasmid-carrying bacteria. Our results provide new insights that establish the basis for developing CS-informed treatment strategies to combat plasmid-mediated ABR.

Distribution of Beta-Lactamase Producing Gram-Negative Bacterial Isolates in Isabela River of Santo Domingo, Dominican Republic

Bacteria carrying antibiotic resistance genes (ARGs) are naturally prevalent in lotic ecosystems such as rivers. Their ability to spread in anthropogenic waters could lead to the emergence of multidrug-resistant bacteria of clinical importance. For this study, three regions of the Isabela river, an important urban river in the city of Santo Domingo, were evaluated for the presence of ARGs. The Isabela river is surrounded by communities that do not have access to proper sewage systems; furthermore, water from this river is consumed daily for many activities, including recreation and sanitation. To assess the state of antibiotic resistance dissemination in the Isabela river, nine samples were collected from these three bluedistinct sites in June 2019 and isolates obtained from these sites were selected based on resistance to beta-lactams. Physico-chemical and microbiological parameters were in accordance with the Dominican legislation. Matrix-assisted laser desorption ionization-time of flight mass spectrometry analyses of ribosomal protein composition revealed a total of 8 different genera. Most common genera were as follows: Acinetobacter (44.6%) and Escherichia (18%). Twenty clinically important bacterial isolates were identified from urban regions of the river; these belonged to genera Escherichia (n = 9), Acinetobacter (n = 8), Enterobacter (n = 2), and Klebsiella (n = 1). Clinically important multi-resistant isolates were not obtained from rural areas. Fifteen isolates were selected for genome sequencing and analysis. Most isolates were resistant to at least three different families of antibiotics. Among beta-lactamase genes encountered, we found the presence of blaTEM, blaOXA, blaSHV, and blaKPC through both deep sequencing and PCR amplification. Bacteria found from genus Klebsiella and Enterobacter demonstrated ample repertoire of antibiotic resistance genes, including resistance from a family of last resort antibiotics reserved for dire infections: carbapenems. Some of the alleles found were KPC-3, OXA-1, OXA-72, OXA-132, CTX-M-55, CTX-M-15, and TEM-1.

National Prevalence of Salmonella enterica Serotype Kentucky ST198 with High-Level Resistance to Ciprofloxacin and Extended-Spectrum Cephalosporins in China, 2013 to 2017

Salmonella enterica serotype Kentucky is frequently associated with high-level fluoroquinolone resistance and has gained epidemiological importance globally. A retrospective screening was performed to understand the national prevalence of ciprofloxacin-resistant S Kentucky in China. S. enterica strains (n = 15,405) were collected within the frame of two national surveillance networks between 2013 and 2017. Thirty-three S. Kentucky strains were detected in 5 of 10 provinces, and 27 were assigned to sequence type 198 (ST198). The 27 isolates were multidrug resistant, with high-level resistance to ciprofloxacin, and 21 isolates were further resistant to extended-spectrum cephalosporins (ESCs). Phylogenomic analysis classified ST198 isolates into two clades (198.1 and 198.2), and recent occurrences of inter-/intraregion and interhost transmission were identified. Phylogenetic reconstruction with a global collection showed that one subclade of clade 198.2 was clustered with historical strains from Egypt, and the other one was clustered with strains from Southeast Asia. Isolates of clade 198.1 were clustered with strains isolated from North America. The various patterns of mutations detected in quinolone resistance-determining regions of GyrA and ParC are accordant with the phylogenetic structure. These findings indicate that our isolates may have various origins. SGI1 was exclusively detected in isolates of clade 198.2 with a highly mosaic structure, which were mainly identified as SGI1-K derivatives. Plasmid-mediated quinolone resistance genes qnrS1 and aac(6')-Ib-cr were identified in three isolates, and bla CTX-M-9 and bla CTX-M-27 were detected in 20 of 21 ESC-resistant isolates. This is the first report of the genetic and epidemiological characterization for the S Kentucky epidemic clone ST198 in China, warranting the necessity of surveillance for the high-risk clone.IMPORTANCE Ciprofloxacin and extended-spectrum cephalosporins are the choice for treatment of severe nontyphoidal S. enterica infections in adults. S. enterica serotype Kentucky ST198 has gained epidemiological importance globally, because the clone is frequently resistant to both of these high-level-resistance drug groups. The genetic and epidemiological characterization of S. Kentucky has been well studied in Western countries; however, the information is unclear for China. To fill in the gap, we here did a retrospective screening on a large collection in China, and ST198 isolates were systematically analyzed by whole-genome sequencing. Our study revealed that multidrug-resistant ST198 has spread in five provinces, and the occurrences of interregion and cross-host clonal disseminations were detected. Of note, phylogenomic analysis suggests that the Chinese isolates may have emerged with diverse origins, including Egypt, Southeast Asia, and North America. This study warrants the necessity of surveillance for the high-risk clone to prevent its further dissemination in China.

mMGE: a database for human metagenomic extrachromosomal mobile genetic elements

Extrachromosomal mobile genetic elements (eMGEs), including phages and plasmids, that can move across different microbes, play important roles in genome evolution and shaping the structure of microbial communities. However, we still know very little about eMGEs, especially their abundances, distributions and putative functions in microbiomes. Thus, a comprehensive description of eMGEs is of great utility. Here we present mMGE, a comprehensive catalog of 517 251 non-redundant eMGEs, including 92 492 plasmids and 424 759 phages, derived from diverse body sites of 66 425 human metagenomic samples. About half the eMGEs could be further grouped into 70 074 clusters using relaxed criteria (referred as to eMGE clusters below). We provide extensive annotations of the identified eMGEs including sequence characteristics, taxonomy affiliation, gene contents and their prokaryotic hosts. We also calculate the prevalence, both within and across samples for each eMGE and eMGE cluster, enabling users to see putative associations of eMGEs with human phenotypes or their distribution preferences. All eMGE records can be browsed or queried in multiple ways, such as eMGE clusters, metagenomic samples and associated hosts. The mMGE is equipped with a user-friendly interface and a BLAST server, facilitating easy access/queries to all its contents easily. mMGE is freely available for academic use at: https://mgedb.comp-sysbio.org.

Complete Genome Sequences for 36 Canadian Salmonella enterica Serovar Typhimurium and I 1,4,[5],12:i:- Isolates from Clinical and Animal Sources

Here, we report the complete genome sequences for 36 Canadian isolates of Salmonella enterica subsp. enterica serovar Typhimurium and its monophasic variant I 1,4,[5]:12:i:– from both clinical and animal sources. These genome sequences will provide useful references for understanding the genetic variation within this prominent serotype.

Here, we report the complete genome sequences for 36 Canadian isolates of Salmonella enterica subsp. enterica serovar Typhimurium and its monophasic variant I 1,4,[5]:12:i:– from both clinical and animal sources. These genome sequences will provide useful references for understanding the genetic variation within this prominent serotype.

Detection of carbapenem-resistance genes in bacteria isolated from wastewater in Ontario

Bacterial carbapenem resistance is a major public health concern since these antimicrobials are often the last resort to treat serious human infections. To evaluate methodologies for detection of carbapenem resistance, carbapenem-tolerant bacteria were isolated from wastewater treatment plants in Toronto, Ottawa, and Arnprior, Ontario. A total of 135 carbapenem-tolerant bacteria were recovered. Polymerase chain reaction (PCR) indicated the presence of carbapenem hydrolysing enzymes KPC ( n = 10), GES ( n = 5), VIM ( n = 7), and IMP ( n = 1), and β-lactamases TEM ( n = 7), PER ( n = 1), and OXA-variants ( n = 16). A subset of 46 isolates were sequenced and analysed using ResFinder and CARD-RGI. Both programs detected carbapenem resistance genes in 35 sequenced isolates and antimicrobial resistance genes (ARGs) conferring resistance to multiple class of other antibiotics. Where β-lactamase resistance genes were not initially identified, lowering the thresholds for ARG detection enabled identification of closely related β-lactamases. However, no known carbapenem resistance genes were found in seven sequenced Pseudomonas spp. isolates. Also of note was a multi-drug-resistant Klebsiella pneumoniae isolate from Ottawa, which harboured resistance to seven antimicrobial classes including β-lactams. These results highlight the diversity of genes encoding carbapenem resistance in Ontario and the utility of whole genome sequencing over PCR for ARG detection where resistance may result from an assortment of genes.

Listeria valentina sp. nov., isolated from a water trough and the faeces of healthy sheep

In the context of a study on the occurrence of Listeria species in an animal farm environment in Valencia, Spain, six Listeria-like isolates could not be assigned to any known species. Phylogenetic analysis based on the 16S rRNA gene and on 231 Listeria core genes grouped these isolates in a monophyletic clade within the genus Listeria, with highest similarity to Listeria thailandensis. Whole-genome sequence analyses based on in silico DNA-DNA hybridization, the average nucleotide blast and the pairwise amino acid identities against all currently known Listeria species confirmed that these isolates constituted a new taxon within the genus Listeria. Phenotypically, these isolates differed from other Listeria species mainly by the production of acid from inositol, the absence of acidification in presence of methyl α-d-glucoside, and the absence of α-mannosidase and nitrate reductase activities. The name Listeria valentina sp. nov. is proposed for this novel species, and the type strain is CLIP 2019/00642^T (=CIP 111799^T=DSM 110544^T).

gplas: a comprehensive tool for plasmid analysis using short-read graphs

Summary

Plasmids can horizontally transmit genetic traits, enabling rapid bacterial adaptation to new environments and hosts. Short-read whole-genome sequencing data are often applied to large-scale bacterial comparative genomics projects but the reconstruction of plasmids from these data is facing severe limitations, such as the inability to distinguish plasmids from each other in a bacterial genome. We developed gplas, a new approach to reliably separate plasmid contigs into discrete components using sequence composition, coverage, assembly graph information and network partitioning based on a pruned network of plasmid unitigs. Gplas facilitates the analysis of large numbers of bacterial isolates and allows a detailed analysis of plasmid epidemiology based solely on short-read sequence data.

Availability and implementation

Gplas is written in R, Bash and uses a Snakemake pipeline as a workflow management system. Gplas is available under the GNU General Public License v3.0 at https://gitlab.com/sirarredondo/gplas.git.

Supplementary information

Supplementary data are available at Bioinformatics online.

Survey and Sequence Characterization of Bovine Mastitis-Associated Escherichia coli in Dairy Herds

Escherichia coli is frequently associated with mastitis in cattle. "Pathogenic" and "commensal" isolates appear to be genetically similar. With a few exceptions, no notable genotypic differences have been found between commensal and mastitis-associated E. coli. In this study, 24 E. coli strains were isolated from dairy cows with clinical mastitis in three geographic regions of Australia (North Queensland, South Queensland, and Victoria), sequenced, then genomically surveyed. There was no observed relationship between sequence type (ST) and region (p = 0.51). The most common Multi Locus Sequence Type was ST10 (38%), then ST4429 (13%). Pangenomic analysis revealed a soft-core genome of 3,463 genes, including genes associated with antibiotic resistance, chemotaxis, motility, adhesion, biofilm formation, and pili. A total of 36 different plasmids were identified and generally found to have local distributions (p = 0.02). Only 2 plasmids contained antibiotic resistance genes, a p1303_5-like plasmid encoding multidrug-resistance (trimethoprim, quaternary ammonium, beta-lactam, streptomycin, sulfonamide, and kanamycin) from two North Queensland isolates on the same farm, while three Victorian isolates from the same farm contained a pCFSAN004177P_01-like plasmid encoding tetracycline-resistance. This pattern is consistent with a local spread of antibiotic resistance through plasmids of bovine mastitis cases. Notably, co-occurrence of plasmids containing virulence factors/antibiotic resistance with putative mobilization was rare, though the multidrug resistant p1303_5-like plasmid was predicted to be conjugative and is of some concern. This survey has provided greater understanding of antibiotic resistance within E. coli-associated bovine mastitis which will allow greater prediction and improved decision making in disease management.

The role of mobile genetic elements in organic micropollutant degradation during biological wastewater treatment

Wastewater treatment plants (WWTPs) are crucial for producing clean effluents from polluting sources such as hospitals, industries, and municipalities. In recent decades, many new organic compounds have ended up in surface waters in concentrations that, while very low, cause (chronic) toxicity to countless organisms. These organic micropollutants (OMPs) are usually quite recalcitrant and not sufficiently removed during wastewater treatment. Microbial degradation plays a pivotal role in OMP conversion. Microorganisms can adapt their metabolism to the use of novel molecules via mutations and rearrangements of existing genes in new clusters. Many catabolic genes have been found adjacent to mobile genetic elements (MGEs), which provide a stable scaffold to host new catabolic pathways and spread these genes in the microbial community. These mobile systems could be engineered to enhance OMP degradation in WWTPs, and this review aims to summarize and better understand the role that MGEs might play in the degradation and wastewater treatment process. Available data about the presence of catabolic MGEs in WWTPs are reviewed, and current methods used to identify and measure MGEs in environmental samples are critically evaluated. Finally, examples of how these MGEs could be used to improve micropollutant degradation in WWTPs are outlined. In the near future, advances in the use of MGEs will hopefully enable us to apply selective augmentation strategies to improve OMP conversion in WWTPs.

Extensive antimicrobial resistance mobilization via multicopy plasmid encapsidation mediated by temperate phages

OBJECTIVES

To investigate the relevance of multicopy plasmids in antimicrobial resistance and assess their mobilization mediated by phage particles.

METHODS

Several databases with complete sequences of plasmids and annotated genes were analysed. The 16S methyltransferase gene armA conferring high-level aminoglycoside resistance was used as a marker in eight different plasmids, from different incompatibility groups, and with differing sizes and plasmid copy numbers. All plasmids were transformed into Escherichia coli bearing one of four different lysogenic phages. Upon induction, encapsidation of armA in phage particles was evaluated using qRT-PCR and Southern blotting.

RESULTS

Multicopy plasmids carry a vast set of emerging clinically important antimicrobial resistance genes. However, 60% of these plasmids do not bear mobility (MOB) genes. When carried on these multicopy plasmids, mobilization of a marker gene armA into phage capsids was up to 10000 times more frequent than when it was encoded by a large plasmid with a low copy number.

CONCLUSIONS

Multicopy plasmids and phages, two major mobile genetic elements (MGE) in bacteria, represent a novel high-efficiency transmission route of antimicrobial resistance genes that deserves further investigation.

Identification of the conjugative and mobilizable plasmid fragments in the plasmidome using sequence signatures

Plasmids are the key element in horizontal gene transfer in the microbial community. Recently, a large number of experimental and computational methods have been developed to obtain the plasmidomes of microbial communities. Distinguishing transmissible plasmid sequences, which are derived from conjugative or at least mobilizable plasmids, from non-transmissible plasmid sequences in the plasmidome is essential for understanding the diversity of plasmids and how they regulate the microbial community. Unfortunately, due to the highly fragmented characteristics of DNA sequences in the plasmidome, effective identification methods are lacking. In this work, we used information entropy from information theory to assess the randomness of synonymous codon usage over 4424 plasmid genomes. The results showed that for all amino acids, the choice of a synonymous codon in conjugative and mobilizable plasmids is more random than that in non-transmissible plasmids, indicating that transmissible plasmids have different sequence signatures from non-transmissible plasmids. Inspired by this phenomenon, we further developed a novel algorithm named PlasTrans. PlasTrans takes the triplet code sequences and base sequences of plasmid DNA fragments as input and uses the convolutional neural network of the deep learning technique to further extract the more complex signatures of the plasmid sequences and identify the conjugative and mobilizable DNA fragments. Tests showed that PlasTrans could achieve an AUC of as high as 84-91%, even though the fragments only contained hundreds of base pairs. To the best of our knowledge, this is the first quantitative analysis of the difference in sequence signatures between transmissible and non-transmissible plasmids, and we developed the first tool to perform transferability annotation for DNA fragments in the plasmidome. We expect that PlasTrans will be a useful tool for researchers who analyse the properties of novel plasmids in the microbial community and horizontal gene transfer, especially the spread of resistance genes and virulence factors associated with plasmids. PlasTrans is freely available via https://github.com/zhenchengfang/PlasTrans.

Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores

Plasmids are extrachromosomal genetic elements that replicate independently of the chromosome and play a vital role in the environmental adaptation of bacteria. Due to potential mobilization or conjugation capabilities, plasmids are important genetic vehicles for antimicrobial resistance genes and virulence factors with huge and increasing clinical implications. They are therefore subject to large genomic studies within the scientific community worldwide. As a result of rapidly improving next-generation sequencing methods, the quantity of sequenced bacterial genomes is constantly increasing, in turn raising the need for specialized tools to (i) extract plasmid sequences from draft assemblies, (ii) derive their origin and distribution, and (iii) further investigate their genetic repertoire. Recently, several bioinformatic methods and tools have emerged to tackle this issue; however, a combination of high sensitivity and specificity in plasmid sequence identification is rarely achieved in a taxon-independent manner. In addition, many software tools are not appropriate for large high-throughput analyses or cannot be included in existing software pipelines due to their technical design or software implementation. In this study, we investigated differences in the replicon distributions of protein-coding genes on a large scale as a new approach to distinguish plasmid-borne from chromosome-borne contigs. We defined and computed statistical discrimination thresholds for a new metric: the replicon distribution score (RDS), which achieved an accuracy of 96.6 %. The final performance was further improved by the combination of the RDS metric with heuristics exploiting several plasmid-specific higher-level contig characterizations. We implemented this workflow in a new high-throughput taxon-independent bioinformatics software tool called Platon for the recruitment and characterization of plasmid-borne contigs from short-read draft assemblies. Compared to PlasFlow, Platon achieved a higher accuracy (97.5 %) and more balanced predictions (F1=82.6 %) tested on a broad range of bacterial taxa and better or equal performance against the targeted tools PlasmidFinder and PlaScope on sequenced Escherichia coli isolates. Platon is available at: http://platon.computational.bio/.

Prediction of Salmonella serovars isolated from clinical and food matrices in Lebanon and genomic-based investigation focusing on Enteritidis serovar

Salmonella enterica subsp. enterica serovars are considered major causes of food poisoning and we performed this study because Salmonella is a burden in Lebanon. The present study investigated the ability of genomic information to predict serovar using a collection of Salmonella isolates from infected humans (n = 24) and contaminated food (n = 63) in Lebanon. Further, the phylogenomic relationships of the serovar the predominated in Lebanon (i.e., S. Enteritidis; n = 25) were investigated in comparison with isolates from other countries (n = 130) based on coregenome single nucleotide polymorphisms (SNPs). Genetic elements, specifically Salmonella pathogenicity islands (SPIs), plasmid replicons, and antibiotic-resistance genes were screened in S. Enteritidis genomes (n = 155). Our results revealed that the Salmonella serovars identification by seroagglutination from the samples isolated in Lebanon (n = 87) was highly correlated with the genomic-based prediction of serovars (80.4-85.0% with SeqSero1 and 93.1-94.2% with SeqSero2). The Salmonella serovars isolated from human and food samples in Lebanon were mainly Enteritidis (28.7%) and Infantis (26%). To a rare extent, other serovars included Amager, Anatum, Bredeney, Chincol, Heidelberg, Hofit, Kentucky, Montevideo, Muenster, Newport, Schwarzengrund, Senftenberg and Typhimurium. In comparison with other countries, S. Enteritidis samples isolated in Lebanon (56 ± 27 intra-group pairwise SNP differences) presented a strong phylogenomic relativeness at the coregenome level with samples, as for example with samples isolated from Syria (65 ± 31 inter-group pairwise SNP differences). Most of the studied S. Enteritidis genomes encoded 10 SPIs involved in survival in immune cells (i.e. SPIs 1, 2, 3, 4, 5, 12, 13, 14, 16 and 17). The plasmid replicons IncFIB (S)_1 and IncFII (S)_1 encoding elements involved in virulence were identified in the majority of the S. Enteritidis genomes (94% and 96%, respectively), the majority exhibiting aminoglycosides (gene aac(6')-Iaa_1). The IncI_1_Alpha replicon responsible for ampicillin-resistance was only detected in 2 of 25 S. Enteritidis Lebanese strains. Genomic-based risk assessment of Salmonella serovars in Lebanon showed that food imported from Syria might be an origin of the S. Enteritidis human cases in Lebanon. The detection of several SPIs involved in the survival, plasmid replicons involved in virulence, and aminoglycoside-resistance genes, emphasizes that S. Enteritidis is of paramount importance for public health in Lebanon and other countries.

Mobility of β-Lactam Resistance Under Bacterial Co-infection and Ampicillin Treatment in a Mouse Model

Ingestion of food- or waterborne antibiotic-resistant bacteria may lead to the dissemination of antibiotic-resistance genes in the gut microbiota and the development of antibiotic-resistant bacterial infection, a significant threat to animal and public health. Food or water may be contaminated with multiple resistant bacteria, but animal models on gene transfer were mainly based on single-strain infections. In this study, we investigated the mobility of β-lactam resistance following infection with single- versus multi-strain of resistant bacteria under ampicillin treatment. We characterized three bacterial strains isolated from food-animal production systems, Escherichia coli O80:H26 and Salmonella enterica serovars Bredeney and Heidelberg. Each strain carries at least one conjugative plasmid that encodes a β-lactamase. We orally infected mice with each or all three bacterial strain(s) in the presence or absence of ampicillin treatment. We assessed plasmid transfer from the three donor bacteria to an introduced E. coli CV601gfp recipient in the mouse gut, and evaluated the impacts of the bacterial infection on gut microbiota and gut health. In the absence of ampicillin treatment, none of the donor or recipient bacteria established in the normal gut microbiota and plasmid transfer was not detected. In contrast, the ampicillin treatment disrupted the gut microbiota and enabled S. Bredeney and Heidelberg to colonize and transfer their plasmids to the E. coli CV601gfp recipient. E. coli O80:H26 on its own failed to colonize the mouse gut. However, during co-infection with the two Salmonella strains, E. coli O80:H26 colonized and transferred its plasmid to the E. coli CV601gfp recipient and a residential E. coli O2:H6 strain. The co-infection significantly increased plasmid transfer frequency, enhanced Proteobacteria expansion and resulted in inflammation in the mouse gut. Our findings suggest that single-strain infection models for evaluating in vivo gene transfer may underrepresent the consequences of multi-strain infections following the consumption of heavily contaminated food or water.

Detection of extended-spectrum beta-lactamase (ESBL) genes and plasmid replicons in Enterobacteriaceae using PlasmidSPAdes assembly of short-read sequence data

Knowledge of the epidemiology of plasmids is essential for understanding the evolution and spread of antimicrobial resistance. PlasmidSPAdes attempts to reconstruct plasmids using short-read sequence data. Accurate detection of extended-spectrum beta-lactamase (ESBL) genes and plasmid replicon genes is a prerequisite for the use of plasmid assembly tools to investigate the role of plasmids in the spread and evolution of ESBL production in Enterobacteriaceae. This study evaluated the performance of PlasmidSPAdes plasmid assembly for Enterobacteriaceae in terms of detection of ESBL-encoding genes, plasmid replicons and chromosomal wgMLST genes, and assessed the effect of k-mer size. Short-read sequence data for 59 ESBL-producing Enterobacteriaceae were assembled with PlasmidSPAdes using different k-mer sizes (21, 33, 55, 77, 99 and 127). For every k-mer size, the presence of ESBL genes, plasmid replicons and a selection of chromosomal wgMLST genes in the plasmid assembly was determined. Out of 241 plasmid replicons and 66 ESBL genes detected by whole-genome assembly, 213 plasmid replicons [88 %; 95 % confidence interval (CI): 83.9-91.9] and 43 ESBL genes (65 %; 95 % CI: 53.1-75.6) were detected in the plasmid assemblies obtained by PlasmidSPAdes. For most ESBL genes (83.3 %) and plasmid replicons (72.0 %), detection results did not differ between the k-mer sizes used in the plasmid assembly. No optimal k-mer size could be defined for the number of ESBL genes and plasmid replicons detected. For most isolates, the number of chromosomal wgMLST genes detected in the plasmid assemblies decreased with increasing k-mer size. Based on our findings, PlasmidSPAdes is not a suitable plasmid assembly tool for short-read sequence data for ESBL-encoding plasmids of Enterobacteriaceae.

Genomic and Proteomic Characterization of the Extended-Spectrum β-Lactamase (ESBL)-Producing Escherichia coli Strain CCUG 73778: A Virulent, Nosocomial Outbreak Strain

Escherichia coli strain CCUG 78773 is a virulent extended-spectrum β-lactamase (ESBL)-producing ST131-O25b type strain isolated during an outbreak at a regional university hospital. The complete and closed genome sequence, comprising one chromosome (5,076,638 bp) and six plasmids (1718–161,372 bp), is presented. Characterization of the genomic features detected the presence of 59 potential antibiotic resistance factors, including three prevalent β-lactamases. Several virulence associated elements were determined, mainly related with adherence, invasion, biofilm formation and antiphagocytosis. Twenty-eight putative type II toxin-antitoxin systems were found. The plasmids were characterized, through in silico analyses, confirming the two β-lactamase-encoding plasmids to be conjugative, while the remaining plasmids were mobilizable. BLAST analysis of the plasmid sequences showed high similarity with plasmids in E. coli from around the world. Expression of many of the described virulence and AMR factors was confirmed by proteomic analyses, using bottom-up, liquid chromatography-tandem mass spectrometry (LC-MS/MS). The detailed characterization of E. coli strain CCUG 78773 provides a reference for the relevance of genetic elements, as well as the characterization of antibiotic resistance and the spread of bacteria harboring ESBL genes in the hospital environment.

Reconstruction of plasmids by shotgun sequencing from environmental DNA: which bioinformatic workflow?

Plasmids play important roles in microbial evolution and also in the spread of antibiotic resistance. Plasmid sequences are extensively studied from clinical isolates but rarely from the environment with a metagenomic approach focused on the plasmid fraction referred to as the plasmidome. A clear challenge in this context is to define a workflow for discriminating plasmids from chromosomal contaminants existing in the plasmidome. For this purpose, we benchmarked existing tools from assembly to detection of the plasmids by reference-free methods (cBar and PlasFlow) and database-guided approaches. Our simulations took into account short-reads alone or combined with moderate long-reads like those actually generated in environmental genomics experiments. This benchmark allowed us to select the best tools for limiting false-positives associated to plasmid prediction tools and a combination of reference-guided methods based on plasmid and bacterial databases.