A genomic overview of the population structure of Salmonella

Antimicrobial Resistance in Salmonella enterica Serovar Paratyphi B Variant Java in Poultry from Europe and Latin America

Salmonella enterica serovar Paratyphi B variant Java sequence type 28 is prevalent in poultry and poultry meat. We investigated the evolutionary relatedness between sequence type 28 strains from Europe and Latin America using time-resolved phylogeny and principal component analysis. We sequenced isolates from Colombia, Guatemala, Costa Rica, and the Netherlands and complemented them with publicly available genomes from Europe, Africa, and the Middle East. Phylogenetic time trees and effective population sizes (N_e) showed separate clustering of strains from Latin America and Europe. The separation is estimated to have occurred during the 1980s. N_e of strains increased sharply in Europe around 1995 and in Latin America around 2005. Principal component analysis on noncore genes showed a clear distinction between strains from Europe and Latin America, whereas the plasmid gene content was similar. Regardless of the evolutionary separation, similar features of resistance to β-lactams and quinolones/fluoroquinolones indicated parallel evolution of antimicrobial resistance in both regions.

Genomic profiling of antimicrobial resistance genes in clinical isolates of Salmonella Typhi from patients infected with Typhoid fever in India

The development of multidrug resistance in Salmonella enterica serovar Typhi currently forms a major roadblock for the treatment of enteric fever. This poses a major health problem in endemic regions and extends to travellers returning from developing countries. The appearance of fluoroquinolone non-susceptible strains has resulted in use of ceftriaxone as drug of choice with azithromycin being recommended for uncomplicated cases of typhoid fever. A recent sporadic instance of decreased susceptibility to the latest drug regime has necessitated a detailed analysis of antimicrobial resistance genes and possible relationships with their phenotypes to facilitate selection of future treatment regimes. Whole genome sequencing (WGS) was conducted for 133 clinical isolates from typhoid patients. Sequence output files were processed for pan-genome analysis and prediction of antimicrobial resistance genes. The WGS analyses disclosed the existence of fluoroquinolone resistance conferring mutations in gyrA, gyrB, parC and parE genes of all strains. Acquired resistance determining mechanisms observed included catA1 genes for chloramphenicol resistance, dfrA7, dfrA15, sul1 and sul2 for trimethoprim-sulfamethoxazole and bla_TEM-116/bla_TEM-1B genes for amoxicillin. No resistance determinants were found for ceftriaxone and cefixime. The genotypes were further correlated with their respective phenotypes for chloramphenicol, ampicillin, co-trimoxazole, ciprofloxacin and ceftriaxone. A high correlation was observed between genotypes and phenotypes in isolates of S. Typhi. The pan-genome analysis revealed that core genes were enriched in metabolic functions and accessory genes were majorly implicated in pathogenesis and antimicrobial resistance. The pan-genome of S. Typhi appears to be closed (B_pan  =  0.09) as analysed by Heap’s law. Simpson’s diversity index of 0.51 showed a lower level of genetic diversity among isolates of S. Typhi. Overall, this study augments the present knowledge that WGS can help predict resistance genotypes and eventual correlation with phenotypes, enabling the chance to spot AMR determinants for fast diagnosis and prioritize antibiotic use directly from sequence.

Precise phylogenetic analysis of microbial isolates and genomes from metagenomes using PhyloPhlAn 3.0

Microbial genomes are available at an ever-increasing pace, as cultivation and sequencing become cheaper and obtaining metagenome-assembled genomes (MAGs) becomes more effective. Phylogenetic placement methods to contextualize hundreds of thousands of genomes must thus be efficiently scalable and sensitive from closely related strains to divergent phyla. We present PhyloPhlAn 3.0, an accurate, rapid, and easy-to-use method for large-scale microbial genome characterization and phylogenetic analysis at multiple levels of resolution. PhyloPhlAn 3.0 can assign genomes from isolate sequencing or MAGs to species-level genome bins built from >230,000 publically available sequences. For individual clades of interest, it reconstructs strain-level phylogenies from among the closest species using clade-specific maximally informative markers. At the other extreme of resolution, it scales to large phylogenies comprising >17,000 microbial species. Examples including Staphylococcus aureus isolates, gut metagenomes, and meta-analyses demonstrate the ability of PhyloPhlAn 3.0 to support genomic and metagenomic analyses.

Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles

For the first time, in the literature review we presents the molecular genetic structure of Salmonella Enteritidis populations in Russia, and particularly, in Siberia and the Far East of the country. Pathogen population in Russia has been compared with Salmonella populations circulating in different countries of the world. It has been shown that the microbial population is heterogeneous, but it is possible to identify the dominant and main genotypes, which determine up to 90% of the total population morbidity. The data were obtained as a result of a 30-year monitoring (1988 to 2018) by studying the microbial plasmid profiles. It was shown that the same S. Enteritidis clones circulate throughout Russia, however, their significance in the population morbidity may vary depending on geographic and temporal characteristics. Population is characterized by heterogeneity and relative stability of the plasmid types' structure. At the same time, the population is also specified by variability, reflected as a simple change of the annual number of constantly detected plasmid types, and the appearance of new ones that can play a significant role in the etiology of Salmonella infection.

CTX-M-55-type ESBL-producing Salmonella enterica are emerging among retail meats in Phnom Penh, Cambodia

Background

Salmonella enterica is a leading cause of human gastroenteritis. S. enterica strains that produce ESBLs (ESBL-Salm) remain rare in Europe and North America, but less is known about their prevalence among animal-derived foods in countries with weaker food safety practices and unregulated veterinary antibiotic use.

Objectives

To examine the prevalence and characteristics of ESBL-Salm from retail meats in Phnom Penh, Cambodia.

Methods

We tested fish, pork and chicken from two markets for ESBL- and carbapenemase-producing Salmonella from September-December 2016, using cefotaxime- and ertapenem-supplemented media, respectively. ESBL-Salm were sequenced and their genomes characterized. We performed plasmid conjugation experiments to assess the co-transferability of ESBL-encoding genes and MDR phenotypes.

Results

Twenty-six of 150 fish and meat samples (17%) were positive for ESBL-Salm, including 10/60 fish (17%), 15/60 pork (25%) and 1/30 chicken (3%). Carbapenemase-producing Salmonella strains were not detected. Pork-origin ESBL-Salm were primarily serotypes Rissen (10/15) or a monophasic variant of Typhimurium 4,5,12:i:- (3/15), whereas Saintpaul (3/10) and Newport (4/10) were more common among fish. Most ESBL enzymes were encoded by blaCTX-M-55 genes (24/26) harboured on conjugative IncA/C2 (n = 14) or IncHI2 (n = 10) plasmids. Resistance to up to six additional drug classes was co-transferred by each plasmid type. ESBL-Salm were resistant to almost every antibiotic recommended for severe salmonellosis treatment.

Conclusions

CTX-M-55-type S. enterica are highly prevalent among pork and fish from Phnom Penh markets and their spread appears to be mediated by MDR IncA/C2 and IncHI2 plasmids. Food safety must be improved and veterinary antibiotic use should be regulated to protect public health.

Genomic Comparison of Salmonella Enteritidis Strains Isolated from Laying Hens and Humans in the Abruzzi Region during 2018

Salmonellosis is a major cause of bacterial foodborne infection. Since 2016, an increased number of cases of gastroenteritis caused by Salmonella enterica serovar Enteritidis linked to eggs produced in Poland has been reported in Europe. In Italy, S. Enteritidis is one of the three most commonly reported serotypes, associated mainly with the consumption of contaminated eggs and derived products. In our work, we analysed 61 strains of S. Enteritidis obtained from humans and farms in the Abruzzi region, Italy, in 2018. We used Multiple-Loci Variable-Number Tandem Repeat (VNTR) analysis (MLVA)-based typing and Whole-Genome Sequencing (WGS) tools to identify closely related strains and perform cluster analysis. We found two clusters of genetically similar strains. The first one was present in the local farms and isolated from human cases and had single-linkage distance of no more than two core genes and less than five Single-Nucleotide Polymorphisms (SNPs). The second cluster contained strains isolated from humans and from a dessert (tiramisù) sample that shared identical core genome and were assigned the same SNP address. Cluster 2 isolates were found to be genetically similar to an S. Enteritidis strain from a multi-country outbreak linked to Polish eggs.

Prevalence, antibiotic resistance, PFGE and MLST characterization of Salmonella in swine mesenteric lymph nodes

