Identifying a list of Salmonella serotypes of concern to target for reducing risk of salmonellosis.

There is an increasing awareness in the field of Salmonella epidemiology that focusing control efforts on those serotypes which cause severe human health outcomes, as opposed to broadly targeting all Salmonella, will likely lead to the greatest advances in decreasing the incidence of salmonellosis. Yet, little guidance exists to support validated, scientific selection of target serotypes. The goal of this perspective is to develop an approach to identifying serotypes of greater concern and present a case study using meat- and poultry-attributed outbreaks to examine challenges in developing a standardized framework for defining target serotypes.

Growth assessment of Salmonella enterica multi-serovar populations in poultry rinsates with commonly used enrichment and plating media.

Isolation of Salmonella from enrichment cultures of food or environmental samples is a complicated process. Numerous factors including fitness in various selective enrichment media, relative starting concentrations in pre-enrichment, and competition among multi-serovar populations and associated natural microflora, come together to determine which serovars are identified from a given sample. A recently developed approach for assessing the relative abundance (RA) of multi-serovar Salmonella populations (CRISPR-SeroSeq or Deep Serotyping, DST) is providing new insight into how these factors impact the serovars observed, especially when different selective enrichment methods are used to identify Salmonella from a primary enrichment sample. To illustrate this, we examined Salmonella-positive poultry pre-enrichment samples through the selective enrichment process in Tetrathionate (TT) and Rappaport Vassiliadis (RVS) broths and assessed recovery of serovars with each medium. We observed the RA of serovars detected post selective enrichment varied depending on the medium used, initial concentration, and competitive fitness factors, all which could result in minority serovars in pre-enrichment becoming dominant serovars post selective enrichment. The data presented provide a greater understanding of culture biases and lays the groundwork for investigations into robust enrichment and plating media combinations for detecting Salmonella serovars of greater concern for human health.

Escherichia coli O157:H7 tir 255 T > A allele strains differ in chromosomal and plasmid composition.

Shiga toxin-producing Escherichia coli (STEC) O157:H7 strains with the T allele in the translocated intimin receptor polymorphism (tir) 255 A > T gene associate with human disease more than strains with an A allele; however, the allele is not thought to be the direct cause of this difference. We sequenced a diverse set of STEC O157:H7 strains (26% A allele, 74% T allele) to identify linked differences that might underlie disease association. The average chromosome and pO157 plasmid size and gene content were significantly greater within the tir 255 A allele strains. Eighteen coding sequences were unique to tir 255 A allele chromosomes, and three were unique to tir 255 T allele chromosomes. There also were non-pO157 plasmids that were unique to each tir 255 allele variant. The overall average number of prophages did not differ between tir 255 allele strains; however, there were different types between the strains. Genomic and mobile element variation linked to the tir 255 polymorphism may account for the increased frequency of the T allele isolates in human disease.

Rates of evolutionary change of resident Escherichia coli O157:H7 differ within the same ecological niche.

BACKGROUND: Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a pathogen known to reside in cattle feedlots. This retrospective study examined 181 STEC O157:H7 strains collected over 23 years from a closed-system feedlot. All strains were subjected to short-read sequencing, with a subset of 36 also subjected to long-read sequencing. RESULTS: Over 96% of the strains fell into four phylogenetically distinct clades. Clade membership was associated with multiple factors including stx composition and the alleles of a well-characterized polymorphism (tir 255 T > A). Small plasmids (2.7 to 40 kb) were found to be primarily clade specific. Within each clade, chromosomal rearrangements were observed along with a core phageome and clade specific phages. Across both core and mobile elements of the genome, multiple SNP alleles were in complete linkage disequilibrium across all strains within specific clades. Clade evolutionary rates varied between 0.9 and 2.8 SNP/genome/year with two tir A allele clades having the lowest evolutionary rates. Investigation into possible causes of the differing rates was not conclusive but revealed a synonymous based mutation in the DNA polymerase III of the fastest evolving clade. Phylogenetic trees generated through our bioinformatic pipeline versus the NCBI's pathogen detection project were similar, with the two tir A allele clades matching individual NCBI SNP clusters, and the two tir T allele clades assigned to multiple closely-related SNP clusters. CONCLUSIONS: In one ecological niche, a diverse STEC O157:H7 population exhibited different rates of evolution that associated with SNP alleles in linkage disequilibrium in the core genome and mobile elements, including tir 255 T > A.

