Cronobacter sakazakii, Cronobacter malonaticus, and Cronobacter dublinensis Genotyping Based on CRISPR Locus Diversity

Cronobacter strains harboring CRISPR-Cas systems are important foodborne pathogens that cause serious neonatal infections. CRISPR typing is a new molecular subtyping method to track the sources of pathogenic bacterial outbreaks and shows a promise in typing Cronobacter, however, this molecular typing procedure using routine PCR method has not been established. Therefore, the purpose of this study was to establish such methodology, 257 isolates of Cronobacter sakazakii, C. malonaticus, and C. dublinensis were used to verify the feasibility of the method. Results showed that 161 C. sakazakii strains could be divided into 129 CRISPR types (CTs), among which CT15 (n = 7) was the most prevalent CT followed by CT6 (n = 4). Further, 65 C. malonaticus strains were divided into 42 CTs and CT23 (n = 8) was the most prevalent followed by CT2, CT3, and CT13 (n = 4). Finally, 31 C. dublinensis strains belonged to 31 CTs. There was also a relationship among CT, sequence type (ST), food types, and serotype. Compared to multi-locus sequence typing (MLST), this new molecular method has greater power to distinguish similar strains and had better accordance with whole genome sequence typing (WGST). More importantly, some lineages were found to harbor conserved ancestral spacers ahead of their divergent specific spacer sequences; this can be exploited to infer the divergent evolution of Cronobacter and provide phylogenetic information reflecting common origins. Compared to WGST, CRISPR typing method is simpler and more affordable, it could be used to identify sources of Cronobacter food-borne outbreaks, from clinical cases to food sources and the production sites.

Introduction

The Cronobacter (formerly Enterobacter sakazakii) genus including C. sakazakii, C. malonaticus, C. dublinensis, C. turicensis, C. universalis, C. muytjensii, and C. condimenti comprises opportunistic foodborne pathogens that can cause rare but life-threatening diseases in neonates and immune-compromised infants, including meningitis, necrotizing enterocolitis, and septicemia (Iversen et al., 2008; Kucerova et al., 2011; Joseph et al., 2012a; Zeng et al., 2018a). Latest study reported an acute gastroenteritis outbreak caused by C. sakazakii in a senior high school of China (Yong et al., 2018). Moreover, this genus has been isolated from the environment, food, and clinical sources (Ueda, 2017; Yong et al., 2018; Zeng et al., 2018a; Li et al., 2019). Further, C. sakazakii, C. malonaticus, and C. dublinensis are three prevalent species and some reports have indicated that the principal sources of these organisms might be soil, water, and vegetables (Ueda, 2017; Zeng et al., 2018b). However, the epidemiology and reservoirs of Cronobacter spp. is still unsure (Holy and Forsythe, 2014).

Some molecular subtyping methods have been developed to study the epidemiology of pathogenic bacteria, including pulsed field gel electrophoresis (PFGE) and multi-locus sequence typing (MLST), but both still have some disadvantages (Ogrodzki and Forsythe, 2017). PFGE is limited for a portion of Cronobacter strains that cannot be typed due to intrinsic DNase activity; moreover, it does not provide the phylogenetic relationship between strains. MLST has been established for Cronobacter genus based on seven housekeeping genes (Joseph et al., 2012b). A curated open access MLST database has been established for the genus with more than 2200 strains and associated metadata¹. This database has enabled the recognition of certain Cronobacter clonal lineages within the genus as pathogenic variants, whereas others are primarily commensal organisms associated with the environment. The discrimination power of MLST is weaker than that of whole genome sequence typing (WGST), and this method lacks information about historical ancestors (Forsythe et al., 2014). WGST is a new method for subtyping bacteria, but its high costs still limit its application (Forsythe et al., 2014; Deng et al., 2015).

CRISPR-Cas system is an adaptive immune system for bacteria, providing bacteria with sequence-specific, acquired defense against phages and plasmids (Barrangou, 2013; Westra et al., 2014). The evolution of CRISPR-Cas has led to the discovery of a diverse set of CRISPR-Cas systems, which can be then classified into distinct classes, types, and subtypes, combined with the analysis of signature protein families and features of cas loci architectures that unambiguously partition most CRISPR–Cas loci (Makarova et al., 2015; Shmakov et al., 2017). The activity of a CRISPR locus occurs in three stages as follows: adaptation through the incorporation of new spacers into the existing repeat-spacer array; expression of the repeat-spacer array and the consequent processing of that array into CRISPR RNAs (crRNAs); interference during which invasive target sequences are recognized and destroyed by the crRNA-effector complex (Barrangou et al., 2007). As new spacers are added to one end of the CRISPR array, polarity exists; specifically, spacers at the leader distal end are more ancient and are often shared among bacterial common ancestors. The acquisition, loss, and duplication of spacers have made CRISPR arrays be the fastest evolving loci in bacteria (Paez-Espino et al., 2013; Shariat and Dudley, 2014).

The first application of CRISPR loci in bacterial genotyping was spacer-oligonucleotide typing (or “spoligotyping”) of Mycobacterium tuberculosis strains (Groenen et al., 1993; Streicher et al., 2007). Its principle is PCR amplification of the CRISPR array with labeled primers that recognize the directed repeat sequences, then hybridization of the PCR products to a membrane containing probes bearing spacer DNA sequences (Streicher et al., 2007). The “next-generation” microbead-spoligotyping approach was an assay termed CRISPOL (for “CRISPR polymorphism”) applied to Salmonella (Fabre et al., 2012). The first application of sequence-based CRISPR typing was group A Streptococcus (GAS) M1 serotype (Hoe et al., 1999). Considered the temporal organization of spacers, the sequencing of CRISPR arrays has been a extremely useful tool to genotype bacteria like Yersinia species, E. coli, and Salmonella enterica (Cui et al., 2008; Fricke et al., 2011; Yin et al., 2013; Li et al., 2014; Bugarel et al., 2018), and it has also been used to investigate bacterial diversity based on metagenomic data (Berg Miller et al., 2012; Sun et al., 2016). Recently, some useful tools to extract spacers and visualize the spacer content with color schemes were developed (Biswas et al., 2016; Couvin et al., 2018; Dion et al., 2018; Nethery and Barrangou, 2019).

