Microevolution and Gain or Loss of Mobile Genetic Elements of Outbreak-Related Listeria monocytogenes in Food Processing Environments Identified by Whole Genome Sequencing Analysis

Persistence of microbiological hazards in food and feed production and processing environments

Listeria monocytogenes (in the meat, fish and seafood, dairy and fruit and vegetable sectors), Salmonella enterica (in the feed, meat, egg and low moisture food sectors) and Cronobacter sakazakii (in the low moisture food sector) were identified as the bacterial food safety hazards most relevant to public health that are associated with persistence in the food and feed processing environment (FFPE). There is a wide range of subtypes of these hazards involved in persistence in the FFPE. While some specific subtypes are more commonly reported as persistent, it is currently not possible to identify universal markers (i.e. genetic determinants) for this trait. Common risk factors for persistence in the FFPE are inadequate zoning and hygiene barriers; lack of hygienic design of equipment and machines; and inadequate cleaning and disinfection. A well-designed environmental sampling and testing programme is the most effective strategy to identify contamination sources and detect potentially persistent hazards. The establishment of hygienic barriers and measures within the food safety management system, during implementation of hazard analysis and critical control points, is key to prevent and/or control bacterial persistence in the FFPE. Once persistence is suspected in a plant, a 'seek-and-destroy' approach is frequently recommended, including intensified monitoring, the introduction of control measures and the continuation of the intensified monitoring. Successful actions triggered by persistence of L. monocytogenes are described, as well as interventions with direct bactericidal activity. These interventions could be efficient if properly validated, correctly applied and verified under industrial conditions. Perspectives are provided for performing a risk assessment for relevant combinations of hazard and food sector to assess the relative public health risk that can be associated with persistence, based on bottom-up and top-down approaches. Knowledge gaps related to bacterial food safety hazards associated with persistence in the FFPE and priorities for future research are provided.

Tetracycline resistance in Listeria monocytogenes and L. innocua from wild black bears (Ursus americanus) in the United States is mediated by novel transposable elements

IMPORTANCE

Listeria monocytogenes causes severe foodborne illness and is the only human pathogen in the genus Listeria. Previous surveys of AMR in Listeria focused on clinical sources and food or food processing environments, with AMR in strains from wildlife and other natural ecosystems remaining under-explored. We analyzed 185 sequenced strains from wild black bears (Ursus americanus) from the United States, including 158 and 27 L. monocytogenes and L. innocua, respectively. Tetracycline resistance was the most prevalent resistance trait. In L. monocytogenes, it was encountered exclusively in serotype 4b strains with the novel Tn916-like element Tn916.1039. In contrast, three distinct, novel tetracycline resistance elements (Tn5801.UAM, Tn5801.551, and Tn6000.205) were identified in L. innocua. Interestingly, Tn5801.551 was identical to elements in L. monocytogenes from a major foodborne outbreak in the United States in 2011. The findings suggest the importance of wildlife and non-pathogenic Listeria species as reservoir for resistance elements in Listeria.

The Slaughterhouse as Hotspot of CC1 and CC6 Listeria monocytogenes Strains with Hypervirulent Profiles in an Integrated Poultry Chain of Italy

In Europe, very few studies are available regarding the diversity of Listeria monocytogenes (L. monocytogenes) clonal complexes (CCs) and sequence types (ST) in poultry and on the related typing of isolates using whole genome sequencing (WGS). In this study, we used a WGS approach to type 122 L. monocytogenes strains isolated from chicken neck skin samples collected in two different slaughterhouses of an integrated Italian poultry company. The studied strains were classified into five CCs: CC1-ST1 (21.3%), CC6-ST6 (22.9%), CC9-ST9 (44.2%), CC121-ST121 (10.6%) and CC193-ST193 (0.8%). CC1 and CC6 strains presented a virulence gene profile composed of 60 virulence genes and including the Listeria Pathogenicity Island 3, aut_IVb, gltA and gltB. According to cgMLST and SNPs analysis, long-term persistent clusters belonging to CC1 and CC6 were found in one of the two slaughterhouses. The reasons mediating the persistence of these CCs (up to 20 months) remain to be elucidated, and may involve the presence and the expression of stress response and environmental adaptation genes including heavy metals resistance genes (cadAC, arsBC, CsoR-copA-copZ), multidrug efflux pumps (mrpABCEF, EmrB, mepA, bmrA, bmr3, norm), cold-shock tolerance (cspD) and biofilm-formation determinants (lmo0673, lmo2504, luxS, recO). These findings indicated a serious risk of poultry finished products contamination with hypervirulent L. monocytogenes clones and raised concern for the consumer health. In addition to the AMR genes norB, mprF, lin and fosX, ubiquitous in L. monocytogenes strains, we also identified parC for quinolones, msrA for macrolides and tetA for tetracyclines. Although the phenotypical expression of these AMR genes was not tested, none of them is known to confer resistance to the primary antibiotics used to treat listeriosis The obtained results increase the data on the L. monocytogenes clones circulating in Italy and in particular in the poultry chain.

