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

Virulence characterization and comparative genomics of Listeria monocytogenes sequence type 155 strains

BACKGROUND

Listeria (L.) monocytogenes strains show a high diversity regarding stress tolerance and virulence potential. Genome studies have mainly focused on specific sequence types (STs) predominantly associated with either food or human listeriosis. This study focused on the prevalent ST155, showing equal distribution among clinical and food isolates. We evaluated the virulence potential of 20 ST155 strains and performed comparative genomic analysis of 130 ST155 strains isolated from food, food processing environments and human listeriosis cases in different countries and years.

RESULTS

The in vitro virulence assays using human intestinal epithelial Caco2 and hepatocytic HEPG2 cells showed an impaired virulence phenotype for six of the 20 selected ST155 strains. Genome analysis revealed no distinct clustering of strains from the same source category (food, food processing environment, and clinical isolates). All strains harbored an intact inlA and inlB locus, except four strains, which had an internal deletion in the inlA gene. All strains harbored LIPI-1, but prfA was present in a longer variant in six strains, all showing impaired virulence. The longer PrfA variant resulted in lower expression of inlA, inlB, and prfA, and no expression of hly and actA. Regarding stress-related gene content, SSI-1 was present, whereas qacH was absent in all strains. 34.6% of the strains harbored a plasmid. All but one ST155 plasmids showed high conservation and harbored cadA2, bcrABC, and a triphenylmethane reductase.

CONCLUSIONS

This study contributes to an enhanced understanding of L. monocytogenes ST155 strains, being equally distributed among isolates from humans, food, and food processing environments. The conservation of the present genetic traits and the absence of unique inherent genetic features makes these types of STs especially interesting since they are apparently equally adapted to the conditions in food processing environments, as well as in food as to the human host environment. However, a ST155-specific mutation resulting in a longer PrfA variant impaired the virulence potential of several ST155 strains.

An outbreak of human listeriosis associated with frozen sweet corn consumption: Investigations in the UK

The use of Whole genome sequencing (WGS) identified a multi-country outbreak of human listeriosis associated with consumption of frozen sweet corn produced in Hungary. The purpose of this report was to summarise information on the cases occurring in the UK which were part of this outbreak and outline investigations on the presence of Listeria monocytogenes in the affected food chain. Prior to the international recall of this product in 2018, 12 UK cases of listeriosis were identified as infected by the outbreak strain between 2015 and 18. Epidemiological and microbiological investigations confirmed these cases as belonging to the outbreak. A further case occurred in 2019 and a contaminated frozen pack from one of the implicated batches of sweet corn was recovered from the patient's domestic freezer. The outbreak strain was also detected in products from a sandwich manufacturer in 2018 which added frozen sweet corn directly to sandwich fillings. The sandwich manufacturer's sweet corn was supplied by a distributor in England which obtained frozen products from the Hungarian manufacturer implicated in the outbreak. Within the distributor's premises, 208 food and environmental samples were taken: L. monocytogenes was detected in 44% of 70 samples of frozen sweet corn and 5% of 79 other foods. The outbreak strain was detected in the frozen sweet corn, in one other frozen food (mixed vegetables) and in the factory environment. The outbreak strain was also recovered from frozen beans on retail sale in the first four months of 2019. Five other L. monocytogenes strains together with two other Listeria species were detected in samples from the importer's premises. One of the L. monocytogenes strains in the importer's factory, which was distinct from the outbreak strain, was also recovered from sweet corn collected from the sandwich manufacturer, sweet corn tested in England in 2013 and 2016 and the blood of two cases of human listeriosis which occurred in England in 2014. This report shows how analysis by WGS provides evidence to understand complex food chains. This report also highlights risks for transmission of human listeriosis from frozen sweet corn and the potential for misuse of this food as a ready-to-eat product.

First Report on the Finding of Listeria mnocytogenes ST121 Strain in a Dolphin Brain

