Integration of metataxonomic data sets into microbial association networks highlights shared bacterial community dynamics in fermented vegetables.

The management of food fermentation is still largely based on empirical knowledge, as the dynamics of microbial communities and the underlying metabolic networks that produce safe and nutritious products remain beyond our understanding. Although these closed ecosystems contain relatively few taxa, they have not yet been thoroughly characterized with respect to how their microbial communities interact and dynamically evolve. However, with the increased availability of metataxonomic data sets on different fermented vegetables, it is now possible to gain a comprehensive understanding of the microbial relationships that structure plant fermentation. In this study, we applied a network-based approach to the integration of public metataxonomic 16S data sets targeting different fermented vegetables throughout time. Specifically, we aimed to explore, compare, and combine public 16S data sets to identify shared associations between amplicon sequence variants (ASVs) obtained from independent studies. The workflow includes steps for searching and selecting public time-series data sets and constructing association networks of ASVs based on co-abundance metrics. Networks for individual data sets are then integrated into a core network, highlighting significant associations. Microbial communities are identified based on the comparison and clustering of ASV networks using the "stochastic block model" method. When we applied this method to 10 public data sets (including a total of 931 samples) targeting five varieties of vegetables with different sampling times, we found that it was able to shed light on the dynamics of vegetable fermentation by characterizing the processes of community succession among different bacterial assemblages. IMPORTANCE: Within the growing body of research on the bacterial communities involved in the fermentation of vegetables, there is particular interest in discovering the species or consortia that drive different fermentation steps. This integrative analysis demonstrates that the reuse and integration of public microbiome data sets can provide new insights into a little-known biotope. Our most important finding is the recurrent but transient appearance, at the beginning of vegetable fermentation, of amplicon sequence variants (ASVs) belonging to Enterobacterales and their associations with ASVs belonging to Lactobacillales. These findings could be applied to the design of new fermented products.

Tell me if you prefer bovine or poultry sectors and I'll tell you who you are: Characterization of Salmonella enterica subsp. enterica serovar Mbandaka in France.

Introduction: In north-western France, Salmonella enterica susp. enterica serovar Mbandaka (S. Mbandaka) is most frequently isolated from bovine and dairy samples. While this serovar most often results in asymptomatic carriage, for a number of years it has caused episodes of abortions, which have serious economic consequences for the sector. Interestingly, this serovar is also isolated from Gallus gallus in the same geographic zone. Despite its prevalence in bovines in north-western France, S. Mbandaka has not been broadly studied at the genomic level, and its prevalence and host adaptation are still not fully understood. Methods: In this study, we analyzed the genomic diversity of 304 strains of S. Mbandaka isolated from the bovine and poultry sectors in this area over a period of 5 years. A phylogenetic analysis was carried out and two approaches were followed to identify conserved genes and mutations related to host associations. The first approach targeted the genes compiled in the MEGARESv2, Resfinder, VFDB and SPI databases. Plasmid and phage contents were also investigated. The second approach refers to an in-house algorithm developed for this study that computes sensitivity, specificity, and accuracy of accessory genes and core variants according to predefined genomes groups. Results and discussion: All the analyzed strains belong to the multi-locus sequence type profile ST413, and the phylogenomic analysis revealed main clustering by host (bovine and poultry), emphasizing the circulation of 12 different major clones, of which seven circulate in poultry and five in the bovine sector in France and a likely food production chain adaptation of these clones. All strains present resistance determinants including heavy metals and biocides that could explain the ability of this serovar to survive and persist in the environment, within herds, and in food processing plants. To explore the wild animal contribution to the spread of this serovar in north-western France, we retrieved S. Mbandaka genomes isolated from wild birds from EnteroBase and included them in the phylogenomic analysis together with our collection. Lastly, screening of accessory genes and major variants allowed us to identify conserved specific mutations characteristic of each major cluster. These mutations could be used to design useful probes for food safety surveillance.

A retrospective and regional approach assessing the genomic diversity of Salmonella Dublin.

