Clinical isolates of Yersinia enterocolitica Biotype 1A represent two phylogenetic lineages with differing pathogenicity-related properties
BMC Microbiology volume 12, Article number: 208 (2012)
Y. enterocolitica biotype (BT) 1A strains are often isolated from human clinical samples but their contribution to disease has remained a controversial topic. Variation and the population structure among the clinical Y. enterocolitica BT 1A isolates have been poorly characterized. We used multi-locus sequence typing (MLST), 16S rRNA gene sequencing, PCR for ystA and ystB, lipopolysaccharide analysis, phage typing, human serum complement killing assay and analysis of the symptoms of the patients to characterize 298 clinical Y. enterocolitica BT 1A isolates in order to evaluate their relatedness and pathogenic potential.
A subset of 71 BT 1A strains, selected based on their varying LPS patterns, were subjected to detailed genetic analyses. The MLST on seven house-keeping genes (adk, argA, aroA, glnA, gyrB, thrA, trpE) conducted on 43 of the strains discriminated them into 39 MLST-types. By Bayesian analysis of the population structure (BAPS) the strains clustered conclusively into two distinct lineages, i.e. Genetic groups 1 and 2. The strains of Genetic group 1 were more closely related (97% similarity) to the pathogenic bio/serotype 4/O:3 strains than Genetic group 2 strains (95% similarity). Further comparison of the 16S rRNA genes of the BT 1A strains indicated that altogether 17 of the 71 strains belong to Genetic group 2. On the 16S rRNA analysis, these 17 strains were only 98% similar to the previously identified subspecies of Y. enterocolitica. The strains of Genetic group 2 were uniform in their pathogenecity-related properties: they lacked the ystB gene, belonged to the same LPS subtype or were of rough type, were all resistant to the five tested yersiniophages, were largely resistant to serum complement and did not ferment fucose. The 54 strains in Genetic group 1 showed much more variation in these properties. The most commonly detected LPS types were similar to the LPS types of reference strains with serotypes O:6,30 and O:6,31 (37%), O:7,8 (19%) and O:5 (15%).
The results of the present study strengthen the assertion that strains classified as Y. enterocolitica BT 1A represent more than one subspecies. Especially the BT 1A strains in our Genetic group 2 commonly showed resistance to human serum complement killing, which may indicate pathogenic potential for these strains. However, their virulence mechanisms remain unknown.
Yersinia enterocolitica species has six biotypes (BTs) of which five (1B, 2, 3, 4, 5) contain pathogenic strains. Y. enterocolitica ssp. enterocolitica consists mainly of the strains of BT 1B, which are considered highly virulent. Low-virulent ssp. palearctica encompasses BTs 2–5 and 1A. Since BT 1A strains lack most of the classical virulence markers, this biotype is often considered non-pathogenic. Nevertheless, BT 1A strains are commonly isolated from patients with diarrhoea. Reports supporting the pathogenicity of some BT 1A strains comprise clinical data[1–7] and cell experiments[8–10]. A virulence marker commonly found in BT 1A strains is the gene ystB encoding heat-stable Yersinia enterotoxin B whereas they usually lack the ystA gene found from Y. enterocolitica 4/O:3 strains. Yersinia enterotoxins A and B are homologues to enterotoxins found in enterotoxigenic E. coli (ETEC) and Vibrio cholerae non-O1 strains. Higher rates of diarrhoea, weight loss, and death have been detected when young rabbits were infected with a Y. enterocolitica strain that produces heat-stable enterotoxin compared to the infection with a knock-out mutant. A majority of the Y. enterocolitica BT 1A strains possess the ystB gene and some excrete heat-stable YstB enterotoxin at 37°C in experimental conditions corresponding to those found in ileum[14, 15].
The BT 1A strains are genetically the most heterogeneous of all the Y. enterocolitica biotypes[16–19]. They belong to numerous serotypes, with at least 17 having been identified. It has been suggested that BT 1A should be separated into its own subspecies based on genetic differences on a DNA microarray against Y. enterocolitica ssp. enterocolitica BT 1B strain 8081. Likewise, a number of other studies utilizing different methods have suggested that Y. enterocolitica BT 1A strains could be divided into two main clusters[16, 21–25]. However, since the studies have been conducted on different sets of strains, it is impossible to know whether all the methods would divide the strains into two clusters similarly. Recently, two genome sequences of BT 1A strains with no evident structural differences were published. Notable differences between an environmental serotype O:36 and a clinical BT 1A/O:5 strains were the presence of a Rtx toxin-like gene cluster and remnants of a P2-like prophage in the clinical BT 1A/O:5 isolate.
BT 1A was the predominant biotype of Y. enterocolitica detected among Yersinia isolates from human clinical stool samples in Finland in 2006, as also in other European countries. Of the Finnish patients with a BT 1A strain, 90% suffered from diarrhoea and abdominal pain, but only 35% had fever. Furthermore, 3% of the patients had reactive arthritis compared to 0.3% of the controls. We hypothesized that certain BT 1A strains might have a higher pathogenic potential than others. In order to study this, the clinical BT 1A isolates were investigated using multilocus sequence typing (MLST), 16S rRNA sequencing, yst-PCR, lipopolysaccharide (LPS) analysis, sensitivity to five yersiniophages and serum killing assay. MLST results were analysed with BAPS (Bayesian Analysis of Population Structure) program, genetic and phenotypic characteristics of the BT 1A strains were compared and statistical analysis was applied to assess their correlation with the symptoms of the patients.
Genetic population structure and phylogeny
In the MLST analysis, a subset of 43 Y. enterocolitica BT 1A strains were discriminated into 39 MLST types and the 10 4/O:3, 3/O:3 or 2/O:9 strains were discriminated into four different MLST types. The genetic diversity indexes for the genes used in MLST were 0.86 (adk), 0.93 (argA), 0.93 (aroA), 0.83 (glnA), 0.82 (gyrB), 0.94 (thrA) and 0.89 (trpE). Bayesian analysis of the MLST sequences divided the BT 1A strains into two distinct genetic clusters, which were clearly separated from the tight cluster formed by the strains of BT’s 2–4 and from the BT 1B strain (8018) (Figure1). One of the BT 1A clusters contained 36 BT 1A and two non-biotypeable strains and was designated as BT 1A Genetic group 1. Another cluster contained five BT 1A strains and was designated as BT 1A Genetic group 2. Ten bio/serotype 3-4/O:3 and 2/O:9 strains clustered closely together, and the single BT 1B strain was located in the vicinity of this cluster. BAPS analysis did not indicate any significant level of mosaicism among the isolates, i.e. no isolates contained variation typical to more than one cluster.
When concatenated MLST sequences (4580 bp) were compared to each other, the BT 1A Genetic group 2 strains were 95–96% similar to BT 1A Genetic group 1, bio/serotype 4/O:3 and 2/O:9, as well as to Y. enterocolitica ssp. enterocolitica strains of biotype 1B (Table1). The BT 1A Genetic group 1 strains were 97% similar to bio/serotype 4/O:3 and 2/O:9 and Y. enterocolitica ssp. enterocolitica strains (Table1). A neighbour-joining tree depicting the relatedness of the selected Yersinia strains and species based on the MLST sequence concatenates is shown in an additional file (Additional file1).
Comparison of the partial 16S rRNA gene sequences (1310 bp) revealed BT 1A Genetic group 2 strains to have over 99% similarity among themselves, 98–99% similarity to Genetic group 1 strains and 98% similarity to Y. enterocolitica BT 2-4/O:3 or O:9 strains (Table2). Actually, the 16S rRNA gene sequences of BT 1A Genetic group 2 were more similar (99%) to Y. intermedia, Y. mollaretii, Y. aldovae and Y. bercovieri than to BT 1A Genetic group 1 (Table2). When the results obtained from representative subsets of 71 strains and analysed using 16S rRNA gene sequencing and MLST were combined, two genetic groups were formed: 17 strains were in Genetic group 2 and 54 in Genetic group 1.
Of all the BT 1A Genetic group 1 strains included in the MLST analysis, none were ystA positive in PCR, but 98% were ystB positive. All five of the BT 1A Genetic group 2 strains were both ystA and ystB negative in PCR. The 4/O:3, 3/O:3 and 2/O:9 strains were all ystA positive and ystB negative in PCR. When also the BT 1A strains that were not included in the MLST analysis were tested for ystA and ystB, 12 further strains were found to be negative in ystB PCR. They were also subjected to 16S rRNA gene sequencing and were found to be part of BT 1A Genetic group 2 (Figure2).
Based on the characteristics of the lipopolysaccrarides (LPS) in silver-stained DOC-PAGE gels, the 298 Y. enterocolitica BT 1A strains were classified into four main LPS types (A-D), with each containing several subtypes (Table3). The subtype characteristics are described in detail in an additional file (Additional file2). Fourteen of the BT 1A Genetic group 2 strains expressed homopolymeric O-PS subtype A2 LPS, similar to LPS of serotype O:10 (Table3). Three of the BT 1A genetic group 2 strains had rough LPS type. The strains of Genetic group 1 and the non-typeable strains expressed a great variety of LPS types and subtypes. All 77 Y. enterocolitica 4/O:3 and 3/O:3 strains included in the LPS analysis expressed homopolymeric subtype A3 O-PS and the five Y. enterocolitica 2/O:9 strains subtype A2 O-PS (Table3). Three of the ystB negative strains of BT 1A Genetic group 1 belonged to LPS group A2, two to C1 and one to B1c.
Phage sensitivity of the strains was tested using five Yersinia specific bacteriophages (Table4). Most of the 63 bio/serotype 3–4/O:3 strains were sensitive to ϕYeO3-12, PY100 and ϕR1-RT, in addition 7 strains were sensitive to ϕ80-81. The single bio/serotype 2/O:9 strain was infected by ϕR1-RT only. The 273 BT 1A and non-biotypeable strains representing different LPS-types showed variable phage sensitivity patterns further demonstrating the heterogeneity of this group of strains. However, all 17 of the BT 1A Genetic subgroup 2 strains were resistant to all the tested phages.
Resistance to human serum complement-mediated killing was most common (99%) in the LPS subtype A3 strains, which included the known pathogenic Y. enterocolitica serotype O:3 strains (Table5). Of the strains in the LPS subtype C2, which included the BT 1A/O:5 isolates, 87% were serum resistant. Serum resistance was also high (67%) among subtype C1 strains, which included BT 1A strains with similar LPS-structure to reference strains of serotypes O:6, O:6,30 and O:6,31. Of the BT 1A LPS subtype A2 (O:10) strains, 72% showed resistance to complement killing. However, 13 of the 14 (93%) BT 1A Genetic group 2 strains among the LPS subtype A2 showed high resistance to complement killing. As a whole, 14 of the 17 (82%) strains of the BT 1A Genetic group 2 were resistant to serum complement killing (Figure2). Among the LPS B-subtypes, which included a number of the BT1A Genetic group1 isolates, complement resistance was rather low or non-existing (Table5).
Statistical analysis of patient symptoms
The symptoms (diarrhoea, vomiting, fever, abdominal pain and blood in stools) of patients with BT 1A did not differ significantly when the statistical analyses were based on the genetic grouping or serum resistance of the BT 1A isolates. The patients with isolates belonging to different LPS-groups were symptomatic, but due to the small amount of patients in analyses, no significant statistical inference could be made.
The strains previously identified by phenotypic tests to belong to Y. enterocolitica BT 1A formed two phylogenetic clusters based on MLST analysis, Genetic groups 1 and 2. BT 1A Genetic group 1 comprised of isolates with related 16S rRNA gene sequences but with great variation in their pathogenicity-associated properties. On the contrary, BT 1A Genetic group 2 was found to be rather uniform and phylogenetically distinct from the other Y. enterocolitica BT 1A strains. The genetic similarity of this group to Genetic group 1 was 95–96% based on the MLST sequences and 98–99% based on the 16S rRNA gene sequences. All the 17 strains determined to belong to Y. enterocolitica BT 1A Genetic group 2 were ystB negative in PCR and were resistant to the five tested yersiniophages. Additionally, none of them fermented fucose, as determined in our previous study. Likewise, pathogenic pYV + yersinia strains do not ferment fucose, whilst 91% of the BT 1A strains other than those of Genetic group 2 do. Of the Genetic group 2 strains 82% were resistant to serum complement killing and 76% belonged to LPS type A2.
Remarkably, the 16S rRNA sequences of BT 1A Genetic group 2 were more similar to Y. intermedia, Y. mollaretii, Y. aldovae and Y. bercovieri than to Y. enterocolitica 16S rRNA sequences. However, a previous study indicated that the use of MLST of house-keeping genes determined genetic relatedness among Yersiniae better than 16S rRNA. Studies using both DNA hybridization and 16S rRNA gene sequence data have illustrated that if two strains show less than 97% 16S rRNA gene sequence similarity, they are separate species. Nevertheless, even 99% similarity of 16S rRNA genes does not guarantee that bacterial strains represent the same species. Howard and colleagues have already suggested that BT 1A strains should be designated as a third subspecies of Y. enterocolitica based on the comparison of whole genomes using DNA microarray. It is likely that the genetic difference between the two phylogenetic groups of Y. enterocolitica BT 1A discovered in the present study may also be high enough to justify designation of different subspecies or even species. Although further analyses would be needed for species designation, our data add insight into the phylogeny of the genus Yersinia, which is continuously evolving: three novel Yersinia species, Y. entomophaga, Y. pekkanenii and Y. nurmii were described as recently as 2010[31–33].
This is the first time that two phylogenetic clusters of Y. enterocolitica BT 1A strains are reported based on the sequence analysis of house-keeping genes, but similar results indicating the existence of two main clusters of BT 1A strains have been obtained with other molecular methods, such as ribotyping and REP-ERIC, gyrB-RFLP, AFPL, MLEE[23, 24] and, most recently, MALDI-TOF mass spectrometry to identify the protein mass patterns. Gulati and Virdi found in their study that sequences of gyrB genes of representative two clonal groups of BT 1A strains were only 97% similar. When we compared these gyrB sequences to our data, sequences DQ140396 and DQ140397 were clustered with BT 1A Genetic groups 1 and 2 of our study, respectively. This is further justification for the separation of BT 1A strains into two phylogenetic lineages. As in our study, the presence of ystB gene correlated with the clonal groups, except in one strain. The lack of the ystB gene in PCR test does not always correlate with the phylogenetic lineages, since our study also found six strains without the ystB gene in BT 1A Genetic group 1. However, only the use of hybridization analysis or sequencing would confirm the PCR results. In a recent study of the whole genome sequences no evident structural difference was found with ystB-positive BT 1A/O:5 and BT 1A/O:36 strains. Therefore, it is likely that the two whole genome sequences represent one of the genetic groups of BT 1A of the present study. Blast searches showed that the sequences we obtained for Genetic group 1 were nearly identical with the ones from the above mentioned whole genome sequences, while for Genetic group 2 no matching sequences were detected.
We used DOC-PAGE based classification of LPS to subtype our Y. enterocolitica strains. This method offered a practical substitute for O-serotyping, since there are no commercial O-specific antisera available for numerous Y. enterocolitica serotypes. The results were consistent with earlier O-serotyping of the BT 1A strains using available commercial antisera which demonstrated that 42 subtype C2 strains were of serotype O:5 and that 56 subtype B2 strains agglutinated with anti-O:8 antiserum indicating that they probably were of the common serotype O:7,8. However, the strains with O:8 antigen, were found in LPS subgroups B2c and B2d which indicates that the classification of subgroups of B2 was tentative and differences could also be inherent to the silver staining procedure. The clinical BT 1A strains showed a wide diversity in their LPS types and this is most likely also reflected in their O-serotypes. The majority of the strains, 37%, had LPS subtype C1 that is similar to that of serotypes O:6,30 and O:6,31, and 15% of the strains had subtype C2, i.e., that of serotype O:5. Globally, the serotypes O:6 and O:5 have been the dominant serotypes of BT 1A associated with diarrhoea. In the present study the strains of LPS subtype C1 and C2 as well as the strains of BT 1A Genetic group 2, demonstrated significant resistance to complement killing, which suggests that the strains of these subgroups may have more pathogenic potential than the other studied strains. Bacterial pathogens have several strategies to resist host defence mechanisms, including resistance to the bactericidal activity of the human serum complement. Pathogenic Y. enterocolitica 4/O:3 strains are able to resist serum killing by YadA- and Ail-mediated binding of the serum complement regulatory proteins factor H and C4 binding protein[36–38]. The BT 1A strains of the present study did not possess the yadA gene and only one strain had the ail gene. Thus the resistance to complement killing of these BT 1A strains must have another, unresolved mechanism.
Although the potential pathogenicity of BT 1A strains remains controversial, there are a few studies that show an association to disease. For instance, BT 1A/O:6,30 was associated with spondyloarthropaties of patients in England and South-Wales. Also, in a study of antibody production, it was found that a patient with symptoms of diarrhoea and reactive arthritis had IgG, IgA and IgM antibodies against the BT 1A/O:6 strain isolated from her fecal sample. We found symptomatic patients with isolates of both BT 1A genetic groups, but did not find statistical differences between the genetic groups and the clinical picture of the symptoms of these patients. It may be that the patients’ genetic or other factors such as gut environment are relevant in the disease caused by BT 1A strains.
The results of our study present strong evidence that strains classified as Y. enterocolitica BT 1A represent more than one subspecies. BT 1A Genetic group 1 consisted of strains with a variety of pathogenicity-related properties, whereas all 17 strains of BT 1A Genetic group 2 lacked the ystB gene, belonged either to the same LPS subtype or were rough, were all resistant to the five tested yersiniophages and were largely resistant to serum complement killing. Furthermore, none of them fermented fucose. Although several studies have been conducted to reveal the significance of the BT 1A strains in causing disease, indisputable results have not been obtained. This study shows, however, that BT 1A is a very heterogenous group of strains, some of which might be potential pathogens. Therefore, better understanding of the genetic and phenotypic variability and clustering of these strains, as achieved in our study, would be crucial in determining the pathogenic role of the strains belonging to the defined clusters.
Altogether 298 BT 1A, 75 bioserotype 4/O:3, two 3/O:3, five 2/O:9 and two non-biotypable Y. enterocolitica strains isolated in 2006 from human samples were utilized in the study. Only one strain per person was included in the study.
MLST analysis was done on 53 Y. enterocolitica strains (43 BT 1A and 10 BT’s 2–4 strains) that represented various LPS patterns. Additionally, two reference strains, NCTC11174 (O:9) and NCTC11176 (O:3), were included in the analysis. Genomic DNA was extracted using Jetflex Genomic DNA purification kit (Genomed, Löhne, Germany). Fragments of seven house-keeping genes (adk, argA, aroA, glnA, gyrB, thrA, trpE) were amplified by PCR. For the adk, argA, aroA, glnA, thrA and trpE genes, the primers available in the MLST database for Y. pseudotuberculosis at the ERI, University College Cork, were used (http://mlst.ucc.ie/mlst/dbs/Ypseudotuberculosis/documents/primersPseudotuberculosis_html). The gyrB gene amplification was done with the primers described earlier. The 25 μl amplification reactions consisted of 0.25 μM of primers, 0.2 mM dNTP, 2.5 U AmpliTaq Gold (Applied Biosystems, Foster City, USA) and 10 × buffer supplied with the enzyme. The thermal cycle consisted of 10 min denaturation at 94°C, followed by 35 cycles of denaturation for 30 s at 94°C, annealing for 30 s at 51°C, and elongation for 30 s at 72°C and finally for 3 min at 72°C. The PCR fragments were sequenced in both directions with an ABI 3730xl DNA Analyzer (Applied Biosystems). The Diversity indexes for each MLST gene were calculated by eBURST v3[40, 41]. The MLST sequences of 53 Y. enterocolitica strains obtained in the study were deposited to EMBL/GenBank database under the accession numbers HE803367- HE803737.
Analysis of the MLST data
Population genetic analyses were performed using BAPS (Bayesian Analysis of Population Structure) software[42–44] with the second-order Markov model and the standard MLST data input option as in, e.g.,[45, 46]. The optimal number of clusters was calculated using 10 runs of the estimation algorithm with the prior upper bound of the number of clusters varying in the range (5,15) over the 10 replicates. All estimation runs resulted in an identical partition of the sequence type data with 4 clusters (estimated P = 1.000). Admixture analysis was done using 100 Monte Carlo replicates for allele frequencies and by generating 100 reference genotypes to calculate p-values. For reference cases we used 10 iterations in estimation according to the guidelines of[44, 47].
Mosaicism is defined as sequence types composed of sequence characteristic of more than one BAPS group. Significance of admixture or mosaicism was determined for each sequence type using the threshold p < 0.05.
Maximum likelihood tree was constructed by using the concatenated sequences under the general time-reversible model as implemented in the MEGA5 software.
16S RNA gene sequencing and tree construction
16S rRNA gene sequencing was obtained for 36 Y. enterocolitica BT 1A strains with the primers FD1mod, pHr, pDf, and pEr in conditions described earlier. The sequences were used to construct a Neighbour-joining tree using Phylip. The 16S rRNA gene sequences of 28 Y. enterocolitica BT 1A strains obtained during this study were deposited to the EMBL/GenBank database under the accession numbers HE803738 - HE803765. Eight of the BT 1A strains were sequenced during our previous studies and have accession numbers FM958217 - FM958223 and FN812721.
ystA and ystB PCR
For the ystA gene specific PCR the forward primer 3-ATC GAC ACC AAT AAC CGC TGAG −5 and reverse primer 3- CCA ATC ACT ACT GAC TTC GGCT −5 were used for 38 Y. enterocolitica BT 1A strains. For the ystB gene specific PCR the forward primer 3-GTA CAT TAG GCC AAG AGA CG −5 and the reverse primer 3-GCA ACA TAC CTC ACA ACA CC −5 were used for 298 Y. enterocolitica BT 1A strains. Chromosomal DNA was used as a template; the conditions for PCR amplification were as described earlier.
DOC-PAGE analysis of LPS
LPS samples of 298 Y. enterocolitica BT 1A strains were prepared by the small scale proteinase K method as described earlier. Briefly, the bacteria were grown for 14–16 h with shaking in 2 ml of LB at 22°C (RT); the OD600 was determined, the bacteria were then pelleted by centrifugation, and the pellet was re-suspended in DOC lysis buffer (2% DOC, 4% 2-mercaptoethanol, 10% glycerol and 0.002% bromophenol blue in 1 M Tris–HCl buffer, pH 6.8) in a volume adjusted according to the density of the culture (i.e., 100 μl / OD600 =1). The suspension was heated to 100°C for 10 min and then 2–4 μl of proteinase K (20 mg/ml) was added and the suspension was incubated overnight at 60°C. An aliquot of 10 μl was loaded on the gel and analysed in 12% DOC-PAGE and the LPS bands were visualized by silver staining as described earlier. The DOC-PAGE-based LPS classification of Y. enterocolitica and Y. enterocolitica –like bacteria has been described elsewhere. Briefly, based on the O-polysaccharide (O-PS) the strains are classified into four main LPS types: (i) type A, LPS with homopolymeric O-PS, (ii) type B, LPS with ladder-forming heteropolymeric O-PS, (iii) type C, LPS with single-length O-PS, and (iv) type D, rough or semi-rough LPS without O-PS or with a lipid A core substituted with a single O-repeat unit, respectively.
Phage sensitivity assay
The following bacteriophages were used in the typing scheme: фR1–37[57, 58] that infects Yersinia expressing the outer core hexasaccharide in LPS; PY100 that infects a broad range of Yersinia strains; фYeO3–12 that uses the Y. enterocolitica serotype O:3 O-PS as receptor[60, 61]; ϕR1-RT that is a bacteriophage originating from the sewage of Turku, Finland and infects Y. enterocolitica serotype O:3 grown at RT (Skurnik, unpublished); and ф80–18 that is a serotype O:8 O-PS specific phage. For altogether 273 Y. enterocolitica BT 1A strains, a 40 μl aliquot from a bacterial culture grown for 14–16 h at RT or 37°C with shaking in LB was mixed with 3.5 ml of molten 0.4% soft agar adjusted to 50°C, mixed briefly and poured on an LA plate. After the soft agar had solidified, 10 μl drops of different phage suspensions (~105 plaque forming units ml-1) were pipetted onto the surface and the plates were incubated at RT or 37°C 14–16 h. Phage sensitivity was scored as a clear lysis zone in the soft agar.
Complement killing assay
Blood was obtained from healthy human donors who were devoid of anti-Yersinia antibodies. Sera were pooled and stored in aliquots at −70°C. The killing assay for 298 Y. enterocolitica BT 1A strains was performed essentially as described previously. Briefly, for bactericidal assay, bacteria were grown to stationary phase overnight in 5 ml of MedECa (MedE: 0.1 g l-1 MgSO4 × 7H2O, 2 g l-1 citric acid, 10 g l-1 K2HPO4, 3.5 g l-1 NaNH4HPO4 × 4H2O and 1 mg l-1 vitamin B1, supplemented with 0.2% glucose, 0.2% casamino acids and 2.5 mM CaCl2) at 37°C without shaking. The cultures were diluted 1:1000–5000 into PBS to obtain a suspension of ca. 105 cfu/ml and 10 μl of the suspension was mixed with 20 μl of normal human serum (NHS) or heat-inactivated serum (HIS, 30 min at 56°C). After 60 min incubation at 37°C, the complement reaction was stopped by transferring the tubes on ice and the addition of 70 μl of ice-cold BHI. Aliquots of 20 μl were cultured on LA-plates and the surviving bacteria were counted after 48 hr incubation at RT. The serum bactericidal effect was calculated as the survival percentage taking the bacterial counts obtained with bacteria incubated in HIS as 100%. The survival was scored as follows: >50% survival, +++; 5–50% survival, ++; and 0.01–5% survival, +; and no colonies, 0.
Statistical analysis of the symptoms of the patients
We compared the symptoms of diarrhoea, vomiting, fever, abdominal pain and blood in stools among 98 patients with a Y. enterocolitica BT 1A isolate, who had answered a questionnaire about the symptoms and had less than six weeks from the onset of symptoms to the sample-taking. Comparisons (Fischer’s exact test) were done among these patients separately for BT 1A genetic groups 1 and 2 (n = 94 and n = 4); for LPS groups: A1-A3 (n = 5), B1-B4 (n = 41), C1 (n = 37), C2 (n = 10), D1 (n = 5); and for serum resistance groups (n = 46 and n = 52). Analyses were done with STATA 9.0.
Informed consent was obtained from the patients who participated in the questionnaire study. The study was approved by the Ethics Committee of National Institute for Health and Welfare (THL). The voluntary healthy blood donors whose sera were used in serum-killing assay gave their verbal consent. They were informed of the details of the study and their blood samples were pooled and used for the study without an individual being identified.
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.
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.
Ratnam S, Mercer E, Picco B, Parsons S, Butler R: A nosocomial outbreak of diarrheal disease due to Yersinia enterocolitica serotype 0:5, biotype 1. J Infect Dis. 1982, 145: 242-247. 10.1093/infdis/145.2.242.
Greenwood MH, Hooper WL: Excretion of Yersinia spp. associated with consumption of pasteurized milk. Epidemiol Infect. 1990, 104: 345-350. 10.1017/S0950268800047361.
Ebringer R, Colthorpe D, Burden G, Hindley C, Ebringer A: Yersinia enterocolitica biotype I. Diarrhoea and episodes of HLA B27 related ocular and rheumatic inflammatory disease in South-East England. Scand J Rheumatol. 1982, 11: 171-176. 10.3109/03009748209098186.
Skurnik M, Nurmi T, Granfors K, Koskela M, Tiilikainen AS: Plasmid associated antibody production against Yersinia enterocolitica in man. Scand J Infect Dis. 1983, 15: 173-177.
Huovinen E, Sihvonen L, Virtanen M, Haukka K, Siitonen A, Kuusi M: Symptoms and sources of Yersinia enterocolitica -infection: a case–control study. BMC Infect Dis. 2010, 10: 122-131. 10.1186/1471-2334-10-122.
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.
McNally A, Dalton T, Ragione RML, 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.
Singh I, Virdi JS: Interaction of Yersinia enterocolitica biotype 1A strains of diverse origin with cultured cells in vitro. Jpn J Infect Dis. 2005, 58: 31-33.
Nair GB, Takeda Y: The heat-stable enterotoxins. Microb Pathog. 1998, 24: 123-131. 10.1006/mpat.1997.0177.
Delor I, Cornelis GR: Role of Yersinia enterocolitica Yst toxin in experimental infection of young rabbits. Infect Immun. 1992, 60: 4269-4277.
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.
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.
Mikulskis AV, Delor I, Thi VH, Cornelis GR: Regulation of the Yersinia enterocolitica enterotoxin Yst gene. Influence of growth phase, temperature, osmolarity, pH and bacterial host factors. Mol Microbiol. 1994, 14: 905-915. 10.1111/j.1365-2958.1994.tb01326.x.
Kuehni-Boghenbor K, On SL, Kokotovic B, Baumgartner A, Wassenaar TM, Wittwer M, Bissig-Choisat B, Frey J: Genotyping of human and porcine Yersinia enterocolitica, Yersinia intermedia, and Yersinia bercovieri strains from Switzerland by amplified fragment length polymorphism analysis. Appl Environ Microbiol. 2006, 72: 4061-4066. 10.1128/AEM.01996-05.
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.
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.
Falcao JP, Falcao DP, Pitondo-Silva A, Malaspina AC, Brocchi M: Molecular typing and virulence markers of Yersinia enterocolitica strains from human, animal and food origins isolated between 1968 and 2000 in Brazil. J Med Microbiol. 2006, 55: 1539-1548. 10.1099/jmm.0.46733-0.
Bhagat N, Virdi JS: The Enigma of Yersinia enterocolitica biovar 1A. Crit Rev Microbiol. 2011, 37: 25-39. 10.3109/1040841X.2010.506429.
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.
Gulati PS, Virdi JS: The rrn locus and gyrB 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.
Mallik S, Virdi JS: Genetic relationships between clinical and non-clinical strains of Yersinia enterocolitica biovar 1A as revealed by multilocus enzyme electrophoresis and multilocus restriction typing. BMC Microbiol. 2010, 10: 158-10.1186/1471-2180-10-158.
Dolina M, Peduzzi R: Population genetics of human, animal, and environmental Yersinia strains. Appl Environ Microbiol. 1993, 59: 442-450.
Stephan R, Cernela N, Ziegler D, Pfluger V, Tonolla M, Ravasi D, Fredriksson-Ahomaa M, Hächler H: Rapid species specific identification and subtyping of Yersinia enterocolitica by MALDI-TOF mass spectrometry. J Microbiol Methods. 2011, 87: 150-153. 10.1016/j.mimet.2011.08.016.
Batzilla J, Heesemann J, Rakin A: The pathogenic potential of Yersinia enterocolitica 1A. Int J Med Microbiol. 2011, 301: 556-571. 10.1016/j.ijmm.2011.05.002.
Sihvonen LM, Haukka K, Kuusi M, Virtanen MJ, Siitonen A: Yersinia enterocolitica and Y. enterocolitica-like species in clinical stool specimens of humans: identification and prevalence of bio/serotypes in Finland. Eur J Clin Microbiol Infect Dis. 2009, 28: 757-765. 10.1007/s10096-008-0696-y.
Fredriksson-Ahomaa M, Cernela N, Hachler H, Stephan R: Yersinia enterocolitica strains associated with human infections in Switzerland 2001–2010. Eur J Clin Microbiol Infect Dis. 2012, 31: 1543-1550. 10.1007/s10096-011-1476-7.
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.
Staley JT: The bacterial species dilemma and the genomic-phylogenetic species concept. Phil Trans Roy Soc Lond B Biol Sci. 2006, 361: 1899-1909. 10.1098/rstb.2006.1914.
Murros-Kontiainen AE, Fredriksson-Ahomaa M, Korkeala H, Johansson P, Rahkila R, Björkroth J: Yersinia nurmii sp. nov. Int J Syst Evol Microbiol. 2011, 61: 2368-2372. 10.1099/ijs.0.024836-0.
Murros-Kontiainen AE, Johansson P, Niskanen T, Fredriksson-Ahomaa M, Korkeala H, Björkroth J: Yersinia pekkanenii sp. nov. Int J Syst Evol Microbiol. 2011, 61: 2363-2367. 10.1099/ijs.0.019984-0.
Hurst MR, Becher SA, Young SD, Nelson TL, Glare TR: Yersinia entomophaga sp. nov. isolated from the New Zealand grass grub Costelytra zealandica. Int J Syst Evol Microbiol. 2011, 61: 844-849. 10.1099/ijs.0.024406-0.
Bhagat N, Virdi J: 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.
Lambris JD, Ricklin D, Geisbrecht BV: Complement evasion by human pathogens. Nat Rev Microbiol. 2008, 6: 132-142. 10.1038/nrmicro1824.
Biedzka-Sarek M, Jarva H, Hyytiainen H, Meri S, Skurnik M: Characterization of complement factor H binding to Yersinia enterocolitica serotype O:3. Infect Immun. 2008, 76: 4100-4109. 10.1128/IAI.00313-08.
Biedzka-Sarek M, Salmenlinna S, Gruber M, Lupas AN, Meri S, Skurnik M: Functional mapping of YadA- and Ail-mediated binding of human factor H to Yersinia enterocolitica serotype O:3. Infect Immun. 2008, 76: 5016-5027. 10.1128/IAI.00314-08.
Kirjavainen V, Jarva H, Biedzka-Sarek M, Blom AM, Skurnik M, Meri S: Yersinia enterocolitica serum resistance proteins YadA and Ail bind the complement regulator C4b-binding protein. PLoS Pathog. 2008, 4: e1000140-10.1371/journal.ppat.1000140.
Sihvonen LM, Hallanvuo S, Haukka K, Skurnik M, Siitonen A: The ail gene is present in some Yersinia enterocolitica biotype 1A strains. Foodborne Pathog Dis. 2011, 8: 455-457. 10.1089/fpd.2010.0747.
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.
Spratt BG, Hanage WP, Li B, Aanensen DM, Feil EJ: Displaying the relatedness among isolates of bacterial species – the eBURST approach. FEMS Microbiol Lett. 2004, 241: 129-134. 10.1016/j.femsle.2004.11.015.
Corander J, Tang J: Bayesian analysis of population structure based on linked molecular information. Math Biosci. 2007, 205: 19-31. 10.1016/j.mbs.2006.09.015.
Corander J, Marttinen P, Siren J, Tang J: Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinforma. 2008, 9: 539-10.1186/1471-2105-9-539.
Tang J, Hanage WP, Fraser C, Corander J: Identifying currents in the gene pool for bacterial populations using an integrative approach. PLoS Comput Biol. 2009, 5: e1000455-10.1371/journal.pcbi.1000455.
Hanage WP, Fraser C, Tang J, Connor TR, Corander J: Hyper-recombination, diversity, and antibiotic resistance in pneumococcus. Science. 2009, 324: 1454-1457. 10.1126/science.1171908.
Corander J, Connor RR, O’Dwyer CA, Kroll JS, Hanage WP: Population structure in the Neisseria, and the biological significance of fuzzy species. J Royal Soc Interface. 2012, 9: 1208-1215. 10.1098/rsif.2011.0601.
Corander J, Marttinen P: Bayesian identification of admixture events using multilocus molecular markers. Mol Ecol. 2006, 15: 2833-2843. 10.1111/j.1365-294X.2006.02994.x.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Kotilainen P, Jalava J, Meurman O, Lehtonen OP, Rintala E, Seppälä OP, Eerola E, Nikkari S: Diagnosis of meningococcal meningitis by broad-range bacterial PCR with cerebrospinal fluid. J Clin Microbiol. 1998, 36: 2205-2209.
Edwards U, Rogall T, Blöcker H, Emde M, Böttger EC: Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA. Nucleic Acids Res. 1989, 17: 7843-7853. 10.1093/nar/17.19.7843.
Felsenstein J: PHYLIP (Phylogeny Inference Package). 3.6a3. 2001, Department of Genome Sciences, University of Washington, Seattle
Thoerner P, Bin Kingombe CI, Bogli-Stuber K, Bissig-Choisat B, Wassenaar TM, Frey J, Jemmi T: PCR detection of virulence genes in Yersinia enterocolitica and Yersinia pseudotuberculosis and investigation of virulence gene distribution. Appl Environ Microbiol. 2003, 69: 1810-1816. 10.1128/AEM.69.3.1810-1816.2003.
Ramamurthy T, Yoshino K, Huang X, Balakrish Nair G, Carniel E, Maruyamad T, Fukushimae H, Takeda T: The novel heat-stable enterotoxin subtype gene (ystB) of Yersinia enterocolitica: nucleotide sequence and distribution of the yst genes. Microb Pathog. 1997, 23: 189-200. 10.1006/mpat.1997.0146.
Bengoechea JA, Zhang L, Toivanen P, Skurnik M: Regulatory network of lipopolysaccharide O-antigen biosynthesis in Yersinia enterocolitica includes cell envelope-dependent signals. Mol Microbiol. 2002, 44: 1045-1062. 10.1046/j.1365-2958.2002.02940.x.
Skurnik M, Venho R, Toivanen P, al-Hendy A: A novel locus of Yersinia enterocolitica serotype O:3 involved in lipopolysaccharide outer core biosynthesis. Mol Microbiol. 1995, 17: 575-594. 10.1111/j.1365-2958.1995.mmi_17030575.x.
Skurnik M, Toivonen S: Identification of distinct lipopolysaccharide patterns among Yersinia enterocolitica and Y. enterocolitica -like bacteria. Biochemistry (Mosc). 2011, 76: 823-831. 10.1134/S0006297911070133.
Kiljunen S, Hakala K, Pinta E, Huttunen S, Pluta P, Gador A, Lönnberg H, Skurnik M: Yersiniophage phiR1-37 is a tailed bacteriophage having a 270 kb DNA genome with thymidine replaced by deoxyuridine. Microbiol. 2005, 151: 4093-4102. 10.1099/mic.0.28265-0.
Skurnik M: Role of YadA in Yersinia-enterocolitica-induced reactive arthritis: a hypothesis. Trends Microbiol. 1995, 3: 318-319. 10.1016/S0966-842X(00)88961-1.
Schwudke D, Ergin A, Michael K, Volkmar S, Appel B, Knabner D, Konietzny A, Strauch E: Broad-host-range Yersinia phage PY100: genome sequence, proteome analysis of virions, and DNA packaging strategy. J Bacteriol. 2008, 190: 332-342. 10.1128/JB.01402-07.
Al-Hendy A, Toivanen P, Skurnik M: The effect of growth temperature on the biosynthesis of Yersinia enterocolitica O:3 lipopolysaccharide: temperature regulates the transcription of the rfb but not of the rfa region. Microb Pathog. 1991, 10: 81-86. 10.1016/0882-4010(91)90068-L.
Pajunen M, Kiljunen S, Skurnik M: Bacteriophage phiYeO3-12, specific for Yersinia enterocolitica serotype O:3, is related to coliphages T3 and T7. J Bacteriol. 2000, 182: 5114-5120. 10.1128/JB.182.18.5114-5120.2000.
Zhang L, Skurnik M: Isolation of an R- M + mutant of Yersinia enterocolitica serotype O:8 and its application in construction of rough mutants utilizing mini-Tn5 derivatives and lipopolysaccharide-specific phage. J Bacteriol. 1994, 176: 1756-1760.
Biedzka-Sarek M, Venho R, Skurnik M: Role of YadA, Ail, and lipopolysaccharide in serum resistance of Yersinia enterocolitica serotype O:3. Infect Immun. 2005, 73: 2232-2244. 10.1128/IAI.73.4.2232-2244.2005.
We wish to acknowledge the excellent technical assistance of Heini Flinck, Tarja Heiskanen, Katriina Mälkönen and Ahmed Mohammed Ahmed. Harri Sihvonen is thanked for assistance with figure preparation. This work was supported by a grant (4850/501/2004) from the Finnish Ministry of Agriculture and Forestry.
The authors’ declare that they have no competing interests.
LMS conducted the MLST work, combined all the results together and drafted the manuscript. KJ contributed to the genomic analyses. ST and MS conducted and analyzed the LPS, serum resistance and phage typing assays. EH and MK analysed the clinical data and JC did the BAPS and phylogenetic analysis of the MLST data. AS and KH participated in planning of the work, analyzing the results and writing the article. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Neighbour-joining tree based on seven concatenated MLST genes (4580 bp). Neighbour-joining bootstrap confidence values over 75% (1000 replicates) are given in the branches. BT 1A strains were ystB positive in PCR and had positive reaction in fucose fermentation unless otherwise indicated. sr=serum resistance; pt= phage type, which encodes reaction to 5 phages (φR1–37, PY100, φYeO3–1, φR1-RT, φ80–81). Strains sequenced in the present study are marked bold. In addition, the following GenBank sequences were used: Y. enterocolitica 8081 (AM286415), Y. aldovae ATCC 35236 (ACCB00000000), Y. kristensenii ATCC 33638: (ACCA00000000), Y. intermedia ATCC 29909 (AALF00000000), Y. frederiksenii ATCC 33641 (AALE00000000), Y. mollaretii ATCC 43969 (AALD00000000), Y. bercovieri ATCC 43970 (AALC00000000), Y. rohdei ATCC 43380 (ACCD00000000) and Y. ruckeri ATCC 29473 (ACCC00000000). (DOCX 649 KB)
Additional 2: Analysis of Y. enterocolitica LPS by DOC-PAGE and silver staining. The picture is compiled of gel images with different LPS types as indicated above the lanes by the LPS type codes that are explained in the text box. Please note that LPS types A2, B1c, B1d, B2a, B2c and B4 are not shown. The gel regions where O-PS and lipid A (LA) plus core migrate are indicated by arrows. (DOCX 230 KB)
About this article
Cite this article
Sihvonen, L.M., Jalkanen, K., Huovinen, E. et al. Clinical isolates of Yersinia enterocolitica Biotype 1A represent two phylogenetic lineages with differing pathogenicity-related properties. BMC Microbiol 12, 208 (2012). https://doi.org/10.1186/1471-2180-12-208