Characterization and evolution of Salmonella CRISPR-Cas systems

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

High-resolution Serotyping Reveals Salmonella Surveillance Challenges in the Turkey Industry

Despite extensive Salmonella controls used at processing, 5.5% of salmonellosis cases are linked to turkey. This study had two objectives: (i) to summarize USDA-FSIS turkey Salmonella verification program data and (ii) to evaluate Salmonella through turkey production and processing of 22 flocks. In objective 1, USDA-FSIS data show the average Salmonella prevalence in ground turkey from 2016 to 2022 was 15.9%, and that the leading serovar changes frequently. For objective 2, bootsocks (n = 22) were collected on-farm right after load-out. At processing, prescald wingtips (n = 6 composites of 10/flock), prechill wingtips (n = 6 composites of 10/flock), mechanically separated turkey (MST; n = 6 bins/flock), and ground turkey (n = 6 bins/flock) were collected. Salmonella prevalence was determined by a commercial qPCR and culture confirmed. In 33.2% of PCR-positive samples, Salmonella was not confirmed by culture, highlighting a discrepancy between molecular and culture detection. On-farm, 8/22 flocks were Salmonella positive, compared to 21 flocks that were positive at one or more processing locations, including 18 flocks that were positive in at least one final product sample. A logistic regression showed higher Salmonella prevalence in prescald (53.8%) than in prechill (18.2%), MST (27.3%) or ground turkey (26.5%). CRISPR-SeroSeq analysis of 148 culture-positive samples detected 18 Salmonella serovars and showed 35.1% of samples contained multiple serovars. In 16 flocks, one or more serovars detected in final products were absent from any upstream samples. Two-thirds of final product samples containing serovar Typhimurium typed as a live-attenuated Typhimurium vaccine strain. Salmonella on-farm and at prescald did not reflect Salmonella observed in final product. These data underscore the complexity of serovar tracking in turkey production and highlight challenges to identify surveillance samples that accurately represent Salmonella in turkey products.

Defence systems and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) is a fundamental process in prokaryotic evolution, contributing significantly to diversification and adaptation. HGT is typically facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages, which often impose fitness costs on their hosts. However, a considerable number of bacterial genes are involved in defence mechanisms that limit the propagation of MGEs, suggesting they may actively restrict HGT. In our study, we investigated whether defence systems limit HGT by examining the relationship between the HGT rate and the presence of 73 defence systems across 12 bacterial species. We discovered that only six defence systems, three of which were different CRISPR-Cas subtypes, were associated with a reduced gene gain rate at the species evolution scale. Hosts of these defence systems tend to have a smaller pangenome size and fewer phage-related genes compared to genomes without these systems. This suggests that these defence mechanisms inhibit HGT by limiting prophage integration. We hypothesize that the restriction of HGT by defence systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and the fitness effect of HGT in bacterial populations.

Multireceptor phage cocktail against Salmonella enterica to circumvent phage resistance

Non-Typhoidal Salmonella (NTS) is one of the most common food-borne pathogens worldwide, with poultry products being the major vehicle for pathogenesis in humans. The use of bacteriophage (phage) cocktails has recently emerged as a novel approach to enhancing food safety. Here, a multireceptor Salmonella phage cocktail of five phages was developed and characterized. The cocktail targets four receptors: O-antigen, BtuB, OmpC, and rough Salmonella strains. Structural analysis indicated that all five phages belong to unique families or subfamilies. Genome analysis of four of the phages showed they were devoid of known virulence or antimicrobial resistance factors, indicating enhanced safety. The phage cocktail broad antimicrobial spectrum against Salmonella, significantly inhibiting the growth of all 66 strains from 20 serovars tested in vitro. The average bacteriophage insensitive mutant (BIM) frequency against the cocktail was 6.22 × 10-6 in S. Enteritidis, significantly lower than that of each of the individual phages. The phage cocktail reduced the load of Salmonella in inoculated chicken skin by 3.5 log10 CFU/cm2 after 48 h at 25°C and 15°C, and 2.5 log10 CFU/cm2 at 4°C. A genome-wide transduction assay was used to investigate the transduction efficiency of the selected phage in the cocktail. Only one of the four phages tested could transduce the kanamycin resistance cassette at a low frequency comparable to that of phage P22. Overall, the results support the potential of cocktails of phage that each target different host receptors to achieve complementary infection and reduce the emergence of phage resistance during biocontrol applications.

Poultry Processing Interventions Reduce Salmonella Serovar Complexity on Postchill Young Chicken Carcasses as Determined by Deep Serotyping

Nearly 20% of salmonellosis cases are attributed to broilers, with renewed efforts to reduce Salmonella during broiler production and processing. A limitation to Salmonella culture is that often a single colony is picked for characterization, favoring isolation of the most abundant serovar found in a sample, while low abundance serovars can remain undetected. We used a deep serotyping approach, CRISPR-SeroSeq (serotyping by sequencing the clustered regularly interspaced palindromic repeats), to assess Salmonella serovar complexity during broiler processing and to determine the impact of antimicrobial interventions upon serovar population dynamics. Paired hot rehang and postchill young chicken carcasses were collected from establishments across the United States from August to November 2022. CRISPR-SeroSeq was performed on Salmonella culture-positive hot rehang (n = 153) and postchill (n = 38) samples, including 31 paired hot rehang and postchill samples. Multiple serovars were detected in 48.4% (74/153) and 7.9% (3/38) of hot rehang and postchill samples, respectively. On average, hot rehang carcasses contained 1.6 serovars, compared to 1.1 serovars at postchill (Mann Whitney U, p = 0.00018). Nineteen serovars were identified with serovar Kentucky the most common at hot rehang (72.5%; 111/153) and postchill (73.7%; 28/38). Serovar Infantis prevalence was higher at hot rehang (39.9%; 61/153) than in postchill (7.9%; 3/38). At hot rehang, serovar Enteritidis was outnumbered by other serovars 81.3% (13/16) of the time but was always the single or most abundant serovar detected when it was present at postchill (n = 5). We observed 98.4% (188/191) concordance between traditional isolation with serotyping and CRISPR-SeroSeq. Deep serotyping was able to explain serovar discrepancies between paired hot rehang and postchill samples when only traditional isolation and serotyping methods were used. These data demonstrate that processing interventions are effective in reducing Salmonella serovar complexity.

Defense systems and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) is a fundamental process in the evolution of prokaryotes, making major contributions to diversification and adaptation. Typically, HGT is facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages that generally impose fitness costs on their hosts. However, a substantial fraction of bacterial genes is involved in defense mechanisms that limit the propagation of MGEs, raising the possibility that they can actively restrict HGT. Here we examine whether defense systems curb HGT by exploring the connections between HGT rate and the presence of 73 defense systems in 12 bacterial species. We found that only 6 defense systems, 3 of which are different CRISPR-Cas subtypes, are associated with the reduced gene gain rate on the scale of species evolution. The hosts of such defense systems tend to have a smaller pangenome size and harbor fewer phage-related genes compared to genomes lacking these systems, suggesting that these defense mechanisms inhibit HGT by limiting the integration of prophages. We hypothesize that restriction of HGT by defense systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and fitness effect of HGT in bacterial populations.

CRISPR genotyping methods: Tracing the evolution from spoligotyping to machine learning

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems provide prokaryotes with adaptive immunity defenses against foreign genetic invaders. The identification of CRISPR-Cas function is among the most impactful discoveries of recent decades that have shaped the development of genome editing in various organisms paving the way for a plethora of promising applications in biotechnology and health. Even before the discovery of CRISPR-Cas biological role, the particular structure of CRISPR loci has been explored for epidemiological genotyping of bacterial pathogens. CRISPR-Cas loci are arranged in CRISPR arrays of mostly identical direct repeats intercalated with invader-derived spacers and an operon of cas genes encoding the Cas protein components. Each small CRISPR RNA (crRNA) encoded within the CRISPR array constitutes a key functional unit of this RNA-based CRISPR-Cas defense system guiding the Cas effector proteins toward the foreign nucleic acids for their destruction. The information acquired from prior invader encounters and stored within CRISPR arrays turns out to be extremely valuable in tracing the microevolution and epidemiology of major bacterial pathogens. We review here the history of CRISPR-based typing strategies highlighting the first PCR-based methods that have set the stage for recent developments of high-throughput sequencing and machine learning-based approaches. A great amount of whole genome sequencing and metagenomic data accumulated in recent years opens up new avenues for combining experimental and computational approaches of high-resolution CRISPR-based typing.

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.

