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

Intraspecies competition among Salmonella enterica isolates in the lettuce leaf apoplast

Multiple Salmonella enterica serovars and strains have been reported to be able to persist inside the foliar tissue of lettuce (Lactuca sativa L.), potentially resisting washing steps and reaching the consumer. Intraspecies variation of the bacterial pathogen and of the plant host can both significantly affect the outcome of foliar colonization. However, current understanding of the mechanisms underlying this phenomenon is still very limited. In this study, we evaluated the foliar fitness of 14 genetically barcoded S. enterica isolates from 10 different serovars, collected from plant and animal sources. The S. enterica isolates were vacuum-infiltrated individually or in pools into the leaves of three- to four-week-old lettuce plants. To estimate the survival capacity of individual isolates, we enumerated the bacterial populations at 0- and 10- days post-inoculation (DPI) and calculated their net growth. The competition of isolates in the lettuce apoplast was assessed through the determination of the relative abundance change of barcode counts of each isolate within pools during the 10 DPI experimental period. Isolates exhibiting varying apoplast fitness phenotypes were used to evaluate their capacity to grow in metabolites extracted from the lettuce apoplast and to elicit the reactive oxygen species burst immune response. Our study revealed that strains of S. enterica can substantially differ in their ability to survive and compete in a co-inhabited lettuce leaf apoplast. The differential foliar fitness observed among these S. enterica isolates might be explained, in part, by their ability to utilize nutrients available in the apoplast and to evade plant immune responses in this niche.

Genomic and phenotypic characterization of Salmonella enterica serovar Kentucky

Non-typhoidal Salmonella are extremely diverse and different serovars can exhibit varied phenotypes, including host adaptation and the ability to cause clinical illness in animals and humans. In the USA, Salmonella enterica serovar Kentucky is infrequently found to cause human illness, despite being the top serovar isolated from broiler chickens. Conversely, in Europe, this serovar falls in the top 10 serovars linked to human salmonellosis. Serovar Kentucky is polyphyletic and has two lineages, Kentucky-I and Kentucky-II; isolates belonging to Kentucky-I are frequently isolated from poultry in the USA, while Kentucky-II isolates tend to be associated with human illness. In this study, we analysed whole-genome sequences and associated metadata deposited in public databases between 2017 and 2021 by federal agencies to determine serovar Kentucky incidence across different animal and human sources. Of 5151 genomes, 90.3 % were from isolates that came from broilers, while 5.9 % were from humans and 3.0 % were from cattle. Kentucky-I isolates were associated with broilers, while isolates belonging to Kentucky-II and a new lineage, Kentucky-III, were more commonly associated with cattle and humans. Very few serovar Kentucky isolates were associated with turkey and swine sources. Phylogenetic analyses showed that Kentucky-III genomes were more closely related to Kentucky-I, and this was confirmed by CRISPR-typing and multilocus sequence typing (MLST). In a macrophage assay, serovar Kentucky-II isolates were able to replicate over eight times better than Kentucky-I isolates. Analysis of virulence factors showed unique patterns across these three groups, and these differences may be linked to their association with different hosts.

Combined Quantification and Deep Serotyping for Salmonella Risk Profiling in Broiler Flocks

Despite a reduction of Salmonella contamination on final poultry products, the level of human salmonellosis cases attributed to poultry has remained unchanged over the last few years. There needs to be improved effort to target serovars which may survive antimicrobial interventions and cause illness, as well as to focus on lessening the amount of contamination entering the processing plant. Advances in molecular enumeration approaches allow for the rapid detection and quantification of Salmonella in pre- and postharvest samples, which can be combined with deep serotyping to properly assess the risk affiliated with a poultry flock. In this study, we collected a total of 160 boot sock samples from 20 broiler farms across four different integrators with different antibiotic management programs. Overall, Salmonella was found in 85% (68/80) of the houses, with each farm having at least one Salmonella-positive house. The average Salmonella quantity across all four complexes was 3.6 log10 CFU/sample. Eleven different serovars were identified through deep serotyping, including all three key performance indicators (KPIs; serovars Enteritidis, Infantis, and Typhimurium) defined by the U.S. Department of Agriculture-Food Safety and Inspection Service (USDA-FSIS). There were eight multidrug resistant isolates identified in this study, and seven which were serovar Infantis. We generated risk scores for each flock based on the presence or absence of KPIs, the relative abundance of each serovar as calculated with CRISPR-SeroSeq (serotyping by sequencing the clustered regularly interspaced palindromic repeats), and the quantity of Salmonella organisms detected. The work presented here provides a framework to develop directed processing approaches and highlights the limitations of conventional Salmonella sampling and culturing methods. IMPORTANCE Nearly one in five foodborne Salmonella illnesses are derived from chicken, making it the largest single food category to cause salmonellosis and indicating a need for effective pathogen mitigation. Although industry has successfully reduced Salmonella incidence in poultry products, there has not been a concurrent reduction in human salmonellosis linked to chicken consumption. New efforts are focused on improved control at preharvest, which requires improved Salmonella surveillance. Here, we present a high-resolution surveillance approach that combines quantity and identity of Salmonella in broiler flocks prior to processing which will further support improved Salmonella controls in poultry. We developed a framework for this approach, indicating that it is possible and important to harness deep serotyping and molecular enumeration to inform on-farm management practices and to minimize risk of cross-contamination between flocks at processing. Additionally, this framework could be adapted to Salmonella surveillance in other food animal production systems.

