Background: Listeria monocytogenesClonal 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 ofL. monocytogenesto 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 FPEsL. monocytogenesisolates, 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 namedListeriagenomic island 3 (LGI3), composed by plant-associated loci and chromosomally integrating cadmium-resistance determinantscadA1C.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 inL. monocytogenesCCs 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 residentL. monocytogenesin FPEs to prevent contamination of RTE seafood.

Many of the bacterial genomic studies exploring evolution processes of the host adaptation focus on the accessory genome describing how the gains and losses of genes can explain the colonization of new habitats. Consequently, we developed a new approach focusing on the coregenome in order to describe the host adaptation of Salmonella serovars. In the present work, we propose bioinformatic tools allowing (i) robust phylogenetic inference based on SNPs and recombination events, (ii) identification of fixed SNPs and InDels distinguishing homoplastic and non-homoplastic coregenome variants, and (iii) gene-ontology enrichment analyses to describe metabolic processes involved in adaptation of Salmonella enterica subsp. enterica to mammalian- (S. Dublin), multi- (S. Enteritidis), and avian- (S. Pullorum and S. Gallinarum) hosts. The ‘VARCall’ workflow produced a robust phylogenetic inference confirming that the monophyletic clade S. Dublin diverged from the polyphyletic clade S. Enteritidis which includes the divergent clades S. Pullorum and S. Gallinarum (i). The scripts ‘phyloFixedVar’ and ‘FixedVar’ detected non-synonymous and non-homoplastic fixed variants supporting the phylogenetic reconstruction (ii). The scripts ‘GetGOxML’ and ‘EveryGO’ identified representative metabolic pathways related to host adaptation using the first gene-ontology enrichment analysis based on bacterial coregenome variants (iii). We propose in the present manuscript a new coregenome approach coupling identification of fixed SNPs and InDels with regards to inferred phylogenetic clades, and gene-ontology enrichment analysis in order to describe the adaptation of Salmonella serovars Dublin (i.e. mammalian-hosts), Enteritidis (i.e. multi-hosts), Pullorum (i.e. avian-hosts) and Gallinarum (i.e. avian-hosts) at the coregenome scale. All these polyvalent Bioinformatic tools can be applied on other bacterial genus without additional developments.

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.

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.

Bacillus thuringiensis (Bt), belonging to the Bacillus cereus (Bc) group, is commonly used as a biopesticide worldwide due to its ability to produce insecticidal crystals during sporulation. The use of Bt, especially subspecies aizawai and kurstaki, to control pests such as Lepidoptera, generally involves spraying mixtures containing spores and crystals on crops intended for human consumption. Recent studies have suggested that the consumption of commercial Bt strains may be responsible for foodborne outbreaks (FBOs). However, its genetic proximity to Bc strains has hindered the development of routine tests to discriminate Bt from other Bc, especially Bacillus cereus sensu stricto (Bc ss), well known for its involvement in FBOs. Here, to develop tools for the detection and the discrimination of Bt in food, we carried out a genome-wide association study (GWAS) on 286 complete genomes of Bc group strains to identify and validate in silico new molecular markers specific to different Bt subtypes. The analyses led to the determination and the in silico validation of 128 molecular markers specific to Bt, its subspecies aizawai, kurstaki and four previously described proximity clusters associated with these subspecies. We developed a command line tool based on a 14-marker workflow, to carry out a computational search for Bt-related markers from a putative Bc genome, thereby facilitating the detection of Bt of interest for food safety, especially in the context of FBOs.

In March 2024, the French genomic surveillance of enteric yersiniosis identified a cluster of lineage 2/3–9b Yersinia enterocolitica isolates, corresponding to bioserotype 2/O:9. An outbreak investigation was conducted to identify the source and implement control measures. A total of 175 confirmed cases were identified in France with sampling dates between 27 January and 23 August. Case interviews and trace-back investigations identified unpasteurised soft goat’s milk cheese from one manufacturer in France as the probable source of the outbreak. Yersinia enterocolitica belonging to the same cluster as the case isolates was isolated from 23 samples from the manufacturer and a goat farm supplying milk to the dairy. Cheeses from the manufacturer were recalled in France. A Rapid Alert System for Food and Feed (RASFF) notification was issued, allowing the identification of cheese distribution from the manufacturer to 29 countries. An Epipulse alert led to the identification of seven additional cases in Belgium, Luxembourg and Norway, illustrating the value of international warning systems. This outbreak demonstrates risks with consumption of unpasteurised cheese and emphasises the need of rigorous good hygiene practices on dairy farms and dairies, especially when milk is processed without pasteurisation.