A Critical Review of Risk Assessment Models for Listeria monocytogenes in Produce

A review of quantitative risk assessment (QRA) models of Listeria monocytogenes in produce was carried out, with the objective of appraising and contrasting the effectiveness of the control strategies placed along the food chains. Despite nine of the thirteen QRA models recovered being focused on fresh or RTE leafy greens, none of them represented important factors or sources of contamination in the primary production, such as the type of cultivation, water, fertilisers or irrigation method/practices. Cross-contamination at processing and during consumer's handling was modelled using transfer rates, which were shown to moderately drive the final risk of listeriosis, therefore highlighting the importance of accurately representing the transfer coefficient parameters. Many QRA models coincided in the fact that temperature fluctuations at retail or temperature abuse at home were key factors contributing to increasing the risk of listeriosis. In addition to a primary module that could help assess current on-farm practices and potential control measures, future QRA models for minimally processed produce should also contain a refined sanitisation module able to estimate the effectiveness of various sanitisers as a function of type, concentration and exposure time. Finally, L. monocytogenes growth in the products down the supply chain should be estimated by using realistic time-temperature trajectories, and validated microbial kinetic parameters, both of them currently available in the literature.

Updated Parameters for Listeria monocytogenes Dose-Response Model Considering Pathogen Virulence and Age and Sex of Consumer

Better knowledge regarding the Listeria monocytogenes dose-response (DR) model is needed to refine the assessment of the risk of foodborne listeriosis. In 2018, the European Food Safety Agency (EFSA) derived a lognormal Poisson DR model for 14 different age-sex sub-groups, marginally to strain virulence. In the present study, new sets of parameters are developed by integrating the EFSA model for these sub-groups together with three classes of strain virulence characteristics ("less virulent", "virulent", and "more virulent"). Considering classes of virulence leads to estimated relative risks (RRs) of listeriosis following the ingestion of 1000 bacteria of "less virulent" vs. "more virulent" strains ranging from 21.6 to 24.1, depending on the sub-group. These relatively low RRs when compared with RRs linked to comorbidities described in the literature suggest that the influence of comorbidity on the occurrence of invasive listeriosis for a given exposure is much more important than the influence of the virulence of the strains. The updated model parameters allow better prediction of the risk of invasive listeriosis across a population of interest, provided the necessary data on population demographics and the proportional contribution of strain virulence classes in food products of interest are available. An R package is made available to facilitate the use of these dose-response models.

A Critical Review of Risk Assessment Models for Listeria monocytogenes in Seafood

Invasive listeriosis, due to its severe nature in susceptible populations, has been the focus of many quantitative risk assessment (QRA) models aiming to provide a valuable guide in future risk management efforts. A review of the published QRA models of Listeria monocytogenes in seafood was performed, with the objective of appraising the effectiveness of the control strategies at different points along the food chain. It is worth noting, however, that the outcomes of a QRA model are context-specific, and influenced by the country and target population, the assumptions that are employed, and the model architecture itself. Studies containing QRA models were retrieved through a literature search using properly connected keywords on Scopus and PubMed®. All 13 QRA models that were recovered were of short scope, covering, at most, the period from the end of processing to consumption; the majority (85%) focused on smoked or gravad fish. Since the modelled pathways commenced with the packaged product, none of the QRA models addressed cross-contamination events. Many models agreed that keeping the product's temperature at 4.0-4.5 °C leads to greater reductions in the final risk of listeriosis than reducing the shelf life by one week and that the effectiveness of both measures can be surpassed by reducing the initial occurrence of L. monocytogenes in the product (at the end of processing). It is, therefore, necessary that future QRA models for RTE seafood contain a processing module that can provide insight into intervention strategies that can retard L. monocytogenes' growth, such as the use of bacteriocins, ad hoc starter cultures and/or organic acids, and other strategies seeking to reduce cross-contamination at the facilities, such as stringent controls for sanitation procedures. Since risk estimates were shown to be moderately driven by growth kinetic parameters, namely, the exponential growth rate, the minimum temperature for growth, and the maximum population density, further work is needed to reduce uncertainties.

A Critical Review of Risk Assessment Models for Listeria monocytogenes in Meat and Meat Products

A review of the published quantitative risk assessment (QRA) models of L. monocytogenes in meat and meat products was performed, with the objective of appraising the intervention strategies deemed suitable for implementation along the food chain as well as their relative effectiveness. A systematic review retrieved 23 QRA models; most of them (87%) focused on ready-to-eat meat products and the majority (78%) also covered short supply chains (end processing/retail to consumption, or consumption only). The processing-to-table scope was the choice of models for processed meats such as chorizo, bulk-cooked meat, fermented sausage and dry-cured pork, in which the effects of processing were simulated. Sensitivity analysis demonstrated the importance of obtaining accurate estimates for lag time, growth rate and maximum microbial density, in particular when affected by growth inhibitors and lactic acid bacteria. In the case of deli meats, QRA models showed that delicatessen meats sliced at retail were associated with a higher risk of listeriosis than manufacture pre-packed deli meats. Many models converged on the fact that (1) controlling cold storage temperature led to greater reductions in the final risk than decreasing the time to consumption and, furthermore, that (2) lower numbers and less prevalence of L. monocytogenes at the end of processing were far more effective than keeping low temperatures and/or short times during retail and/or home storage. Therefore, future listeriosis QRA models for meat products should encompass a processing module in order to assess the intervention strategies that lead to lower numbers and prevalence, such as the use of bio-preservation and novel technologies. Future models should be built upon accurate microbial kinetic parameters, and should realistically represent cross-contamination events along the food chain.

Comparison of Listeria monocytogenes alternative detection methods for food microbiology official controls in Europe

Listeriosis remains one of the most severe foodborne diseases regarding mortality rate. L. monocytogenes can grow under stressful conditions and contaminate various food categories. Regulation (EC) No 2073/2005 modified on microbiological criteria for foodstuffs includes either qualitative or quantitative food safety criteria on L. monocytogenes and mentions EN ISO 11290-1 Standard as the reference detection method. The use of alternative certified methods is allowed at European level for own checks and by certain European countries in the frame of official controls. Here we report the comparison of alternative certified methods for the detection of L. monocytogenes, according to the criteria expressed by the National Reference Laboratories (NRLs) for L. monocytogenes network in Europe, through an enquiry conducted in 2022 by the European Union Reference Laboratory (EURL) for L. monocytogenes. Emerging detection methods make it possible to diversify the panel of available detection methods. Each method has advantages and limitations and the choice of a method depends on the needs, organization and objectives of the laboratory. Our study focused on a comparison by category in order to provide an overview of the differences between the methods as complete and simple as possible, mainly by examining the certification files of the alternative methods, which constitute a reliable and publicly accessible source of information about performances of the methods.