This study evaluated 250 animals from 25 different processing lots, processed in four slaughterhouses in São Paulo state, Brazil for the presence of Salmonella in the mesenteric lymph nodes (10 g sample of each animal) and characterized the antibiotics resistance profile, the Pulsed Field Gel Electrophoresis - PFGE and Multi Locus Sequence Typing - MLST profiles of selected strains. The pathogen was present in 36.4% (n = 91, CL 95% 30.4-43.4) of samples and 72% (n = 18, CL 95% 50.6-87.9%) of the analyzed lots. The main serovars were S. Typhimurium (n = 23), Salmonella enterica subsp. enterica 1.4,5,12:i:- (n = 17), followed by S. Infantis (n = 12) and S. Havana (n = 11). Twenty-eight strains (30%) were classified as other serovars. Sixty-eight percent of the strains were resistant to Streptomycin and tetracycline, followed by ampicillin and sulphonamides (62.6%), chloramphenicol (56.0%), trimethoprim-sulfamethoxazole (41.8%) and nalidixic acid (40.7%). The antibiotics with lower resistance rates were cephalothin and aztreonam (both with 3.3% resistant), and ceftriaxone and cefepime (both with 7.7%). Multidrug-resistant strains (MDR) accounted for 70.3% of the isolates. Eight strains were submitted to MLST: four S. Typhimurium and one S.1.4,5,12:i:-, all belonging to the ST 19, two Salmonella Infantis, belonging to the ST 32 and one S. Derby, belonging to ST 40. Twenty-one isolates with different antibiotics resistance profiles from the most prevalent serovars were selected for PFGE analysis. Serovar S. Typhimurium (n = 11) revealed 4 pulsotypes and 1 cluster and S. 1.4,5,12:i:- (n = 10) revealed 5 pulsotypes and 4 clusters. The high prevalence of the pathogen, with its high rates of antibiotics resistance and belonging to genetic groups that are often associated with disease in humans, shows that the production chain of pork is a potential source of infection in salmonellosis cases. Therefore, effective preventive measures for pathogen control are needed to reduce the risk of foodborne diseases.

Outbreak of Salmonella Newport associated with internationally distributed raw goats' milk cheese, France, 2018

Raw milk cheeses are commonly consumed in France and are also a common source of foodborne outbreaks (FBOs). Both an FBO surveillance system and a laboratory-based surveillance system aim to detect Salmonella outbreaks. In early August 2018, five familial FBOs due to Salmonella spp. were reported to a regional health authority. Investigation identified common exposure to a raw goats' milk cheese, from which Salmonella spp. were also isolated, leading to an international product recall. Three weeks later, on 22 August, a national increase in Salmonella Newport ST118 was detected through laboratory surveillance. Concomitantly isolates from the earlier familial clusters were confirmed as S. Newport ST118. Interviews with a selection of the laboratory-identified cases revealed exposure to the same cheese, including exposure to batches not included in the previous recall, leading to an expansion of the recall. The outbreak affected 153 cases, including six cases in Scotland. S. Newport was detected in the cheese and in the milk of one of the producer's goats. The difference in the two alerts generated by this outbreak highlight the timeliness of the FBO system and the precision of the laboratory-based surveillance system. It is also a reminder of the risks associated with raw milk cheeses.

An African perspective on the prevalence, fate and effects of carbapenem resistance genes in hospital effluents and wastewater treatment plant (WWTP) final effluents: A critical review

This article provides an overview of the antibiotic era and discovery of earliest antibiotics until the present day state of affairs, coupled with the emergence of carbapenem-resistant bacteria. The ways of response to challenges of antibiotic resistance (AR) such as the development of novel strategies in the search of new antibiotics, designing more effective preventive measures as well as the ecology of AR have been discussed. The applications of plant extract and chemical compounds like nanomaterials which are based on recent developments in the field of antimicrobials, antimicrobial resistance (AMR), and chemotherapy were briefly discussed. The agencies responsible for environmental protection have a role to play in dealing with the climate crisis which poses an existential threat to the planet, and contributes to ecological support towards pathogenic microorganisms. The environment serves as a reservoir and also a vehicle for transmission of antimicrobial resistance genes hence, as dominant inhabitants we have to gain a competitive advantage in the battle against AMR.

A genoserotyping system for a fast and objective identification of Salmonella serotypes commonly isolated from poultry and pork food sectors in Belgium

Humans are mostly contaminated by Salmonella through the consumption of pork- and poultry-derived food products. Therefore, a strict monitoring of Salmonella serotypes in food-producing animals is needed to limit the transmission of the pathogen to humans. Additionally, Salmonella can lead to economic loss in the food sector. Previously, a genoserotyping method using the MOL-PCR and Luminex technology was developed for the identification of the 6 Salmonella serotypes, and their variants, subjected to an official control in the Belgian food sector. In this study, 3 additional assays using the same technology were developed for the rapid and cost-effective detection of 13 dangerous highly invasive serotypes or other serotypes frequently isolated from the Belgian poultry and pork sector, i.e. Agona, Anatum, Brandenburg, Choleraesuis, Derby, Enteritidis vaccine strains, Gallinarum var. Gallinarum/Pullorum, Livingstone, Mbandaka, Minnesota, Ohio, Rissen and Senftenberg. Moreover, the previously developed first MOL-PCR assay was improved for S. Paratyphi B and serogroup O:3 detection. Finally, a Decision Support System hosted by a web application was created for an automatic and objective interpretation of the Luminex raw data. The 3 new assays and the modifications of the first assay were validated with a 100% accuracy, using 553 Salmonella and non-Salmonella strains in total.

Hierarchically Labeled Database Indexing Allows Scalable Characterization of Microbiomes

Increasingly available microbial reference data allow interpreting the composition and function of previously uncharacterized microbial communities in detail, via high-throughput sequencing analysis. However, efficient methods for read classification are required when the best database matches for short sequence reads are often shared among multiple reference sequences. Here, we take advantage of the fact that microbial sequences can be annotated relative to established tree structures, and we develop a highly scalable read classifier, PRROMenade, by enhancing the generalized Burrows-Wheeler transform with a labeling step to directly assign reads to the corresponding lowest taxonomic unit in an annotation tree. PRROMenade solves the multi-matching problem while allowing fast variable-size sequence classification for phylogenetic or functional annotation. Our simulations with 5% added differences from reference indicated only 1.5% error rate for PRROMenade functional classification. On metatranscriptomic data PRROMenade highlighted biologically relevant functional pathways related to diet-induced changes in the human gut microbiome.

A guide to machine learning for bacterial host attribution using genome sequence data

With the ever-expanding number of available sequences from bacterial genomes, and the expectation that this data type will be the primary one generated from both diagnostic and research laboratories for the foreseeable future, then there is both an opportunity and a need to evaluate how effectively computational approaches can be used within bacterial genomics to predict and understand complex phenotypes, such as pathogenic potential and host source. This article applied various quantitative methods such as diversity indexes, pangenome-wide association studies (GWAS) and dimensionality reduction techniques to better understand the data and then compared how well unsupervised and supervised machine learning (ML) methods could predict the source host of the isolates. The study uses the example of the pangenomes of 1203 Salmonella enterica serovar Typhimurium isolates in order to predict 'host of isolation' using these different methods. The article is aimed as a review of recent applications of ML in infection biology, but also, by working through this specific dataset, it allows discussion of the advantages and drawbacks of the different techniques. As with all such sub-population studies, the biological relevance will be dependent on the quality and diversity of the input data. Given this major caveat, we show that supervised ML has the potential to add real value to interpretation of bacterial genomic data, as it can provide probabilistic outcomes for important phenotypes, something that is very difficult to achieve with the other methods.

Importation of Extensively Drug-Resistant Salmonella enterica Serovar Typhi Cases in Ontario, Canada

A strain of extensively drug-resistant (XDR) Salmonella enterica serovar Typhi has caused a large ongoing outbreak in Pakistan since 2016. In Ontario, Canada, 10 cases of mainly bloodstream infections (n = 9) were identified in patients who traveled to Pakistan. Whole-genome sequencing showed that Canadian cases were genetically related to the Pakistan outbreak strain. The appearance of XDR typhoid cases in Ontario prompted a provincial wide alert to physicians to recommend treatment with carbapenems or azithromycin in suspected typhoid cases with travel history to Pakistan.

Genomic Characterization of Escherichia coli Isolates Belonging to a New Hybrid aEPEC/ExPEC Pathotype O153:H10-A-ST10 eae-beta1 Occurred in Meat, Poultry, Wildlife and Human Diarrheagenic Samples

Different surveillance studies (2005–2015) in northwest Spain revealed the presence of eae-positive isolates of Escherichia coli O153:H10 in meat for human consumption, poultry farm, wildlife and human diarrheagenic samples. The aim of this study was to explore the genetic and genomic relatedness between human and animal/meat isolates, as well as the mechanism of its persistence. We also wanted to know whether it was a geographically restricted lineage, or whether it was also reported elsewhere. Conventional typing showed that 32 isolates were O153:H10-A-ST10 fimH54, fimAv_MT78, traT and eae-beta1. Amongst these, 21 were CTX-M-32 or SHV-12 producers. The PFGE XbaI-macrorestriction comparison showed high similarity (>85%). The plasmidome analysis revealed a stable combination of IncF (F2:A-:B-), IncI1 (STunknown) and IncX1 plasmid types, together with non-conjugative Col-like plasmids. The core genome investigation based on the cgMLST scheme from EnteroBase proved close relatedness between isolates of human and animal origin. Our results demonstrate that a hybrid MDR aEPEC/ExPEC of the clonal group O153:H10-A-ST10 (CH11-54) is circulating in our region within different hosts, including wildlife. It seems implicated in human diarrhea via meat transmission, and in the spreading of ESBL genes (mainly of CTX-M-32 type). We found genomic evidence of a related hybrid aEPEC/ExPEC in at least one other country.

