High-Resolution Genomic Comparisons within Salmonella enterica Serotypes Derived from Beef Feedlot Cattle: Parsing the Roles of Cattle Source, Pen, Animal, Sample Type, and Production Period.

Salmonella enterica is a major foodborne pathogen, and contaminated beef products have been identified as one of the primary sources of Salmonella-related outbreaks. Pathogenicity and antibiotic resistance of Salmonella are highly serotype and subpopulation specific, which makes it essential to understand high-resolution Salmonella population dynamics in cattle. Time of year, source of cattle, pen, and sample type (i.e., feces, hide, or lymph nodes) have previously been identified as important factors influencing the serotype distribution of Salmonella (e.g., Anatum, Lubbock, Cerro, Montevideo, Kentucky, Newport, and Norwich) that were isolated from a longitudinal sampling design in a research feedlot. In this study, we performed high-resolution genomic comparisons of Salmonella isolates within each serotype using both single-nucleotide polymorphism-based maximum-likelihood phylogeny and hierarchical clustering of core-genome multilocus sequence typing. The importance of the aforementioned features in clonal Salmonella expansion was further explored using a supervised machine learning algorithm. In addition, we identified and compared the resistance genes, plasmids, and pathogenicity island profiles of the isolates within each subpopulation. Our findings indicate that clonal expansion of Salmonella strains in cattle was mainly influenced by the randomization of block and pen, as well as the origin/source of the cattle, i.e., regardless of sampling time and sample type (i.e., feces, lymph node, or hide). Further research is needed concerning the role of the feedlot pen environment prior to cattle placement to better understand carryover contributions of existing strains of Salmonella and their bacteriophages. IMPORTANCESalmonella serotypes isolated from outbreaks in humans can also be found in beef cattle and feedlots. Virulence factors and antibiotic resistance are among the primary defense mechanisms of Salmonella, and are often associated with clonal expansion. This makes understanding the subpopulation dynamics of Salmonella in cattle critical for effective mitigation. There remains a gap in the literature concerning subpopulation dynamics within Salmonella serotypes in feedlot cattle from the beginning of feeding up until slaughter. Here, we explore Salmonella population dynamics within each serotype using core-genome phylogeny and hierarchical classifications. We used machine learning to quantitatively parse the relative importance of both hierarchical and longitudinal clustering among cattle host samples. Our results reveal that Salmonella populations in cattle are highly clonal over a 6-month study period and that clonal dissemination of Salmonella in cattle is mainly influenced spatially by experimental block and pen, as well by the geographical origin of the cattle.

Salmonella enterica and Escherichia coli in Wheat Flour: Detection and Serotyping by a Quasimetagenomic Approach Assisted by Magnetic Capture, Multiple-Displacement Amplification, and Real-Time Sequencing.

Food safety is a new area for novel applications of metagenomics analysis, which not only can detect and subtype foodborne pathogens in a single workflow but may also produce additional information with in-depth analysis capabilities. In this study, we applied a quasimetagenomic approach by combining short-term enrichment, immunomagnetic separation (IMS), multiple-displacement amplification (MDA), and nanopore sequencing real-time analysis for simultaneous detection of Salmonella and Escherichia coli in wheat flour. Tryptic soy broth was selected for the 12-h enrichment of samples at 42°C. Enrichments were subjected to IMS using beads capable of capturing both Salmonella and E. coli MDA was performed on harvested beads, and amplified DNA fragments were subjected to DNA library preparation for sequencing. Sequencing was performed on a portable device with real-time basecalling adaptability, and resulting sequences were subjected to two parallel pipelines for further analysis. After 1 h of sequencing, the quasimetagenomic approach could detect all targets inoculated at approximately 1 CFU/g flour to the species level. Discriminatory power was determined by simultaneous detection of dual inoculums of Salmonella and E. coli, absence of detection in control samples, and consistency in microbial flora composition of the same flour samples over several rounds of experiments. The total turnaround time for detection was approximately 20 h. Longer sequencing for up to 15 h enabled serotyping for many of the samples with more than 99% genome coverage, which could be subjected to other appropriate genetic analysis pipelines in less than a total of 36 h.IMPORTANCE Enterohemorrhagic Escherichia coli (EHEC) and Salmonella are of serious concern in low-moisture foods, including wheat flour and its related products, causing illnesses, outbreaks, and recalls. The development of advanced detection methods based on molecular principles of analysis is essential to incorporate into interventions intended to reduce the risk from these pathogens. In this work, a quasimetagenomic method based on real-time sequencing analysis and assisted by magnetic capture and DNA amplification was developed. This protocol is capable of detecting multiple Salmonella and/or E. coli organisms in the sample within less than a day, and it can also generate sufficient whole-genome sequences of the target organisms suitable for subsequent bioinformatics analysis. Multiplex detection and identification were accomplished in less than 20 h and additional whole-genome analyses of different nature were attained within 36 h, in contrast to the several days required in previous sequencing pipelines.

