The strains previously identified by phenotypic tests to belong to
Y. enterocolitica BT 1A
[
27] 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
[
27]. 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
[
29]. 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
[
30]. Nevertheless, even 99% similarity of 16S rRNA genes does not guarantee that bacterial strains represent the same species. Howard and colleagues
[
17] 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
[
21],
gyrB-RFLP
[
22], AFPL
[
16], MLEE
[
23,
24] and, most recently, MALDI-TOF mass spectrometry to identify the protein mass patterns
[
25]. Gulati and Virdi
[
22] 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
[
22] 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
[
34]. 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
[
26]. 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
[
27] 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
[
20]. 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
[
35]. 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
[
39]. 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
[
5]. 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
[
6]. 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.