Emergence of a Novel Salmonella enterica Serotype Reading Clonal Group Is Linked to Its Expansion in Commercial Turkey Production, Resulting in Unanticipated Human Illness in North America

Two separate human outbreaks of Salmonella enterica serotype Reading occurred between 2017 and 2019 in the United States and Canada, and both outbreaks were linked to the consumption of raw turkey products. In this study, a comprehensive genomic investigation was conducted to reconstruct the evolutionary history of S. Reading from turkeys and to determine the genomic context of outbreaks involving this infrequently isolated Salmonella serotype. A total of 988 isolates of U.S. origin were examined using whole-genome-based approaches, including current and historical isolates from humans, meat, and live food animals. Broadly, isolates clustered into three major clades, with one apparently highly adapted turkey clade. Within the turkey clade, isolates clustered into three subclades, including an "emergent" clade that contained only isolates dated 2016 or later, with many of the isolates from these outbreaks. Genomic differences were identified between emergent and other turkey subclades, suggesting that the apparent success of currently circulating subclades is, in part, attributable to plasmid acquisitions conferring antimicrobial resistance, gain of phage-like sequences with cargo virulence factors, and mutations in systems that may be involved in beta-glucuronidase activity and resistance towards colicins. U.S. and Canadian outbreak isolates were found interspersed throughout the emergent subclade and the other circulating subclade. The emergence of a novel S Reading turkey subclade, coinciding temporally with expansion in commercial turkey production and with U.S. and Canadian human outbreaks, indicates that emergent strains with higher potential for niche success were likely vertically transferred and rapidly disseminated from a common source.IMPORTANCE Increasingly, outbreak investigations involving foodborne pathogens are difficult due to the interconnectedness of food animal production and distribution, and homogeneous nature of industry integration, necessitating high-resolution genomic investigations to determine their basis. Fortunately, surveillance and whole-genome sequencing, combined with the public availability of these data, enable comprehensive queries to determine underlying causes of such outbreaks. Utilizing this pipeline, it was determined that a novel clone of Salmonella Reading has emerged that coincided with increased abundance in raw turkey products and two outbreaks of human illness in North America. The rapid dissemination of this highly adapted and conserved clone indicates that it was likely obtained from a common source and rapidly disseminated across turkey production. Key genomic changes may have contributed to its apparent continued success in commercial turkeys and ability to cause illness in humans.

Molecular epidemiology of Salmonella Infantis in Europe: insights into the success of the bacterial host and its parasitic pESI-like megaplasmid

Salmonella Infantis is one of the five serovars most frequently causing human salmonellosis in Europe, mainly associated with poultry. A clone harbouring a conjugative plasmid of emerging S. Infantis (pESI)-like megaplasmid, carrying multidrug resistant (MDR) and extended-spectrum beta-lactamases (ESBL) genes, has spread in the Italian broiler chicken industry also causing human illness. This work is aimed at elucidating the molecular epidemiology of S. Infantis and pESI-like in Europe using whole-genome sequencing and bioinformatics analysis, and to investigate the genetic relatedness of S. Infantis clones and pESI-like from animals, meat, feed and humans provided by institutions of nine European countries. Two genotyping approaches were used: chromosome or plasmid SNP-based analysis and the minimum spanning tree (MST) algorithm based on core-genome multilocus sequence typing (cgMLST). The European S. Infantis population appeared heterogeneous, with different genetic clusters defined at core-genome level. However, pESI-like variants present in 64.1 % of the isolates were more genetically homogeneous and capable of infecting different clonal lineages in most of the countries. Two different pESI-like with ESBL genes (n=82) were observed: bla_CTX-M-1-positive in European isolates and bla_CTX-M-65-positive in American isolates (study outgroup). Both variants had toxin-antitoxin systems, resistance genes towards tetracyclines, trimethoprim, sulphonamides and aminoglycosides, heavy metals (merA) and disinfectants (qacEΔ). Worryingly, 66 % of the total isolates studied presented different gyrA chromosomal point mutations associated with (fluoro)quinolone resistance (MIC range 0.125–0.5 mg/L), while 18 % displayed transferable macrolide resistance mediated by mph, mef and erm(B) genes. Proper intervention strategies are needed to prevent further dissemination/transmission of MDR S. Infantis and pESI-like along the food chain in Europe.

Metagenomic Analysis of Acquired Antibiotic Resistance Determinants in the Gut Microbiota of Wild Boars (Sus Scrofa) - Preliminary Results

Introduction

Land application of manure that contains antibiotics and resistant bacteria may facilitate the establishment of an environmental reservoir of antibiotic-resistant microbes, promoting their dissemination into agricultural and natural habitats. The main objective of this study was to search for acquired antibiotic resistance determinants in the gut microbiota of wild boar populations living in natural habitats.

Material and Methods

Gastrointestinal samples of free-living wild boars were collected in the Zemplén Mountains in Hungary and were characterised by culture-based, metagenomic, and molecular microbiological methods. Bioinformatic analysis of the faecal microbiome of a hunted wild boar from Japan was used for comparative studies. Also, shotgun metagenomic sequencing data of two untreated sewage wastewater samples from North Pest (Hungary) from 2016 were analysed by bioinformatic methods. Minimum spanning tree diagrams for seven-gene MLST profiles of 104 E. coli strains isolated in Europe from wild boars and domestic pigs were generated in Enterobase.

Results

In the ileum of a diarrhoeic boar, a dominant E. coli O112ab:H2 strain with intermediate resistance to gentamicin, tobramycin, and amikacin was identified, displaying sequence type ST388 and harbouring the EAST1 toxin astA gene. Metagenomic analyses of the colon and rectum digesta revealed the presence of the tetQ, tetW, tetO, and mefA antibiotic resistance genes that were also detected in the gut microbiome of four other wild boars from the mountains. Furthermore, the tetQ and cfxA genes were identified in the faecal microbiome of a hunted wild boar from Japan.

Conclusion

The gastrointestinal microbiota of the free-living wild boars examined in this study carried acquired antibiotic resistance determinants that are highly prevalent among domestic livestock populations.

Alignment- and reference-free phylogenomics with colored de Bruijn graphs

BACKGROUND

The increasing amount of available genome sequence data enables large-scale comparative studies. A common task is the inference of phylogenies-a challenging task if close reference sequences are not available, genome sequences are incompletely assembled, or the high number of genomes precludes multiple sequence alignment in reasonable time.

RESULTS

We present a new whole-genome based approach to infer phylogenies that is alignment- and reference-free. In contrast to other methods, it does not rely on pairwise comparisons to determine distances to infer edges in a tree. Instead, a colored de Bruijn graph is constructed, and information on common subsequences is extracted to infer phylogenetic splits.

CONCLUSIONS

The introduced new methodology for large-scale phylogenomics shows high potential. Application to different datasets confirms robustness of the approach. A comparison to other state-of-the-art whole-genome based methods indicates comparable or higher accuracy and efficiency.

First detection of a plasmid located carbapenem resistant blaVIM-1 gene in E. coli isolated from meat products at retail in Belgium in 2015

Carbapenemase-producing Enterobacteriaceae (CPE) confer resistance to antibiotics that are of critical importance to human medicine. There have only been a few reported cases of CPEs in the European food chain. We report the first detection of a carbapenemase-producing Escherichia coli (ST 5869) in the Belgian food chain. Our aim was to characterize the origin of the carbapenem resistance in the E. coli isolate. The isolate was detected during the screening of 178 minced pork samples and was shown to contain the carbapenemase gene bla_VIM-1 by PCR and Sanger sequencing. Whole genome short and long read sequencing (MiSeq and MinION) was performed to characterize the isolate. With a hybrid assembly we reconstructed a 190,205 bp IncA/C2 plasmid containing bla_VIM-1 (S15FP06257_p), in addition to other critically important resistance genes. This plasmid showed only low similarity to plasmids containing bla_VIM-1 previously reported in Germany. Moreover, no sequences existed in the NCBI nucleotide database that completely covered S15FP06257_p. Analysis of the bla_VIM-1 gene cassette demonstrated that it likely originated from an integron of a Klebsiella plasmid reported previously in a clinical isolate in Europe, suggesting that the meat could have been contaminated by human handling in one of the steps of the food chain. This study shows the relevance of fully reconstructing plasmids to characterize their genetic content and to allow source attribution. This is especially important in view of the potential risk of antimicrobial resistance gene transmission through mobile elements as was reported here for the of public health concern bla_VIM-1.

The Not so Good, the Bad and the Ugly: Differential Bacterial Adhesion and Invasion Mediated by Salmonella PagN Allelic Variants