Listeria monocytogenes (Lm) is a ubiquitous bacterium that causes the foodborne illness, listeriosis. Clonal complexes (CC), such as CC121, are overrepresented in the food production industry, and are rarely reported in animals and the environment. Working within a European-wide project, we investigated the routes by which strains are transmitted from environments and animals to food and the food production environment (FPE). In this context, we report, for the first time, the occurrence of a ST121 (CC121) strain isolated from a dolphin brain. The genome was compared with the genomes of 376 CC121 strains. Genomic comparisons showed that 16 strains isolated from food were the closest to the dolphin strain. Like most of the food strains analyzed here, the dolphin strain included genomic features (transposon Tn6188, plasmid pLM6179), both described as being associated with the strain’s adaptation to the FPE. Like all 376 strains, the dolphin strain contained a truncated actA gene and inlA gene, both described as being associated with attenuated virulence. Despite this fact, the strain was able to cross blood-brain barrier in immunosuppressed dolphin exposed polychlorinated biphenyl and invaded by parasites. Our data suggest that the dolphin was infected by a food-related strain released into the Mediterranean Sea.

In-Depth Longitudinal Study of Listeria monocytogenes ST9 Isolates from the Meat Processing Industry: Resolving Diversity and Transmission Patterns Using Whole-Genome Sequencing

Listeria monocytogenes is a deadly foodborne pathogen that is widespread in the environment, and certain types can be established in food factories. The sequence type ST9 dominates in meat processing environments, and this work was undertaken to obtain data needed for the tracking of this subtype. By using whole-genome sequencing (WGS), we revealed the presence of cross-contamination routes between meat factories as well as within a single factory, including the spread from different reservoirs within the same room. It was also possible to estimate the time frame of persistence in the factory, as well as when and how new clones had entered. The present work contributes valuable information about the diversity of ST9 and exemplifies the potential power of WGS in food safety management, allowing the determination of relationships between strains both in an international context and locally between and within factories.

Listeria monocytogenes is a pathogen mostly associated with the consumption of ready-to-eat foods and can cause severe disease and death. It can be introduced into food chains from raw materials, but often the contamination source is the food production environment, where certain clones can persist for years. In the meat chain, ST9 is one of the most commonly encountered L. monocytogenes sequence types, and for effective source tracking, the divergence and spread of ST9 must be understood. In this study, whole-genome sequencing (WGS) was used to characterize and track 252 L. monocytogenes ST9 isolates collected from four Norwegian meat processing plants between 2009 and 2017. The isolates formed distinct clusters relative to genomes found in public databases, and all but three isolates clustered into two major clonal populations. Different contamination patterns were revealed, e.g., evidence of contamination of two factories with a clone that diverged from its ancestor in the late 1990s through a common source of raw materials; breach of hygienic barriers within a factory, leading to repeated detection of two clones in the high-risk zone during a 4- to 6-year period; entry through the purchase and installation of second-hand equipment harboring a previously established clonal population; and spreading and diversification of two clones from two reservoirs within the same production room over a 9-year period. The present work provides data on the diversity of ST9, which is crucial for epidemiological investigations and highlights how WGS can be used for source tracking within food processing factories.

IMPORTANCEListeria monocytogenes is a deadly foodborne pathogen that is widespread in the environment, and certain types can be established in food factories. The sequence type ST9 dominates in meat processing environments, and this work was undertaken to obtain data needed for the tracking of this subtype. By using whole-genome sequencing (WGS), we revealed the presence of cross-contamination routes between meat factories as well as within a single factory, including the spread from different reservoirs within the same room. It was also possible to estimate the time frame of persistence in the factory, as well as when and how new clones had entered. The present work contributes valuable information about the diversity of ST9 and exemplifies the potential power of WGS in food safety management, allowing the determination of relationships between strains both in an international context and locally between and within factories.

Combining Whole-Genome Sequencing and Multimodel Phenotyping To Identify Genetic Predictors of Salmonella Virulence