Investigation of a Listeria monocytogenes Chromosomal Immigration Control Region Reveals Diverse Restriction Modification Systems with Complete Sequence Type Conservation

Listeria monocytogenes is a Gram-positive pathogen responsible for the severe foodborne disease listeriosis. A chromosomal hotspot between lmo0301 and lmo0305 has been noted to harbor diverse restriction modification (RM) systems. Here, we analyzed 872 L. monocytogenes genomes to better understand the prevalence and types of RM systems in this region, designated the immigration control region (ICR). Type I, II, III and IV RM systems were found in 86.1% of strains inside the ICR and in 22.5% of strains flanking the ICR. ICR content was completely conserved within the same multilocus sequence typing-based sequence type (ST), but the same RM system could be identified in diverse STs. The intra-ST conservation of ICR content suggests that this region may drive the emergence of new STs and promote clone stability. Sau3AI-like, LmoJ2 and LmoJ3 type II RM systems as well as type I EcoKI-like, and type IV AspBHI-like and mcrB-like systems accounted for all RM systems in the ICR. A Sau3AI-like type II RM system with specificity for GATC was harbored in the ICR of many STs, including all strains of the ancient, ubiquitous ST1. The extreme paucity of GATC recognition sites in lytic phages may reflect ancient adaptation of these phages to preempt resistance associated with the widely distributed Sau3AI-like systems. These findings indicate that the ICR has a high propensity for RM systems which are intraclonaly conserved and may impact bacteriophage susceptibility as well as ST emergence and stability.

Horizontal Gene Transfer and Loss of Serotype-Specific Genes in Listeria monocytogenes Can Lead to Incorrect Serotype Designations with a Commonly-Employed Molecular Serotyping Scheme

Listeria monocytogenes is a Gram-positive, facultative intracellular foodborne pathogen capable of causing severe, invasive illness (listeriosis). Three serotypes, 1/2a, 1/2b, and 4b, are leading contributors to human listeriosis, with 4b including the major hypervirulent clones. The multiplex PCR scheme developed by Doumith and collaborators employs primers targeting specific lineages (e.g., lineage II-specific lmo0737, lineage I-specific LMOf2365_2059) or serotypes (e.g., serotype 4b-specific LMOf2365_1900). The Doumith scheme (DS) is extensively employed for molecular serotyping of L. monocytogenes due to its high accuracy, relative ease, and affordability. However, for certain strains, the DS serotype designations are in conflict with those relying on antibody-based schemes or whole-genome sequence (WGS) analysis. In the current study, all 27 tested serotype 4b strains with sequence type 782 (ST782) within the hypervirulent clonal complex 2 (CC2) were designated 1/2b/3b using the DS. These strains lacked the serotype 4b-specific gene LMOf2365_1900, while retaining LMOf2365_2059, which, together with prs, yields the DS 1/2b/3b profile. Furthermore, 15 serotype 1/2a strains of four STs, mostly from water, were designated 1/2b/3b using the DS. These strains lacked the lmo0737 cassette but harbored genomic islands with LMOf2365_2059, thus yielding the DS 1/2b/3b profile. Lastly, we investigated a novel, dual 1/2a-1/2b profile obtained using the DS with 21 serotype 1/2a strains of four STs harboring both the lmo0737 cassette and genomic islands with LMOf2365_2059. The findings suggest that for certain strains and clones of L. monocytogenes the DS designations should be viewed with caution and complemented with alternative tools, e.g., traditional serotyping or WGS analysis. IMPORTANCE Listeria monocytogenes is a foodborne pathogen responsible for severe illness (listeriosis), especially in pregnant women and their fetuses, immunocompromised individuals, and the elderly. Three serotypes, 1/2a, 1/2b, and 4b, account for most human listeriosis, with certain serotype 4b clonal complexes (CCs) overrepresented in human disease. Serotyping remains extensively employed in Listeria epidemiologic investigations, and a multiplex PCR-based serotyping scheme is widely used. However, the PCR gene targets can be lost or gained via horizontal gene transfer, leading to novel PCR profiles without known serotype designations or to incorrect serotype assignments. Thus, an entire serotype 4b clone of the hypervirulent CC2 would be misidentified as serotype 1/2b, and several strains of serotype 1/2a would be identified as serotype 1/2b. Such challenges are especially common in novel clones from underexplored habitats, e.g., wildlife and surface water. The findings suggest caution in application of molecular serotyping, while highlighting Listeria's diversity and potential for horizontal gene transfer.