While advances in genomic sequencing have highlighted significant strain variability between and within Salmonella serovars, only a few protein variants have been directly related to evolutionary adaptation for survival, such as host specificity or differential virulence. The current study investigated whether allelic variation of the Salmonella adhesin/invasin PagN influences bacterial interaction with their receptors. The Salmonella enterica, subspecies enterica serovar Typhi (S. Typhi) allelic variant of PagN was found to bind significantly better to different enterocytes as well as to the extracellular matrix protein laminin than did the major Salmonella enterica, subspecies enterica serovar Typhimurium (S. Typhimurium) allele. The two alleles differed at amino acid residues 49 and 109 in two of the four predicted PagN surface loops, and residue substitution analysis revealed that a glutamic acid at residue 49 increased the adhesive and invasive properties of S. Typhi PagN. PagN sequence comparisons from 542 Salmonella strains for six representative S. enterica serovars and S. diarizonae further supported the role of glutamic acid at residues 49 and 109 in optimizing adhesion to cells and laminin, as well as for cell invasion. In summary, this study characterized unique residues in allelic variants of a virulence factor that participates in the colonization and invasive properties of different Salmonella stains, subspecies and serovars.

Large scale automated phylogenomic analysis of bacterial isolates and the Evergreen Online platform

Public health authorities whole-genome sequence thousands of isolates each month for microbial diagnostics and surveillance of pathogenic bacteria. The computational methods have not kept up with the deluge of data and the need for real-time results. We have therefore created a bioinformatics pipeline for rapid subtyping and continuous phylogenomic analysis of bacterial samples, suited for large-scale surveillance. The data is divided into sets by mapping to reference genomes, then consensus sequences are generated. Nucleotide based genetic distance is calculated between the sequences in each set, and isolates are clustered together at 10 single-nucleotide polymorphisms. Phylogenetic trees are inferred from the non-redundant sequences and the clustered isolates are added back. The method is accurate at grouping outbreak strains together, while discriminating them from non-outbreak strains. The pipeline is applied in Evergreen Online, which processes publicly available sequencing data from foodborne bacterial pathogens on a daily basis, updating phylogenetic trees as needed.

Genome-based Salmonella serotyping as the new gold standard

Salmonella enterica is the second most reported bacterial cause of food-borne infections in Europe. Therefore molecular surveillance activities based on pathogen subtyping are an important measure of controlling Salmonellosis by public health agencies. In Germany, at the federal level, this work is carried out by the National Reference Center for Salmonella and other Bacterial Enteric Pathogens (NRC). With rise of next generation sequencing techniques, the NRC has introduced whole-genome-based typing methods for S. enterica in 2016. In this study we report on the feasibility of genome-based in silico serotyping in the German setting using raw sequence reads. We found that SeqSero and seven gene MLST showed 98% and 95% concordance, respectively, with classical serotyping for the here evaluated serotypes, including the most common German serotypes S. Enteritidis and S. Typhimurium as well as less frequently found serotypes. The level of concordance increased to >99% when the results of both in silico methods were combined. However, both tools exhibited misidentification of monophasic variants, in particular monophasic S. Typhimurium and therefore need to be fine-tuned for reliable detection of this epidemiologically important variant. We conclude that with adjustments Salmonella genome-based serotyping might become the new gold standard.

WhatsGNU: a tool for identifying proteomic novelty

To understand diversity in enormous collections of genome sequences, we need computationally scalable tools that can quickly contextualize individual genomes based on their similarities and identify features of each genome that make them unique. We present WhatsGNU, a tool based on exact match proteomic compression that, in seconds, classifies any new genome and provides a detailed report of protein alleles that may have novel functional differences. We use this technique to characterize the total allelic diversity (panallelome) of Salmonella enterica, Mycobacterium tuberculosis, Pseudomonas aeruginosa, and Staphylococcus aureus. It could be extended to others. WhatsGNU is available from https://github.com/ahmedmagds/WhatsGNU.

Comparative genomic analysis reveals high intra-serovar plasticity within Salmonella Napoli isolated in 2005-2017

BACKGROUND

Salmonella enterica subsp. enterica serovar Napoli (S. Napoli) is among the top serovars causing human infections in Italy, although it is relatively uncommon in other European countries; it is mainly isolated from humans and the environment, but neither the reservoir nor its route of infection are clearly defined. This serovar is characterized by high genomic diversity, and molecular evidences revealed important similarities with typhoidal serovars.

RESULTS

179 S. Napoli genomes as well as 239 genomes of typhoidal and non-typhoidal serovars were analyzed in a comparative genomic study. Phylogenetic analysis and draft genome characterization in terms of Multi Locus Sequence Typing (MLST), plasmid replicons, Salmonella Pathogenicity Islands (SPIs), antimicrobial resistance genes (ARGs), phages, biocide and metal-tolerance genes confirm the high genetic variability of S. Napoli, also revealing a within-serovar phylogenetic structure more complex than previously known. Our work also confirms genomic similarity of S. Napoli to typhoidal serovars (S. Typhi and S. Paratyphi A), with S. Napoli samples clustering primarily according to ST, each being characterized by specific genomic traits. Moreover, two major subclades of S. Napoli can be clearly identified, with ST-474 being biphyletic. All STs span among isolation sources and years of isolation, highlighting the challenge this serovar poses to define its epidemiology and evolution. Altogether, S. Napoli strains carry less SPIs and less ARGs than other non-typhoidal serovars and seldom acquire plasmids. However, we here report the second case of an extended-spectrum β-lactamases (ESBLs) producing S. Napoli strain and the first cases of multidrug resistant (MDR) S. Napoli strains, all isolated from humans.

CONCLUSIONS

Our results provide evidence of genomic plasticity of S. Napoli, highlighting genomic similarity with typhoidal serovars and genomic features typical of non-typhoidal serovars, supporting the possibility of survival in different niches, both enteric and non-enteric. Presence of horizontally acquired ARGs and MDR profiles rises concerns regarding possible selective pressure exerted by human environment on this pathogen.

Four European Salmonella Typhimurium datasets collected to develop WGS-based source attribution methods

Zoonotic Salmonella causes millions of human salmonellosis infections worldwide each year. Information about the source of the bacteria guides risk managers on control and preventive strategies. Source attribution is the effort to quantify the number of sporadic human cases of a specific illness to specific sources and animal reservoirs. Source attribution methods for Salmonella have so far been based on traditional wet-lab typing methods. With the change to whole genome sequencing there is a need to develop new methods for source attribution based on sequencing data. Four European datasets collected in Denmark (DK), Germany (DE), the United Kingdom (UK) and France (FR) are presented in this descriptor. The datasets contain sequenced samples of Salmonella Typhimurium and its monophasic variants isolated from human, food, animal and the environment. The objective of the datasets was either to attribute the human salmonellosis cases to animal reservoirs or to investigate contamination of the environment by attributing the environmental isolates to different animal reservoirs.

Phylogenetic analysis revealed that Salmonella Typhimurium ST313 isolated from humans and food in Brazil presented a high genomic similarity

Salmonella Typhimurium sequence type 313 (S. Typhimurium ST313) has caused invasive disease mainly in sub-Saharan Africa. In Brazil, ST313 strains have been recently described, and there is a lack of studies that assessed by whole genome sequencing (WGS)-the relationship of these strains. The aims of this work were to study the phylogenetic relationship of 70 S. Typhimurium genomes comparing strains of ST313 (n = 9) isolated from humans and food in Brazil among themselves, with other STs isolated in this country (n = 31) and in other parts of the globe (n = 30) by 16S rRNA sequences, the Gegenees software, whole genome multilocus sequence typing (wgMLST), and average nucleotide identity (ANI) for the genomes of ST313. Additionally, pangenome analysis was performed to verify the heterogeneity of these genomes. The phylogenetic analyses showed that the ST313 genomes were very similar among themselves. However, the ST313 genomes were usually clustered more distantly to other STs of strains isolated in Brazil and in other parts of the world. By pangenome calculation, the core genome was 2,880 CDSs and 4,171 CDSs singletons for all the 70 S. Typhimurium genomes studied. Considering the 10 ST313 genomes analyzed the core genome was 4,112 CDSs and 76 CDSs singletons. In conclusion, the ST313 genomes from Brazil showed a high similarity among them which information might eventually help in the development of vaccines and antibiotics. The pangenome analysis showed that the S. Typhimurium genomes studied presented an open pangenome, but specifically tending to become close for the ST313 strains.

Concordance of SNP- and allele-based typing workflows in the context of a large-scale international Salmonella Enteritidis outbreak investigation

A large European multi-country Salmonella enterica serovar Enteritidis outbreak associated with Polish eggs was characterized by whole-genome sequencing (WGS)-based analysis, with various European institutes using different analysis workflows to identify isolates potentially related to the outbreak. The objective of our study was to compare the output of six of these different typing workflows (distance matrices of either SNP-based or allele-based workflows) in terms of cluster detection and concordance. To this end, we analysed a set of 180 isolates coming from confirmed and probable outbreak cases, which were representative of the genetic variation within the outbreak, supplemented with 22 unrelated contemporaneous S. enterica serovar Enteritidis isolates. Since the definition of a cluster cut-off based on genetic distance requires prior knowledge on the evolutionary processes that govern the bacterial populations in question, we used a variety of hierarchical clustering methods (single, average and complete) and selected the optimal number of clusters based on the consensus of the silhouette, Dunn2, and McClain-Rao internal validation indices. External validation was done by calculating the concordance with the WGS-based case definition (SNP-address) for this outbreak using the Fowlkes-Mallows index. Our analysis indicates that with complete-linkage hierarchical clustering combined with the optimal number of clusters, as defined by three internal validity indices, the six different allele- and SNP-based typing workflows generate clusters with similar compositions. Furthermore, we show that even in the absence of coordinated typing procedures, but by using an unsupervised machine learning methodology for cluster delineation, the various workflows that are currently in use by six European public-health authorities can identify concordant clusters of genetically related S. enterica serovar Enteritidis isolates; thus, providing public-health researchers with comparable tools for detection of infectious-disease outbreaks.