The Number and Type of Chaperone-Usher Fimbriae Reflect Phylogenetic Clade Rather than Host Range in Salmonella

Salmonella is one of the most successful foodborne pathogens worldwide, owing in part to its ability to colonize or infect a wide range of hosts. Salmonella serovars are known to encode a variety of different fimbriae (hairlike organelles that facilitate binding to surfaces); however, the distribution, number, and sequence diversity of fimbriae encoded across different lineages of Salmonella were unknown. We queried whole-genome sequence (WGS) data for 242 Salmonella enterica subsp. enterica (subspecies enterica) isolates from the top 217 serovars associated with isolation from humans and agricultural animals; this effort identified 2,894 chaperone-usher (CU)-type fimbrial usher sequences, representing the most conserved component of CU fimbriae. On average, isolates encoded 12 different CU fimbrial ushers (6 to 18 per genome), although the distribution varied significantly (P = 1.328E-08) by phylogenetic clade, with isolates in section Typhi having significantly fewer fimbrial ushers than isolates in clade A2 (medians = 10 and 12 ushers, respectively). Characterization of fimbriae in additional non-enterica subspecies genomes suggested that 8 fimbrial ushers were classified as being unique to subspecies enterica isolates, suggesting that the majority of fimbriae were most likely acquired prior to the divergence of subspecies enterica. Characterization of mobile elements suggested that plasmids represent an important vehicle facilitating the acquisition of a wide range of fimbrial ushers, particularly for the acquisition of fimbriae from other Gram-negative genera. Overall, our results suggest that differences in the number and type of fimbriae encoded most likely reflect differences in phylogenetic clade rather than differences in host range. IMPORTANCE Fimbriae of the CU assembly pathway represent important organelles that mediate Salmonella's interactions with host tissues and abiotic surfaces. Our analyses provide a comprehensive overview of the diversity of CU fimbriae in Salmonella spp., highlighting that the majority of CU fimbriae are distributed broadly across multiple subspecies and suggesting that acquisition most likely occurred prior to the divergence of subspecies enterica. Our data also suggest that plasmids represent the primary vehicles facilitating the horizontal transfer of diverse CU fimbriae in Salmonella. Finally, the observed high sequence similarity between some ushers suggests that different names may have been assigned to closely related fimbrial ushers that likely should be represented by a single designation. This highlights the need to establish standard criteria for fimbria classification and nomenclature, which will also facilitate future studies seeking to associate virulence factors with adaptation to or differences in the likelihood of causing disease in a given host.

Application of a High-Throughput Targeted Sequence AmpliSeq Procedure to Assess the Presence and Variants of Virulence Genes in Salmonella