Zoonotic Source Attribution of Salmonella enterica Serotype Typhimurium Using Genomic Surveillance Data, United States.

Increasingly, routine surveillance and monitoring of foodborne pathogens using whole-genome sequencing is creating opportunities to study foodborne illness epidemiology beyond routine outbreak investigations and case-control studies. Using a global phylogeny of Salmonella enterica serotype Typhimurium, we found that major livestock sources of the pathogen in the United States can be predicted through whole-genome sequencing data. Relatively steady rates of sequence divergence in livestock lineages enabled the inference of their recent origins. Elevated accumulation of lineage-specific pseudogenes after divergence from generalist populations and possible metabolic acclimation in a representative swine isolate indicates possible emergence of host adaptation. We developed and retrospectively applied a machine learning Random Forest classifier for genomic source prediction of Salmonella Typhimurium that correctly attributed 7 of 8 major zoonotic outbreaks in the United States during 1998-2013. We further identified 50 key genetic features that were sufficient for robust livestock source prediction.

Comparative Genomics Reveals the Diversity of Restriction-Modification Systems and DNA Methylation Sites in Listeria monocytogenes.

Listeria monocytogenes is a bacterial pathogen that is found in a wide variety of anthropogenic and natural environments. Genome sequencing technologies are rapidly becoming a powerful tool in facilitating our understanding of how genotype, classification phenotypes, and virulence phenotypes interact to predict the health risks of individual bacterial isolates. Currently, 57 closed L. monocytogenes genomes are publicly available, representing three of the four phylogenetic lineages, and they suggest that L. monocytogenes has high genomic synteny. This study contributes an additional 15 closed L. monocytogenes genomes that were used to determine the associations between the genome and methylome with host invasion magnitude. In contrast to previous findings, large chromosomal inversions and rearrangements were detected in five isolates at the chromosome terminus and within rRNA genes, including a previously undescribed inversion within rRNA-encoding regions. Each isolate's epigenome contained highly diverse methyltransferase recognition sites, even within the same serotype and methylation pattern. Eleven strains contained a single chromosomally encoded methyltransferase, one strain contained two methylation systems (one system on a plasmid), and three strains exhibited no methylation, despite the occurrence of methyltransferase genes. In three isolates a new, unknown DNA modification was observed in addition to diverse methylation patterns, accompanied by a novel methylation system. Neither chromosome rearrangement nor strain-specific patterns of epigenome modification observed within virulence genes were correlated with serotype designation, clonal complex, or in vitro infectivity. These data suggest that genome diversity is larger than previously considered in L. monocytogenes and that as more genomes are sequenced, additional structure and methylation novelty will be observed in this organism. IMPORTANCE: Listeria monocytogenes is the causative agent of listeriosis, a disease which manifests as gastroenteritis, meningoencephalitis, and abortion. Among Salmonella, Escherichia coli, Campylobacter, and Listeria-causing the most prevalent foodborne illnesses-infection by L. monocytogenes carries the highest mortality rate. The ability of L. monocytogenes to regulate its response to various harsh environments enables its persistence and transmission. Small-scale comparisons of L. monocytogenes focusing solely on genome contents reveal a highly syntenic genome yet fail to address the observed diversity in phenotypic regulation. This study provides a large-scale comparison of 302 L. monocytogenes isolates, revealing the importance of the epigenome and restriction-modification systems as major determinants of L. monocytogenes phylogenetic grouping and subsequent phenotypic expression. Further examination of virulence genes of select outbreak strains reveals an unprecedented diversity in methylation statuses despite high degrees of genome conservation.