A Critical Review of Risk Assessment Models for Listeria monocytogenes in Dairy Products

A review of the published quantitative risk assessment (QRA) models of L. monocytogenes in dairy products was undertaken in order to identify and appraise the relative effectiveness of control measures and intervention strategies implemented at primary production, processing, retail, and consumer practices. A systematic literature search retrieved 18 QRA models, most of them (9) investigated raw and pasteurized milk cheeses, with the majority covering long supply chains (4 farm-to-table and 3 processing-to-table scopes). On-farm contamination sources, either from shedding animals or from the broad environment, have been demonstrated by different QRA models to impact the risk of listeriosis, in particular for raw milk cheeses. Through scenarios and sensitivity analysis, QRA models demonstrated the importance of the modeled growth rate and lag phase duration and showed that the risk contribution of consumers' practices is greater than in retail conditions. Storage temperature was proven to be more determinant of the final risk than storage time. Despite the pathogen's known ability to reside in damp spots or niches, re-contamination and/or cross-contamination were modeled in only two QRA studies. Future QRA models in dairy products should entail the full farm-to-table scope, should represent cross-contamination and the use of novel technologies, and should estimate L. monocytogenes growth more accurately by means of better-informed kinetic parameters and realistic time-temperature trajectories.

An agent-based model to simulate SARS-CoV-2 contamination of surfaces and meat cuts in processing plants

At the beginning of the COVID-19 pandemic, several contamination clusters were reported in food-processing plants in France and several countries worldwide. Therefore, a need arose to better understand viral transmission in such occupational environments from multiple perspectives: the protection of workers in hotspots of viral circulation; the prevention of supply disruption due to the closure of plants; and the prevention of cluster expansion due to exports of food products contaminated by the virus to other locations. This paper outlines a simulation-based approach (using agent-based models) to study the effects of measures taken to prevent the contamination of workers, surfaces, and food products. The model includes user-defined parameters to integrate characteristics relating to SARS-CoV-2 (variant of concern to be considered, symptom onset…), food-processing plants (dimensions, ventilation…), and other sociodemographic transmission factors based on laboratory experiments as well as industrial and epidemiological investigations. Simulations were performed for a typical meat-processing plant in different scenarios for illustration purposes. The results suggested that increasing the mask-wearing ratio led to great reductions in the probability of observing clusters of more than 25 infections. In the case of clusters, masks being worn by all workers limited the presence of contamination (defined as levels of at least 5 log10 viral RNA copies) on meat cuts at less than 0.05 % and maintained the production capacity of the plant at optimal levels. Increasing the average distance between two workers from less than 1 m to more than 2 m decreased the cluster-occurrence probability by up to 15 % as well as contamination of food products during cluster situations. The developed approach can open up several perspectives in terms of potential communication-support tools for the agri-food sector and further reuses or adaptations for other hazards and occupational environments.

Assessment of the relationship between the MLST genetic diversity of Listeria monocytogenes and growth under selective and non-selective conditions

Listeria monocytogenes can grow under stressful conditions and contaminate various food categories. Progresss in DNA sequencing-based identification methods, such as multi-locus sequence typing (MLST) now allow for more accurate characterization of pathogens. L. monocytogenes MLST genetic diversity is reflected by the different prevalence of the "clonal complexes" (CCs) in foods or infections. Better understanding of the growth potentials of L. monocytogenes is essential for quantitative risk assessment and efficient detection across CCs genetic diversity. Using optical density measurements taken with an automated spectrophotometer, we compared the maximal growth rate and lag phase of 39 strains from 13 different CCs and various food origins, in 3 broths mimicking stresful food conditions (8 °C, aw 0.95 and pH5) and in ISO Standard enrichment broths (Half Fraser and Fraser). This is important as growth could influence risk through pathogen multiplication in food. Besides, enrichment problems could lead to a lack of detection of some CCs. Despite small differences highlighting natural intraspecific variability, our results show that growth performances of L. monocytogenes strains under the conditions tested in selective and non-selective broth do not appear to be strongly correlated to CCs and cannot explain higher CC "virulence" or prevalence.

Harmonized approach to estimate the burden of disease of dietary exposure to four chemical contaminants - A French study

Exposure to chemical contaminants found in foods has been associated with various adverse health effects. Burden of disease studies are increasingly used to estimate the public health impact of such exposures. The aims of this study were to estimate the burden of disease due to dietary exposure to four chemicals in France in 2019 (lead (Pb), cadmium (Cd), methylmercury (MeHg), and inorganic arsenic (i-As)), and to develop harmonized methods that can be applied for other chemicals and countries. We used national food consumption data from the third French national food consumption survey, chemical food monitoring data from the Second French Total Diet Study (TDS), dose-response data and disability weights from scientific literature, and disease incidence and demographics from national statistics. We adopted a risk assessment approach to estimate disease burden, incidence, mortality, and Disability-Adjusted Life Years (DALYs) attributable to dietary exposure to the chemicals. In all models, we harmonized food classification and exposure assessment. We propagated uncertainty through the calculations using Monte Carlo simulation. We estimated that, among these chemicals, i-As and Pb were responsible for the highest disease burden. i-As was estimated to cause 820 DALYs, or approximately 1.25 DALYs per 100,000 inhabitants. The estimated burden of Pb was 1834 to 5936 DALYs, or 2.7 (lower bound) to 8.96 (upper bound) DALYs/100,000. The burden of MeHg (192 DALYs), and Cd (0 DALY) was substantially lower. The foods contributing most to disease burden was drinks (30 %), "other foods" (mostly composite dishes) (19 %), and fish and seafood (7 %). Interpretation of estimates needs to consider all underlying uncertainties, linked with data and knowledge gaps. The harmonized models are the first to make use of data from TDS, which are available in several other countries. Thus, they can be applied to estimate the burden and to rank food-associated chemicals at national level.

Physical Treatments to Control Clostridium botulinum Hazards in Food

