Skip to main content


Genetic relationships between clinical and non-clinical strains of Yersinia enterocolitica biovar 1A as revealed by multilocus enzyme electrophoresis and multilocus restriction typing

Article metrics



Genetic relationships among 81 strains of Y. enterocolitica biovar 1A isolated from clinical and non-clinical sources were discerned by multilocus enzyme electrophoresis (MLEE) and multilocus restriction typing (MLRT) using six loci each. Such studies may reveal associations between the genotypes of the strains and their sources of isolation.


All loci were polymorphic and generated 62 electrophoretic types (ETs) and 12 restriction types (RTs). The mean genetic diversity (H) of the strains by MLEE and MLRT was 0.566 and 0.441 respectively. MLEE (DI = 0.98) was more discriminatory and clustered Y. enterocolitica biovar 1A strains into four groups, while MLRT (DI = 0.77) identified two distinct groups. BURST (Based Upon Related Sequence Types) analysis of the MLRT data suggested aquatic serotype O:6,30-6,31 isolates to be the ancestral strains from which, clinical O:6,30-6,31 strains might have originated by host adaptation and genetic change.


MLEE revealed greater genetic diversity among strains of Y. enterocolitica biovar 1A and clustered strains in four groups, while MLRT grouped the strains into two groups. BURST analysis of MLRT data nevertheless provided newer insights into the probable evolution of clinical strains from aquatic strains.


Yersinia enterocolitica is an important food- and water-borne gastrointestinal agent. It is known to cause a variety of syndromes ranging from mild gastroenteritis to more invasive diseases like terminal ileitis and mesenteric lymphadenitis mimicking appendicitis [1]. Blood transfusion associated septicaemia due to Y. enterocolitica has been reported to have high mortality [2]. Post infectious sequelae include reactive arthritis and erythema nodosum [1].

Y. enterocolitica is classified into six biovars (1A, 1B, 2, 3, 4 and 5) and more than 50 serotypes [3]. On the basis of pathogenicity, it has been grouped into highly pathogenic (biovar 1B), moderately pathogenic (biovars 2-5) and the so called non-pathogenic (biovar 1A) biovars. Recently, using comparative phylogenomics, Howard et al [4] suggested that these groups might represent three subspecies of Y. enterocolitica. The biovar 1A strains are quite heterogeneous serologically and have been isolated from a variety of sources viz. stools of diarrheic humans, animals, food and aquatic sources [5]. The biovar 1A strains are thought to be non-pathogenic as they lack pYV (plasmid for Yersinia virulence) plasmid and major chromosomal virulence determinants [1]. However, some biovar 1A strains are known to produce symptoms indistinguishable from that produced by the pathogenic biovars [6, 7]. Y. enterocolitica biovar 1A has also been implicated in nosocomial [8] and food-borne [9] outbreaks. A serotype O:6,30 (biovar 1A) strain was reported to cause placentitis and abortion in pregnant ewes [10]. Y. enterocolitica biovar 1A was the most predominant biovar isolated from both livestock and humans during a survey in Great Britain in 1999-2000 and surely needs to be studied further [11]. Several recent studies suggest that these strains might possess novel, as yet unidentified, virulence determinants [1216].

Serological heterogeneity notwithstanding, Y. enterocolitica biovar 1A has only limited genetic heterogeneity as revealed by different genotyping methods such as repetitive elements sequence-based PCR (rep-PCR) fingerprinting [17], 16S-23S intergenic spacer (IGS) region and gyrB restriction fragment length polymorphism [18], and multilocus variable number tandem repeat analysis (MLVA) [19]. Overall, these studies revealed presence of two clonal groups among biovar 1A strains. These studies also showed that clinical and non-clinical serotype O:6,30-6,31 (biovar 1A) strains clustered into two separate groups but failed to reveal any unequivocal associations between genotypes and the source of isolation.

Multilocus enzyme electrophoresis (MLEE) is an important tool used to study genetic relationships where allelic variations in housekeeping genes are indexed using electrophoretic mobilities of corresponding enzymes [20, 21]. The technique has been used to study epidemiology of several pathogenic bacteria [2226]. Multilocus restriction typing (MLRT), a recently developed tool, analyses restriction fragment length polymorphism of several housekeeping genes [2729].

The objective of this study was to use MLEE and MLRT to gain further insight into the genetic heterogeneity and relationships among clinical and non-clinical strains of Y. enterocolitica biovar 1A.


Bacterial strains

Eighty one strains of Y. enterocolitica biovar 1A were examined in this study. Of these, sixty-five were isolated from clinical and non-clinical sources in India viz. diarrheic human patients (35), wastewater (18), swine (7) and pork (5) [3032]. All isolates have been authenticated, and deposited with Yersinia National Reference Laboratory and WHO Collaborating Centre, Institut Pasteur, Paris (France). Of the remaining 16 isolates, ten were obtained from Elisabeth Carniel (Yersinia National Reference Laboratory and WHO Collaborating Centre, Institut Pasteur, Paris, France) and six from Jürgen Heesemann (Max von Pattenkofer Institute, Munich, Germany). Y. enterocolitica 8081 (biovar 1B, serotype O:8), kindly provided by Mikael Skurnik (Haartman Institute, Finland) was used as the reference strain for both MLEE and MLRT.