Whole-Genome Sequencing of Drug-Resistant Salmonella enterica Isolates from Dairy Cattle and Humans in New York and Washington States Reveals Source and Geographic Associations.

Multidrug-resistant (MDR) Salmonella enterica can be spread from cattle to humans through direct contact with animals shedding Salmonella as well as through the food chain, making MDR Salmonella a serious threat to human health. The objective of this study was to use whole-genome sequencing to compare antimicrobial-resistant (AMR) Salmonella enterica serovars Typhimurium, Newport, and Dublin isolated from dairy cattle and humans in Washington State and New York State at the genotypic and phenotypic levels. A total of 90 isolates were selected for the study (37 S Typhimurium, 32 S Newport, and 21 S Dublin isolates). All isolates were tested for phenotypic antibiotic resistance to 12 drugs using Kirby-Bauer disk diffusion. AMR genes were detected in the assembled genome of each isolate using nucleotide BLAST and ARG-ANNOT. Genotypic prediction of phenotypic resistance resulted in a mean sensitivity of 97.2 and specificity of 85.2. Sulfamethoxazole-trimethoprim resistance was observed only in human isolates (P < 0.05), while resistance to quinolones and fluoroquinolones was observed only in 6 S Typhimurium isolates from humans in Washington State. S Newport isolates showed a high degree of AMR profile similarity, regardless of source. S Dublin isolates from New York State differed from those from Washington State based on the presence/absence of plasmid replicons, as well as phenotypic AMR susceptibility/nonsusceptibility (P < 0.05). The results of this study suggest that distinct factors may contribute to the emergence and dispersal of AMR S. enterica in humans and farm animals in different regions.IMPORTANCE The use of antibiotics in food-producing animals has been hypothesized to select for AMR Salmonella enterica and associated AMR determinants, which can be transferred to humans through different routes. Previous studies have sought to assess the degree to which AMR livestock- and human-associated Salmonella strains overlap, as well as the spatial distribution of Salmonella's associated AMR determinants, but have often been limited by the degree of resolution at which isolates can be compared. Here, a comparative genomics study of livestock- and human-associated Salmonella strains from different regions of the United States shows that while many AMR genes and phenotypes were confined to human isolates, overlaps between the resistomes of bovine and human-associated Salmonella isolates were observed on numerous occasions, particularly for S Newport. We have also shown that whole-genome sequencing can be used to reliably predict phenotypic resistance across Salmonella isolated from bovine sources.

Targeted Amplicon Sequencing for Single-Nucleotide-Polymorphism Genotyping of Attaching and Effacing Escherichia coli O26:H11 Cattle Strains via a High-Throughput Library Preparation Technique.

Enterohemorrhagic Escherichia coli (EHEC) O26:H11, a serotype within Shiga toxin-producing E. coli (STEC) that causes severe human disease, has been considered to have evolved from attaching and effacing E. coli (AEEC) O26:H11 through the acquisition of a Shiga toxin-encoding gene. Targeted amplicon sequencing using next-generation sequencing technology of 48 phylogenetically informative single-nucleotide polymorphisms (SNPs) and three SNPs differentiating Shiga toxin-positive (stx-positive) strains from Shiga toxin-negative (stx-negative) strains were used to infer the phylogenetic relationships of 178 E. coli O26:H11 strains (6 stx-positive strains and 172 stx-negative AEEC strains) from cattle feces to 7 publically available genomes of human clinical strains. The AEEC cattle strains displayed synonymous SNP genotypes with stx2-positive sequence type 29 (ST29) human O26:H11 strains, while stx1 ST21 human and cattle strains clustered separately, demonstrating the close phylogenetic relatedness of these Shiga toxin-negative AEEC cattle strains and human clinical strains. With the exception of seven stx-negative strains, five of which contained espK, three stx-related SNPs differentiated the STEC strains from non-STEC strains, supporting the hypothesis that these AEEC cattle strains could serve as a potential reservoir for new or existing pathogenic human strains. Our results support the idea that targeted amplicon sequencing for SNP genotyping expedites strain identification and genetic characterization of E. coli O26:H11, which is important for food safety and public health.