Clostridium botulinum produces Botulinum neurotoxins (BoNTs), causing a rare but potentially deadly type of food poisoning called foodborne botulism. This review aims to provide information on the bacterium, spores, toxins, and botulisms, and describe the use of physical treatments (e.g., heating, pressure, irradiation, and other emerging technologies) to control this biological hazard in food. As the spores of this bacterium can resist various harsh environmental conditions, such as high temperatures, the thermal inactivation of 12-log of C. botulinum type A spores remains the standard for the commercial sterilization of food products. However, recent advancements in non-thermal physical treatments present an alternative to thermal sterilization with some limitations. Low- (<2 kGy) and medium (3-5 kGy)-dose ionizing irradiations are effective for a log reduction of vegetative cells and spores, respectively; however, very high doses (>10 kGy) are required to inactivate BoNTs. High-pressure processing (HPP), even at 1.5 GPa, does not inactivate the spores and requires heat combination to achieve its goal. Other emerging technologies have also shown some promise against vegetative cells and spores; however, their application to C. botulinum is very limited. Various factors related to bacteria (e.g., vegetative stage, growth conditions, injury status, type of bacteria, etc.) food matrix (e.g., compositions, state, pH, temperature, aw, etc.), and the method (e.g., power, energy, frequency, distance from the source to target, etc.) influence the efficacy of these treatments against C. botulinum. Moreover, the mode of action of different physical technologies is different, which provides an opportunity to combine different physical treatment methods in order to achieve additive and/or synergistic effects. This review is intended to guide the decision-makers, researchers, and educators in using physical treatments to control C. botulinum hazards.

Clostridium botulinum type C, D, C/D, and D/C: An update

Clostridium botulinum is the main causative agent of botulism, a neurological disease encountered in humans as well as animals. Nine types of botulinum neurotoxins (BoNTs) have been described so far. Amongst these "toxinotypes," the A, the B and E are the most frequently encountered in humans while the C, D, C/D and D/C are mostly affecting domestic and wild birds as well as cattle. In France for instance, many cases and outbreaks are reported in these animal species every year. However, underestimation is very likely at least for avifauna species where the detection of dead animals can be challenging. Knowledge about BoNTs C, D, C/D, and D/C and the diseases they cause in animals and humans is still scarce and unclear. Specifically, the potential role of animal botulism outbreaks in cattle and poultry as a source of human illness needs to be further assessed. In this narrative review, we present the current knowledge about toxinotypes C, D, C/D, and D/C in cattle and poultry with, amongst various other aspects, their epidemiological cycles. We also discuss the zoonotic potential of these toxinotypes and some possible ways of risk mitigation. An adapted and effective management of botulism outbreaks in livestock also requires a better understanding of these less common and known toxinotypes.

Risk of Monkeypox virus (MPXV) transmission through the handling and consumption of food

Monkeypox (MPX) is a zoonotic infectious disease caused by Monkeypox virus (MPXV), an enveloped DNA virus belonging to the Poxviridae family and the Orthopoxvirus genus. Since early May 2022, a growing number of human cases of Monkeypox have been reported in non-endemic countries, with no history of contact with animals imported from endemic and enzootic areas, or travel to an area where the virus usually circulated before May 2022. This qualitative risk assessment aimed to investigate the probability that MPXV transmission occurs through food during its handling and consumption. The risk assessment used "top-down" (based on epidemiological data) and "bottom-up" (following the agent through the food chain to assess the risk of foodborne transmission to human) approaches, which were combined. The "top-down" approach first concluded that bushmeat was the only food suspected as a source of contamination in recorded cases of MPXV, by contact or ingestion. The "bottom-up" approach then evaluated the chain of events required for a human to become ill after handling or consuming food. This approach involves several conditions: (i) the food must be contaminated with MPXV (naturally, by an infected handler or after contact with a contaminated surface); (ii) the food must contain viable virus when it reaches the handler or consumer; (iii) the person must be exposed to the virus and; (iv) the person must be infected after exposure. Throughout the risk assessment, some data gaps were identified and highlighted. The conclusions of the top-down and bottom-up approaches are consistent and suggest that the risk of transmission of MPXV through food is hypothetical and that such an occurrence was never reported. In case of contamination, cooking (e.g., 12 min at 70°C) could be considered effective in inactivating Poxviridae in foods. Recommendations for risk management are proposed. To our knowledge, this is the first risk assessment performed on foodborne transmission of MPXV.

Human and animal botulism surveillance in France from 2008 to 2019

Botulism is a human and animal neurological disease caused by the action of bacterial neurotoxins (botulinum toxins) produced by bacteria from the genus Clostridium. This disease induces flaccid paralysis that can result in respiratory paralysis and heart failure. Due to its serious potential impact on public health, botulism is a closely monitored notifiable disease in France through a case-based passive surveillance system. In humans, this disease is rare, with an average of 10 outbreaks reported each year, mainly due to the consumption of contaminated foods. Type B and to a lesser extend type A are responsible for the majority of cases of foodborne botulism. Each year, an average of 30 outbreaks are recorded on poultry farms, about 20 cases in wild birds and about 10 outbreaks in cattle, involving a large number of animals. Mosaic forms C/D and D/C in birds and cattle, respectively, are the predominant types in animals in France. Types C and D have also been observed to a lesser extent in animals. With the exception of botulinum toxin E, which was exceptionally detected throughout the period in wild birds, the types of botulism found in animal outbreaks are different from those identified in human outbreaks over the last ten years in France and no human botulism outbreaks investigated have been linked to animal botulism. In line with the One Health concept, we present the first integrative approach to the routine surveillance of botulism in humans and animals in France.

Database of SARS-CoV-2 and coronaviruses kinetics relevant for assessing persistence in food processing plants

SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2), a virus causing severe acute respiratory disease in humans, emerged in late 2019. This respiratory virus can spread via aerosols, fomites, contaminated hands or surfaces as for other coronaviruses. Studying their persistence under different environmental conditions represents a key step for better understanding the virus transmission. This work aimed to present a reproducible procedure for collecting data of stability and inactivation kinetics from the scientific literature. The aim was to identify data useful for characterizing the persistence of viruses in the food production plants. As a result, a large dataset related to persistence on matrices or in liquid media under different environmental conditions is presented. This procedure, combining bibliographic survey, data digitalization techniques and predictive microbiological modelling, identified 65 research articles providing 455 coronaviruses kinetics. A ranking step as well as a technical validation with a Gage Repeatability & Reproducibility process were performed to check the quality of the kinetics. All data were deposited in public repositories for future uses by other researchers.