In previous studies, six CRISPR arrays were detected in conserved regions of the Cronobacter genomes; among these, CRISPR1 and CRISPR2 neighbor the I-E type of the “complete” cas gene cluster, whereas CRISPR3 and CRISPR6 integrate with the I-F type of the “complete” cas gene cluster comprising subtype I-E and I-F CRISPR-Cas systems, respectively. Two CRISPR-Cas systems (Subtype I-E and I-F) were found only in C. sakazakii, C. malonaticus, and C. dublinensis isolates, specifically. Unlike subtype I-E, which was commonly detected among Cronobacter strains, subtype I-F was found to be significantly more prevalent in the plant-associated species C. dublinensis than in the human virulence-related species C. sakazakii and C. malonaticus. However, C. condimenti lacked intact CRISPR-Cas system (Zeng et al., 2017, 2018b). At the same time, significantly higher CRISPR activity was also observed in the plant-associated species C. dublinensis than in the virulence-related species C. sakazakii and C. malonaticus (Zeng et al., 2018b). Similar CRISPR array spacers have been rarely detected among species, indicating intensive changes through adaptive acquisition and loss. Thus, differentiated CRISPR activity appears to be the product of environmental selective pressure and might contribute to the bidirectional divergence and speciation of Cronobacter (Zeng et al., 2018b).

CRISPR arrays will be a promising typing method compared with MLST (Ogrodzki and Forsythe, 2016, 2017; Zeng et al., 2018b). However, the identification of CRISPR arrays in Cronobacter was based on whole genome sequences by next generation sequencing, which is associated with high cost. It is necessary to establish this new molecular subtyping method using routine PCR and define the nomenclature system. Among six types of CRISPR arrays, CRISPR1 and CRISPR2 were found in almost all Cronobacter strains, whereas CRISPR3 and CRISPR6 were also found to be preserved in many Cronobacter strains, and all of them have a high diversity among different isolates (Zeng et al., 2018b). However, CRISPR4 and CRISPR5 were not suitable for genotyping, as they were only found in a few C. sakazakii isolates and lacked spacer diversity (Zeng et al., 2017). The diversity of four CRISPR arrays in isolates of this genus could therefore provide a powerful tool to track the origin of genetically similar strains within an outbreak. In this study, we established a CRISPR-based subtyping method for C. sakazakii, C. malonaticus, and C. dublinensis using routine PCR and examined the relationship between CRISPR profiles and other genetic factors.

Materials and Methods

Bacterial Isolates

A total of 257 Cronobacter isolates used in this study were collected from four types of food (powdered infant milk, ready-to-eat food, vegetables, and edible mushroom) in China, including 161 C. sakazakii, 65 C. malonaticus, and 31 C. dublinensis strains. All strains belonged to the large-scale and systematic investigation on the prevalence of Cronobacter spp. in food in China, and detailed information about these strains, including O serotypes, STs, and antibiotic-resistance profiles, are provided in Supplementary Data Sheets 1–3 (Zeng et al., 2017, 2018b; Ling et al., 2018; Li et al., 2019).

CRISPR PCR Amplification and Sequencing

The primers location for the amplification CRISPR1, CRISPR2, CRISPR3, and CRISPR6 are shown in Figure 1, and are in accordance with the genomic sequences encoding CRISPR-Cas systems in Cronobacter reported previously (Zeng et al., 2018b). The sequences of primers for the amplification and sequencing of CRISPR1, CRISPR2, CRISPR3, and CRISPR6 loci are listed in Table 1. PCR reaction was performed using a 50-μL volume, which contained 0.5 μL of PrimeSTAR^®HS DNA Polymerase (2.5 U/μL; Takara, Dalian, Japan), 4 μL of 2.5 mM dNTPs, 0.5 μL of each 10 mM primer, 10 μL of 5 × PrimeSTAR Buffer, and 1 or 2 μL of bacterial DNA template for CRISPR1 and CRISPR2 or CRISPR3 and CRISPR6, respectively; the remaining volume consisted of sterile water. The PCR conditions were as follows: initial denaturation at 98°C for 1 min; 30 cycles of 98°C for 10 s, 58°C for CRISPR1 and 57°C for CRISPR2, CRISPR3, and CRISPR6 for 5 s, and 72°C for 4 min; a final extension at 72°C for 5 min. After identification by electrophoresis, the PCR products were subjected to DNA sequencing directly (Beijing Genomics Institute, Guangzhou, China). All PCR products were sequenced using amplification primers in both the forward and reverse directions to obtain a double-stranded sequence. PCR products larger than 2-kb sometimes required the design of additional primers based on Sanger sequencing results by amplification primers to obtain the complete sequences.

FIGURE 1

