CRISPR Typing Increases the Discriminatory Power of Streptococcus agalactiae Typing Methods

We explored the relevance of a Clustered regularly interspaced short palindromic repeats (CRISPR)-based genotyping tool for Streptococcus agalactiae typing and we compared this method to current molecular methods [multi locus sequence typing (MLST) and capsular typing]. To this effect, we developed two CRISPR marker schemes (using 94 or 25 markers, respectively). Among the 255 S. agalactiae isolates tested, 229 CRISPR profiles were obtained. The 94 and 25 markers made it possible to efficiently separate isolates with a high diversity index (0.9947 and 0.9267, respectively), highlighting a high discriminatory power, superior to that of both capsular typing and MLST (diversity index of 0.9017 for MLST). This method has the advantage of being correlated with MLST [through analysis of the terminal direct repeat (TDR) and ancestral spacers] and to possess a high discriminatory power (through analysis of the leader-end spacers recently acquired, which are the witnesses of genetic mobile elements encountered by the bacteria). Furthermore, this “one-shot” approach presents the benefit of much-reduced time and cost in comparison with MLST. On the basis of these data, we propose that this method could become a reference method for group B Streptococcus (GBS) typing.

Introduction

Streptococcus agalactiae or group B Streptococcus (GBS) is the leading cause of neonatal infections and an emerging pathogen in adults, particularly in elderly and immunocompromised patients (Phares et al., 2008; Skoff et al., 2009; Slotved and Hoffmann, 2020; Vuillemin et al., 2021). It is also a commensal bacteria that colonize the gastrointestinal and genitourinary tracts of 10–30% of healthy humans (Regan et al., 1991; Campbell et al., 2000; Khalil et al., 2017; Cho et al., 2019). Prevention of GBS-related neonatal infections involves the detection of vaginal carriage in pregnant women, followed by intrapartum antibiotic prophylaxis for those colonized (ANAES, 2001; Schrag et al., 2002).

To characterize GBS isolates, two main typing methods are widely used: serotyping and MLST (multi locus sequence typing). Serotyping methods, which are based on differences in capsular polysaccharides (CPS), were first evaluated by phenotypic methods and then by PCR-based molecular methods (Imperi et al., 2010). To date, 10 serotypes (or capsular types) have been identified (Ia, Ib, and II-IX) (Cieslewicz et al., 2005; Slotved et al., 2007). The prevalence and distribution of serotypes are known to differ between geographical regions, ethnic populations, and clinical presentations (Schuchat, 1998; Kwatra et al., 2016). MLST, a sequence-based method, is now widely used to investigate the population structure and genetic lineage of GBS and is currently the reference method for GBS typing. This approach is based on the combination of alleles for seven housekeeping genes (Jones et al., 2003). Unique combinations of the alleles at each locus define the allelic profiles, or sequence types (STs). ST can be clustered in clonal complexes (CCs) when six of the seven alleles are in common (Feil et al., 2004). These methods have highlighted the involvement of serotype III and ST17, in causing more invasive neonate diseases (Manning et al., 2009). Among adults, serotype V isolates belonging to CC1 have been associated with invasive disease (Bergseng et al., 2008; Skoff et al., 2009). More recently, serotypes III, Ia, and IV have gained relevance in this context (Tazi et al., 2011; Lamagni et al., 2013; Teatero et al., 2014, 2015). Despite the clonality observed within GBS populations causing invasive disease, increasing diversity has appeared, hence the need to investigate the population structure. However, between these two typing methods, while serotyping is relatively easy to perform in the laboratory, it is insufficiently discriminating to compare isolates. MLST presents significant disadvantages: it is relatively costly, time-consuming, and labor-intensive. Moreover, as the sequences targeted evolve slowly, MLST is not highly discriminant for epidemiological studies and local surveillance, and isolates are not easy to distinguish at the ST level (Radtke et al., 2010; Haguenoer et al., 2011; Sabat et al., 2013).

For typing, whole genome sequencing (WGS) provides an ideal resolution and accuracy (Schürch et al., 2018; Kayansamruaj et al., 2019; Beauruelle et al., 2020; Perme et al., 2020). This method shows higher discriminatory power than capsular typing and MLST. However, although WGS is increasingly accessible, it is still an expensive technology requiring experience and skill to be used, including bioinformatics analysis. Indeed, WGS is not yet useful for epidemiological studies or local surveillance. An accessible typing method for GBS strains is definitely needed and useful at this moment, especially for low-to-middle-income country and in the view of the advanced development stage of GBS vaccines, which will require robust and continuous surveillance worldwide. At the same time, other approaches have been developed such as the analysis of patterns of virulence gene, prophages content, or CRISPR-Cas [Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR associated sequences (Cas)] analysis (Van der Mee-Marquet et al., 2006; Lopez−Sanchez et al., 2012; Lier et al., 2015; Crestani et al., 2020).

CRISPR-Cas is an adaptive and vertically transmitted immune system present in a large proportion of prokaryotic genomes (Barrangou et al., 2007). CRISPR arrays are made up of a succession of highly conserved repeated sequences, called direct repeats or DRs, interspaced by sequences of a similar length, called spacers. Most CRISPR arrays are flanked on one side by the leader sequence and on the other side by a trailer sequence, which is preceded by the terminal direct repeat (TDR) of the CRISPR array, corresponding to a degenerated or truncated DR. DRs are highly conserved within a locus, whereas spacers vary widely as they derive from foreign genetic elements. In this bacterial immune system, each spacer is acquired in response to mobile genetic element (MGE) mostly considered as invasive elements by the bacteria. Spacers are integrated in a linear, time-oriented manner. Each new spacer is incorporated at the leader end of the CRISPR array, concomitantly to the duplication of a DR. Because of the polarized acquisition of spacers deriving from encounters with various MGE over time, CRISPR arrays constitute a chronological archive of past encounters. Thus, spacers located at the leader-end extremity are recently acquired and represent recent contacts with MGE, whereas spacers located at the trailer-end extremity represent ancestral spacers. Leader-end spacers enable the differentiation of closely related strains and trailer-end spacers reflect broader phylogenetic relationships. Indeed, several studies have shown potential for CRISPR-based typing (Barrangou and Dudley, 2016).