Measurement(s)	Decimal reduction time
Technology Type(s)	Modelling of kinetics
Factor Type(s)	Temperature • Relative humidity • Genus • Strain
Sample Characteristic - Organism	Coronaviridae

Genomic elements located in the accessory repertoire drive the adaptation to biocides in Listeria monocytogenes strains from different ecological niches

In response to the massive use of biocides for controlling Listeria monocytogenes (hereafter Lm) contaminations along the food chain, strains showing biocide tolerance emerged. Here, accessory genomic elements were associated with biocide tolerance through pangenome-wide associations performed on 197 Lm strains from different lineages, ecological, geographical and temporal origins. Mobile elements, including prophage-related loci, the Tn6188_qacH transposon and pLMST6_emrC plasmid, were widespread across lineage I and II food strains and associated with tolerance to benzalkonium-chloride (BC), a quaternary ammonium compound (QAC) widely used in food processing. The pLMST6_emrC was also associated with tolerance to another QAC, the didecyldimethylammonium-chloride, displaying a pleiotropic effect. While no associations were detected for chemically reactive biocides (alcohols and chlorines), genes encoding for cell-surface proteins were associated with BC or polymeric biguanide tolerance. The latter was restricted to lineage I strains from animal and the environment. In conclusion, different genetic markers, with polygenic nature or not, appear to have driven the Lm adaptation to biocide, especially in food strains but also from animal and the environment. These markers could aid to monitor and predict the spread of biocide tolerant Lm genotypes across different ecological niches, finally reducing the risk of such strains in food industrial settings.

Genetic, physiological, and cellular heterogeneities of bacterial pathogens in food matrices: Consequences for food safety

In complex food systems, bacteria live in heterogeneous microstructures, and the population displays phenotypic heterogeneities at the single-cell level. This review provides an overview of spatiotemporal drivers of phenotypic heterogeneity of bacterial pathogens in food matrices at three levels. The first level is the genotypic heterogeneity due to the possibility for various strains of a given species to contaminate food, each of them having specific genetic features. Then, physiological heterogeneities are induced within the same strain, due to specific microenvironments and heterogeneous adaptative responses to the food microstructure. The third level of phenotypic heterogeneity is related to cellular heterogeneity of the same strain in a specific microenvironment. Finally, we consider how these phenotypic heterogeneities at the single-cell level could be implemented in mathematical models to predict bacterial behavior and help ensure microbiological food safety.

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.

Listeria monocytogenes: Investigation of Fitness in Soil Does Not Support the Relevance of Ecotypes

Listeria monocytogenes (Lm) is a ubiquitous bacterium that causes the serious foodborne illness listeriosis. Although soil is a primary reservoir and a central habitat for Lm, little information is available on the genetic features underlying the fitness of Lm strains in this complex habitat. The aim of this study was to identify (i) correlations between the strains fitness in soil, their origin and their phylogenetic position (ii) identify genetic markers allowing Lm to survive in the soil. To this end, we assembled a balanced panel of 216 Lm strains isolated from three major ecological compartments (outdoor environment, animal hosts, and food) and from 33 clonal complexes occurring worldwide. The ability of the 216 strains to survive in soil was tested phenotypically. Hierarchical clustering identified three phenotypic groups according to the survival rate (SR): phenotype 1 “poor survivors” (SR &lt; 2%), phenotype 2 “moderate survivors” (2% &lt; SR &lt; 5%) and phenotype 3 “good survivors” (SR &gt; 5%). Survival in soil depended neither on strains’ origin nor on their phylogenetic position. Genome-wide-association studies demonstrated that a greater number of genes specifically associated with a good survival in soil was found in lineage II strains (57 genes) than in lineage I strains (28 genes). Soil fitness was mainly associated with variations in genes (i) coding membrane proteins, transcription regulators, and stress resistance genes in both lineages (ii) coding proteins related to motility and (iii) of the category “phage-related genes.” The cumulative effect of these small genomic variations resulted in significant increase of soil fitness.

Salmonella enterica subsp. enterica Welikade: guideline for phylogenetic analysis of serovars rarely involved in foodborne outbreaks

BACKGROUND

Salmonella spp. is a major foodborne pathogen with a wide variety of serovars associated with human cases and food sources. Nevertheless, in Europe a panel of ten serovars is responsible for up to 80% of confirmed human cases. Clustering studies by single nucleotide polymorphism (SNP) core-genome phylogenetic analysis of outbreaks due to these major serovars are simplified by the availability of many complete genomes in the free access databases. This is not the case for outbreaks due to less common serovars, such as Welikade, for which no reference genomes are available. In this study, we propose a method to solve this problem. We propose to perform a core genome MLST (cgMLST) analysis based on hierarchical clustering using the free-access EnteroBase to select the most suitable genome to use as a reference for SNP phylogenetic analysis. In this study, we applied this protocol to a retrospective analysis of a Salmonella enterica serovar Welikade (S. Welikade) foodborne outbreak that occurred in France in 2016. Finally, we compared the cgMLST and SNP analyses. SNP phylogenetic reconstruction was carried out considering the effect of recombination events identified by the ClonalFrameML tool. The accessory genome was also explored by phage content and virulome analyses.

RESULTS

Our findings revealed high clustering concordance using cgMLST and SNP analyses. Nevertheless, SNP analysis allowed for better assessment of the genetic distance among strains. The results revealed epidemic clones of S. Welikade circulating within the poultry and dairy sectors in France, responsible for sporadic and non-sporadic human cases between 2012 and 2019.

CONCLUSIONS

This study increases knowledge on this poorly described serovar and enriches public genome databases with 42 genomes from human and non-human S. Welikade strains, including the isolate collected in 1956 in Sri Lanka, which gave the name to this serovar. This is the first genomic analysis of an outbreak due to S. Welikade described to date.

Review: On-farm and processing factors affecting bovine carcass and meat quality