Comparative genomics of Salmonella enterica serovar Enteritidis ST-11 isolated in Uruguay reveals lineages associated with particular epidemiological traits

Salmonella enterica serovar Enteritidis is a major cause of foodborne disease in Uruguay since 1995. We used a genomic approach to study a set of isolates from different sources and years. Whole genome phylogeny showed that most of the strains are distributed in two major lineages (E1 and E2), both belonging to MLST sequence type 11 the major ST among serovar Enteritidis. Strikingly, E2 isolates are over-represented in periods of outbreak abundance in Uruguay, while E1 span all epidemic periods. Both lineages circulate in neighbor countries at the same timescale as in Uruguay, and are present in minor numbers in distant countries. We identified allelic variants associated with each lineage. Three genes, ycdX, pduD and hsdM, have distinctive variants in E1 that may result in defective products. Another four genes (ybiO, yiaN, aas, aceA) present variants specific for the E2 lineage. Overall this work shows that S. enterica serovar Enteritidis strains circulating in Uruguay have the same phylogenetic profile than strains circulating in the region, as well as in more distant countries. Based on these results we hypothesize that the E2 lineage, which is more prevalent during epidemics, exhibits a combination of allelic variants that could be associated with its epidemic ability.

The Virulence of Escherichia coli O157:H7 Isolates in Mice Depends on Shiga Toxin Type 2a (Stx2a)-Induction and High Levels of Stx2a in Stool

In this study we compared nine Shiga toxin (Stx)-producing Escherichia coli O157:H7 patient isolates for Stx levels, stx-phage insertion site(s), and pathogenicity in a streptomycin (Str)-treated mouse model. The strains encoded stx_2a, stx_1a and stx_2a, or stx_2a and stx_2c. All of the strains elaborated 10⁵-10⁶ cytotoxic doses 50% (CD₅₀) into the supernatant after growth in vitro as measured on Vero cells, and showed variable levels of increased toxin production after growth with sub-inhibitory levels of ciprofloxacin (Cip). The stx_2a+stx_2c+ isolates were 90–100% lethal for Str-treated BALB/c mice, though one isolate, JH2013, had a delayed time-to-death. The stx_2a+ isolate was avirulent. Both an stx_2a and a recA deletion mutant of one of the stx_2a+stx_2c+ strains, JH2010, exhibited at least a three-log decrease in cytotoxicity in vitro and both were avirulent in the mice. Stool from Str-treated mice infected with the highly virulent isolates were 10- to 100-fold more cytotoxic than feces from mice infected with the clinical isolate, JH2012, that made only Stx2a. Taken together these findings demonstrate that the stx_2a-phage from JH2010 induces to higher levels in vivo than does the phage from JH2012. The stx_1a+stx_2a+ clinical isolates were avirulent and neutralization of Stx1 in stool from mice infected with those strains indicated that the toxin produced in vivo was primarily Stx1a. Treatment of mice infected with Stx1a+Stx2a+ isolates with Cip resulted in an increase in Stx2a production in vivo and lethality in the mice. Our data suggest that high levels of Stx2a in stool are predictive of virulence in mice.

Emergence of human-adapted Salmonella enterica is linked to the Neolithization process

It has been hypothesized that the Neolithic transition towards an agricultural and pastoralist economy facilitated the emergence of human adapted pathogens. Here, we recovered eight Salmonella enterica subsp. enterica genomes from human skeletons of transitional foragers, pastoralists, and agro-pastoralists in western Eurasia that were up to 6,500 years old. Despite the high genetic diversity of S. enterica all ancient bacterial genomes clustered in a single previously uncharacterized branch that contains S. enterica adapted to multiple mammalian species. All ancient bacterial genomes from prehistoric (agro-)pastoralists fall within a part of this branch that also includes the human-specific S. enterica Paratyphi C, illustrating the evolution of a human pathogen over a period of five thousand years. Bacterial genomic comparisons suggest that the earlier ancient strains were not host specific, differed in pathogenic potential, and experienced convergent pseudogenization that accompanied their downstream host adaptation. These observations support the concept that the emergence of human adapted S. enterica is linked to human cultural transformations.

Understanding the challenges of non-food industrial product contamination

Preventing microbial contamination of non-food products is a major area of industrial microbiology where preservatives are used to stop microbial growth. However, microorganisms occasionally overcome product preservation, causing recalls and the implementation of multiple procedures to prevent further contamination. Correct reporting of microbial contamination in non-food industrial products is vital, especially if spoilage organisms are antimicrobial resistant and pose a health threat. Gram-negative bacteria such as Pseudomonas, Burkholderia and Enterobacteriaceae are frequently reported as non-food product contaminants, including species that overlap current antimicrobial resistance priorities. Historical analysis of recall databases highlighted that for greater than 15% of contamination incidents, the causative microbial agents are reported as unidentified. Here we review the current antimicrobial resistant bacterial species associated with non-food product contamination and evaluate recall reporting in Europe from 2005 to 2018. Our review shows that 49% of microbial contaminants are reported as unidentified despite frequent detection of antimicrobial resistant pathogens; in contrast, 98% of food-related microbial contaminants are classified. Recommendations to fill this microbial identification gap in non-food product recalls are made. Overall, reporting standards for microbial contamination in non-food products must be improved to enable surveillance and for understanding the risks associated with antimicrobial resistant microorganisms.

Escherichia coli ST302: Genomic Analysis of Virulence Potential and Antimicrobial Resistance Mediated by Mobile Genetic Elements

aEPEC are associated with persistent diarrhea, and diarrheal outbreaks in both humans and animals worldwide. They are differentiated from typical EPEC by the lack of bundle-forming pili, and from EHEC by the lack of phage-mediated stx toxins. However, phylogenetic analyses often associate aEPEC with EHEC, promoting the hypothesis that aEPEC are the progenitors of EHEC, which is supported by aEPEC conversion to EHEC by stx-carrying phages. While aEPEC can cause disease outright, the potential to acquire stx, one of the most potent bacterial toxins known, merits close monitoring. Escherichia coli ST302 (O108:H9, O182:H9, O45:H9) are aEPEC that have been isolated from diarrheic human, pig and rabbit hosts, as well as in healthy pigs, however, no study to date has focused on E. coli ST302 strains. Through WGS and hybrid assembly we present the first closed chromosome, and two circularized plasmids of an ST302 strain - F2_18C, isolated from a healthy pig in Australia. A phylogenetic analysis placed E. coli ST302 strains in proximity to EHEC ST32 (O145:H28) strains. Public databases were interrogated for WGSs of E. coli ST302 strains and short-read gene screens were used to compare their virulence-associated gene (VAG) and antimicrobial resistance gene (ARG) cargo. E. coli ST302 strains carry diverse VAGs, including those that typically associated with extraintestinal pathogenic E. coli (ExPEC). Plasmid comparisons showed that pF2_18C_FIB shared homology with EHEC virulence plasmids such as pO103 while pF2_18C_HI2 is a large multidrug resistance IncHI2:ST3 plasmid. A comparison of 33 HI2:ST3 plasmids demonstrated that those of Australian origin have not acquired resistances to extended-spectrum beta-lactams, colistin, fosfomycin or rifampicin, unlike those originating from Asia. F2_18C was shown to carry two additional pathogenicity islands – ETT2, and the STEC-associated PAI_CL3, plasmid-associated heavy metal resistance genes, as well as several unoccupied stx-phage attachment sites. This study sheds light on the virulence and AMR potential of E. coli ST302 strains and informs AMR genomic surveillance.

Epidemiology of Salmonella enterica Serovar Dublin in Cattle and Humans in Denmark, 1996 to 2016: a Retrospective Whole-Genome-Based Study