Salmonella comprises more than 2,600 serovars. Very few environmental and uncommon serovars have been characterized for their potential role in virulence and human infections. A complementary in vitro and in vivo systematic high-throughput analysis of virulence was used to elucidate the association between genetic and phenotypic variations across Salmonella isolates. The goal was to develop a strategy for the classification of isolates as a benchmark and predict virulence levels of isolates. Thirty-five phylogenetically distant strains of unknown virulence were selected from the Salmonella Foodborne Syst-OMICS (SalFoS) collection, representing 34 different serovars isolated from various sources. Isolates were evaluated for virulence in 4 complementary models of infection to compare virulence traits with the genomics data, including interactions with human intestinal epithelial cells, human macrophages, and amoeba. In vivo testing was conducted using the mouse model of Salmonella systemic infection. Significant correlations were identified between the different models. We identified a collection of novel hypothetical and conserved proteins associated with isolates that generate a high burden. We also showed that blind prediction of virulence of 33 additional strains based on the pan-genome was high in the mouse model of systemic infection (82% agreement) and in the human epithelial cell model (74% agreement). These complementary approaches enabled us to define virulence potential in different isolates and present a novel strategy for risk assessment of specific strains and for better monitoring and source tracking during outbreaks.IMPORTANCE Salmonella species are bacteria that are a major source of foodborne disease through contamination of a diversity of foods, including meat, eggs, fruits, nuts, and vegetables. More than 2,600 different Salmonella enterica serovars have been identified, and only a few of them are associated with illness in humans. Despite the fact that they are genetically closely related, there is enormous variation in the virulence of different isolates of Salmonella enterica Identification of foodborne pathogens is a lengthy process based on microbiological, biochemical, and immunological methods. Here, we worked toward new ways of integrating whole-genome sequencing (WGS) approaches into food safety practices. We used WGS to build associations between virulence and genetic diversity within 83 Salmonella isolates representing 77 different Salmonella serovars. Our work demonstrates the potential of combining a genomics approach and virulence tests to improve the diagnostics and assess risk of human illness associated with specific Salmonella isolates.

Predominance of Distinct Listeria Innocua and Listeria Monocytogenes in Recurrent Contamination Events at Dairy Processing Facilities

: The genus Listeria now comprises up to now 21 recognized species and six subspecies, with L. monocytogenes and L. innocua as the most prevalent sensu stricto associated species. Reports focusing on the challenges in Listeria detection and confirmation are available, especially from food-associated environmental samples. L. innocua is more prevalent in the food processing environment (FPE) than L. monocytogenes and has been shown to have a growth advantage in selective enrichment and agar media. Until now, the adaptive nature of L. innocua in FPEs has not been fully elucidated and potential persistence in the FPE has not been observed. Therefore, the aim of this study is to characterize L. innocua (n = 139) and L. monocytogenes (n = 81) isolated from FPEs and cheese products collected at five dairy processing facilities (A-E) at geno- and phenotypic levels. Biochemical profiling was conducted for all L. monocytogenes and the majority of L. innocua (n = 124) isolates and included a rhamnose positive reaction. L. monocytogenes isolates were most frequently confirmed as PCR-serogroups 1/2a, 3a (95%). Pulsed-field gel electrophoresis (PFGE)-typing, applying the restriction enzymes AscI, revealed 33 distinct Listeria PFGE profiles with a Simpson's Index of Diversity of 0.75. Multi-locus sequence typing (MLST) resulted in 27 STs with seven new L. innocua local STs (ST1595 to ST1601). L. innocua ST1597 and ST603 and L. monocytogenes ST121 and ST14 were the most abundant genotypes in dairy processing facilities A-E over time. Either SSI-1 (ST14) or SSI-2 (ST121, all L. innocua) were present in successfully FPE-adapted strains. We identified housekeeping genes common in Listeria isolates and L. monocytogenes genetic lineage III. Wherever there are long-term contamination events of L. monocytogenes and other Listeria species, subtyping methods are helpful tools to identify niches of high risk.

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.