Determination of Evolutionary Relationships of Outbreak-Associated Listeria monocytogenes Strains of Serotypes 1/2a and 1/2b by Whole-Genome Sequencing.

We used whole-genome sequencing to determine evolutionary relationships among 20 outbreak-associated clinical isolates of Listeria monocytogenes serotypes 1/2a and 1/2b. Isolates from 6 of 11 outbreaks fell outside the clonal groups or "epidemic clones" that have been previously associated with outbreaks, suggesting that epidemic potential may be widespread in L. monocytogenes and is not limited to the recognized epidemic clones. Pairwise comparisons between epidemiologically related isolates within clonal complexes showed that genome-level variation differed by 2 orders of magnitude between different comparisons, and the distribution of point mutations (core versus accessory genome) also varied. In addition, genetic divergence between one closely related pair of isolates from a single outbreak was driven primarily by changes in phage regions. The evolutionary analysis showed that the changes could be attributed to horizontal gene transfer; members of the diverse bacterial community found in the production facility could have served as the source of novel genetic material at some point in the production chain. The results raise the question of how to best utilize information contained within the accessory genome in outbreak investigations. The full magnitude and complexity of genetic changes revealed by genome sequencing could not be discerned from traditional subtyping methods, and the results demonstrate the challenges of interpreting genetic variation among isolates recovered from a single outbreak. Epidemiological information remains critical for proper interpretation of nucleotide and structural diversity among isolates recovered during outbreaks and will remain so until we understand more about how various population histories influence genetic variation.

Whole-Genome Sequencing Allows for Improved Identification of Persistent Listeria monocytogenes in Food-Associated Environments.

While the food-borne pathogen Listeria monocytogenes can persist in food associated environments, there are no whole-genome sequence (WGS) based methods to differentiate persistent from sporadic strains. Whole-genome sequencing of 188 isolates from a longitudinal study of L. monocytogenes in retail delis was used to (i) apply single-nucleotide polymorphism (SNP)-based phylogenetics for subtyping of L. monocytogenes, (ii) use SNP counts to differentiate persistent from repeatedly reintroduced strains, and (iii) identify genetic determinants of L. monocytogenes persistence. WGS analysis revealed three prophage regions that explained differences between three pairs of phylogenetically similar populations with pulsed-field gel electrophoresis types that differed by ≤3 bands. WGS-SNP-based phylogenetics found that putatively persistent L. monocytogenes represent SNP patterns (i) unique to a single retail deli, supporting persistence within the deli (11 clades), (ii) unique to a single state, supporting clonal spread within a state (7 clades), or (iii) spanning multiple states (5 clades). Isolates that formed one of 11 deli-specific clades differed by a median of 10 SNPs or fewer. Isolates from 12 putative persistence events had significantly fewer SNPs (median, 2 to 22 SNPs) than between isolates of the same subtype from other delis (median up to 77 SNPs), supporting persistence of the strain. In 13 events, nearly indistinguishable isolates (0 to 1 SNP) were found across multiple delis. No individual genes were enriched among persistent isolates compared to sporadic isolates. Our data show that WGS analysis improves food-borne pathogen subtyping and identification of persistent bacterial pathogens in food associated environments.

Selection and Characterization of Phage-Resistant Mutant Strains of Listeria monocytogenes Reveal Host Genes Linked to Phage Adsorption.

