A genomic overview of the population structure of Salmonella

Integrated Whole-Genome Sequencing Infrastructure for Outbreak Detection and Source Tracing of Salmonella enterica Serotype Enteritidis

As an important foodborne pathogen, Salmonella enterica serotype Enteritidis is recognized as one of the most common causes of human salmonellosis globally. Outbreak detection for this highly homogenous serotype, however, has remained challenging. Rapid advances in sequencing technologies have presented whole-genome sequencing (WGS) as a significant advancement for source tracing and molecular typing of foodborne pathogens. A retrospective analysis was conducted using Salmonella Enteritidis isolates (n = 65) from 11 epidemiologically confirmed outbreaks and a collection of contemporaneous sporadic isolates (n = 258) during 2007-2017 to evaluate the performance of WGS in delineating outbreak-associated isolates. Whole-genome single-nucleotide polymorphism (SNP)-based phylogenetic analysis revealed well-supported clades in concordance with epidemiological evidence and pairwise distances of ≤3 SNPs for all outbreaks. WGS-based framework of outbreak detection was thus proposed and applied prospectively to investigate isolates (n = 66) from nine outbreaks during 2018-2019. We further demonstrated the superior discriminatory power and accuracy of WGS to resolve and delineate outbreaks for pragmatic food source tracing. The proposed integrated WGS framework is the first in China for Salmonella Enteritidis and has the potential to serve as a paradigm for outbreak detection and source tracing of Salmonella throughout the stages of food production, as well as expanded to other foodborne pathogens.

Identification, subtyping, and tracking of dairy spoilage-associated Pseudomonas by sequencing the ileS gene

Pseudomonas spp. are important spoilage bacteria that negatively affect the quality of refrigerated fluid milk and uncultured cheese by generating unwanted odors, flavors, and pigments. They are frequently found in dairy plant environments and enter dairy products predominantly as postpasteurization contaminants. Current subtyping and characterization methods for dairy-associated Pseudomonas are often labor-intensive and expensive or provide limited and possibly unreliable classification information (e.g., to the species level). Our goal was to identify a single-copy gene that could be analyzed in dairy spoilage-associated Pseudomonas for preliminary species-level identification, subtyping, and phenotype prediction. We tested 7 genes previously targeted in a Pseudomonas fluorescens multilocus sequence typing scheme for their individual suitability in this application using a set of 113 Pseudomonas spp. isolates representing the diversity of typical pasteurized milk contamination. For each of the 7 candidate genes, we determined the success rate of PCR and sequencing for these 113 isolates as well as the level of discrimination for species identification and subtyping that the sequence data provided. Using these metrics, we selected a single gene, isoleucyl tRNA synthetase (ileS), which had the most suitable traits for simple and affordable single-gene Pseudomonas characterization. This was based on the number of isolates successfully sequenced for ileS (113/113), the number of unique allelic types assigned (83, compared with 50 for 16S rDNA), nucleotide and sequence diversity measures (e.g., number of unique SNP and Simpson index), and tests for genetic recombination. The discriminatory ability of ileS sequencing was confirmed by separation of 99 additional dairy Pseudomonas spp. isolates, which were indistinguishable by 16S rDNA sequencing, into 28 different ileS allelic types. Further, we used whole-genome sequencing data to demonstrate the similarities in ileS-based phylogenetic clustering to whole-genome-based clustering for 27 closely related dairy-associated Pseudomonas spp. isolates and for 178 Pseudomonas type strains. We also found that dairy-associated Pseudomonas within an ileS cluster typically shared the same proteolytic and lipolytic activities. Use of ileS sequencing provides a promising strategy for affordable initial characterization of Pseudomonas isolates, which will help the dairy industry identify, characterize, and track Pseudomonas in their facilities and products.

BlastFrost: fast querying of 100,000s of bacterial genomes in Bifrost graphs

BlastFrost is a highly efficient method for querying 100,000s of genome assemblies, building on Bifrost, a dynamic data structure for compacted and colored de Bruijn graphs. BlastFrost queries a Bifrost data structure for sequences of interest and extracts local subgraphs, enabling the identification of the presence or absence of individual genes or single nucleotide sequence variants. We show two examples using Salmonella genomes: finding within minutes the presence of genes in the SPI-2 pathogenicity island in a collection of 926 genomes and identifying single nucleotide polymorphisms associated with fluoroquinolone resistance in three genes among 190,209 genomes. BlastFrost is available at https://github.com/nluhmann/BlastFrost/tree/master/data .

Models in parasite and pathogen evolution: Genomic analysis reveals predominant clonality and progressive evolution at all evolutionary scales in parasitic protozoa, yeasts and bacteria

The predominant clonal evolution (PCE) model of pathogenic microorganisms postulates that the impact of genetic recombination in those pathogens' natural populations is not enough to erase a persistent phylogenetic signal at all evolutionary scales from microevolution till geological times in the whole ecogeographical range of the species considered. We have tested this model with a set of representative parasitic protozoa, yeasts and bacteria in the light of the most recent genomic data. All surveyed species, including those that were considered as highly recombining, exhibit similar PCE patterns above and under the species level, from macro- to micro-evolutionary scales (Russian doll pattern), suggesting gradual evolution. To our knowledge, it is the first time that such a strong common evolutionary feature among very diverse pathogens has been evidenced. The implications of this model for basic biology and applied research are exposed. These implications include our knowledge on the pathogens' reproductive mode, their population structure, the possibility to type strain and to follow up epidemics (molecular epidemiology) and to revisit pathogens' taxonomy through a flexible use of the phylogenetic species concept (Cracraft, 1983).

Phylogenetic Methods for Genome-Wide Association Studies in Bacteria

Genome-wide association studies in bacteria have great potential to deliver a better understanding of the genetic basis of many biologically important phenotypes, including antibiotic resistance, pathogenicity, and host adaptation. Such studies need however to account for the specificities of bacterial genomics, especially in terms of population structure, homologous recombination, and genomic plasticity. A powerful way to tackle this challenge is to use a phylogenetic approach, which is based on long-standing methodology for the evolutionary analysis of bacterial genomic data. Here we present both the theoretical and practical aspects involved in the use of phylogenetic methods for bacterial genome-wide association studies.

Epidemiological and Bacteriological Investigations Using Whole-Genome Sequencing in a Recurrent Outbreak of Pullorum Disease on a Quail Farm in France

Simple Summary

Although pullorum disease is endemic in many parts of the world, this avian disease responsible for high economic and commercial losses has been eliminated from organized poultry production in Europe and North America. However, it may still remain in backyards and reappear sporadically on conventional poultry farms. This study presents in detail a recurrent outbreak of pullorum disease on a quail farm. In this case report, we present how epidemiological and bacteriological investigations using molecular sequencing tools were carried out in order to identify the source of contamination. Finally, we identify high-risk sanitary practices, and propose recommendations to manage and control this poultry disease, which has become rare in European countries. Given the development of outdoor farms and the increase in self-consumption family farms, a resurgence of pullorum disease cannot be excluded in the coming years. It is essential to develop effective sanitary barriers to prevent transmission between the two coexisting populations of commercial and non-commercial poultry and to raise awareness among all those involved in the poultry industry to be able to detect any outbreak quickly.

Abstract

An outbreak of pullorum disease causing septicemia and high mortality was diagnosed in 2019 on a quail farm in western France. An initial episode had been detected in another building at the same site eight months earlier. Given the exceptional nature and the extent of the potential economic consequences of pullorum disease, epidemiological and bacteriological investigations using molecular sequencing tools were carried out. Salmonella Gallinarum and Salmonella Infantis were isolated (using the NF U 47-101 reference method) from samples taken from birds at the infected site. A resurgence of the initial episode by horizontal transmission of S. Gallinarum is the most likely hypothesis, supported by whole-genome sequencing (WGS) of the strains isolated during the two episodes. Risk health practices have been identified, including the rearing of animals of different ages and species on the same site. Recurrence is explained by the probable persistence of reservoirs of the pathogen on the site (manure, lesser mealworm beetles). The article also highlights the importance of decontamination measures, including pest control, as a key element in the success of the disease control protocol.

Multilevel genome typing: genomics-guided scalable resolution typing of microbial pathogens

Background

Both long- and short-term epidemiology are fundamental to disease control and require accurate bacterial typing. Genomic data resulting from implementation of whole genome sequencing in many public health laboratories can potentially provide highly sensitive and accurate descriptions of strain relatedness. Previous typing efforts using these data have mainly focussed on outbreak detection.

Aim

We aimed to develop multilevel genome typing (MGT), using consecutive multilocus sequence typing (MLST) schemes of increasing sizes, stepping up from seven-gene MLST to core genome MLST, to allow examination of genetic relatedness at multiple resolution levels.

Methods

