Listeria monocytogenes is frequently found in foods and processing facilities, where it can persist, creating concerns for the food industry. Its ability to survive under a wide range of environmental conditions enhances the potential for cross-contamination of the final food products, leading to possible outbreaks of listeriosis. In this study, whole-genome sequencing (WGS) was applied as a tool to characterize and track 100 L. monocytogenes isolates collected from three food processing environments. These WGS data from environmental and food isolates were analyzed to (i) assess the genomic diversity of L. monocytogenes, (ii) identify possible source(s) of contamination, cross-contamination routes, and persistence, (iii) detect absence/presence of antimicrobial resistance-encoding genes, (iv) assess virulence genotypes, and (v) explore in vivo pathogenicity of selected L. monocytogenes isolates carrying different virulence genotypes. The predominant L. monocytogenes sublineages (SLs) identified were SL101 (21%), SL9 (17%), SL121 (12%), and SL5 (12%). Benzalkonium chloride (BC) tolerance-encoding genes were found in 62% of these isolates, a value that increased to 73% among putative persistent subgroups. The most prevalent gene was emrC followed by bcrABC, qacH-Tn6188, and qacC. The L. monocytogenes major virulence factor inlA was truncated in 31% of the isolates, and only one environmental isolate (L. monocytogenes CFS086) harbored all major virulence factors, including Listeria pathogenicity island 4 (LIPI-4), which has been shown to confer hypervirulence. A zebrafish embryo infection model showed a low (3%) embryo survival rate for all putatively hypervirulent L. monocytogenes isolates assayed. Higher embryo survival rates were observed following infection with unknown virulence potential (20%) and putatively hypovirulent (53 to 83%) L. monocytogenes isolates showing predicted pathogenic phenotypes inferred from virulence genotypes. IMPORTANCE This study extends current understanding of the genetic diversity among L. monocytogenes from various food products and food processing environments. Application of WGS-based strategies facilitated tracking of this pathogen of importance to human health along the production chain while providing insights into the pathogenic potential for some of the L. monocytogenes isolates recovered. These analyses enabled the grouping of selected isolates into three putative virulence categories according to their genotypes along with informing selection for phenotypic assessment of their pathogenicity using the zebrafish embryo infection model. It has also facilitated the identification of those isolates with genes conferring tolerance to commercially used biocides. Findings from this study highlight the potential for the application of WGS as a proactive tool to support food safety controls as applied to L. monocytogenes.

BackgroundAlthough serotype O157:H7 is the predominant enterohemorrhagic Escherichia coli (EHEC), outbreaks of non-O157 EHEC that cause severe foodborne illness, including hemolytic uremic syndrome have increased worldwide. In fact, non-O157 serotypes are now estimated to cause over half of all the Shiga toxin-producing Escherichia coli (STEC) cases, and outbreaks of non-O157 EHEC infections are frequently associated with serotypes O26, O45, O103, O111, O121, and O145. Currently, there are no complete genomes for O145 in public databases.ResultsWe determined the complete genome sequences of two O145 strains (EcO145), one linked to a US lettuce-associated outbreak (RM13514) and one to a Belgium ice-cream-associated outbreak (RM13516). Both strains contain one chromosome and two large plasmids, with genome sizes of 5,737,294 bp for RM13514 and 5,559,008 bp for RM13516. Comparative analysis of the two EcO145 genomes revealed a large core (5,173 genes) and a considerable amount of strain-specific genes. Additionally, the two EcO145 genomes display distinct chromosomal architecture, virulence gene profile, phylogenetic origin of Stx2a prophage, and methylation profile (methylome). Comparative analysis of EcO145 genomes to other completely sequenced STEC and other E. coli and Shigella genomes revealed that, unlike any other known non-O157 EHEC strain, EcO145 ascended from a common lineage with EcO157/EcO55. This evolutionary relationship was further supported by the pangenome analysis of the 10 EHEC str ains. Of the 4,192 EHEC core genes, EcO145 shares more genes with EcO157 than with the any other non-O157 EHEC strains.ConclusionsOur data provide evidence that EcO145 and EcO157 evolved from a common lineage, but ultimately each serotype evolves via a lineage-independent nature to EHEC by acquisition of the core set of EHEC virulence factors, including the genes encoding Shiga toxin and the large virulence plasmid. The large variation between the two EcO145 genomes suggests a distinctive evolutionary path between the two outbreak strains. The distinct methylome between the two EcO145 strains is likely due to the presence of a BsuBI/PstI methyltransferase gene cassette in the Stx2a prophage of the strain RM13514, suggesting a role of horizontal gene transfer-mediated epigenetic alteration in the evolution of individual EHEC strains.

The Campylobacter lari group is a phylogenetic clade within the epsilon subdivision of the Proteobacteria and is part of the thermotolerant Campylobacter spp., a division within the genus that includes the human pathogen Campylobacter jejuni. The C. lari group is currently composed of five species (C. lari, Campylobacter insulaenigrae, Campylobacter volucris, Campylobacter subantarcticus, and Campylobacter peloridis), as well as a group of strains termed the urease-positive thermophilic Campylobacter (UPTC) and other C. lari-like strains. Here we present the complete genome sequences of 11 C. lari group strains, including the five C. lari group species, four UPTC strains, and a lari-like strain isolated in this study. The genome of C. lari subsp. lari strain RM2100 was described previously. Analysis of the C. lari group genomes indicates that this group is highly related at the genome level. Furthermore, these genomes are strongly syntenic with minor rearrangements occurring only in 4 of the 12 genomes studied. The C. lari group can be bifurcated, based on the flagella and flagellar modification genes. Genomic analysis of the UPTC strains indicated that these organisms are variable but highly similar, closely related to but distinct from C. lari. Additionally, the C. lari group contains multiple genes encoding hemagglutination domain proteins, which are either contingency genes or linked to conserved contingency genes. Many of the features identified in strain RM2100, such as major deficiencies in amino acid biosynthesis and energy metabolism, are conserved across all 12 genomes, suggesting that these common features may play a role in the association of the C. lari group with coastal environments and watersheds.