# Comment on Tanmoy et al. CRISPR-Cas Diversity in Clinical Salmonella enterica Serovar Typhi Isolates from South Asian Countries. Genes 2020, 11, 1365

Laetitia Fabre, Elisabeth Njamkepo, François-Xavier Weill

► To cite this version:

Laetitia Fabre, Elisabeth Njamkepo, François-Xavier Weill. Comment on Tanmoy et al. CRISPR-Cas Diversity in Clinical Salmonella enterica Serovar Typhi Isolates from South Asian Countries. Genes 2020, 11, 1365. Genes, 2021, Special Issue Omics Research of Pathogenic Microorganisms, 12 (8), pp.1142. 10.3390/genes12081142. pasteur-03705909

HAL Id: pasteur-03705909

https://pasteur.hal.science/pasteur-03705909

Submitted on 27 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Comment

# Comment on Tanmoy et al. CRISPR-Cas Diversity in Clinical *Salmonella enterica* Serovar Typhi Isolates from South Asian Countries. *Genes* 2020, *11*, 1365

Laetitia Fabre, Elisabeth Njamkepo and François-Xavier Weill * ID

Institut Pasteur, Unité des Bactéries Pathogènes Entériques, 75015 Paris, France; laetitia.fabre@pasteur.fr (L.F.);
elisabeth.njamkepo-nguemkam@pasteur.fr (E.N.)

* Correspondence: francois-xavier.weill@pasteur.fr

Tanmoy et al. [1] report new findings relating to CRISPR locus organization and composition in *Salmonella enterica* serovar Typhi (hereafter referred to as S. Typhi). They reported that S. Typhi isolates can carry up to five different CRISPR loci and about 19% of the tested genomes had three or more CRISPR loci, whereas previous studies reported only two loci [2,3], suggesting that these studies were incomplete due to the use of too small a set of S. Typhi genomes.

As first described by Jansen et al. [4], clustered regularly interspaced short palindromic repeats (CRISPRs) are a family of repeated DNA sequences present in prokaryotes, and they are characterized by 24–47 bp DNA direct repeats (DRs), separated by variable 21–72 bp sequences called “spacers” [5,6]. An A-T-rich “leader sequence” and *cas* (CRISPR-associated sequence) genes are often identified adjacent to the CRISPR locus. Based on the 16 to 39 assembled genomes from *Salmonella* available at the time, including two from S. Typhi (CT18 and Ty2), the *Salmonella* CRISPR region was characterized by our group and others as two loci, CRISPR1 and CRISPR2, separated by less than 20 kb [2,7], with type I-E *cas* genes in the interval between the CRISPR loci [8] (Figure 1).

We also investigated the polymorphism of CRISPR1 and CRISPR2 by performing a PCR (polymerase chain reaction)–Sanger sequencing analysis on 744 *Salmonella* reference strains and isolates of 130 serovars (including 18 clinical S. Typhi isolates). Over 3800 unique spacer sequences were identified, stored and can be queried online at: https://galaxy.pasteur.fr/root?tool_id=toolshed.pasteur.fr/repos/khillion/salmonella_crisprTyping/salmonella_crisprTyping/1.0.0 (accessed on 5 July 2021) [2]. Mean spacer length was 32 bp (29–74 bp), and a 29 bp (26–30 bp) DR consensus sequence was identified; however, single nucleotide polymorphism (SNP) variants were also observed. The strong correlation between spacer content and serovar/multilocus sequence type [2] led to a patent describing a new *Salmonella*-subtyping method [9]. In the two S. Typhi genomes (Ty2 and CT18) and 18 S. Typhi isolates from diverse genotypes (13 different haplotypes [10]), geographic origins (9 countries) and time periods (1918–2006) studied, six CRISPR1 spacer sequences (Typhi1 to Typhi6), one CRISPR2 spacer sequence (EntB0var1) and seven different combined CRISPR1/CRISPR2 profiles were identified [9]. No CRISPR1 spacer was common to all 20 S. Typhi genomes and isolates studied. However, a spacer, EntB0var1, was found in all CRISPR2 sequences. This spacer sequence was then used to develop a S. Typhi serovar-specific PCR assay, which was validated on 188 S. Typhi reference strains and isolates of diverse genotypes (65 different haplotypes [10]), geographic origins (40 countries) and time periods (1918–2009). This CRISPR2 target was then used in a multiplex PCR for the detection of S. Typhi in blood samples from Bangladesh [11].

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Fabre, L.; Njamkepo, E.; Weill, F.-X. Comment on Tanmoy et al. CRISPR-Cas Diversity in Clinical *Salmonella enterica* Serovar Typhi Isolates from South Asian Countries. *Genes* 2020, *11*, 1365. *Genes* 2021, *12*, 1142. https://doi.org/10.3390/genes12081142