The system was applied to Salmonellaenterica serovar Typhimurium. The MLST scheme used at each step (MGT level), defined a given MGT-level specific sequence type (ST). The list of STs generated from all of these increasing MGT levels, was named a genome type (GT). Using MGT, we typed 9,096 previously characterised isolates with publicly available data.

Results

Our approach could identify previously described S. Typhimurium populations, such as the DT104 multidrug resistance lineage (GT 19-2-11) and two invasive lineages of African isolates (GT 313-2-3 and 313-2-752). Further, we showed that MGT-derived clusters can accurately distinguish five outbreaks from each other and five background isolates.

Conclusion

MGT provides a universal and stable nomenclature at multiple resolutions for S. Typhimurium strains and could be implemented as an internationally standardised strain identification system. While established so far only for S. Typhimurium, the results here suggest that MGT could form the basis for typing systems in other similar microorganisms.

Genotyping Study of Salmonella 4,[5],12:i:- Monophasic Variant of Serovar Typhimurium and Characterization of the Second-Phase Flagellar Deletion by Whole Genome Sequencing

After Salmonella Enteritidis and S. Typhimurium, S. 4,[5],12:i:- is the most reported serovar in human clinical cases. During the past 20 years, many tools have been used for its typing and second-phase flagellar deletion characterization. Currently, whole genome sequencing (WGS) and different bioinformatic programs have shown the potential to be more accurate than earlier tools. To assess this potential, we analyzed by WGS and in silico typing a selection of 42 isolates of S. 4,[5],12:i:- and S. Typhimurium with different in vitro characteristics. Comparative analysis showed that SeqSero2 does not differentiate fljB-positive S. 4,[5],12:i:- strains from those of serovar Typhimurium. Our results proved that the strains selected for this work were non-clonal S. 4,[5],12:i:- strains circulating in Spain. Using WGS data, we identified 13 different deletion types of the second-phase flagellar genomic region. Most of the deletions were generated by IS26 insertions, showing orientation-dependent conserved deletion ends. In addition, we detected S. 4,[5],12:i:- strains of the American clonal line that would give rise to the Southern European clone in Spain. Our results suggest that new S. 4,[5],12:i:- strains are continuously emerging from different S. Typhimurium strains via different genetic events, at least in swine products.

Carriage and Gene Content Variability of the pESI-Like Plasmid Associated with Salmonella Infantis Recently Established in United States Poultry Production

Salmonella Infantis carrying extended spectrum β-lactamase bla_CTX-M-65 on a pESI-like megaplasmid has recently emerged in United States poultry. In order to determine the carriage rate and gene content variability of this plasmid in U.S. Salmonella Infantis, whole genome sequences of Salmonella isolates from humans and animals in the U.S. and internationally containing the pESI-like plasmid were analyzed. The U.S. Department of Agriculture Food Safety and Inspection Service (FSIS) identified 654 product sampling isolates containing pESI-like plasmids through hazard analysis and critical control point (HACCP) verification testing in 2017 and 2018. The Centers for Disease Control and Prevention identified 55 isolates with pESI-like plasmids in 2016–2018 through the National Antimicrobial Resistance Monitoring System. Approximately 49% of pESI-like plasmids from FSIS verification isolates and 71% from CDC NARMS contained bla_CTX-M-65. Pan-plasmid genome analysis was also performed. All plasmids contained traN and more than 95% contained 172 other conserved genes; 61% contained bla_CTX-M-65. In a hierarchical clustering analysis, some plasmids from U.S. animal sources clustered together and some plasmids from South America clustered together, possibly indicating multiple plasmid lineages. However, most plasmids contained similar genes regardless of origin. Carriage of the pESI-like plasmid in U.S. appears to be limited to Salmonella Infantis and carriage rates increased from 2017 to 2018.

Multilocus Sequence Typing of Staphylococcus aureus

Multilocus sequence typing (MLST) has been successfully used to differentiate and trace the bacterial species and pathogens that cause outbreaks or epidemics of infectious diseases. MLST provides a powerful solution for molecular epidemiological characterization of bacterial strains, including Staphylococcus aureus, by using the sequences of the internal region of seven housekeeping genes. In previous studies, we utilized MLST to analyze the genotypes of S. aureus isolates from pediatric patients with cystic fibrosis and revealed three prevalent ST types, including ST5, ST30, and ST8 in these isolates. In this chapter, we describe a detailed procedure of MLST for genotyping S. aureus.

Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany

Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates-obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)-were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.

Genomic characterization and phylogenetic analysis of Salmonella enterica serovar Javiana

Salmonella enterica serovar Javiana is the fourth most reported serovar of laboratory-confirmed human Salmonella infections in the U.S. and in Tennessee (TN). Although Salmonella ser. Javiana is a common cause of human infection, the majority of cases are sporadic in nature rather than outbreak-associated. To better understand Salmonella ser. Javiana microbial population structure in TN, we completed a phylogenetic analysis of 111 Salmonella ser. Javiana clinical isolates from TN collected from Jan. 2017 to Oct. 2018. We identified mobile genetic elements and genes known to confer antibiotic resistance present in the isolates, and performed a pan-genome-wide association study (pan-GWAS) to compare gene content between clades identified in this study. The population structure of TN Salmonella ser. Javiana clinical isolates consisted of three genetic clades: TN clade I (n = 54), TN clade II (n = 4), and TN clade III (n = 48). Using a 5, 10, and 25 hqSNP distance threshold for cluster identification, nine, 12, and 10 potential epidemiologically-relevant clusters were identified, respectively. The majority of genes that were found to be over-represented in specific clades were located in mobile genetic element (MGE) regions, including genes encoding integrases and phage structures (91.5%). Additionally, a large portion of the over-represented genes from TN clade II (44.9%) were located on an 87.5 kb plasmid containing genes encoding a toxin/antitoxin system (ccdAB). Additionally, we completed phylogenetic analyses of global Salmonella ser. Javiana datasets to gain a broader insight into the population structure of this serovar. We found that the global phylogeny consisted of three major clades (one of which all of the TN isolates belonged to) and two cgMLST eBurstGroups (ceBGs) and that the branch length between the two Salmonella ser. Javiana ceBGs (1,423 allelic differences) was comparable to those from other serovars that have been reported as polyphyletic (929–2,850 allelic differences). This study demonstrates the population structure of TN and global Salmonella ser. Javiana isolates, a clinically important Salmonella serovar and can provide guidance for phylogenetic cluster analyses for public health surveillance and response.

Whole-Genome Sequence Analysis of an Extensively Drug-Resistant Salmonella enterica Serovar Agona Isolate from an Australian Silver Gull (Chroicocephalus novaehollandiae) Reveals the Acquisition of Multidrug Resistance Plasmids

Although most of the approximately 94 million annual human cases of gastroenteritis due to Salmonella enterica resolve without medical intervention, antimicrobial therapy is recommended for patients with severe disease. Wild birds can be natural hosts of Salmonella that pose a threat to human health; however, multiple-drug-resistant serovars of S. enterica have rarely been described. In 2012, silver gull (Chroicocephalus novaehollandiae) chicks at a major breeding colony were shown to host Salmonella, most isolates of which were susceptible to antibiotics. However, multiple-drug-resistant (MDR) Escherichia coli with resistance to carbapenems, ceftazidime, and fluoroquinolones was reported from this breeding colony. In this paper, we describe a novel MDR Salmonella strain subsequently isolated from the same breeding colony. SG17-135, an isolate of S. enterica with phenotypic resistance to 12 individual antibiotics but only nine antibiotic classes including penicillins, cephalosporins, monobactams, macrolides, fluoroquinolones, aminoglycosides, dihydrofolate reductase inhibitors (trimethoprim), sulfonamides, and glycylcyclines was recovered from a gull chick in 2017. Whole-genome sequence (WGS) analysis of SG17-135 identified it as Salmonella enterica serovar Agona (S Agona) with a chromosome comprising 4,813,284 bp, an IncHI2 ST2 plasmid (pSG17-135-HI2) of 311,615 bp, and an IncX1 plasmid (pSG17-135-X) of 27,511 bp. pSG17-135-HI2 housed a complex resistance region comprising 16 antimicrobial resistance genes including bla _CTX-M-55 The acquisition of MDR plasmids by S. enterica described here poses a serious threat to human health. Our study highlights the importance of taking a One Health approach to identify environmental reservoirs of drug-resistant pathogens and MDR plasmids.IMPORTANCE Defining environmental reservoirs hosting mobile genetic elements that shuttle critically important antibiotic resistance genes is key to understanding antimicrobial resistance (AMR) from a One Health perspective. Gulls frequent public amenities, parklands, and sewage and other waste disposal sites and carry drug-resistant Escherichia coli Here, we report on SG17-135, a strain of Salmonella enterica serovar Agona isolated from the cloaca of a silver gull chick nesting on an island in geographic proximity to the greater metropolitan area of Sydney, Australia. SG17-135 is closely related to pathogenic strains of S Agona, displays resistance to nine antimicrobial classes, and carries important virulence gene cargo. Most of the antibiotic resistance genes hosted by SG17-135 are clustered on a large IncHI2 plasmid and are flanked by copies of IS26 Wild birds represent an important link in the evolution and transmission of resistance plasmids, and an understanding of their behavior is needed to expose the interplay between clinical and environmental microbial communities.