Figure 1. Outline of the new molecular typing method based on four CRISPR arrays of Cronobacter. The locations of PCR primers used to amplify CRISPR loci are shown. Compared to that in C. sakazakii and C. malonaticus, there was a 1-kb plus nucleotide sequence region including one hypothetical protein between the hypothetical protein used for the design of primer E-1F and CRISPR1 in C. dublinensis isolates. The orientation of CRISPR arrays and extraction of spacers were completed by CRISPRDetect. The specific CRISPR type was determined by a combination of sequenced incorporated spacers in CRISPR arrays.

TABLE 1

Table 1. Primers for CRISPR1, CRISPR2, CRISPR3, and CRISPR6 loci.

CRISPR Typing and Cluster Analysis

The orientation of CRISPR spacers were determined by CRISPRDetect and the spacers were extracted using CRISPRCasFinder (Biswas et al., 2016; Couvin et al., 2018). A similarity search of the identified spacer sequences (84% similarity) and the establishment of a unique spacer library were performed as described previously (Zeng et al., 2017). A comparison of these unique spacers to previously studied elements within the ACLAME database (Leplae et al., 2010) was performed to identify potential targets. Every unique spacer among different CRISPR arrays of one species was assigned a single number beginning with 1 from the leader distal end, and lists of CRISPR spacer sequences for C. sakazakii, C. malonaticus, and C. dublinensis are provided as Supplementary Data Sheets 4–6, respectively. Then, every CRISPR array with multiple spacers was assigned a number as a spacer code. CRISPR typing was performed by combining CRISPR1, CRISPR2, CRISPR3, and CRISPR6 into one allele and displayed this as an arrangement of CRISPR spacers. The CRISPR type (CT) of each isolate was defined using a specific number to reflect its unique allelic type. The discrimination index (D) was calculated based on the Simpson’s index of diversity with the equation as previously defined (Hunter and Gaston, 1988). To depict the clustering of subtypes determined by CRISPR diversity, the binary distribution (presence as “1” or absence as “0”) of every spacer in each CRISPR locus was profiled for each strain. The binary distribution patterns of all strains were then combined and used to create a minimum spanning tree, developed utilizing BioNumerics version 7.6.3 (Applied Maths, Belgium). To explore the genetic relationships between CRISPR sequence variability and food type or serotype and CTs and food type, CTs and serotypes were displayed according to the results of the cluster analysis, respectively. Differences in CRISPR spacers comparing antibiotic-resistant and susceptible isolates were also examined. The spacer comparison and conversions to HEX color code were performed using CRISPRstudio software (Dion et al., 2018).

Core Genome Phylogenetic Analyses

Among 257 Cronobacter isolates, whole genome sequences of 117 isolates were established based on core genome analyses and reported in our previous study (Zeng et al., 2017, 2018b). Next, a core genome ML phylogenetic tree was generated based on 287,220 nucleotides from concatenated 563 single-copy core genes sequences using FastTree (Price et al., 2009). The display and annotation of phylogenetic trees were performed using iTOL (Letunic and Bork, 2016).

Results

CRISPR Types

Two hundred and fifty-seven isolates of C. sakazakii, C. malonaticus, and C. dublinensis were used to establish CRISPR subtyping method. The design of primers (sequences were listed in Table 1) and procedure of Cronobacter CRISPR typing method were shown in Figure 1. As shown in Figure 2, CRISPR1 and CRISPR2 loci were more conserved and active than others in all species; moreover, CRISPR2 had the largest average number of spacers. These results were in accordance with our previous study (Zeng et al., 2017, 2018b), and similar spacers were rarely detected among species, indicating intensive changes through adaptive acquisition and loss. In this study, the incidences of CRISPR1, CRISPR2, CRISPR3, and CRISPR6 in 161 C. sakazakii isolates were 99.4% (160/161), 99.4% (160/161), 57.1% (92/161), and 10.6% (17/161), respectively (Supplementary Data Sheet 1). Moreover, 1706 unique spacers were identified in C. sakazakii strains, and these were divided into 129 CTs; CT15 (n = 7) was the most prevalent followed by CT6 (n = 4). Regarding 65 C. malonaticus isolates (Supplementary Data Sheet 2), the incidences of CRISPR1, CRISPR2, CRISPR3, and CRISPR6 were 90.8% (59/65), 100% (65/65), 12.3% (8/65), and 0% (0/65), respectively. For this species, 487 unique spacers were identified in C. malonaticus strains, and they were divided into 42 CTs with CT23 (n = 8) being the most prevalent CT followed by CT2 (n = 4), CT3 (n = 4), and CT13 (n = 4). In C. dublinensis (Supplementary Data Sheet 3), the frequencies of CRISPR1, CRISPR2, CRISPR3, and CRISPR6 were 74.2% (23/31), 100% (31/31), 19.4% (6/31), and 4.9% (2/31), respectively. Further, 1361 unique spacers were identified in 31 C. dublinensis strains, and these belonged to 31 CTs. In these C. sakazakii, C. malonaticus, and C. dublinensis isolates, the discriminatory powers (a single numerical index of discrimination [D]) of CT were 0.9957, 0.9736, and 1.0000, respectively, indicating that there should be a 99.6, 97.4, and 100.0% probability that two unrelated isolates can be separated using the CT scheme. This method has comparable power to distinguish these species. The discriminatory powers of MLST for these species were 0.9669, 0.8986, and 0.9892, respectively, among these isolates. Thus, the CRISPR typing method showed better discriminatory power than MLST.

FIGURE 2

Figure 2. Number of spacers from four CRISPR arrays in C. sakazakii, C. malonaticus, and C. dublinensis strains.

Relationship Between CRISPR Sequence Variability and Food Type and Serotype and Antibiotic Resistance