Identification and characterization of two CRISPR/Cas systems associated with the mosquito microbiome

The microbiome profoundly influences many traits in medically relevant vectors such as mosquitoes, and a greater functional understanding of host-microbe interactions may be exploited for novel microbial-based approaches to control mosquito-borne disease. Here, we characterized two novel clustered regularly interspaced short palindromic repeats (CRISPR)/Cas systems in Serratia sp. Ag1, which was isolated from the gut of an Anopheles gambiae mosquito. Two distinct CRISPR/Cas systems were identified in Serratia Ag1, CRISPR1 and CRISPR2. Based on cas gene composition, CRISPR1 is classified as a type I-E CRISPR/Cas system and has a single array, CRISPR1. CRISPR2 is a type I-F system with two arrays, CRISPR2.1 and CRISPR2.2. RT-PCR analyses show that all cas genes from both systems are expressed during logarithmic growth in culture media. The direct repeat sequences of CRISPRs 2.1 and 2.2 are identical and found in the arrays of other Serratia spp., including S. marcescens and S. fonticola , whereas CRISPR1 is not. We searched for potential spacer targets and revealed an interesting difference between the two systems: only 9 % of CRISPR1 (type I-E) targets are in phage sequences and 91 % are in plasmid sequences. Conversely, ~66 % of CRISPR2 (type I-F) targets are found within phage genomes. Our results highlight the presence of CRISPR loci in gut-associated bacteria of mosquitoes and indicate interplay between symbionts and invasive mobile genetic elements over evolutionary time.

Genomic analysis of Anderson typing phages of Salmonella Typhimrium: towards understanding the basis of bacteria-phage interaction

The Anderson phage typing scheme has been successfully used worldwide for epidemiological surveillance of Salmonella enterica serovar Typhimurium. Although the scheme is being replaced by whole genome sequence subtyping methods, it can provide a valuable model system for study of phage-host interaction. The phage typing scheme distinguishes more than 300 definitive types of Salmonella Typhimurium based on their patterns of lysis to a unique collection of 30 specific Salmonella phages. In this study, we sequenced the genomes of 28 Anderson typing phages of Salmonella Typhimurium to begin to characterize the genetic determinants that are responsible for the differences in these phage type profiles. Genomic analysis of typing phages reveals that Anderson phages can be classified into three different groups, the P22-like, ES18-like and SETP3-like clusters. Most Anderson phages are short tailed P22-like viruses (genus Lederbergvirus); but phages STMP8 and STMP18 are very closely related to the lambdoid long tailed phage ES18, and phages STMP12 and STMP13 are related to the long noncontractile tailed, virulent phage SETP3. Most of these typing phages have complex genome relationships, but interestingly, two phage pairs STMP5 and STMP16 as well as STMP12 and STMP13 differ by a single nucleotide. The former affects a P22-like protein involved in DNA passage through the periplasm during its injection, and the latter affects a gene whose function is unknown. Using the Anderson phage typing scheme would provide insights into phage biology and the development of phage therapy for the treatment of antibiotic resistant bacterial infections.

Analysis of CRISPR/Cas Genetic Structure, Spacer Content and Molecular Epidemiology in Brazilian Acinetobacter baumannii Clinical Isolates

CRISPR/Cas is a molecular mechanism to prevent predatory viruses from invading bacteria via the insertion of small viral sequences (spacers) in its repetitive locus. The nature of spacer incorporation and the viral origins of spacers provide an overview of the genetic evolution of bacteria, their natural viral predators, and the mechanisms that prokaryotes may use to protect themselves, or to acquire mobile genetic elements such as plasmids. Here, we report on the CRISPR/Cas genetic structure, its spacer content, and strain epidemiology through MLST and CRISPR typing in Acinetobacter baumannii, an opportunistic pathogen intimately related to hospital infections and antimicrobial resistance. Results show distinct genetic characteristics, such as polymorphisms specific to ancestor direct repeats, a well-defined degenerate repeat, and a conserved leader sequence, as well as showing most spacers as targeting bacteriophages, and several self-targeting spacers, directed at prophages. There was a particular relationship between CRISPR/Cas and CC113 in the study of Brazilian isolates, and CRISPR-related typing techniques are interesting for subtyping strains with the same MLST profile. We want to emphasize the significance of descriptive genetic research on CRISPR loci, and we argue that spacer or CRISPR typing are helpful for small-scale investigations, preferably in conjunction with other molecular typing techniques such as MLST.

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.

Bioinformatic survey of CRISPR loci across 15 Serratia species

The Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR-associated proteins (CRISPR-Cas) system of prokaryotes is an adaptative immune defense mechanism to protect themselves from invading genetic elements (e.g., phages and plasmids). Studies that describe the genetic organization of these prokaryotic systems have mainly reported on the Enterobacteriaceae family (now reorganized within the order of Enterobacterales). For some genera, data on CRISPR-Cas systems remain poor, as in the case of Serratia (now part of the Yersiniaceae family) where data are limited to a few genomes of the species marcescens. This study describes the detection, in silico, of CRISPR loci in 146 Serratia complete genomes and 336 high-quality assemblies available for the species ficaria, fonticola, grimesii, inhibens, liquefaciens, marcescens, nematodiphila, odorifera, oryzae, plymuthica, proteomaculans, quinivorans, rubidaea, symbiotica, and ureilytica. Apart from subtypes I-E and I-F1 which had previously been identified in marcescens, we report that of I-C and the I-E unique locus 1, I-E*, and I-F1 unique locus 1. Analysis of the genomic contexts for CRISPR loci revealed mdtN-phnP as the region mostly shared (grimesii, inhibens, marcescens, nematodiphila, plymuthica, rubidaea, and Serratia sp.). Three new contexts detected in genomes of rubidaea and fonticola (puu genes-mnmA) and rubidaea (osmE-soxG and ampC-yebZ) were also found. The plasmid and/or phage origin of spacers was also established.

Correlation between CRISPR Loci Diversity in Three Enterobacterial Taxa

CRISPR-Cas is an adaptive immunity system of prokaryotes, composed of CRISPR arrays and the associated proteins. The successive addition of spacer sequences in the CRISPR array has made the system a valuable molecular marker, with multiple applications. Due to the high degree of polymorphism of the CRISPR loci, their comparison in bacteria from various sources may provide insights into the evolution and spread of the CRISPR-Cas systems. The aim of this study was to establish a correlation between the enterobacterial CRISPR loci, the sequence of direct repeats (DR), and the number of spacer units, along with the geographical origin and collection source. For this purpose, 3474 genomes containing CRISPR loci from the CRISPRCasdb of Salmonella enterica, Escherichia coli, and Klebsiella pneumoniae were analyzed, and the information regarding the isolates was recorded from the NCBI database. The most prevalent was the I-E CRISPR-Cas system in all three studied taxa. E. coli also presents the I-F type, but in a much lesser percentage. The systems found in K. pneumoniae can be classified into I-E and I-E*. The I-E and I-F systems have two CRISPR loci, while I-E* has only one locus upstream of the Cas cluster. PCR primers have been developed in this study for each CRISPR locus. Distinct clustering was not evident, but statistically significant relationships occurred between the different CRISPR loci and the number of spacer units. For each of the queried taxa, the number of spacers was significantly different (p < 0.01) by origin (Africa, Asia, Australia and Oceania, Europe, North America, and South America) but was not linked to the isolation source type (human, animal, plant, food, or laboratory strains).

Clustered regularly interspaced short palindromic repeats-Cas system: diversity and regulation in Enterobacteriaceae

Insights into the arms race between bacteria and invading mobile genetic elements have revealed the intricacies of the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system and the counter-defenses of bacteriophages. Incredible spacer diversity but significant spacer conservation among species/subspecies dictates the specificity of the CRISPR-Cas system. Researchers have exploited this feature to type/subtype the bacterial strains, devise targeted antimicrobials and regulate gene expression. This review focuses on the nuances of the CRISPR-Cas systems in Enterobacteriaceae that predominantly harbor type I-E and I-F CRISPR systems. We discuss the systems' regulation by the global regulators, H-NS, LeuO, LRP, cAMP receptor protein and other regulators in response to environmental stress. We further discuss the regulation of noncanonical functions like DNA repair pathways, biofilm formation, quorum sensing and virulence by the CRISPR-Cas system. The review comprehends multiple facets of the CRISPR-Cas system in Enterobacteriaceae including its diverse attributes, association with genetic features, regulation and gene regulatory mechanisms.

