Skip to main content


We’d like to understand how you use our websites in order to improve them. Register your interest.

    We’re sorry, something doesn't seem to be working properly.

    Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.

  • Published:

A multi locus variable number of tandem repeat analysis (MLVA) scheme for Streptococcus agalactiae genotyping



Multilocus sequence typing (MLST) is currently the reference method for genotyping Streptococcus agalactiae strains, the leading cause of infectious disease in newborns and a major cause of disease in immunocompromised children and adults. We describe here a genotyping method based on multiple locus variable number of tandem repeat (VNTR) analysis (MLVA) applied to a population of S. agalactiae strains of various origins characterized by MLST and serotyping.


We studied a collection of 186 strains isolated from humans and cattle and three reference strains (A909, NEM316 and 2603 V/R). Among 34 VNTRs, 6 polymorphic VNTRs loci were selected for use in genotyping of the bacterial population. The MLVA profile consists of a series of allele numbers, corresponding to the number of repeats at each VNTR locus. 98 MLVA genotypes were obtained compared to 51 sequences types generated by MLST. The MLVA scheme generated clusters which corresponded well to the main clonal complexes obtained by MLST. However it provided a higher discriminatory power. The diversity index obtained with MLVA was 0.960 compared to 0.881 with MLST for this population of strains.


The MLVA scheme proposed here is a rapid, cheap and easy genotyping method generating results suitable for exchange and comparison between different laboratories and for the epidemiologic surveillance of S. agalactiae and analyses of outbreaks.


Streptococcus agalactiae, one of the group B streptococci (GBS), is a leading cause of bovine mastitis [1] and has been implicated in cases of invasive disease in humans since the 1960s and 1970s [2]. GBS have emerged as major pathogens in neonates [3] and in elderly adults, in whom they cause invasive infections, such as meningitis, soft tissue infections, endocarditis and osteoarticular infections [4, 5]. There is a considerable body of evidence to suggest a genetic link between bovine isolates and the emerging human isolates [6, 7].

GBS isolates were initially distinguished on the basis of differences in capsule polysaccharides, giving rise to 10 different serotypes [8, 9]. Serotype III has been identified as a marker of late-onset neonatal disease isolates [10], but serotyping does not have sufficient discriminatory power to distinguish between isolates. Molecular methods have therefore been developed to determine the genetic relationships between isolates: multilocus enzyme electrophoresis [11], ribotyping [12], random amplified polymorphism DNA (RAPD) [13, 14] and pulsed-field gel electrophoresis (PFGE) [15]. These methods make it possible to compare isolates and to define particular bacterial genogroups associated with invasive isolates in neonates. These findings were confirmed by multilocus sequence typing, as described by Jones et al. [16]. Other studies have shown that sequence type 17 (ST-17) isolates are associated with invasive behavior [17, 18]. Two methods are currently used to explore the genetic links between isolates: PFGE for epidemiological studies, and MLST for both epidemiological and phylogenetic studies.

Analyses of fully sequenced bacterial genomes have revealed the existence of tandemly repeated sequences varying in size, location and the type of repetition [19]. Tandem repeats (TR) consist of a direct repetition of between one and more than 200 nucleotides, which may or may not be perfectly identical, located within or between genes. Depending on the size of the unit, the TR may be defined as a microsatellite (up to 9 bp) or a minisatellite (more than 9 bp) [19]. A fraction of these repeated sequences display intraspecies polymorphism and are described as VNTRs (variable number of tandem repeats). The proportion of VNTRs in the genome varies between bacterial species. Indeed, variation in the number of repeats at particular loci is used by some bacteria as a means of rapid genomic and phenotypic adaptation to the environment [20].

A molecular typing method based on VNTRs variability has recently been developed and applied to the typing of several bacterial pathogens [19]. Multiple locus VNTR analysis, or MLVA, is a PCR-based method that was originally developed for the typing of Haemophilus influenzae [21], Mycobacterium tuberculosis [22] and two bacterial species with potential for use in bioterrorism, Bacillus anthracis and Yersinia pestis [23, 24]. This method has since been shown to be useful for the genotyping of several other bacterial species causing disease in humans, including Streptococcus pneumoniae [25], Legionella pneumophila [26], Brucella [27, 28], Pseudomonas aeruginosa [29] and Staphylococcus aureus [30]. This technique has several advantages. For example, in bacterial species with high levels of genetic diversity, the study of six to eight markers is sufficient for accurate discrimination between strains [26]. Highly monomorphic species, such as B. anthracis, may be typed by MLVA, but this requires the use of a larger number of markers (25 VNTRs for B. anthracis) [31]. The discriminatory power of MLVA may also be increased by adding extra panels of more polymorphic markers [28] or by sequencing repeated sequences displaying internal variability [26]. Conversely, the evaluation of differences in the number of repeats only, on the basis of MLVA, is a cheap and rapid method that is not technically demanding. The work of Radtke et al. showed relevance of MLVA for S. agalactiae genotyping [32].

Our aim in this study was to develop a MLVA scheme for the genotyping of a population of S. agalactiae strains of various origins previously characterized by MLST.