We have developed a targeted, amplicon-based next-generation sequencing method to detect and analyze 227 virulence genes (VG) of Salmonella (AmpliSeq_{Salm_227VG}) for assessing the pathogenicity potential of Salmonella. The procedure was developed using 80 reference genomes representing 75 epidemiologically-relevant serovars associated with human salmonellosis. We applied the AmpliSeq_{Salm_227VG} assay to (a) 35 previously characterized field strains of Salmonella consisting of serovars commonly incriminated in foodborne illnesses and (b) 34 Salmonella strains with undisclosed serological or virulence attributes, and were able to divide Salmonella VGs into two groups: core VGs and variable VGs. The commonest serovars causing foodborne illnesses such as Enteritidis, Typhimurium, Heidelberg and Newport had a high number of VGs (217–227). In contrast, serovars of subspecies not commonly associated with human illnesses, such as houtenae, arizonae and salame, tended to have fewer VGs (177–195). Variable VGs were not only infrequent but, when present, displayed considerable sequence variation: safC, sseL, sseD, sseE, ssaK and stdB showed the highest variation and were linked to strain pathogenicity. In a chicken infection model, VGs belonging to rfb and sse operons showed differences and were linked with pathogenicity. The high-throughput, targeted NGS-based AmpliSeq_{Salm_227VG} procedure provided previously unknown information about variation in select virulence genes that can now be applied to a much larger population of Salmonella for evaluating pathogenicity of various serovars of Salmonella and for risk assessment of foodborne salmonellosis.

Differences in biofilm formation of Salmonella serovars on two surfaces under two temperature conditions

AIMS

Salmonella is extremely diverse, with >2500 serovars that are genetically and phenotypically diverse. The aim of this study was to build a collection of Salmonella isolates that are genetically diverse and to evaluate their ability to form biofilm under different conditions relevant to a processing environment.

METHODS AND RESULTS

Twenty Salmonella isolates representative of 10 serovars were subtyped using Clustered regularly interspaced short palindromic repeats (CRISPR)-typing to assess the genetic diversity between isolates of each serovar. Biofilm formation of the isolates on both plastic and stainless-steel surfaces at 25 and 15°C was assessed. At 25°C, 8/20 isolates each produced strong and moderate biofilm on plastic surface compared to stainless-steel (3/20 and 13/20 respectively). At 15°C, 5/20 produced strong biofilm on plastic surface and none on stainless-steel. Several isolates produced weak biofilm on plastic (11/20) and stainless-steel (16/20) surfaces. Serovar Schwarzengrund consistently produced strong biofilm while serovars Heidelberg and Newport produced weak biofilm.

CONCLUSION

These results suggest that Salmonellae differ in their attachment depending on the surface and temperature conditions encountered, which may influence persistence in the processing environment.

SIGNIFICANCE AND IMPACT OF STUDY

These differences in biofilm formation could provide useful information for mitigation of Salmonella in processing environments.

Phylogeographic Clustering Suggests that Distinct Clades of Salmonella enterica Serovar Mississippi Are Endemic in Australia, the United Kingdom, and the United States

Salmonella enterica serovar Mississippi is the 2nd and 14th leading cause of human clinical salmonellosis in the Australian island state of Tasmania and the United States, respectively. Despite its public health relevance, relatively little is known about this serovar. Comparison of whole-genome sequence (WGS) data of S. Mississippi isolates with WGS data for 317 additional S. enterica serovars placed one clade of S. Mississippi within S. enterica clade B (“clade B Mississippi”) and the other within section Typhi in S. enterica clade A (“clade A Mississippi”), suggesting that these clades evolved from different ancestors. Phylogenetic analysis of 364 S. Mississippi isolates from Australia, the United Kingdom, and the United States suggested that the isolates cluster geographically, with U.S. and Australian isolates representing different subclades (Ai and Aii, respectively) within clade A Mississippi and clade B isolates representing the predominant S. Mississippi isolates in the United Kingdom. Intraclade comparisons suggested that different mobile elements, some of which encode virulence factors, are responsible for the observed differences in gene content among isolates within these clades. Specifically, genetic differences among clade A isolates reflect differences in prophage contents, while differences among clade B isolates are due to the acquisition of a 47.1-kb integrative conjugative element (ICE). Phylogenies inferred from antigenic components (fliC, fljB, and O-antigen-processing genes) support that clade A and B Mississippi isolates acquired these loci from different ancestral serovars. Overall, these data support that different S. Mississippi phylogenetic clades are endemic in Australia, the United Kingdom, and the United States.

IMPORTANCE The number of known so-called “polyphyletic” serovars (i.e., phylogenetically distinct clades with the same O and H antigenic formulas) continues to increase as additional Salmonella isolates are sequenced. While serotyping remains a valuable tool for reporting and monitoring Salmonella, more discriminatory analyses for classifying polyphyletic serovars may improve surveillance efforts for these serovars, as we found that for S. Mississippi, distinct genotypes predominate at different geographic locations. Our results suggest that the acquisition of genes encoding O and H antigens from different ancestors led to the emergence of two Mississippi clades. Furthermore, our results suggest that different mobile elements contribute to the microevolution and diversification of isolates within these two clades, which has implications for the acquisition of novel adaptations, such as virulence factors.