Salmonella enterica serovar Dublin is a cattle-adapted S. enterica serovar causing both intestinal and systemic infection in its bovine host, and it is also a serious threat to human health. The present study aimed to determine the population structure of S Dublin isolates obtained from Danish cattle herds and to investigate how cattle isolates relate to Danish human isolates, as well as to non-Danish human and bovine isolates. Phylogenetic analysis of 197 Danish cattle isolates from 1996 to 2016 identified three major clades corresponding to distinct geographical regions of cattle herds. Persistence of closely related isolates within the same herd and their circulation between epidemiologically linked herds for a period of more than 20 years were demonstrated. These findings suggest that a lack of internal biosecurity and, to some extent, also a lack of external biosecurity in the herds have played an important role in the long-term persistence of S Dublin in Danish cattle herds in the period investigated. Global population analysis revealed that Danish cattle isolates clustered separately from bovine isolates from other countries, whereas human isolates were geographically spread. Resistance genes were not commonly demonstrated in Danish bovine isolates; only the isolates within one Danish clade were found to often harbor two plasmids of IncFII/IncFIB and IncN types, the latter plasmid carrying bla _TEM-1, tetA, strA, and strB antibiotic resistance genes.IMPORTANCE S Dublin causes economic losses in cattle production, and the bacterium is a public health concern. A surveillance and control program has been in place in Denmark since 2002 with the ultimate goal to eradicate S Dublin from Danish cattle herds; however, a small proportion of herds have remained positive for many years. In this study, we demonstrate that herds with persistent infection often were infected with the same strain for many years, indicating that internal biosecurity has to be improved to curb the infection. Further, domestic cases of S Dublin infection in humans were found to be caused both by Danish cattle isolates and by isolates acquired abroad. This study shows the strength of whole-genome sequencing to obtain detailed information on epidemiology of S Dublin and allows us to suggest internal biosecurity as a main way to control this bacterium in Danish cattle herds.

Distinct but Intertwined Evolutionary Histories of Multiple Salmonella enterica Subspecies

S. enterica is a major foodborne pathogen, which can be transmitted via several distinct routes from animals and environmental sources to human hosts. Multiple subspecies and serotypes of S. enterica exhibit considerable differences in virulence, host specificity, and colonization. This study provides detailed insights into the dynamics of recombination and its contributions to S. enterica subspecies evolution. Widespread recombination within the species means that new adaptations arising in one lineage can be rapidly transferred to another lineage. We therefore predict that recombination has been an important factor in the emergence of several major disease-causing strains from diverse genomic backgrounds and their ability to adapt to disparate environments.

Salmonella is responsible for many nontyphoidal foodborne infections and enteric (typhoid) fever in humans. Of the two Salmonella species, Salmonella enterica is highly diverse and includes 10 known subspecies and approximately 2,600 serotypes. Understanding the evolutionary processes that generate the tremendous diversity in Salmonella is important in reducing and controlling the incidence of disease outbreaks and the emergence of virulent strains. In this study, we aim to elucidate the impact of homologous recombination in the diversification of S. enterica subspecies. Using a data set of previously published 926 Salmonella genomes representing the 10 S. enterica subspecies and Salmonella bongori, we calculated a genus-wide pan-genome composed of 84,041 genes and the S. enterica pan-genome of 81,371 genes. The size of the accessory genomes varies between 12,429 genes in S. enterica subsp. arizonae (subsp. IIIa) to 33,257 genes in S. enterica subsp. enterica (subsp. I). A total of 12,136 genes in the Salmonella pan-genome show evidence of recombination, representing 14.44% of the pan-genome. We identified genomic hot spots of recombination that include genes associated with flagellin and the synthesis of methionine and thiamine pyrophosphate, which are known to influence host adaptation and virulence. Last, we uncovered within-species heterogeneity in rates of recombination and preferential genetic exchange between certain donor and recipient strains. Frequent but biased recombination within a bacterial species may suggest that lineages vary in their response to environmental selection pressure. Certain lineages, such as the more uncommon non-enterica subspecies (non-S. enterica subsp. enterica), may also act as a major reservoir of genetic diversity for the wider population.

IMPORTANCES. enterica is a major foodborne pathogen, which can be transmitted via several distinct routes from animals and environmental sources to human hosts. Multiple subspecies and serotypes of S. enterica exhibit considerable differences in virulence, host specificity, and colonization. This study provides detailed insights into the dynamics of recombination and its contributions to S. enterica subspecies evolution. Widespread recombination within the species means that new adaptations arising in one lineage can be rapidly transferred to another lineage. We therefore predict that recombination has been an important factor in the emergence of several major disease-causing strains from diverse genomic backgrounds and their ability to adapt to disparate environments.

The EnteroBase user's guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity

EnteroBase is an integrated software environment that supports the identification of global population structures within several bacterial genera that include pathogens. Here, we provide an overview of how EnteroBase works, what it can do, and its future prospects. EnteroBase has currently assembled more than 300,000 genomes from Illumina short reads from Salmonella, Escherichia, Yersinia, Clostridioides, Helicobacter, Vibrio, and Moraxella and genotyped those assemblies by core genome multilocus sequence typing (cgMLST). Hierarchical clustering of cgMLST sequence types allows mapping a new bacterial strain to predefined population structures at multiple levels of resolution within a few hours after uploading its short reads. Case Study 1 illustrates this process for local transmissions of Salmonella enterica serovar Agama between neighboring social groups of badgers and humans. EnteroBase also supports single nucleotide polymorphism (SNP) calls from both genomic assemblies and after extraction from metagenomic sequences, as illustrated by Case Study 2 which summarizes the microevolution of Yersinia pestis over the last 5000 years of pandemic plague. EnteroBase can also provide a global overview of the genomic diversity within an entire genus, as illustrated by Case Study 3, which presents a novel, global overview of the population structure of all of the species, subspecies, and clades within Escherichia.

A proposed core genome scheme for analyses of the Salmonella genus

The salmonellae are found in a wide range of animal hosts and many food products for human consumption. Most cases of human disease are caused by S. enterica subspecies I; however as opportunistic pathogens the other subspecies (II-VI) and S. bongori are capable of causing disease. Loci that were not consistently present in all of the species and subspecies were removed from a previously proposed core genome scheme (EBcgMLSTv2.0), the removal of these 252 loci resulted in a core genus scheme (SalmcgMLSTv1.0). SalmcgMLSTv1.0 clustered isolates from the same subspecies more rapidly and more accurately grouped isolates from different subspecies when compared with EBcgMLSTv2.0. All loci within the EBcgMLSTv2.0 scheme were present in over 98% of S. enterica subspecies I isolates and should, therefore, continue to be used for subspecies I analyses, while the SalmcgMLSTv1.0 scheme is more appropriate for cross genus investigations.

Coresistance to quaternary ammonium compounds in extended-spectrum beta-lactamase-producing Escherichia coli

Background and Aim: Extended-spectrum β-lactamases (ESBL) in Escherichia coli constitutes one of the major threats to modern medicine, and the increasing pollution with quaternary ammonium compounds (QACs) has been suspected to contribute to the spread of ESBL-producing bacteria. The aim of the study was to investigate ESBLA and ESBLM-C-producing E. coli isolates for their coresistance to QACs and their phylogeny isolated from a Swedish University Hospital. Materials and Methods: Coresistance in E. coli with production of ESBL enzymes of the type blaCTX-M (n=23) was compared to E. coli producing AmpC type ESBL enzymes blaCMY and blaDHA (n=27). All isolates were tested for susceptibility to antibiotics and QACs, and high-quality whole-genome sequences were analyzed for resistance determinants. Results: The plasmid-borne small multidrug resistance (SMR) efflux pump sugE(p) was solely present in blaCMY-producing E. coli (n=9), within the same genetic environment blaCMY–blc–sugE(p). Other small multidrug efflux pumps were found without association for ESBL-types: emrE (n=5) and the truncated qacEΔ1 (n=18). Conclusion: Coresistance of ESBL enzymes and SMR efflux pumps in E. coli was common and might indicate that other substances than antibiotics contribute to the spread and emergence of antibiotic resistance.

Decreasing prevalence of contamination with extended-spectrum beta-lactamase-producing Enterobacteriaceae (ESBL-E) in retail chicken meat in the Netherlands

Retail chicken meat is a potential source of extended-spectrum beta-lactamase-producing Enterobacteriaceae (ESBL-E). In the past decade, vast national efforts were undertaken to decrease the antibiotic use in the veterinary sector, resulting in a 58% decrease in antibiotic sales in the sector between 2009 and 2014. This decrease in antibiotic use was followed by a decrease in ESBL-E prevalence in broilers. The current study investigates the prevalence of contamination with ESBL-E in retail chicken meat purchased in the Netherlands between December 2013 and August 2015. It looks at associations between the prevalence of contamination with ESBL-E and sample characteristics such as method of farming (free-range or conventional), supermarket chain of purchase and year of purchase. In the current study, 352 chicken meat samples were investigated for the presence of ESBL-E using selective culture methods. Six samples were excluded due to missing isolates or problems obtaining a good quality sequence leaving 346 samples for further analyses. Of these 346 samples, 188 (54.3%) were positive for ESBL-E, yielding 216 ESBL-E isolates (Escherichia coli (n = 204), Klebsiella pneumoniae (n = 11) and Escherichia fergusonii (n = 1)). All ESBL-E isolates were analysed using whole-genome sequencing. The prevalence of contamination with ESBL-E in retail chicken meat decreased from 68.3% in 2014 to 44.6% in 2015, absolute risk difference 23.7% (95% confidence interval (CI): 12.6% - 34.1%). The ESBL-E prevalence was lower in free-range chicken meat (36.4%) compared with conventional chicken meat (61.5%), absolute risk difference 25.2% (95% CI: 12.9% - 36.5%). The prevalence of contamination with ESBL-E varied between supermarket chains, the highest prevalence of contamination was found in supermarket chain 4 (76.5%) and the lowest in supermarket chain 1 (37.8%). Pairwise isolate comparisons using whole-genome multilocus sequence typing (wgMLST) showed that clustering of isolates occurs more frequently within supermarket chains than between supermarket chains. In conclusion, the prevalence of contamination with ESBL-E in retail chicken in the Netherlands decreased over time; nevertheless, it remains substantial and as such a potential source for ESBL-E in humans.