Horizontal gene transfer enables programmable gene stability in synthetic microbiota

The functions of many microbial communities exhibit remarkable stability despite fluctuations in the compositions of these communities. To date, a mechanistic understanding of this function-composition decoupling is lacking. Statistical mechanisms have been commonly hypothesized to explain such decoupling. Here, we proposed that dynamic mechanisms, mediated by horizontal gene transfer (HGT), also enable the independence of functions from the compositions of microbial communities. We combined theoretical analysis with numerical simulations to illustrate that HGT rates can determine the stability of gene abundance in microbial communities. We further validated these predictions using engineered microbial consortia of different complexities transferring one or more than a dozen clinically isolated plasmids, as well as through the reanalysis of data from the literature. Our results demonstrate a generalizable strategy to program the gene stability of microbial communities.

Atypical Serogroup IVb-v1 of Listeria monocytogenes Assigned to New ST2801, Widely Spread and Persistent in the Environment of a Pork-Meat Producing Plant of Central Italy

In this study, we characterized 84 Listeria monocytogenes (Lm) strains having an atypical IVb-v1 profile and isolated in a meat producing plant of Central Italy. They were assigned to the new MLST type ST2801 (CC218). The new ST was widespread in the food-producing environment where it was able to persist for over a year even after cleaning and sanitation. Cluster analysis identified three main clusters genetically close to each other (0–22 allelic differences and 0–28 SNPs) from two different cgMLST types, suggesting a common source. The coexistence of closely related clusters over time could be the result of a different evolution path starting from a common ancestor first introduced in the plant and/or the consequence of the repetitive reintroduction of closely related clones probably by raw materials. All the strains presented several determinants for heavy metals resistance, stress response, biofilm production, and multidrug efflux pumps with no significant differences among the clusters. A total of 53 strains carried pLI100 and the j1776 plasmids, while in one strain, the pLM33 was found in addition to pLI100. Only the strains carrying plasmids presented cadA and cadC for cadmium resistance and the mco gene encoding a multicopper oxidase and gerN for an additional Na+/H+-K+ antiporter. All the strains presented a virulence profile including a full-length inlA gene and the additional LIPI-3. The isolation of a new ST with a large pattern of stress-adaptation genes and able to persist is an important contribution to deepening the current knowledge on the uncommon IVb-v1 and in general on the genomic diversity of Lm.

Prevalence and Clonal Diversity of over 1,200 Listeria monocytogenes Isolates Collected from Public Access Waters near Produce Production Areas on the Central California Coast during 2011 to 2016