Interference of the acid stress on the expression of llsX by Listeria monocytogenes pathogenic island 3 (LIPI-3) variants

Listeria monocytogenes harbor different virulence factors, with a highly heterogeneous distribution between distinct lineages and serotypes. The Listeria Pathogenicity Island 3 (LIPI-3), mainly described in lineage I, encodes for Listeriolysin S (LLS), a virulence factor expressed by L. monocytogenes in the gastrointestinal tract during in vivo infections. The aim of this study was to carry out a comparative genotypic analysis of LIPI-3 identified in L. monocytogenes isolates obtained in Brazil and subjected to whole genomic sequencing (WGS). In addition, transcription of llsX expression under different acid stress conditions was evaluated by RT-PCR. Homologues of the eight LIPI-3 genes (llsAGHXBYDP) were identified in 15 isolates (all from lineage I) representative of different sequence types: ST1 (n = 3), ST3 (n = 6), ST218 (n = 5) and ST288 (n = 1). Single nucleotide polymorphism (SNP) analysis revealed that genetic variation resulted in modification of the final peptide LLS for ST218 (serogroup IVb-v1) and ST288 (serogroup IIb). Selected strains from ST3 and ST288 were subjected to acid stress conditions and the expression of llsX, a LIPI-3 gene, was observed: only F2365 (4b/ST1) presented llsX expression after six hours of acid stress, indicating relevant differences when compared to isolates IIb (ST3 and 288). The results highlight the presence of genomic variations on LIPI-3 and llsX expression under acid stress conditions, demanding further studies to evaluate if these mutations have an impact on L. monocytogenes virulence in vivo.

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.

Exploring Listeria monocytogenes Transcriptomes in Correlation with Divergence of Lineages and Virulence as Measured in Galleria mellonella

As for many opportunistic pathogens, the virulence potential of Listeria monocytogenes is highly heterogeneous between isolates and correlated, to some extent, with phylogeny and gene repertoires. In sharp contrast with copious data on intraspecies genome diversity, little is known about transcriptome diversity despite the role of complex genetic regulation in pathogenicity. The current study implemented RNA sequencing to characterize the transcriptome profiles of 33 isolates under optimal in vitro growth conditions. Transcript levels of conserved single-copy genes were comprehensively explored from several perspectives, including phylogeny, in silico-predicted virulence category based on epidemiological multilocus sequence typing (MLST) data, and in vivo virulence phenotype assessed in Galleria mellonella Comparing baseline transcriptomes between isolates was intrinsically more complex than standard genome comparison because of the inherent plasticity of gene expression in response to environmental conditions. We show that the relevance of correlation analyses and their statistical power can be enhanced by using principal-component analysis to remove the first level of irrelevant, highly coordinated changes linked to growth phase. Our results highlight the major contribution of transcription factors with key roles in virulence to the diversity of transcriptomes. Divergence in the basal transcript levels of a substantial fraction of the transcriptome was observed between lineages I and II, echoing previously reported epidemiological differences. Correlation analysis with in vivo virulence identified numerous sugar metabolism-related genes, suggesting that specific pathways might play roles in the onset of infection in G. mellonella IMPORTANCE Listeria monocytogenes is a multifaceted bacterium able to proliferate in a wide range of environments from soil to mammalian host cells. The accumulated genomic data underscore the contribution of intraspecies variations in gene repertoire to differential adaptation strategies between strains, including infection and stress resistance. It seems very likely that the fine-tuning of the transcriptional regulatory network is also a key component of the phenotypic diversity, albeit more difficult to investigate than genome content. Some studies reported incongruity in the basal transcriptome between isolates, suggesting a putative relationship with phenotypes, but small isolate numbers hampered proper correlation analyses with respect to their characteristics. The present study is the embodiment of the promising approach that consists of analyzing correlations between transcriptomes and various isolate characteristics. Statistically significant correlations were found with phylogenetic groups, epidemiological evidence of virulence potential, and virulence in Galleria mellonella larvae used as an in vivo model.