Antimicrobial Resistance and Genomic Characterization of Two mcr-1-Harboring Foodborne Salmonella Isolates Recovered in China, 2016

The mcr-1 gene mediating mobile colistin resistance in Escherichia coli was first reported in China in 2016 followed by reports among different species worldwide, especially in E. coli and Klebsiella. However, data on its transmission in Salmonella are still lacking. This study analyzed the antimicrobial resistance (AMR) profiles and the mcr-1 gene presence in 755 foodborne Salmonella from 26 provinces of mainland, China in 2016. Genomic features of two mcr-1-carrying isolates, genome sequencing, serotypes and further resistance profiles were studied. Among the 755 Salmonella tested, 72.6% were found to be resistant to at least one antimicrobial agent and 10% were defined as multi-drug resistant (MDR). Salmonella Derby CFSA231 and Salmonella Typhimurium CFSA629 were mcr-1-harboring isolates. Both expressed an MDR phenotype and included a single circular chromosome and one plasmid. Among the 22 AMR genes identified in S. Derby CFSA231, only the mcr-1 gene was localized on the IncX4 type plasmid pCFSA231 while 20 chromosomal AMR genes, including four plasmid-mediated quinolone resistance (PMQR) genes, were mapped within a 64 kb Salmonella genomic island (SGI) like region. S. Typhimurium CFSA629 possessed 11 resistance genes including an mcr-1.19 variant and two ESBL genes. Two IS26-flanked composite-like transposons were identified. Additionally, 153 and 152 virulence factors were separately identified in these two isolates with secretion system and fimbrial adherence determinants as the dominant virulence classes. Our study extends our concern on mcr-1-carrying Salmonella in regards to antimicrobial resistance and virulence factors, and highlight the importance of surveillance to mitigate dissemination of mcr-encoding genes among foodborne Salmonella.

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. IMPORTANCE Salmonella 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.

Foodborne pathogens in the omics era

Outbreaks and deaths related to Foodborne Diseases (FBD) occur constantly in the world, as a result of the consumption of contaminated foodstuffs with pathogens such as Listeria monocytogenes, Escherichia coli, Staphylococcus aureus, Salmonella spp, Clostridium spp. and Campylobacter spp. The purpose of this review is to discuss the main omic techniques applied in foodborne pathogen and to demonstrate their functionalities through the food chain and to guarantee the food safety. The main techniques presented are genomic, transcriptomic, secretomic, proteomic, and metabolomic, which together, in the field of food and nutrition, are known as "Foodomics." This review had highlighted the potential of omics to integrate variables that contribute to food safety and to enable us to understand their application on foodborne diseases. The appropriate use of these techniques had driven the definition of critical parameters to achieve successful results in the improvement of consumers health, costs and to obtain safe and high-quality products.

Phylogeny of Salmonella enterica subspecies arizonae by whole-genome sequencing reveals high incidence of polyphyly and low phase 1 H antigen variability

Salmonella enterica subspecies arizonae is frequently associated with animal reservoirs, particularly reptiles, and can cause illness in some mammals, including humans. Using whole-genome sequencing data, core genome phylogenetic analyses were performed using 112 S. enterica subsp. arizonae isolates, representing 46 of 102 described serovars. Nearly one-third of these are polyphyletic, including two serovars that appear in four and five distinct evolutionary lineages. Subspecies arizonae has a monophasic H antigen. Among the 46 serovars investigated, only 8 phase 1 H antigens were identified, demonstrating high conservation for this antigen. Prophages and plasmids were found throughout this subspecies, including five novel prophages. Polyphyly was also reflected in prophage content, although some clade-specific enrichment for some phages was observed. IncFII(S) was the most frequent plasmid replicon identified and was found in a quarter of S. enterica subsp. arizonae genomes. Salmonella pathogenicity islands (SPIs) 1 and 2 are present across all Salmonella, including this subspecies, although effectors sipA, sptP and arvA in SPI-1 and sseG and ssaI in SPI-2 appear to be lost in this lineage. SPI-20, encoding a type VI secretion system, is exclusive to this subspecies and is well maintained in all genomes sampled. A number of fimbral operons were identified, including the sas operon that appears to be a synapomorphy for this subspecies, while others exhibited more clade-specific patterns. This work reveals evolutionary patterns in S. enterica subsp. arizonae that make this subspecies a unique lineage within this very diverse species.