Minimum spanning trees were generated using BioNumerics software to analyze the distribution of CTs among different types of food and their relationship with serotypes (Figure 3). C. sakazakii, C. malonaticus, and C. dublinensis strains isolated from vegetables showed higher CRISPR diversity than those from other types of food, and this was especially true for C. dublinensis. This was in accordance with previous studies showing a higher frequency and diversity of Cronobacter in vegetables compared to that in other types of food, supporting the contention that this species is plant-associated (Ueda, 2017; Ling et al., 2018; Silva et al., 2019). There was also a relationship between CT and serotype. As shown in Figure 2, when the maximum distance between nodes in the same partition was set to 10, CT6-, CT64-, CT15-, CT41-, and CT85-associated partitions were the five main partitions in C. sakazakii (Figure 3A). C. sakazakii serotype O2 was found among all strains of the CT85-associated partition and most strains of partition CT6; moreover, serotype O1 predominated the CT64- and CT48-associated partitions and most strains in the CT15-associated partition were serotype O4 (Figure 3B). For C. malonaticus, CT13-, CT23-, and CT2-associated partitions were the three major partitions (Figure 3C). All strains in CT23- and CT3-associated partitions were serotype O1, whereas serotype O2 was predominant in the CT13-associated partition (Figure 3D). Based on the limited number of isolates and high diversity of CRISPR sequences in C. dublinensis, there were no major partitions reported in this study, whereas O1 was the predominant serotype (Figure 3E). In accordance with previous studies, 96.9% (249/257) of isolates were resistant or intermediate to cephalothin, whereas most were susceptible to other antibiotics (Brandao et al., 2017; Ling et al., 2018). In total, there were three isolates resistant to two or more antibiotics in this study. Comparisons of CRISPR sequence variability between the resistant strains and other strains were also performed (Supplementary Data Sheets 1–3), there was no significant relationship between antibiotic resistance and CRISPR variability in C. sakazakii, C. malonaticus, and C. dublinensis.

FIGURE 3

Figure 3. Minimum spanning tree of CRISPR data from 161 C. sakazakii, 65 C. malonaticus, and 31 C. dublinensis isolates. Minimum spanning tree of C. sakazakii (A), C. malonaticus (C), and C. dublinensis (E) isolates with color corresponding to each type of food indicated in the legend on the right side of (A). Minimum spanning tree of C. sakazakii (B), C. malonaticus (D), and C. dublinensis (F) with color corresponding to each serotype indicated in the legend on the right side. Each circle represents one CRISPR type (CT), and the area of the circle corresponds to the number of isolates. The maximum distance between nodes in the same partition was set to 10.

Accordance Among CRISPR Typing, MLST, and WGST

A core genome ML phylogenetic tree based on whole genome sequences of 117 strains was generated to evaluate the consistency between CRISPR typing and WGST. As shown in Figure 3, CRISPR profiles were conserved among phylogenetically related strains and these had a close relationship with ST types. At the same time, the strains with different STs but belonging to the same clonal complex (CC) also had similar CRISPR profiles and belonged to the same partition. C. sakazakii CC4, C. sakazakii CC8, and C. malonaticus CC7 were major pathogenic CCs in previous studies, all the strains in these CCs formed distinct clusters in the phylogenetic tree, and belonged to C. sakazakii CT6-, C. sakazakii CT64-, and C. malonaticus CT13-associated partitions, respectively. Moreover, this approach was found to distinguish the same ST into smaller units (Figures 4, 5). For example, seven ST64 isolates formed a small lineage in this phylogenetic tree, and three ST64 strains within the CT89-associated partition were more closely related than other strains of different CTs. At the same time, the phylogenetic distance between other ST64 strains was also in accordance with the differences in CRISPR spacer composition (Supplementary Data Sheet 1). The same phenomenon was observed for the ST23 strain (Figures 4, 5A). In C. malonaticus CC7, C. malonaticus ST7 isolates typed as CT12, CT14, CT13, and CT15 were more closely phylogenetically related to C. malonaticus ST211 isolates typed as CT11 than other ST7 isolates themselves (Figures 4, 5B), implying better accordance between the CRISPR typing method and WGST. However, there were also few inconsistent results, for example, C. sakazakii ST4 isolate cro7 and ST267 isolate cro1511C1 were C. sakazakii CT2, but cro7 have more closely phylogenetic relationship with another C. sakazakii ST4 isolate 7G. Combined with all the results, the CRISPR typing method shows better discriminatory power than MLST and has a better accordance with WGST.

FIGURE 4

Figure 4. Phylogeny of 118 C. sakazakii, C malonaticus, and C. dublinensis strains inferred by whole genome sequences types (WGSTs). The STs and CRISPR types (CTs) of each isolate were listed on the right side and the CRISPR profiles of clonal complex 4 (CC4), CC8, CC7, ST148, ST60, and ST77 strains were also shown.

FIGURE 5

Figure 5. CRISPR spacer overview. Organization of spacer content of CRISPR alleles identified in 20 C. sakazakii isolates (A) and 33 C. malonaticus isolates (B). Repeats were not shown in this figure, and only spacers were displayed. Color schemes were provided at the spacer level to visualize differences among isolates based on the software CRISPRStudio. Spacers are shown in the order of predicted acquisition in the locus (right, ancestral spacers; left, newly acquired spacers).

The Phylogenetic Information Inferred by Sequence Diversity in CRISPR Arrays