A 5-year survey of public access surface waters in an agricultural region of the Central California Coast was done to assess the prevalence of the foodborne pathogen Listeria monocytogenes. In nature, L. monocytogenes lives as a saprophyte in soil and water, which are reservoirs for contamination of preharvest produce. Moore swabs were deployed biweekly in lakes, ponds, streams, and rivers during 2011 to 2016. L. monocytogenes was recovered in 1,224 of 2,922 samples, resulting in 41.9% prevalence. Multiple subtypes were isolated from 97 samples, resulting in 1,323 L. monocytogenes isolates. Prevalence was higher in winter and spring and after rain events in some waterways. Over 84% of the isolates were serotype 4b. Whole-genome sequencing was done on 1,248 isolates, and in silico multilocus sequence typing revealed 74 different sequence types (STs) and 39 clonal complexes (CCs). The clones most isolated, CC639, CC183, and CC1, made up 27%, 19%, and 13%, respectively, of the sequenced isolates. Other types were CC663, CC6, CC842, CC4, CC2, CC5, and CC217. All sequenced isolates contained intact copies of core L. monocytogenes virulence genes, and pathogenicity islands LIPI-3 and LIPI-4 were identified in 73% and 63%, respectively, of the sequenced isolates. The virulence factor internalin A was predicted to be intact in all but four isolates, while genes important for sanitizer and heavy metal resistance were found in <5% of the isolates. These waters are not used for crop irrigation directly, but they are available to wildlife and can flood fields during heavy rains. IMPORTANCE Listeria monocytogenes serotype 4b and 1/2a strains are implicated in most listeriosis, and hypervirulent listeriosis stems from strains containing pathogenicity islands LIPI-3 and LIPI-4. The waters and sediments in the Central California Coast agricultural region contain widespread and diverse L. monocytogenes populations, and all the isolates contain intact virulence genes. Emerging clones CC183 and CC639 were the most abundant clones, and major clones CC1, CC4, and CC6 were well represented. CC183 was responsible for three produce-related outbreaks in the last 7 years. Most of the isolates in the survey differ from those of lesser virulence that are often isolated from foods and food processing plants because they contain genes encoding an intact virulence factor, internalin A, and most did not contain genes for sanitizer and heavy metal resistance. This isolate collection is important for understanding L. monocytogenes populations in agricultural and natural regions.

Genetic Diversity of Listeria monocytogenes Isolated From Three Commercial Tree Fruit Packinghouses and Evidence of Persistent and Transient Contamination