Hash-Based Core Genome Multilocus Sequence Typing for Clostridium difficile

Pathogen whole-genome sequencing has huge potential as a tool to better understand infection transmission. However, rapidly identifying closely related genomes among a background of thousands of other genomes is challenging. Here, we describe a refinement to core genome multilocus sequence typing (cgMLST) in which alleles at each gene are reproducibly converted to a unique hash, or short string of letters (hash-cgMLST). This avoids the resource-intensive need for a single centralized database of sequentially numbered alleles. We test the reproducibility and discriminatory power of cgMLST/hash-cgMLST compared to those of mapping-based approaches in Clostridium difficile, using repeated sequencing of the same isolates (replicates) and data from consecutive infection isolates from six English hospitals. Hash-cgMLST provided the same results as standard cgMLST, with minimal performance penalty. Comparing 272 replicate sequence pairs using reference-based mapping, there were 0, 1, or 2 single-nucleotide polymorphisms (SNPs) between 262 (96%), 5 (2%), and 1 (<1%) of the pairs, respectively. Using hash-cgMLST, 218 (80%) of replicate pairs assembled with SPAdes had zero gene differences, and 31 (11%), 5 (2%), and 18 (7%) pairs had 1, 2, and >2 differences, respectively. False gene differences were clustered in specific genes and associated with fragmented assemblies, but were reduced using the SKESA assembler. Considering 412 pairs of infections with ≤2 SNPS, i.e., consistent with recent transmission, 376 (91%) had ≤2 gene differences and 16 (4%) had ≥4. Comparing a genome to 100,000 others took <1 min using hash-cgMLST. Hash-cgMLST is an effective surveillance tool for rapidly identifying clusters of related genomes. However, cgMLST/hash-cgMLST generate more false variants than mapping-based approaches. Follow-up mapping-based analyses are likely required to precisely define close genetic relationships.

The Occurrence and Characterization of Extended-Spectrum-Beta-Lactamase-Producing Escherichia coli Isolated from Clinical Diagnostic Specimens of Equine Origin

Escherichia coli isolates were recovered from clinical specimens of equine patients admitted to the Texas A&M Veterinary Medical Teaching Hospital over a five-year period. Ceftiofur resistance was used as a marker for potential extended-spectrum beta-lactamase (ESBL)-activity, and of the 48 ceftiofur-resistant E. coli isolates, 27.08% (n = 13) were phenotypically ESBL-positive. Conventional PCR analysis followed by the _large-scalebla Finder multiplex PCR detected the ESBL genes, CTX-M-1 and SHV, in seven out of the 13 isolates. Moreover, beta-lactamase genes of TEM-1-type, BER-type (AmpC), and OXA-type were also identified. Sequencing of these genes resulted in identification of a novel TEM-1-type gene, called bla_TEM-233, and a study is currently underway to determine if this gene confers the ESBL phenotype. Furthermore, this report is the first to have found E. coli ST1308 in horses. This subtype, which has been reported in other herbivores, harbored the SHV-type ESBL gene. Finally, one out of 13 E. coli isolates was PCR-positive for the carbapenemase gene, bla_IMP-1 despite the lack of phenotypically proven resistance to imipenem. With the identification of novel ESBL gene variant and the demonstrated expansion of E. coli sequence types in equine patients, this study underscores the need for more investigation of equines as reservoirs for ESBL-producing pathogens.

Genotypic characterization and antimicrobial resistance profile of Salmonella isolated from chicken, pork and the environment at abattoirs and supermarkets in Chongqing, China

Background

Salmonella is one of the most important foodborne pathogens, causing outbreaks of human salmonellosis worldwide. Owing to large scales of consumption markets, pork and poultry that contaminated by Salmonella could pose a tremendous threat to public health. The aim of this study was to investigate the contamination of Salmonella from chicken, pork and the environment in slaughtering and retail processes in Chongqing, China.

Results

A total of 115 Salmonella isolates were recovered from 1112 samples collected from pork, chicken and the environment. Compared with the isolation rate of samples from chicken (9.50%) and the environment (6.23%), samples from pork had a significant higher isolation rate (44.00%). The isolation rates in slaughterhouses (10.76%) and in supermarkets (10.07%) showed no statistical difference. Thirty different serotypes were identified among all the isolates. S. Derby (n = 26), S. London (n = 16) and S. Rissen (n = 12) were the dominant serotypes. Antimicrobial susceptibility testing revealed that 73.04% isolates were resistant to tetracycline, followed by 66.96% to ampicillin and 59.13% to doxycycline. More than half (50.43%) of the isolates were multidrug resistant (MDR), and most of the MDR isolates were from supermarkets. Multilocus sequence typing results showed 24 out of 115 isolates were ST40, which was the most prevalent. Furthermore, isolates from supermarkets had 20 different sequence types while isolates from slaughterhouses only had 8 different sequence types.

Conclusion

Our study highlighted that Salmonella was more frequently isolated in pork production chain than that in chicken. Compared with isolates from slaughterhouses, isolates from supermarkets had more MDR profiles and represented a wider range of serotypes and sequence types, indicating that the retail process had more diverse sources of Salmonella contamination than that of slaughtering process.

Multidrug-Resistant Salmonella enterica Serovar Rissen Clusters Detected in Azores Archipelago, Portugal

Gastrointestinal infections caused by nontyphoidal Salmonella (NTS) remain one of the main causes of foodborne illness worldwide. Within the multiple existing Salmonella enterica serovars, the serovar Rissen is rarely reported, particularly as a cause of human salmonellosis. Between 2015 and 2017, the Portuguese National Reference Laboratory of Gastrointestinal Infections observed an increase in the number of clinical cases caused by multidrug-resistant (MDR) S. enterica serovar Rissen, particularly from the Azores archipelago. In the present study, we analyzed by whole genome sequencing (WGS) all clinical, animal, food, and environmental isolates received up to 2017 in the Portuguese Reference Laboratories. As such, through a wgMLST-based gene-by-gene analysis, we aimed to identify potential epidemiological clusters linking clinical and samples from multiple sources, while gaining insight into the genetic diversity of S. enterica serovar Rissen. We also investigated the genetic basis driving the observed multidrug resistance. By integrating 60 novel genomes with all publicly available serovar Rissen genomes, we observed a low degree of genetic diversity within this serovar. Nevertheless, the majority of Portuguese isolates showed high degree of genetic relatedness and a potential link to pork production. An in-depth analysis of these isolates revealed the existence of two major clusters from the Azores archipelago composed of MDR isolates, most of which were resistant to at least five antimicrobials. Considering the well-known spread of MDR between gastrointestinal bacteria, the identification of MDR circulating clones should constitute an alert to public health authorities. Finally, this study constitutes the starting point for the implementation of the "One Health" approach for Salmonella surveillance in Portugal.

Prevalence and genomic characteristics of zoonotic gastro-intestinal pathogens and ESBL/pAmpC producing Enterobacteriaceae among Swedish corvid birds

Introduction: Wild birds pose a potential threat to animal and human health by spreading infectious diseases. In the present study, we studied the occurrence of bacterial zoonotic pathogens as well as enterobacteria with transferrable antimicrobial resistance genes among Swedish corvids. Materials and methods: Intestines from 66 jackdaws, crows, rooks and magpies from the vicinity of livestock farms at 14 locations in 7 counties were analysed by direct culture or PCR screening followed by culture. Isolates were investigated by whole-genome sequencing. Results and discussion: Campylobacter jejuni were detected in 82% and Yersinia in 3% of the birds. ESBL-producing E. coli were found in one sample (2%) and carried bla CTX-M-55. No Enterobacteriaceae with transferable carbapenem resistance were identified. No Salmonella or E. coli O157:H7 were found, but PCR analysis for enterohaemorrhagic E. coli virulence genes revealed 35% positive samples for intimin, 9% for verotoxin 1 and 17% for verotoxin 2. C. jejuni isolates from corvids were compared to previously published isolates from Swedish sources by multi-locus sequence typing based on genome sequences. All corvid C. jejuni isolates formed a cluster, intermingled with human and chicken isolates. Our results indicate that C. jejuni is ubiquitous among Swedish corvid birds, with sporadic transmission to poultry and humans.

Multiplex ligation reaction based on probe melting curve analysis: a pragmatic approach for the identification of 30 common Salmonella serovars

BACKGROUND

While Salmonella serotyping is of paramount importance for the disease intervention of salmonellosis, a fast and easy-to-operate molecular serotyping solution remains elusive. We have developed a multiplex ligation reaction based on probe melting curve analysis (MLMA) for the identification of 30 common Salmonella serovars.

METHODS

Serovar-specific primers and probes were designed based on a comparison of gene targets (wzx and wzy encoding for somatic antigen biosynthesis; fliC and fljB for flagellar antigens) from 5868 Salmonella genomes. The ssaR gene, a type III secretion system component, was included for the confirmation of Salmonella.