Prophage Diversity Across Salmonella and Verotoxin-Producing Escherichia coli in Agricultural Niches of British Columbia, Canada

Prophages have long been regarded as an important contributor to the evolution of Salmonella and Verotoxin-producing E. coli (VTEC), members of the Enterobacteriaceae that cause millions of cases of foodborne illness in North America. In S. Typhimurium, prophages provide many of the genes required for invasion; similarly, in VTEC, the Verotoxin-encoding genes are located in cryptic prophages. The ability of prophages to quickly acquire and lose genes have driven their rapid evolution, leading to highly diversified populations of phages that can infect distantly-related bacterial hosts. To defend against foreign genetic materials (i.e., phages), bacteria have evolved Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) immunity, consisting of variable spacer regions that match short nucleic acid sequences of invaders previously encountered. The number of spacer regions varies widely amongst Enterobacteriaceae, and there is currently no clear consensus if the accumulation of spacers is linked to genomic prophage abundance. Given the immense prophage diversity and contribution to bacterial host phenotypes, we analyzed the prophage sequences within 118 strains of Salmonella and VTEC, 117 of which are of agricultural origin. Overall, 130 unique prophage sequences were identified and they were found to be remarkably diverse with <50% nucleotide similarity, particularly with the Gifsy-1 group which was identified in several Salmonella serovars and interestingly, a strain of VTEC. Additionally, we identified a novel plasmid-like phage that carried antibiotic resistance and bacteriocin resistance genes. The strains analyzed carried at least six distinct spacers which did not possess homology to prophages identified in the same genome. In fact, only a fraction of all identified spacers (14%) possessed significant homology to known prophages. Regression models did not discern a correlation between spacer and prophage abundance in our strains, although the relatively high number of spacers in our strains (an average of 27 in Salmonella and 19 in VTEC) suggest that high rates of infection may occur in agricultural niches and be a contributing driver in bacterial evolution. Cumulatively, these results shed insight into prophage diversity of Salmonella and VTEC, which will have further implications when informing development of phage therapies against these foodborne pathogens.

Evaluation of the use of CRISPR loci for discrimination of Salmonella enterica subsp. enterica serovar Enteritidis strains recovered in Canada and comparison with other subtyping methods

Salmonella enterica subsp. enterica serovar Enteritidis remains one of the most important foodborne pathogens worldwide. To minimise its public health impact when outbreaks of the disease occur, timely investigation to identify and recall the contaminated food source is necessary. Central to this approach is the need for rapid and accurate identification of the bacterial subtype epidemiologically linked to the outbreak. While traditional methods of S. Enteritidis subtyping, such as pulsed field gel electrophoresis (PFGE) and phage typing (PT), have played an important role, the clonal nature of this organism has spurred efforts to improve subtyping resolution and timeliness through molecular based approaches. This study uses a cohort of 92 samples, recovered from a variety of sources, to compare these two traditional methods for S. Enteritidis subtyping with recently developed molecular techniques. These latter methods include the characterisation of two clustered regularly interspaced short palindromic repeats (CRISPR) loci, either in isolation or together with sequence analysis of virulence genes such as fimH. For comparison, another molecular technique developed in this laboratory involved the scoring of 60 informative single nucleotide polymorphisms (SNPs) distributed throughout the genome. Based on both the number of subtypes identified and Simpson's index of diversity, the CRISPR method was the least discriminatory and not significantly improved with the inclusion of fimH gene sequencing. While PT analysis identified the most subtypes, the SNP-PCR process generated the greatest index of diversity value. Combining methods consistently improved the number of subtypes identified, with the SNP/CRISPR typing scheme generating a level of diversity comparable with that of PT/PFGE. While these molecular methods, when combined, may have significant utility in real-world situations, this study suggests that CRISPR analysis alone lacks the discriminatory capability required to support investigations of foodborne disease outbreaks.

The CRISPR-Cas System Differentially Regulates Surface-Attached and Pellicle Biofilm in Salmonella enterica Serovar Typhimurium

The CRISPR-Cas mediated regulation of biofilm by Salmonella enterica serovar Typhimurium was investigated by deleting CRISPR-Cas components ΔcrisprI, ΔcrisprII, ΔΔcrisprI crisprII, and Δcas op. We determined that the system positively regulates surface biofilm while inhibiting pellicle biofilm formation. Results of real-time PCR suggest that the flagellar (fliC, flgK) and curli (csgA) genes were repressed in knockout strains, causing reduced surface biofilm. The mutants displayed altered pellicle biofilm architecture. They exhibited bacterial multilayers and a denser extracellular matrix with enhanced cellulose and less curli, ergo weaker pellicles than those of the wild type. The cellulose secretion was more in the knockout strains due to the upregulation of bcsC, which is necessary for cellulose export. We hypothesized that the secreted cellulose quickly integrates into the pellicle, leading to enhanced pellicular cellulose in the knockout strains. We determined that crp is upregulated in the knockout strains, thereby inhibiting the expression of csgD and, hence, also of csgA and bcsA. The conflicting upregulation of bcsC, the last gene of the bcsABZC operon, could be caused by independent regulation by the CRISPR-Cas system owing to a partial match between the CRISPR spacers and bcsC gene. The cAMP-regulated protein (CRP)-mediated regulation of the flagellar genes in the knockout strains was probably circumvented through the regulation of yddx governing the availability of the sigma factor σ²⁸ that further regulates class 3 flagellar genes (fliC, fljB, and flgK). Additionally, the variations in the lipopolysaccharide (LPS) profile and expression of LPS-related genes (rfaC, rfbG, and rfbI) in knockout strains could also contribute to the altered pellicle architecture. Collectively, we establish that the CRISPR-Cas system differentially regulates the formation of surface-attached and pellicle biofilm.

IMPORTANCE In addition to being implicated in bacterial immunity and genome editing, the CRISPR-Cas system has recently been demonstrated to regulate endogenous gene expression and biofilm formation. While the function of individual cas genes in controlling Salmonella biofilm has been explored, the regulatory role of CRISPR arrays in biofilm is less studied. Moreover, studies have focused on the effects of the CRISPR-Cas system on surface-associated biofilms, and comprehensive studies on the impact of the system on pellicle biofilm remain an unexplored niche. We demonstrate that the CRISPR array and cas genes modulate the expression of various biofilm genes in Salmonella, whereby surface and pellicle biofilm formation is distinctively regulated.

Clustered Regularly Interspaced Short Palindromic Repeats Genotyping of Multidrug-Resistant Salmonella Heidelberg Strains Isolated From the Poultry Production Chain Across Brazil