Listeria-infecting phages are readily isolated from Listeria-containing environments, yet little is known about the selective forces they exert on their host. Here, we identified that two virulent phages, LP-048 and LP-125, adsorb to the surface of Listeria monocytogenes strain 10403S through different mechanisms. We isolated and sequenced, using whole-genome sequencing, 69 spontaneous mutant strains of 10403S that were resistant to either one or both phages. Mutations from 56 phage-resistant mutant strains with only a single mutation mapped to 10 genes representing five loci on the 10403S chromosome. An additional 12 mutant strains showed two mutations, and one mutant strain showed three mutations. Two of the loci, containing seven of the genes, accumulated the majority (n = 64) of the mutations. A representative mutant strain for each of the 10 genes was shown to resist phage infection through mechanisms of adsorption inhibition. Complementation of mutant strains with the associated wild-type allele was able to rescue phage susceptibility for 6 out of the 10 representative mutant strains. Wheat germ agglutinin, which specifically binds to N-acetylglucosamine, bound to 10403S and mutant strains resistant to LP-048 but did not bind to mutant strains resistant to only LP-125. We conclude that mutant strains resistant to only LP-125 lack terminal N-acetylglucosamine in their wall teichoic acid (WTA), whereas mutant strains resistant to both phages have disruptive mutations in their rhamnose biosynthesis operon but still possess N-acetylglucosamine in their WTA.

Transcriptomic Analysis of the Adaptation of Listeria monocytogenes to Growth on Vacuum-Packed Cold Smoked Salmon.

The foodborne pathogen Listeria monocytogenes is able to survive and grow in ready-to-eat foods, in which it is likely to experience a number of environmental stresses due to refrigerated storage and the physicochemical properties of the food. Little is known about the specific molecular mechanisms underlying survival and growth of L. monocytogenes under different complex conditions on/in specific food matrices. Transcriptome sequencing (RNA-seq) was used to understand the transcriptional landscape of L. monocytogenes strain H7858 grown on cold smoked salmon (CSS; water phase salt, 4.65%; pH 6.1) relative to that in modified brain heart infusion broth (MBHIB; water phase salt, 4.65%; pH 6.1) at 7°C. Significant differential transcription of 149 genes was observed (false-discovery rate [FDR], <0.05; fold change, ≥2.5), and 88 and 61 genes were up- and downregulated, respectively, in H7858 grown on CSS relative to the genes in H7858 grown in MBHIB. In spite of these differences in transcriptomes under these two conditions, growth parameters for L. monocytogenes were not significantly different between CSS and MBHIB, indicating that the transcriptomic differences reflect how L. monocytogenes is able to facilitate growth under these different conditions. Differential expression analysis and Gene Ontology enrichment analysis indicated that genes encoding proteins involved in cobalamin biosynthesis as well as ethanolamine and 1,2-propanediol utilization have significantly higher transcript levels in H7858 grown on CSS than in that grown in MBHIB. Our data identify specific transcriptional profiles of L. monocytogenes growing on vacuum-packaged CSS, which may provide targets for the development of novel and improved strategies to control L. monocytogenes growth on this ready-to-eat food.

Molecular evolution patterns reveal life history features of mycoplasma-related endobacteria associated with arbuscular mycorrhizal fungi.

The mycoplasma-related endobacteria (MRE), representing a recently discovered lineage of Mollicutes, are widely distributed across arbuscular mycorrhizal fungi (AMF, Glomeromycota). AMF colonize roots of most terrestrial plants and improve plant mineral nutrient uptake in return for plant-assimilated carbon. The role of MRE in the biology of their fungal hosts is unknown. To start characterizing this association, we assessed partitioning of MRE genetic diversity within AMF individuals and across the AMF phylogeographic range. We further used molecular evolution patterns to make inferences about MRE codivergence with AMF, their lifestyle and antiquity of the Glomeromycota-MRE association. While we did not detect differentiation between MRE derived from different continents, high levels of diversity were apparent in MRE populations within AMF host individuals. MRE exhibited significant codiversification with AMF over ecological time and the absence of codivergence over evolutionary time. Moreover, genetic recombination was evident in MRE. These patterns indicate that, while MRE transmission is predominantly vertical, their complex intrahost populations are likely generated by horizontal transmission and recombination. Based on predictions of evolutionary theory, we interpreted these observations as a suggestion that MRE may be antagonists of AMF. Finally, we detected a marginally significant signature of codivergence of MRE with Glomeromycota and the Endogone lineage of Mucoromycotina, implying that the symbiosis between MRE and fungi may predate the divergence between these two groups of fungi.

Listeria booriae sp. nov. and Listeria newyorkensis sp. nov., from food processing environments in the USA.