From a historically rare serotype, Salmonella enterica subsp. enterica Dublin slowly became one of the most prevalent Salmonella in cattle and raw milk cheese in some regions of France. We present a retrospective genomic analysis of 480 S. Dublin isolates to address the context, evolutionary dynamics, local diversity and the genesis processes of regional S. Dublin outbreaks events between 2015 and 2017. Samples were clustered and assessed for correlation against metadata including isolation date, isolation matrices, geographical origin and epidemiological hypotheses. Significant findings can be drawn from this work. We found that the geographical distance was a major factor explaining genetic groups in the early stages of the cheese production processes (animals, farms) while down-the-line transformation steps were more likely to host genomic diversity. This supports the hypothesis of a generalised local persistence of strains from animal to finished products, with occasional migration. We also observed that the bacterial surveillance is representative of diversity, while targeted investigations without genomics evidence often included unrelated isolates. Combining both approaches in phylogeography methods allows a better representation of the dynamics, of outbreaks.

The Invasin and Complement-Resistance Protein Rck of Salmonella is More Widely Distributed than Previously Expected.

The rck open reading frame (ORF) on the pefI-srgC operon encodes an outer membrane protein responsible for invasion of nonphagocytic cell lines and resistance to complement-mediated killing. Until now, the rck ORF was only detected on the virulence plasmids of three serovars of Salmonella subsp. enterica (i.e., Bovismorbificans, Enteritidis, and Typhimurium). The increasing number of Salmonella genome sequences allowed us to use a combination of reference sequences and whole-genome multilocus sequence typing (wgMLST) data analysis to probe the presence of the operon and of rck in a wide array of isolates belonging to all Salmonella species and subspecies. We established the presence of partial or complete operons in 61 subsp. enterica serovars as well as in 4 other subspecies with various syntenies and frequencies. The rck ORF itself was retrieved in 36 subsp. enterica serovars and in two subspecies with either chromosomal or plasmid-borne localization. It displays high conservation of its sequence within the genus, and we demonstrated that most of the allelic variations identified did not alter the virulence properties of the protein. However, we demonstrated the importance of the residue at position 38 (at the level of the first extracellular loop of the protein) in the invasin function of Rck. Altogether, our results highlight that rck is not restricted to the three formerly identified serovars and could therefore have a more important role in virulence than previously expected. Moreover, this work raises questions about the mechanisms involved in rck acquisition and about virulence plasmid distribution and evolution. IMPORTANCE The foodborne pathogen Salmonella is responsible for a wide variety of pathologies depending on the infected host, the infecting serovars, and its set of virulence factors. However, the implication of each of these virulence factors and their role in the specific host-pathogen interplay are not fully understood. The significance of our research is in determining the distribution of one of these factors, the virulence plasmid-encoded invasin and resistance to complement killing protein Rck. In addition to providing elements of reflection concerning the mechanisms of acquisition of specific virulence genes in certain serotypes, this work will help to understand the role of Rck in the pathogenesis of Salmonella.

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

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

Comparative phenotypic, genotypic and genomic analyses of Bacillus thuringiensis associated with foodborne outbreaks in France.

Bacillus thuringiensis (Bt) belongs to the Bacillus cereus (Bc) group, well known as an etiological agent of foodborne outbreaks (FBOs). Bt distinguishes itself from other Bc by its ability to synthesize insecticidal crystals. However, the search for these crystals is not routinely performed in food safety or clinical investigation, and the actual involvement of Bt in the occurrence of FBOs is not known. In the present study, we reveal that Bt was detected in the context of 49 FBOs declared in France between 2007 and 2017. In 19 of these FBOs, Bt was the only microorganism detected, making it the most likely causal agent. Searching for its putative origin of contamination, we noticed that more than 50% of Bt isolates were collected from dishes containing raw vegetables, in particular tomatoes (48%). Moreover, the genomic characterization of isolates showed that most FBO-associated Bt isolates exhibited a quantified genomic proximity to Bt strains, used as biopesticides, especially those from subspecies aizawai and kurstaki. Taken together, these results strengthen the hypothesis of an agricultural origin for the Bt contamination and call for further investigations on Bt pesticides.

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.

NAuRA: Genomic Tool to Identify Staphylococcal Enterotoxins in Staphylococcus aureus Strains Responsible for FoodBorne Outbreaks.