Our collection consisted of 186 epidemiologically unrelated S. agalactiae strains, isolated from humans and cattle between 1966 and 2004 in France. Five of the 152 human strains were isolated from the gastric fluid of neonates, 71 were isolated from cases of vaginal carriage, 59 were isolated from cerebrospinal fluid and 17 were isolated from cultures of blood from adults presenting confirmed endocarditis according to the modified Duke criteria [33]. The 34 bovine strains were isolated from cattle presenting clinical signs of mastitis. We also studied three reference strains: NEM316, A909 and 2603 V/R. Each strain had previously been identified on the basis of Gram-staining, colony morphology, beta-hemolysis and Lancefield group antigen determination (Slidex Strepto Kit®, bioMérieux, Marcy l'Etoile, France). The capsular serotype was identified with the Pastorex® rapid latex agglutination test (Bio-Rad, Hercules, USA) and by molecular serotyping, as described by Manning et al. [34]. We were unable to determine the serotype for 20 strains.

DNA extraction

The bacteria were lysed mechanically with glass beads and their genomic DNA was extracted with an Invisorb® Spin Cell Mini kit (Invitek, Berlin, Germany).

MLST and assignment to clonal clusters

MLST was carried out as previously described [16]. Briefly, PCR was used to amplify small (≈ 500 bp) fragments from seven housekeeping genes (adhP, pheS, atr, glnA, sdhA, glcK and tkt) chosen on the basis of their chromosomal location and sequence diversity. The seven PCR products were purified and sequenced and an allele number was assigned to each fragment on the basis of its sequence. A sequence type (ST), based on the allelic profile of the seven amplicons, was assigned to each strain. The sequences of all new alleles and the composition of the new STs identified are available from Strains were grouped into clonal complexes (CCs) with eBURST software [35]. An eBURST clonal complex (CC) was defined as all allelic profiles sharing six identical alleles with at least one other member of the group. The term "singleton ST" refers to a ST that did not cluster into a CC.

Identification of VNTR loci

Tandem repeats were identified in the sequenced genomes of the three reference strains, NEM316, A909 and 2603 V/R, with the Microbial Tandem Repeats Database[36] and the Tandem Repeats Finder program [37]. We determined the size of the repeat sequence and the number of repeat units for the three reference strains. BLAST analysis was carried out to determine whether the repeats were located within or between genes and to identify a hypothetical function for the open reading frame involved. The TR locus name was defined according to the following nomenclature: common name_size of the repeat sequence_size of the amplicon for the reference strain_corresponding number of repeats (Table 1). The primers used for amplification targeted the 5' and 3' flanking regions of selected loci and matched the sequences present at these positions in the genomes of strains NEM316, A909 and 2603 V/R. We initially selected and evaluated 34 tandem repeats with repeat units of more than 9 bp in length. Some TRs were not present in all the strains, some were present in all strains and displayed no polymorphism, and others were too large for amplification in standard conditions. Six TRs were retained for this study, selected on the basis of their greater stability and discriminatory power for four of the six (Table 1).

Table 1 Characteristics of the 6 VNTR loci selected for MLVA scheme to genotype the 186 strains of S. agalactiae

Multiple locus VNTR analysis (MLVA)

The primers used for the VNTRs amplification are presented in Table 2. Three loci have already been described by Radtke et al. in a contemporary study but were amplified here with other primers [32] (Table 2). For the SAG7 locus, no amplification was observed with primers directly flanking the TR for 14% (26/189) of the strains. A second primer pair targeting larger consensual flanking regions was designed to confirm the absence of the locus. PCR was performed in a final volume of 25 μl containing 10 ng DNA, 1 × PCR Reaction Buffer, 2 mM MgCl2 (Applied Biosystems), 5% DMSO (dimethyl sulfoxide), 1 unit of Taq DNA polymerase (Applied Biosystems), 200 μM of each dNTP and 0.5 μM of each flanking primer (Eurogentec, Belgium). Amplification was performed in a 2720 Thermal Cycler (Applied Biosystems) under the following conditions: initial denaturation for 5 min at 94°C, followed by 30 cycles of denaturation for 30 s at 94°C, annealing for 30 s at 50°C and elongation for 60 s at 72°C plus a final elongation step for 7 min at 72°C. We separated 10 μl of PCR product by electrophoresis in a 2% agarose gel (Eurogentec, Belgium), which was also loaded with a 100 bp DNA size ladder (New England BioLabs). Electrophoresis was performed in 20 cm-long gels, in 1× TBE buffer (89 mM Tris-Borate, 2.5 mM EDTA) containing 1 μg/ml ethidium bromide run at 10 V/cm. In each run, at least one lane was loaded with PCR product from one of the reference strains, NEM316, A909 or 2603 V/R. The gels were photographed under ultraviolet illumination, with Vision-Capt® Software (Vilber-Lourmat, Marne la Vallée, France). The number of repeats for each VNTR was deduced from amplicon size, by comparison with the reference strain, for which the number of repeats was known. The allele number corresponded to the number of repeats. For the SAG7 locus, the lack of a VNTR was revealed by the absence of amplification with the first primer pair and the amplification of a fragment of the expected size with the second primer pair, which targeted larger consensual flanking regions. In this case, an allele number of 0 was given. For the SAG21 locus, a 117 bp PCR product was obtained, demonstrating deletion of the inserted sequence and, thus, the absence of a VNTR. An allele number of 0 was also assigned in this case. The MLVA genotype of a strain was expressed as its allelic profile, corresponding to the number of repeats at each VNTR, listed in the order SAG2, SAG3, SAG4, SAG7, SAG21, SAG22.

Table 2 Primers used in the MLVA scheme

Data analysis

The polymorphism index of individual or combined VNTR loci was calculated with the Hunter-Gaston diversity index [38], an application of Simpson's index of diversity [39]. Confidence intervals (CI) were calculated as described by Grundmann et al. [40]. The categorical coefficient (also called Hamming's distance) and unweighted pair group method with arithmetic mean (UPGMA) clustering approaches were run within BioNumerics. A cutoff value of 50% similarity was applied to define MLVA clusters. The minimum spanning tree (MST) was generated with BioNumerics. Each circle represents an MLVA genotype and its size is proportional to the number of strains. A logarithmic scale was used when drawing branches. The thicker branches link the MLVA genotypes differing by only one allele, the thinner branches link MLVA genotypes differing by more than one allele.