Twenty-Four-Month Longitudinal Study Suggests Little to No Horizontal Gene Transfer In Situ between Third-Generation Cephalosporin-Resistant Salmonella and Third-Generation Cephalosporin-Resistant Escherichia coli in a Beef Cattle Feedyard.

ABSTRACT: Third-generation cephalosporins (3GCs) are preferred treatments for serious human Salmonella enterica infections. Beef cattle are suspected to contribute to human 3GC-resistant Salmonella infections. Commensal 3GC-resistant Escherichia coli are thought to act as reservoirs of 3GC resistance because these strains are isolated more frequently than are 3GC-resistant Salmonella strains at beef cattle feedyards. During each of 24 consecutive months, four samples of pen surface material were obtained from five pens (N = 480) at a Nebraska feedyard to determine to the contribution of 3GC-resistant E. coli to the occurrence of 3GC-resistant Salmonella. Illumina whole genome sequencing was performed, and susceptibility to 14 antimicrobial agents was determined for 121 3GC-susceptible Salmonella, 121 3GC-resistant Salmonella, and 203 3GC-resistant E. coli isolates. 3GC-susceptible Salmonella isolates were predominantly from serotypes Muenchen (70.2%) and Montevideo clade 1 (23.1%). 3GC-resistant Salmonella isolates were predominantly from serotypes Montevideo clade 2 (84.3%). One bla gene type (blaCMY-2) and the IncC plasmid replicon were present in 100 and 97.5% of the 3GC-resistant Salmonella, respectively. Eleven bla gene types were detected in the 3GC-resistant E. coli, which were distributed across 42 multilocus sequence types. The blaCMY-2 gene and IncC plasmid replicon were present in 37.9 and 9.9% of the 3GC-resistant E. coli, respectively. These results suggest that 3GC resistance in Salmonella was primarily due the persistence of Salmonella Montevideo clade 2 with very minimal or no contribution from 3GC-resistant E. coli via horizontal gene transfer and that 3GC-resistant E. coli may not be a useful indicator for 3GC-resistant Salmonella in beef cattle production environments.

A Conserved Histophilus somni 23S Intervening Sequence Yields Functional, Fragmented 23S rRNA.

Histophilus somni is a Gram-negative bacterial organism that acts as an opportunistic pathogen and is a fastidious member of the Pasteurellaceae family associated with diseases of respiratory, reproductive, cardiac, and other tissues of ruminants. We identified an intervening sequence (IVS) embedded in all five copies of the 23S rRNA gene in the closed genome sequence of the H. somni isolate USDA-ARS-USMARC-63250 that may play an important role in affecting the biology of the organism. Sequencing the RNA from this isolate shows that most of the IVS is cleaved from the transcript, resulting in independent fragments of this structural rRNA that remain functional within the bacterial ribosome. The IVS lies between positions 1170 and 1278 bp of the 3,017-bp gene and exhibits self-complementarity between its 5' and 3' ends that predicts a stem-loop structure interrupting helix-45 in the transcribed 23S rRNA. Excision removes a 94-nucleotide (nt) stem-loop structure that displays an unusual 1-nt 3' end overhang instead of the more typical 2-nt overhang commonly observed at the ends of other excised IVS stem-loops. A comparison with genomes of other H. somni isolates indicates that this IVS is highly conserved, with 31 of 32 complete genomes having similar interruptions of canonical 23S rRNA genes. The potential biological effects of either the released IVS or the fragmentation of the functional 23S rRNA are unknown, but fragmentation may enhance rRNA degradation in ways that contribute to the regulation of gene expression. IMPORTANCE The genome biology underlying H. somni virulence, pathogenicity, environmental adaptability, and broad tissue tropism is understood poorly. We identified a novel H. somni 109-nt IVS stem-loop structure, of which the central portion is excised from the 23S rRNA transcript, resulting in the fragmentation of this rRNA in the H. somni isolate USDA-ARS-USMARC-63250 and the release of a 94-nt structured RNA of unknown function. We determined that this peculiar rRNA biology is widespread among sequenced H. somni isolates, suggesting it has importance to organism biology. The fragmented 23S rRNA molecules remain functional in the ribosome, given that the isolate grows in culture. The structured excised portion of the IVS, presumably due to the action of the endoribonuclease III, has an unusual 1-nt 3' end overhang. This newly discovered H. somni 23S rRNA fragmentation may enhance rRNA degradation providing a previously unrecognized avenue for regulating H. somni biological processes.