Food contamination by staphylococcal enterotoxins (SEs) is responsible for many food poisoning outbreaks (FPOs) each year, and they represent the third leading cause of FPOs in Europe. SEs constitute a protein family with 27 proteins. However, enzyme immunoassays can only detect directly in food the five classical SEs (SEA-SEE). Thus, molecular characterization methods of strains found in food are now used for FPO investigations. Here, we describe the development and implementation of a genomic analysis tool called NAuRA (Nice automatic Research of alleles) that can detect the presence of 27 SEs genes in just one analysis- and create a database of allelic data and protein variants for harmonizing analyses. This tool uses genome assembly data and the 27 protein sequences of SEs. To include the different divergence levels between SE-coding genes, parameters of coverage and identity were generated from 10,000 simulations and a dataset of 244 assembled genomes from strains responsible for outbreaks in Europe as well as the RefSeq reference database. Based on phylogenetic inference performed using maximum-likelihood on the core genomes of the strains in this collection, we demonstrated that strains responsible for FPOs are distributed throughout the phylogenetic tree. Moreover, 71 toxin profiles were obtained using the NAuRA pipeline and these profiles do not follow the evolutionary history of strains. This study presents a pioneering method to investigate strains isolated from food at the genomic level and to analyze the diversity of all 27 SE-coding genes together.

Source Attribution Study of Sporadic Salmonella Derby Cases in France.

Salmonella enterica subsp. enterica serovar Derby is one of the most frequent causes of gastroenteritis in humans. In Europe, this pathogen is one of the top five most commonly reported serovars in human cases. In France, S. Derby has been among the ten most frequently isolated serovars in humans since the year 2000. The main animal hosts of this serovar are pigs and poultry, and white meat is the main source of human contamination. We have previously shown that this serovar is polyphyletic and that three distinct genetic lineages of S. Derby cohabit in France. Two of them are associated with pork and one with poultry. In this study, we conducted a source attribution study based on single nucleotide polymorphism analysis of a large collection of 440 S. Derby human and non-human isolates collected in 2014-2015, to determine the contribution of each lineage to human contamination. In France, the two lineages associated with pork strains, and corresponding to the multilocus sequence typing (MLST) profiles ST39-ST40 and ST682 were responsible for 94% of human contaminations. Interestingly, the ST40 profile is responsible for the majority of human cases (71%). An analysis of epidemiologic data and the structure of the pork sector in France allowed us to explain the spread and the sporadic pattern of human cases that occurred in the studied period.

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

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

Dynamics of mobile genetic elements of Listeria monocytogenes persisting in ready-to-eat seafood processing plants in France.

BACKGROUND: Listeria monocytogenes Clonal Complexes (CCs) have been epidemiologically associated with foods, especially ready-to-eat (RTE) products for which the most likely source of contamination depends on the occurrence of persisting clones in food-processing environments (FPEs). As the ability of L. monocytogenes to adapt to environmental stressors met in the food chain challenges the efforts to its eradication from FPEs, the threat of persistent strains to the food industry and public health authorities continues to rise. In this study, 94 food and FPEs L. monocytogenes isolates, representing persistent subtypes contaminating three French seafood facilities over 2-6 years, were whole-genome sequenced to characterize their genetic diversity and determine the biomarkers associated with long-term survival in FPEs. RESULTS: Food and FPEs isolates belonged to five CCs, comprising long-term intra- and inter-plant persisting clones. Mobile genetic elements (MGEs) such as plasmids, prophages and transposons were highly conserved within CCs, some of which harboured genes for resistance to chemical compounds and biocides used in the processing plants. Some of these genes were found in a 90.8 kbp plasmid, predicted to be" mobilizable", identical in isolates from CC204 and CC155, and highly similar to an 81.6 kbp plasmid from isolates belonging to CC7. These similarities suggest horizontal transfer between isolates, accompanied by deletion and homologous recombination in isolates from CC7. Prophage profiles characterized persistent clonal strains and several prophage-loci were plant-associated. Notably, a persistent clone from CC101 harboured a novel 31.5 kbp genomic island that we named Listeria genomic island 3 (LGI3), composed by plant-associated loci and chromosomally integrating cadmium-resistance determinants cadA1C. CONCLUSIONS: Genome-wide analysis indicated that inter- and intra-plant persisting clones harbour conserved MGEs, likely acquired in FPEs and maintained by selective pressures. The presence of closely related plasmids in L. monocytogenes CCs supports the hypothesis of horizontal gene transfer conferring enhanced survival to FPE-associated stressors, especially in hard-to-clean harbourage sites. Investigating the MGEs evolutionary and transmission dynamics provides additional resolution to trace-back potentially persistent clones. The biomarkers herein discovered provide new tools for better designing effective strategies for the removal or reduction of resident L. monocytogenes in FPEs to prevent contamination of RTE seafood.