Academic Editors: J. Peter, W. Young and Thierry Wirth

Received: 23 March 2021

Accepted: 6 July 2021

Published: 28 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The CRISPR/Cas system structure of *S. enterica* serovar Typhi. The structure shown is that for the representative *S. enterica* serovar Typhi strain Ty2 (GenBank accession no. AE014613.1). Two CRISPR loci (CRISPR1 and CRISPR2) are present. The CRISPR-associated (*cas*) genes *cas2*, *cas1*, *cse3*, *cas5e*, *cse4*, *cse2*, *cse1* and *cas3* genes of the I-E type are located between the CRISPR loci. Diamonds represent direct repeats (DRs), with colored rectangles indicating spacers. The primers used to extract the two CRISPR loci from genomic sequences are shown as red horizontal arrows. The coordinates of the region (based on AE014613.1) are indicated.</figcaption></figure>

For comparison of the results reported by Tanmoy et al. [1] with those from our previous studies [2,3], the 1059 genomes described by the authors were downloaded from EBI-ENA (https://www.ebi.ac.uk/ena/browser/home, accessed on 24 November 2020) and assembled with SPAdes [12], according to the authors' parameters. The metrics of the assemblies (N50, genome size and N contigs) revealed evidence of the contamination of some genomes (ERR2663487, ERR2663542, ERR2663589, ERR2663887 and ERR2663969) with other *Salmonella* serovars (Enteritidis, Paratyphi A and Worthington), which was confirmed by molecular serotyping and/or multilocus sequence typing (Supplementary Materials Table S1, "Comment" column). The genotypi program [13] was used to check the genotypes of the downloaded genomes. Intriguingly, for more than 500 of the 1059 genomes, discrepancies were observed between the results presented by the authors and our analysis (Supplementary Materials Table S2). The discrepancies related to genome sequences from a previous publication by the authors [14]. We found that the strain name/accession code pairs shown in the authors' Dataset S1 [14] did not match those available from EBI-ENA (Supplementary Materials Table S3). We suspect that a single lane shift occurred between the accession codes and associated strain names after strain 311189_226186 (ERR2663487), either during construction of the spreadsheet or after its submission to EBI-ENA. We therefore used only the EBI-ENA accession codes as identifiers for the genomes in our review. This major issue also precluded the use of metadata, other than the country of origin (Bangladesh) and of the data published in the authors' papers [1,14]. We analyzed the CRISPR region with both CRISPRCasFinder [15] and a script developed in-house and based on *in silico* PCR for identification of the sequences of CRISPR1, CRISPR2 and spacer/DR sequences with the SalmCRISPRtyping script [2] (Supplementary Materials Table S1).

All the *S. Typhi* genomes belonged to the type I-E CRISPR-Cas system. No new CRISPR spacer or DR sequences were identified. Eleven combined CRISPR1/CRISPR2 profiles were identified among the 1059 genomes tested (Table 1 and Supplementary Materials Table S1). The four new combined CRISPR1/CRISPR2 profiles identified since the previous analysis of 20 *S. Typhi* genomes [2,3] were due entirely to variations of the number of DR64 at the 3'-end of CRISPR1. No clear association between the combined CRISPR1/CRISPR2 profiles and genotype was observed (Table 1).
<table><caption>**Table 1.** CRISPR arrays, CRISPR profiles and genotypes of the 1059 *S. enterica* serovar Typhi genomes described by Tanmoy et al. [1].</caption><thead><tr><th>CRISPR1</th><th>CRISPR2</th><th>Combined CRISPR Profile</th><th>No. of Strains</th><th>Genotypes (n)</th></tr></thead><tbody><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P1</td><td>905</td><td><b>4.3.1.1 (294), 4.3.1.2 (212), 4.3.1.3 (109), 4.3.1.1.P1 (88), 4.3.1 (50), 3.3.2 (47), 3.3.2.Bd2 (24), 2.0.0 (21), 2.2.0 (11), 3.3.0 (9), 2.0.1 (7), 2.1.7 (7), 3.3.1 (6), 2.3.3 (3), 3.0.0 (3), 3.0.2 (3), ND (3), 1.2.1 (2), 2.2.2 (2), 2.2.4 (1), 2.4.0 (1), 3.1.0 (1), 3.1.2 (1), 4.1.0 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P2</td><td>76</td><td><b>3.2.2 (73), 3.2.1 (2), 3.0.0 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR</td><td>DR27-EntB0var1-DR</td><td>P3</td><td>40</td><td><b>3.3.2.Bd1 (21), 2.5.0 (7), 2.2.0 (4), 3.3.2 (3), 2.2.2 (2), 3.0.1 (1), 3.3.0 (1), ND (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR64-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P4</td><td>19</td><td><b>2.3.3 (17), 2.2.1 (1), 2.3.4 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P5</td><td>5</td><td><b>4.3.1.1 (4), 4.1.0 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR</td><td>DR27-EntB0var1-DR</td><td>P6</td><td>3</td><td><b>3.3.1 (3)</b></td></tr><tr><td>DR1A1-Typhi1-DR</td><td>DR27-EntB0var1-DR</td><td>P7</td><td>2</td><td><b>4.3.1.1 (1), 2.0.0 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi6-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P8</td><td>2</td><td><b>4.1.0 (2)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi5-DR15-Typhi6-DR</td><td>DR27-EntB0var1-DR</td><td>P9</td><td>2</td><td><b>3.0.1 (2)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR15-Typhi6-DR64-DR64-DR64-DR64-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P10</td><td>1</td><td><b>2.3.4 (1)</b></td></tr><tr><td>DR1A1-Typhi1-DR-Typhi2-DR-Typhi3-DR-Typhi4-DR-Typhi5-DR-Typhi6-DR64-DR</td><td>DR27-EntB0var1-DR</td><td>P11</td><td>1</td><td><b>4.3.1.1 (1)</b></td></tr></tbody></table>

As the authors surprisingly found no correspondence between the CRISPR spacer sequences they identified and those from our previous studies (except for Ts32v), a further analysis of the sequences listed in their Dataset S3 [1] was conducted. As illustrated in Supplementary Materials Table S4, six of their most frequent “Group-A” spacers (Ts32h, Ts32c, Ts32l, Ts32e, Ts32i and Ts32g) corresponded to our previously described CRISPR1 spacers (Typhi1, Typhi2, Typhi3, Typhi4, Typhi5 and Typhi6, respectively). Another six of their “Group-A” spacers (Ts34f, Ts36a, Ts34a, Ts34b, Ts34c and Ts34e) were actually the CRISPR1 spacers described above, but with the inclusion of two to four nucleotides from the upstream or downstream DR sequences (Supplementary Materials Table S4). Seventeen of their identified S. Typhi “Group-A” spacers (Ts32t, Ts32u, Ts32q, Ts32s, Ts32p, Ts32k, Ts32a, Ts32d, Ts32v, Ts32f, Ts32n, Ts32w, Ts33a, Ts32b, Ts32r, Ts32o and Ts32m) actually correspond to spacer sequences from S. Enteritidis (Supplementary Materials
Table S4). Eight of these spacers (Ts32t, Ts32u, Ts32q, Ts32s, Ts32p, Ts32k, Ts32a and Ts32d) were identified in only one genome (ERR2663969, corresponding to ERR2663968 in Reference [1]), which was a mix-up between S. Typhi and S. Enteritidis (Supplementary Materials Table S1). Surprisingly, S. Enteritidis spacers Ts32f, Ts32n, Ts32w, Ts33a, Ts32b, Ts32r, Ts32o and Ts32m could not be identified in any of the 1059 genomes, including ERR2663969. Regarding the DR sequences, the “Group-A” Td29a corresponded to our DR, the most frequently found direct repeat in S. Typhi (Supplementary Materials Table S4). We confirmed the presence of DR/Td29a in all 1059 S. Typhi genomes, whereas the authors did not find it in 5/1059. The frequent DR variants DR15 (SNP variant) and DR64 (26 bp long variant) we previously found in S. Typhi CRISPR1 were not identified by Tanmoy et al. [1]. Two of the “Group-B” DR sequences (see below), Td39a and Td29b, were the exact reverse complemented sequences of another two DRs, Td39b and Td29c, respectively [1].

We took a closer look at the different CRISPR loci to understand the discrepancies between the authors' study and this review. First, our CRISPR2 locus was not identified by the authors in any of the 1059 genomes studied. This is not particularly surprising, as we had already reported in 2012 that CRISPRCasFinder software [15] was unable to detect the short CRISPR2 locus of S. Typhi strains Ty2 and CT18, which have a unique spacer (EntB0var1) between two DRs (DR27 and DR), one of which is degenerate (identity of 20/29 bp). Surprisingly, this undetected CRISPR2 target was previously used by the first and senior authors of the Tanmoy et al. [1] article, for the detection of S. Typhi by PCR in clinical samples [11]. The number of loci we identified (Supplementary Materials Table S1) with CRISPRCasFinder software (including loci with a low evidence score) did not really match the number of loci reported by the authors. The nucleotide sequences of each locus defined in the authors' Table 4 were then reconstructed according to the DR-spacer sequences described in their Dataset S3 [1], and the blastn algorithm of blast+ v.2.6.0 was then used to check for these sequences. For “Group-A” loci, patterns a1 and a3 (described in their Table 4) corresponded to CRISPR1 and CRISPR2 of S. Enteritidis, respectively. Therefore, only patterns a2, and a4 to a7 consist of S. Typhi CRISPR1 spacer sequences. One of these patterns, a6, found in genomes ERR2663783 and ERR2663776, is actually identical to the 421 bp length variant of a2 and should therefore be withdrawn. In comparison, our analysis identified 11 alleles for CRISPR1 (Table 1). Some of our profiles containing a limited number of CRISPR1 spacers, such as P1 (one spacer) and P6 (two spacers), were not identified by the authors (Supplementary Materials Table S5). Our profiles P5 and P8 correspond to their patterns a4 and a7, respectively. For the other profiles, there was no good correlation with their patterns. In particular, their prevalent pattern, a2, could be broken down into eight different profiles according to our analysis, due to the presence of DR64, sometimes repeated, in various profiles (Table 1). DR64 was not identified by the authors.

Regarding the new CRISPR loci identified by the authors (corresponding to “Group-B” loci), we found that the b4, and b21/b22 patterns corresponded to the normal CRISPR1 array found downstream from the *iap* gene [2] (“group-A” locus according to Reference [1]) and not to new CRISPR loci. Pattern b21, in particular, corresponds to our profile P6 (see above). The other “Group-B” loci were scored as “low evidence” by CRISPRCasFinder (levels 0 and 1 versus level 4 for CRISPR1), and no *cas* genes were detected in their vicinity. Furthermore, most of these new CRISPR loci consist of a minimal array (DR-spacer-DR). Investigation of the most frequent “Group-B” patterns, namely b1, b9 and b10–b13 (Supplementary Materials Figure S1) identified by the authors revealed that these “Group-A” loci consisted of repeated sequences, some being genuine variable number of tandem repeats (VNTRs). Hence, in the region defining “Group-B” loci patterns b10–b13, a large sequence of 93 bp (fusion of Ts54a/b and Td39a), repeated up to six times, depending on the genome, was misinterpreted as a CRISPR locus. Consequently, single loci consisting of limited numbers of repeated sequences with no nearby CRISPR-*cas* machinery, such as the “Group-B” loci described by the authors, should not be considered to be CRISPR loci.
This analysis confirms our previous results [2,3,9] and those of other groups [7,8] by showing that the genetically homogeneous *S. Typhi* population contains a single CRISPR-cas system (type I-E), with two adjacent CRISPR loci. Both CRISPR loci contain a limited number of spacers (1–6 in CRISPR1 and one in CRISPR2), as observed in other host-adapted *Salmonella* serovars with altered *cas* genes [2].

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/genes12081142/s1. Figure S1: Analysis of various “Group-B” CRISPR loci patterns described by Tanmoy et al. [1]. Table S1: Main characteristics of the 1059 *S. enterica* serovar Typhi genomes studied. Table S2: Discrepancies between genotype data from the study by Tanmoy et al. [1] and our study. Table S3: Strains described in Tanmoy et al. [14] and their different accession codes. Table S4: Correspondence between the spacer and DR sequences described by Tanmoy et al. [1] and those from Fabre et al. [2], and prevalence of these sequences among the 1059 *S. enterica* serovar Typhi genomes studied. Table S5: Correspondence between the “Group-A” CRISPR locus patterns from Tanmoy et al. [1] and our CRISPR profiles.

**Author Contributions:** Conceptualization, L.F., E.N. and F.-X.W.; formal analysis, L.F., E.N.; investigation, F.-X.W.; data curation, L.F.; writing—original draft preparation, L.F., F.-X.W.; writing—review and editing, E.N.; supervision, F.-X.W.; funding acquisition, F.-X.W. All authors have read and agreed to the published version of the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Tanmoy, A.M.; Saha, C.; Sajib, M.S.I.; Saha, S.; Komurian-Pradel, F.; van Belkum, A.; Louwen, R.; Saha, S.K.; Endtz, H.P. CRISPR-Cas Diversity in clinical *Salmonella enterica* serovar Typhi isolates from South Asian countries. *Genes* **2020**, *11*, 1365. [CrossRef] [PubMed]

2. Fabre, L.; Zhang, J.; Guigon, G.; Le Hello, S.; Guibert, V.; Accou-Demartin, M.; de Romans, S.; Lim, C.; Roux, C.; Passet, V.; et al. CRISPR typing and subtyping for improved laboratory surveillance of *Salmonella* infections. *PLoS ONE* **2012**, *7*, e36995. [CrossRef]

3. Fabre, L.; Le Hello, S.; Roux, C.; Issenhuth-Jeanjean, S.; Weill, F.X. CRISPR is an optimal target for the design of specific PCR assays for *Salmonella enterica* serotypes Typhi and Paratyphi A. *PLoS Negl. Trop. Dis.* **2014**, *8*, e2671. [CrossRef] [PubMed]

4. Jansen, R.; Embden, J.D.; Gaastra, W.; Schouls, L.M. Identification of genes that are associated with DNA repeats in prokaryotes. *Mol. Microbiol.* **2002**, *43*, 1565–1575. [CrossRef] [PubMed]

5. Grissa, I.; Vergnaud, G.; Pourcel, C. The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats. *BMC Bioinform.* **2007**, *8*, 172. [CrossRef] [PubMed]

6. Horvath, P.; Barrangou, R. CRISPR/Cas, the Immune System of Bacteria and Archaea. *Science* **2010**, *327*, 167–170. [CrossRef] [PubMed]

7. Touchon, M.; Rocha, E. The Small, Slow and Specialized CRISPR and Anti-CRISPR of Escherichia and Salmonella. *PLoS ONE* **2010**, *5*, e11126. [CrossRef] [PubMed]

8. Makarova, K.S.; Haft, D.H.; Barrangou, R.; Brouns, S.J.J.; Charpentier, E.; Horvath, P.; Moineau, S.; Mojica, F.J.M.; Wolf, Y.I.; Yakunin, A.F.; et al. Evolution and classification of the CRISPR-Cas systems. *Nat. Rev. Microbiol. Genet.* **2011**, *9*, 467–477. [CrossRef] [PubMed]

9. Weill, F.X.; Fabre-Berland, L.; Guibert, V.; Diancourt, L.; Brisse, S. Molecular Typing and Subtyping of Salmonella by Identification of the Variable Nucleotide Sequences of the CRISPR Loci. French Patent FR07/09188, 28 December 2007.

10. Roumagnac, P.; Weill, F.-X.; Dolecek, C.; Baker, S.; Brisse, S.; Chinh, N.T.; Le, T.A.H.; Acosta, C.J.; Farrar, J.; Dougan, G.; et al. Evolutionary History of *Salmonella* Typhi. *Science* **2006**, *314*, 1301–1304. [CrossRef] [PubMed]

11. Pouzol, S.; Tanmoy, A.M.; Ahmed, D.; Khanam, F.; Brooks, W.A.; Bhuyan, G.S.; Fabre, L.; Bryant, J.E.; Gustin, M.-P.; Vanhems, P.; et al. Clinical Evaluation of a Multiplex PCR for the Detection of *Salmonella enterica* Serovars Typhi and Paratyphi A from Blood Specimens in a High-Endemic Setting. *Am. J. Trop. Med. Hyg.* **2019**, *101*, 513–520. [CrossRef] [PubMed]

12. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. *J. Comput. Biol.* **2012**, *19*, 455–477. [CrossRef] [PubMed]

13. Wong, V.K.; Baker, S.; Connor, T.R.; Pickard, D.; Page, A.J.; Dave, J.; Murphy, N.; Holliman, R.; Sefton, A.; Millar, M.; et al. An extended genotyping framework for *Salmonella enterica* serovar Typhi, the cause of human typhoid. *Nat. Commun.* **2016**, *7*, 12827. [CrossRef] [PubMed]
14. Tanmoy, A.M.; Westeel, E.; De Bruyne, K.; Goris, J.; Rajoharison, A.; Sajib, M.S.I.; Van Belkum, A.; Saha, S.K.; Pradel, F.; Endtz, H.P. *Salmonella enterica* Serovar Typhi in Bangladesh: Exploration of Genomic Diversity and Antimicrobial Resistance. *mBio* **2018**, *9*, e02112-18. [CrossRef] [PubMed]

15. Couvin, D.; Bernheim, A.; Toffano-Nioche, C.; Touchon, M.; Michalik, J.; Néron, B.; Rocha, E.P.C.; Vergnaud, G.; Gautheret, D.; Pourcel, C. CRISPRCasFinder, an update of CRISFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. *Nucleic Acids Res.* **2018**, *46*, W246–W251. [CrossRef] [PubMed]