In GBS, two CRISPR-Cas systems have been characterized, a type I-C system named CRISPR2, which is rare and most often incomplete, suggesting little or no activity, and a type II-A system, named CRISPR1, which is ubiquitous and functional (Lopez−Sanchez et al., 2012). The CRISPR1 array is made up of highly conserved DR of 36 bp, separated by spacers of approximately 30 bp, which are extremely diverse in sequence and in number across GBS strains studied. Similarities between CRISPR1 spacer sequences and MGEs have been previously reported (Lopez−Sanchez et al., 2012; Lier et al., 2015). Comparative sequence analysis across numerous GBS isolates emphasized that CRISPR1 array is extremely diverse and evolves in vivo, demonstrating the dynamics of the system (Lopez−Sanchez et al., 2012; Lier et al., 2015; Beauruelle et al., 2017, 2018).

A typing method based on CRISPR array analysis has been shown to be a useful tool for comparing GBS isolates (Lopez−Sanchez et al., 2012; Lier et al., 2015; Beauruelle et al., 2017, 2018; Gajic et al., 2019). This method presents the advantage of being linked to the genetic lineage defined by MLST through analysis of the TDR and ancestral spacers. Moreover, exploration of the recent evolution of the isolate, especially encounters with MGE, is possible through the analysis of the leader-end spacers. Based on the spacers’ variability, it seems to be in congruence with MLST with a greater discriminating power. However, these advantages still have to be confirmed.

The aim of our study was to evaluate CRISPR1 analysis as a high-resolution S. agalactiae typing method and to compare this method to capsular typing and MLST.

Materials and Methods

Bacterial Isolates

Two hundred and fifty-five S. agalactiae isolates were analyzed. A total of 224 isolates were collected at the University Hospital of Tours, France, between 2002 and 2016; 26 were collected by the National Reference Center for Streptococci in 2012 and five were reference strains. Among the clinical isolates, 198 were isolated from adults and 52 were isolated from neonates. Among them, 164 were isolated in a non-invasive setting from cutaneous lesions (n = 2), urinary tract (n = 5), gastric aspirate (n = 7), respiratory tract (n = 9), genital tract (n = 119), or intestinal tract (n = 23), and 86 were isolated in an invasive setting from blood (n = 49), CSF (n = 33), or joint fluid cultures (n = 2). Isolates were grown on GBS selective media or blood agar plates. Identification was performed using MALDI-TOF MS (Vitek MS, BioMerieux France, or Bruker Daltonics Germany). The five reference strains considered were 2,603 V/R (capsular type V, ST110), A909 (capsular type Ia, ST7), NEM316 (capsular type III, ST23), COH1 (capsular type III, ST17), and BM110 (capsular type III, ST17) (Glaser et al., 2002; Tettelin et al., 2005; Da Cunha et al., 2014; Supplementary Table 1).

DNA Extraction

Genomic DNA was extracted following enzymatic lysis with mutanolysin (Sigma). A bacterial suspension of 1.5 McFarland was prepared in 500 μl of water containing 50 U of mutanolysin. The suspension was incubated for 1 h at 56°C, followed by 10 min at 100°C, leading to cell lysis. Lysates were centrifuged for 3 min at 1,500 × g and the supernatants containing DNA were collected.

Capsular Typing

Capsular typing was performed by PCR-based methods based on previously described methods (Imperi et al., 2010). The assay identified each serotype (Ia to IX) by analysis of band patterns on agarose gel.

Multilocus Sequence Typing

Multi locus sequence typing was carried out as previously described (Jones et al., 2003). Allelic profiles and ST were assigned using the international MLST database¹. CCs were defined using the stringent group definition (6/7 shared alleles) and eBURST analysis² applied to the 250 isolates of the study.

CRISPR1 Array Amplification and Sequencing

CRISPR1 array amplification was performed using CRISPR1 PCR-F and CRISPR1 PCR-R primers that target the CRISPR1-flanking regions as previously described (Lopez−Sanchez et al., 2012; Lier et al., 2015; Beauruelle et al., 2017). PCR products were sequenced using the internal sequencing primers CRISPR1 SEQ-F and CRISPR1 SEQ-R. For CRISPR1 regions exceeding 1.3 kb, primers targeting internal spacers were used to complete the sequencing. Primer sequences are presented in Supplementary Table 1.

CRISPR1 Array Analysis

Spacers, repeats, and flanking regions for each sequence were identified using a macro-enabled Excel tool (P. Horvath, DuPont). This tool allows the identification and extraction of CRISPR features from nucleotide sequences, and the graphic representation of spacers as colored cells in Excel spreadsheets. The macro-enabled Excel tool is optional and sequence obtained could be analyzed without it by the identification of the DR sequences (Supplementary Table 2) separating each spacer sequence. Spacer sequences were compared to the dictionary of spacers established earlier and expanded previously (Lopez−Sanchez et al., 2012; Beauruelle et al., 2017, 2018). New spacers identified in this study further expanded the dictionary and were numbered incrementally. The original contributions presented in the study are publicly available. These data can be found at http://crispr.i2bc.paris-saclay.fr/CRISPRcompar/Dict/Dict.php.

Spacers and TDR Selection

Graphic representation of CRISPR1 arrays made it possible to separate isolates according to their CRISPR1 array composition. Isolates were clustered according to their ancestral spacers and TDR composition. Among each cluster, specific markers, corresponding to spacers and TDR, were selected in view of separating isolates according to their CRISPR1 array similarity. Markers were selected according to their frequency among CRISPR array and their specificity to each phylogenetic group. Markers were first selected visually thanks to the macro-enabled Excel tool used. Markers selected were then analyzed using a binary code (Supplementary Table 3) to (i) evaluate the frequency of the marker selected among isolates and among each group, (ii) evaluate the absence of redundancy between markers selected, and (iii) evaluate the specificity of each marker for a group. Markers selected had to be widely present and specific of each CRISPR array belonging to a given phylogenetic group. Markers selected were present in at least three-quarters of the isolates of the group, were non-redundant, and were specific of each group considered. Our aim was to obtain the best compromise between marker number and discriminatory index (DI).