Metagenomics of the modern and historical human oral microbiome with phylogenetic studies on Streptococcus mutans and Streptococcus sobrinus

We have recently developed bioinformatic tools to accurately assign metagenomic sequence reads to microbial taxa: SPARSE for probabilistic, taxonomic classification of sequence reads; EToKi for assembling and polishing genomes from short-read sequences; and GrapeTree, a graphic visualizer of genetic distances between large numbers of genomes. Together, these methods support comparative analyses of genomes from ancient skeletons and modern humans. Here, we illustrate these capabilities with 784 samples from historical dental calculus, modern saliva and modern dental plaque. The analyses revealed 1591 microbial species within the oral microbiome. We anticipated that the oral complexes of Socransky et al., which were defined in 1998, would predominate among taxa whose frequencies differed by source. However, although some species discriminated between sources, we could not confirm the existence of the complexes. The results also illustrate further functionality of our pipelines with two species that are associated with dental caries, Streptococcus mutans and Streptococcus sobrinus. They were rare in historical dental calculus but common in modern plaque, and even more common in saliva. Reconstructed draft genomes of these two species from metagenomic samples in which they were abundant were combined with modern public genomes to provide a detailed overview of their core genomic diversity. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Genomic Serotyping, Clinical Manifestations, and Antimicrobial Resistance of Nontyphoidal Salmonella Gastroenteritis in Hospitalized Children in Ho Chi Minh City, Vietnam

Nontyphoidal Salmonella (NTS) are among the most common etiological agents of diarrheal diseases worldwide and have become the most commonly detected bacterial pathogen in children hospitalized with diarrhea in Vietnam. Aiming to better understand the epidemiology, serovar distribution, antimicrobial resistance (AMR), and clinical manifestation of NTS gastroenteritis in Vietnam, we conducted a clinical genomics investigation of NTS isolated from diarrheal children admitted to one of three tertiary hospitals in Ho Chi Minh City.

Nontyphoidal Salmonella (NTS) are among the most common etiological agents of diarrheal diseases worldwide and have become the most commonly detected bacterial pathogen in children hospitalized with diarrhea in Vietnam. Aiming to better understand the epidemiology, serovar distribution, antimicrobial resistance (AMR), and clinical manifestation of NTS gastroenteritis in Vietnam, we conducted a clinical genomics investigation of NTS isolated from diarrheal children admitted to one of three tertiary hospitals in Ho Chi Minh City. Between May 2014 and April 2016, 3,166 children hospitalized with dysentery were recruited into the study; 478 (∼15%) children were found to be infected with NTS by stool culture. Molecular serotyping of the 450 generated genomes identified a diverse collection of serogroups (B, C1, C2 to C3, D1, E1, G, I, K, N, O, and Q); however, Salmonella enterica serovar Typhimurium was the most predominant serovar, accounting for 41.8% (188/450) of NTS isolates. We observed a high prevalence of AMR to first-line treatments recommended by WHO, and more than half (53.8%; 242/450) of NTS isolates were multidrug resistant (MDR; resistant to ≥3 antimicrobial classes). AMR gene detection positively correlated with phenotypic AMR testing, and resistance to empirical antimicrobials was associated with a significantly longer hospitalization (0.91 days; P = 0.04). Our work shows that genome sequencing is a powerful epidemiological tool to characterize the serovar diversity and AMR profiles in NTS. We propose a revaluation of empirical antimicrobials for dysenteric diarrhea and endorse the use of whole-genome sequencing for sustained surveillance of NTS internationally.

Accurate reconstruction of bacterial pan- and core genomes with PEPPAN

Bacterial genomes can contain traces of a complex evolutionary history, including extensive homologous recombination, gene loss, gene duplications, and horizontal gene transfer. To reconstruct the phylogenetic and population history of a set of multiple bacteria, it is necessary to examine their pangenome, the composite of all the genes in the set. Here we introduce PEPPAN, a novel pipeline that can reliably construct pangenomes from thousands of genetically diverse bacterial genomes that represent the diversity of an entire genus. PEPPAN outperforms existing pangenome methods by providing consistent gene and pseudogene annotations extended by similarity-based gene predictions, and identifying and excluding paralogs by combining tree- and synteny-based approaches. The PEPPAN package additionally includes PEPPAN_parser, which implements additional downstream analyses, including the calculation of trees based on accessory gene content or allelic differences between core genes. To test the accuracy of PEPPAN, we implemented SimPan, a novel pipeline for simulating the evolution of bacterial pangenomes. We compared the accuracy and speed of PEPPAN with four state-of-the-art pangenome pipelines using both empirical and simulated data sets. PEPPAN was more accurate and more specific than any of the other pipelines and was almost as fast as any of them. As a case study, we used PEPPAN to construct a pangenome of approximately 40,000 genes from 3052 representative genomes spanning at least 80 species of Streptococcus The resulting gene and allelic trees provide an unprecedented overview of the genomic diversity of the entire Streptococcus genus.

Salmonella identified in pigs in Kenya and Malawi reveals the potential for zoonotic transmission in emerging pork markets

Salmonella is a major cause of foodborne disease globally. Pigs can carry and shed non-typhoidal Salmonella (NTS) asymptomatically, representing a significant reservoir for these pathogens. To investigate Salmonella carriage by African domestic pigs, faecal and mesenteric lymph node samples were taken at slaughter in Nairobi, Busia (Kenya) and Chikwawa (Malawi) between October 2016 and May 2017. Selective culture, antisera testing and whole genome sequencing were performed on samples from 647 pigs; the prevalence of NTS carriage was 12.7% in Busia, 9.1% in Nairobi and 24.6% in Chikwawa. Two isolates of S. Typhimurium ST313 were isolated, but were more closely related to ST313 isolates associated with gastroenteritis in the UK than bloodstream infection in Africa. The discovery of porcine NTS carriage in Kenya and Malawi reveals potential for zoonotic transmission of diarrhoeal strains to humans in these countries, but not for transmission of clades specifically associated with invasive NTS disease in Africa.

Virulence and resistance genes profiles and clonal relationships of non-typhoidal food-borne Salmonella strains isolated in Tunisia by whole genome sequencing