MLST genotyping

MLST was performed on the 189 S. agalactiae strains, identifying a total of 51 individual STs. Eburst analysis clustered the STs into five clonal complexes (CC17, CC19, CC10, CC23 and CC7), two groups with only two STs and six singletons (Table 3). Two of the CCs -- CC17 (73 strains) and CC19 (63 strains) -- accounted for 72% (136/189) of the strains. CC23 accounted for 8% (15/189) of the strains. The various serotypes of S. agalactiae were distributed between multiple CCs and singleton STs. STs were characterized by a predominant serotype: serotype V in ST-1, serotype III in ST-17 and ST-19, serotype Ib in ST-10 and ST-12. ST-23 contained two serotypes (serotype Ia and III; Table 3). The population was therefore representative of S. agalactiae diversity in terms of anatomic origin, serotypes and clonal complexes (Table 3).

Table 3 Distribution of the 186 S. agalactiae strains studied and the 3 reference strains (NEM316, A909 and 2603 V/R), as a function of serotype and origin, within MLST clonal complexes

Description of the MLVA scheme

The six VNTRs were amplified from all 189 strains. MLVA was carried out with individual PCRs and agarose gel electrophoresis of the amplicons, as shown in Figure 1, for a subset of VNTRs. The repeat unit size of the six VNTRs was between 18 bp and 159 bp, making it straightforward to estimate the size of amplicons on agarose gels. For SAG2, SAG3, SAG4 and SAG7, amplicons were between 114 and 573 bp in size and were readily resolved by 2% agarose gel electrophoresis (Table 1). For SAG21 (48 bp repeat unit) and SAG22 (159 bp repeat unit), few amplicons exceeded 1,000 bp and extensive electrophoretic separation was required for precise estimations of size. For SAG21, three strains gave rise to amplicons of more than 1500 bp in size. This made it difficult to assess the number of repeats with any degree of precision, and an arbitrary allele number of > 30 was assigned in these cases. For SAG7, no amplification with the first primer pair was observed for 14% of strains. This locus is part of a genomic island and a second primer pair targeting consensual flanking regions beyond the borders of this genomic island was designed to confirm the deletion of the VNTR locus. The number of alleles was between two for SAG3 and 26 for SAG21. Thus, this MLVA method combined markers with a low discriminatory power (Hunter and Gaston's index of diversity or HGDI < 0.5) with highly discriminant markers, such as SAG21. With the exception of SAG2, the VNTRs used in this MLVA method were located within open reading frames (Table 1). SAG2 is located upstream from the gene encoding the ribosomal protein S10; SAG3 is located within dnaJ, encoding a co-chaperone protein (Hsp40). SAG21 is located within fbsA, encoding a protein involved in adhesion. SAG4, SAG7 and SAG22 are located in a "predicted coding region" of unknown function.

Figure 1

Polymorphism of four VNTRs. The polymorphism of VNTRs (SAG2, SAG3, SAG4 and SAG22) is shown by agarose gel electrophoresis of PCR products. The first strain on each gel is the reference strain and the PCR products were loaded alongside a 100 bp DNA size ladder (the sizes in base pairs are shown on the left side of the first panel). The allele number, corresponding to the number of repeats, is indicated under the band.

MLVA genotyping and clustering

The MLVA scheme resolved 98 genotypes among the 189 strains (Table 4). Five MLVA genotypes were represented by more than five strains: genotype 46 (n = 32), genotype 47 (n = 13), genotype 33 (n = 11), genotype 57 (n = 7) and genotype 51 (n = 6). Seventy-five MLVA genotypes were represented by only one strain (Table 4). S. agalactiae strains of different origins were spread among a number of MLVA genotypes. However 66% (39/59) of the strains isolated from cerebrospinal fluid were confined to four MLVA genotypes (genotypes 46, 47, 51 and 57). An MLVA cluster was defined by a cutoff value of 50% similarity with the UPGMA algorithm (Figures 2 and 3). Nine MLVA clusters, each containing more than four strains, were identified (MLVA clusters 1 to 9) (Figures 2 and 3 and Figure 4A). All clusters other than cluster 1 were congruent with the two algorithms, UPGMA and MST.

Table 4 MLVA genotypes resolved by the MLVA-6 scheme
Figure 2

MLVA clustering of the 189 strains of S. agalactiae by the UPGMA method, run in BioNumerics. The names of strains (Strain Id), MLST clonal complex (CC), MLST sequence type (ST), MLST profile, serotype and the origin of strains are shown on the right. A cutoff value of 50% similarity was applied to define MLVA clusters (named MLVA cluster 1 to MLVA cluster 9). The colors used are based on MLVA clusters.

Figure 3

MLVA clustering of the 189 strains of S. agalactiae by the UPGMA method, run in BioNumerics. The names of strains (Strain Id), MLST clonal complex (CC), MLST sequence type (ST), MLST profile, serotype and the origin of strains are shown on the right. A cutoff value of 50% similarity was applied to define MLVA clusters (named MLVA cluster 1 to MLVA cluster 9). The colors used are based on MLVA clusters.

Figure 4

Minimum spanning tree (MST) representation of the MLVA clustering. The colors used in figure 4A are based on MLVA clusters. The colors used in figure 4B are based on MLST clonal complexes. White circles correspond to genotypes not clustered by MLVA or MLST. The MLVA data for 189 strains, including 3 reference strains, were analyzed in BioNumerics. Each circle represents an MLVA genotype and its size is proportional to the number of strains. A logarithmic scale was used when drawing branches. The thicker branches link the MLVA genotypes differing by only one allele, the thinner branches link MLVA genotypes differing by more than one allele.

Comparison of MLVA and MLST clustering

MLVA clustering showed a clonal distribution of the population similar to that obtained by MLST (Figure 4). All human strains of MLST CC17 clustered together in MLVA cluster 9 and the bovine strains of MLST CC17 belonged to several MLVA clusters, suggesting greater heterogeneity of this population (Figure 4). With the exception of 3 strains, the MLST CC19 strains clustered into 2 linked MLVA clusters, MLVA cluster 6 and MLVA cluster 7. The MLST CC23 strains of serotype III and the MLST CC10 strains clustered into MLVA cluster 2. The strains from MLST CC23 serotype Ia also formed a separate group, the MLVA cluster 8.

Discrimination of S. agalactiaestrains by MLVA

The diversity index obtained with MLVA was 0.960 (95% CI [0.943 - 0.978]), which is greater than that obtained with MLST (0.881). For the population studied, MLVA distinguished 98 genotypes, whereas MLST distinguished 51 different STs. A much higher level of diversity was observed with MLVA, particularly within the major CCs. For example, the 73 CC17 strains were separated into 12 STs by MLST and 22 MLVA genotypes; the 63 CC19 strains were separated into 15 STs by MLST and 35 MLVA genotypes and the 15 CC23 strains were separated into 6 STs by MLST and 15 MLVA genotypes. Nevertheless, two genotypes (46 and 47) accounted for 76% (45/59) of CC17 strains of human origin. For this particular genogroup, the discriminatory power of the MLVA method was greater than that of MLST, although it remained low.


In this study, we applied the multi locus VNTR analysis (MLVA) typing method to S. agalactiae. VNTR analysis, a method based on tandem repeat polymorphisms at multiple loci, has been successfully applied to many other bacterial species [30, 41]. We investigated the relevance of this tool for the genotyping of S. agalactiae, by testing this method on six VNTR loci in 189 strains previously characterized by MLST and serotyping. The MLVA-6 scheme is inexpensive and can be carried out with the equipment routinely used for PCR amplification and agarose gel electrophoresis. For the six VNTR loci, amplification was achieved with all the strains tested. For SAG7, a second PCR targeting a larger flanking region was required for 14% of the strains, which did not have a 16 kb genomic island encompassing the VNTR. The repeat sizes of the six VNTRs were sufficiently large for evaluation of the number of repeats on agarose gels. Moreover, the conversion of results into allelic profiles should make it possible to construct databases for exchange between laboratories. The MLVA-6 scheme includes a set of markers with different diversity indices, making it suitable for epidemiological studies. Markers with a moderate diversity and small number of alleles (presumably reflecting their slow rate of evolution) define clusters, whereas markers displaying more rapid evolution reflect variability within clusters. The MLVA-6 method described here is a rapid, reproducible and epidemiologically meaningful typing tool.

Three loci studied in the present MLVA scheme are in common with the MLVA scheme proposed by Radtke et al. [32]. The 3 additional loci studied here provide more weight to clusters while maintaining a high discrimination power. Moreover, in the MLVA scheme proposed here, only one locus (SAG7) was missing in some strains (14%), and another primer pair targeting larger consensual flanking region confirmed the absence of this locus with a specific amplification. Unlike Radtke et al., we sought to develop a MLVA scheme in which a PCR product was amplified in all strains whether the VNTR was present or absent. In fact, negative amplification may result from the lack of a VNTR locus or modification of the flanking regions, especially as some VNTRs are close to transposases or insertion sequences such as SAG4 (alias SATR1) which is close to IS1381. Thus, the possibility of negative amplification for 3 out of 5 VNTR loci in the Radtke et al. MLVA analysis could be a real problem in terms of resolution and reproducibility of the genotyping method. Nevertheless, cumulative works allow to define the best set of VNTR loci, as has already been done for other bacterial species such as Mycobacterium tuberculosis [22, 4246] and Staphylococcus aureus [30, 4749]. Finally, the study of 34 isolates of bovine origin provided information about their distribution, especially those belonging to MLST CC17.

Population analysis by MLVA revealed a clonal distribution of the strains similar to that obtained by MLST. The greater discriminatory index of MLVA (0.96) made it possible to distinguish between strains within the clonal complexes defined by MLST. Thus, MLVA divided CC23 into two groups: one associated with serotype III and the other associated with serotype Ia. Moreover, MLVA also separated CC17 into two groups: one corresponding to strains of human origin and the other, containing several related STs (ST-61, ST-64, ST-301 etc.), corresponding to strains of animal origin only. A previous study analyzing 75 strains of S. agalactiae of human and animal origin by whole-genome DNA-array hybridization also separated ST-23 strains into two clusters, one associated with serotype III and the other with serotype Ia [50]. Each of these two clusters was associated with a particular pattern of surface protein expression. This previous study also separated the bovine and human CC17 strains [50]. These results are consistent with an ancient divergence of these clusters, whereas other observations based on MLST analysis suggest that ST-17 strains may have arisen from a bovine ancestor [6]. The lack of a strict correlation between the results of MLST and MLVA may be accounted for by differences in the markers used for MLST (targeting housekeeping genes) and MLVA (targeting a set of diverse regions that may or may not be conserved). Unlike MLST, MLVA targets several types of markers: genes involved in metabolism, genes associated with virulence and a genomic island. Indeed, SAG2 is located upstream from the gene encoding the ribosomal protein S10 which is involved in transcription and translation, and SAG3 is located within dnaJ, which encodes a member of the Hsp70 family, a co-chaperone protein (Hsp40). The SAG21 locus encodes a surface protein involved in virulence, FbsA. The SAG7 locus is located on a genomic island and belongs to a gene encoding a hypothetical protein whose function has not yet been identified, like most of the genes of genomic islands [51]. Clustering based on MLVA data was almost identical with the UPGMA and MST algorithms except for cluster 1. The differences in mathematical calculation between the two methods may account for the observed differences in strain clustering. This phenomenom has been previously observed in MLVA studies [52].

Some VNTRs for the alpha C protein have already been described in S. agalactiae [41, 53, 54]. One of these VNTRs is involved in regulating gene expression: a pentanucleotide repeat located upstream from the promoter regulates expression in vitro by phase variation. Another is an intragenic VNTR that modifies the size of the alpha C protein, thereby altering its antigenicity and strain virulence [53]. These two VNTR loci were not included in the MLVA method proposed here, in one case because the small size of the repeat unit (5 bp) complicates the mode of PCR fragment size assessment [19]. The amplicons of the second VNTR locus not included were more than 2000 bp in size, again making it difficult to evaluate repeat number. Tandem repeats were also found in the gene encoding another surface protein, FbsA, which interacts with epithelial cells and is involved in invasion of the central nervous system of colonized neonates. Its ability to bind to fibrinogen depends on the number of repeats of a unit of 16 amino acids present at its N-terminus [55]. A particular number of repeats is associated with the greater potential of the ST-17 strains implicated in neonatal meningitis to adhere to fibrinogen [56]. This major marker was included in our MLVA method and corresponds to SAG21.


The MLVA method proposed here is a simple genotyping method producing results that can be exchanged between laboratories. MLVA generated major clusters that corresponded well to the main clonal complexes obtained by MLST. However its discriminatory power provided was greater that that of MLST. MLVA could also therefore be used as an epidemiological tool, given its high discriminatory power, making it possible to distinguish between strains of homogenous lineages. The specificities of the VNTRs for each phylogenetic lineage raise questions about the role of VNTRs in the adaptation of S. agalactiae to its environment and in virulence. Further studies are required to clarify these issues.


  1. 1.

    Keefe GP: Streptococcus agalactiae mastitis: a review. Can Vet J. 1997, 38: 429-437.

  2. 2.

    Schuchat A: Group B streptococcal disease: from trials and tribulations to triumph and trepidation. Clin Infect Dis. 2001, 33: 751-756. 10.1086/322697.

  3. 3.

    Bohnsack JF, Whiting A, Gottschalk M, Dunn DM, Weiss R, Azimi PH, Philips JB, Weisman LE, Rhoads GG, Lin F-YC: Population structure of invasive and colonizing strains of Streptococcus agalactiae from neonates of six U.S. Academic Centers from 1995 to 1999. J Clin Microbiol. 2008, 46: 1285-1291. 10.1128/JCM.02105-07.

  4. 4.

    Edwards MS, Rench MA, Palazzi DL, Baker CJ: Group B streptococcal colonization and serotype-specific immunity in healthy elderly persons. Clin Infect Dis. 2005, 40: 352-357. 10.1086/426820.

  5. 5.

    Farley MM: Group B streptococcal disease in nonpregnant adults. Clin Infect Dis. 2001, 33: 556-561. 10.1086/322696.

  6. 6.

    Bisharat N, Crook DW, Leigh J, Harding RM, Ward PN, Coffey TJ, Maiden MC, Peto T, Jones N: Hyperinvasive neonatal group B streptococcus has arisen from a bovine ancestor. J Clin Microbiol. 2004, 42: 2161-2167. 10.1128/JCM.42.5.2161-2167.2004.

  7. 7.

    Héry-Arnaud G, Bruant G, Lanotte P, Brun S, Picard B, Rosenau A, van der Mee-Marquet N, Rainard P, Quentin R, Mereghetti L: Mobile genetic elements provide evidence for a bovine origin of clonal complex 17 of Streptococcus agalactiae. Appl Environ Microbiol. 2007, 73: 4668-4672. 10.1128/AEM.02604-06.

  8. 8.

    Lindahl G, Stålhammar-Carlemalm M, Areschoug T: Surface proteins of Streptococcus agalactiae and related proteins in other bacterial pathogens. Clin Microbiol Rev. 2005, 18: 102-127. 10.1128/CMR.18.1.102-127.2005.

  9. 9.

    Slotved H-C, Kong F, Lambertsen L, Sauer S, Gilbert GL: Serotype IX, a proposed new Streptococcus agalactiae serotype. J Clin Microbiol. 2007, 45: 2929-2936. 10.1128/JCM.00117-07.

  10. 10.

    Musser JM, Mattingly SJ, Quentin R, Goudeau A, Selander RK: Identification of a high-virulence clone of type III Streptococcus agalactiae (group B Streptococcus) causing invasive neonatal disease. Proc Natl Acad Sci USA. 1989, 86: 4731-4735. 10.1073/pnas.86.12.4731.

  11. 11.

    Quentin R, Huet H, Wang FS, Geslin P, Goudeau A, Selander RK: Characterization of Streptococcus agalactiae strains by multilocus enzyme genotype and serotype: identification of multiple virulent clone families that cause invasive neonatal disease. J Clin Microbiol. 1995, 33: 2576-2581.

  12. 12.

    Blumberg HM, Stephens DS, Licitra C, Pigott N, Facklam R, Swaminathan B, Wachsmuth IK: Molecular epidemiology of group B streptococcal infections: use of restriction endonuclease analysis of chromosomal DNA and DNA restriction fragment length polymorphisms of ribosomal RNA genes (ribotyping). J Infect Dis. 1992, 166: 574-579. 10.1093/infdis/166.3.574.

  13. 13.

    Chatellier S, Huet H, Kenzi S, Rosenau A, Geslin P, Quentin R: Genetic diversity of rRNA operons of unrelated Streptococcus agalactiae strains isolated from cerebrospinal fluid of neonates suffering from meningitis. J Clin Microbiol. 1996, 34: 2741-2747.

  14. 14.

    Chatellier S, Ramanantsoa C, Harriau P, Rolland K, Rosenau A, Quentin R: Characterization of Streptococcus agalactiae strains by randomly amplified polymorphic DNA analysis. J Clin Microbiol. 1997, 35: 2573-2579.

  15. 15.

    Rolland K, Marois C, Siquier V, Cattier B, Quentin R: Genetic features of Streptococcus agalactiae strains causing severe neonatal infections, as revealed by pulsed-field gel electrophoresis and hylB gene analysis. J Clin Microbiol. 1999, 37: 1892-1898.

  16. 16.

    Jones N, Bohnsack JF, Takahashi S, Oliver KA, Chan M-S, Kunst F, Glaser P, Rusniok C, Crook DWM, Harding RM, Bisharat N, Spratt BG: Multilocus sequence typing system for group B streptococcus. J Clin Microbiol. 2003, 41: 2530-2536. 10.1128/JCM.41.6.2530-2536.2003.

  17. 17.

    Lamy M-C, Dramsi S, Billoët A, Réglier-Poupet H, Tazi A, Raymond J, Guérin F, Couvé E, Kunst F, Glaser P, Trieu-Cuot P, Poyart C: Rapid detection of the "highly virulent" group B Streptococcus ST-17 clone. Microbes Infect. 2006, 8: 1714-1722. 10.1016/j.micinf.2006.02.008.

  18. 18.

    Luan S-L, Granlund M, Sellin M, Lagergård T, Spratt BG, Norgren M: Multilocus sequence typing of Swedish invasive group B streptococcus isolates indicates a neonatally associated genetic lineage and capsule switching. J Clin Microbiol. 2005, 43: 3727-3733. 10.1128/JCM.43.8.3727-3733.2005.

  19. 19.

    Lindstedt B-A: Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis. 2005, 26: 2567-2582. 10.1002/elps.200500096.

  20. 20.

    Martin P, van de Ven T, Mouchel N, Jeffries AC, Hood DW, Moxon ER: Experimentally revised repertoire of putative contingency loci in Neisseria meningitidis strain MC58: evidence for a novel mechanism of phase variation. Mol Microbiol. 2003, 50: 245-257. 10.1046/j.1365-2958.2003.03678.x.

  21. 21.

    Van Belkum A, Melchers WJ, Ijsseldijk C, Nohlmans L, Verbrugh H, Meis JF: Outbreak of amoxicillin-resistant Haemophilus influenzae type b: variable number of tandem repeats as novel molecular markers. J Clin Microbiol. 1997, 35: 1517-1520.

  22. 22.

    Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C: Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol. 2000, 36: 762-771.

  23. 23.

    Keim P, Price LB, Klevytska AM, Smith KL, Schupp JM, Okinaka R, Jackson PJ, Hugh-Jones ME: Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis. J Bacteriol. 2000, 182: 2928-2936. 10.1128/JB.182.10.2928-2936.2000.

  24. 24.

    Le Flèche P, Hauck Y, Onteniente L, Prieur A, Denoeud F, Ramisse V, Sylvestre P, Benson G, Ramisse F, Vergnaud G: A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. 2001, 1: 2-10.1186/1471-2180-1-2.

  25. 25.

    Koeck J-L, Njanpop-Lafourcade B-M, Cade S, Varon E, Sangare L, Valjevac S, Vergnaud G, Pourcel C: Evaluation and selection of tandem repeat loci for Streptococcus pneumoniae MLVA strain typing. BMC Microbiol. 2005, 5: 66-10.1186/1471-2180-5-66.

  26. 26.

    Pourcel C, Visca P, Afshar B, D'Arezzo S, Vergnaud G, Fry NK: Identification of variable-number tandem-repeat (VNTR) sequences in Legionella pneumophila and development of an optimized multiple-locus VNTR analysis typing scheme. J Clin Microbiol. 2007, 45: 1190-1199. 10.1128/JCM.02078-06.

  27. 27.

    Al Dahouk S, Flèche PL, Nöckler K, Jacques I, Grayon M, Scholz HC, Tomaso H, Vergnaud G, Neubauer H: Evaluation of Brucella MLVA typing for human brucellosis. J Microbiol Methods. 2007, 69: 137-145. 10.1016/j.mimet.2006.12.015.

  28. 28.

    Le Flèche P, Jacques I, Grayon M, Al Dahouk S, Bouchon P, Denoeud F, Nöckler K, Neubauer H, Guilloteau LA, Vergnaud G: Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiol. 2006, 6: 1471-1484.

  29. 29.

    Vu-Thien H, Corbineau G, Hormigos K, Fauroux B, Corvol H, Clément A, Vergnaud G, Pourcel C: Multiple-locus variable-number tandem-repeat analysis for longitudinal survey of sources of Pseudomonas aeruginosa infection in cystic fibrosis patients. J Clin Microbiol. 2007, 45: 3175-3183. 10.1128/JCM.00702-07.

  30. 30.

    Pourcel C, Hormigos K, Onteniente L, Sakwinska O, Deurenberg RH, Vergnaud G: Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value. J Clin Microbiol. 2009, 47: 3121-3128. 10.1128/JCM.00267-09.

  31. 31.

    Lista F, Faggioni G, Valjevac S, Ciammaruconi A, Vaissaire J, le Doujet C, Gorgé O, De Santis R, Carattoli A, Ciervo A, Fasanella A, Orsini F, D'Amelio R, Pourcel C, Cassone A, Vergnaud G: Genotyping of Bacillus anthracis strains based on automated capillary 25-loci multiple locus variable-number tandem repeats analysis. BMC Microbiol. 2006, 6: 33-10.1186/1471-2180-6-33.

  32. 32.

    Radtke A, Lindstedt B-A, Afset JE, Bergh K: Rapid multiple-locus variant-repeat assay (MLVA) for genotyping of Streptococcus agalactiae. J Clin Microbiol. 2010, 48: 2502-2508. 10.1128/JCM.00234-10.

  33. 33.

    Li JS, Sexton DJ, Mick N, Nettles R, Fowler VG, Ryan T, Bashore T, Corey GR: Proposed modifications to the Duke criteria for the diagnosis of infective endocarditis. Clin Infect Dis. 2000, 30: 633-638. 10.1086/313753.

  34. 34.

    Manning SD, Lacher DW, Davies HD, Foxman B, Whittam TS: DNA polymorphism and molecular subtyping of the capsular gene cluster of group B streptococcus. J Clin Microbiol. 2005, 43: 6113-6116. 10.1128/JCM.43.12.6113-6116.2005.

  35. 35.

    Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004, 186: 1518-1530. 10.1128/JB.186.5.1518-1530.2004.

  36. 36.

    Denoeud F, Vergnaud G: Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource. BMC Bioinformatics. 2004, 5: 4-10.1186/1471-2105-5-4.

  37. 37.

    Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.

  38. 38.

    Hunter PR, Gaston MA: Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol. 1988, 26: 2465-2466.

  39. 39.

    Simpson EH: Measurement of diversity. Nature. 1949, 163: 688-10.1038/163688a0.

  40. 40.

    Grundmann H, Hori S, Tanner G: Determining confidence intervals when measuring genetic diversity and the discriminatory abilities of typing methods for microorganisms. J Clin Microbiol. 2001, 39: 4190-4192. 10.1128/JCM.39.11.4190-4192.2001.

  41. 41.

    Puopolo KM, Madoff LC: Upstream short sequence repeats regulate expression of the alpha C protein of group B Streptococcus. Mol Microbiol. 2003, 50: 977-991. 10.1046/j.1365-2958.2003.03745.x.

  42. 42.

    Frothingham R, Meeker-O'Connell WA: Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. Microbiology. 1998, 144: 1189-1196. 10.1099/00221287-144-5-1189.

  43. 43.

    Supply P, Lesjean S, Savine E, Kremer K, van Soolingen D, Locht C: Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units. J Clin Microbiol. 2001, 39: 3563-3571. 10.1128/JCM.39.10.3563-3571.2001.

  44. 44.

    Mazars E, Lesjean S, Banuls AL, Gilbert M, Vincent V, Gicquel B, Tibayrenc M, Locht C, Supply P: High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology. Proc Natl Acad Sci USA. 2001, 98: 1901-1906. 10.1073/pnas.98.4.1901.

  45. 45.

    Le Flèche P, Fabre M, Denoeud F, Koeck J-L, Vergnaud G: High resolution, on-line identification of strains from the Mycobacterium tuberculosis complex based on tandem repeat typing. BMC Microbiol. 2002, 2: 37-10.1186/1471-2180-2-37.

  46. 46.

    Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rüsch-Gerdes S, Willery E, Savine E, de Haas P, van Deutekom H, Roring S, Bifani P, Kurepina N, Kreiswirth B, Sola C, Rastogi N, Vatin V, Gutierrez MC, Fauville M, Niemann S, Skuce R, Kremer K, Locht C, van Soolingen D: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006, 44: 4498-4510. 10.1128/JCM.01392-06.

  47. 47.

    Sabat A, Krzyszton-Russjan J, Strzalka W, Filipek R, Kosowska K, Hryniewicz W, Travis J, Potempa J: New method for typing Staphylococcus aureus strains: multiple-locus variable-number tandem repeat analysis of polymorphism and genetic relationships of clinical isolates. J Clin Microbiol. 2003, 41: 1801-1804. 10.1128/JCM.41.4.1801-1804.2003.

  48. 48.

    Francois P, Huyghe A, Charbonnier Y, Bento M, Herzig S, Topolski I, Fleury B, Lew D, Vaudaux P, Harbarth S, van Leeuwen W, van Belkum A, Blanc DS, Pittet D, Schrenzel J: Use of an automated multiple-locus, variable-number tandem repeat-based method for rapid and high-throughput genotyping of Staphylococcus aureus isolates. J Clin Microbiol. 2005, 43: 3346-3355. 10.1128/JCM.43.7.3346-3355.2005.

  49. 49.

    Hardy KJ, Ussery DW, Oppenheim BA, Hawkey PM: Distribution and characterization of staphylococcal interspersed repeat units (SIRUs) and potential use for strain differentiation. Microbiology. 2004, 150: 4045-4052. 10.1099/mic.0.27413-0.

  50. 50.

    Brochet M, Couvé E, Zouine M, Vallaeys T, Rusniok C, Lamy M-C, Buchrieser C, Trieu-Cuot P, Kunst F, Poyart C, Glaser P: Genomic diversity and evolution within the species Streptococcus agalactiae. Microbes Infect. 2006, 8: 1227-1243. 10.1016/j.micinf.2005.11.010.

  51. 51.

    Tettelin H, et al: Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome". Proc Natl Acad Sci USA. 2005, 102: 13950-13955. 10.1073/pnas.0506758102.

  52. 52.

    Dauchy FA, Degrange S, Charron A, Dupon M, Xin Y, Bebear C, Maugein J: Variable-number tandem-repeat markers for typing Mycobacterium intracellulare strains isolated in humans. BMC Microbiol. 2010, 10: 93-10.1186/1471-2180-10-93.

  53. 53.

    Gravekamp C, Kasper DL, Michel JL, Kling DE, Carey V, Madoff LC: Immunogenicity and protective efficacy of the alpha C protein of group B streptococci are inversely related to the number of repeats. Infect Immun. 1997, 65: 5216-5221.

  54. 54.

    Madoff LC, Michel JL, Gong EW, Kling DE, Kasper DL: Group B streptococci escape host immunity by deletion of tandem repeat elements of the alpha C protein. Proc Natl Acad Sci USA. 1996, 93: 4131-4136. 10.1073/pnas.93.9.4131.

  55. 55.

    Schubert A, Zakikhany K, Schreiner M, Frank R, Spellerberg B, Eikmanns BJ, Reinscheid DJ: A fibrinogen receptor from group B Streptococcus interacts with fibrinogen by repetitive units with novel ligand binding sites. Mol Microbiol. 2002, 46: 557-569. 10.1046/j.1365-2958.2002.03177.x.

  56. 56.

    Rosenau A, Martins K, Amor S, Gannier F, Lanotte P, van der Mee-Marquet N, Mereghetti L, Quentin R: Evaluation of the ability of Streptococcus agalactiae strains isolated from genital and neonatal specimens to bind to human fibrinogen and correlation with characteristics of the fbsA and fbsB genes. Infect Immun. 2007, 75: 1310-1317. 10.1128/IAI.00996-06.

Download references


This work was presented in part at the 20 European Congress of Clinical Microbiology and Infectious Diseases (ECCMID) in Vienna, April 2010 (poster No P 1698). We thank Nicolas Bery for the initial trials and Mazen Salloum.

Author information



Corresponding author

Correspondence to Philippe Lanotte.

Additional information

Authors' contributions

EH and GB carried out the molecular genetic studies by MLST and MLVA. CP performed BioNumerics analysis of data and helped to draft the manuscript. MFL and ASD contributed to MLST analysis. AR and RQ participated in the design of the study. LM participated in the design of the study and helped to draft the manuscript. EH and PL conceived the study and draft the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Haguenoer, E., Baty, G., Pourcel, C. et al. A multi locus variable number of tandem repeat analysis (MLVA) scheme for Streptococcus agalactiae genotyping. BMC Microbiol 11, 171 (2011).

Download citation


  • Mastitis
  • Genomic Island
  • Clonal Complex
  • Unweighted Pair Group Method With Arithmetic
  • Streptococcus Agalactiae