The serotypes, sources of isolation, country of origin and reference laboratory accession numbers of these strains have been reported previously [17]. All strains were maintained as glycerol stocks at -40°C.

Multilocus enzyme electrophoresis (MLEE)

The enzyme extracts were prepared as per the method described by Selander et al [20]. Briefly, cultures grown overnight in tryptone soy broth (TSB) were harvested by centrifugation at 10,000 g for 10 min at 4°C. The cells were washed twice in potassium phosphate buffer (0.15 M, pH 7.0) and the pellet was resuspended in 2 ml of buffer (10 mM Tris-HCl, 1 mM EDTA and 0.5 mM NADP, pH 6.8). The bacteria were lysed by sonication (Sonics) on ice and centrifuged at 13,000 g for 30 min at 4°C to obtain the supernatant (enzyme extract), which was stored in aliquots of 200 μl each at -40°C until use.

The enzyme extracts were subjected to horizontal gel electrophoresis in 0.9% (w/v) agarose and stained for specific enzyme activities according to the procedures described by Selander et al [20]. The enzymes studied were: malate dehydrogenase (MDH; EC, malic enzyme (ME; EC, glucose-6-phosphate dehydrogenase (G6P; EC, isocitrate dehydrogenase (IDH; EC, alpha esterase (EST-A; EC and glutamate dehydrogenase (GD2; EC The enzymes MDH, ME, G6P and IDH were electrophoresed in Tris citrate buffer (pH 8.0). For EST-A, potassium phosphate buffer (gel buffer, pH 7.0; electrode buffer, pH 6.7) was used and GD2 was electrophoresed in a lithium hydroxide buffer (gel buffer, pH 8.3; electrode buffer, pH 8.1).

Replicate samples from reference strain were run on each gel, which facilitated comparison of the gels. The mobilities of the enzymes from different samples on the same gel were compared. For each enzyme, the distinct mobility variants were designated as electromorphs and numbered in order of decreasing rate of anodal migration. The electromorphs of an enzyme were equated with alleles at the corresponding structural gene locus. Each strain was characterized on the basis of combination of its electromorphs obtained for the six enzymes. Distinct profiles of electromorphs corresponding to multilocus genotypes were designated as electrophoretic types (ETs).

Statistical analyses

Computer programs written by Prof T. S. Whittam were used to analyze the ET data and calculation of genetic diversity [20]. Genetic diversity (h) at an enzyme locus (i.e., the probability that two isolates differ at the j locus) was calculated from the allele frequencies as h j = n (1 - Σx i 2)/n - 1), where x i is the frequency of the i th allele at the j locus and n is the number of isolates [33]. Mean genetic diversity per locus (H) was calculated as the arithmetic average of h values for all loci. The genetic distances between pairs of ETs were calculated as the proportions of loci at which dissimilar electromorphs occurred. Clustering of data was performed from a matrix of pairwise genetic distances by the average-linkage method (unweighted pair group method using arithmetic averages or UPGMA).

Multilocus restriction typing (MLRT)

Genomic DNA was extracted using DNeasy tissue kit (Qiagen) as per the manufacturer's instructions. The six genes encoding housekeeping enzymes: malate dehydrogenase (mdh), adenylate cyclase (cya), glutamine synthetase (glnA), glucose-6-phosphate dehydrogenase (zwf), isocitrate dehydrogenase (icdA) and glutamate dehydrogenase (gdhA) were selected. For amplification of these genes, Yersinia consensus primers were designed using nucleotide sequences from Y. enterocolitica 8081 (biovar 1B, AM286415), Y. pestis (AE009952) and Y. pseudotuberculosis (BX936398) available at EMBL and GenBank databases, after pairwise alignment of the sequences using ClustalW Primers were designed with PRIMER SELECT software (DNAStar), and synthesized from Microsynth. The details of the primers are given in Table 1.

Table 1 Details of primers and restriction enzymes used for multilocus restriction typing (MLRT) of Y. enterocolitica biovar 1A

Polymerase chain reactions were performed in 25 μl of reaction mixture containing 1 × PCR buffer (10 mM Tris-HCl pH 8.8, 50 mM KCl, 0.1% Triton X-100, 1.5 mM MgCl2), 200 μM of each dNTP (MBI Fermentas), 20 pmoles each of forward and reverse primers, 2 U DyNAzyme™ II DNA polymerase (Finnzymes) and 100 ng of template DNA. All amplifications were performed in a PTC-100™ thermal cycler (MJ Research) according to the following cycling conditions: initial denaturation for 5 min at 94°C, 30 amplification cycles each consisting of 1 min denaturation at 94°C, annealing for 45 s at the temperatures as given in Table 1, and 1 min elongation at 72°C. The final extension was carried out at 72°C for 10 min. 5 μl of the PCR product was electrophoresed in 1% (w/v) agarose gel containing 0.5 μg ml-1 ethidium bromide (EtBr) at 80 V for 1 h in 1 × Tris-acetate EDTA buffer (1 × TAE: 40 mM Tris acetate, 1 mM EDTA, pH 8.0). The 100 bp DNA ladder (New England Biolabs) served as the molecular size marker.

The restriction enzymes for MLRT were selected by an in silico restriction analysis of respective gene sequences of Y. enterocolitica 8081 (biovar 1B) available in GenBank using MapDraw (DNAStar) such that polymorphism in the restriction sites was revealed. The PCR amplicons of six genes for all the 81 strains were digested with enzymes as shown in Table 1. Restriction digestion was carried out overnight at 37°C in 25 μl reaction mixture containing 8 μl of the PCR amplicon, 2.5 μl of 10 × buffer and 2 U of restriction enzyme (New England Biolabs). Restriction digests were analyzed by agarose gel electrophoresis (2.5% gel containing 0.5 μg ml-1 EtBr in 1 × TAE buffer). Gels were run at 60 V and photographed under UV transillumination. The 50 bp and 100 bp DNA ladders (New England Biolabs or MBI Fermentas) served as the molecular weight standards.

The restriction patterns for all the isolates were analyzed using Diversity Database Software (version 2, Bio-Rad). Distinct restriction patterns for each locus were considered to represent separate alleles, and each allele was assigned a numeral. As with MLEE, the combination of alleles at each of the six loci gave a restriction type (RT). Strains were considered different if the allele of any of the six loci differed. The genetic diversity h was calculated as described for MLEE. The restriction profile for each isolate was entered into a database and used to construct a phylogenetic tree based on unweighted-pair group method with average (UPGMA) linkage of distance, using the START (Sequence Type Analysis and Recombination Tests) software package In addition, clonal complexes within 81 biovar 1A strains were investigated using the BURST (Based Upon Related Sequence Types) algorithm of START software package.

DNA sequencing and analysis

For each allele identified for the six genes used in MLRT, one amplicon was sequenced to confirm its identity. PCR products were purified with the QIAquick gel extraction kit (Qiagen) and DNA sequencing was performed by the Big-Dye terminator kit using an automated DNA sequencer (ABI PRISM 3730 genetic analyzer).

Linkage disequilibrium analysis

Linkage disequilibrium for MLEE and MLRT data was calculated on the basis of the distribution of allelic mismatches between pairs of bacterial isolates among all the loci examined. The ratio of the variance observed (VO) in mismatches to the variance expected (VE) at linkage equilibrium provides a measure of multilocus linkage disequilibrium and can be expressed as the index of association (IA) as: IA = VO/VE - 1 [34, 35]. For populations in linkage equilibrium, VO = VE and IA is not significantly different from zero, whereas values of IA significantly greater than zero indicate that recombination has been rare or absent. To determine whether VO was significantly different from VE in any sample, a Monte Carlo procedure was iterated, wherein alleles are repeatedly scrambled to eliminate any effect of linkage disequilibrium [36]. The LIAN version 3.5 software program [37] was used to calculate IA and standardized IA (ISA) values and perform Monte Carlo procedure.

Calculation of diversity index

Simpson's diversity index (DI), a measure of the discriminatory ability of a given typing method, was calculated for MLEE and MLRT as described by Hunter and Gaston [38].

Nucleotide sequence accession number

The nucleotide sequence data of six genes used in MLRT study reported in this paper have been deposited in GenBank database under the accession numbers FJ899547-FJ899554 and GQ229153-GQ229162.


Electrophoretic types (ETs) and genetic diversity

Activities of six enzymes were detected in all 81 strains of Y. enterocolitica biovar 1A and the reference strain Y. enterocolitica 8081. All enzyme loci studied were polymorphic and the number of alleles ranged from three (isocitrate dehydrogenase) to fifteen (glucose-6-phosphate dehydrogenase) (Table 2). The mean number of alleles per locus was 7.5. Esterase was the most polymorphic (h = 0.827), while glutamate dehydrogenase was the least polymorphic locus (h = 0.250). The mean genetic diversity (H) of all strains was 0.566 ± 0.088. Among the 81 Y. enterocolitica biovar 1A strains, 62 ETs (electrophoretic types) were identified. The reference strain Y enterocolitica 8081 formed a distinct ET, ET63 (Table 3). Fifty seven ETs were represented only once in the data set. The ETs which were represented more than once were ET1 (ten isolates), ET8 (six isolates), ET6 (three isolates), and ETs 20, 36 and 42 (two isolates each).

Table 2 Genetic diversity at six enzyme loci in Y. enterocolitica biovar 1A
Table 3 Details of electrophoretic types (ETs) and restriction types (RTs) of Y. enterocolitica biovar 1A

The genetic relationships among strains of Y. enterocolitica biovar 1A as revealed by cluster analysis using UPGMA are shown in Figure 1. The shortest genetic distance (0.167) between the ETs corresponded to a single locus difference. The strains were grouped into 4 groups (I to IV) diverging at genetic distance of 0.76. The group I comprising 38 ETs (ET1-20, 22-24, 26, 29-35, 42-44, 49, 60-62) was the largest with 56 isolates belonging to different serotypes and sources. This group was highly diverse with several subclusters. This group also contained the most common ET, ET1 which was represented by 9 clinical isolates belonging to serotypes O:6,30-6,31 (3 isolates) and O:6,30 (6 isolates), and one pork isolate of serotype O:7,8-8-8,19. Another ET, ET18 was also predominant and contained 6 Indian strains which included three wastewater serotype O:6,30-6,31 isolates, one wastewater serotype O:10-34 isolate and two NAG isolates one each of aquatic and clinical source. Group II included 4 ETs (ET56-59) containing one pig throat isolate and 3 clinical isolates. Group III was formed by 18 isolates representing 17 ETs (ET 21, 25, 27, 28, 36-41, 48, 50-55). These strains belonged to diverse serotypes and sources from India (15 isolates) and France (3 isolates). The three French isolates formed a separate subgroup at a genetic distance of 0.64. Group IV included three European clinical serotype O:6,30 isolates representing ETs 45-47. MLEE dendrogram revealed that ET1 and ET36 represented by multiple isolates showed close association (linkage distance = 0.0) between isolates from pork/pig throat and human.

Figure 1

UPGMA dendrogram showing genetic relationships among 62 electrophoretic types (ETs) of Y. enterocolitica biovar 1A. NAG: non-agglutinable, ND: not determined, NK: not known.

Multilocus restriction typing

PCR amplicons were obtained for all six loci using primers and PCR conditions given in Table 1. For each of the six loci, PCR amplicons of respective sizes were obtained for all the 81 strains of Y. enterocolitica biovar 1A. The amplicons were digested with restriction enzymes as shown in Table 1. The RFLP profiles for each of the six loci are given in Additional file 1. Collating the PCR-RFLP data for six loci in 81 strains, 12 restriction types (RTs) were identified (Table 3). Reference strain Y. enterocolitica 8081 (biovar 1B, serotype O:8) was represented by a distinct RT, RT13. RT1 was the most common restriction type and was present among 31 (37%) isolates. The second commonest type was RT2, represented by 20 (25%) isolates while RT3 was the third commonest (15 isolates, 19%) restriction type. Reproducibility of MLRT was checked by repeating RFLP using selected isolates. Same allelic profiles were obtained indicating reproducibility of MLRT.

The number of alleles present per locus and genetic diversity among 81 strains of Y. enterocolitica biovar 1A as determined by MLRT are given in Table 2. Glucose-6-phosphate dehydrogenase (zwf) locus was the most diverse (h = 0.644) while isocitrate dehydrogenase (icdA) was least diverse (h = 0.336). The mean genetic diversity (H) of all isolates was 0.441 ± 0.048.

The genetic relationships among strains analyzed by cluster analysis using UPGMA are shown in Figure 2. MLRT clustered biovar 1A strains into two clonal groups (A and B) while the reference strain (Y. enterocolitica 8081, biovar 1B) formed a separate group, at the linkage distance of 0.78. The group A comprising most (64 of 81) of the isolates was represented by 9 different RTs. Within the group A, two subgroups were identified namely, A-I and A-II. In subgroup A-I, all wastewater serotype O:6,30-6,31 isolates, human NAG and European O:6,30 isolates were present. Subgroup A-II comprised of all clinical O:6,30-6,31 isolates, most clinical O:6,30 isolates, three pork and pig throat isolates each, and five wastewater isolates belonging to different serotypes. The most common RT, RT1 representing 31 isolates was present in this subgroup. The group B comprised of 15 isolates belonging to RT3 and a single isolate each of RT8 and RT11. Genotypically, this group was quite homogeneous despite belonging to different serotypes, sources and geographic origin.

Figure 2

Dendrogram showing relationships of Y. enterocolitica biovar 1A strains based on analysis of restriction types (RTs) generated by MLRT. The dendrogram was constructed using UPGMA algorithm available in the START software package. NAG: non-agglutinable, ND: not determined, NK: not known.

The analysis of MLRT data by BURST program identified two clonal complexes (Figure 3) corresponding to the clonal groups identified above. The clonal complex A comprising 9 RTs (64 strains) revealed that wastewater serotype O:6,30-6,31 isolates represented by RT2 were present in the innermost circle as ancestral strains. The clinical serotype O:6,30-6,31 strains represented by RT1 and RT12 were present in the outer circle as single locus variants (Figure 3a) The double locus variants (RT5 and RT9) and the satellite RTs (RT6 and RT10) were represented by serotypes which are relatively not common. However, not much information could be inferred from clonal complex B (Figure 3b).

Figure 3

Clonal complexes identified among 81 strains of Y. enterocolitica biovar 1A by BURST analysis of MLRT data. a) Clonal complex A, b) Clonal complex B. Each number denotes a restriction type (RT; refer to Figure 2). Radial distribution shows divergent RTs. Ancestral RT is shown in the innermost circle. Single locus variants (SLV) are shown in the second circle and double locus variants (DLV) are represented in the outermost circle. Satellite RTs (RTs present outside the outermost circle) vary by more than two loci from the ancestral type. Lines indicate whether the RT is SLV (solid line) or DLV (dashed line).

Sequencing of amplicons from representative strains confirmed the identity of the genes. Analysis of the sequences also confirmed the restriction patterns observed for each of the six genes. This is the first report on MLRT of Y. enterocolitica.

Analysis of linkage disequilibrium and discriminatory indices

The frequency of recombination in natural populations can be estimated by calculating index of association (IA) between loci [35]. The results of the analysis of multilocus linkage disequilibrium in Y. enterocolitica are summarized in Table 4. The IA and ISA values for the 81 strains studied by MLEE were 0.613 and 0.128 respectively, which differed significantly (p < 0.001) from zero indicating that the strains were in linkage disequilibrium. Similarly, significant level of linkage disequilibrium was observed on analysis of MLRT data. The IA and ISA values were 3.357 and 0.672 respectively, and differed significantly (p < 0.001) from zero. Simpson's diversity index (DI) for MLEE and MLRT was 0.98 and 0.77 respectively.

Table 4 Multilocus linkage disequilibrium analysis of Y. enterocolitica biovar 1A strains


Indexing allelic variations in sets of housekeeping genes provides a good measure of overall genetic heterogeneity in populations of microorganisms [21]. Methods based on this principle such as MLEE, MLRT and MLST (multilocus sequence typing) provide good insight into the genetic relationships among strains. In the present study, we used MLEE and MLRT to assess the genetic relationships among 81 strains of Y. enterocolitica biovar 1A isolated from clinical and non-clinical sources.

MLEE clustered Y. enterocolitica biovar 1A into four groups. A close analysis of data presented by Dolina and Peduzzi [23] who studied human, animal and aquatic strains of Y. enterocolitica isolated from Switzerland by MLEE, revealed that 51 biovar 1A strains clustered into two major groups, although minor clusters having one and six isolates each were also observed. Another study that used fluorescent amplified fragment length polymorphism (FAFLP) also clustered biovar 1A strains into two groups: one group comprised of biovar 1A strains; while a few biovar 1A strains clustered with atypical pathogenic biovars constituting the second group [39]. Further study by comparative genomic DNA microarray however showed that these biovar 1A strains constituted a single group [4]. Other studies using rep-PCR genotyping [17], 16S-23S IGS and gyr B RFLP [18], and MLVA [19] have also clustered biovar 1A strains into two clonal groups. MLEE revealed a total of 62 electrophoretic types (ETs) among 81 biovar 1A strains and showed high degree of discrimination (DI = 0.98). Studies of allelic variation by MLEE also revealed sufficient genetic diversity (H = 0.566) among strains of Y. enterocolitica biovar 1A. Similar genetic diversity was also reported in previous MLEE studies on Y. enterocolitica[22, 23].

In the present study however, based on the number of distinct ETs generated, the clinical serotype O:6,30 and O:6,30-6,31 isolates were shown to be heterogeneous with mean genetic diversities (H) of 0.514 ± 0.112 and 0.442 ± 0.078 respectively. Previous studies in which other techniques namely rep-PCR [17], 16S-23S IGS and gyr B RFLP [18], and MLVA [19] were used to type these strains did not reveal this heterogeneity. Fearnley et al [39] also reported heterogeneity among serotype O:6,30 strains wherein seven AFLP types were identified among eight strains.

In the MLEE dendrogram, two ETs showed some pork and pig strains to be identical to the strains isolated from diarrheic human subjects suggesting that like pathogenic biovars [11, 22, 40], pigs may be the source of biovar 1A strains isolated from human patients. No such grouping of human and pork/pig isolates was evident from earlier studies [17, 18]. However, this observation needs to be explored further by making use of a larger number of pig/pork isolates belonging to biovar 1A.

Multilocus restriction typing (MLRT) has recently been used to discern phylogenetic relationships among strains of Streptococcus pneumoniae[41], Neisseria meningitidis[28, 42], Burkholderia cepacia[27, 43], Staphylococcus aureus[44] and Escherichia coli[29]. MLRT has been reported to show good correlation with PFGE [27, 29] and has been advocated as a cost effective alternative to MLST, which is relatively an expensive technique [28, 42]. In the present study, MLRT divided 81 strains of Y. enterocolitica biovar 1A into 12 RTs based on a combination criteria of number of alleles and restriction patterns observed at each of the six loci examined. Cluster analysis of MLRT data revealed two clonal groups - A and B. The reference strain Y. enterocolitica 8081 (biovar 1B) formed a distinct RT. Although MLRT profiles showed good reproducibility, the method failed to rival the discriminatory ability of MLEE. In the context of Y. enterocolitica biovar 1A, the discriminatory ability of MLRT (DI = 0.77) was lower than even rep-PCR (DI = 0.84) [17] and MLVA (DI = 0.87) [19].

Two clonal complexes were identified following BURST analysis of MLRT data. The primary clonal complex contained all but 3 RTs, representing 78% of the isolates. The other complex contained the remaining strains. The approach used in the BURST analysis specifically examines the relationships between closely related genotypes in the clonal complexes [45]. This analysis revealed that in the primary clonal complex, wastewater serotype O:6,30-6,31 isolates represented the ancestral strains while, clinical serotype O:6,30-6,31 strains occupied radial position as single locus variants. This observation corroborates the recent findings obtained from the study of VNTR loci which also suggested that the clinical serotype O:6,30-6,31 strains probably originated from the wastewater strains, by host adaptation and genetic change [19].

The analysis of linkage disequilibrium indicated clonal structure for Y. enterocolitica biovar 1A as values of IA and ISA were found to be significantly different from zero for both MLEE and MLRT data. Other genera, which have been reported to have clonal structure, include Salmonella enterica serovar Paratyphi B [46], Mycobacterium spp. [47], Vibrio cholerae[24] and Pseudomonas stutzeri[25].

Both MLEE and MLRT showed European strains to be more heterogeneous than the Indian strains. MLEE revealed that each of the 15 strains from France and Germany had distinct electrophoretic profiles indicating their heterogeneity. MLRT also revealed that the European strains, which displayed 5 RTs were more heterogeneous compared to Indian isolates. Genetic heterogeneity of European biovar 1A strains has been reported earlier using PFGE [48] and FAFLP [39]. A previous study using multilocus variable number tandem repeat analysis also identified 13 MLVA types among 15 European biovar 1A strains [19]. This suggests that European and Indian strains may constitute separate groups and might be evolving independently in two different settings. It would be interesting to explore these evolutionary aspects by comparative whole genome sequencing or multilocus sequence typing of Indian and European strains. It was also observed that strains with different serotypes (O antigen) types produced identical ETs or RTs and were closely related genetically. Also, in some cases, same O antigen was shared by strains that were different genotypically. These observations indicate O antigen switching in strains of Y. enterocolitica as suggested recently by MLST [49]. Such observations have however been reported in other bacteria also [24, 41, 50]. Thus, given the enormous discriminatory power of genotyping techniques such observations also emphasize the need to discuss threadbare, the question of suitability of widely used typing techniques like serotyping.


More diversity was observed among clinical and non-clinical strains of Y. enterocolitica biovar 1A when MLEE was used. Sixty-two electrophoretic types were identified among 81 strains, which clustered into four distinct groups. MLRT identified 12 restriction types and was distinctly less discriminatory, clustering the strains into two groups. The BURST analysis of the MLRT data nevertheless provided newer insights into the probable evolution of clinical strains from those present in the aquatic environments.


  1. 1.

    Bottone EJ: Yersinia: enterocolitica overview and epidemiologic correlates. Microbes Infect. 1999, 1: 323-333. 10.1016/S1286-4579(99)80028-8.

  2. 2.

    Leclercq A, Martin L, Vergnes ML, Ounnoughene N, Laran JF, Giraud P, Carniel E: Fatal Yersinia enterocolitica biotype 4 serovar O:3 sepsis after red blood cell transfusion. Transfusion. 2005, 45: 814-818. 10.1111/j.1537-2995.2005.04363.x.

  3. 3.

    Cornelis G, Laroche Y, Balligand G, Sory MP, Wauters G: Yersinia enterocolitica a primary model for bacterial invasiveness. Rev Infect Dis. 1987, 9: 64-87.

  4. 4.

    Howard SL, Gaunt MW, Hinds J, Witney AA, Stabler R, Wren BW: Application of comparative phylogenomics to study the evolution of Yersinia enterocolitica and to identify genetic differences relating to pathogenicity. J Bacteriol. 2006, 188: 3645-3653. 10.1128/JB.188.10.3645-3653.2006.

  5. 5.

    Tennant SM, Grant TH, Robins-Browne RM: Pathogenicity of Yersinia enterocolitica biotype 1A. FEMS Immunol Med Microbiol. 2003, 38: 127-137. 10.1016/S0928-8244(03)00180-9.

  6. 6.

    Morris JG, Prado V, Ferreccio C, Robins-Browne RM, Bordun AM, Cayazzo M, Kay BA, Levine MM: Yersinia enterocolitica isolated from two cohorts of young children in Santiago Chile: incidence of and lack of correlation between illness and proposed virulence factors. J Clin Microbiol. 1991, 29: 2784-2788.

  7. 7.

    Burnens AP, Frey A, Nicolet J: Association between clinical presentation biogroups and virulence attributes of Yersinia enterocolitica strains in human diarrhoeal disease. Epidemiol Infect. 1996, 116: 27-34. 10.1017/S0950268800058921.

  8. 8.

    Ratnam S, Mercer E, Picco B, Parsons S, Butler R: A nosocomial outbreak of diarrheal disease due to Yersinia enterocolitica serotype O:5, biotype 1. J Infect Dis. 1982, 145: 242-247.

  9. 9.

    Greenwood MH, Hooper WL: Excretion of Yersinia spp associated with consumption of pasteurized milk. Epidemiol Infect. 1990, 104: 345-350. 10.1017/S0950268800047361.

  10. 10.

    Corbel MJ, Ellis B, Richardson C, Bradley R: Experimental Yersinia enterocolitica placentitis in sheep. Br Vet J. 1992, 148: 339-349.

  11. 11.

    McNally A, Cheasty T, Fearnley C, Dalziel RW, Paiba GA, Manning G, Newell DG: Comparison of the biotypes of Yersinia enterocolitica isolated from pigs, cattle and sheep at slaughter and from humans with yersiniosis in Great Britain during 1999-2000. Lett Appl Microbiol. 2004, 39: 103-108. 10.1111/j.1472-765X.2004.01548.x.

  12. 12.

    Grant T, Bennett-Wood V, Robins-Browne RM: Characterization of the interaction between Yersinia enterocolitica biotype 1A and phagocytes and epithelial cells in vitro. Infect Immun. 1999, 67: 4367-4375.

  13. 13.

    Singh I, Virdi JS: Production of Yersinia stable toxin (YST) and distribution of yst genes in biotype 1A strains of Yersinia enterocolitica. J Med Microbiol. 2004, 53: 1065-1068. 10.1099/jmm.0.45527-0.

  14. 14.

    McNally A, Dalton T, La Ragione RM, Stapleton K, Manning G, Newell DG: Yersinia enterocolitica isolates of differing biotypes from humans and animals are adherent, invasive and persist in macrophages, but differ in cytokine secretion profiles in vitro. J Med Microbiol. 2006, 55: 1725-1734. 10.1099/jmm.0.46726-0.

  15. 15.

    McNally A, La Ragione RM, Best A, Manning G, Newell DG: An aflagellate mutant Yersinia enterocolitica biotype 1A strain displays altered invasion of epithelial cells, persistence in macrophages, and cytokine secretion profiles in vitro. Microbiology. 2007, 153: 1339-1349. 10.1099/mic.0.2006/000919-0.

  16. 16.

    Bhagat N, Virdi JS: Distribution of virulence-associated genes in Yersinia enterocolitica biovar 1A correlates with clonal groups and not the source of isolation. FEMS Microbiol Lett. 2007, 266: 177-183. 10.1111/j.1574-6968.2006.00524.x.

  17. 17.

    Sachdeva P, Virdi JS: Repetitive elements sequence (REP/ERIC)-PCR based genotyping of clinical and environmental strains of Yersinia enterocolitica biotype 1A reveal existence of limited number of clonal groups. FEMS Microbiol Lett. 2004, 240: 193-201. 10.1016/j.femsle.2004.09.029.

  18. 18.

    Gulati PS, Virdi JS: The rrn locus and gyr B genotyping confirm the existence of two clonal groups in strains of Yersinia enterocolitica subspecies palearctica biovar 1A. Res Microbiol. 2007, 158: 236-243. 10.1016/j.resmic.2006.11.011.

  19. 19.

    Gulati P, Varshney RK, Virdi JS: Multilocus variable number tandem repeat analysis as a tool to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A. J Appl Microbiol. 2009, 107: 875-884. 10.1111/j.1365-2672.2009.04267.x.

  20. 20.

    Selander RK, Caugant DA, Gilmour MN, Whittam TS: Methods of multilocus enzyme electrophoresis for bacterial population genetics and systematic. Appl Environ Microbiol. 1986, 51: 873-884.

  21. 21.

    Musser JM: Molecular population genetic analysis of emerged bacterial pathogens: selected insights. Emerg Infect Dis. 1996, 2: 1-17. 10.3201/eid0201.960101.

  22. 22.

    Caugant DA, Aleksic S, Mollaret HH, Selander RK, Kapperud G: Clonal diversity and relationship among strains of Yersinia enterocolitica. J Clin Microbiol. 1989, 27: 2678-2683.

  23. 23.

    Dolina M, Peduzzi R: Population genetics of human, animal, and environmental Yersinia strains. Appl Environ Microbiol. 1993, 59: 442-450.

  24. 24.

    Farfán M, Miñana D, Fusté MC, Lorén JG: Genetic relationships between clinical and environmental Vibrio cholerae isolates based on multilocus enzyme electrophoresis. Microbiology. 2000, 146: 2613-2626.

  25. 25.

    Rius N, Fuste MC, Guasp C, Lalucat J, Loren JG: Clonal population structure of Pseudomonas stutzeri a species with exceptional genetic diversity. J Bacteriol. 2001, 183: 736-744. 10.1128/JB.183.2.736-744.2001.

  26. 26.

    Scortichini M, Natalini E, Angelucci L: Clonal population structure of Pseudomonas avellanae strains of different origin based on multilocus enzyme electrophoresis. Microbiology. 2003, 149: 2891-2900. 10.1099/mic.0.26380-0.

  27. 27.

    Coenye T, LiPuma JJ: Multilocus restriction typing: A novel tool for studying global epidemiology of Burkholderia cepacia complex infection in cystic fibrosis. J Infect Dis. 2002, 185: 1454-1462. 10.1086/340279.

  28. 28.

    Bennett DE, Cafferkey MT: Multilocus restriction typing: A tool for Neisseria meningitidis strain discrimination. J Med Microbiol. 2003, 52: 781-787. 10.1099/jmm.0.05225-0.

  29. 29.

    Helgerson AF, Sharma V, Dow AM, Schroeder R, Post K, Cornick NA: Edema disease caused by a clone of Escherichia coli O147. J Clin Microbiol. 2006, 44: 3074-3077. 10.1128/JCM.00617-06.

  30. 30.

    Singh I, Virdi JS: Isolation biochemical characterization and in vitro tests of pathogenicity of Yersinia enterocolitica isolated from pork. Curr Sci. 1999, 77: 1019-1021.

  31. 31.

    Sinha I, Choudhary I, Virdi JS: Isolation of Yersinia enterocolitica and Yersinia intermedia from wastewaters and their biochemical and serological characteristics. Curr Sci. 2000, 79: 510-513.

  32. 32.

    Singh I, Bhatnagar S, Virdi JS: Isolation and characterization of Yersinia enterocolitica from diarrheic human subjects and other sources. Curr Sci. 2003, 84: 1353-1355.

  33. 33.

    Nei M: Estimation of average heterozygosity and genetic distance from a small sample of individuals. Genetics. 1978, 89: 583-590.

  34. 34.

    Brown AH, Feldman MW, Nevo E: Multilocus structure of natural populations of Hordeum spontaneum. Genetics. 1980, 96: 523-536.

  35. 35.

    Maynard Smith J, Smith NH, O'Rourke M, Spratt BG: How clonal are bacteria?. Proc Nat Acad Sci USA. 1993, 90: 4384-4388. 10.1073/pnas.90.10.4384.

  36. 36.

    Souza V, Nguyen TT, Hudson RR, Piñero D, Lenski RE: Hierarchical analysis of linkage disequilibrium in Rhizobium populations: Evidence for sex?. Proc Natl Acad Sci USA. 1992, 89: 8389-8393. 10.1073/pnas.89.17.8389.

  37. 37.

    Haubold H, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics. 2000, 16: 847-848. 10.1093/bioinformatics/16.9.847.

  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.

    Fearnley C, On SLW, Kokotovic B, Manning G, Cheasty T, Newell DG: Application of fluorescent amplified fragment length polymorphism for comparison of human and animal isolates of Yersinia enterocolitica. Appl Environ Microbiol. 2005, 71: 4960-4965. 10.1128/AEM.71.9.4960-4965.2005.

  40. 40.

    Tauxe RV, Vandepitte J, Wauters G, Martin SM, Goossens V, DeMol P, Van Noyen R, Thiers G: Yersinia enterocolitica infections and pork: the missing link. Lancet. 1987, 1: 1129-1132. 10.1016/S0140-6736(87)91683-7.

  41. 41.

    Muller-Graf CDM, Whatmore AM, King SJ, Trzcinski K, Pickerill AP, Doherty N, Paul J, Griffiths D, Crook D, Dowson CG: Population biology of Streptococcus pneumoniae isolated from oropharyngeal carriage and invasive disease. Microbiology. 1999, 145: 3283-3293.

  42. 42.

    Dyet KH, Simmonds RS, Martin DR: Multilocus restriction typing method to predict the sequence type of meningococci. J Clin Microbiol. 2004, 42: 1742-1745. 10.1128/JCM.42.4.1742-1745.2004.

  43. 43.

    Coenye T, Spilker T, Martin A, LiPuma JJ: Comparative assessment of genotyping methods for epidemiologic study of Burkholderia cepacia genomovar III. J Clin Microbiol. 2002, 40: 3300-3307. 10.1128/JCM.40.9.3300-3307.2002.

  44. 44.

    Diep BA, Perdreau-Remington F, Sensabaugh GF: Clonal characterization of Staphylococcus aureus by multilocus restriction fragment typing, a rapid screening approach for molecular epidemiology. J Clin Microbiol. 2003, 41: 4559-4564. 10.1128/JCM.41.10.4559-4564.2003.

  45. 45.

    Jolley KA, Feil EJ, Chan MS, Maiden MC: Sequence type analysis and recombinational tests (START). Bioinformatics. 2001, 17: 1230-1231. 10.1093/bioinformatics/17.12.1230.

  46. 46.

    Selander RK, Beltran P, Smith NH, Barker RM, Crichton PB, Old DC, Musser JM, Whittam TS: Genetic population structure, clonal phylogeny, and pathogenicity of Salmonella paratyphi B. Infect Immun. 1990, 58: 1891-1901.

  47. 47.

    Feizabadi MM, Robertson ID, Cousins DV, Dawson DJ, Hampson DJ: Use of multilocus enzyme electrophoresis to examine genetic relationships amongst isolates of Mycobacterium intracellulare and related species. Microbiology. 1997, 143: 1461-1469. 10.1099/00221287-143-4-1461.

  48. 48.

    Najdenski H, Iteman I, Carniel E: Efficient subtyping of pathogenic Yersinia enterocolitica strains by pulsed-field gel electrophoresis. J Clin Microbiol. 1994, 32: 2913-2920.

  49. 49.

    Kotetishvili M, Kreger A, Wauters G, Morris JG, Sulakvelidze A, Stine OC: Multilocus sequence typing for studying genetic relationships among Yersinia species. J Clin Microbiol. 2005, 43: 2674-2684. 10.1128/JCM.43.6.2674-2684.2005.

  50. 50.

    Beltrán P, Delgado G, Navarro A, Trujillo F, Selander RK, Cravioto A: Genetic diversity and population structure of Vibrio cholerae. J Clin Microbiol. 1999, 37: 581-590.

Download references


SM acknowledges Senior Research Fellowship from Council for Scientific and Industrial Research, New Delhi, India. The research grants to JSV from Department of Biotechnology, Indian Council of Medical Research and University of Delhi to strengthen R & D doctoral research programme are acknowledged gratefully.

Author information

Correspondence to Jugsharan S Virdi.

Additional information

Authors' contributions

SM carried out the experimental part of the study. JSV conceived and supervised the work. Both authors participated in interpretation of data and preparation of the final manuscript.

Electronic supplementary material

Additional file 1:Representative restriction profiles of six genes of Y. enterocolitica biovar 1A. (a) Malate dehydrogenase (mdh) digested with Hae III and Rsa I; (b) adenylate cyclase (cya) digested with Hae III and Sau96 I; (c) gluamine synthetase (gln A) digested with Sau96 I and Bgl I; (d) glucose-6-phosphate dehydrogenase (zwf) digested with Rsa I and BstN I; (e) isocitrate dehydrogenase (icdA) digested with Rsa I and Hae III; (f) glutamate dehydrogenase (gdhA) digested with Hae III and Mse I. Numbers above lanes represent the name of strain used to obtain the restriction pattern. Digestion products were compared to 100 bp (M) or 50 bp (M') DNA ladder. (JPEG 535 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article


  • Clonal Complex
  • Yersinia Enterocolitica
  • Clonal Group
  • Variable Number Tandem Repeat Analysis
  • Multilocus Enzyme Electrophoresis