In addition to the good discriminatory power of CRISPR arrays in distinguishing Cronobacter strains, the phylogeny information conserved in the iterative spacer acquisition process can be used to infer common ancestry. The CRISPR alleles of C. sakazakii ST23, ST264, ST148, and ST566 strains are shown in Figure 5A. These ST isolates belonged to different CCs and no apparent close phylogenetic relationship among these STs was observed in Figure 3. In contrast to the high diversity of CRISPR spacers among these strains at four CRISPR loci, it was interesting to note that all of these strains harbored some conserved ancestral spacers in CRISPR1. There were seven ancestral spacers conserved in CRISPR1, and some ST23, ST264, and ST148 strains preserved all of these ancestral spacers. Moreover, one additional spacer inserted in the fourth and fifth ancestral spacers was detected in all ST148 strains. Unlike the first seven conserved spacer sequences, the newly incorporated spacers showed lineage specificity. The ancestral spacers might be important proof of lineage divergence. As shown in Figure 5B, C. malonaticus CC7, ST211, ST11, and ST139 isolates preserved one ancestral sequence in CRISPR2. ST139 strains had three to four ancestral spacers that were common with CC7 isolates, which suggests a closer phylogenetic relationship between these lineages.

Discussion

The Cronobacter genus including seven species are opportunistic foodborne human pathogens that can cause rare but serious diseases in neonates and immune-compromised infants (Iversen et al., 2008; Kucerova et al., 2011; Joseph et al., 2012a). In our previous studies, C. sakazakii, C. malonaticus, and C. dublinensis were three prevalent species in food, however, we have never isolated C. universalis and C. condimenti strains (Xu et al., 2015; Ling et al., 2018; Li et al., 2019). In total, five C. turicensis and one C. muytjensii strains were isolated from vegetables and ready-to-eat foods (Xu et al., 2015; Ling et al., 2018) and we have successfully performed CRISPR arrays on these isolates using the same primers used for C. sakazakii. However, for the limited number of strains, whether these primers will be suitable for these species is unknown. Thus, we only constructed a CRISPR typing method for C. sakazakii, C. malonaticus, and C. dublinensis in this study.

In this study, CRISPR arrays were detected in all Cronobacter isolates; moreover, 1706, 487, and 1361 unique spacers were identified in 161 C. sakazakii, 65 C. malonaticus, and 31 C. dublinensis isolates. In accordance with a previous study (Zeng et al., 2018b), the number of CRISPR spacers in C. dublinensis isolates was greater than that in C. sakazakii and C. malonaticus. CRISPR1 and CRISPR2 were preserved in all three species and more active than other CRISPR loci; further, CRISPR3 was found in some strains of these species; however, CRISPR6 was only detected in some C. sakazakii and C. dublinensis strains (Supplementary Data Sheets 1–3). Whether there is a need to use four CRISPR loci for C. malonaticus CRISPR typing should be examined in the future using more isolates. Moreover, in these C. sakazakii, C. malonaticus, and C. dublinensis isolates, the discriminatory powers of the CRISPR typing method for all three species were comparable. According to our results, the CRISPR typing method shows better discriminatory power than MLST and has a better accordance with WGST. The largest outbreak of C. sakazakii occurred in a neonatal intensive care unit in France (1994), lasting over 3 months and claiming the lives of three neonates. A recent study used whole genome sequencing data of 26 isolates obtained from this outbreak to reveal relatedness (Masood et al., 2015). To examine the accuracy of CRISPR typing for the identification of pathogens in the Cronobacter outbreak, we downloaded these genome sequences and extracted CRISPR arrays for molecular typing. All C. sakazakii ST4, ST12, and ST13 strains belonged to CT2, CT50, and CT52, respectively. This was in accordance with the data obtained from the outbreak but had weaker discriminatory power compared to whole genome SNP analyses. In this study, 19 C. sakazakii ST4 strains isolated from several types of food in China were divided into 14 CTs including CT2, whereas four C. sakazakii ST13 isolates were divided into four CTs, but without CT52. Thus, the better discriminatory power of CRISPR typing could make it more useful than MLST to differentiate potential sources of Cronobacter outbreaks in the future.

Polarity exists as new spacers are always added to the proximal end of the CRISPR array; in addition, spacers at the leader distal end were found to be more ancient and were shared among phylogenetically related Cronobacter isolates. Spacer loss and gain make CRISPR elements the fastest evolving loci in Cronobacter, supporting previous speculation that CRISPR-Cas systems have an important impact on the evolution of this genus (Zeng et al., 2018b). CRISPR spacer variability in Cronobacter can divide an ST into smaller units and has better accordance with WGST than MLST. The CRISPR1 and CRISPR2 spacers in three species were more active, as shown in Figure 5, and some phylogenetically distant lineages were found to preserve some ancestral spacers at CRISPR1 or CRISPR2, respectively, although no similar spacers existed in other CRISPR loci. These ancestral spacers are important proof of lineage divergence; thus, CRISPR1 and CRISPR2 in Cronobacter can provide phylogenetic anchors reflecting common origins. Unfortunately, despite the extremely high variation in CRISPR spacer sequences of Cronobacter, many lineages had a unique CRISPR pattern, and no common ancestral spacers were found among these different clonal isolates. In summary, CRISPR diversity can be used to unfold a complete evolutionary story of strain divergence and relatedness, showing unique advantages compared to other genotyping methods.

The advantages of CRISPR-based genotyping methods have been demonstrated for some bacteria widely found in the food supply chain such as Streptococcus thermophilus (Horvath et al., 2008) and Lactobacillus buchneri (Briner and Barrangou, 2014). It can also be used for pathogenic strains like Escherichia coli (Yin et al., 2013; Barrangou and Dudley, 2016), Salmonella (Fabre et al., 2012; Li et al., 2014), Clostridium difficile (Andersen et al., 2016), and Mycobacterium tuberculosis (Streicher et al., 2007; Zhang et al., 2010). Finally, we also found a relationship between CT, ST, food types, and serotypes among Cronobacter isolates, and this phenomenon has also been found in other foodborne pathogens (Li et al., 2014, 2018; Bugarel et al., 2018).