This paper reviews the current state of knowledge on beef carcass and meat quality, with particular emphasis on on-farm and processing factors associated with its high and inconsistent variability. The diversity of livestock systems comes from the diversity of breeds (dairy or beef), ages and gender (bulls, steers, heifers, cull cows) used to produce either mainly beef or beef and milk. In addition, there are factors linked to farming practices (including diet, especially grazing) which significantly influence the sensory, nutritional, technological and extrinsic (such as image) quality attributes of meat. These can become factors of positive differentiation when controlled by the application and certification of technical specifications. Finally, preslaughter (such as stress), slaughter (such as the chilling and hanging method of carcasses) and postslaughter (such as ageing, packaging and cooking) conditions have a strong influence on the microbiological, sensory, technological and image quality attributes of beef. In this review, potential synergisms or antagonisms between the different quality attributes are highlighted. For example, finishing cattle on grass, compared to indoor fattening on a high concentrate diet, has the advantage of producing leaner meat with a higher proportion of omega-3 fatty acids while exhibiting superior oxidative stability, but with the consequence of a darker meat colour and lower productivity, as well as higher seasonality and land surface requirements. Moreover, the control of on-farm factors is often guided by productivity (growth rate, feed conversion ratio) and carcass quality attributes (weight, conformation and fatness). Genetic selection has often been oriented in this direction, without taking other quality attributes into account. Finally, the interactions between all these factors (and especially between on-farm and slaughter or processing factors) are not considered in the quality grading schemes in European countries. This means that positive efforts at farm level may be mitigated or even eliminated by poor slaughtering or processing conditions. All these considerations explain why between-animal variability in quality can be high, even when animals come from the same farming system. The ability to predict the sensory and nutritional properties of meat according to production factors has become a major objective of the supply chain.

Review: Production factors affecting poultry carcass and meat quality attributes

Poultry meat mainly comes from standard production system using high growth rate strains reared under indoor intensive conditions. However, it is possible to find also different alternative systems using outdoor extensive rearing conditions and slow-growing lines. These different production systems can affect carcass and meat quality. In this review, quality has been broken down into six properties: commercial, organoleptic, nutritional, technological, sanitary and image, the latter covering the ethical, cultural and environmental dimensions associated with the way the meat is produced, as well as its origin and being particularly valued in many quality labels. The quality of meat is built and can deteriorate along the continuum from the conception of the animal to the fork. Our review details the different factors implicated in the determinism of poultry meat properties and pinpoints critical periods, such as the preslaughter and slaughter periods, and key factors, such as the feeding regimen, via its direct effect on the fatty acid profile, the antioxidant and volatile compound contents, and indirect effects mediated via the growth rate of the bird. Our review also highlights potential antagonisms between different dimensions of quality. The genetic selection for breast meat yield, for example, has been effective in producing carcasses with higher meat yield, but resulting since a decade in the increased occurrence of quality defects and myopathies (white striping, wooden breast, spaghetti meat and deep pectoral disease). Outdoor access has positive effects on the image and nutritional properties (through its effect on the fatty acid profile of meat lipids), but it increases the exposition risk to environmental contaminants and pathologies (parasites, virus, bacteria); it also increases the variability in meat quality linked to the variability of animal performance and slaughter age. The orientation towards more agro-ecological low-input farming systems may present benefits for the image and nutritional properties, but also risks for the commercial (low carcass weight and low breast yield, irregularity in supply), organoleptic (stronger flavour, less tender and darker colour of the meat) and in terms of variability of the different properties that constitute quality. Efforts should be made in the future to better take into account the various dimensions of quality, in consumer information, payment to farmers and genetic selection.

Review: Factors affecting sheep carcass and meat quality attributes

Sheep meat comes from a wide variety of farming systems utilising outdoor extensive to indoor intensive with animals of various ages at slaughter. In Europe, slaughter may occur from 4 weeks of age in suckling light lambs to adult ages. More than any other animal species used for meat production, there are strong country-specific preferences for sheep meat quality linked to production system characteristics such as dairy or grassland-based systems. This article critically reviews the current state of knowledge on factors affecting sheep carcass and meat quality. Quality has been broken down into six core attributes: commercial, organoleptic, nutritional, technological, safety and image, the latter covering aspects of ethics, culture and environment associated with the way the meat is produced and its origin, which are particularly valued in the many quality labels in Europe. The quality of meat is built but can also deteriorate along the continuum from the conception of the animal to the consumer. Our review pinpoints critical periods, such as the gestation and the preslaughter and slaughter periods, and key factors, such as the animal diet, via its direct effect on the fatty acid profile, the antioxidant and volatile content, and indirect effects mediated via the age of the animal. It also pinpoints methodological difficulties in predicting organoleptic attributes, particularly odour and flavour. Potential antagonisms between different dimensions of quality are highlighted. For example, pasture-feeding has positive effects on the image and nutritional attributes (through its effect on the fatty acid profile of meat lipids), but it increases the risk of off-odours and off-flavours for sensitive consumersand the variability in meat quality linked to variability of animal age at slaughter. The orientation towards more agro-ecological, low-input farming systems may therefore present benefits for the image and nutritional properties of the meat, but also risks for the commercial (insufficient carcass fatness, feed deficiencies at key periods of the production cycle, irregularity in supply), organoleptic (stronger flavour and darker colour of the meat) and variability of sheep carcass and meat quality. Furthermore, the genetic selection for lean meat yield has been effective in producing carcasses that yield more meat, but at a penalty to the intramuscular fat content and eating quality of the meat, and making it more difficult to finish lambs on grass. Various tools to assess and predict quality are in development to better consider the various dimensions of quality in consumer information, payment to farmers and genetic selection.

Development and validation of the Standard method EN ISO 19020 - microbiology of the food chain - Horizontal method for the immunoenzymatic detection of staphylococcal enterotoxins in foodstuffs

In the frame of the CEN Mandate M/381 from the European Commission to CEN (European Committee for Standardization), a method for the detection of staphylococcal enterotoxins in foodstuffs has been developed, validated and standardized. An extraction procedure based on dialysis concentration followed by an immuno-enzymatic detection has been defined. In addition, performance criteria (minimum values of sensitivity, specificity and level of detection) to be achieved by the commercially available immuno-enzymatic kits that could be used to detect staphylococcal enterotoxins in food matrices, were developed. A 2-stage validation study was conducted: The first stage aimed at selecting the commercial kits to be included in the second stage, which consisted in an interlaboratory study, using eight matrices covering five food categories (ready-to-eat food, meat products, milk products, dessert and fish). Results showed that two detection kits included in the study met the pre-defined performance criteria. The implementation of dialysis concentration step increased significantly the sensitivity of the method. The method developed allowed to achieve the Benchmark Dose lower limit (BMD₁₀) estimated at 6.1 ng. In 2019, finally, the European Commission recognized this standard as the European Union reference method for the detection of staphylococcal enterotoxins in food.

A microscopy-based approach for determining growth probability and lag time of individual bacterial cells