Data Analysis

The DI described by Hunter and Gaston was used as a numerical index for the discriminatory power of each typing method (Hunter and Gaston, 1988). The categorical coefficient, unweighted pair group with arithmetic mean (UPGMA), and the minimum spanning tree (MST) were run using BioNumerics 7.6.2 software (Applied Maths, Sint-Martens-Latem, Belgium). In MST, each circle represents a CRISPR genotype or a CC and its size is proportional to the number of isolates. The thicker branches link the genotypes differing by only one spacer, the thinner branches link genotypes differing by more than one spacer. Congruence between CRISPR1 typing, MLST, and serotyping was calculated using Rand Index (BioNumerics software).

Results

Capsular Typing and Multi Locus Sequence Typing

Capsular Typing

Among the 255 S. agalactiae isolates analyzed, six capsular types were represented. The main capsular type was type III (n = 93), followed by type Ia (n = 56), V (n = 54), Ib (n = 23), II (n = 17), and IV (n = 12). Capsular types VI, VII, VIII, and IX were not found among our isolates (Table 1 and Supplementary Table 4).

TABLE 1

Table 1. Characteristics of the 255 GBS isolates tested classified by sequence type.

MLST

A total of 48 different STs were found among the 255 isolates. The most common was ST17 (n = 56), followed by ST23 (n = 37), ST1 (n = 32), ST19 (n = 26), ST8 (n = 12), and ST28 (n = 10). Other STs were represented by fewer isolates (<10) (Table 1). Invasive isolates belonged mainly to ST17 followed by ST23 and ST1, while non-invasive isolates were distributed mainly among ST23, ST17, ST1, and ST19 isolates. eBURST analysis clustered the ST into seven CCs [CC17 (grouping ST17, ST291, and ST1004), CC1 (grouping ST1, ST196, ST2, ST1005, ST459, ST136, ST370, ST297, and ST173), CC23 (grouping ST23, ST220, ST1002, ST481, ST385, ST88, ST144, ST366, and ST391), CC19 (grouping ST19, ST28, ST27, ST389, ST182, and ST386), CC8 (grouping ST8, ST12, ST10, and ST390), CC6 (grouping ST6, ST7, ST255, and ST41), and CC4 (grouping ST4, ST3, and ST243)], one group with two singletons ST (ST26 and ST388), and five other singletons (ST22, ST24, ST130, ST327, and ST569) (Supplementary Table 5). All CCs except CC17 contained mainly non-invasive isolates [CC19 (n = 31, 77.5%), CC1 (n = 38, 74.5%), CC23 (n = 34, 70.8%), CC6 (n = 4, 66.7%), CC8 (n = 17, 63%), CC4 (n = 3, 60%), and CC17 (n = 22, 39.3%)] (Supplementary Table 5). Using UPGMA and MST, the 48 different profiles (corresponding to the 48 STs) were grouped into 21 clusters with more than one isolate (up to 54 isolates by cluster) (Supplementary Figure 5).

CRISPR1 Array Analysis

We generated a complete CRISPR1 sequence for all the 255 isolates. The number of spacers ranged from three to 29 per isolate, corresponding to a CRISPR array size of 266–1,800 bp. Among the 255 isolates, a specific CRISPR1 array was observed for 92% of isolates (n = 229). For the other 8%, the same CRISPR1 array was common for two isolates (n = 22) or for four isolates (n = 4). Using the macro-enabled Excel tool, isolates could be distributed according to their CRISPR1 array homology (Supplementary Figure 1). First, isolates were grouped into six clusters according to their TDR and ancestral spacers. These six clusters showed a correlation between CC defined by MLST and eBURST analysis: a CC23 cluster; a cluster grouping CC17 and ST130 isolates; a cluster grouping CC1, CC4, and CC19 isolates; a cluster grouping CC6 and CC8 isolates; a cluster grouping ST26–ST388 isolates; and a ST22 cluster (Figure 1). CRISPR1 clusters were then divided into subgroups according to more recently acquired spacers. In this way, a total of 14 clusters and subgroups were defined, namely, three subgroups in cluster CC1–CC4–CC19, allowing separation between CC19 (two subgroups) and the group CC1–CC4 (two subgroups), five subgroups in the CC17 cluster (COH1 type, BM110 type, ST130 type, and two other subgroups), and two subgroups in the CC23 cluster (Figure 1 and Supplementary Figure 1). This CRISPR genotyping approach can be used to rank some singletons defined by eBURST analysis into clusters: ST24 (CRISPR1 cluster CC23), ST130 (CRISPR1 cluster CC17), ST569 (CRISPR1 cluster CC6–CC8), and ST327 (CRISPR1 cluster CC19) (Supplementary Table 5 and Supplementary Figure 1).

FIGURE 1