Conclusion

In conclusion, we developed a CRISPR typing method for C. sakazakii, C. malonaticus, and C. dublinensis. Compared to MLST, this new molecular method has greater power to distinguish similar strains and had better accordance with WGST. Compared to WGST, CRISPR typing is simpler and more affordable, and it could be useful for the identification of sources of Cronobacter outbreaks, in addition to performing microbial risk assessment during food processing. More importantly, CRISPR diversity can be used to infer the divergent evolution of Cronobacter and provide phylogenetic anchors reflecting common origins. In the future, it would be meaningful to generate a comprehensive Cronobacter complex database of CRISPR spacers for the global application of CRISPR typing, and pool the results of different research groups to explore the epidemiology and reservoirs of Cronobacter spp.

Data Availability

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

HZ and QW conceived and designed the experiments. HZ, CL, WH, MC, TL, HW, and NL performed the experiments. HZ, JZ, SC, JW, and YD analyzed the data. HZ drafted the manuscript. QW supervised the project. All authors read and approved the final manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (31601571), the National Key R&D Program of China (2017YFC1601201), Local Innovative and Research Teams Project of Guangdong PEARL River Talents Program (2017BT01S174), the Natural Science Foundation of Guangdong Province (2016A030310315), Pearl River S&T Nova Program of Guangzhou (201806010062), and GDAS’ Special Project of Science and Technology Development (2017GDASCX-0201).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01989/full#supplementary-material

Footnotes

1. ^ https://pubmlst.org/cronobacter/

References

Andersen, J. M., Shoup, M., Robinson, C., Britton, R., Olsen, K. E., and Barrangou, R. (2016). CRISPR diversity and microevolution in clostridium difficile. Genome Biol. Evol. 8, 2841–2855. doi: 10.1093/gbe/evw203

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrangou, R. (2013). CRISPR-Cas systems and RNA-guided interference. Wiley Interdiscip. Rev. RNA 4, 267–278. doi: 10.1002/wrna.1159

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrangou, R., and Dudley, E. G. (2016). CRISPR-Based typing and next-generation tracking technologies. Annu. Rev. Food Sci. Technol. 7, 395–411. doi: 10.1146/annurev-food-022814-015729

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrangou, R., Fremaux, C., Deveau, H., Richards, M., Boyaval, P., Moineau, S., et al. (2007). CRISPR provides acquired resistance against viruses in prokaryotes. Science 315, 1709–1712. doi: 10.1126/science.1138140

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg Miller, M. E., Yeoman, C. J., Chia, N., Tringe, S. G., Angly, F. E., Edwards, R. A., et al. (2012). Phage-bacteria relationships and CRISPR elements revealed by a metagenomic survey of the rumen microbiome. Environ. Microbiol. 14, 207–227. doi: 10.1111/j.1462-2920.2011.02593.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Biswas, A., Staals, R. H., Morales, S. E., Fineran, P. C., and Brown, C. M. (2016). CRISPRDetect: a flexible algorithm to define CRISPR arrays. BMC Genomics 17:356. doi: 10.1186/s12864-016-2627-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandao, M. L., Umeda, N. S., Jackson, E., Forsythe, S. J., and de Filippis, I. (2017). Isolation, molecular and phenotypic characterization, and antibiotic susceptibility of Cronobacter spp. from Brazilian retail foods. Food Microbiol 63, 129–138. doi: 10.1016/j.fm.2016.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Briner, A. E., and Barrangou, R. (2014). Lactobacillus buchneri genotyping on the basis of clustered regularly interspaced short palindromic repeat (CRISPR) locus diversity. Appl. Environ. Microbiol. 80, 994–1001. doi: 10.1128/AEM.03015-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Bugarel, M., Bakker, H. D., Grout, J., Vignaud, M. L., Loneragan, G. H., Fach, P., et al. (2018). CRISPR-based assay for the molecular identification of highly prevalent Salmonella serotypes. Food Microbiol. 71, 8–16. doi: 10.1016/j.fm.2017.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Couvin, D., Bernheim, A., Toffano-Nioche, C., Touchon, M., Michalik, J., Neron, B., 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. doi: 10.1093/nar/gky425

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Y., Li, Y., Gorge, O., Platonov, M. E., Yan, Y., Guo, Z., et al. (2008). Insight into microevolution of yersinia pestis by clustered regularly interspaced short palindromic repeats. PLoS One 3:e2652. doi: 10.1371/journal.pone.0002652

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, X., Shariat, N., Driebe, E. M., Roe, C. C., Tolar, B., Trees, E., 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. doi: 10.1128/JCM.02332-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Dion, M. B., Labrie, S. J., Shah, S. A., and Moineau, S. (2018). CRISPRStudio: a user-friendly software for rapid CRISPR array visualization. Viruses 10:E602. doi: 10.3390/v10110602

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabre, L., Zhang, J., Guigon, G., Le Hello, S., Guibert, V., Accou-Demartin, M., et al. (2012). CRISPR typing and subtyping for improved laboratory surveillance of Salmonella infections. PLoS One 7:e36995. doi: 10.1371/journal.pone.0036995

PubMed Abstract | CrossRef Full Text | Google Scholar

Forsythe, S. J., Dickins, B., and Jolley, K. A. (2014). Cronobacter, the emergent bacterial pathogen Enterobacter sakazakii comes of age; MLST and whole genome sequence analysis. BMC Genomics 15:1121. doi: 10.1186/1471-2164-15-1121