Sampling of seafood and dairy processing facilities in the north-eastern USA produced 18 isolates of Listeria spp. that could not be identified at the species-level using traditional phenotypic and genotypic identification methods. Results of phenotypic and genotypic analyses suggested that the isolates represent two novel species with an average nucleotide blast identity of less than 92% with previously described species of the genus Listeria. Phylogenetic analyses based on whole genome sequences, 16S rRNA gene and sigB gene sequences confirmed that the isolates represented by type strain FSL M6-0635(T) and FSL A5-0209 cluster phylogenetically with Listeria cornellensis. Phylogenetic analyses also showed that the isolates represented by type strain FSL A5-0281(T) cluster phylogenetically with Listeria riparia. The name Listeria booriae sp. nov. is proposed for the species represented by type strain FSL A5-0281(T) ( =DSM 28860(T) =LMG 28311(T)), and the name Listeria newyorkensis sp. nov. is proposed for the species represented by type strain FSL M6-0635(T) ( =DSM 28861(T) =LMG 28310(T)). Phenotypic and genotypic analyses suggest that neither species is pathogenic.

Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States.

BACKGROUND: Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar. RESULTS: In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here "NY bovine clade"), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche. CONCLUSIONS: Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.

Genomic comparison of sporeforming bacilli isolated from milk.

BACKGROUND: Sporeformers in the order Bacillales are important contributors to spoilage of pasteurized milk. While only a few Bacillus and Viridibacillus strains can grow in milk at 6°C, the majority of Paenibacillus isolated from pasteurized fluid milk can grow under these conditions. To gain a better understanding of genomic features of these important spoilage organisms and to identify candidate genomic features that may facilitate cold growth in milk, we performed a comparative genomic analysis of selected dairy associated sporeformers representing isolates that can and cannot grow in milk at 6°C. RESULTS: The genomes for seven Paenibacillus spp., two Bacillus spp., and one Viridibacillus sp. isolates were sequenced. Across the genomes sequenced, we identified numerous genes encoding antimicrobial resistance mechanisms, bacteriocins, and pathways for synthesis of non-ribosomal peptide antibiotics. Phylogenetic analysis placed genomes representing Bacillus, Paenibacillus and Viridibacillus into three distinct well supported clades and further classified the Paenibacillus strains characterized here into three distinct clades, including (i) clade I, which contains one strain able to grow at 6°C in skim milk broth and one strain not able to grow under these conditions, (ii) clade II, which contains three strains able to grow at 6°C in skim milk broth, and (iii) clade III, which contains two strains unable to grow under these conditions. While all Paenibacillus genomes were found to include multiple copies of genes encoding ß-galactosidases, clade II strains showed significantly higher numbers of genes encoding these enzymes as compared to clade III strains. Genome comparison of strains able to grow at 6°C and strains unable to grow at this temperature identified numerous genes encoding features that might facilitate the growth of Paenibacillus in milk at 6°C, including peptidases with cold-adapted features (flexibility and disorder regions in the protein structure) and cold-adaptation related proteins (DEAD-box helicases, chaperone DnaJ). CONCLUSIONS: Through a comparative genomics approach we identified a number of genomic features that may relate to the ability of selected Paenibacillus strains to cause spoilage of refrigerated fluid milk. With additional experimental evidence, these data will facilitate identification of targets to detect and control Gram positive spore formers in fluid milk.

Rapid whole-genome sequencing for surveillance of Salmonella enterica serovar enteritidis.

For Salmonella enterica serovar Enteritidis, 85% of isolates can be classified into 5 pulsed-field gel electrophoresis (PFGE) types. However, PFGE has limited discriminatory power for outbreak detection. Although whole-genome sequencing has been found to improve discrimination of outbreak clusters, whether this procedure can be used in real-time in a public health laboratory is not known. Therefore, we conducted a retrospective and prospective analysis. The retrospective study investigated isolates from 1 confirmed outbreak. Additional cases could be attributed to the outbreak strain on the basis of whole-genome data. The prospective study included 58 isolates obtained in 2012, including isolates from 1 epidemiologically defined outbreak. Whole-genome sequencing identified additional isolates that could be attributed to the outbreak, but which differed from the outbreak-associated PFGE type. Additional putative outbreak clusters were detected in the retrospective and prospective analyses. This study demonstrates the practicality of implementing this approach for outbreak surveillance in a state public health laboratory.