RESULTS

All gene targets were detected and gave expected Tm values during assay evaluation. Cross reactions were not demonstrated between the 30 serovars (n = 211), or with an additional 120 serovars (n = 120) and other Enterobacteriaceae (n = 3). The limit of identification for all targets ranged from using 1.2 ng/μL to 1.56 ng/μL of DNA. The intra- and inter-assay standard deviations and the coefficients of variation were no more than 0.5 °C and less than 1% respectively, indicating high reproducibility. From consecutive outpatient stool samples (n = 3590) collected over a 10-month period at 11 sentinel hospitals in Shenzhen, China, we conducted a multicenter study using the traditional Salmonella identification workflow and the MLMA assay workflow in parallel. From Salmonella isolates (n = 496, 13.8%) derived by both workflows, total agreement (kappa = 1.0) between the MLMA assay and conventional serotyping was demonstrated.

CONCLUSIONS

With an assay time of 2.5 h, this simple assay has shown promising potential to provide rapid and high-throughput identification of Salmonella serovars for clinical and public health laboratories to facilitate timely surveillance of salmonellosis.

Koutsoumanis K, Allende A, Alvarez-Ordóñez A, Bolton D, Bover-Cid S, Chemaly M, Davies R, De Cesare A, Hilbert F, Lindqvist R, Nauta M, Peixe L, Ru G, Simmons M, Skandamis P, Suffredini E, Jenkins C, Malorny B, Ribeiro Duarte AS, Torpdahl M, da Silva Felício MT, Guerra B, Rossi M, Herman L. Whole genome sequencing and metagenomics for outbreak investigation, source attribution and risk assessment of food-borne microorganisms

This Opinion considers the application of whole genome sequencing (WGS) and metagenomics for outbreak investigation, source attribution and risk assessment of food‐borne pathogens. WGS offers the highest level of bacterial strain discrimination for food‐borne outbreak investigation and source‐attribution as well as potential for more precise hazard identification, thereby facilitating more targeted risk assessment and risk management. WGS improves linking of sporadic cases associated with different food products and geographical regions to a point source outbreak and can facilitate epidemiological investigations, allowing also the use of previously sequenced genomes. Source attribution may be favoured by improved identification of transmission pathways, through the integration of spatial‐temporal factors and the detection of multidirectional transmission and pathogen–host interactions. Metagenomics has potential, especially in relation to the detection and characterisation of non‐culturable, difficult‐to‐culture or slow‐growing microorganisms, for tracking of hazard‐related genetic determinants and the dynamic evaluation of the composition and functionality of complex microbial communities. A SWOT analysis is provided on the use of WGS and metagenomics for Salmonella and Shigatoxin‐producing Escherichia coli (STEC) serotyping and the identification of antimicrobial resistance determinants in bacteria. Close agreement between phenotypic and WGS‐based genotyping data has been observed. WGS provides additional information on the nature and localisation of antimicrobial resistance determinants and on their dissemination potential by horizontal gene transfer, as well as on genes relating to virulence and biological fitness. Interoperable data will play a major role in the future use of WGS and metagenomic data. Capacity building based on harmonised, quality controlled operational systems within European laboratories and worldwide is essential for the investigation of cross‐border outbreaks and for the development of international standardised risk assessments of food‐borne microorganisms.

Multidrug-Resistant and Clinically Relevant Gram-Negative Bacteria Are Present in German Surface Waters

Water is considered to play a role in the dissemination of antibiotic-resistant Gram-negative bacteria including those encoding Extended-spectrum beta-lactamases (ESBL) and carbapenemases. To investigate the role of water for their spread in more detail, we characterized ESBL/Carbapenemase-producing bacteria from surface water and sediment samples using phenotypic and genotypic approaches. ESBL/Carbapenemase-producing isolates were obtained from water/sediment samples. Species and antibiotic resistance were determined. A subset of these isolates (n = 33) was whole-genome-sequenced and analyzed for the presence of antibiotic resistance genes and virulence determinants. Their relatedness to isolates associated with human infections was investigated using multilocus sequence type and cgMLST-based analysis. Eighty-nine percent of the isolates comprised of clinically relevant species. Fifty-eight percent exhibited a multidrug-resistance phenotype. Two isolates harbored the mobile colistin resistance gene mcr-1. One carbapenemase-producing isolate identified as Enterobacter kobei harbored bla VIM- 1. Two Escherichia coli isolates had sequence types (ST) associated with human infections (ST131 and ST1485) and a Klebsiella pneumoniae isolate was classified as hypervirulent. A multidrug-resistant (MDR) Pseudomonas aeruginosa isolate encoding known virulence genes associated with severe lung infections in cystic fibrosis patients was also detected. The presence of MDR and clinically relevant isolates in recreational and surface water underlines the role of aquatic environments as both reservoirs and hot spots for MDR bacteria. Future assessment of water quality should include the examination of the multidrug resistance of clinically relevant bacterial species and thus provide an important link regarding the spread of MDR bacteria in a One Health context.

Vertical transmission of Salmonella Enteritidis with heterogeneous antimicrobial resistance from breeding chickens to commercial chickens in China

Human salmonellosis caused by the consumption of eggs and chicken meat contaminated with Salmonella Enteritidis has become a continuing public health concern worldwide. In this study we adopted whole genome sequencing (WGS) to determine the genetic relationship and antimicrobial resistance of S. enterica strains isolated from a poultry breeding enterprise that consists of one breeding chicken farm, one egg hatchery and one commercial chicken farm. A total of 148 S. enterica including 147 S. Enteritidis strains were isolated from 2100 fecal swab samples, with 16 (5.3 %, 16/300) from breeding chicken farm, 38 (4.2 %, 38/900) from egg hatchery and 94 (10.4 %, 94/900) from commercial chicken farm. WGS revealed that all 147 S. Enteritidis strains belonged to ST11, and further divided into 4 different ribosomal STs and 64 core genome STs. Single nucleotide polymorphism typing suggested the presence of the vertical transmission of S. Enteritidis from breeding chicken to commercial chicken. Three different antimicrobial-resistant plasmids including one bla_CTX-M-14-carrying plasmid and two virulence-resistance plasmids were characterized, resulting in the heterogeneous antimicrobial resistance of clonally related S. Enteritidis strains. Routine surveillance in breeding chicken farms is conducive to the control of S. Enteritidis from farm to fork.

Implications of Mobile Genetic Elements for Salmonella enterica Single-Nucleotide Polymorphism Subtyping and Source Tracking Investigations

Single-nucleotide polymorphisms (SNPs) are widely used for whole-genome sequencing (WGS)-based subtyping of foodborne pathogens in outbreak and source tracking investigations. Mobile genetic elements (MGEs) are commonly present in bacterial genomes and may affect SNP subtyping results if their evolutionary history and dynamics differ from that of the bacterial chromosomes. Using Salmonella enterica as a model organism, we surveyed major categories of MGEs, including plasmids, phages, insertion sequences, integrons, and integrative and conjugative elements (ICEs), in 990 genomes representing 21 major serotypes of S. enterica We evaluated whether plasmids and chromosomal MGEs affect SNP subtyping with 9 outbreak clusters of different serotypes found in the United States in 2018. The median total length of chromosomal MGEs accounted for 2.5% of a typical S. enterica chromosome. Of the 990 analyzed S. enterica isolates, 68.9% contained at least one assembled plasmid sequence. The median total length of assembled plasmids in these isolates was 93,671 bp. Plasmids that carry high densities of SNPs were found to substantially affect both SNP phylogenies and SNP distances among closely related isolates if they were present in the reference genome for SNP subtyping. In comparison, chromosomal MGEs were found to have limited impact on SNP subtyping. We recommend the identification of plasmid sequences in the reference genome and the exclusion of plasmid-borne SNPs from SNP subtyping analysis.IMPORTANCE Despite increasingly routine use of WGS and SNP subtyping in outbreak and source tracking investigations, whether and how MGEs affect SNP subtyping has not been thoroughly investigated. Besides chromosomal MGEs, plasmids are frequently entangled in draft genome assemblies and yet to be assessed for their impact on SNP subtyping. This study provides evidence-based guidance on the treatment of MGEs in SNP analysis for Salmonella to infer phylogenetic relationship and SNP distance between isolates.

[Bioinformatics of next generation sequencing in clinical microbiology diagnosis]

Massive parallel sequencing (High-Throughput Sequencing [HTS]) allows to read millions or billions of DNA sequences or fragments (reads) in parallel and is revolutionizing microbiology research, moving from laboratory methods to computed-assisted analyses, with the compelling use of Bioinformatics. The time and cost reduction in studies on the microbiota, microbiome and metagenome, allows to rapidly progress in diagnosis, taxonomy, epidemiology, comparative genomics, virulence, discovery of genes or variants of interest and the association of microorganisms with traditionally considered non-microbial diseases. In this review, the terminology, the sequencing technologies and their applications are described for microbial analysis using open-source bioinformatics software, analysis pipelines, databases and web platforms that allow a user-friendly bioinformatics approach affordable by the clinical microbiologist and infectious disease practitioners.