PubMed Abstract | CrossRef Full Text | Google Scholar

Fricke, W. F., Mammel, M. K., McDermott, P. F., Tartera, C., White, D. G., Leclerc, J. E., et al. (2011). Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPR-mediated adaptive sublineage evolution. J. Bacteriol. 193, 3556–3568. doi: 10.1128/JB.00297-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Groenen, P. M., Bunschoten, A. E., van Soolingen, D., and van Embden, J. D. (1993). Nature of DNA polymorphism in the direct repeat cluster of Mycobacterium tuberculosis; application for strain differentiation by a novel typing method. Mol. Microbiol. 10, 1057–1065. doi: 10.1111/j.1365-2958.1993.tb00976.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoe, N., Nakashima, K., Grigsby, D., Pan, X., Dou, S. J., Naidich, S., et al. (1999). Rapid molecular genetic subtyping of serotype M1 group A Streptococcus strains. Emerg. Infect. Dis. 5, 254–263. doi: 10.3201/eid0502.990210

PubMed Abstract | CrossRef Full Text | Google Scholar

Holy, O., and Forsythe, S. (2014). Cronobacter spp. as emerging causes of healthcare-associated infection. J. Hosp. Infect. 86, 169–177. doi: 10.1016/j.jhin.2013.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Horvath, P., Romero, D. A., Coute-Monvoisin, A. C., Richards, M., Deveau, H., Moineau, S., et al. (2008). Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus. J. Bacteriol. 190, 1401–1412. doi: 10.1128/JB.01415-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunter, P. R., and Gaston, M. A. (1988). Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J. Clin. Microbiol. 26, 2465–2466.

PubMed Abstract | Google Scholar

Iversen, C., Mullane, N., McCardell, B., Tall, B. D., Lehner, A., Fanning, S., et al. (2008). Cronobacter gen. nov., a new genus to accommodate the biogroups of Enterobacter sakazakii, and proposal of Cronobacter sakazakii gen. nov., comb. nov., Cronobacter malonaticus sp. nov., Cronobacter turicensis sp. nov., Cronobacter muytjensii sp. nov., Cronobacter dublinensis sp. nov., Cronobacter genomospecies 1, and of three subspecies, Cronobacter dublinensis subsp. dublinensis subsp. nov., Cronobacter dublinensis subsp. lausannensis subsp. nov. and Cronobacter dublinensis subsp. lactaridi subsp. nov. Int. J. Syst. Evol. Microbiol. 58, 1442–1447. doi: 10.1099/ijs.0.65577-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Joseph, S., Desai, P., Ji, Y., Cummings, C. A., Shih, R., Degoricija, L., et al. (2012a). Comparative analysis of genome sequences covering the seven Cronobacter species. PLoS One 7:e49455. doi: 10.1371/journal.pone.0049455

PubMed Abstract | CrossRef Full Text | Google Scholar

Joseph, S., Sonbol, H., Hariri, S., Desai, P., McClelland, M., and Forsythe, S. J. (2012b). Diversity of the Cronobacter genus as revealed by multilocus sequence typing. J. Clin. Microbiol. 50, 3031–3039. doi: 10.1128/JCM.00905-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Kucerova, E., Joseph, S., and Forsythe, S. (2011). The Cronobacter genus: ubiquity and diversity. Q. Assur. Safety Crops Foods 3, 104–122. doi: 10.1111/j.1757-837x.2011.00104.x

CrossRef Full Text | Google Scholar

Leplae, R., Lima-Mendez, G., and Toussaint, A. (2010). ACLAME: a CLAssification of mobile genetic elements, update 2010. Nucleic Acids Res. 38, D57–D61. doi: 10.1093/nar/gkp938

PubMed Abstract | CrossRef Full Text | Google Scholar

Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245. doi: 10.1093/nar/gkw290

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Zeng, H., Zhang, J., He, W., Ling, N., Chen, M., et al. (2019). Prevalence, antibiotic susceptibility, and molecular characterization of Cronobacter spp. isolated from edible mushrooms in china. Front. Microbiol. 10:283. doi: 10.3389/fmicb.2019.00283

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Li, P., Xie, J., Yi, S., Yang, C., Wang, J., et al. (2014). New clustered regularly interspaced short palindromic repeat locus spacer pair typing method based on the newly incorporated spacer for Salmonella enterica. J. Clin. Microbiol. 52, 2955–2962. doi: 10.1128/JCM.00696-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Wang, X., Yin, K., Hu, Y., Xu, H., Xie, X., 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. doi: 10.1016/j.ijfoodmicro.2017.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Ling, N., Li, C., Zhang, J., Wu, Q., Zeng, H., He, W., et al. (2018). Prevalence and molecular and antimicrobial characteristics of Cronobacter spp. isolated from raw vegetables in china. Front. Microbiol. 9:1149. doi: 10.3389/fmicb.2018.01149

PubMed Abstract | CrossRef Full Text | Google Scholar

Makarova, K. S., Wolf, Y. I., Alkhnbashi, O. S., Costa, F., Shah, S. A., Saunders, S. J., et al. (2015). An updated evolutionary classification of CRISPR-Cas systems. Nat. Rev. Microbiol. 13, 722–736. doi: 10.1038/nrmicro3569

PubMed Abstract | CrossRef Full Text | Google Scholar