The development of relevant predictive models for single-cell lag time and growth probability near growth limits is of critical importance for predicting pathogen behavior in foods. The classical methods for data acquisition in this field are based on turbidity measurements of culture media in microplate wells inoculated with approximately one bacterial cell per well. Yet, these methods are labour intensive and would benefit from higher throughput. In this study, we developed a quantitative experimental method using automated microscopy to determine the single-cell growth probability and lag time. The developed method consists of the use of direct cell observation with phase-contrast microscopy equipped with a 100× objective and a high-resolution device camera. The method is not a time-lapse method but is based on the observation of high numbers of colonies for a given time. Automation of image acquisition and image analysis was used to reach a high throughput. The single-cell growth probabilities and lag times of four strains of Listeria monocytogenes were determined at 4 °C. The microscopic method was shown to be a promising method for the determination of individual lag times and growth probability at the single-cell level.

Modeling the Inactivation of Viruses from the Coronaviridae Family in Response to Temperature and Relative Humidity in Suspensions or on Surfaces

Temperature and relative humidity are major factors determining virus inactivation in the environment. This article reviews inactivation data regarding coronaviruses on surfaces and in liquids from published studies and develops secondary models to predict coronaviruses inactivation as a function of temperature and relative humidity. A total of 102 D values (i.e., the time to obtain a log₁₀ reduction of virus infectivity), including values for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), were collected from 26 published studies. The values obtained from the different coronaviruses and studies were found to be generally consistent. Five different models were fitted to the global data set of D values. The most appropriate model considered temperature and relative humidity. A spreadsheet predicting the inactivation of coronaviruses and the associated uncertainty is presented and can be used to predict virus inactivation for untested temperatures, time points, or any coronavirus strains belonging to Alphacoronavirus and Betacoronavirus genera.IMPORTANCE The prediction of the persistence of SARS-CoV-2 on fomites is essential in investigating the importance of contact transmission. This study collects available information on inactivation kinetics of coronaviruses in both solid and liquid fomites and creates a mathematical model for the impact of temperature and relative humidity on virus persistence. The predictions of the model can support more robust decision-making and could be useful in various public health contexts. A calculator for the natural clearance of SARS-CoV-2 depending on temperature and relative humidity could be a valuable operational tool for public authorities.

Modeling the Inactivation of Viruses from the Coronaviridae Family in Response to Temperature and Relative Humidity in Suspensions or on Surfaces

Temperature and relative humidity are major factors determining virus inactivation in the environment. This article reviews inactivation data regarding coronaviruses on surfaces and in liquids from published studies and develops secondary models to predict coronaviruses inactivation as a function of temperature and relative humidity. A total of 102 D values (i.e., the time to obtain a log10 reduction of virus infectivity), including values for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), were collected from 26 published studies. The values obtained from the different coronaviruses and studies were found to be generally consistent. Five different models were fitted to the global data set of D values. The most appropriate model considered temperature and relative humidity. A spreadsheet predicting the inactivation of coronaviruses and the associated uncertainty is presented and can be used to predict virus inactivation for untested temperatures, time points, or any coronavirus strains belonging to Alphacoronavirus and Betacoronavirus genera.IMPORTANCE The prediction of the persistence of SARS-CoV-2 on fomites is essential in investigating the importance of contact transmission. This study collects available information on inactivation kinetics of coronaviruses in both solid and liquid fomites and creates a mathematical model for the impact of temperature and relative humidity on virus persistence. The predictions of the model can support more robust decision-making and could be useful in various public health contexts. A calculator for the natural clearance of SARS-CoV-2 depending on temperature and relative humidity could be a valuable operational tool for public authorities.

Rapid assessment and prediction of the efficiency of two preservatives against S. aureus in cosmetic products using High Content Screening-Confocal Laser Scanning Microscopy

Most cosmetic products are susceptible to microbiological spoilage due to contaminations that could happen during fabrication or by consumer’s repetitive manipulation. The composition of cosmetic products must guarantee efficient bacterial inactivation all along with the product shelf life, which is usually assessed by challenge-tests. A challenge-test consists in inoculating specific bacteria, i.e. Staphylococcus aureus, in the formula and then investigating the bacterial log reduction over time. The main limitation of this method is relative to the time-consuming protocol, where 30 days are needed to obtain results. In this study, we have proposed a rapid alternative method coupling High Content Screening—Confocal Laser Scanning Microscopy (HCS-CLSM), image analysis and modeling. It consists in acquiring real-time S. aureus inactivation kinetics on short-time periods (typically 4h) and in predicting the efficiency of preservatives on longer scale periods (up to 7 days). The action of two preservatives, chlorphenesin and benzyl alcohol, was evaluated against S. aureus at several concentrations in a cosmetic matrix. From these datasets, we compared two secondary models to determine the logarithm reduction time (Dc) for each preservative concentration. Afterwards, we used two primary inactivation models to predict log reductions for up to 7 days and we compared them to observed log reductions. The IQ model better fits datasets and the Q value gives information about the matrix level of interference.

AB_SA: Accessory genes-Based Source Attribution - tracing the source of Salmonella enterica Typhimurium environmental strains

The partitioning of pathogenic strains isolated in environmental or human cases to their sources is challenging. The pathogens usually colonize multiple animal hosts, including livestock, which contaminate the food-production chain and the environment (e.g. soil and water), posing an additional public-health burden and major challenges in the identification of the source. Genomic data opens up new opportunities for the development of statistical models aiming to indicate the likely source of pathogen contamination. Here, we propose a computationally fast and efficient multinomial logistic regression source-attribution classifier to predict the animal source of bacterial isolates based on 'source-enriched' loci extracted from the accessory-genome profiles of a pangenomic dataset. Depending on the accuracy of the model's self-attribution step, the modeller selects the number of candidate accessory genes that best fit the model for calculating the likelihood of (source) category membership. The Accessory genes-Based Source Attribution (AB_SA) method was applied to a dataset of strains of Salmonella enterica Typhimurium and its monophasic variant (S. enterica 1,4,[5],12:i:-). The model was trained on 69 strains with known animal-source categories (i.e. poultry, ruminant and pig). The AB_SA method helped to identify 8 genes as predictors among the 2802 accessory genes. The self-attribution accuracy was 80 %. The AB_SA model was then able to classify 25 of the 29 S. enterica Typhimurium and S. enterica 1,4,[5],12:i:- isolates collected from the environment (considered to be of unknown source) into a specific category (i.e. animal source), with more than 85 % of probability. The AB_SA method herein described provides a user-friendly and valuable tool for performing source-attribution studies in only a few steps. AB_SA is written in R and freely available at https://github.com/lguillier/AB_SA.

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.

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

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

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.