Locus of Heat Resistance (LHR) in Meat-Borne Escherichia coli: Screening and Genetic Characterization.

Microbial resistance to processing treatments poses a food safety concern, as treatment tolerant pathogens can emerge. Occasional foodborne outbreaks caused by pathogenic Escherichia coli have led to human and economic losses. Therefore, this study screened for the extreme heat resistance (XHR) phenotype as well as one known genetic marker, the locus of heat resistance (LHR), in 4,123 E. coli isolates from diverse meat animals at different processing stages. The prevalences of XHR and LHR among the meat-borne E. coli were found to be 10.3% and 11.4%, respectively, with 19% agreement between the two. Finished meat products showed the highest LHR prevalence (24.3%) compared to other processing stages (0 to 0.6%). None of the LHR+E. coli in this study would be considered pathogens based on screening for virulence genes. Four high-quality genomes were generated by whole-genome sequencing of representative LHR+ isolates. Nine horizontally acquired LHRs were identified and characterized, four plasmid-borne and five chromosomal. Nine newly identified LHRs belong to ClpK1 LHR or ClpK2 LHR variants sharing 61 to 68% nucleotide sequence identity, while one LHR appears to be a hybrid. Our observations suggest positive correlation between the number of LHR regions present in isolates and the extent of heat resistance. The isolate exhibiting the highest degree of heat resistance possessed four LHRs belonging to three different variant groups. Maintenance of as many as four LHRs in a single genome emphasizes the benefits of the LHR in bacterial physiology and stress response.IMPORTANCE Currently, a "multiple-hurdle" approach based on a combination of different antimicrobial interventions, including heat, is being utilized during meat processing to control the burden of spoilage and pathogenic bacteria. Our recent study (M. Guragain, G. E. Smith, D. A. King, and J. M. Bosilevac, J Food Prot 83:1438-1443, 2020, https://doi.org/10.4315/JFP-20-103) suggests that U.S. beef cattle harbor Escherichia coli that possess the locus of heat resistance (LHR). LHR seemingly contributes to the global stress tolerance in bacteria and hence poses a food safety concern. Therefore, it is important to understand the distribution of the LHRs among meat-borne bacteria identified at different stages of different meat processing systems. Complete genome sequencing and comparative analysis of selected heat-resistant bacteria provide a clearer understanding of stress and heat resistance mechanisms. Further, sequencing data may offer a platform to gain further insights into the genetic background that provides optimal bacterial tolerance against heat and other processing treatments.

Rapid estimation of Salmonella enterica contamination level in ground beef - Application of the time-to-positivity method using a combination of molecular detection and direct plating.

Little progress has been made in decreasing the incidence rate of salmonellosis in the US over the past decade. Mitigating the contribution of contaminated raw meat to the salmonellosis incidence rate requires rapid methods for quantifying Salmonella, so that highly contaminated products can be removed before entering the food chain. Here we evaluated the use of Time-to-Positivity (TTP) as a rapid, semi-quantitative approach for estimating Salmonella contamination levels in ground beef. Growth rates of 14 Salmonella strains (inoculated at log 1 to -2 CFU/g) were characterized in lean ground beef mTSB enrichments and time-to-detection was determined using culture and molecular detection methods. Enrichments were sampled at five timepoints and results were used to construct a prediction model of estimated contamination level by TTP (superscript indicates time in hours) defined as TTP4: ≥5 CFU/g; TTP6: ≤5, ≥1 CFU/g; TTP8: ≤1, ≥0.01 CFU/g; with samples negative at 8 h estimated ≤0.01 CFU/g. Model performance measures showed high sensitivity (100%) and specificity (83% and 93% for two detection methods) for samples with a TTP4, with false negative rates of 0%.

Consecutive Treatments with a Multicomponent Sanitizer Inactivate Biofilms Formed by Escherichia coli O157:H7 and Salmonella enterica and Remove Biofilm Matrix.

ABSTRACT: Many foodborne pathogens, including Escherichia coli O157:H7 and Salmonella enterica, can develop biofilms on contact surfaces at meat processing plants. Owing to the high tolerance of the biofilm cells associated with the three-dimensional biofilm structure and the well-expressed bacterial extracellular polymeric substances, it is a real challenge to completely inactivate and remove mature biofilms, as well as further prevent biofilm reoccurrence and pathogen survival. In the present study, we evaluated the effectiveness of consecutive treatments (10 to 120 min per treatment) by repeatedly applying a multicomponent sanitizer, based on a functional mechanism by synergistic combination of hydrogen peroxide and quaternary ammonia compounds, against biofilms formed by E. coli O157:H7 and S. enterica strains. Biofilms on stainless steel surfaces were treated with 2.5, 5, or 10% (recommended working concentration) of the sanitizer applied as a foam or liquid solution. Our results showed that the multicomponent sanitizer significantly (P < 0.05) reduced the amount of viable biofilm cells at all concentrations, as enumerable bacteria were only detected after low-concentration treatments (2.5 or 5%) with short exposure periods (10 or 20 min per treatment). Treatments with high concentrations (5 or 10%) of the sanitizer, multiple consecutive treatments (2 or 3 treatments), and sufficient exposure time (>60 min per treatment) effectively controlled pathogen survival postsanitization. Examination with a scanning electron microscope showed that treatment with the sanitizer at 5% strength significantly dissolved the connecting extracellular polysaccharide matrix and removed the majority of the biofilm matrix. No intact biofilm structure was detected after the 10% sanitizer treatment; instead, scattered individual bacteria with visibly altered cell morphology were observed. The treated bacteria exhibited indented and distorted shapes with shortened cell length and increased surface roughness, indicating severe cell injury and death. Our observations indicated that consecutive treatments with the multicomponent sanitizer was effective in inactivating E. coli O157:H7 and S. enterica biofilms and preventing pathogen reoccurrence.

Differences between predicted outer membrane proteins of genotype 1 and 2 Mannheimia haemolytica.

BACKGROUND: Mannheimia haemolytica strains isolated from North American cattle have been classified into two genotypes (1 and 2). Although members of both genotypes have been isolated from the upper and lower respiratory tracts of cattle with or without bovine respiratory disease (BRD), genotype 2 strains are much more frequently isolated from diseased lungs than genotype 1 strains. The mechanisms behind the increased association of genotype 2 M. haemolytica with BRD are not fully understood. To address that, and to search for interventions against genotype 2 M. haemolytica, complete, closed chromosome assemblies for 35 genotype 1 and 34 genotype 2 strains were generated and compared. Searches were conducted for the pan genome, core genes shared between the genotypes, and for genes specific to either genotype. Additionally, genes encoding outer membrane proteins (OMPs) specific to genotype 2 M. haemolytica were identified, and the diversity of their protein isoforms was characterized with predominantly unassembled, short-read genomic sequences for up to 1075 additional strains. RESULTS: The pan genome of the 69 sequenced M. haemolytica strains consisted of 3111 genes, of which 1880 comprised a shared core between the genotypes. A core of 112 and 179 genes or gene variants were specific to genotype 1 and 2, respectively. Seven genes encoding predicted OMPs; a peptidase S6, a ligand-gated channel, an autotransporter outer membrane beta-barrel domain-containing protein (AOMB-BD-CP), a porin, and three different trimeric autotransporter adhesins were specific to genotype 2 as their genotype 1 homologs were either pseudogenes, or not detected. The AOMB-BD-CP gene, however, appeared to be truncated across all examined genotype 2 strains and to likely encode dysfunctional protein. Homologous gene sequences from additional M. haemolytica strains confirmed the specificity of the remaining six genotype 2 OMP genes and revealed they encoded low isoform diversity at the population level. CONCLUSION: Genotype 2 M. haemolytica possess genes encoding conserved OMPs not found intact in more commensally prone genotype 1 strains. Some of the genotype 2 specific genes identified in this study are likely to have important biological roles in the pathogenicity of genotype 2 M. haemolytica, which is the primary bacterial cause of BRD.

Classification of 16S rRNA reads is improved using a niche-specific database constructed by near-full length sequencing.

Surveys of microbial populations in environmental niches of interest often utilize sequence variation in the gene encoding the ribosomal small subunit (the 16S rRNA gene). Generally, these surveys target the 16S genes using semi-degenerate primers to amplify portions of a subset of bacterial species, sequence the amplicons in bulk, and assign to putative taxonomic categories by comparison to databases purporting to connect specific sequences in the main variable regions of the gene to specific organisms. Due to sequence length constraints of the most popular bulk sequencing platforms, the primers selected amplify one to three of the nine variable regions, and taxonomic assignment is based on relatively short stretches of sequence (150-500 bases). We demonstrate that taxonomic assignment is improved through reduced unassigned reads by including a survey of near-full-length sequences specific to the target environment, using a niche of interest represented by the upper respiratory tract (URT) of cattle. We created a custom Bovine URT database from these longer sequences for assignment of shorter, less expensive reads in comparisons of the upper respiratory tract among individual animals. This process improves the ability to detect changes in the microbial populations of a given environment, and the accuracy of defining the content of that environment at increasingly higher taxonomic resolution.

Effectiveness and Functional Mechanism of a Multicomponent Sanitizer against Biofilms Formed by Escherichia coli O157:H7 and Five Salmonella Serotypes Prevalent in the Meat Industry.

ABSTRACT: Biofilm formation by Escherichia coli O157:H7 and Salmonella enterica at meat processing plants poses a potential risk of meat product contamination. Many common sanitizers are unable to completely eradicate biofilms formed by these foodborne pathogens because of the three-dimensional biofilm structure and the presence of bacterial extracellular polymeric substances (EPSs). A novel multifaceted approach combining multiple chemical reagents with various functional mechanisms was used to enhance the effectiveness of biofilm control. We tested a multicomponent sanitizer consisting of a quaternary ammonium compound (QAC), hydrogen peroxide, and the accelerator diacetin for its effectiveness in inactivating and removing Escherichia coli O157:H7 and Salmonella enterica biofilms under meat processing conditions. E. coli O157:H7 and Salmonella biofilms on common contact surfaces were treated with 10, 20, or 100% concentrations of the multicomponent sanitizer solution for 10 min, 1 h, or 6 h, and log reductions in biofilm mass were measured. Scanning electron microscopy (SEM) was used to directly observe the effect of sanitizer treatment on biofilm removal and bacterial morphology. After treatment with the multicomponent sanitizer, viable E. coli O157:H7 and Salmonella biofilm cells were below the limit of detection, and the prevalence of both pathogens was low. After treatment with a QAC-based control sanitizer, surviving bacterial cells were countable, and pathogen prevalence was higher. SEM analysis of water-treated control samples revealed the three-dimensional biofilm structure with a strong EPS matrix connecting bacteria and the contact surface. Treatment with 20% multicomponent sanitizer for 10 min significantly reduced biofilm mass and weakened the EPS connection. The majority of the bacterial cells had altered morphology and compromised membrane integrity. Treatment with 100% multicomponent sanitizer for 10 min dissolved the EPS matrix, and no intact biofilm structure was observed; instead, scattered clusters of bacterial aggregates were detected, indicating the loss of cell viability and biofilm removal. These results indicate that the multicomponent sanitizer is effective, even after short exposure with dilute concentrations, against E. coli O157:H7 and S. enterica biofilms.

Publisher Correction: A Computational Method to Quantify the Effects of Slipped Strand Mispairing on Bacterial Tetranucleotide Repeats.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A Computational Method to Quantify the Effects of Slipped Strand Mispairing on Bacterial Tetranucleotide Repeats.

The virulence and pathogenicity of bacterial pathogens are related to their adaptability to changing environments. One process enabling adaptation is based on minor changes in genome sequence, as small as a few base pairs, within segments of genome called simple sequence repeats (SSRs) that consist of multiple copies of a short sequence (from one to several nucleotides), repeated in series. SSRs are found in eukaryotes as well as prokaryotes, and length variation in them occurs at frequencies up to a million-fold higher than bacterial point mutations through the process of slipped strand mispairing (SSM) by DNA polymerase during replication. The characterization of SSR length by standard sequencing methods is complicated by the appearance of length variation introduced during the sequencing process that obscures the lower abundance repeat number variants in a population. Here we report a computational approach to correct for sequencing process-induced artifacts, validated for tetranucleotide repeats by use of synthetic constructs of fixed, known length. We apply this method to a laboratory culture of Histophilus somni, prepared from a single colony, and demonstrate that the culture consists of populations of distinct sequence phase and length variants at individual tetranucleotide SSR loci.

Complete Genome Sequence of a Salmonella enterica subsp. enterica Serovar Fresno Isolate Recovered from a Bovine Lymph Node.

Salmonella enterica serovar Fresno is an infrequently isolated serovar whose ecology and genomic characteristics have not yet been described. To further understand the genomic characteristics of this serovar, we sequenced the complete genome of a single isolate recovered from a bovine lymph node at harvest.

Complete Genome Sequences of Four Salmonella enterica Strains (Including Those of Serotypes Montevideo, Mbandaka, and Lubbock) Isolated from Peripheral Lymph Nodes of Healthy Cattle.

Salmonella enterica serotype Lubbock emerged most likely from a Salmonella enterica serotype Mbandaka ancestor that acquired by recombination the fliC operon from Salmonella enterica serotype Montevideo. Here, we report the complete genome sequence of two S. Lubbock, one S. Montevideo, and one S. Mbandaka strain isolated from bovine lymph nodes.

Closed Genome Sequences and Antibiograms of 16 Pasteurella multocida Isolates from Bovine Respiratory Disease Complex Cases and Apparently Healthy Controls.

Pasteurella multocida is an animal-associated Gram-negative member of the Pasteurellaceae family. It is an opportunistic pathogen and is one of the principal bacterial species contributing to bovine respiratory disease complex (BRDC) in feedlot cattle. We present 16 closed genome sequences and antibiograms of isolates cultured from calves exhibiting clinical signs of BRDC and from control calves not showing signs of BRDC.

Complete Closed Genome Sequences of Three Salmonella enterica subsp. enterica Serovar Dublin Strains Isolated from Cattle at Harvest.

Salmonella enterica subsp. enterica serovar Dublin is a host-adapted pathogen for cattle that can cause invasive disease in humans. To facilitate genomic comparisons characterizing virulence determinants of this pathogen, we present the complete genome sequences of three S. Dublin strains isolated from bovine sources at harvest.

Whole-Genome Sequence Analysis of CTX-M Containing Escherichia coli Isolates from Retail Meats and Cattle in the United States.

In recent years, there have been increased reports on the detection of extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli and Salmonella strains from food-producing animals and animal products in the United States. We characterized 18 ESBL E. coli isolates from cattle (n = 5), chicken breast (n = 5), ground turkey (n = 6), ground beef (n = 1), and pork chops (n = 1) that were collected by the National Antimicrobial Resistance Monitoring System (NARMS) between 2011 and 2015. In vitro antimicrobial susceptibility testing was done against a panel of 14 antimicrobials followed by a secondary panel of 9 ß-lactam agents. Whole-genome sequencing was used to characterize the resistome, plasmids, and the genetic structures of the ESBL genes. All ESBL-producing E. coli isolates were resistant to at least three antimicrobial classes and carried various blaCTX-M genes. Most of the cattle and ground turkey isolates carried blaCTX-M-27. In chicken breast isolates, blaCTX-M-1 was present as part of an ISEcp1 transposition unit carried on a plasmid that shares sequence similarity with the backbone structure of the IncI plasmid. Isolates carrying the blaCTX-M-14 and blaCTX-M-15 genes, widely distributed in human clinical isolates, were also isolated. To our knowledge, this is the first report of the widely distributed blaCTX-M-14 and blaCTX-M-15 in E. coli isolates from retail meat samples in the United States. Different insertional sequences were identified upstream of these blaCTX-Ms, including ISEcp1, IS26, and IS903-D. CTX-M in E. coli from food animals and retail chicken breast were often present on plasmids with other resistance genes. Other resistance genes identified included aadA, strA, strB, aac(3)-IId, aac(3)-VIa, aph(3')-Ic, blaTEM, blaHERA-3, floR, sul1, sul2, catA1, tetA, tetB, dfrA, and qacE. These data describe the emergence of CTX-M-carrying E. coli isolates in food animals and animal products monitored by NARMS program.

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association.

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.