Masood, N., Moore, K., Farbos, A., Paszkiewicz, K., Dickins, B., McNally, A., et al. (2015). Genomic dissection of the 1994 Cronobacter sakazakii outbreak in a french neonatal intensive care unit. BMC Genomics 16:750. doi: 10.1186/s12864-015-1961-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Nethery, M. A., and Barrangou, R. (2019). CRISPR Visualizer: rapid identification and visualization of CRISPR loci via an automated high-throughput processing pipeline. RNA Biol. 16, 577–584. doi: 10.1080/15476286.2018.1493332

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogrodzki, P., and Forsythe, S. J. (2016). CRISPR-cas loci profiling of Cronobacter sakazakii pathovars. Future Microbiol. 11, 1507–1519. doi: 10.2217/fmb-2016-0070

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogrodzki, P., and Forsythe, S. J. (2017). DNA-Sequence based typing of the Cronobacter genus using mlst, CRISPR-Cas array and capsular profiling. Front. Microbiol. 8:1875. doi: 10.3389/fmicb.2017.01875

PubMed Abstract | CrossRef Full Text | Google Scholar

Paez-Espino, D., Morovic, W., Sun, C. L., Thomas, B. C., Ueda, K., Stahl, B., et al. (2013). Strong bias in the bacterial CRISPR elements that confer immunity to phage. Nat. Commun. 4:1430. doi: 10.1038/ncomms2440

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | CrossRef Full Text | Google Scholar

Shariat, N., and Dudley, E. G. (2014). CRISPRs: molecular signatures used for pathogen subtyping. Appl. Environ. Microbiol. 80, 430–439. doi: 10.1128/AEM.02790-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Shmakov, S., Smargon, A., Scott, D., Cox, D., Pyzocha, N., Yan, W., et al. (2017). Diversity and evolution of class 2 CRISPR-Cas systems. Nat. Rev. Microbiol. 15, 169–182. doi: 10.1038/nrmicro.2016.184

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, J. N., Vasconcellos, L., Forsythe, S. J., de Filippis, I., and Luiz Lima Brandao, M. (2019). Molecular and phenotypical characterization of Cronobacter species isolated with high occurrence from oats and linseeds. FEMS Microbiol. Lett. 366:fny289. doi: 10.1093/femsle/fny289

PubMed Abstract | CrossRef Full Text | Google Scholar

Streicher, E. M., Victor, T. C., van der Spuy, G., Sola, C., Rastogi, N., van Helden, P. D., et al. (2007). Spoligotype signatures in the Mycobacterium tuberculosis complex. J. Clin. Microbiol. 45, 237–240. doi: 10.1128/JCM.01429-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, C. L., Thomas, B. C., Barrangou, R., and Banfield, J. F. (2016). Metagenomic reconstructions of bacterial CRISPR loci constrain population histories. ISME J. 10, 858–870. doi: 10.1038/ismej.2015.162

PubMed Abstract | CrossRef Full Text | Google Scholar

Ueda, S. (2017). Occurrence of Cronobacter spp. in dried foods, fresh vegetables and soil. Biocontrol Sci. 22, 55–59. doi: 10.4265/bio.22.55

PubMed Abstract | CrossRef Full Text | Google Scholar

Westra, E. R., Buckling, A., and Fineran, P. C. (2014). CRISPR-Cas systems: beyond adaptive immunity. Nat. Rev. Microbiol. 12, 317–326. doi: 10.1038/nrmicro3241

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Li, C., Wu, Q., Zhang, J., Huang, J., and Yang, G. (2015). Prevalence, molecular characterization, and antibiotic susceptibility of Cronobacter spp. in chinese ready-to-eat foods. Int. J. Food Microbiol. 204, 17–23. doi: 10.1016/j.ijfoodmicro.2015.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, S., Jensen, M. A., Bai, J., Debroy, C., Barrangou, R., and Dudley, E. G. (2013). The evolutionary divergence of Shiga toxin-producing Escherichia coli is reflected in clustered regularly interspaced short palindromic repeat (CRISPR) spacer composition. Appl. Environ. Microbiol. 79, 5710–5720. doi: 10.1128/AEM.00950-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Yong, W., Guo, B., Shi, X., Cheng, T., Chen, M., Jiang, X., et al. (2018). An investigation of an acute gastroenteritis outbreak: Cronobacter sakazakii, a potential cause of food-borne illness. Front. Microbiol. 9:2549. doi: 10.3389/fmicb.2018.02549

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, H., Lei, T., He, W., Zhang, J., Liang, B., Li, C., et al. (2018a). Novel multidrug-resistant Cronobacter sakazakii causing meningitis in neonate, China, 2015. Emerg. Infect. Dis. 24, 2121–2124. doi: 10.3201/eid2411.180718

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, H., Zhang, J., Wu, Q., He, W., Wu, H., Ye, Y., et al. (2018b). Reconstituting the history of Cronobacter evolution driven by differentiated CRISPR activity. Appl. Environ. Microbiol. 84:e00267-18. doi: 10.1128/AEM.00267-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, H., Zhang, J., Li, C., Xie, T., Ling, N., Wu, Q., et al. (2017). The driving force of prophages and CRISPR-Cas system in the evolution of Cronobacter sakazakii. Sci. Rep. 7:40206. doi: 10.1038/srep40206

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Abadia, E., Refregier, G., Tafaj, S., Boschiroli, M. L., Guillard, B., et al. (2010). Mycobacterium tuberculosis complex CRISPR genotyping: improving efficiency, throughput and discriminative power of ‘spoligotyping’ with new spacers and a microbead-based hybridization assay. J. Med. Microbiol. 59(Pt 3), 285–294. doi: 10.1099/jmm.0.016949-0

PubMed Abstract | CrossRef Full Text | Google Scholar