Whole genome sequencing (WGS) has made impressive progress in the field of molecular biology. Its most common application for public health is in the area of surveillance of food-borne diseases. WGS has the potential for providing a large amount of information, such as the identification of the strain type, the characterization of antibiotic resistance and virulence, and phylogeny. In our study, thirty-nine non-typhoidal Salmonella strains were isolated from diverse sources in Tunisia. Non-typhoidal Salmonella are among the most common pathogens contaminating food animals. The presence of virulence and antimicrobial resistance determinants in those strains were investigated using whole genome sequencing (WGS) and appropriate data analysis. The genomes were screened for several Salmonella virulence genes using the Virulence Factor Database VFDB. Twelve different virulence profiles, which correspond to the 12 identified serovars, were recognized. Several antimicrobial resistance genes were also detected: aac (6')-Iaa, sul1, tetA, bla-TEM and qnrS genes. Phylogenetic relationships among the strains were further assessed by a cgMLST analysis. The resulting phylogenetic tree consisted of several clusters consistently with the in silico multilocus sequence typing (MLST) and serotyping. Our findings demonstrated that WGS and subsequent data analysis provided an accurate tool for genetic characterization of bacterial strains compared to usual molecular typing techniques. To the best of our knowledge, this is the first report of an application of WGS for genetic characterization of food-borne Tunisian strains.

Microbe hunting in the modern era: reflecting on a decade of microbial genomic epidemiology

Since the first recognition that infectious microbes serve as the causes of many human diseases, physicians and scientists have sought to understand and control their spread. For the past 150+ years, these ‘microbe hunters’ have learned to combine epidemiological information with knowledge of the infectious agent(s). In this essay, I reflect on the evolution of microbe hunting, beginning with the history of pre-germ theory epidemiological studies, through the microbiological and molecular eras. Now in the genomic age, modern-day microbe hunters are combining pathogen whole-genome sequencing with epidemiological data to enhance epidemiological investigations, advance our understanding of the natural history of pathogens and drivers of disease, and ultimately reshape our plans and priorities for global disease control and eradication. Indeed, as we have seen during the ongoing Covid-19 pandemic, the role of microbe hunters is now more important than ever. Despite the advances already made by microbial genomic epidemiology, the field is still maturing, with many more exciting developments on the horizon.

Analysis of Salmonella enterica serovar Enteritidis isolates from chickens and chicken meat products in Malaysia using PFGE, and MLST

BACKGROUND

Salmonella is a very important foodborne pathogen causing illness in humans. The emergence of drug-resistant strains also constitutes a serious worry to global health and livestock productivity. This study investigated Salmonella isolates from chicken and chicken meat products using the phenotypic antimicrobial screening as well as the molecular characteristics of Salmonella isolates. Upon serotyping of the isolates, the antimicrobial susceptibility profiling using a panel of 9 commonly used antimicrobials was done. Subsequently, the molecular profiles of all the isolates were further determined using Pulsed Field Gel Electrophoresis (PFGE) and the Whole Genome Multi-Locus Sequence Type (wgMLST) analysis in order to obtain the sequence types.

RESULTS

The PFGE data was input into FPQuest software, and the dendrogram generated was studied for possible genetic relatedness among the isolates. All the isolates were found to belong to the Salmonella Enteritidis serotype with notable resistance to tetracycline, gentamycin, streptomycin, and sulfadimidine. The S. Enteritidis isolates tested predominantly subtyped into the ST11 and ST1925, which was found to be a single cell variant of ST11. The STs were found to occur in chicken meats, foods, and live chicken cloacal swabs, which may indicate the persistence of the bacteria in multiple foci.

CONCLUSION

The data demonstrate the presence of S. Enteritidis among chickens, indicating its preference and reservoir status for enteric Salmonella pathogens.

Clinical and molecular characteristics of carbapenem non-susceptible Escherichia coli: A nationwide survey from Oman

The prevalence of carbapenem-resistant Enterobacterales (CRE) in the Arabian Peninsula is predicted to be high, as suggested from published case reports. Of particular concern, is carbapenem-resistant E. coli (CR-EC), due to the importance of this species as a community pathogen. Herein, we conducted a comprehensive molecular characterization of putative CR-EC strains from Oman. We aim to establish a baseline for future molecular monitoring. We performed whole-genome sequencing (WGS) for 35 putative CR-EC. Isolates were obtained from patients at multiple centers in 2015. Genetic relatedness was investigated using several typing approaches such as MLST, SNP calling, phylogroup and CRISPR typing. Maxiuium likelihood SNP-tree was performed by RAxML after variant calling and removal of recombination regions with Snippy and Gubbins, respectively. Resistance genes, plasmid replicon types, virulence genes, and prophage were also characterised. The online databases CGE, CRISPRcasFinder, Phaster and EnteroBase were used for the in silico analyses. Screening for mutations in genes regulating the expression of porins and efflux pump as well as mutations lead to fluoroquinolones resistance were performed with CLC Genomics Workbench. The genetic diversity suggests a polyclonal population structure with 21 sequence types (ST), of which ST38 being the most prevalent (11%). SNPs analysis revealed possible transmission episodes. Whereas, CRISPR typing helped to spot outlier strains belonged to phylogroups other than B2 which was CRISPR-free. The virulent phylogroups B2 and D were detected in 4 and 9 isolates, respectively. In some strains bacteriophages acted as vectors for virulence genes. Regarding resistance to β-lactam, 22 were carbapenemase producers, 3 carbapenem non-susceptible but carbapenemase-negative, 9 resistant to expanded-spectrum cephalosporins, and one isolate with susceptibility to cephalosporins and carbapenems. Thirteen out of the 22 (59%) carbapenemase-producing isolates were NDM and 7 (23%) were OXA-48-like which mirrors the situation in Indian subcontinent. Two isolates co-produced NDM and OXA-48-like enzymes. In total, 80% (28/35) were CTX-M-15 producers and 23% (8/35) featured AmpC. The high-risk subclones ST131-H30Rx/C2, ST410-H24RxC and ST1193-H64RxC were detected, the latter associated with NDM. To our knowledge, this is the first report of ST1193-H64Rx subclone with NDM. In conclusion, strains showed polyclonal population structure with OXA-48 and NDM as the only carbapenemases in CR-EC from Oman. We detected the high-risk subclone ST131-H30Rx/C2, ST410-H24RxC and ST1193-H64RxC. The latter was reported with carbapenemase gene for the first time here.

Comparison of conventional molecular and whole-genome sequencing methods for subtyping Salmonella enterica serovar Enteritidis strains from Tunisia

We sought to determine the relative value of conventional molecular methods and whole-genome sequencing (WGS) for subtyping Salmonella enterica serovar Enteritidis recovered from 2000 to 2015 in Tunisia and to investigate the genetic diversity of this serotype. A total of 175 Salmonella Enteritidis isolates were recovered from human, animal, and foodborne outbreak samples. Pulsed-field gel electrophoresis (PFGE), multiple locus variable-number tandem repeat analysis (MLVA), and whole-genome sequencing were performed. Eight pulsotypes were detected for all isolates with PFGE (DI = 0.518). Forty-five Salmonella Enteritidis isolates were selected for the MLVA and WGS techniques. Eighteen MLVA profiles were identified and classified into two major clusters (DI = 0.889). Core genome multilocus typing (cgMLST) analysis revealed 16 profiles (DI = 0.785). Whole-genome analysis indicated 660 single-nucleotide polymorphism (SNP) divergences dividing these isolates into 43 haplotypes (DI = 0.997). The phylogenetic tree supported the classification of Salmonella Enteritidis isolates into two distinct lineages subdivided into five clades and seven subclades. Pairwise SNP differences between the isolates ranged between 302 and 350. We observed about 311 SNP differences between the two foodborne outbreaks, while only less or equal to 4 SNP differences within each outbreak. SNP-based WGS typing showed an excellent discriminatory power comparing with the conventional methods such as PFGE and MLVA. Besides, we demonstrate the added value of WGS as a complementary subtyping method to discriminate outbreak from non-outbreak isolates belonging to common subtypes. It is important to continue the survey of Salmonella Enteritidis lineages in Tunisia using WGS.

High Prevalence of Human-Associated Escherichia coli in Wetlands Located in Eastern France

Escherichia coli that are present in the rivers are mostly brought by human and animal feces. Contamination occurs mostly through wastewater treatment plant (WWTP) outflows and field amendment with sewage sludge or manure. However, the survival of these isolates in river-associated wetlands remains unknown. Here, we assessed E. coli population structure in low-anthropized wetlands located along three floodplains to identify the major source of contamination of wetlands, whose functioning is different from the rivers. We retrieved 179 E. coli in water samples collected monthly from 19 sites located in eastern France over 1 year. Phylogroups B1 and B2 were dominant in the E. coli population, while phylogroup A was dominant in isolates resistant to third-generation cephalosporins, which harbored the extended-spectrum β-lactamase (ESBL) encoding genes bla_CTX–M–15 and bla_CTX–M–27 in half of the cases. The high proportion of isolates from human source can be attributed to WWTP outflows and the spread of sewage sludge. We analyzed the distribution of the isolates belonging to the most human-associated phylogroups (B2 and D) on a phylogenetic tree of the whole species and compared it with that of isolates retrieved from patients and from WWTP outflows. The distribution of the three E. coli populations was similar, suggesting the absence of a specific population in the environment. Our results suggest that a high proportion of E. coli isolates that reach and survive in low-anthropized environments such as wetlands are from human source. To the best of our knowledge, this is the first study assessing E. coli contamination and resistance genes in natural freshwater wetlands.

A database for risk assessment and comparative genomic analysis of foodborne Vibrio parahaemolyticus in China

Vibrio parahaemolyticus is a major foodborne pathogen worldwide. The increasing number of cases of V. parahaemolyticus infections in China indicates an urgent need to evaluate the prevalence and genetic diversity of this pathogenic bacterium. In this paper, we introduce the Foodborne Vibrio parahaemolyticus genome database (FVPGD), the first scientific database of foodborne V. parahaemolyticus distribution and genomic data in China, based on our previous investigations of V. parahaemolyticus contamination in different kinds of food samples across China from 2011 to 2016. The dataset includes records of 2,499 food samples and 643 V. parahaemolyticus strains from supermarkets and marketplaces distributed over 39 cities in China; 268 whole-genome sequences have been deposited in this database. A spatial view on the risk situations of V. parahaemolyticus contamination in different food types is provided. Additionally, the database provides a functional interface of sequence BLAST, core genome multilocus sequence typing, and phylogenetic analysis. The database will become a powerful tool for risk assessment and outbreak investigations of foodborne pathogens in China.

MLST based serotype prediction for the accurate identification of non typhoidal Salmonella serovars

Traditional serotyping based on the phenotypic variation of O- and H-antigen remains as the gold-standard for the identification and classification of Salmonella isolates for last 70 years. Although this classification is a globally recognized nomenclature, huge diversity of Salmonella serotypes have made the serovar identification to be very complex. Seven gene multilocus sequence typing (MLST) on the other hand can provide serovar prediction as well as the evolutionary origin between the serovars. In this study non typhoidal Salmonella (NTS) strains (n = 45) isolated from clinical samples (blood, faeces and pus) were identified by traditional phenotypic serotyping and biochemical testing. All the tested Salmonella isolates were designated as serovar Typhimurium based on phenotyping. However, by MLST 60% (27/45) of the isolates were S. Typhimurium, 35.5% (16/45) were S. Agona (ST13), 2.2% (1/45) were S. Kentucky (ST198) and 2.2% (1/45) were S. Saintpaul (ST27). MLST analysis assigned S. Typhimurium isolates as ST36 (18/127), ST19 (7/27) and ST313 (2/27). Mismatches in serovar designation between MLST database and phenotypic serotyping can be due to the misinterpretation of phenotypic serotyping as the antigenic structures of S. Typhimurium, S. Agona differs by a surface antigen. MLST based phylogeny of study isolates showed clustering according to sequence types. Concordance between MLST based sequence type and phenotypic serotype is important to provide insights into genetic population structure of Salmonella.

Correction: Multilocus Sequence Typing as a Replacement for Serotyping in Salmonella enterica

[This corrects the article DOI: 10.1371/journal.ppat.1002776.].

Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

Background: Most publicly available genomes of Salmonella enterica are from human disease in the US and the UK, or from domesticated animals in the US. Methods: Here we describe a historical collection of 10,000 strains isolated between 1891-2010 in 73 different countries. They encompass a broad range of sources, ranging from rivers through reptiles to the diversity of all S. enterica isolated on the island of Ireland between 2000 and 2005. Genomic DNA was isolated, and sequenced by Illumina short read sequencing. Results: The short reads are publicly available in the Short Reads Archive. They were also uploaded to EnteroBase, which assembled and annotated draft genomes. 9769 draft genomes which passed quality control were genotyped with multiple levels of multilocus sequence typing, and used to predict serovars. Genomes were assigned to hierarchical clusters on the basis of numbers of pair-wise allelic differences in core genes, which were mapped to genetic Lineages within phylogenetic trees. Conclusions: The University of Warwick/University College Cork (UoWUCC) project greatly extends the geographic sources, dates and core genomic diversity of publicly available S. enterica genomes. We illustrate these features by an overview of core genomic Lineages within 33,000 publicly available Salmonella genomes whose strains were isolated before 2011. We also present detailed examinations of HC400, HC900 and HC2000 hierarchical clusters within exemplar Lineages, including serovars Typhimurium, Enteritidis and Mbandaka. These analyses confirm the polyphyletic nature of multiple serovars while showing that discrete clusters with geographical specificity can be reliably recognized by hierarchical clustering approaches. The results also demonstrate that the genomes sequenced here provide an important counterbalance to the sampling bias which is so dominant in current genomic sequencing.

Genome-wide Identification and Characterization of a Superfamily of Bacterial Extracellular Contractile Injection Systems

Several phage-tail-like nanomachines were shown to play an important role in the interactions between bacteria and their eukaryotic hosts. These apparatuses appear to represent a new injection paradigm. Here, with three verified extracellular contractile injection systems (eCISs), a protein profile and genomic context-based iterative approach was applied to identify 631 eCIS-like loci from the 11,699 publicly available complete bacterial genomes. The eCIS superfamily, which is phylogenetically diverse and sub-divided into six families, is distributed among Gram-negative and -positive bacteria in addition to archaea. Our results show that very few bacteria are seen to possess intact operons of both eCIS and type VI secretion systems (T6SSs). An open access online database of all detected eCIS-like loci is presented to facilitate future studies. The presence of this bacterial injection machine in a multitude of organisms suggests that it may play an important ecological role in the life cycles of many bacteria.

•

eCIS loci are widely distributed among bacteria genomes

•

eCIS loci encode phage-tail-like proteinaceous machines

•

eCIS superfamily is grouped into six families with distinct genetic features

Twentieth-century emergence of antimicrobial resistant human- and bovine-associated Salmonella enterica serotype Typhimurium lineages in New York State

Salmonella enterica serotype Typhimurium (S. Typhimurium) boasts a broad host range and can be transmitted between livestock and humans. While members of this serotype can acquire resistance to antimicrobials, the temporal dynamics of this acquisition is not well understood. Using New York State (NYS) and its dairy cattle farms as a model system, 87 S. Typhimurium strains isolated from 1999 to 2016 from either human clinical or bovine-associated sources in NYS were characterized using whole-genome sequencing. More than 91% of isolates were classified into one of four major lineages, two of which were largely susceptible to antimicrobials but showed sporadic antimicrobial resistance (AMR) gene acquisition, and two that were largely multidrug-resistant (MDR). All four lineages clustered by presence and absence of elements in the pan-genome. The two MDR lineages, one of which resembled S. Typhimurium DT104, were predicted to have emerged circa 1960 and 1972. The two largely susceptible lineages emerged earlier, but showcased sporadic AMR determinant acquisition largely after 1960, including acquisition of cephalosporin resistance-conferring genes after 1985. These results confine the majority of AMR acquisition events in NYS S. Typhimurium to the twentieth century, largely within the era of antibiotic usage.

Application of Whole-Genome Sequences and Machine Learning in Source Attribution of Salmonella Typhimurium

Prevention of the emergence and spread of foodborne diseases is an important prerequisite for the improvement of public health. Source attribution models link sporadic human cases of a specific illness to food sources and animal reservoirs. With the next generation sequencing technology, it is possible to develop novel source attribution models. We investigated the potential of machine learning to predict the animal reservoir from which a bacterial strain isolated from a human salmonellosis case originated based on whole-genome sequencing. Machine learning methods recognize patterns in large and complex data sets and use this knowledge to build models. The model learns patterns associated with genetic variations in bacteria isolated from the different animal reservoirs. We selected different machine learning algorithms to predict sources of human salmonellosis cases and trained the model with Danish Salmonella Typhimurium isolates sampled from broilers (n = 34), cattle (n = 2), ducks (n = 11), layers (n = 4), and pigs (n = 159). Using cgMLST as input features, the model yielded an average accuracy of 0.783 (95% CI: 0.77-0.80) in the source prediction for the random forest and 0.933 (95% CI: 0.92-0.94) for the logit boost algorithm. Logit boost algorithm was most accurate (valid accuracy: 92%, CI: 0.8706-0.9579) and predicted the origin of 81% of the domestic sporadic human salmonellosis cases. The most important source was Danish produced pigs (53%) followed by imported pigs (16%), imported broilers (6%), imported ducks (2%), Danish produced layers (2%), Danish produced cattle and imported cattle (<1%) while 18% was not predicted. Machine learning has potential for improving source attribution modeling based on sequence data. Results of such models can inform risk managers to identify and prioritize food safety interventions.

Low Dose Infection of Hens in Lay with Salmonella enterica Serovar Enteritidis from Different Genomic Clades

Salmonella enterica serovar Enteritidis is the leading cause of salmonellosis in people, and modeling of infections in chickens is used to identify intervention strategies. A review of 80 manuscripts encompassing 119 experiments indicated that the mean dose of infection was 10⁸ CFU per bird. Experiments of less than 10⁶ CFU were primarily conducted in immature birds. To address a lack of information on the impact of low dosages on the hen at lay, two experiments were conducted in triplicate. Experiment A addressed issues associated with vaccination; thus, hens were infected intramuscularly at 10³, 10⁵, and 10⁷ CFU. For Experiment B, which was focused more on colonization and invasion, hens were infected orally with 5 × 10³ CFU with 4 strains from different genomic clades. Samples from liver, spleen, ovarian pedicle, and paired ceca in both experiments were cultured 5, 6, 7, and 8 days postinfection. Eggshell microbiome taxa were assessed in Experiment B. Results indicated that dosages of 10³ CFU in both experiments produced enough positive samples to be used within models. The intramuscular route resulted in approximately twice as many positive samples as the oral route. The kinetics of infection appeared to differ between low and high dosages suggestive of a J-curve response. These results could impact risk assessments if the hen at lay has a nonlinear response to infectious dose.

Atypical Non-H2S-Producing Monophasic Salmonella Typhimurium ST3478 Strains from Chicken Meat at Processing Stage Are Adapted to Diverse Stresses

Poultry products are still an important cause of Salmonella infections worldwide, with an increasingly reported expansion of less-frequent serotypes or atypical strains that are frequently multidrug-resistant. Nevertheless, the ability of Salmonella to survive antimicrobials promoted in the context of antibiotic reducing/replacing and farming rethinking (e.g., organic acids and copper in feed/biocides) has been scarcely explored. We investigated Salmonella occurrence (conventional and molecular assays) among chicken meat at the processing stage (n = 53 batches/29 farms) and characterized their tolerance to diverse stress factors (antibiotics, copper, acid pH, and peracetic acid). Whole-genome sequencing was used to assess adaptive features and to perform comparative analysis. We found a low Salmonella occurrence (4%) and identified S. Enteritidis/ST11 plus atypical non-H2S-producing S. 1,4,[5],12:i:-/ST3478. The ST3478 presented the ability to grow under diverse stresses (antibiotics, copper, and acid-pH). Comparative genomics among ST3478 isolates showed similar antibiotic/metal resistance gene repertoires and identical nonsense phsA thiosulfate reductase mutations (related to H2S-negative phenotype), besides their close phylogenetic relationship by cgMLST and SNPs. This study alerts for the ongoing national and international spread of an emerging monophasic Salmonella Typhimurium clonal lineage with an enlarged ability to survive to antimicrobials/biocides commonly used in poultry production, being unnoticed by conventional Salmonella detection approaches due to an atypical non-H2S-producing phenotype.

Municipal Wastewater Surveillance Revealed a High Community Disease Burden of a Rarely Reported and Possibly Subclinical Salmonella enterica Serovar Derby Strain

Clinical surveillance of enteric pathogens like Salmonella is integral to track outbreaks and endemic disease trends. However, clinic-centered disease monitoring biases toward detection of severe cases and underestimates the incidence of self-limiting gastroenteritis and asymptomatic strains. Monitoring pathogen loads and diversity in municipal wastewater (MW) can provide insight into asymptomatic or subclinical infections which are not reflected in clinical cases. Subclinical infection patterns may explain the unusual observation from a year-long sampling campaign in Hawaii: Salmonella enterica serovar Derby was the most abundant pulsotype in MW but was detected infrequently in clinics over the sampling period. Using whole-genome sequencing data of Salmonella isolates from MW and public databases, we demonstrate that the Derby serovar has lower virulence potential than other clinical serovars, particularly based on its reduced profile of genes linked with immune evasion and symptom production, suggesting its potential as a subclinical salmonellosis agent. Furthermore, MW had high abundance of a rare Derby sequence type (ST), ST-72 (rather than the more common ST-40). ST-72 isolates had higher frequencies of fimbrial adherence genes than ST-40 isolates; these are key virulence factors involved in colonization and persistence of infections. However, ST-72 isolates lack the Derby-specific Salmonella pathogenicity island 23 (SPI-23), which invokes host immune responses. In combination, ST-72's genetic features may lead to appreciable infection rates without obvious symptom production, allowing for subclinical persistence in the community. This study demonstrated wastewater's capability to provide community infectious disease information-such as background infection rates of subclinical enteric illness-which is otherwise inaccessible through clinical approaches.IMPORTANCE Wastewater-based epidemiology (WBE) has been conventionally used to analyze community health via the detection of chemicals, such as legal and illicit drugs; however, municipal wastewater contains microbiological determinants of health and disease as well, including enteric pathogens. Here, we demonstrate that WBE can be used to examine subclinical community salmonellosis patterns. Derby was the most abundant Salmonella serovar detected in Hawaii wastewater over a year-long sampling study, with few corresponding clinical cases. Comparative genomics analyses indicate that the normally rare strain of S Derby found in wastewater has a unique combination of genes which allow it to persist as a subclinical infection without producing symptoms of severe gastroenteritis. This study shows that WBE can be used to explore trends in community infectious disease patterns which may not be reflected in clinical monitoring, shedding light on overall enteric disease burden and rates of asymptomatic cases.

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.

Whole-genome sequencing analysis and CRISPR genotyping of rare antibiotic-resistant Salmonella enterica serovars isolated from food and related sources

For decades, Salmonella Typhimurium and Salmonella Enteritidis have prevailed in several countries as agents of salmonellosis outbreaks. In Brazil, the largest exporter of poultry meat, relatively little attention has been paid to infrequent serovars. Here, we report the emergence and characterization of rare serovars isolated from food and related sources collected between 2014 and 2016 in Brazil. Twenty-two Salmonella enterica isolates were analyzed through the use of whole-genome sequencing (WGS) and clustered regularly interspaced short palindromic repeats (CRISPR) genotyping. These isolates were classified into 10 infrequent serovars, including S. Abony, S. Isangi, S. Rochdale, S. Saphra, S. Orion, S. Ouakam, S. Grumpensis, S. Carrau, S. Abaetetuba, and S. Idikan. The presence of six antimicrobial resistance (AMR) genes, qnrB19, bla_CMY-2, tetA, aac(6')-Iaa, sul2 and fosA7, which encode resistance to quinolones, third-generation cephalosporin, tetracycline, aminoglycoside, sulfonamide and fosfomycin, respectively, were confirmed by WGS. All S. Isangi harbored qnrB19 with conserved genomic context across strains, while S. Abony harbored bla_CMY-2. Twelve (54.5%) strains displayed chromosomal mutations in parC (Thr57→Ser). Most serovars were classified as independent lineages, except S. Abony and S. Abaetetuba, which phylogenetically nested with Salmonella strains from different countries. CRISPR analysis revealed that the spacer content was strongly correlated with serovar and multi-locus sequence type for all strains, independently confirming the observed phylogenetic patterns, and highlighting the value of CRISPR-based genotyping for Salmonella. These findings add valuable information to the epidemiology of S. enterica in Brazil, where the emergency of antibiotic-resistant Salmonella continues to evolve.

A publicly accessible database for Clostridioides difficile genome sequences supports tracing of transmission chains and epidemics

Clostridioides difficile is the primary infectious cause of antibiotic-associated diarrhea. Local transmissions and international outbreaks of this pathogen have been previously elucidated by bacterial whole-genome sequencing, but comparative genomic analyses at the global scale were hampered by the lack of specific bioinformatic tools. Here we introduce a publicly accessible database within EnteroBase (http://enterobase.warwick.ac.uk) that automatically retrieves and assembles C. difficile short-reads from the public domain, and calls alleles for core-genome multilocus sequence typing (cgMLST). We demonstrate that comparable levels of resolution and precision are attained by EnteroBase cgMLST and single-nucleotide polymorphism analysis. EnteroBase currently contains 18 254 quality-controlled C. difficile genomes, which have been assigned to hierarchical sets of single-linkage clusters by cgMLST distances. This hierarchical clustering is used to identify and name populations of C. difficile at all epidemiological levels, from recent transmission chains through to epidemic and endemic strains. Moreover, it puts newly collected isolates into phylogenetic and epidemiological context by identifying related strains among all previously published genome data. For example, HC2 clusters (i.e. chains of genomes with pairwise distances of up to two cgMLST alleles) were statistically associated with specific hospitals (P<10-4) or single wards (P=0.01) within hospitals, indicating they represented local transmission clusters. We also detected several HC2 clusters spanning more than one hospital that by retrospective epidemiological analysis were confirmed to be associated with inter-hospital patient transfers. In contrast, clustering at level HC150 correlated with k-mer-based classification and was largely compatible with PCR ribotyping, thus enabling comparisons to earlier surveillance data. EnteroBase enables contextual interpretation of a growing collection of assembled, quality-controlled C. difficile genome sequences and their associated metadata. Hierarchical clustering rapidly identifies database entries that are related at multiple levels of genetic distance, facilitating communication among researchers, clinicians and public-health officials who are combatting disease caused by C. difficile.

Whole-Genome-Based Survey for Polyphyletic Serovars of Salmonella enterica subsp. enterica Provides New Insights into Public Health Surveillance

Serotyping has traditionally been considered the basis for surveillance of Salmonella, but it cannot distinguish distinct lineages sharing the same serovar that vary in host range, pathogenicity and epidemiology. However, polyphyletic serovars have not been extensively investigated. Public health microbiology is currently being transformed by whole-genome sequencing (WGS) data, which promote the lineage determination using a more powerful and accurate technique than serotyping. The focus in this study is to survey and analyze putative polyphyletic serovars. The multi-locus sequence typing (MLST) phylogenetic analysis identified four putative polyphyletic serovars, namely, Montevideo, Bareilly, Saintpaul, and Muenchen. Whole-genome-based phylogeny and population structure highlighted the polyphyletic nature of Bareilly and Saintpaul and the multi-lineage nature of Montevideo and Muenchen. The population of these serovars was defined by extensive genetic diversity, the open pan genome and the small core genome. Source niche metadata revealed putative existence of lineage-specific niche adaptation (host-preference and environmental-preference), exhibited by lineage-specific genomic contents associated with metabolism and transport. Meanwhile, differences in genetic profiles relating to virulence and antimicrobial resistance within each lineage may contribute to pathogenicity and epidemiology. The results also showed that recombination events occurring at the H1-antigen loci may be an important reason for polyphyly. The results presented here provide the genomic basis of simple, rapid, and accurate identification of phylogenetic lineages of these serovars, which could have important implications for public health.

Extended-Spectrum-β-Lactamase- and Plasmid AmpC-Producing Escherichia coli Causing Community-Onset Bloodstream Infection: Association of Bacterial Clones and Virulence Genes with Septic Shock, Source of Infection, and Recurrence

Invasive infections due to extended-spectrum-β-lactamase- and pAmpC-producing Escherichia coli (ESBL/pAmpC-EC) are an important cause of morbidity, often caused by the high-risk clone sequence type (ST131) and isolates classified as extraintestinal pathogenic E. coli (ExPEC). The relative influence of host immunocompetence versus microbiological virulence factors in the acquisition and outcome of bloodstream infections (BSI) is poorly understood. Herein, we used whole-genome sequencing on 278 blood culture isolates of ESBL/pAmpC-EC from 260 patients with community-onset BSI collected from 2012 to 2015 in Stockholm to study the association of virulence genes, sequence types, and antimicrobial resistance with severity of disease, infection source, ESBL/pAmpC-EC BSI low-risk patients, and patients with repeated episodes. ST131 subclade C2 comprised 29% of all patients. Factors associated with septic shock in multivariable analysis were patient host factors (hematologic cancer or transplantation and reduced daily living activity), presence of the E. coli virulence factor iss (increased serum survival), absence of phenotypic multidrug resistance, and absence of the genes pap and hsp Adhesins, particularly pap, were associated with urinary tract infection (UTI) source, while isolates from post-prostate biopsy sepsis had a low overall number of virulence operons, including adhesins, and commonly belonged to ST131 clades A, B, and subclade C1, ST1193, and ST648. ST131 was associated with recurrent episodes. In conclusion, the most interesting finding is the association of iss with septic shock. Adhesins are important for UTI pathogenesis, while otherwise low-pathogenic isolates from the microbiota can cause post-prostate biopsy sepsis.

Emergence of Multidrug-Resistant Salmonella enterica Subspecies enterica Serovar Infantis of Multilocus Sequence Type 2283 in German Broiler Farms

During the last decade, Salmonella enterica subspecies enterica serovar Infantis (S. Infantis) has become more prevalent across Europe with an increased capability to persist in broiler farms. In this study, we aimed to identify potential genetic causes for the increased emergence and longer persistence of S. Infantis in German poultry farms by high-throughput-sequencing. Broiler derived S. Infantis strains from two decades, the 1990s (n = 12) and the 2010s (n = 18), were examined phenotypically and genotypically to detect potential differences responsible for increased prevalence and persistence. S. Infantis organisms were characterized by serotyping and determining antimicrobial susceptibility using the microdilution method. Genotypic characteristics were analyzed by whole genome sequencing (WGS) to detect antimicrobial resistance and virulence genes as well as plasmids. To detect possible clonal relatedness within S. Infantis organisms, 17 accessible genomes from previous studies about emergent S. Infantis were downloaded and analyzed using complete genome sequence of SI119944 from Israel as reference. In contrast to the broiler derived antibiotic-sensitive S. Infantis strains from the 1990s, the majority of strains from the 2010s (15 out of 18) revealed a multidrug-resistance (MDR) phenotype that encodes for at least three antimicrobials families: aminoglycosides [ant(3“)-Ia], sulfonamides (sul1), and tetracyclines [tet(A)]. Moreover, these MDR strains carry a virulence gene pattern missing in strains from the 1990s. It includes genes encoding for fimbriae clusters, the yersiniabactin siderophore, mercury and disinfectants resistance and toxin/antitoxin complexes. In depth genomic analysis confirmed that the 15 MDR strains from the 2010s carry a pESI-like megaplasmid with resistance and virulence gene patterns detected in the emerged S. Infantis strain SI119944 from Israel and clones inside and outside Europe. Genotyping analysis revealed two sequence types (STs) among the resistant strains from the 2010s, ST2283 (n = 13) and ST32 (n = 2). The sensitive strains from the 1990s, belong to sequence type ST32 (n = 10) and ST1032 (n = 2). Therefore, this study confirms the emergence of a MDR S. Infantis pESI-like clone of ST2283 in German broiler farms with presumably high tendency of dissemination. Further studies on the epidemiology and control of S. Infantis in broilers are needed to prevent the transfer from poultry into the human food chain.

Corynebacterium silvaticum sp. nov., a unique group of NTTB corynebacteria in wild boar and roe deer

A total of 34 Corynebacterium sp. strains were isolated from caseous lymph node abscesses of wild boar and roe deer in different regions of Germany. They showed slow growth on Columbia sheep blood agar and sparse growth on Hoyle's tellurite agar. Cellular fatty acid analysis allocated them in the C. diphtheriae group of genus Corynebacterium. MALDI-TOF MS using specific database extensions and rpoB sequencing resulted in classification as C. ulcerans. Their quinone system is similar to C. ulcerans, with major menaquinone MK-8(H2). Their complex polar lipid profile includes major lipids phosphatidylinositol, phosphatidylinositol-mannoside, diphosphatidylglycerol, but also unidentified glycolipids, distinguishing them clearly from C. ulcerans. They ferment glucose, ribose and maltose (like C. ulcerans), but do not utilise d-xylose, mannitol, lactose, sucrose and glycogen (like C. pseudotuberculosis). They showed activity of catalase, urease and phospholipase D, but variable results for alkaline phosphatase and alpha-glucosidase. All were non-toxigenic, tox gene bearing and susceptible to clindamycin, penicillin and erythromycin. In 16SrRNA gene and RpoB protein phylogenies the strains formed distinct brancheswith C. ulcerans as nearest relative.Whole genome sequencing revealed the unique sequence type 578, a distinctbranch in pangenomic core genome MLST, average nucleotide identities <91%, enhancedgenome sizes (2.55 Mbp) and G/C content (54.4 mol%) compared to related species.These results suggest that the strains represent a novel species, for which wepropose the name Corynebactriumsilvaticum sp. nov., based on their first isolation from forest-dwellinggame animals. The type strain isKL0182^T (= CVUAS 4292^T = DSM 109166^T = LMG 31313^T= CIP 111 672^T).

The phenotypic and molecular characteristics of antimicrobial resistance of Salmonella enterica subsp. enterica serovar Typhimurium in Henan Province, China

Background

Salmonella enterica subsp. enterica serovar Typhimurium infections continue to be a significant public health threat worldwide. The aim of this study was to investigate antibiotic resistance among 147 S. Typhimurium isolates collected from patients in Henan, China from 2006 to 2015.

Methods

147 S. Typhimurium isolates were collected from March 2006 to November 2015 in Henan Province, China. Antimicrobial susceptibility testing was performed, and the resistant genes of ciprofloxacin, cephalosporins (ceftriaxone and cefoxitin) and azithromycin were detected and sequenced. Clonal relationships were assessed by multilocus sequence typing (MLST) and pulsed-field gel electrophoresis (PFGE).

Results

Of the 147 isolates, 91.1% were multidrug resistant (MDR), with 4.1% being resistant to all antibiotic classes tested. Of concern, 13 MDR isolates were co-resistant to the first-line treatments cephalosporins and ciprofloxacin, while three were also resistant to azithromycin. Seven PFGE patterns were identified among the 13 isolates. All of the isolates could be assigned to one of four main groups, with a similarity value of 89%. MLST assigned the 147 isolates into five STs, including two dominant STs (ST19 and ST34). Of the 43 ciprofloxacin-resistant isolates, 39 carried double gyrA mutations (Ser83Phe, Asp87Asn/Tyr/Gly) and a single parC (Ser80Arg) mutation, including 1 isolate with four mutations (gyrA: Ser83Phe, Asp87Gly; parC: Ser80Arg; parE: Ser458Pro). In addition, 12 isolates not only carried mutations in gyrA and parC but also had at least one plasmid-mediated quinolone resistance (PMQR) gene. Among the 32 cephalosporin-resistant isolates, the most common extended-spectrum β-lactamase (ESBL) gene was bla_OXA-1, followed by bla_CTX-M, bla_TEM-1, and bla_CMY-2. Moreover, the mphA gene was identified in 5 of the 15 azithromycin-resistant isolates. Four MDR isolates contained ESBL and PMQR genes, and one of them also carried mphA in addition.

Conclusion

The high level of antibiotic resistance observed in S. Typhimurium poses a great danger to public health, so continuous surveillance of changes in antibiotic resistance is necessary.

Outbreak of Salmonella enterica serotype Poona in infants linked to persistent Salmonella contamination in an infant formula manufacturing facility, France, August 2018 to February 2019

We describe a Salmonella Poona outbreak involving 31 infant cases in France. Following outbreak detection on 18 January 2019, consumption of rice-based infant formula manufactured at a facility in Spain was identified as the probable cause, leading to a recall on 24 January. Whole genome sequencing analysis linked present outbreak isolates to a 2010–11 S. Poona outbreak in Spain associated with formula manufactured in the same facility, indicating a persistent source of contamination.

Genetic diversity and delineation of Salmonella Agona outbreak strains by next generation sequencing, Bavaria, Germany, 1993 to 2018

Background

In 2017, a food-borne Salmonella Agona outbreak caused by infant milk products from a French supplier occurred in Europe. Simultaneously, S. Agona was detected in animal feed samples in Bavaria.

Aim

Using next generation sequencing (NGS) and three data analysis methods, this study’s objectives were to verify clonality of the Bavarian feed strains, rule out their connection to the outbreak, explore the genetic diversity of Bavarian S. Agona isolates from 1993 to 2018 and compare the analysis approaches employed, for practicality and ability to delineate outbreaks caused by the genetically monomorphic Agona serovar.

Methods

In this observational retrospective study, three 2017 Bavarian feed isolates were compared to a French outbreak isolate and 48 S. Agona isolates from our strain collections. The later included human, food, feed, veterinary and environmental isolates, of which 28 were epidemiologically outbreak related. All isolates were subjected to NGS and analysed by: (i) a publicly available species-specific core genome multilocus sequence typing (cgMLST) scheme, (ii) single nucleotide polymorphism phylogeny and (iii) an in-house serovar-specific cgMLST scheme. Using additional international S. Agona outbreak NGS data, the cluster resolution capacity of the two cgMLST schemes was assessed.

Results

We could prove clonality of the feed isolates and exclude their relation to the French outbreak. All approaches confirmed former Bavarian epidemiological clusters.

Conclusion

Even for S. Agona, species-level cgMLST can produce reasonable resolution, being standardisable by public health laboratories. For single samples or homogeneous sample sets, higher resolution by serovar-specific cgMLST or SNP genotyping can facilitate outbreak investigations.

High-resolution characterisation of ESBL/pAmpC-producing Escherichia coli isolated from the broiler production pyramid

The presence of extended-spectrum β-lactamase (ESBL) or plasmid-mediated AmpC β-lactamase (pAmpC)-producing Escherichia coli (ESBL/pAmpC-EC) in livestock is a public health risk given the likelihood of their transmission to humans via the food chain. We conducted whole genome sequencing on 100 ESBL/pAmpC-EC isolated from the broiler production to explore their resistance and virulence gene repertoire, characterise their plasmids and identify transmission events derived from their phylogeny. Sequenced isolates carried resistance genes to four antimicrobial classes in addition to cephalosporins. Virulence gene analysis assigned the majority of ESBL/pAmpC-EC to defined pathotypes. In the complex genetic background of ESBL/pAmpC-EC, clusters of closely related isolates from various production stages were identified and indicated clonal transmission. Phylogenetic comparison with publicly available genomes suggested that previously uncommon ESBL/pAmpC-EC lineages could emerge in poultry, while others might contribute to the maintenance and dissemination of ESBL/pAmpC genes in broilers. The majority of isolates from diverse E. coli lineages shared four dominant plasmids (IncK2, IncI1, IncX3 and IncFIB/FII) with identical ESBL/pAmpC gene insertion sites. These plasmids have been previously reported in diverse hosts, including humans. Our findings underline the importance of specific plasmid groups in the dissemination of cephalosporin resistance genes within the broiler industry and across different reservoirs.

Easy phylotyping of Escherichia coli via the EzClermont web app and command-line tool

The Clermont PCR method for phylotyping Escherichia coli remains a useful classification scheme even though genome sequencing is now routine, and higher-resolution sequence typing schemes are now available. Relating present-day whole-genome E. coli classifications to legacy phylotyping is essential for harmonizing the historical literature and understanding of this important organism. Therefore, we present EzClermont – a novel in silico Clermont PCR phylotyping tool to enable ready application of this phylotyping scheme to whole-genome assemblies. We evaluate this tool against phylogenomic classifications, and an alternative software implementation of Clermont typing. EzClermont is available as a web app at www.ezclermont.org, and as a command-line tool at https://nickp60.github.io/EzClermont/.

Network Approach to Source Attribution of Salmonella enterica Serovar Typhimurium and Its Monophasic Variant

Salmonella enterica subspecies enterica serovar Typhimurium and its monophasic variant are among the most common Salmonella serovars associated with human salmonellosis each year. Related infections are often due to consumption of contaminated meat of pig, cattle, and poultry origin. In order to evaluate novel microbial subtyping methods for source attribution, an approach based on weighted networks was applied on 141 human and 210 food and animal isolates of pigs, broilers, layers, ducks, and cattle collected in Denmark from 2013 to 2014. A whole-genome SNP calling was performed along with cgMLST and wgMLST. Based on these genomic input data, pairwise distance matrices were built and used as input for construction of a weighted network where nodes represent genomes and links to distances. Analyzing food and animal Typhimurium genomes, the coherence of source clustering ranged from 89 to 90% for animal source, from 84 to 85% for country, and from 63 to 65% for year of isolation and was equal to 82% for serotype, suggesting animal source as the first driver of clustering formation. Adding human isolate genomes to the network, a percentage between 93.6 and 97.2% clustered with the existing component and only a percentage between 2.8 and 6.4% appeared as not attributed to any animal sources. The majority of human genomes were attributed to pigs with probabilities ranging from 83.9 to 84.5%, followed by broilers, ducks, cattle, and layers in descending order. In conclusion, a weighted network approach based on pairwise SNPs, cgMLST, and wgMLST matrices showed promising results for source attribution studies.

Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation

Salmonella enterica serotype Typhimurium (S. Typhimurium) is a leading cause of gastroenteritis and bacteraemia worldwide, and a model organism for the study of host-pathogen interactions. Two S. Typhimurium strains (SL1344 and ATCC14028) are widely used to study host-pathogen interactions, yet genotypic variation results in strains with diverse host range, pathogenicity and risk to food safety. The population structure of diverse strains of S. Typhimurium revealed a major phylogroup of predominantly sequence type 19 (ST19) and a minor phylogroup of ST36. The major phylogroup had a population structure with two high order clades (α and β) and multiple subclades on extended internal branches, that exhibited distinct signatures of host adaptation and anthropogenic selection. Clade α contained a number of subclades composed of strains from well characterized epidemics in domesticated animals, while clade β contained multiple subclades associated with wild avian species. The contrasting epidemiology of strains in clade α and β was reflected by the distinct distribution of antimicrobial resistance (AMR) genes, accumulation of hypothetically disrupted coding sequences (HDCS), and signatures of functional diversification. These observations were consistent with elevated anthropogenic selection of clade α lineages from adaptation to circulation in populations of domesticated livestock, and the predisposition of clade β lineages to undergo adaptation to an invasive lifestyle by a process of convergent evolution with of host adapted Salmonella serotypes. Gene flux was predominantly driven by acquisition and recombination of prophage and associated cargo genes, with only occasional loss of these elements. The acquisition of large chromosomally-encoded genetic islands was limited, but notably, a feature of two recent pandemic clones (DT104 and monophasic S. Typhimurium ST34) of clade α (SGI-1 and SGI-4).

Characterization of Multidrug-Resistant Isolates of Salmonella enterica Serovars Heidelberg and Minnesota from Fresh Poultry Meat Imported to Portugal

Salmonella enterica serovars Heidelberg and Minnesota frequently display several genetic mobile elements making them potential spreaders of resistance genes. Here, we phenotypically determined the antibiotic resistance profile and subsequently performed whole-genome sequencing on 36 isolates recovered from samples of fresh poultry meat, within the Portuguese Official Inspection Plan for Imported Foodstuffs. Several isolates of both serovars showed high genetic relatedness either with isolates from raw poultry meat imported to the Netherlands from Brazil or with isolates from samples from the broiler production chain in Brazil. The multidrug-resistant (MDR) character was common to the vast majority (94.4%) of isolates from both serovars, and several isolates carried the plasmid IncA/C2 containing the β-lactamase gene bla_CMY-2 and IncX1 containing a type IV secretion system. These results somehow mirror the scenario observed in the Netherlands, showing the introduction, through fresh imported poultry meat in compliance with European legislation, of MDR Salmonella enterica serovars Heidelberg and Minnesota in Europe, with the potential spread of resistance markers. These data suggest the need to revise the hygiene criteria for foodstuffs monitoring before its placement on the market, with the determination of the resistome being an invaluable contribute to limit the dissemination of resistance markers.

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.