Biofilm Formation of Listeria monocytogenes Strains Under Food Processing Environments and Pan-Genome-Wide Association Study

Concerns about food contamination by Listeria monocytogenes are on the rise with increasing consumption of ready-to-eat foods. Biofilm production of L. monocytogenes is presumed to be one of the ways that confer its increased resistance and persistence in the food chain. In this study, a collection of isolates from foods and food processing environments (FPEs) representing persistent, prevalent, and rarely detected genotypes was evaluated for biofilm forming capacities including adhesion and sessile biomass production under diverse environmental conditions. The quantity of sessile biomass varied according to growth conditions, lineage, serotype as well as genotype but association of clonal complex (CC) 26 genotype with biofilm production was evidenced under cold temperature. In general, relative biofilm productivity of each strain varied inconsistently across growth conditions. Under our experimental conditions, there were no clear associations between biofilm formation efficiency and persistent or prevalent genotypes. Distinct extrinsic factors affected specific steps of biofilm formation. Sudden nutrient deprivation enhanced cellular adhesion while a prolonged nutrient deficiency impeded biofilm maturation. Salt addition increased biofilm production, moreover, nutrient limitation supplemented by salt significantly stimulated biofilm formation. Pan-genome-wide association study (Pan-GWAS) assessed genetic composition with regard to biofilm phenotypes for the first time. The number of reported genes differed depending on the growth conditions and the number of common genes was low. However, a broad overview of the ontology contents revealed similar patterns regardless of the conditions. Functional analysis showed that functions related to transformation/competence and surface proteins including Internalins were highly enriched.

Critical Orientation in the Jungle of Currently Available Methods and Types of Data for Source Attribution of Foodborne Diseases

With increased interest in source attribution of foodborne pathogens, there is a need to sort and assess the applicability of currently available methods. Herewith we reviewed the most frequently applied methods for source attribution of foodborne diseases, discussing their main strengths and weaknesses to be considered when choosing the most appropriate methods based on the type, quality, and quantity of data available, the research questions to be addressed, and the (epidemiological and microbiological) characteristics of the pathogens in question. A variety of source attribution approaches have been applied in recent years. These methods can be defined as top–down, bottom–up, or combined. Top–down approaches assign the human cases back to their sources of infection based on epidemiological (e.g., outbreak data analysis, case-control/cohort studies, etc.), microbiological (i.e., microbial subtyping), or combined (e.g., the so-called ‘source-assigned case-control study’ design) methods. Methods based on microbial subtyping are further differentiable according to the modeling framework adopted as frequency-matching (e.g., the Dutch and Danish models) or population genetics (e.g., Asymmetric Island Models and STRUCTURE) models, relying on the modeling of either phenotyping or genotyping data of pathogen strains from human cases and putative sources. Conversely, bottom–up approaches like comparative exposure assessment start from the level of contamination (prevalence and concentration) of a given pathogen in each source, and then go upwards in the transmission chain incorporating factors related to human exposure to these sources and dose-response relationships. Other approaches are intervention studies, including ‘natural experiments,’ and expert elicitations. A number of methodological challenges concerning all these approaches are discussed. In absence of an universally agreed upon ‘gold’ standard, i.e., a single method that satisfies all situations and needs for all pathogens, combining different approaches or applying them in a comparative fashion seems to be a promising way forward.

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.

Combining Quantitative Risk Assessment of Human Health, Food Waste, and Energy Consumption: The Next Step in the Development of the Food Cold Chain?

The preservation of perishable food via refrigeration in the supply chain is essential to extend shelf life and provide consumers with safe food. However, electricity consumed in refrigeration processes has an economical and an environmental impact. This study focuses on the cold chain of cooked ham, including transport, cold room in supermarket, display cabinet, transport by consumer, and domestic refrigerator, and aims to predict the risk for human health associated with Listeria monocytogenes, the amount of food wasted due to the growth of spoilage bacteria, and the electrical consumption to maintain product temperature through the cold chain. A set of eight intervention actions were tested to evaluate their impact on the three criteria. Results show that the modification of the thermostat of the domestic refrigerator has a high impact on food safety and food waste and a limited impact on the electrical consumption. Inversely, the modification of the airflow rate in the display cabinet has a high impact on electrical consumption and a limited impact on food safety and food waste. A cost-benefit analysis approach and two multicriteria decision analysis methods were used to rank the intervention actions. These three methodologies show that setting the thermostat of the domestic refrigerator to 4 °C presents the best compromise between the three criteria. The impact of decisionmaker preferences (criteria weight) and limitations of these three approaches are discussed. The approaches proposed by this study may be useful in decision making to evaluate global impact of intervention actions in issues involving conflicting outputs.

LiSEQ - whole-genome sequencing of a cross-sectional survey of Listeria monocytogenes in ready-to-eat foods and human clinical cases in Europe

We present the LiSEQ (Listeria SEQuencing) project, funded by the European Food Safety Agency (EFSA) to compare Listeria monocytogenes isolates collected in the European Union from ready-to-eat foods, compartments along the food chain (e.g. food-producing animals, food-processing environments) and humans. In this article, we report the molecular characterization of a selection of this data set employing whole-genome sequencing analysis. We present an overview of the strain diversity observed in different sampled sources, and characterize the isolates based on their virulence and resistance profile. We integrate into our analysis the global L. monocytogenes genome collection described by Moura and colleagues in 2016 to assess the representativeness of the LiSEQ collection in the context of known L. monocytogenes strain diversity.

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.

Development of a Real-Time Cell Analysis (RTCA) Method as a Fast and Accurate Method for Detecting Infectious Particles of the Adapted Strain of Hepatitis A Virus