A Simple and Robust Statistical Method to Define Genetic Relatedness of Samples Related to Outbreaks at the Genomic Scale - Application to Retrospective Salmonella Foodborne Outbreak Investigations.

The investigation of foodborne outbreaks (FBOs) from genomic data typically relies on inspecting the relatedness of samples through a phylogenomic tree computed on either SNPs, genes, kmers, or alleles (i.e., cgMLST and wgMLST). The phylogenomic reconstruction is often time-consuming, computation-intensive and depends on hidden assumptions, pipelines implementation and their parameterization. In the context of FBO investigations, robust links between isolates are required in a timely manner to trigger appropriate management actions. Here, we propose a non-parametric statistical method to assert the relatedness of samples (i.e., outbreak cases) or whether to reject them (i.e., non-outbreak cases). With typical computation running within minutes on a desktop computer, we benchmarked the ability of three non-parametric statistical tests (i.e., Wilcoxon rank-sum, Kolmogorov-Smirnov and Kruskal-Wallis) on six different genomic features (i.e., SNPs, SNPs excluding recombination events, genes, kmers, cgMLST alleles, and wgMLST alleles) to discriminate outbreak cases (i.e., positive control: C+) from non-outbreak cases (i.e., negative control: C-). We leveraged four well-characterized and retrospectively investigated FBOs of Salmonella Typhimurium and its monophasic variant S. 1,4,[5],12:i:- from France, setting positive and negative controls in all the assays. We show that the approaches relying on pairwise SNP differences distinguished all four considered outbreaks in contrast to the other tested genomic features (i.e., genes, kmers, cgMLST alleles, and wgMLST alleles). The freely available non-parametric method written in R has been designed to be independent of both the phylogenomic reconstruction and the detection methods of genomic features (i.e., SNPs, genes, kmers, or alleles), making it widely and easily usable to anybody working on genomic data from suspected samples.

Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale.

BACKGROUND: Salmonella enterica subsp. enterica is a public health issue related to food safety, and its adaptation to animal sources remains poorly described at the pangenome scale. Firstly, serovars presenting potential mono- and multi-animal sources were selected from a curated and synthetized subset of Enterobase. The corresponding sequencing reads were downloaded from the European Nucleotide Archive (ENA) providing a balanced dataset of 440 Salmonella genomes in terms of serovars and sources (i). Secondly, the coregenome variants and accessory genes were detected (ii). Thirdly, single nucleotide polymorphisms and small insertions/deletions from the coregenome, as well as the accessory genes were associated to animal sources based on a microbial Genome Wide Association Study (GWAS) integrating an advanced correction of the population structure (iii). Lastly, a Gene Ontology Enrichment Analysis (GOEA) was applied to emphasize metabolic pathways mainly impacted by the pangenomic mutations associated to animal sources (iv). RESULTS: Based on a genome dataset including Salmonella serovars from mono- and multi-animal sources (i), 19,130 accessory genes and 178,351 coregenome variants were identified (ii). Among these pangenomic mutations, 52 genomic signatures (iii) and 9 over-enriched metabolic signatures (iv) were associated to avian, bovine, swine and fish sources by GWAS and GOEA, respectively. CONCLUSIONS: Our results suggest that the genetic and metabolic determinants of Salmonella adaptation to animal sources may have been driven by the natural feeding environment of the animal, distinct livestock diets modified by human, environmental stimuli, physiological properties of the animal itself, and work habits for health protection of livestock.