Genomic epidemiology of Salmonella enterica serotype Enteritidis based on population structure of prevalent lineages.

Salmonella enterica serotype Enteritidis is one of the most commonly reported causes of human salmonellosis. Its low genetic diversity, measured by fingerprinting methods, has made subtyping a challenge. We used whole-genome sequencing to characterize 125 S. enterica Enteritidis and 3 S. enterica serotype Nitra strains. Single-nucleotide polymorphisms were filtered to identify 4,887 reliable loci that distinguished all isolates from each other. Our whole-genome single-nucleotide polymorphism typing approach was robust for S. enterica Enteritidis subtyping with combined data for different strains from 2 different sequencing platforms. Five major genetic lineages were recognized, which revealed possible patterns of geographic and epidemiologic distribution. Analyses on the population dynamics and evolutionary history estimated that major lineages emerged during the 17th-18th centuries and diversified during the 1920s and 1950s.

Comparative genomic and morphological analyses of Listeria phages isolated from farm environments.

The genus Listeria is ubiquitous in the environment and includes the globally important food-borne pathogen Listeria monocytogenes. While the genomic diversity of Listeria has been well studied, considerably less is known about the genomic and morphological diversity of Listeria bacteriophages. In this study, we sequenced and analyzed the genomes of 14 Listeria phages isolated mostly from New York dairy farm environments as well as one related Enterococcus faecalis phage to obtain information on genome characteristics and diversity. We also examined 12 of the phages by electron microscopy to characterize their morphology. These Listeria phages, based on gene orthology and morphology, together with previously sequenced Listeria phages could be classified into five orthoclusters, including one novel orthocluster. One orthocluster (orthocluster I) consists of large genome (~135-kb) myoviruses belonging to the genus "Twort-like viruses," three orthoclusters (orthoclusters II to IV) contain small-genome (36- to 43-kb) siphoviruses with icosahedral heads, and the novel orthocluster V contains medium-sized-genome (~66-kb) siphoviruses with elongated heads. A novel orthocluster (orthocluster VI) of E. faecalis phages, with medium-sized genomes (~56 kb), was identified, which grouped together and shares morphological features with the novel Listeria phage orthocluster V. This new group of phages (i.e., orthoclusters V and VI) is composed of putative lytic phages that may prove to be useful in phage-based applications for biocontrol, detection, and therapeutic purposes.

Listeria floridensis sp. nov., Listeria aquatica sp. nov., Listeria cornellensis sp. nov., Listeria riparia sp. nov. and Listeria grandensis sp. nov., from agricultural and natural environments.

Sampling of agricultural and natural environments in two US states (Colorado and Florida) yielded 18 Listeria-like isolates that could not be assigned to previously described species using traditional methods. Using whole-genome sequencing and traditional phenotypic methods, we identified five novel species, each with a genome-wide average BLAST nucleotide identity (ANIb) of less than 85% to currently described species. Phylogenetic analysis based on 16S rRNA gene sequences and amino acid sequences of 31 conserved loci showed the existence of four well-supported clades within the genus Listeria; (i) a clade representing Listeria monocytogenes, L. marthii, L. innocua, L. welshimeri, L. seeligeri and L. ivanovii, which we refer to as Listeria sensu stricto, (ii) a clade consisting of Listeria fleischmannii and two newly described species, Listeria aquatica sp. nov. (type strain FSL S10-1188(T) = DSM 26686(T) = LMG 28120(T) = BEI NR-42633(T)) and Listeria floridensis sp. nov. (type strain FSL S10-1187(T) = DSM 26687(T) = LMG 28121(T) = BEI NR-42632(T)), (iii) a clade consisting of Listeria rocourtiae, L. weihenstephanensis and three novel species, Listeria cornellensis sp. nov. (type strain TTU A1-0210(T) = FSL F6-0969(T) = DSM 26689(T) = LMG 28123(T) = BEI NR-42630(T)), Listeria grandensis sp. nov. (type strain TTU A1-0212(T) = FSL F6-0971(T) = DSM 26688(T) = LMG 28122(T) = BEI NR-42631(T)) and Listeria riparia sp. nov. (type strain FSL S10-1204(T) = DSM 26685(T) = LMG 28119(T) = BEI NR- 42634(T)) and (iv) a clade containing Listeria grayi. Genomic and phenotypic data suggest that the novel species are non-pathogenic.

