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.
Citation: Susan Nadin-Davis, Louise Pope, John Devenish, Ray Allain, Dele Ogunremi. 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[J]. AIMS Microbiology, 2022, 8(3): 300-317. doi: 10.3934/microbiol.2022022
[1] | Alexey V. Rakov, Natalya A. Kuznetsova, Anatoly A. Yakovlev . Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles. AIMS Microbiology, 2020, 6(2): 106-120. doi: 10.3934/microbiol.2020007 |
[2] | Ogueri Nwaiwu, Chiugo Claret Aduba . An in silico analysis of acquired antimicrobial resistance genes in Aeromonas plasmids. AIMS Microbiology, 2020, 6(1): 75-91. doi: 10.3934/microbiol.2020005 |
[3] | Zeling Xu, Shuzhen Chen, Weiyan Wu, Yongqi Wen, Huiluo Cao . Type I CRISPR-Cas-mediated microbial gene editing and regulation. AIMS Microbiology, 2023, 9(4): 780-800. doi: 10.3934/microbiol.2023040 |
[4] | Mohamed-Yousif Ibrahim Mohamed, Ihab Habib, Hazim O. Khalifa . Salmonella in the food chain within the Gulf Cooperation Council countries. AIMS Microbiology, 2024, 10(3): 468-488. doi: 10.3934/microbiol.2024023 |
[5] | Ashrafus Safa, Jinath Sultana Jime, Farishta Shahel . Cholera toxin phage: structural and functional diversity between Vibrio cholerae biotypes. AIMS Microbiology, 2020, 6(2): 144-151. doi: 10.3934/microbiol.2020009 |
[6] | Bahram Nikmanesh, Kazem Ahmadikia, Muhammad Ibrahim Getso, Sanaz Aghaei Gharehbolagh, Shima Aboutalebian, Hossein Mirhendi, Shahram Mahmoudi . Candida africana and Candida dubliniensis as causes of pediatric candiduria: A study using HWP1 gene size polymorphism. AIMS Microbiology, 2020, 6(3): 272-279. doi: 10.3934/microbiol.2020017 |
[7] | Oluwafolajimi Adesanya, Tolulope Oduselu, Oluwawapelumi Akin-Ajani, Olubusuyi M. Adewumi, Olusegun G. Ademowo . An exegesis of bacteriophage therapy: An emerging player in the fight against anti-microbial resistance. AIMS Microbiology, 2020, 6(3): 204-230. doi: 10.3934/microbiol.2020014 |
[8] | Maria Parapouli, Anastasios Vasileiadis, Amalia-Sofia Afendra, Efstathios Hatziloukas . Saccharomyces cerevisiae and its industrial applications. AIMS Microbiology, 2020, 6(1): 1-31. doi: 10.3934/microbiol.2020001 |
[9] | Alaa Fathalla, Amal Abd el-mageed . Salt tolerance enhancement Of wheat (Triticum Asativium L) genotypes by selected plant growth promoting bacteria. AIMS Microbiology, 2020, 6(3): 250-271. doi: 10.3934/microbiol.2020016 |
[10] | Moh. A'inurrofiqin, Endang Sutriswati Rahayu, Dian Anggraini Suroto, Tyas Utami, Yunika Mayangsari . Safety assessment of the indigenous probiotic strain Lactiplantibacillus plantarum subsp. plantarum Kita-3 using Sprague–Dawley rats as a model. AIMS Microbiology, 2022, 8(4): 403-421. doi: 10.3934/microbiol.2022028 |
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.
Foodborne gastroenteritis is a major public health concern globally, resulting in significant morbidity, mortality and economic losses [1], with certain serovars of Salmonella enterica subsp. enterica figuring prominently [2] due to their wide dispersal in water sources used in agriculture [3]. In many developed countries, including Canada and the USA, S. Enteritidis is responsible for a significant proportion of all laboratory-confirmed cases of human salmonellosis [4]–[6] and has been associated with several large disease outbreaks resulting from contamination of fresh produce and poultry products [1],[7]–[10]. Indeed, poultry meat and egg products have been recognised as a significant exposure risk [8],[11] and have been the subject of major food recalls [12]. National programmes aim to prevent contamination of products at the source [13], but when these fail, timely identification of the food source of the disease-causing bacteria and removal of the product from the food supply system are instrumental in limiting the impact on public health.
A key requirement for the successful implementation of an investigation into foodborne salmonellosis is the accurate identification of the bacterial strain epidemiologically linked with both food contamination and human disease. To this end, several subtyping methods have been developed with the aim of devising cost-effective protocols that are timely and reproducible in a variety of laboratory settings [14]. One of the earliest subtyping schemes was phage typing, in which strains are evaluated for susceptibility to a panel of bacteriophages [15]. However, this method could be undertaken by only a limited number of laboratories due to the need for access to a panel of well characterised phages; furthermore, it has been found to lack robust epidemiological capability [16]. The development of pulsed field gel electrophoresis (PFGE), which characterises whole genomes based on the band patterns produced upon treatment of DNA with selected restriction endonucleases, provided a useful alternative that, until very recently was the gold standard technique, widely applied to Salmonella in food outbreak investigations [17],[18]. However, it has limited capability to differentiate S. Enteritidis strains. An alternative subtyping technique, multiple-locus variable number tandem repeat analysis (MLVA), which scores the number of short tandem repeats present at several different defined loci throughout the genome, has been developed into a standardised protocol employed for epidemiological investigations of S. Enteritidis worldwide, often in conjunction with PFGE [19],[20]. Multi-locus sequence typing (MLST), which involves determination of the nucleotide sequences for a set of housekeeping genes, was proposed as a novel means of discriminating between closely related micro-organisms [21]. Its use as an alternative to serotyping for identification of S. enterica strains was proposed [22], but it was insufficiently discriminatory for use in foodborne outbreak investigations.
Since the more widespread availability of whole genome sequencing (WGS), there has been a shift towards replacement of these traditional tools by in silico analyses of WGS data to enable both serovar identification [23] and highly sensitive strain subtyping through methods such as core genome MLST [24]. While WGS provides for the most comprehensive genetic characterisation of a strain, the technology remains reasonably expensive per strain, requires considerable bioinformatics, human and material resources, and may take many days to complete. Methods which use WGS data to devise simpler and faster tools for strain subtyping have thus also been developed. These include an SNP-based approach in which nucleotide polymorphisms at several specific locations across the genome are evaluated by a PCR-based assay to provide a dataset used for phylogenetic assessment [25].
Yet another subtyping approach which has gained interest for application to many bacterial pathogens involves sequence analysis of clustered regularly interspaced short palindromic repeats (CRISPR) loci [26]. This sequence element is considered part of an innate immune system utilised by many bacteria, including Salmonella, to prevent invasion by mobile genetic elements and phages through the identification and subsequent degradation of foreign DNA [27]. Each CRISPR locus consists of several direct repeat sequences which are interspersed by spacer sequences derived from foreign DNA. A set of CRISPR-associated (Cas) genes is responsible for mediating this protective mechanism through processing of non-native DNA into short segments that are integrated as spacer sequences into the CRISPR locus. Following transcription of this locus, one or more cas proteins cleave the RNA at the spacer sequences, to generate small interfering RNA/protein complexes that identify complementary unintegrated DNA sequences and elicit their subsequent cleavage and degradation. This process thereby prevents persistence and propagation of foreign genetic material within the cell [28]. The acquisition or loss of spacer sequences in the CRISPR locus reflects the challenge history of the strain, and thus characterisation of this sequence element would be expected to provide a sensitive subtyping target. Indeed, comparative genomics of S. enterica strains, which harbour two CRISPR loci, suggested the importance of CRISPR-mediated immunity in regulating gene acquisition from mobile genetic elements and thereby influencing lineage evolution and diversification [29]. Studies showing that CRISPR polymorphism correlated well with Salmonella serotype and MLST designation also suggested the utility of these loci for serovar subtyping [30],[31]. A method designated CRISPR-including multi-virulence-locus sequence typing (CRISPR-MLVST), which includes sequence characterisation of both CRISPR loci and two virulence genes (fimH and sseL), was described for subtyping of several Salmonella serovars, including S. Enteritidis [32],[33]. Similar protocols have been employed to subtype other serovars of significant public health importance, including S. Typhimurium and S. Heidelberg [31],[34]. An Australian study found that the CRISPR and prophage profiles of a collection of S. Typhimurium strains correlated with core genome evolution [35]. Subtyping of S. Enteritidis isolates obtained in China, by CRISPR analysis only, identified several subtypes recovered from different poultry production areas [36] and suggested the pathogen's transmission from swine and poultry to humans [37]. In contrast to these reports, other studies have concluded that CRISPR loci of many bacteria, including those of Salmonella, exhibit quite limited sequence variation and may no longer be actively accumulating spacer sequences, thus making these loci relatively poor epidemiological markers [38],[39].
To explore this issue further, this study used a collection of S. Enteritidis strains, recovered from a variety of environmental and food sources across Canada, to determine the value of strain subtyping by both CRISPR locus and fimH gene characterisation compared to the traditional methods of PFGE, PT and an SNP-PCR approach previously developed by our group [25]. CRISPR analysis alone was found to be the least sensitive of all the subtyping methods analysed, though when combined with other methods such as SNP-based typing, high levels of strain discrimination could be achieved.
This study included a total of 89 S. Enteritidis field strains and three reference strains, as detailed in Table 1. In addition, the well characterised S. Enteritidis reference strain P125109 (GenBank accession NC_011294) was included in the CRISPR analysis to confirm that our analyses yielded results in complete concordance with previously reported CRISPR sequences. Many of these strains were recovered from environmental swabs taken in poultry production facilities as part of an animal health programme to monitor these operations for Salmonella contamination. Culturing of these swabs and identification of S. Enteritidis were undertaken as previously described [40]. A few others were isolated from food samples or food processing facilities using standard Salmonella isolation procedures, comprising both pre-enrichment and enrichment steps and plating on selective agar plates [41]. Confirmation of serotype and phage typing were undertaken at the Salmonella reference laboratory of the Public Health Agency of Canada (PHAC) in Guelph, Ontario.
Pure isolates of S. Enteritidis were grown on tryptic soy agar plates and incubated at 36 °C for 18–22 hrs. Cells were scraped off the agar surface and resuspended in phosphate buffered saline prior to total DNA extraction using a Wizard genomic DNA isolation kit as per the supplier's directions (Promega, Madison, Wisconsin). Purified DNA solutions were quantified spectroscopically using a Nanovue instrument (GE Biosciences) and stored at −20 °C.
The two CRISPR loci, CRISPR-1 and CRISPR-2, were amplified separately by PCR using the primer pairs described by Liu and colleagues [33], as detailed in Table 2. Each 50 µL PCR reaction contained 1X PCR buffer, 1.5 mM MgCl2, 0.2 mM dNTP, 2.5 U Taq DNA polymerase (all supplied by Invitrogen Life Technologies), 1 ng of DNA template and 0.5 µM forward and reverse primers (synthesised by Integrated DNA Technologies, Coralville, Iowa). Amplifications were performed on a GeneAmp 9700 thermocycler (Applied Biosystems, ThermoFisher) using the following cycling conditions: an initial denaturation at 94 °C for 2 min followed by 30 cycles of 94 °C for 1 min, 52 °C for 1 min, 72 °C for 1 min 15 sec and a final extension at 72 °C for 5 min. Product generation was verified by standard gel electrophoresis, and amplicons were purified using a Wizard PCR purification kit (Promega, Madison, Wisconsin) as per the manufacturer's directions.
PCR products were sequenced in both directions using the PCR primers with a BDTv3.1 cycle sequencing kit (Applied Biosystems) and a 9700 thermocycler as per kit instructions. Sequencing products were purified using a BDXterminator purification kit (Applied Biosystems) and analysed on a 3500xl genetic analyser (Applied Biosystems). Forward and reverse sequence reads were assembled into a consensus sequence using the Lasergene v.11 software package (DNASTAR Inc., Madison, Wisconsin).
Sequence alignments for each of the two loci were performed using MEGA X [42]. CRISPR loci were analysed using the CRISPRFinder tool [43] available at https://crisprcas.i2bc.paris-saclay.fr to identify the conserved direct repeat (DR) motif and the intervening spacer sequences, each of which was identified numerically. For each strain, the complete CRISPR locus was assigned a single letter according to the spacer sequences it retained, and the combination of scores for both CRISPR loci yielded the final double letter strain type. For some samples, a number was also included to identify minor sequence variations to type. The sequence data from both loci were assembled into a concatenated alignment for generation of a phylogenetic tree using BioNumerics v. 6.01 software (Applied Maths, Sint-Martens-Latem, Belgium).
To explore if additional subtype information could be realised through sequence analysis of the virulence-associated fimH gene, this locus was amplified from genomic DNA prepared from 70 samples. The 1087 bp PCR product was generated using the fimH-F and fimH-R primer pair, which straddles the 1008 bp ORF [32] (Table 2), and procedures similar to those used for the CRISPR loci. Amplicons were sequenced from both DNA strands using fimH PCR and internal primers (fimH-Fint and fimH-Rint), and the reads were assembled using DNASTAR Lasergene v.11 software. An alignment of all assembled sequences was generated in MEGA X to identify all fimH gene SNPs using the P125109 strain sequence as the reference, which was designated as allele A. All other alleles thus identified were assigned a single letter designation from B to D (Table 3). Table 1 summarises the fimH types of all 70 samples thus analysed.
fimH allele | Sample(s) | Location and nature of SNP in fimH ORF | CRISPR subtype |
B | BUR-FI-2009 OLF10012-1 |
A49 to G49 G112 to A112 C259 to T259 C292 to T292 A730 to G730 C770 to T770 T794 to G794 |
DG |
C | 27655 | T466 to C466 | AG-1 |
D | BUR-FI-2005 | T878 to C878 | AE |
Differences from the sequence of the reference strain P125109 (positions 588156 to 589163), which was assigned as allele A, are indicated. The CRISPR subtype for each of the four samples is also indicated.
PFGE was performed using a standard method [44] modified as described [40]. Samples were electrophoresed for 20 hrs on the CHEF Mapper (Bio-Rad Laboratories, Mississauga, ON), and the data were analysed using the BioNumerics v6.01 software. PFGE patterns were assigned by PulseNet Canada (National Microbiology Laboratory, Winnipeg, Canada).
For each strain, a set of 60 base positions scattered throughout the genome was determined as described [25] and used to generate a fasta file. An alignment of these fasta files was used to construct a phylogenetic tree using the unweighted pair group method with arithmetic mean (UPGMA), with the MEGA X software employing 1000 bootstrap replicates.
The discriminatory capabilities of individual and combined methods were assessed using Simpson's index of diversity [45], calculated using the formula
where Ds is Simpson's index of diversity, N is the total number of samples, and ni is the number of samples in the ith group.
All 92 strains examined generated both CRISPR 1 and CRISPR 2 loci for sequencing. The sequence of the 29 base DR was highly conserved within the CRISPR 1 locus, while that of the CRISPR 2 locus exhibited much more variation, with base substitutions observed in several samples (Table 4). The CRISPR 1 locus comprised up to 11 spacer sequences, while the CRISPR 2 locus comprised up to 12, with significant sequence variations in spacers 3 and 4 for individual strains such that they are identified as 3′ and 4′. Full details of all spacer sequences, including all variations thereof, are presented in Table S1, while Figure 1 summarises the combinations of spacers that defined each CRISPR subtype. A total of 7 spacer combinations, labelled with a single letter code (A to G), were observed for each CRISPR locus (Figure 1). Combining the letter codes for each CRISPR locus generated a two-letter binary code, with additional small variations due to base substitutions indicated by a numerical code (1, 2, etc.) (Table S1). When both CRISPR loci sequences were combined, they identified a total of 16 distinct S. Enteritidis subtypes. In addition, the P125109 reference strain yielded sequences identical to those published previously for the two CRISPR loci and was assigned to the HB subtype based on our classification scheme.
CRISPR | ||
1 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | GTGTTTATCCCCGCTGACGCGGGGAACAC | |
CRISPR | ||
2 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | ACGGCTATCCTTGTTGGCGCGGGGAACAC | |
variation (single isolate) | ACGGCTATCCTGGTTGGCGCGGGGAACAC | |
Internal repeat variants | GGGTTTATYCCCGCTGGCGCGGGGAACAA | |
GGGTTTATCCCCGCTGGCGCGGGGAACAC | ||
CGGTTTATCCCCGCTGGCGAGGGGAACAC | ||
CGGTTTATCCCCGATGGCGCGGGGAACAC | ||
CGGTTTATCTCCGCTGGGGCGGGGAACAC | ||
CGGTCTATCCCCGCTGGCGCGGGGAACAC | ||
CGCTTTATCCCCGCTGGCGCGGGGAACAC | ||
Internal and 3′ terminal repeat | CAGTTTATCCCCGCTGGCGCGGGGAACAC | |
Bases different from those of each consensus sequence are underlined |
These 92 strains yielded 19 different PFGE patterns, with the majority falling into three groups: SENXAI.0003_ SENBNI.0003 (n = 35), SENXAI.0006_ SENBNI.0007 (n = 14) and SENXAI.0038_ SENBNI.0016 (n = 15). Six samples comprised the PFGE type SENXAI.0214_ SENBNI.0225, while the remaining 15 subtypes were distributed among 22 strains.
Of the 26 different phage types (PTs) identified amongst this collection, the most common were PT8 (n = 23), PT13 (n = 13), PT13a (n = 9) and PT23 (n = 13). Eleven atypical PTs were observed, two of which were unattributed to a group, and one sample could not be typed. Four strains appeared to be either of mixed type or to have modified their PT during assessment; for the purpose of Simpson's index of diversity calculation, as detailed below, the rarer type was employed.
SNP-PCR typing scored the base sequence of 60 positions scattered throughout the genome and subsequently used these concatenated data for phylogenetic tree generation. Based on hundreds of S. Enteritidis samples, this process has been reported to identify many distinct clades [25], of which 20 were identified in the sample collection (Figure 2).
Simpson's index of diversity was calculated for each subtyping method, as summarised in Table 5; the full calculation summary of these values is provided (Table S2). This clearly shows that, for the sample set examined in this study, the CRISPR analysis was the least discriminatory, identifying just 16 subtypes with an index of diversity of 0.7367. Although PT identified the most subtypes (n = 25), the SNP-PCR subtyping, which identified 20 subtypes, had an index of diversity of 0.9259 and was scored as the most discriminatory compared to the traditional methods of PFGE and PT.
When subtyping methods were paired, not surprisingly, their discriminatory capabilities increased; the combined discriminatory capabilities of the two traditional methods versus the two molecular methods are summarised (Table 5). Comparison of the subtyping assignments for individual samples did reveal some level of concordance between methods for some subtypes, as well as significant variations in subtype structure for many samples (Table 1 and Table S3). This is illustrated by the phylogenetic tree (Figure 3) representing the CRISPR subtypes for all samples, to which has been added the corresponding SNP-PCR data. For example, all eight strains with a CC CRISPR subtype were classified as SNP-PCR subtype 4, but most other CRISPR subtypes were scattered amongst several SNP-PCR subtypes. The two most dominant CRISPR subtypes, namely AA-1 (n = 39) and AG-1 (n = 26), were each classified into multiple separate SNP-PCR subtypes. The AA-1 CRISPR subtype included strains belonging to SNP-PCR subtypes 11, 14, 15, 16, 20, 22, 24, 25 and 26, whereas the AG-1 CRISPR subtype was composed of SNP-PCR subtypes 6, 7, 11, 13, 20, 28, 29 and 30. Conversely, some SNP-PCR subtypes are restricted to a specific CRISPR subtype (e.g., SNP-PCR subtypes 14, 15, 16, 22, 24 and 26 all correspond to CRISPR subtype AA-1), but this is not always the case. Often, each SNP-PCR subtype is composed of two or three CRISPR subtypes, e.g., SNP-PCR 5 (AE and DG CRISPR subtypes), SNP-PCR 6 (AB, AG-1 and GG CRISPR subtypes) and SNP-PCR 25 (AA-1, AA-2 and AG-3 CRISPR subtypes). As a result, when both typing methods were combined, the number of distinct subtypes increased to 31, with an index of diversity of 0.9439, which was the highest observed in this study. Combining PFGE and PT analysis yielded the greatest number of distinct subtypes at 45, yet the increased index of diversity of 0.9422 observed for the combined traditional methods was marginally lower than for the combined molecular methods of CRISPR/SNP-PCR subtyping.
As the sequencing of additional genes in concert with CRISPR analysis has been reported to improve S. Enteritidis subtype discrimination [33], an initial screen of 70 strains employed in this study, including the P125109 strain used as the reference, was undertaken to explore the utility of including the sequence of the 1008 bp fimH gene in the subtyping scheme. Based upon the fimH gene sequence of the reference, which was scored as allele A, only three other alleles of this gene were identified based on SNPs at nine positions (Tables 1 and 3). However, seven of these SNPs were associated with the two distinctive strains with a DG CRISPR profile, one was found in the single strain with an AE CRISPR profile, and the last SNP was found in a single strain of CRISPR type AG-1. As a result, this additional sequence analysis added only a single subtype to this set of 70 samples.
The primary application of bacterial typing in disease control is to reveal epidemiological links between isolates responsible for human illness and the source of the causative agent. For outbreaks of foodborne gastroenteritis, identifying the food product contaminated with a specific bacterial pathogen enables its recall and limitation of the disease outbreak. Such investigations require timely action on the part of public health officials, thereby limiting the extent of the outbreak; this in turn requires efficient laboratory identification of the strain of the bacterium involved. Given the importance of S. Enteritidis in human foodborne disease, this study has explored sequence typing of the two CRISPR loci as a rapid method of potential utility for such investigations.
However, CRISPR analysis of the 92 strains included in this study was less discriminatory compared to the traditional typing methods of PFGE and PT. While CRISPR typing did clearly distinguish strains recovered from distinct food sources such as seafood (clams), a high number of isolates fell into two CRISPR groups, AA-1 (n = 39) and AG-1 (n = 26), representing 42.4% and 28.3%, respectively, of all samples. With few exceptions, these isolates were all strongly associated with poultry production facilities or poultry food products, often with multiple samples recovered from single submissions. Given that these are the products primarily responsible for many human infections, this clearly limits the ability of CRISPR-based typing alone to locate the precise source of many food contaminants. Clearly the interpretation of the value of this approach will depend on the range of isolates requiring discrimination.
Comparison of the subtyping assignments by these different methods revealed very limited concordance in their grouping patterns, except for a set of eight samples that were clearly discriminated by their CRISPR CC subtype, a SENBNI.0225 PFGE profile, a PT of 9b or a closely related atypical subtype and SNP subtype 4. The close evolutionary relationship of these samples, which all originated from duck producing facilities and which harbour an unusual virulence plasmid, has been reported previously [46]. The lack of similar concordance in the subtyping profiles for most other samples reflects the independent nature of the features being targeted, and, while this does not limit their utility for strain identification, it suggests that in general many of these targets may not reveal the evolutionary paths of the strains tested. Nevertheless, the SNP-PCR method, when applied to a large population of S. Enteritidis (n = 1,227), provides insights into the genetic structure of the organism [25].
Some studies, in which CRISPR sequence analysis has been combined with that of two additional genes, fimH and sseL (CRISPR-MVLST), have reported satisfactory discriminatory capability [33],[34]. Indeed, given that the CRISPR-MVLST typing scheme for S. Enteritidis has been reported to be more discriminatory than PFGE [47], the value of including additional sequence information for these two genes was explored. However, neither target was found to significantly improve discriminatory capability for this sample cohort. The fimH locus identified just a single additional type out of 70 of the 92 isolates of this study, while review of whole genome sequence data previously reported for some of these isolates [48] indicated that the sseL gene was highly conserved. Given the additional effort involved in the amplification and sequencing of these targets, their inclusion in this subtyping scheme was not considered worthwhile.
Combination of CRISPR-MVLST with PFGE has been reported to improve discriminatory capability for clinical S. Enteritidis isolates [49]. Indeed, in this study the CRISPR/SNP-PCR combination, which identified 31 groups, yielded a marginally higher index of diversity than the PFGE/PT subtyping combination, even though the latter identified a much higher number of groups. However, the significant technical challenges posed by the PT method and concerns over the epidemiological significance of this subtyping tool, due to reports that some Salmonella serovars can change PT [50],[51], undermine its value for timely food borne illness investigations. Indeed, this problem was identified in four samples of this study for which a single distinct PT result was not obtained. This suggests a high number of groups, in itself, may have limited biological value, and other measures of sample diversity may be more meaningful.
This study provides a cautionary note on the use of CRISPR characterisation for S. Enteritidis subtyping. While CRISPR analysis alone did provide some discriminatory capability, in isolation this method was inferior to the other subtyping methods examined in this study. Even with the addition of fimH gene sequencing, it was not found to be a sufficiently sensitive subtyping tool. Indeed, this study reinforces prior observations that these loci are relatively conserved across this serovar and thus do not appear to reflect recent acquisition of spacer sequences [38], thereby diminishing their use for high resolution epidemiological studies. Combining SNP-PCR and CRISPR sequence determination, analyses more readily performed by the standard microbiology laboratory and in a more timely fashion than the traditional techniques, would offer an alternative to those laboratories lacking the necessary hardware or bioinformatics support for whole genome sequence analysis.
[1] |
Carstens CK, Salazar JK, Darkoh C (2019) Multistate outbreaks of foodborne illness in the United States associated with fresh produce from 2010 to 2017. Front Microbiol 10: 2667-2667. https://doi.org/10.3389/fmicb.2019.02667 ![]() |
[2] |
Majowicz S, Musto J, Scallan E, et al. (2010) The global burden of nontyphoidal Salmonella. Clin Infect Dis 50: 882-889. https://doi.org/10.1086/650733 ![]() |
[3] |
Liu H, Whitehouse CA, Li B (2018) Presence and persistence of Salmonella in water: the impact on microbial quality of water and food safety. Front Public Health 6: 159-159. https://doi.org/10.3389/fpubh.2018.00159 ![]() |
[4] | (2018) PHACNational Enteric Surveillance Program Annual Summary 2016. Guelph, Ontario: Public Health Agency, Government of Canada. |
[5] |
Chai S, White P, Lathrop S, et al. (2012) Salmonella enterica serotype Enteritidis: increasing incidence of domestically acquired infections. Clin Infect Dis 54: S488-497. https://doi.org/10.1093/cid/cis231 ![]() |
[6] |
Lane C, LeBaigue S, Esan O, et al. (2014) Salmonella enterica serovar Enteritidis, England and Wales, 1945–2011. Emerging Infect Dis 20: 1097-1104. https://doi.org/10.3201/eid2007.121850 ![]() |
[7] |
Kozak GK, MacDonald D, Landry L, et al. (2013) Foodborne outbreaks in Canada linked to produce: 2001 through 2009. J Food Prot 76: 173-183. https://doi.org/10.4315/0362-028X.JFP-12-126 ![]() |
[8] |
Chai SJ, Cole D, Nisler A, et al. (2017) Poultry: the most common food in outbreaks with known pathogens, United States, 1998–2012. Epidemiol Infect 145: 316-325. https://doi.org/10.1017/S0950268816002375 ![]() |
[9] |
Middleton D, Savage R, Tighe M, et al. (2014) Risk factors for sporadic domestically acquired Salmonella serovar Enteritidis infections: a case-control study in Ontario, Canada, 2011. Epidemiol Infect 142: 1411-1421. https://doi.org/10.1017/S0950268813001945 ![]() |
[10] |
Nesbitt A, Ravel A, Murray R, et al. (2012) Integrated surveillance and potential sources of Salmonella Enteritidis in human cases in Canada from 2003 to 2009. Epidemiol Infect 140: 1757-1772. https://doi.org/10.1017/S0950268811002548 ![]() |
[11] |
Taylor M, Leslie M, Ritson M, et al. (2012) Investigation of the concurrent emergence of Salmonella enteritidis in humans and poultry in British Columbia, Canada, 2008-2010. Zoonoses Public Health 59: 584-892. https://doi.org/10.1111/j.1863-2378.2012.01500.x ![]() |
[12] |
Kuehn B (2010) Salmonella cases traced to egg producers. J Am Medl Assoc 304: 1316. https://doi.org/10.1001/jama.2010.1330 ![]() |
[13] |
Chousalkar K, Gast R, Martelli F, et al. (2018) Review of egg-related salmonellosis and reduction strategies in United States, Australia, United Kingdom and New Zealand. Crit Rev Microbiol 44: 290-303. https://doi.org/10.1080/1040841X.2017.1368998 ![]() |
[14] |
Tang S, Orsi RH, Luo H, et al. (2019) Assessment and comparison of molecular subtyping and characterization methods for Salmonella. Front Microbiol 10: 1591. https://doi.org/10.3389/fmicb.2019.01591 ![]() |
[15] |
Ward LR, de Sa JDH, Rowe B (1987) A phage-typing scheme for Salmonella enteritidis. Epidemiol Infect 99: 291-294. https://doi.org/10.1017/S0950268800067765 ![]() |
[16] |
Carrique-Mas JJ, Papadopoulou C, Evans SJ, et al. (2008) Trends in phage types and antimicrobial resistance of Salmonella enterica serovar Enteritidis isolated from animals in Great Britain from 1990 to 2005. Vet Rec 162: 541-546. https://doi.org/10.1136/vr.162.17.541 ![]() |
[17] |
Peters TM, Berghold C, Brown D, et al. (2007) Relationship of pulsed-field profiles with key phage types of Salmonella enterica serotype Enteritidis in Europe: results of an international multi-centre study. Epidemiol Infect 135: 1274-1281. https://doi.org/10.1017/S0950268807008102 ![]() |
[18] |
Ribot EM, Fair MA, Gautom R, et al. (2006) Standardization of pulsed-field gel electrophoresis protocols for the subtyping of Escherichia coli O157:H7, Salmonella, and Shigella for PulseNet. Foodborne Pathog Dis 3: 59-67. https://doi.org/10.1089/fpd.2006.3.59 ![]() |
[19] |
Nadon CA, Trees E, Ng LK, et al. (2013) Development and application of MLVA methods as a tool for inter-laboratory surveillance. Eurosurveillance 18: 20565. https://doi.org/10.2807/1560-7917.ES2013.18.35.20565 ![]() |
[20] |
Bertrand S, De Lamine de Bex G, Wildemauwe C, et al. (2015) Multi Locus Variable-Number Tandem Repeat (MLVA) typing tools improved the surveillance of Salmonella Enteritidis: A 6 Years Retrospective Study. PLoS One 10: e0117950. https://doi.org/10.1371/journal.pone.0117950 ![]() |
[21] |
Maiden MCJ, Bygraves JA, Feil E, et al. (1998) Multilocus sequence typing: A portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci 95: 3140-3145. https://doi.org/10.1073/pnas.95.6.3140 ![]() |
[22] |
Achtman M, Wain J, Weill F-X, et al. (2012) Multilocus sequence typing as a replacement for serotyping in Salmonella enterica. PLoS Pathog 8: e1002776. https://doi.org/10.1371/journal.ppat.1002776 ![]() |
[23] |
Yoshida C, Kruczkiewicz P, Laing C, et al. (2016) The Salmonella In Silico Typing Resource (SISTR): an open web-accessible tool for rapidly typing and subtyping draft Salmonella genome assemblies. PLoS One 11: e0147101. https://doi.org/10.1371/journal.pone.0147101 ![]() |
[24] |
Pearce ME, Alikhan NF, Dallman TJ, et al. (2018) Comparative analysis of core genome MLST and SNP typing within a European Salmonella serovar Enteritidis outbreak. Int J Food Microbiol 274: 1-11. https://doi.org/10.1016/j.ijfoodmicro.2018.02.023 ![]() |
[25] |
Ogunremi D, Kelly H, Dupras AA, et al. (2014) Development of a new molecular subtyping tool for Salmonella enterica serovar Enteritidis based on single nucleotide polymorphism genotyping using PCR. J Clin Microbiol 52: 4275-4285. https://doi.org/10.1128/JCM.01410-14 ![]() |
[26] |
Shariat N, Dudley EG (2014) CRISPRs: molecular signatures used for pathogen subtyping. Appl Environ Microbiol 80: 430-439. https://doi.org/10.1128/AEM.02790-13 ![]() |
[27] |
Horvath P, Barrangou R (2010) CRISPR/Cas, the immune system of bacteria and archaea. Science 327: 167-170. https://doi.org/10.1126/science.1179555 ![]() |
[28] |
Burmistrz M, Krakowski K, Krawczyk-Balska A (2020) RNA-targeting CRISPR-Cas systems and their applications. Int J Mol Sci 21: 1122. https://doi.org/10.3390/ijms21031122 ![]() |
[29] |
Fricke WF, Mammel MK, McDermott PF, et al. (2011) Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPR-mediated adaptive sublineage evolution. J Bacteriol 193: 3556-3568. https://doi.org/10.1128/JB.00297-11 ![]() |
[30] |
Fabre L, Le Hello S, Roux C, et al. (2014) CRISPR Is an optimal target for the design of specific PCR assays for Salmonella enterica serotypes Typhi and Paratyphi A. PLoS Neglected Trop Dis 8: e2671. https://doi.org/10.1371/journal.pntd.0002671 ![]() |
[31] |
Fabre L, Zhang J, Guigon G, et al. (2012) CRISPR typing and subtyping for improved laboratory surveillance of Salmonella infections. PLoS ONE 7: e36995. https://doi.org/10.1371/journal.pone.0036995 ![]() |
[32] |
Liu F, Barrangou R, Gerner-Smidt P, et al. (2011) Novel virulence gene and Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) multilocus sequence typing scheme for subtyping of the major serovars of Salmonella enterica subsp. enterica. Appl Environ Microbiol 77: 1946-1956. https://doi.org/10.1128/AEM.02625-10 ![]() |
[33] |
Liu F, Kariyawasam S, Jayarao BM, et al. (2011) Subtyping Salmonella enterica serovar Enteritidis isolates from different sources by using sequence typing based on virulence genes and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs). Appl Environ Microbiol 77: 4520-4526. https://doi.org/10.1128/AEM.00468-11 ![]() |
[34] |
Shariat N, Sandt CH, DiMarzio MJ, et al. (2013) CRISPR-MVLST subtyping of Salmonella enterica subsp. enterica serovars Typhimurium and Heidelberg and application in identifying outbreak isolates. BMC Microbiol 13: 254-254. https://doi.org/10.1186/1471-2180-13-254 ![]() |
[35] |
Fu S, Hiley L, Octavia S, et al. (2017) Comparative genomics of Australian and international isolates of Salmonella Typhimurium: correlation of core genome evolution with CRISPR and prophage profiles. Sci Rep 7: 9733-9733. https://doi.org/10.1038/s41598-017-06079-1 ![]() |
[36] |
Fei X, He X, Guo R, et al. (2017) Analysis of prevalence and CRISPR typing reveals persistent antimicrobial-resistant Salmonella infection across chicken breeder farm production stages. Food Control 77: 102-109. https://doi.org/10.1016/j.foodcont.2017.01.023 ![]() |
[37] |
Li Q, Wang X, Yin K, et al. (2018) Genetic analysis and CRISPR typing of Salmonella enterica serovar Enteritidis from different sources revealed potential transmission from poultry and pig to human. Int J Food Microbiol 266: 119-125. https://doi.org/10.1016/j.ijfoodmicro.2017.11.025 ![]() |
[38] |
Shariat N, Timme RE, Pettengill JB, et al. (2015) Characterization and evolution of Salmonella CRISPR-Cas systems. Microbiology 161: 374-386. https://doi.org/10.1099/mic.0.000005 ![]() |
[39] |
Touchon M, Rocha EPC (2010) The small, slow and specialized CRISPR and anti-CRISPR of Escherichia and Salmonella. PLoS ONE 5: e11126. https://doi.org/10.1371/journal.pone.0011126 ![]() |
[40] |
Nadin-Davis S, Pope L, Ogunremi D, et al. (2019) A real-time PCR regimen for testing environmental samples for Salmonella enterica subsp. enterica serovars of concern to the poultry industry, with special focus on Salmonella Enteritidis. Can J Microbiol 65: 162-173. https://doi.org/10.1139/cjm-2018-0417 ![]() |
[41] | Reid A Isolation and identification of Salmonella from food and environmental samples: Standard operating procedure MFHPB 20. Microbiological Methods Committee, Bureau of Microbial Hazards, Health Canada (2009). |
[42] |
Kumar S, Stecher G, Li M, et al. (2018) MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol 35: 1547-1549. https://doi.org/10.1093/molbev/msy096 ![]() |
[43] |
Couvin D, Bernheim A, Toffano-Nioche C, et al. (2018) CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res 46: W246-W251. https://doi.org/10.1093/nar/gky425 ![]() |
[44] | PulseNetInternationalStandard Operating Procedure for PulseNet PFGE of Escherichia coli O157:H7, Escherichia coli non-O157 (STEC), Salmonella serotypes, Shigella sonnei and Shigella flexneri. (2013). |
[45] |
Hunter PR, Gaston MA (1988) Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol 26: 2465-2466. https://doi.org/10.1128/jcm.26.11.2465-2466.1988 ![]() |
[46] |
Nadin-Davis S, Pope L, Chmara J, et al. (2020) An unusual Salmonella Enteritidis strain carrying a modified virulence plasmid lacking the prot6e gene represents a geographically widely distributed lineage. Front Microbiol 11: 1322. https://doi.org/10.3389/fmicb.2020.01322 ![]() |
[47] |
Deng X, Shariat N, Driebe EM, et al. (2015) Comparative analysis of subtyping methods against a whole-genome-sequencing standard for Salmonella enterica serotype Enteritidis. J Clin Microbiol 53: 212-218. https://doi.org/10.1128/JCM.02332-14 ![]() |
[48] |
Ogunremi D, Devenish J, Amoako K, et al. (2014) High resolution assembly and characterization of genomes of Canadian isolates of Salmonella Enteritidis. BMC Genomics 15: 713. https://doi.org/10.1186/1471-2164-15-713 ![]() |
[49] |
Shariat N, DiMarzio MJ, Yin S, et al. (2013) The combination of CRISPR-MVLST and PFGE provides increased discriminatory power for differentiating human clinical isolates of Salmonella enterica subsp. enterica serovar Enteritidis. Food Microbiol 34: 164-173. https://doi.org/10.1016/j.fm.2012.11.012 ![]() |
[50] |
Hiley L, Graham RMA, Jennison AV (2021) Characterisation of IncI1 plasmids associated with change of phage type in isolates of Salmonella enterica serovar Typhimurium. BMC Microbiol 21: 92-92. https://doi.org/10.1186/s12866-021-02151-z ![]() |
[51] |
Tankouo-Sandjong B, Kinde H, Wallace I (2012) Development of a sequence typing scheme for differentiation of Salmonella Enteritidis strains. FEMS Microbiol Lett 331: 165-175. https://doi.org/10.1111/j.1574-6968.2012.02568.x ![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
1. | Ruimin Gao, Jasmine Rae Frost, 2023, 10.5772/intechopen.109712 |
fimH allele | Sample(s) | Location and nature of SNP in fimH ORF | CRISPR subtype |
B | BUR-FI-2009 OLF10012-1 |
A49 to G49 G112 to A112 C259 to T259 C292 to T292 A730 to G730 C770 to T770 T794 to G794 |
DG |
C | 27655 | T466 to C466 | AG-1 |
D | BUR-FI-2005 | T878 to C878 | AE |
Differences from the sequence of the reference strain P125109 (positions 588156 to 589163), which was assigned as allele A, are indicated. The CRISPR subtype for each of the four samples is also indicated.
CRISPR | ||
1 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | GTGTTTATCCCCGCTGACGCGGGGAACAC | |
CRISPR | ||
2 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | ACGGCTATCCTTGTTGGCGCGGGGAACAC | |
variation (single isolate) | ACGGCTATCCTGGTTGGCGCGGGGAACAC | |
Internal repeat variants | GGGTTTATYCCCGCTGGCGCGGGGAACAA | |
GGGTTTATCCCCGCTGGCGCGGGGAACAC | ||
CGGTTTATCCCCGCTGGCGAGGGGAACAC | ||
CGGTTTATCCCCGATGGCGCGGGGAACAC | ||
CGGTTTATCTCCGCTGGGGCGGGGAACAC | ||
CGGTCTATCCCCGCTGGCGCGGGGAACAC | ||
CGCTTTATCCCCGCTGGCGCGGGGAACAC | ||
Internal and 3′ terminal repeat | CAGTTTATCCCCGCTGGCGCGGGGAACAC | |
Bases different from those of each consensus sequence are underlined |
fimH allele | Sample(s) | Location and nature of SNP in fimH ORF | CRISPR subtype |
B | BUR-FI-2009 OLF10012-1 |
A49 to G49 G112 to A112 C259 to T259 C292 to T292 A730 to G730 C770 to T770 T794 to G794 |
DG |
C | 27655 | T466 to C466 | AG-1 |
D | BUR-FI-2005 | T878 to C878 | AE |
CRISPR | ||
1 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | GTGTTTATCCCCGCTGACGCGGGGAACAC | |
CRISPR | ||
2 | Consensus repeat | CGGTTTATCCCCGCTGGCGCGGGGAACAC |
5′ terminal repeat | ACGGCTATCCTTGTTGGCGCGGGGAACAC | |
variation (single isolate) | ACGGCTATCCTGGTTGGCGCGGGGAACAC | |
Internal repeat variants | GGGTTTATYCCCGCTGGCGCGGGGAACAA | |
GGGTTTATCCCCGCTGGCGCGGGGAACAC | ||
CGGTTTATCCCCGCTGGCGAGGGGAACAC | ||
CGGTTTATCCCCGATGGCGCGGGGAACAC | ||
CGGTTTATCTCCGCTGGGGCGGGGAACAC | ||
CGGTCTATCCCCGCTGGCGCGGGGAACAC | ||
CGCTTTATCCCCGCTGGCGCGGGGAACAC | ||
Internal and 3′ terminal repeat | CAGTTTATCCCCGCTGGCGCGGGGAACAC | |
Bases different from those of each consensus sequence are underlined |