Large-Scale Genomic Analyses and Toxinotyping of Clostridium perfringens Implicated in Foodborne Outbreaks in France.

Clostridium perfringens is both an ubiquitous environmental bacterium and the fourth most common causative agent of foodborne outbreaks (FBOs) in France and Europe. These outbreaks are known to be caused by C. perfringens enterotoxin (CPE) encoded by the cpe gene. However, additional information on the toxin/virulence gene content of C. perfringens has become available in the last few years. Therefore, to understand the enteropathogenicity of this bacterium, we need to describe the toxin and virulence genes content of strains involved in FBOs. In this study, we used a new real-time PCR typing technique based on a comprehensive set of 17 genes encoding virulence factors. The analysis was performed on a collection of 141 strains involved in 42 FBOs in the Paris region. It was combined with whole genome sequence (WGS) phylogenomic reconstruction, based on the coregenome single nucleotide polymorphisms (SNPs) of 58 isolates, representatives of the identified virulence gene profiles. Two or three different virulence gene profiles were detected in 10 FBOs, demonstrating that C. perfringens FBOs may be associated with heterogeneous strains. cpe-positive strains were isolated in 23 outbreaks, confirming the prominent role of CPE in pathogenicity. However, while C. perfringens was the sole pathogen isolated from the incriminated food, the cpe gene was not detected in strains related to 13 outbreaks. This result indicates either that the standard method was not able to isolate cpe+ strains or that the cpe gene may not be the only determinant of the enterotoxigenic potential of C. perfringens strains. Using phylogenomic reconstruction, we identified two clades distinguishing chromosomal cpe-positive from cpe-negative and plasmid-borne cpe. Important epidemiological information was also garnered from this phylogenomic reconstruction that revealed unexpected links between different outbreaks associated with closely related strains (seven SNP differences) and having common virulence gene profiles. This study provides new insight into the characterization of foodborne C. perfringens and highlights the potential of WGS for the investigation of FBOs.

Genetic Diversity of Salmonella Derby from the Poultry Sector in Europe.

Salmonella Derby (S. Derby) is emerging in Europe as a predominant serovar in fattening turkey flocks. This serovar was recorded as being predominant in the turkey sector in 2014 in the United Kingdom (UK). Only two years later, in 2016, it was also recorded in the turkey and broiler sectors in Ireland and Spain. These S. Derby isolates were characterised as members of the multilocus sequence type (MLST) profile 71 (ST71). For the first time, we characterise by whole genome sequencing (WGS) analysis a panel of 90 S. Derby ST71 genomes to understand the routes of transmission of this emerging pathogen within the poultry/turkey food trade. Selected panel included strains isolated as early as 2010 in five leading European g countries for turkey meat production. Twenty-one of the 90 genomes were extracted from a public database-Enterobase. Five of these originated from the United States (n=3), China (n=1) and Taiwan (n=1) isolated between 1986 and 2016. A phylogenomic analysis at the core-genome level revealed the presence of three groups. The largest group contained 97.5% of the European strains and included both, turkey and human isolates that were genetically related by an average of 35 ± 15 single nucleotide polymorphism substitutions (SNPs). To illustrate the diversity, the presence of antimicrobial resistance genes and phages were characteised in 30, S. Derby ST71 genomes, including 11 belonging to this study This study revealed an emergent turkey-related S. Derby ST71 clone circulating in at least five European countries (the UK, Germany, Poland, Italy, and France) since 2010 that causes human gastroenteritis. A matter of concern is the identification of a gyrA mutation involved in resistance to quinolone, present in the Italian genomes. Interestingly, the diversity of phages seems to be related to the geographic origins. These results constitute a baseline for following the spread of this emerging pathogen and identifying appropriate monitoring and prevention measures.

Correction for Sévellec et al., "Complete Genome Sequence of Salmonella enterica subsp. enterica Serotype Derby, Associated with the Pork Sector in France".

[This corrects the article DOI: 10.1128/MRA.01027-18.].

Machine Learning Methods as a Tool for Predicting Risk of Illness Applying Next-Generation Sequencing Data.