The phylogenomics of CRISPR-Cas system and revelation of its features in Salmonella

Salmonellae display intricate evolutionary patterns comprising over 2500 serovars having diverse pathogenic profiles. The acquisition and/or exchange of various virulence factors influences the evolutionary framework. To gain insights into evolution of Salmonella in association with the CRISPR-Cas genes we performed phylogenetic surveillance across strains of 22 Salmonella serovars. The strains differed in their CRISPR1-leader and cas operon features assorting into two main clades, CRISPR1-STY/cas-STY and CRISPR1-STM/cas-STM, comprising majorly typhoidal and non-typhoidal Salmonella serovars respectively. Serovars of these two clades displayed better relatedness, concerning CRISPR1-leader and cas operon, across genera than between themselves. This signifies the acquisition of CRISPR1/Cas region could be through a horizontal gene transfer event owing to the presence of mobile genetic elements flanking CRISPR1 array. Comparison of CRISPR and cas phenograms with that of multilocus sequence typing (MLST) suggests differential evolution of CRISPR/Cas system. As opposed to broad-host-range, the host-specific serovars harbor fewer spacers. Mapping of protospacer sources suggested a partial correlation of spacer content with habitat diversity of the serovars. Some serovars like serovar Enteritidis and Typhimurium that inhabit similar environment/infect similar hosts hardly shared their protospacer sources.

Iron-Uptake Systems of Chicken-Associated Salmonella Serovars and Their Role in Colonizing the Avian Host

Iron is an essential micronutrient for most bacteria. Salmonella enterica strains, representing human and animal pathogens, have adopted several mechanisms to sequester iron from the environment depending on availability and source. Chickens act as a major reservoir for Salmonella enterica strains which can lead to outbreaks of human salmonellosis. In this review article we summarize the current understanding of the contribution of iron-uptake systems to the virulence of non-typhoidal S. enterica strains in colonizing chickens. We aim to address the gap in knowledge in this field, to help understand and define the interactions between S. enterica and these important hosts, in comparison to mammalian models.

Whole-Genome-Based Survey for Polyphyletic Serovars of Salmonella enterica subsp. enterica Provides New Insights into Public Health Surveillance

Serotyping has traditionally been considered the basis for surveillance of Salmonella, but it cannot distinguish distinct lineages sharing the same serovar that vary in host range, pathogenicity and epidemiology. However, polyphyletic serovars have not been extensively investigated. Public health microbiology is currently being transformed by whole-genome sequencing (WGS) data, which promote the lineage determination using a more powerful and accurate technique than serotyping. The focus in this study is to survey and analyze putative polyphyletic serovars. The multi-locus sequence typing (MLST) phylogenetic analysis identified four putative polyphyletic serovars, namely, Montevideo, Bareilly, Saintpaul, and Muenchen. Whole-genome-based phylogeny and population structure highlighted the polyphyletic nature of Bareilly and Saintpaul and the multi-lineage nature of Montevideo and Muenchen. The population of these serovars was defined by extensive genetic diversity, the open pan genome and the small core genome. Source niche metadata revealed putative existence of lineage-specific niche adaptation (host-preference and environmental-preference), exhibited by lineage-specific genomic contents associated with metabolism and transport. Meanwhile, differences in genetic profiles relating to virulence and antimicrobial resistance within each lineage may contribute to pathogenicity and epidemiology. The results also showed that recombination events occurring at the H1-antigen loci may be an important reason for polyphyly. The results presented here provide the genomic basis of simple, rapid, and accurate identification of phylogenetic lineages of these serovars, which could have important implications for public health.

Complete Genome Sequences of Four Salmonella enterica Strains Associated with Pistachios Assembled Using a Combination of Short- and Long-Read Sequencing

Here, we report the genomes of two Salmonella enterica subsp. enterica serovar Montevideo strains (CFSAN005645 and FCC0123) and two Salmonella enterica subsp. enterica serovar Senftenberg strains (FSW0104 and CFSAN087304) isolated from pistachios. The genomes were closed using a hybrid assembly method using short- and long-read sequencing technology.

Here, we report the genomes of two Salmonella enterica subsp. enterica serovar Montevideo strains (CFSAN005645 and FCC0123) and two Salmonella enterica subsp. enterica serovar Senftenberg strains (FSW0104 and CFSAN087304) isolated from pistachios. The genomes were closed using a hybrid assembly method using short- and long-read sequencing technology.