Analysis of Microbial Composition of Edible Insect Products Available for Human Consumption within the United States Using Traditional Microbiological Methods and Whole Genome Sequencing.

Edible insects offer a promising protein source for humans, but their food safety risks have not been previously investigated within the United States. Therefore, the aim of this study was to investigate the microbial content of processed edible insect products. A total of eight different types of edible insect products, including diving beetles, silkworms, grasshoppers, Jamaican crickets, mealworms, mole crickets, whole roasted crickets, and 100% pure cricket powder, were purchased from a large online retailer for the analysis. All the products were purchased in August 2022 and examined between August 2022 and November 2022. Traditional microbiological methods were employed to determine microbial counts for each product type using three replicates (total number of samples = 24). This included assessing aerobic bacterial spore, lactic acid bacteria, Enterobacteriaceae, total viable counts, and the presence of Salmonella. Additionally, whole genome sequencing was employed to further characterize selected colonies (n = 96). Microbial counts data were statistically analyzed using one-way ANOVA, while sequence data were taxonomically classified using Sepia.Bacilluscereusgroup isolates underwent additional characterization with Btyper3. Product type significantly influenced total viable counts, bacterial spore counts, and lactic acid bacteria counts (P = 0.00391, P = 0.0065, and P < 0.001, respectively), with counts ranging from < 1.70 to 6.01 Log10 CFU/g, <1.70 to 5.25 Log10 CFU/g, and < 1.70 to 4.86 Log10 CFU/g, respectively. Enterobacteriaceae were only detected in mole crickets (<2.30 Log10 CFU/g) and house cricket powder (<2.15 Log10 CFU/g). All samples were negative for Salmonella. Whole genome sequencing revealed the presence of 12 different bacterial genera among the analyzed isolates, with a majority belonging to the Bacillus genus. Some of the isolates of Bacillus cereus group were identified as biovar Emeticus. Overall, although edible insects offer a promising food alternative, the presence of Bacillus cereus group in some products could raise concerns regarding food safety.

Genomic characterization provides new insight into Salmonella phage diversity.

BACKGROUND: Salmonella is a widely distributed foodborne pathogen that causes tens of millions of salmonellosis cases globally every year. While the genomic diversity of Salmonella is increasingly well studied, our knowledge of Salmonella phage genomic diversity is still rather limited, despite the contributions of both lysogenic and lytic phages to Salmonella virulence, diversity and ecology (e.g., through horizontal gene transfer and Salmonella lysis). To gain a better understanding of phage diversity in a specific ecological niche, we sequenced 22 Salmonella phages isolated from a number of dairy farms from New York State (United States) and analyzed them using a comparative genomics approach. RESULTS: Classification of the 22 phages according to the presence/absence of orthologous genes allowed for classification into 8 well supported clusters. In addition to two phage clusters that represent novel virulent Salmonella phages, we also identified four phage clusters that each contained previously characterized phages from multiple continents. Our analyses also identified two clusters of phages that carry putative virulence (e.g., adhesins) and antimicrobial resistance (tellurite and bicyclomycin) genes as well as virulent and temperate transducing phages. Insights into phage evolution from our analyses include (i) identification of DNA metabolism genes that may facilitate nucleotide synthesis in phages with a G+C % distinct from Salmonella, and (ii) evidence of Salmonella phage tailspike and fiber diversity due to both single nucleotide polymorphisms and major re-arrangements, which may affect the host specificity of Salmonella phages. CONCLUSIONS: Genomics-based characterization of 22 Salmonella phages isolated from dairy farms allowed for identification of a number of novel Salmonella phages. While the comparative genomics analyses of these phages provide a number of new insights in the evolution and diversity of Salmonella phages, they only represent a first glimpse into the diversity of Salmonella phages that is likely to be discovered when phages from different environments are characterized.