Figure 1. CRISPR clustering according to TDR and ancestral spacers. The CRISPR1 arrays are represented using a macro-enabled Excel tool, whereby spacers are converted into two-color symbols based on spacer sequence. Repeats are not shown except terminal direct repeats (TDRs), which are represented by different colored borders according to their sequence. Spacers were identified by a number attributed following the spacer dictionary (http://crispr.i2bc.paris-saclay.fr/CRISPR compar/Dict/Dict.php). Arrays are oriented with respect to the leader sequence located on the left. A total of 14 clusters were obtained. The six groups were defined by TDR and ancestral spacers. These groups were then divided into subgroups, based on more recent spacers.

Among each cluster, spacers and TDR were selected to define phylogenetic lineages. Each spacer selected was cluster- and subgroup-specific. We selected two groups of markers. The first involved 94 markers, namely, the five different TDRs and 89 spacers (Supplementary Table 3A). The second involved 25 markers, namely, the five different TDRs and 20 spacers (Supplementary Table 3B). Among the 255 isolates, using the UPGMA algorithm, the 94-marker scheme defined a total of 172 different profiles and 45 groups with more than one isolate (up to eight isolates) (Supplementary Figures 3, 4). Likewise, the 25-marker selection defined a total of 42 different profiles and 29 clusters with more than one isolate (up to 40 isolates) (Figure 2 and Supplementary Figure 4). Using these two-marker selections, groups previously defined were divided in different clusters (Figure 2 and Supplementary Figure 2). Thus, CC17 isolates were divided into 30 groups and five groups, respectively, using these two-marker selections (94 and 25 markers, respectively).

FIGURE 2

Figure 2. Minimum spanning tree (MST) representation of the 25 CRISPR1 markers scheme clustering. Each circle represents a CRISPR1 genotype and its size is proportional to the number of strains. Each color represents CC or singleton defined by MLST (e.g., yellow for CC17). A high level of correlation between this marker selection and MLST type was observed; circles (representing CRISPR1 genotype) are characterized mostly by a same color, especially for CC17 and CC23, whereas isolates belonging to CC1 and CC19 are more dispersed.

Differences Among Clusters

Some groups displayed lower spacer composition diversity. Among the 26 isolates sharing a mutual CRISPR1 profile, 17 belonged to CC17, and the nine others belonged to CC8 (n = 2), CC23 (n = 3), CC19 (n = 2), and ST388 (n = 2). Indeed, the majority of CRISPR1 arrays sharing a mutual profile belonged to CC17 isolates and represented 31% (17/54) of them. Among these CC17 isolates with common CRISPR1 array, the main isolates (76%, 13/17) belonged to the first CRISPR1 subgroup (COH1 type) and represented 41% of isolates of the COH1 type subgroup (13/32). By contrast, isolates belonging to CC1 and CC19 shared a high degree of polymorphism and could not be easily clustered into their specific CC defined by MLST and eBURST (Supplementary Figure 1).

Discrimination of Isolates by CRISPR Analysis

CRISPR1 array analysis makes it possible to separate isolates within a same ST or capsular type. Among the 255 isolates, six capsular types and 48 different STs were observed, while CRISPR1 analysis separated them into 229 different CRISPR1 profiles, and the two marker schemes into 172 (94-marker scheme) and 42 (25-marker scheme) profiles. The diversity index was compared for the two CRISPR1 specific marker schemes as well as for MLST. The diversity index was 0.9947 for the CRISPR1-specific 94-marker scheme, 0.9267 for the CRISPR1-specific 25-marker scheme, and 0.9017 for the MLST method on this population (Table 2).

TABLE 2

Table 2. Hunter and Gaston diversity index of CRISPR typing scheme with 94 markers or 25 markers and MLST.

Congruence Between CRISPR1 Typing, MLST, and Capsular Typing

Congruence between the two CRISPR1-specific markers, MLST and capsular typing, was analyzed using the Rand Index (Supplementary Table 6). A strong congruence was observed between the 25-marker scheme and the 94-marker scheme (Rand Index: 0,932). The Rand Index between CRISPR1 marker selections and MLST was 0.904 and 0.902 for the 25-marker scheme and the 94-marker scheme, respectively. The Rand Index between CRISPR1 marker selections and capsular typing was 0.797 and 0.766 for the 25-marker scheme and the 94-marker scheme, respectively. The Rand Index between MLST and capsular typing was 0.84.

Discussion

The aim of our study was to evaluate CRISPR1 analysis as a high-resolution S. agalactiae typing method. We explored the relevance of the CRISPR1-based genotyping tool, and we compared this method to current molecular standards. We developed two schemes of markers selection (94 or 25 selected markers corresponding to spacers and TDR), and we simultaneously characterized these isolates by MLST and capsular typing. We used a representative library of the species by selecting a wide variety of isolates that differed by their anatomic and geographic origin, as well as by their phylogenetic origin.

Markers were selected in view of clustered isolates according to their CRISPR1 array similarity. TDR and ancestral spacers allow the isolates to be linked to their common ancestor. More recent spacers allow separating isolates in the subgroup within this first selection. We defined two schemes of marker selection. The largest selection allowed the clustering of the majority of isolates. To facilitate data mining, we tried to reduce the number of markers as much as possible. A number of 25 markers appears sufficient to discriminate isolates successfully. The discriminatory power of this CRISPR1 approach is superior to that of both capsular typing and MLST (diversity index of 0.9947 and 0.9267 for the two marker schemes vs. 0.9017 for MLST). The CRISPR1-typing approach and the two-CRISPR1 marker selection, including the simplest, allow one to efficiently separate isolates and to successfully discriminate isolates that were considered indistinguishable by MLST. Furthermore, compared to publish data, the CRISPR1 discriminatory power appears superior to that of pulsed-field gel electrophoresis, which is considered highly discriminatory (Simpson’s diversity index = 0.92) (Pillai et al., 2009). Aside from its high discriminatory power, this CRISPR1 typing approach assesses the phylogenetic structure of the S. agalactiae population, with the two schemes of CRISPR1 marker selection showing a clonal distribution of the population similar to that obtained by MLST. The CRISPR1 clustering approach generated major clusters that corresponded well to the main CC obtained by MLST and eBURST analysis. Moreover, this approach presents the major advantage of being highly discriminant thanks to the variability of the CRISPR1 array. The high DI of the CRISPR1 approach made it possible to distinguish between isolates within the same ST or CC. This approach delivers a unique DNA fingerprint and makes it possible to separate isolates, even clonal bacteria such as S. agalactiae. A strong congruence was observed between the 25-marker scheme and the 94-marker scheme (Rand Index: 0.932), as well as between the two marker schemes and MLST (Rand Index: 0.904 and 0.902 for the 25-marker scheme and the 94-marker scheme, respectively). The congruence between CRISPR1 marker selection and capsular type was lower (0.797 and 0.766 for the 25-marker scheme and the 94-marker scheme, respectively). This could probably be explained by the high discriminatory power of the CRISPR typing technique compared to serotyping.

Previous studies highlight the interest of CRISPR for GBS typing (Lopez−Sanchez et al., 2012; Lier et al., 2015; Beauruelle et al., 2017, 2018; Gajic et al., 2019). In particular, these studies highlight the correlation between CRISPR and MLST and the strong discriminatory power of CRISPR typing. In the present study, we confirm the interest of GBS CRISPR typing by specifying the added value of this method compared to other commonly used, including those based on the calculation of diversity and congruence index. Moreover, we show that this discriminatory power of CRISPR1 typing is also applicable with a limited number of markers, offering a good compromise between discriminatory power, phylogenetic data, and simplicity.

We then challenged our marker selection with other groups of isolates based on CRISPR array publicly available from two publications (Lopez−Sanchez et al., 2012; Gajic et al., 2019). Gajic et al. (2019) characterized 87 CRISPR1 array from GBS isolates from Serbia (invasive and non-invasive human isolates). Following our spacer selection, all of the isolates from Gajic et al. (2019) study could be clustered in the groups and 90% of them in the subgroups. All except one (Serbie 41 isolate) could not be clustering using our two marker schemes. We also evaluated our spacer selection with the isolates from the Lopez−Sanchez et al. (2012). study. This publication presents the great advantage of including CRISPR1 array of GBS isolated from animals and from several geographic regions. Among the 351 GBS isolates, 87% (n = 306) could be clustered in the groups and 80% in the subgroups. Marker selections allow clustering in subgroups 97% and 86% of isolates with the 94- and 25-marker scheme, respectively. The 45 isolates, almost exclusively animal isolates, which could not be clustered in a group, belonged to CC340, CC260, CC61–67, and CC103. In our study, these CCs were not present and could not be classified using our marker selections. Further studies are needed to enlarge the diversity of GBS isolates studied to be able to classify isolates belonging to less frequent or peculiar lineages.

Within CRISPR1 locus analysis, CRISPR1 diversity differed among clusters. CC17 isolates shared a low degree of polymorphism compared to other CCs, as noted previously (Lier et al., 2015; Beauruelle et al., 2017). This moderate CRISPR1 diversity could be explained by the slow rate of evolution of ST17 isolates, characterized by a low rate of recombination, which is known for contributing to CRISPR array diversity (Da Cunha et al., 2014). Another explanation might be that CC17 isolates encode a specific repertoire of surface proteins, suggesting a specific colonization site, and therefore a genetic isolation (Da Cunha et al., 2014). Conversely, some CCs shared a high degree of polymorphism and were difficult to cluster using CRISPR1 array. This relates in particular to isolates belonging to CC1, CC4, and CC19. Interestingly, CC1 isolates are increasing in adult invasive infections since the 1990s, suggesting the evolution of this CC, especially ST1 isolates (Bergseng et al., 2008; Skoff et al., 2009). We can hypothesize that this CRISPR1 diversity could be due to (i) a more rapid evolution of these CC isolates, (ii) a more active CRISPR-Cas system, or (iii) a higher diversity of MGE attacking this CC leading to enhanced CRISPR immunization of this phylogenetic lineage. Nevertheless, lineages defined by CRISPR1 analysis, including in less variable CCs, correspond to different lineages previously highlighted by WGS, which is undoubtedly the gold standard to compare isolates (Da Cunha et al., 2014). Indeed, according to TDR and ancestral spacers, CC17 isolates could be divided into four subgroups and CC23 isolates could be divided into two subgroups (Figure 1) as previously described by SNP analysis (Da Cunha et al., 2014). Furthermore, WGS has previously confirmed CRISPR1-based clustering for the S. agalactiae population (Kayansamruaj et al., 2019; Beauruelle et al., 2020). Conversely, some CCs shared a high degree of polymorphism and were difficult to cluster using CRISPR1 array. This relates in particular to isolates belonging to CC1, CC4, and CC19, suggesting a greater heterogeneity of this population. This may be due to the high recombination rate in some CCs as previously described (Da Cunha et al., 2014).

CRISPR1 array analysis is easy to perform with just one array of limited length (266 to 1,800 bp) to be analyzed (unlike MLST, which requires seven loci to be sequenced) and offered supplementary information than that obtained by MLST and capsular typing combined. Moreover, this “one-shot” approach generates data (spacers dictionary) that enable comparison of isolates between different laboratories. Considering these data, CRISPR1 typing appears as an ad hoc tool to compare isolates and analyzed GBS transmission. Nowadays, although incidence of late onset neonatal GBS disease and adult GBS disease is increasing (Bekker et al., 2014; Ballard et al., 2016), their transmission routes are poorly understood. This CRISPR1 based approach could be a useful tool to explore source of transmission of these GBS infections. Moreover, this approach has the ability to assess the genetic relatedness among these isolates and to provide a better understanding of the physiopathology of these infections. Whereas the analysis of ancestral spacer and TDR allows typing, the analysis of spacers at the leader end (recently acquired spacers) allows subtyping and provides specific evidence on the recent evolution of isolates, especially encounters with MGEs. MGEs are key factors for the evolution of bacteria, including for GBS as highlighted by the insertion of integrative and conjugative elements that caused the expansion of few clones in human, particularly adapted to their host (Da Cunha et al., 2014). Indeed, following GBS–MGEs contact, CRISPR analysis gives us valuable clues. Similarly, GBS possesses a broad animal host spectrum and studies proved that some GBS genotypes can cause human invasive diseases through animal sources as food-borne zoonotic infections (Tan et al., 2016). A deep phylogenetic analysis such as CRISPR typing appears useful to analyze the circulation of different GBS genotypes in humans and animals in different countries and to monitor potential emerging zoonotic GBS clones. CRISPR1-based typing could be used to explore the genetic relatedness among humans and animal isolates such as cattle and fish, which are also an important source of GBS infection (Radtke et al., 2010; Liu et al., 2016). In the present work, we selected a wide variety of human isolates but we did not analyze animal isolates, which is one limitation of this work.

Another interest to this CRISPR1-based typing method was the clinical evaluation of GBS vaccine. The advanced development stage of GBS vaccine requires robust and continuous surveillance worldwide, including in low-income countries. This CRISPR approach has already proven its effectiveness to evaluate evolution and diversity of GBS vaginal carriage (Beauruelle et al., 2017, 2018). Indeed, this low-cost and easy-to-use method appears useful to evaluate the diversity of the species, including in low-income countries where there is little information available about the GBS isolate characteristics. The majority of vaccine development was based on polysaccharide conjugate vaccine, and evaluation of GBS diversity is usually based on the capsular type (Berner, 2021; World Health Organization, 2021). However, serotype replacement should be kept in mind (Meehan et al., 2014). Moreover, GBS vaccine efficacy on different GBS genotypes as well as GBS carriage evolution over time (same or new isolate) could be evaluated with this tool.

Although WGS is increasingly accessible, it is still an expensive technology requiring experience and skill to be used, including bioinformatics analysis. Nowadays, whole-genome sequencing still does not appear to be a routine method for genotyping, thus rendering CRISPR-based typing technology useful during the transition from the current molecular typing method to the omic level. Given these data, we assume that this method could become an actual reference method for phylogenetic GBS typing.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: http://crispr.i2bc.paris-saclay.fr/CRISPRcompar/Dict/Dict.php.

Author Contributions

CB and PL designed the study, developed methods, and wrote the manuscript. CB, LT, CL, PG, and CP selected bacterial isolates and performed the research. TC analyzed data and performed bioinformatic analysis. AP and LM provided critical feedback of the manuscript. All authors analyzed data, contributed to manuscript revision, read, and approved the submitted version.

Conflict of Interest

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.

Acknowledgments

We would like to thank P. Horvath for providing the macro-enabled Excel tool for CRISPR analysis and for helpful discussion and to the Institut de Biologie Integrative de la cellule (I2BC) for maintaining the web server of CRISPR (http://crispr.i2bc.paris-saclay.fr).

Supplementary Material

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

Footnotes

1. ^ http://pubmlst.org/sagalactiae/
2. ^ http://eburst.mlst.net/v3/mlst_datasets/

References

ANAES (2001). [Antenatal prevention of early neonatal bacterial infection (September 2001)]. J. Gynecol. Obstet. Biol. Reprod. 32, 68–74.

Google Scholar

Ballard, M. S., Schønheyder, H. C., Knudsen, J. D., Lyytikäinen, O., Dryden, M., Kennedy, K. J., et al. (2016). The changing epidemiology of group B Streptococcus bloodstream infection: a multi-national population-based assessment. Infect. Dis. 48, 386–391. doi: 10.3109/23744235.2015.1131330

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrangou, R., and Dudley, E. G. (2016). CRISPR-Based Typing and Next-Generation Tracking Technologies. Ann. 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

Beauruelle, C., Branger, M., Cochard, T., Pastuszka, A., Biet, F., and Lanotte, P. (2020). Whole-genome sequencing confirms the coexistence of different colonizing Group B Streptococcus isolates underscored by CRISPR typing. Microbiol. Resour. Announc. 9, e1359–e1319. doi: 10.1128/MRA.01359-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Beauruelle, C., Pastuszka, A., Horvath, P., Perrotin, F., Mereghetti, L., and Lanotte, P. (2017). CRISPR: A useful genetic feature to follow vaginal carriage of Group B Streptococcus. Front. Microbiol. 8:1981. doi: 10.3389/fmicb.2017.01981

PubMed Abstract | CrossRef Full Text | Google Scholar

Beauruelle, C., Pastuszka, A., Mereghetti, L., and Lanotte, P. (2018). Group B Streptococcus vaginal carriage in pregnant women deciphered by CRISPR analysis. J. Clin. Microbiol. 56, e1949–e1917. doi: 10.1128/JCM.01949-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Bekker, V., Bijlsma, M. W., van de Beek, D., Kuijpers, T. W., and van der Ende, A. (2014). Incidence of invasive group B Streptococcal disease and pathogen genotype distribution in newborn babies in the Netherlands over 25 years: a nationwide surveillance study. Lancet Infect. Dis. 14, 1083–1089. doi: 10.1016/s1473-3099(14)70919-3

CrossRef Full Text | Google Scholar

Bergseng, H., Rygg, M., Bevanger, L., and Bergh, K. (2008). Invasive Group B Streptococcus (GBS) disease in Norway 1996–2006. Eur. J. Clin. Microbiol. Infect. Dis. 27, 1193–1199. doi: 10.1007/s10096-008-0565-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Berner, R. (2021). Group B streptococcus vaccines: one step further. Lancet Infect Dis. 21, 158–160. doi: 10.1016/S1473-3099(20)30451-5

CrossRef Full Text | Google Scholar

Campbell, J. R., Hillier, S. L., Krohn, M. A., Ferrieri, P., Zaleznik, D. F., and Baker, C. J. (2000). Group B Streptococcal colonization and serotype-specific immunity in pregnant women at delivery. Obstet. Gynecol. 96, 498–503. doi: 10.1016/s0029-7844(00)00977-7

CrossRef Full Text | Google Scholar

Cho, C.-Y., Tang, Y.-H., Chen, Y.-H., Wang, S.-Y., Yang, Y.-H., Wang, T.-H., et al. (2019). Group B Streptococcal infection in neonates and colonization in pregnant women: An epidemiological retrospective analysis. J. Microbiol. Immunol. Infect. 52, 265–272. doi: 10.1016/j.jmii.2017.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cieslewicz, M. J., Chaffin, D., Glusman, G., Kasper, D., Madan, A., Rodrigues, S., et al. (2005). Structural and genetic diversity of Group B Streptococcus capsular polysaccharides. Infect. Immunit. 73, 3096–3103. doi: 10.1128/IAI.73.5.3096-3103.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Crestani, C., Forde, T. L., and Zadoks, R. N. (2020). Development and Application of a Prophage Integrase typing scheme for Group B Streptococcus. Front. Microbiol. 31:1993. doi: 10.3389/fmicb.2020.01993

PubMed Abstract | CrossRef Full Text | Google Scholar

Da Cunha, V., Davies, M. R., Douarre, P.-E., Rosinski-Chupin, I., Margarit, I., Spinali, S., et al. (2014). Streptococcus agalactiae clones infecting humans were selected and fixed through the extensive use of tetracycline. Nat. Comm. 4:4544. doi: 10.1038/ncomms5544

PubMed Abstract | CrossRef Full Text | Google Scholar

Feil, E. J., Li, B. C., Aanensen, D. M., Hanage, W. P., and Spratt, B. G. (2004). eBURST: Inferring patterns of evolutionary descent among clusters of related bacterial genotypes from Multilocus Sequence Typing Data. J. Bacteriol. 186, 1518–1530. doi: 10.1128/JB.186.5.1518-1530.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Gajic, I., Plainvert, C., Kekic, D., Dmytruk, N., Mijac, V., Tazi, A., et al. (2019). Molecular epidemiology of invasive and non-invasive group B Streptococcus circulating in Serbia. Int. J. Med. Microbiol. 309, 19–25. doi: 10.1016/j.ijmm.2018.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Glaser, P., Rusniok, C., Buchrieser, C., Chevalier, F., Frangeul, L., Msadek, T., et al. (2002). Genome sequence of Streptococcus agalactiae, a pathogen causing invasive neonatal disease. Mol. Microbiol. 45, 1499–1513. doi: 10.1046/j.1365-2958.2002.03126.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haguenoer, E., Baty, G., Pourcel, C., Lartigue, M.-F., Domelier, A.-S., Rosenau, A., et al. (2011). A multi locus variable number of tandem repeat analysis (MLVA) scheme for Streptococcus agalactiae genotyping. BMC Microbiol. 11:171. doi: 10.1186/1471-2180-11-171

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. doi: 10.1128/JCM.26.11.2465-2466.1988

PubMed Abstract | CrossRef Full Text | Google Scholar

Imperi, M., Pataracchia, M., Alfarone, G., Baldassarri, L., Orefici, G., and Creti, R. (2010). A multiplex PCR assay for the direct identification of the capsular type (Ia to IX) of Streptococcus agalactiae. J. Microbiol. Methods 80, 212–214. doi: 10.1016/j.mimet.2009.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, N., Bohnsack, J. F., Takahashi, S., Oliver, K. A., Chan, M.-S., Kunst, F., et al. (2003). Multilocus sequence typing system for group B Streptococcus. J. Clin. Microbiol. 41, 2530–2536. doi: 10.1128/jcm.41.6.2530-2536.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kayansamruaj, P., Soontara, C., Unajak, S., Dong, H. T., Rodkhum, C., Kondo, H., et al. (2019). Comparative genomics inferred two distinct populations of piscine pathogenic Streptococcus agalactiae, serotype Ia ST7 and serotype III ST283, in Thailand and Vietnam. Genomics 111, 1657–1667. doi: 10.1016/j.ygeno.2018.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalil, M. R., Uldbjerg, N., Thorsen, P. B., and Møller, J. K. (2017). Intrapartum PCR assay versus antepartum culture for assessment of vaginal carriage of group B Streptococci in a Danish cohort at birth. PLoS One 12:0180262. doi: 10.1371/journal.pone.0180262

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwatra, G., Cunnington, M. C., Merrall, E., Adrian, P. V., Ip, M., Klugman, K. P., et al. (2016). Prevalence of maternal colonisation with group B Streptococcus: a systematic review and meta-analysis. Lancet Infect. Dis. 16, 1076–1084. doi: 10.1016/S1473-3099(16)30055-X

CrossRef Full Text | Google Scholar

Lamagni, T. L., Keshishian, C., Efstratiou, A., Guy, R., Henderson, K. L., Broughton, K., et al. (2013). Emerging trends in the epidemiology of invasive Group B Streptococcal disease in England and Wales, 1991-2010. Clin. Infect. Dis. 57, 682–688. doi: 10.1093/cid/cit337

PubMed Abstract | CrossRef Full Text | Google Scholar

Lier, C., Baticle, E., Horvath, P., Haguenoer, E., Valentin, A.-S., Glaser, P., et al. (2015). Analysis of the type II-A CRISPR-Cas system of Streptococcus agalactiae reveals distinctive features according to genetic lineages. Front. Genet. 6:214. doi: 10.3389/fgene.2015.00214

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G., Zhu, J., Chen, K., Gao, T., Yao, H., Liu, Y., et al. (2016). Development of Streptococcus agalactiae vaccines for tilapia. Dis. Aquat. Org. 122, 163–170. doi: 10.3354/dao03084

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez−Sanchez, M.-J., Sauvage, E., Cunha, V. D., Clermont, D., Hariniaina, E. R., Gonzalez−Zorn, B., et al. (2012). The highly dynamic CRISPR1 system of Streptococcus agalactiae controls the diversity of its mobilome. Mole. Microb. 85, 1057–1071. doi: 10.1111/j.1365-2958.2012.08172.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Manning, S. D., Springman, A. C., Lehotzky, E., Lewis, M. A., Whittam, T. S., and Davies, H. D. (2009). Multilocus Sequence Types associated with neonatal group B Streptococcal sepsis and meningitis in Canada. J. Clin. Microb. 47, 1143–1148. doi: 10.1128/JCM.01424-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Meehan, M., Cunney, R., and Cafferkey, M. (2014). Molecular epidemiology of group B streptococci in Ireland reveals a diverse population with evidence of capsular switching. Eur. J. Clin. Microbiol. Infect. Dis. 33, 1155–1162. doi: 10.1007/s10096-014-2055-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Perme, T., Golparian, D., Bombek Ihan, M., Rojnik, A., Luèovnik, M., Kornhauser Cerar, L., et al. (2020). Genomic and phenotypic characterisation of invasive neonatal and colonising group B Streptococcus isolates from Slovenia, 2001–2018. BMC Infect. Dis. 20:5599. doi: 10.1186/s12879-020-05599-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Phares, C. R., Lynfield, R., Farley, M. M., Mohle-Boetani, J., Harrison, L. H., Petit, S., et al. (2008). Epidemiology of invasive group B Streptococcal disease in the United States, 1999-2005. JAMA 299, 2056–2065. doi: 10.1001/jama.299.17.2056

PubMed Abstract | CrossRef Full Text | Google Scholar

Pillai, P., Srinivasan, U., Zhang, L., Borchardt, S. M., Debusscher, J., Marrs, C. F., et al. (2009). Streptococcus agalactiae pulsed-field gel electrophoresis patterns cross capsular types. Epidemiol. Infect. 137:1420. doi: 10.1017/S0950268809002167

PubMed Abstract | CrossRef Full Text | Google Scholar

Radtke, A., Lindstedt, B.-A., Afset, J. E., and Bergh, K. (2010). Rapid Multiple-Locus Variant-Repeat Assay (MLVA) for Genotyping of Streptococcus agalactiae. J. Clin. Microbiol. 48, 2502–2508. doi: 10.1128/JCM.00234-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Regan, J. A., Klebanoff, M. A., and Nugent, R. P. (1991). The epidemiology of group B Streptococcal colonization in pregnancy. Vaginal Infections and Prematurity Study Group. Obstet. Gynecol. 77, 604–610.

Google Scholar

Sabat, A. J., Budimir, A., Nashev, D., Sá-Leão, R., Van Dijl, J. M., Laurent, F., et al. (2013). Overview of molecular typing methods for outbreak detection and epidemiological surveillance. Euro Surveill 18, 20380.

Google Scholar

Schrag, S., Gorwitz, R., Fultz-Butts, K., and Schuchat, A. (2002). Prevention of perinatal group B Streptococcal disease. Revised guidelines from CDC. MMWR Recomm. Rep. 51, 1–22.

Google Scholar

Schuchat, A. (1998). Epidemiology of group B Streptococcal disease in the United States: shifting paradigms. Clin. Microbiol. Rev. 11, 497–513. doi: 10.1128/cmr.11.3.497

PubMed Abstract | CrossRef Full Text | Google Scholar

Schürch, A. C., Arredondo-Alonso, S., Willems, R. J. L., and Goering, R. V. (2018). Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene–based approaches. Clin. Microbiol. Infect. 24, 350–354. doi: 10.1016/j.cmi.2017.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Skoff, T. H., Farley, M. M., Petit, S., Craig, A. S., Schaffner, W., Gershman, K., et al. (2009). Increasing burden of invasive Group B Streptococcal disease in nonpregnant adults, 1990–2007. Clin. Infect. Dis. 49, 85–92. doi: 10.1086/599369

PubMed Abstract | CrossRef Full Text | Google Scholar

Slotved, H.-C., and Hoffmann, S. (2020). The epidemiology of invasive Group B Streptococcus in Denmark From 2005 to 2018. Front. Public Health 8:40. doi: 10.3389/fpubh.2020.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Slotved, H.-C., Kong, F., Lambertsen, L., Sauer, S., and Gilbert, G. L. (2007). Serotype IX, a Proposed New Streptococcus agalactiae Serotype. J. Clin. Microbiol. 45, 2929–2936. doi: 10.1128/JCM.00117-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, S., Lin, Y., Foo, K., Koh, H. F., Tow, C., Zhang, Y., et al. (2016). Group B Streptococcus Serotype III Sequence Type 283 Bacteremia Associated with Consumption of Raw Fish. Emerg. Infect. Dis. 22, 1970–1973. doi: 10.3201/eid2211.160210

PubMed Abstract | CrossRef Full Text | Google Scholar

Tazi, A., Morand, P. C., Réglier-Poupet, H., Dmytruk, N., Billoët, A., Antona, D., et al. (2011). Invasive group B Streptococcal infections in adults, France (2007–2010). Clin. Microbiol. Infect. 17, 1587–1589. doi: 10.1111/j.1469-0691.2011.03628.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Teatero, S., Athey, T. B. T., Van Caeseele, P., Horsman, G., Alexander, D. C., Melano, R. G., et al. (2015). Emergence of Serotype IV Group B Streptococcus adult invasive disease in Manitoba and Saskatchewan, Canada, Is Driven by Clonal Sequence Type 459 Strains. J. Clin. Microbiol. 53, 2919–2926. doi: 10.1128/JCM.01128-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Teatero, S., McGeer, A., Low, D. E., Li, A., Demczuk, W., Martin, I., et al. (2014). Characterization of invasive Group B Streptococcus strains from the greater Toronto Area, Canada. J. Clin. Microbiol. 52, 1441–1447. doi: 10.1128/JCM.03554-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Tettelin, H., Masignani, V., Cieslewicz, M. J., Donati, C., Medini, D., Ward, N. L., et al. (2005). Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial “pan-genome.”. Proc. Natl. Acad. Sci. USA 102, 13950–13955. doi: 10.1073/pnas.0506758102

PubMed Abstract | CrossRef Full Text | Google Scholar

Van der Mee-Marquet, N., Domelier, A.-S., Mereghetti, L., Lanotte, P., Rosenau, A., van Leeuwen, W., et al. (2006). Prophagic DNA Fragments in Streptococcus agalactiae Strains and Association with Neonatal Meningitis. J. Clin. Microbiol. 44, 1049–1058. doi: 10.1128/JCM.44.3.1049-1058.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuillemin, X., Hays, C., Plainvert, C., Dmytruk, N., Louis, M., Touak, G., et al. (2021). Invasive group B Streptococcus infections in non-pregnant adults: a retrospective study, France, 2007–2019. Clin. Microbiol. Infect. 27, .e1–.e129. doi: 10.1016/j.cmi.2020.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2021). https://www.who.int/immunization/documents/research/who_ivb_17.09/en/ (accessed April 29, 2021).

Google Scholar