Salmonella enterica subsp. enterica serovar Heidelberg has been associated with a broad host range, such as poultry, dairy calves, swine, wild birds, environment, and humans. The continuous evolution of S. Heidelberg raises a public health concern since there is a global dispersal of lineages harboring a wide resistome and virulome on a global scale. Here, we characterized the resistome, phylogenetic structure and clustered regularly interspaced short palindromic repeats (CRISPR) array composition of 81 S. Heidelberg strains isolated from broiler farms (n = 16), transport and lairage (n = 5), slaughterhouse (n = 22), and retail market (n = 38) of the poultry production chain in Brazil, between 2015 and 2016 using high-resolution approaches including whole-genome sequencing (WGS) and WGS-derived CRISPR genotyping. More than 91% of the S. Heidelberg strains were multidrug-resistant. The total antimicrobial resistance (AMR) gene abundances did not vary significantly across regions and sources suggesting the widespread distribution of antibiotic-resistant strains from farm to market. The highest AMR gene abundance was observed for fosA7, aac(6')-Iaa, sul2, tet(A), gyrA, and parC for 100% of the isolates, followed by 88.8% for bla CMY-2. The β-lactam resistance was essentially driven by the presence of the plasmid-mediated AmpC (pAmpC) bla CMY-2 gene, given the isolates which did not carry this gene were susceptible to cefoxitin (FOX). Most S. Heidelberg strains were classified within international lineages, which were phylogenetically nested with Salmonella strains from European countries; while CRISPR genotyping analysis revealed that the spacer content was overall highly conserved, but distributed into 13 distinct groups. In summary, our findings underscore the potential role of S. Heidelberg as a key pathogen disseminated from farm to fork in Brazil and reinforce the importance of CRISPR-based genotyping for salmonellae. Hence, we emphasized the need for continuous mitigation programs to monitor the dissemination of this high-priority pathogen.

Theragnostic application of nanoparticle and CRISPR against food-borne multi-drug resistant pathogens

Foodborne infection is one of the leading sources of infections spreading across the world. Foodborne pathogens are recognized as multidrug-resistant (MDR) pathogens posing a significant problem in the food industry and healthy consumers resulting in enhanced economic burden, and nosocomial infections. The continued search for enhanced microbial detection tools has piqued the interest of the CRISPR-Cas system and Nanoparticles. CRISPR-Cas system is present in the bacterial genome of some prokaryotes and is repurposed as a theragnostic tool against MDR pathogens. Nanoparticles and composites have also emerged as an efficient tool in theragnostic applications against MDR pathogens. The diagnostic limitations of the CRISPR-Cas system are believed to be overcome by a synergistic combination of the nanoparticles system and CRISPR-Cas using nanoparticles as vehicles. In this review, we have discussed the diagnostic application of CRISPR-Cas technologies along with their potential usage in applications like phage resistance, phage vaccination, strain typing, genome editing, and antimicrobial. we have also elucidated the antimicrobial and detection role of nanoparticles against foodborne MDR pathogens. Moreover, the novel combinatorial approach of CRISPR-Cas and nanoparticles for their synergistic effects in pathogen clearance and drug delivery vehicles has also been discussed.

Assessing Salmonella Prevalence and Complexity Through Processing Using Different Culture Methods

Conventional Salmonella surveillance requires a week for isolation, confirmation, and subsequent serotyping. We previously showed that this could be reduced by 24 h by combining the pre-enrichment and enrichment steps into a single selective pre-enrichment step and was tested on directly after picking. The goal of this study was 2-fold: 1) to evaluate the use of selective pre-enrichment through each step of processing, including postintervention when the Salmonella load is reduced, and 2) to assess any changes in serovar populations in Salmonella positive samples. Duplicate carcass drip samples, each representative of 500 broiler carcasses, were collected by catching processing water drip under moving carcass shackle lines in each of three commercial broiler slaughter plants. Samples were collected post-pick, post-inside-outside bird wash (IOBW), and post-chill; duplicate wing rinses were performed pre- and post-antimicrobial parts dip. Each processing plant was sampled 6 times for a total of 180 samples collected. The number of Salmonella positives identified with selective pre-enrichment conditions (48/180) was similar to traditional selective enrichment culture conditions (52/180), showed good concordance in recovery rate between the 2 culture methods (Fisher's exact test, P = 0.72). We also found that the incidence of Salmonella reduced dramatically after antimicrobial intervention (post-pick 66.7% vs. post chill 8.3%). When serovar populations were evaluated in Salmonella positive samples using CRISPR-SeroSeq, we detected four different Salmonella serovars, Kentucky, Infantis, Schwarzengrund, and Typhimurium, and their incidence rose between post-pick and post-IOBW. The relative abundance of Infantis within individual samples increased between post-pick and post-IOBW while the relative abundance of the other 3 serovars decreased. These results suggest that a selective pre-enrichment step reduces the time required for Salmonella isolation without negatively affecting detection and serovar profiles in culture positive samples were not altered between culture conditions used.

Regional Salmonella Differences in United States Broiler Production from 2016 to 2020 and the Contribution of Multiserovar Populations to Salmonella Surveillance

Poultry remains a considerable source of foodborne salmonellosis despite significant reduction of Salmonella incidence during processing. There are multiple entry points for Salmonella during production that can lead to contamination during slaughter, and it is important to distinguish the serovars present between the different stages to enact appropriate controls. National Salmonella data from the U.S. Department of Agriculture-Food Safety Inspection Service (USDA-FSIS) monitoring of poultry processing was analyzed from 2016 to 2020. The overall Salmonella incidence at processing in broiler carcasses and intact parts (parts) decreased from 9.00 to 6.57% over this period. The incidence in parts was higher (11.15%) than in carcasses (4.78%). Regional differences include higher proportions of serovars Infantis and Typhimurium in the Atlantic and higher proportion of serovar Schwarzengrund in the Southeast. For Georgia, the largest broiler-producing state, USDA-FSIS data were compared to Salmonella monitoring data from breeder flocks over the same period, revealing serovar Kentucky as the major serovar in breeders (67.91%) during production but not at processing, suggesting that it is more effectively removed during antimicrobial interventions. CRISPR-SeroSeq was performed on breeder samples collected between 2020 and 2021 to explain the incongruence between pre- and postharvest and showed that 32% of samples contain multiple serovars, with up to 11 serovars found in a single flock. High-resolution sequencing identifies serovar patterns at the population level and can provide insight to develop targeted controls. The work presented may apply to other food production systems where Salmonella is a concern, since it overcomes limitations associated with conventional culture. IMPORTANCE Salmonella is a leading cause of bacterial foodborne illness in the United States, with poultry as a significant Salmonella reservoir. We show the relative decrease in Salmonella over a 5-year period from 2016 to 2020 in processed chicken parts and highlight regional differences with respect to the prevalence of clinically important Salmonella serovars. Our results show that the discrepancy between Salmonella serovars found in pre- and postharvest poultry during surveillance are due in part by the limited detection depth offered by traditional culture techniques. Despite the reduction of Salmonella at processing, the number of human salmonellosis cases has remained stable, which may be attributed to differences in virulence among serovars and their associated risk. When monitoring for Salmonella, it is imperative to identify all serovars present to appropriately assess public health risk and to implement the most effective Salmonella controls.

CRISPR-Cas System: An Adaptive Immune System’s Association with Antibiotic Resistance in Salmonella enterica Serovar Enteritidis

Several factors are involved in the emergence of antibiotic-resistant bacteria and pose a serious threat to public health safety. Among them, clustered regularly interspaced short palindromic repeat- (CRISPR-) Cas system, an adaptive immune system, is thought to be involved in the development of antibiotic resistance in bacteria. The current study was aimed at determining not only the presence of antibiotic resistance and CRISPR-Cas system but also their association with each other in Salmonella enteritidis isolated from the commercial poultry. A total of 139 samples were collected from poultry birds sold at the live bird markets of Lahore City, and both phenotypic and genotypic methods were used to determine antimicrobial resistance. The presence of the CRISPR-Cas system was determined by PCR, followed by sequencing. All isolates of S. enteritidis (100%) were resistant to nalidixic acid, whereas 95% of isolates were resistant to ampicillin. Five multidrug-resistant isolates (MDR) such as S. enteritidis isolate (S. E1, S. E2, S. E4, S. E5, and S. E8) were found in the present study. The CRISPR-Cas system was detected in all of these MDR isolates, and eight spacers were detected within the CRISPR array. In addition, an increased expression of CRISPR-related genes was observed in the standard strain and MDR S. enteritidis isolates. The association of the CRISPSR-Cas system with multiple drug resistance highlights the exogenous acquisition of genes by horizontal transfer. The information could be used further to combat antibiotic resistance in pathogens like Salmonella.

The genetic basis of phage susceptibility, cross-resistance and host-range in Salmonella

Though bacteriophages (phages) are known to play a crucial role in bacterial fitness and virulence, our knowledge about the genetic basis of their interaction, cross-resistance and host-range is sparse. Here, we employed genome-wide screens in Salmonella enterica serovar Typhimurium to discover host determinants involved in resistance to eleven diverse lytic phages including four new phages isolated from a therapeutic phage cocktail. We uncovered 301 diverse host factors essential in phage infection, many of which are shared between multiple phages demonstrating potential cross-resistance mechanisms. We validate many of these novel findings and uncover the intricate interplay between RpoS, the virulence-associated general stress response sigma factor and RpoN, the nitrogen starvation sigma factor in phage cross-resistance. Finally, the infectivity pattern of eleven phages across a panel of 23 genome sequenced Salmonella strains indicates that additional constraints and interactions beyond the host factors uncovered here define the phage host range.

Comparative Genomic Analyses and CRISPR-Cas Characterization of Cutibacterium acnes Provide Insights Into Genetic Diversity and Typing Applications

Cutibacterium acnes is an important member of the human skin microbiome and plays a critical role in skin health and disease. C. acnes encompasses different phylotypes that have been found to be associated with different skin phenotypes, suggesting a genetic basis for their impact on skin health. Here, we present a comprehensive comparative analysis of 255 C. acnes genomes to provide insights into the species genetic diversity and identify unique features that define various phylotypes. Results revealed a relatively small and open pan genome (6,240 genes) with a large core genome (1,194 genes), and three distinct phylogenetic clades, with multiple robust sub-clades. Furthermore, we identified several unique gene families driving differences between distinct C. acnes clades. Carbohydrate transporters, stress response mechanisms and potential virulence factors, potentially involved in competitive growth and host colonization, were detected in type I strains, which are presumably responsible for acne. Diverse type I-E CRISPR-Cas systems and prophage sequences were detected in select clades, providing insights into strain divergence and adaptive differentiation. Collectively, these results enable to elucidate the fundamental differences among C. acnes phylotypes, characterize genetic elements that potentially contribute to type I-associated dominance and disease, and other key factors that drive the differentiation among clades and sub-clades. These results enable the use of comparative genomics analyses as a robust method to differentiate among the C. acnes genotypes present in the skin microbiome, opening new avenues for the development of biotherapeutics to manipulate the skin microbiota.

PHIDA: A High Throughput Turbidimetric Data Analytic Tool to Compare Host Range Profiles of Bacteriophages Isolated Using Different Enrichment Methods

Bacteriophages are viruses that infect bacteria and are present in niches where bacteria thrive. In recent years, the suggested application areas of lytic bacteriophage have been expanded to include therapy, biocontrol, detection, sanitation, and remediation. However, phage application is constrained by the phage's host range-the range of bacterial hosts sensitive to the phage and the degree of infection. Even though phage isolation and enrichment techniques are straightforward protocols, the correlation between the enrichment technique and host range profile has not been evaluated. Agar-based methods such as spotting assay and efficiency of plaquing (EOP) are the most used methods to determine the phage host range. These methods, aside from being labor intensive, can lead to subjective and incomplete results as they rely on qualitative observations of the lysis/plaques, do not reflect the lytic activity in liquid culture, and can overestimate the host range. In this study, phages against three bacterial genera were isolated using three different enrichment methods. Host range profiles of the isolated phages were quantitatively determined using a high throughput turbidimetric protocol and the data were analyzed with an accessible analytic tool "PHIDA". Using this tool, the host ranges of 9 Listeria, 14 Salmonella, and 20 Pseudomonas phages isolated with different enrichment methods were quantitatively compared. A high variability in the host range index (HRi) ranging from 0.86-0.63, 0.07-0.24, and 0.00-0.67 for Listeria, Salmonella, and Pseudomonas phages, respectively, was observed. Overall, no direct correlation was found between the phage host range breadth and the enrichment method in any of the three target bacterial genera. The high throughput method and analytics tool developed in this study can be easily adapted to any phage study and can provide a consensus for phage host range determination.

Shining light on a deep-sea bacterial symbiont population structure with CRISPR

Many foundation species in chemosynthesis-based ecosystems rely on environmentally acquired symbiotic bacteria for their survival. Hence, understanding the biogeographic distributions of these symbionts at regional scales is key to understanding patterns of connectivity and predicting resilience of their host populations (and thus whole communities). However, such assessments are challenging because they necessitate measuring bacterial genetic diversity at fine resolutions. For this purpose, the recently discovered clustered regularly interspaced short palindromic repeats (CRISPR) constitutes a promising new genetic marker. These DNA sequences harboured by about half of bacteria hold their viral immune memory, and as such, might allow discrimination of different lineages or strains of otherwise indistinguishable bacteria. In this study, we assessed the potential of CRISPR as a hypervariable phylogenetic marker in the context of a population genetic study of an uncultured bacterial species. We used high-throughput CRISPR-based typing along with multi-locus sequence analysis (MLSA) to characterize the regional population structure of the obligate but environmentally acquired symbiont species Candidatus Endoriftia persephone on the Juan de Fuca Ridge. Mixed symbiont populations of Ca. Endoriftia persephone were sampled across individual Ridgeia piscesae hosts from contrasting habitats in order to determine if environmental conditions rather than barriers to connectivity are more important drivers of symbiont diversity. We showed that CRISPR revealed a much higher symbiont genetic diversity than the other housekeeping genes. Several lines of evidence imply this diversity is indicative of environmental strains. Finally, we found with both CRISPR and gene markers that local symbiont populations are strongly differentiated across sites known to be isolated by deep-sea circulation patterns. This research showed the high power of CRISPR to resolve the genetic structure of uncultured bacterial populations and represents a step towards making keystone microbial species an integral part of conservation policies for upcoming mining operations on the seafloor.

Antimicrobial Resistance Hidden within Multiserovar Salmonella Populations

Salmonella enterica can exist in food animals as multiserovar populations, and different serovars can harbor diverse antimicrobial resistance (AMR) profiles. Conventional Salmonella isolation assesses AMR only in the most abundant members of a multiserovar population, which typically reflects their relative abundance in the initial sample. Therefore, AMR in underlying serovars is an undetected reservoir that can readily be expanded upon antimicrobial use. CRISPR-SeroSeq profiling demonstrated that 60% of cattle fecal samples harbored multiple serovars, including low levels of Salmonella serovar Reading in 11% of samples, which were not found by culture-based Salmonella isolation. An in vitro challenge revealed that Salmonella serovar Reading was tetracycline resistant, while more abundant serovars were susceptible. This study highlights the importance of AMR surveillance in multiserovar populations.

Characterisation of Phage Susceptibility Variation in Salmonellaenterica Serovar Typhimurium DT104 and DT104b

The surge in mortality and morbidity rates caused by multidrug-resistant (MDR) bacteria prompted a renewal of interest in bacteriophages (phages) as clinical therapeutics and natural biocontrol agents. Nevertheless, bacteria and phages are continually under the pressure of the evolutionary phage-host arms race for survival, which is mediated by co-evolving resistance mechanisms. In Anderson phage typing scheme of Salmonella Typhimurium, the epidemiologically related definitive phage types, DT104 and DT104b, display significantly different phage susceptibility profiles. This study aimed to characterise phage resistance mechanisms and genomic differences that may be responsible for the divergent phage reaction patterns in S. Typhimurium DT104 and DT104b using whole genome sequencing (WGS). The analysis of intact prophages, restriction-modification systems (RMS), plasmids and clustered regularly interspaced short palindromic repeats (CRISPRs), as well as CRISPR-associated proteins, revealed no unique genetic determinants that might explain the variation in phage susceptibility among the two phage types. Moreover, analysis of genes coding for potential phage receptors revealed no differences among DT104 and DT104b strains. However, the findings propose the need for experimental assessment of phage-specific receptors on the bacterial cell surface and analysis of bacterial transcriptome using RNA sequencing which will explain the differences in bacterial susceptibility to phages. Using Anderson phage typing scheme of Salmonella Typhimurium for the study of bacteria-phage interaction will help improving our understanding of host-phage interactions which will ultimately lead to the development of phage-based technologies, enabling effective infection control.

The CRISPR-Cas System Is Involved in OmpR Genetic Regulation for Outer Membrane Protein Synthesis in Salmonella Typhi

The CRISPR-Cas cluster is found in many prokaryotic genomes including those of the Enterobacteriaceae family. Salmonella enterica serovar Typhi (S. Typhi) harbors a Type I-E CRISPR-Cas locus composed of cas3, cse1, cse2, cas7, cas5, cas6e, cas1, cas2, and a CRISPR1 array. In this work, it was determined that, in the absence of cas5 or cas2, the amount of the OmpC porin decreased substantially, whereas in individual cse2, cas6e, cas1, or cas3 null mutants, the OmpF porin was not observed in an electrophoretic profile of outer membrane proteins. Furthermore, the LysR-type transcriptional regulator LeuO was unable to positively regulate the expression of the quiescent OmpS2 porin, in individual S. Typhi cse2, cas5, cas6e, cas1, cas2, and cas3 mutants. Remarkably, the expression of the master porin regulator OmpR was dependent on the Cse2, Cas5, Cas6e, Cas1, Cas2, and Cas3 proteins. Therefore, the data suggest that the CRISPR-Cas system acts hierarchically on OmpR to control the synthesis of outer membrane proteins in S. Typhi.

Exploiting heterologous and endogenous CRISPR-Cas systems for genome editing in the probiotic Clostridium butyricum

Clostridium butyricum has been widely used as a probiotic for humans and food animals. However, the mechanisms of beneficial effects of C. butyricum on the host remain poorly understood, largely due to the lack of high-throughput genome engineering tools. Here, we report the exploitation of heterologous Type II CRISPR-Cas9 system and endogenous Type I-B CRISPR-Cas system in probiotic C. butyricum for seamless genome engineering. Although successful genome editing was achieved in C. butyricum when CRISPR-Cas9 system was employed, the expression of toxic cas9 gene result in really poor transformation, spurring us to develop an easy-applicable and high-efficient genome editing tool. Therefore, the endogenous Type I-B CRISPR-Cas machinery located on the megaplasmid of C. butyricum was co-opted for genome editing. In vivo plasmid interference assays identified that ACA and TAA were functional protospacer adjacent motif sequences needed for site-specific CRISPR attacking. Using the customized endogenous CRISPR-Cas system, we successfully deleted spo0A and aldh genes in C. butyricum, yielding an efficiency of up to 100%. Moreover, the conjugation efficiency of endogenous CRISPR-Cas system was dramatically enhanced due to the precluding expression of cas9. Altogether, the two approaches developed herein remarkably expand the existing genetic toolbox available for investigation of C. butyricum.

Survival of Salmonella Under Heat Stress is Associated with the Presence/Absence of CRISPR Cas Genes and Iron Levels

Clustered regularly interspaced short palindromic repeats (CRISPR) cas genes have been linked to stress response in Salmonella. Our aim was to identify the presence of CRISPR cas in Salmonella and its response to heat in the presence of iron. Whole genomes of Salmonella (n = 50) of seven serovars were compared to identify the presence of CRISPR cas genes, direct-repeats and spacers. All Salmonella genomes had all cas genes present except S. Newport 2393 which lacked these genes. Gene-specific primers were used to confirm the absence of these genes in S. Newport 2393. The presence/absence of CRISPR cas genes was further investigated among 469 S. Newport genomes from PATRIC with 283 genomes selected for pan-genome analysis. The response of eleven Salmonella strains of various serovars to gradual heat in ferrous and ferric forms of iron was investigated. A total of 32/283 S. Newport genomes that lacked all CRISPR cas genes clustered together. S. Newport 2393 was the most heat-sensitive strain at higher iron levels (200 and 220 pm) in ferrous and ferric forms of iron. The absence of CRISPR cas genes in S. Newport 2393 may contribute to its increase in heat sensitivity and iron may play a role in this. The high reduction in numbers of most Salmonella strains exposed to heat makes it unfeasible to extract RNA and conduct transcription studies. Further studies should be conducted to validate the survival of Salmonella when exposed to heat in the presence/absence of CRISPR cas genes and different iron levels.

Salmonella Genomics and Population Analyses Reveal High Inter- and Intraserovar Diversity in Freshwater

Freshwater can support the survival of the enteric pathogen Salmonella, though temporal Salmonella diversity in a large watershed has not been assessed. At 28 locations within the Susquehanna River basin, 10-liter samples were assessed in spring and summer over 2 years. Salmonella prevalence was 49%, and increased river discharge was the main driver of Salmonella presence. The amplicon-based sequencing tool, CRISPR-SeroSeq, was used to determine serovar population diversity and detected 25 different Salmonella serovars, including up to 10 serovars from a single water sample. On average, there were three serovars per sample, and 80% of Salmonella-positive samples contained more than one serovar. Serovars Give, Typhimurium, Thompson, and Infantis were identified throughout the watershed and over multiple collections. Seasonal differences were evident: serovar Give was abundant in the spring, whereas serovar Infantis was more frequently identified in the summer. Eight of the ten serovars most commonly associated with human illness were detected in this study. Crucially, six of these serovars often existed in the background, where they were masked by a more abundant serovar(s) in a sample. Serovars Enteritidis and Typhimurium, especially, were masked in 71 and 78% of samples where they were detected, respectively. Whole-genome sequencing-based phylogeny demonstrated that strains within the same serovar collected throughout the watershed were also very diverse. The Susquehanna River basin is the largest system where Salmonella prevalence and serovar diversity have been temporally and spatially investigated, and this study reveals an extraordinary level of inter- and intraserovar diversity.IMPORTANCESalmonella is a leading cause of bacterial foodborne illness in the United States, and outbreaks linked to fresh produce are increasing. Understanding Salmonella ecology in freshwater is of importance, especially where irrigation practices or recreational use occur. As the third largest river in the United States east of the Mississippi, the Susquehanna River is the largest freshwater contributor to the Chesapeake Bay, and it is the largest river system where Salmonella diversity has been studied. Rainfall and subsequent high river discharge rates were the greatest indicators of Salmonella presence in the Susquehanna and its tributaries. Several Salmonella serovars were identified, including eight commonly associated with foodborne illness. Many clinically important serovars were present at a low frequency within individual samples and so could not be detected by conventional culture methods. The technologies employed here reveal an average of three serovars in a 10-liter sample of water and up to 10 serovars in a single sample.

Genomic signatures of adaptation to natural settings in non-typhoidal Salmonella enterica Serovars Saintpaul, Thompson and Weltevreden

Salmonella enterica is a pathogenic bacterium responsible for intestinal illness and systemic diseases such as typhoid and paratyphoid fevers. Among clinical manifestation classification, non-typhoidal Salmonella is mainly known as foodborne pathogen associated with the consumption of fecal contaminated food and water. Even though Salmonella hosts include humans and warm-blooded animals, it has been found in non-host environments as river water where the bacteria use different strategies to fitness the environment persisting and establishment. Now with the availability of WGS and bioinformatics tools, we can explore bacterial genomes with higher resolution to increase our understanding of specific genetic signatures among environmental and clinical isolates, being the goal of this work. Pangenome construction allowed the detection of specific environmental and clinical gene clusters related to metabolism and secretion systems as the main signature respectively. Specifically, D-galactonate degradation pathway was observed mainly in environmental genomes while T3SS and flagellum genes were detected for all clinical but not for all environmental isolates. Gene duplication and pseudogenes accumulation were detected as the main adaptation strategy for environmental isolates; thus, isolation source may play an important role in genome plasticity, conferring a selective advantage to survive and persist for environmental Salmonella isolates. Intact prophage sequences with cargo genes were observable for both isolation sources playing an important role in virulence contribution.

CRISPR Typing of Salmonella Isolates

Polymerase chain reaction (PCR) and sequencing-based subtyping tools are useful for rapid analyses of Salmonella isolates. Here we describe the process of clustered regularly interspaced short palindromic repeat-multiple virulence locus sequence typing (CRISPR-MVLST) for Salmonella subtyping.

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.

Incidence of Salmonella serovars isolated from commercial animal feed mills in the United States and serovar identification using CRISPR analysis

AIMS

In this study, we sought to determine the incidence and diversity of Salmonella in a broad collection of commercial animal feeds collected from animal feed mills across the United States over an 11-month period and utilize CRISPR analysis to identify individual serovars.

METHODS AND RESULTS

Over two independent trials, 387 feed samples from 135 different animal feed mills in the United States were screened for Salmonella. A total of 6·2% (24/387) of samples were contaminated with Salmonella, which is concordant with similar studies. Clustered regularly interspaced short palindromic repeats (CRISPR)-typing was used to serotype Salmonella isolates, and serovars Infantis and Tennessee were the most common.

CONCLUSIONS

Serogroups O:4 and O:7 were enriched in the feed samples, suggesting that these serogroups are better adapted to surviving in low moisture animal feeds. The study supports the utility of CRISPR to determine serovar type since most of the serovars identified in this study have been also isolated and identified in earlier studies using more classical serotyping methods.

SIGNIFICANCE AND IMPACT OF THE STUDY

This work contributes to a growing body of literature concerning the Salmonella prevalence in animal feeds and highlights the need to effectively mitigate pathogens in livestock and poultry feed.

CRISPR-Cas Diversity in Clinical Salmonella enterica Serovar Typhi Isolates from South Asian Countries

Typhoid fever, caused by Salmonella enterica serovar Typhi (S. Typhi), is a global health concern and its treatment is problematic due to the rise in antimicrobial resistance (AMR). Rapid detection of patients infected with AMR positive S. Typhi is, therefore, crucial to prevent further spreading. Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR-associated genes (CRISPR-Cas), is an adaptive immune system that initially was used for typing purposes. Later, it was discovered to play a role in defense against phages and plasmids, including ones that carry AMR genes, and, at present, it is being explored for its usage in diagnostics. Despite the availability of whole-genome sequences (WGS), very few studied the CRISPR-Cas system of S. Typhi, let alone in typing purposes or relation to AMR. In the present study, we analyzed the CRISPR-Cas system of S. Typhi using WGS data of 1059 isolates obtained from Bangladesh, India, Nepal, and Pakistan in combination with demographic data and AMR status. Our results reveal that the S. Typhi CRISPR loci can be classified into two groups: A (evidence level >2) and B (evidence level ≤2), in which we identified a total of 47 unique spacers and 15 unique direct repeats. Further analysis of the identified spacers and repeats demonstrated specific patterns that harbored significant associations with genotype, demographic characteristics, and AMR status, thus raising the possibility of their usage as biomarkers. Potential spacer targets were identified and, interestingly, the phage-targeting spacers belonged to the group-A and plasmid-targeting spacers to the group-B CRISPR loci. Further analyses of the spacer targets led to the identification of an S. Typhi protospacer adjacent motif (PAM) sequence, TTTCA/T. New cas-genes known as DinG, DEDDh, and WYL were also discovered in the S. Typhi genome. However, a specific variant of the WYL gene was only identified in the extensively drug-resistant (XDR) lineage from Pakistan and ciprofloxacin-resistant lineage from Bangladesh. From this work, we conclude that there are strong correlations between variations identified in the S. Typhi CRISPR-Cas system and endemic AMR positive S. Typhi isolates.

Transcription termination and antitermination of bacterial CRISPR arrays

A hallmark of CRISPR-Cas immunity systems is the CRISPR array, a genomic locus consisting of short, repeated sequences ('repeats') interspersed with short, variable sequences ('spacers'). CRISPR arrays are transcribed and processed into individual CRISPR RNAs that each include a single spacer, and direct Cas proteins to complementary sequences in invading nucleic acid. Most bacterial CRISPR array transcripts are unusually long for untranslated RNA, suggesting the existence of mechanisms to prevent premature transcription termination by Rho, a conserved bacterial transcription termination factor that rapidly terminates untranslated RNA. We show that Rho can prematurely terminate transcription of bacterial CRISPR arrays, and we identify a widespread antitermination mechanism that antagonizes Rho to facilitate complete transcription of CRISPR arrays. Thus, our data highlight the importance of transcription termination and antitermination in the evolution of bacterial CRISPR-Cas systems.

Molecular Epidemiology of Acinetobacter calcoaceticus-Acinetobacter baumannii Complex Isolated From Children at the Hospital Infantil de México Federico Gómez

The Acinetobacter calcoaceticus-baumannii (Acb) complex is regarded as a group of phenotypically indistinguishable opportunistic pathogens responsible for mainly causing hospital-acquired pneumonia and bacteremia. The aim of this study was to determine the frequency of isolation of the species that constitute the Acb complex, as well as their susceptibility to antibiotics, and their distribution at the Hospital Infantil de Mexico Federico Gomez (HIMFG). A total of 88 strains previously identified by Vitek 2®, 40 as Acinetobacter baumannii and 48 as Acb complex were isolated from 52 children from 07, January 2015 to 28, September 2017. A. baumannii accounted for 89.77% (79/88) of the strains; Acinetobacter pittii, 6.82% (6/88); and Acinetobacter nosocomialis, 3.40% (3/88). Most strains were recovered mainly from patients in the intensive care unit (ICU) and emergency wards. Blood cultures (BC) provided 44.32% (39/88) of strains. The 13.63% (12/88) of strains were associated with primary bacteremia, 3.4% (3/88) with secondary bacteremia, and 2.3% (2/88) with pneumonia. In addition, 44.32% (39/88) were multidrug-resistant (MDR) strains and, 11.36% (10/88) were extensively drug-resistant (XDR). All strains amplified the bla_OXA-51 gene; 51.13% (45/88), the bla_OXA-23 gene; 4.54% (4/88), the bla_OXA-24 gene; and 2.27% (2/88), the bla_OXA-58 gene. Plasmid profiles showed that the strains had 1–6 plasmids. The strains were distributed in 52 pulsotypes, and 24 showed identical restriction patterns, with a correlation coefficient of 1.0. Notably, some strains with the same pulsotype were isolated from different patients, wards, or years, suggesting the persistence of more than one clone. Twenty-seven sequence types (STs) were determined for the strains based on a Pasteur multilocus sequence typing (MLST) scheme using massive sequencing; the most prevalent was ST 156 (27.27%, 24/88). The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas I-Fb system provided amplification in A. baumannii and A. pittii strains (22.73%, 20/88). This study identified an increased number of MDR strains and the relationship among strains through molecular typing. The data suggest that more than one strain could be causing an infection in some patient. The implementation of molecular epidemiology allowed the characterization of a set of strains and identification of different attributes associated with its distribution in a specific environment.

Whole-genome sequencing analysis and CRISPR genotyping of rare antibiotic-resistant Salmonella enterica serovars isolated from food and related sources

For decades, Salmonella Typhimurium and Salmonella Enteritidis have prevailed in several countries as agents of salmonellosis outbreaks. In Brazil, the largest exporter of poultry meat, relatively little attention has been paid to infrequent serovars. Here, we report the emergence and characterization of rare serovars isolated from food and related sources collected between 2014 and 2016 in Brazil. Twenty-two Salmonella enterica isolates were analyzed through the use of whole-genome sequencing (WGS) and clustered regularly interspaced short palindromic repeats (CRISPR) genotyping. These isolates were classified into 10 infrequent serovars, including S. Abony, S. Isangi, S. Rochdale, S. Saphra, S. Orion, S. Ouakam, S. Grumpensis, S. Carrau, S. Abaetetuba, and S. Idikan. The presence of six antimicrobial resistance (AMR) genes, qnrB19, bla_CMY-2, tetA, aac(6')-Iaa, sul2 and fosA7, which encode resistance to quinolones, third-generation cephalosporin, tetracycline, aminoglycoside, sulfonamide and fosfomycin, respectively, were confirmed by WGS. All S. Isangi harbored qnrB19 with conserved genomic context across strains, while S. Abony harbored bla_CMY-2. Twelve (54.5%) strains displayed chromosomal mutations in parC (Thr57→Ser). Most serovars were classified as independent lineages, except S. Abony and S. Abaetetuba, which phylogenetically nested with Salmonella strains from different countries. CRISPR analysis revealed that the spacer content was strongly correlated with serovar and multi-locus sequence type for all strains, independently confirming the observed phylogenetic patterns, and highlighting the value of CRISPR-based genotyping for Salmonella. These findings add valuable information to the epidemiology of S. enterica in Brazil, where the emergency of antibiotic-resistant Salmonella continues to evolve.

Conserved CRISPR arrays in Salmonella enterica serovar Infantis can serve as qPCR targets to detect Infantis in mixed serovar populations

Salmonellosis is a leading bacterial cause of foodborne illness, and numerous Salmonella enterica serovars have been responsible for foodborne outbreaks. In the United States outbreaks are often linked to poultry and poultry-related products. The prevalence of Salmonella serovar Infantis has been increasing in poultry processing facilities over the past few years and in 2018 was identified as the causative agent for a large multistate outbreak linked to raw chicken. CRISPR-typing is a subtyping approach based on PCR and the sequencing of two Salmonella loci, CRISPR1 and CRISPR2. CRISPR-typing was used to interrogate 138 recent (2018-2019) isolates and genomes of ser. Infantis. Results show that the CRISPR elements are remarkably conserved in this serovar. The most conserved spacers, and those also unique to ser. Infantis, were used as targets to develop a ser. Infantis-specific qPCR assay. This assay was able to detect ser. Infantis in mixed serovar cultures of Salmonella, down to 0·1% of the population, highlighting the utility of this molecular approach in improving surveillance sensitivity for this important food safety pathogen. SIGNIFICANCE AND IMPACT OF THE STUDY: The incidence of human salmonellosis cases caused by Salmonella enterica serovar Infantis (ser. Infantis) has been increasing, as has its prevalence in broiler chickens, which are a frequent reservoir of Salmonella. A cluster of ser. Infantis genetically linked to an outbreak strain have been identified in numerous processing facilities. A qPCR assay targeting CRISPR elements that are unique to ser. Infantis has been developed and can detect this serovar directly from mixed cultures. This assay is sensitive enough to reveal ser. Infantis within a mixed Salmonella population where it constitutes only 0·1% of the population. The rapid nature of qPCR lends this assay to high-throughput screening of poultry samples to detect this important pathogen.

Source Tracking Based on Core Genome SNV and CRISPR Typing of Salmonella enterica Serovar Heidelberg Isolates Involved in Foodborne Outbreaks in Québec, 2012

Whole-genome sequencing (WGS) is the method of choice for bacterial subtyping and it is rapidly replacing the more traditional methods such as pulsed-field gel electrophoresis (PFGE). Here we used the high-resolution core genome single nucleotide variant (cgSNV) typing method to characterize clinical and food from Salmonella enterica serovar Heidelberg isolates in the context of source attribution. Additionally, clustered regularly interspaced short palindromic repeats (CRISPR) analysis was included to further support this method. Our results revealed that cgSNV was highly discriminatory and separated the outbreak isolates into distinct clusters (0-4 SNVs). CRISPR analysis was also able to distinguish outbreak strains from epidemiologically unrelated isolates. Specifically, our data clearly demonstrated the strength of these two methods to determine the probable source(s) of a 2012 epidemiologically characterized outbreak of S. Heidelberg. Using molecular cut-off of 0-10 SNVs, the cgSNV analysis of 246 clinical and food isolates of S. Heidelberg collected in Québec, in the same year of the outbreak event, revealed that retail and abattoir chicken isolates likely represent an important source of human infection to S. Heidelberg. Interestingly, the isolates genetically related by cgSNV also harbored the same CRISPR as outbreak isolates and clusters. This indicates that CRISPR profiles can be useful as a complementary approach to determine source attribution in foodborne outbreaks. Use of the genomic analysis also allowed to identify a large number of cases that were missed by PFGE, indicating that most outbreaks are probably underestimated. Although epidemiological information must still support WGS-based results, cgSNV method is a highly discriminatory method for the resolution of outbreak events and the attribution of these events to their respective sources. CRISPR typing can serve as a complimentary tool to this analysis during source tracking.

The Role of the Salmonella spvB IncF Plasmid and Its Resident Entry Exclusion Gene traS on Plasmid Exclusion

Salmonella enterica cause significant illnesses worldwide. There has been a marked increase in resistance to fluoroquinolones and β-lactams/cephalosporins, antibiotics commonly used to treat salmonellosis. However, S. enterica serovars vary in their resistance to these and other antibiotics. The systemic virulence of some Salmonella serovars is due to a low copy number, IncF plasmid (65-100 kb) that contains the ADP-ribosylating toxin, SpvB. This virulence plasmid is present in only nine Salmonella serovars. It is possible that the spvB-virulence plasmid excludes other plasmids and may explain why antibiotic resistance is slow to develop in certain Salmonella serovars such as S. Enteritidis. The distribution of plasmid entry exclusion genes traS/traT and traY/excA are variable in Salmonella IncF and IncI plasmids, respectively and may account for differences in emergent antimicrobial resistance for some Salmonella serovars. The goal of this study is to determine the contribution of the Salmonella spvB-virulence plasmid in F-plasmid exclusion. From conjugation experiments, S. Typhimurium exhibited lower conjugation frequency with incFI and incFII plasmids when the spvB-virulence plasmid is present. Furthermore, introduction of cloned incFI traS into a "plasmidless" S. Typhimurium LT2 strain and Escherichia coli DH5α excluded incFI plasmid. However, deletion of the virulence plasmid traS did not affect plasmid exclusion significantly compared to a spvB control deletion. In addition, differences in F plasmid conjugation in natural Salmonella isolates did not correlate with IncF or SpvB-virulence plasmid genotype. There appear to be other plasmid or chromosomal genes at play in plasmid exclusion that may be responsible for the slow development of antibiotic resistance in certain serovars.

Genetic Diversity of Non-O157 Shiga Toxin-Producing Escherichia coli Recovered From Patients in Michigan and Connecticut

Shiga toxin-producing Escherichia coli (STEC) are important foodborne pathogens and non-O157 serotypes have been gradually increasing in frequency. The non-O157 STEC population is diverse and is often characterized using serotyping and/or multilocus sequence typing (MLST). Although spacers within clustered regularly interspaced repeat (CRISPR) regions were shown to comprise horizontally acquired DNA elements, this region does not actively acquire spacers in STEC. Hence, it is useful for further characterizing non-O157 STEC and examining relationships between strains. Our study goal was to evaluate the genetic relatedness of 41 clinical non-O157 isolates identified in Michigan between 2001 and 2005 while comparing to 114 isolates from Connecticut during an overlapping time period. Whole genome sequencing (WGS) was performed, and sequences were extracted for serotyping, MLST and CRISPR analysis. Phylogenetic analysis of MLST and CRISPR data was performed using the Neighbor joining and unweighted pair group method with arithmetic mean (UPGMA) algorithms, respectively. In all, 29 serogroups were identified; eight were unique to Michigan and 13 to Connecticut. “Big-six” serogroup frequencies were similar by state (Michigan: 73.2%, Connecticut: 81.6%), though STEC O121 was not found in Michigan. The distribution of sequence types (STs) and CRISPR profiles was also similar across states. Interestingly, big-six serogroups such as O103 and O26, grouped into different STs located on distinct branches of the phylogeny, further confirming that serotyping alone is not adequate for evaluating strain relatedness. Comparatively, the CRISPR analysis identified 361 unique spacers that grouped into 80 different CRISPR profiles. CRISPR spacers 231 and 317 were isolated from 79.2% (n = 118) and 59.1% (n = 88) of strains, respectively, regardless of serogroup and ST. Spacer profiles clustered according to the MLST analysis, though some discrepancies were noted. Indeed, use of both MLST and CRISPR typing enhanced the discriminatory power when compared to the use of each tool separately. These data highlight the genetic diversity of clinical STEC from different locations and show that CRISPR profiling can be used alongside MLST to discriminate related strains. Use of targeted sequencing approaches are particularly helpful for sites without WGS capabilities and can help define which strains require additional characterization using more discriminatory methods.

CRISPR-cas3 of Salmonella Upregulates Bacterial Biofilm Formation and Virulence to Host Cells by Targeting Quorum-Sensing Systems

Salmonella is recognized as one of the most common microbial pathogens worldwide. The bacterium contains the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) systems, providing adaptive immunity against invading foreign nucleic acids. Previous studies suggested that certain bacteria employ the Cas proteins of CRISPR-Cas systems to target their own genes, which also alters the virulence during invasion of mammals. However, whether CRISPR-Cas systems in Salmonella have similar functions during bacterial invasion of host cells remains unknown. Here, we systematically analyzed the genes that are regulated by Cas3 in a type I-E CRISPR-Cas system and the virulence changes due to the deletion of cas3 in Salmonella enterica serovar Enteritidis. Compared to the cas3 gene wild-type (cas3 WT) Salmonella strain, cas3 deletion upregulated the lsrFGBE genes in lsr (luxS regulated) operon related to quorum sensing (QS) and downregulated biofilm-forming-related genes and Salmonella pathogenicity island 1 (SPI-1) genes related to the type three secretion system (T3SS). Consistently, the biofilm formation ability was downregulated in the cas3 deletion mutant (Δcas3). The bacterial invasive and intracellular capacity of Δcas3 to host cells was also reduced, thereby increasing the survival of infected host cells and live chickens. By the transcriptome-wide screen (RNA-Seq), we found that the cas3 gene impacts a series of genes related to QS, the flagellum, and SPI-1-T3SS system, thereby altering the virulence phenotypes. As QS SPI-1-T3SS and CRISPR-Cas systems are widely distributed in the bacteria kingdom, our findings extend our understanding of virulence regulation and pathogenicity in mammalian hosts for Salmonella and potentially other bacteria.

Whole-genome sequencing reveals resistome of highly drug-resistant retail meat and human Salmonella Dublin

Non-typhoidal Salmonella (NTS) are a significant source of foodborne illness worldwide, with disease symptoms most often presenting as self-limiting gastroenteritis; however, occasionally the infection spreads and becomes invasive, frequently requiring anti-microbial treatment. The cattle-adapted Dublin serovar of NTS has commonly been associated with invasive illness and anti-microbial resistance (AMR). Here, the enhanced resolution conferred by whole-genome sequencing was utilized to elucidate and compare the resistome and genetic relatedness of 14 multidrug-resistant (MDR) and one pan-susceptible S. Dublin, isolated primarily in Pennsylvania, from fresh retail meat (one isolate) and humans (14 isolates). Twelve different genetic AMR determinants, including both acquired and chromosomal, were identified. Furthermore, comparative plasmid analysis indicated that AMR was primarily conferred by a putative IncA/C2 plasmid. A single pan-susceptible S. Dublin isolate, collected from the same timeframe and geographical region as the MDR isolates, did not carry an IncA/C2 replicon sequence within its genome. Moreover, the pan-susceptible isolate was genetically distinct from its MDR counterparts, as it was separated by ≥267 single nucleotide polymorphisms (SNPs), whereas there was a ≤38 SNP distance between the MDR isolates. Collectively, this data set advances our understanding of the genetic basis of the highly drug-resistant nature of S. Dublin, a serovar with significant public health implications.