Next-generation sequencing (NGS) data present an untapped potential to improve microbial risk assessment (MRA) through increased specificity and redefinition of the hazard. Most of the MRA models do not account for differences in survivability and virulence among strains. The potential of machine learning algorithms for predicting the risk/health burden at the population level while inputting large and complex NGS data was explored with Listeria monocytogenes as a case study. Listeria data consisted of a percentage similarity matrix from genome assemblies of 38 and 207 strains of clinical and food origin, respectively. Basic Local Alignment (BLAST) was used to align the assemblies against a database of 136 virulence and stress resistance genes. The outcome variable was frequency of illness, which is the percentage of reported cases associated with each strain. These frequency data were discretized into seven ordinal outcome categories and used for supervised machine learning and model selection from five ensemble algorithms. There was no significant difference in accuracy between the models, and support vector machine with linear kernel was chosen for further inference (accuracy of 89% [95% CI: 68%, 97%]). The virulence genes FAM002725, FAM002728, FAM002729, InlF, InlJ, Inlk, IisY, IisD, IisX, IisH, IisB, lmo2026, and FAM003296 were important predictors of higher frequency of illness. InlF was uniquely truncated in the sequence type 121 strains. Most important risk predictor genes occurred at highest prevalence among strains from ready-to-eat, dairy, and composite foods. We foresee that the findings and approaches described offer the potential for rethinking the current approaches in MRA.

Insights from genome-wide approaches to identify variants associated to phenotypes at pan-genome scale: Application to L. monocytogenes' ability to grow in cold conditions.

Intraspecific variability of the behavior of most foodborne pathogens is well described and taken into account in Quantitative Microbial Risk Assessment (QMRA), but factors (strain origin, serotype, …) explaining these differences are scarce or contradictory between studies. Nowadays, Whole Genome Sequencing (WGS) offers new opportunities to explain intraspecific variability of food pathogens, based on various recently published bioinformatics tools. The objective of this study is to get a better insight into different existing bioinformatics approaches to associate bacterial phenotype(s) and genotype(s). Therefore, a dataset of 51 L. monocytogenes strains, isolated from multiple sources (i.e. different food matrices and environments) and belonging to 17 clonal complexes (CC), were selected to represent large population diversity. Furthermore, the phenotypic variability of growth at low temperature was determined (i.e. qualitative phenotype), and the whole genomes of selected strains were sequenced. The almost exhaustive gene content, as well as the core genome SNPs based phylogenetic reconstruction, were derived from the whole sequenced genomes. A Bayesian inference method was applied to identify the branches on which the phenotype distribution evolves within sub-lineages. Two different Genome Wide Association Studies (i.e. gene- and SNP-based GWAS) were independently performed in order to link genetic mutations to the phenotype of interest. The genomic analyses presented in this study were successfully applied on the selected dataset. The Bayesian phylogenetic approach emphasized an association with "slow" growth ability at 2 °C of the lineage I, as well as CC9 of the lineage II. Moreover, both gene- and SNP-GWAS approaches displayed significant statistical associations with the tested phenotype. A list of 114 significantly associated genes, including genes already known to be involved in the cold adaption mechanism of L. monocytogenes and genes associated to mobile genetic elements (MGE), resulted from the gene-GWAS. On the other hand, a group of 184 highly associated SNPs were highlighted by SNP-GWAS, including SNPs detected in genes which were already likely involved in cold adaption; hypothetical proteins; and intergenic regions where for example promotors and regulators can be located. The successful application of combined bioinformatics approaches associating WGS-genotypes and specific phenotypes, could contribute to improve prediction of microbial behaviors in food. The implementation of this information in hazard identification and exposure assessment processes will open new possibilities to feed QMRA-models.

Complete Genome Sequence of Salmonella enterica subsp. enterica Serotype Derby, Associated with the Pork Sector in France.

In the European Union, Salmonella enterica subsp. enterica serovar Derby is the most abundant serotype isolated from pork. Recent studies have shown that this serotype is polyphyletic. However, one main genomic lineage, characterized by sequence type 40 (ST40), the presence of the Salmonella pathogenicity island 23, and showing resistance to streptomycin, sulphonamides, and tetracycline (STR-SSS-TET), is pork associated. Here, we describe the complete genome sequence of a strain from this lineage isolated in France.