Hepatitis A virus (HAV) is one of the most common agents causing acute liver disease worldwide. HAV has been increasingly reported as the cause of foodborne disease outbreaks. The standard method currently available for detection of the genome of HAV in vulnerable foodstuffs is by RT-qPCR (ISO 15216). Despite its usefulness in the investigation of foodborne viruses, the use of RT-qPCR in food virology has been shown to overestimate the quantity of infectious virus or to highly underestimate the effect of the treatment on virus inactivation. The gold standard methods currently used for evaluating the efficacy of inactivation treatments on the adapted strain of HAV (HM175/18f) are either the plaque assay or the end-point dilution assay (TCID₅₀). However, both assays are labor-intensive and time-consuming. The aim of this study was to evaluate the use of the xCELLigence real-time cell analysis (RTCA) system for detecting the infectivity of the adapted strain of HAV. Kinetics of cell impedance showed that HAV induced a decrease in cell index (CI) correlated with the onset of HAV-induced cell death. In addition, the time to which the HAV-induced CI drop occurred was dependent on the viral concentration. An inverse linear relation could be established over a range of 5 log₁₀ between the concentration of HAV and the time to reach 50% of CI decrease (TCI₅₀), showing that the RTCA assay could be used as a titration method for HAV. In addition, the RTCA-based assay could be performed in less than 6 days instead of 12 to 14 days with the gold standard methods. Therefore, the RTCA-based titration method is a powerful and suitable tool for high-throughput screening of anti-viral treatments. Its usefulness in HAV inactivation studies will improve the assessment of viral risk in food virology, as controlling transmission of viruses through their removal from foodstuffs is also an important challenge in reducing the burden of viral foodborne illnesses.

Survival of Cronobacter in powdered infant formula and their variation in biofilm formation

Cronobacter is a ubiquitous Gram-negative pathogen bacterium capable of surviving in low water activity environments, in particular powdered infant formula (PIF). Seven Cronobacter strains representing four different species (C. sakazakii, n = 4; C. malonaticus, n = 1; C. muytjensii, n = 1; C. turicensis, n = 1) were subjected to dry stress and stored in PIF at room temperature. The resulting survivor curves showed that Cronobacter sp. can survive for extended periods of at least 3 months with a significant, but moderate, variability regarding the level of resistance between species; however, no correlation was evident regarding the origin of strains. These results are evaluated with regard to other key characteristics, including genomic profiles and biofilm formation capacities of the strains.

SIGNIFICANCE AND IMPACT OF THE STUDY

Cronobacter can survive extended periods of at least 3 months in PIF, with moderately significant interspecific variability in desiccation resistance. Results are evaluated with regard to genomic profiles and biofilm formation capacities of the strains, and contribute to an improved understanding of the environmental persistence of Cronobacter in contaminated PIF, and subsequent risk to infant exposure.

Population Genetic Structure of Listeria monocytogenes Strains Isolated From the Pig and Pork Production Chain in France

Listeria monocytogenes is an ubiquitous pathogenic bacterium, transmissible to humans through the consumption of contaminated food. The pork production sector has been hit hard by a series of L. monocytogenes-related food poisoning outbreaks in France. An overview of the diversity of strains circulating at all levels of the pork production chain, from pig farming (PF) to finished food products (FFP), is needed to identify the contamination routes and improve food safety. Until now, no typing data has been available on strains isolated across the entire pig and pork production chain. Here, we analyzed the population genetic structure of 687 L. monocytogenes strains isolated over the last 20 years in virtually all the French départements from three compartments of this production sector: PF, the food processing environment (FPE), and FFP. The genetic structure was described based on Multilocus sequence typing (MLST) clonal complexes (CCs). The CCs were obtained by mapping the PFGE profiles of the strains. The distribution of CCs was compared firstly between the three compartments and then with CCs obtained from 1106 strains isolated from other food production sectors in France. The predominant CCs of pig and pork strains were not equally distributed among the three compartments: the CC37, CC59, and CC77 strains, rarely found in FPE and FFP, were prevalent in PF. The two most prevalent CCs in the FPE and FFP compartments, CC9 and CC121, were rarely or never detected in PF. No CC was exclusively associated with the pork sector. Three CCs (CC5, CC6, and CC2) were considered ubiquitous, because they were observed in comparable proportions in all food production sectors. The two most prevalent CCs in all sectors were CC9 and CC121, but their distribution was disparate. CC9 was associated with meat products and food products combining several food categories, whereas CC121 was not associated with any given sector. Based on these results, CC121 is likely able to colonize a larger diversity of food products than CC9. Both CCs being associated with the food production suggests, that certain processing steps, such as slaughtering or stabilization treatments, favor their settlement and the recontamination of the food produced.

Listeria monocytogenes Sequence Types 121 and 14 Repeatedly Isolated Within One Year of Sampling in a Rabbit Meat Processing Plant: Persistence and Ecophysiology

Listeria monocytogenes is a foodborne pathogen adapted to survive and persist in multiple environments. Following two previous studies on prevalence and virulence of L. monocytogenes ST121 and ST14 repeatedly collected in a the same rabbit-meat processing plant, the research questions of the present study were to: (1) assess persistence of L. monocytogenes isolates from the rabbit-plant; (2) select genes associated to physiological adaptation to the food-processing environment; (3) compare presence/absence/truncation of these genes in newly sequenced and publicly available ST121 and ST14 genomes. A total of 273 draft genomes including ST121 and ST14 newly sequenced and publicly available draft genomes were analyzed. Whole-genome Single Nucleotide Polymorfism (wgSNP) analysis was performed separately on the assemblies of ST121 and ST14 draft genomes. SNPs alignments were used to infer phylogeny. A dataset of L. monocytogenes ecophysiology genes was built based on a comprehensive literature review. The 94 selected genes were screened on the assemblies of all ST121 and ST14 draft genomes. Significant gene enrichments were evaluated by statistical analyses. A persistent ST14 clone, including 23 out of 27 newly sequenced genomes, was circulating in the rabbit-meat plant along with two not persistent clones. A significant enrichment was observed in ST121 genomes concerning stress survival islet 2 (SSI-2) (alkaline and oxidative stress), qacH gene (resistance to benzalkonium chloride), cadA1C gene cassette (resistance to 70 mg/l of cadmium chloride) and a truncated version of actA gene (biofilm formation). Conversely, ST14 draft genomes were enriched with a full-length version of actA gene along with the Listeria Genomic Island 2 (LGI 2) including the ars operon (arsenic resistance) and the cadA4C gene cassette (resistance to 35 mg/l of cadmium chloride). Phenotypic tests confirmed ST121 as a weak biofilm producer in comparison to ST14. In conclusion, ST121 carried the qacH gene and was phenotypically resistant to quaternary ammonium compounds. This property might contribute to the high prevalence of ST121 in food processing plants. ST14 showed greater ability to form biofilms, which might contribute to the occasional colonization and persistence on harborage sites where sanitizing procedures are difficult to display.