Whole genome analysis was performed on 501 isolates obtained from a previous survey which recovered 139 positive environmental sponge samples (i.e., up to 4 isolates per sample) from a total of 719 samples collected at 40 standardized sites in 3 commercial apple packinghouse facilities (i.e., P1, P2, and P3) over 3 successive seasons in a single production year. After excluding duplicated isolates, the data from 156 isolates revealed the clonal diversity of L. monocytogenes and allowed the detection of transient contamination, persistent contamination, and cross-area transmission events. Facility P2 with the poorest sanitary conditions had the least diversity (Shannon's index of 0.38). P2 contained a Clonal Complex (CC) 554, serogroup IVb-v1 strain that persisted throughout the year and spread across the entire facility, a singleton Sequence Type (ST) 1003, lineage III strain that persisted through two seasons and spread across two areas of the facility, and 3 other clones from transient contaminations. P1 and P3, facilities with better sanitary conditions, had much higher diversity (i.e., 15 clones with a Shannon's index of 2.49 and 10 clones with a Shannon's index of 2.10, respectively) that were the result of transient contamination. Facilities P1 and P3 had the highest incidence (43.1%) of lineage III isolates, followed by lineage I (31.3%) and lineage II (25.5%) isolates. Only 1 isolate in the three facilities contained a premature stop codon in virulence gene inlA. Fourteen samples yielded 2-3 clones per sample, demonstrating the importance of choosing appropriate methodologies and selecting a sufficient number of isolates per sample for studying L. monocytogenes diversity. Only 1 isolate, belonging to CC5 and from facility P3, contained a known plasmid, and this was also the only isolate containing benzalkonium chloride tolerance genes. The persistent CC554 strain did not exhibit stronger sanitizer resistance than other isolates and did not contain any confirmed molecular determinants of L. monocytogenes stress resistance that were differentially present in other isolates, such as genes involved in sanitizer tolerance, heavy metal resistance, biofilm-forming, stress survival islet 1 (SSI-1), stress survival islet 2 (SSI-2) or Listeria genomic island (LGI2).

Hyper-Virulent Listeria monocytogenes Strains Associated With Respiratory Infections in Central Italy

Listeria monocytogenes (Lm) is a foodborne pathogen causing listeriosis. Invasive forms of the disease mainly manifest as septicaemia, meningitis and maternal-neonatal infections. Lm-associated respiratory infections are very rare and little known. We reported two Lm respiratory infection cases occurred in Central Italy during the summer of 2020, in the midst of the SARS-CoV2 pandemic. In addition to collect the epidemiological and clinical characteristics of the patients, we used Whole Genome Sequencing to study the genomes of the Lm isolates investigating their virulence and antimicrobial profiles and the presence of genetic mobile elements. Both the strains belonged to hypervirulent MLST clonal complexes (CC). In addition to the Listeria Pathogenicity Island 1 (LIPI-1), the CC1 strain also carried LIPI-3 and the CC4 both LIPI-3 and LIPI-4. Genetic determinants for antimicrobial and disinfectants resistance were found. The CC1 genome presented prophage sequences but they did not interrupt the comK gene, involved in the phagosomal escape of Lm. None of the strains carried plasmids. Lm is an important, although rare, opportunistic pathogen for respiratory tract and lung infections. To avoid dangerous diagnostic delays of these severe clinical forms, it is important to sensitize hospital laboratories to this rare manifestation of listeriosis considering Lm in the differential diagnosis of respiratory infections.

Harnessing Whole Genome Sequence Data for Facility-Specific Signatures for Listeria monocytogenes: A Case Study With Turkey Processing Plants in the United States

Listeria monocytogenes is a Gram-positive foodborne pathogen responsible for the severe disease listeriosis and notorious for its ability to persist in food processing plants, leading to contamination of processed, ready-to-eat foods. L. monocytogenes persistence in various food processing environments (FPEs) has been extensively investigated by various subtyping tools, with increasing use of whole genome sequencing (WGS). However, major knowledge gaps remain. There is a need for facility-specific molecular signatures not only for adequate attribution of L. monocytogenes to a specific FPE but also for improved understanding of the ecology and evolution of L. monocytogenes in the food processing ecosystem. Furthermore, multiple strains can be recovered from a single FPE sample, but their diversity can be underestimated with common molecular subtyping tools. In this study we investigated a panel of 54 L. monocytogenes strains from four turkey processing plants in the United States. A combination of WGS and phenotypic assays was employed to assess strain persistence as well as identify facility-specific molecular signatures. Comparative analysis of allelic variation across the whole genome revealed that allelic profiles have the potential to be specific to individual processing plants. Certain allelic profiles remained associated with individual plants even when closely-related strains from other sources were included in the analysis. Furthermore, for certain sequence types (STs) based on the seven-locus multilocus sequence typing scheme, presence and location of premature stop codons in inlA, inlB length, prophage sequences, and the sequence content of a genomic hotspot could serve as plant-specific signatures. Interestingly, the analysis of different isolates from the same environmental sample revealed major differences not only in serotype and ST, but even in the sequence content of strains of the same ST. This study highlights the potential for WGS data to be deployed for identification of facility-specific signatures, thus facilitating the tracking of strain movement through the food chain. Furthermore, deployment of WGS for intra-sample strain analysis allows for a more complete environmental surveillance of L. monocytogenes in food processing facilities, reducing the risk of failing to detect strains that may be clinically relevant and potentially novel.

Characterization of Mobile Genetic Elements Using Long-Read Sequencing for Tracking Listeria monocytogenes from Food Processing Environments

Recently developed nanopore sequencing technologies offer a unique opportunity to rapidly close the genome and to identify complete sequences of mobile genetic elements (MGEs). In this study, 17 isolates of Listeria monocytogenes (Lm) epidemic clone II (ECII) from seven ready-to-eat meat or poultry processing facilities, not known to be associated with outbreaks, were shotgun sequenced, and among them, five isolates were further subjected to long-read sequencing. Additionally, 26 genomes of Lm ECII isolates associated with three listeriosis outbreaks in the U.S. and South Africa were obtained from the National Center for Biotechnology Information (NCBI) database and analyzed to evaluate if MGEs may be used as a high-resolution genetic marker for identifying and sourcing the origin of Lm. The analyses identified four comK prophages in 11 non-outbreak isolates from four facilities and three comK prophages in 20 isolates associated with two outbreaks that occurred in the U.S. In addition, three different plasmids were identified among 10 non-outbreak isolates and 14 outbreak isolates. Each comK prophage and plasmid was conserved among the isolates sharing it. Different prophages from different facilities or outbreaks had significant genetic variations, possibly due to horizontal gene transfer. Phylogenetic analysis showed that isolates from the same facility or the same outbreak always closely clustered. The time of most recent common ancestor of the Lm ECII isolates was estimated to be in March 1816 with the average nucleotide substitution rate of 3.1 × 10⁻⁷ substitutions per site per year. This study showed that complete MGE sequences provide a good signal to determine the genetic relatedness of Lm isolates, to identify persistence or repeated contamination that occurred within food processing environment, and to study the evolutionary history among closely related isolates.