# Patent

## Title
Combination of marker genes for characterizing a Lactobacillus sakei strain

## Abstract

The present invention relates to a new combination of marker genes for characterizing a Lactobacillus sakei strain. In particular, the present invention concerns the use of a pattern of presence or absence of marker genes in the genome of the strain to be characterized for classifying and identifying said strain.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is a U.S. National Stage application of international application PCT/EP2009/052075, filed in English on Feb. 20, 2009, which designated the United States, and which claims the benefit of European patent application 08300118.0, filed Feb. 27, 2008 and U.S. provisional application 61/031,763, filed Feb. 27, 2008. Each of these applications is incorporated by reference herein in its entirety.
FIELD OF THE INVENTION
The present invention relates to a combination of marker genes for characterizing a Lactobacillus sakei strain. In particular, the present invention concerns the use of a pattern of presence or absence of marker genes in the genome of the strain to be characterized for classifying and identifying said strain.
BACKGROUND OF THE INVENTION
Fresh meat and fish are nutritious but highly perishable foods. During production and storage they are exposed to unavoidable microbial contamination from the processing environment. Such contamination may include spoilage organisms and pathogens. It is therefore a priority for food processors to restrict the growth of contaminants so that they do not develop to potentially dangerous levels. One of the methods employed is the use of safe bacteria to curb growth of spoilage and disease-causing micro-organisms. The meat-borne lactic acid bacterium Lactobacillus sakei shows in this view excellent properties.
L. sakei has the ability to survive and grow on fresh meat, forming the dominant population when selective techniques are applied. Some strains are widely used in Europe for artisanal and larger scale manufacture of fermented sausages because of their useful preservative properties. But they could also be used as meat biopreservative by preventing the growth of unwanted bacteria (Vermeiren et al. 2004 Int. J. Food Microbiol. 96:149-164).
Strains of L. sakei can display an important variability in phenotypic traits and have for long been considered difficult to classify. Previous studies have disclosed approaches for L. sakei strains classification. Studies using numerical analysis of RAPD patterns (Berthier and Ehrlich 1999 Int. J. Syst. Bacteriol. 49:997-1007) or SDS-PAGE soluble protein content patterns have suggested the division of the strains into two subgroups, although weakly defined and not comparable according to the studies.
Thus, L. sakei strains classification obtained through these approaches remains unsatisfactory. An aim of the present invention is to provide better methods of classification, characterization and identification of L. sakei strains, which could in particular be used to identify L. sakei strains present on food or in a cocktail of bacteria used as biopreservative.
The inventors of the present invention have herein identified, based on an in silico study of the L. sakei genome and on experimental validation on a large collection of L. sakei strains, a combination of markers for characterizing and detecting these strains. In particular, 29 marker genes have been isolated, whose combination provides an optimal way of characterizing L. sakei strains.
SUMMARY OF THE INVENTION
The present invention thus relates to a combination of markers enabling for characterizing a L. sakei strain which comprises at least two marker genes selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29).
In a preferred embodiment, said combination comprises all of the marker genes SEQ ID NO:1 to SEQ ID NO:29.
The present invention also relates to a method for characterizing a L. sakei strain which comprises the step consisting of determining the presence or absence of at least one marker gene selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29) in said L. sakei strain.
In a preferred embodiment, the presence or absence of all of said marker genes SEQ ID NO: 1 to SEQ ID NO: 29 is determined.
Preferably, the presence or absence of said marker gene(s) is determined by amplification, or by hybridization with probes specific of said marker gene(s).
In a preferred embodiment, a classification of said L. sakei strain is performed by analyzing the pattern of presence or absence of the above marker gene(s) and calculating a Jaccard index or Dice Coefficient or any binary distance matrix with regard to a set of reference L. sakei strains.
The present invention also concerns a DNA array which comprises a combination of markers as defined above.
Another aspect of the present invention concerns a method of comparing at least two L. sakei strains comprising the steps consisting of

    
    
        a) determining the presence or absence of at least one marker gene selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29), in a first L. sakei strain, and
        b) determining the presence or absence of said at least one marker gene in a second L. sakei strain,

wherein, if the pattern of presence or absence of said at least one marker gene is different between said first and said second L. sakei strains, then the L. sakei strains are different.

    
    


In a preferred embodiment, the presence or absence of all of said marker genes SEQ ID NO: 1 to SEQ ID NO: 29 is determined.
In a particular embodiment, the above method of comparing at least two L. sakei strains further comprises the step of analysing the pattern of presence or absence of said marker genes by calculating a Jaccard index, or a Dice coefficient or any binary distance matrix.
Another aspect of the present invention relates to a method of quantifying a specific L. sakei strain in a sample comprising the steps consisting of:
a) identifying and/or characterizing the L. sakei strains present in said sample by a method as defined above,
b) determining the marker gene(s) which is(are) differently present in said specific L. strain with respect to other strains present in said sample, and
c) amplifying said gene(s) by quantitative PCR.



DESCRIPTION OF THE FIGURES
FIG. 1: L. sakei genomic diversity. Distribution of the 60 genetic markers among the 73 L. sakei isolates and the 2 L. curvatus isolates taken as outgroup reference. Genes are ordered on the left of the table by their position in the 23K chromosome, i.e. LSA0088, LSA0118, LSA0157, LSA0165, LSA0172, LSA0178, LSA0212, LSA0216, LSA0217, LSA0218, LSA0219_b, LSA0306, LSA0439, LSA509, LSA510, LSA0530, LSA0564_ac, LSA0565, LSA0567, LSA0572, LSA724, LSA727, LSA1006, LSA1182/3, LSA1220, LSA1222, LSA1227, LSA1232, LSA1283, LSA1509, LSA1510, LSA1512, LSA1510_a, LSA1510_g, LSA1572, LSA1579/80, LSA1581, LSA1584, LSA1641, LSA1640, LSA1640, LSA1720, LSA1724, LSA1730, LSA1731, LSA1806, LSA1809, LSA1874, sspT, sspA, FGP21-0001, DrsA, FGP332-0001, FGP332-0002, FGP332-0003, FGP332-0006, FGP332-0005, lacC, lacG. Colours: white, absent; black, present. The dendrogram showing estimates of genomic relationships of the strains was constructed by complete linkage hierarchical analysis. The scale represents the distance at each node. A coefficient of 1 would denote complete independence, and zero would indicate absolute identity. P-values at node indicate confidence of the clustering by multiscale bootstrap resampling using the PVCLUST program. Strains were grouped on the basis of cluster branches which confidence was above 90% and with a maximum distance of 0.5 between isolates. Main groups of strains are indicated with their respective percentage of the total population. The subspecies type strains and the reference strain L. sakei 23K are indicated by vertical rectangle.
FIG. 2: (A) Histogram distribution of genome size among L. sakei isolates. Strains are represented by rug bars above the genome size axis. Histogram bars are representing the number of strains within a genome size window of 30 kb (average standard deviation of PFGE measurements). The Gaussian probability distribution of genome size in the population (estimated density on the right axis) is shown by the grey smooth line. The subspecies type strains and strain 23K are indicated by asterisks. (B) Boxplot showing genome size distribution of L. sakei isolates according to their genotypic clustering. The horizontal black line is showing the average genome size of 2,020 kb. Clusters E and F showing few variations between them were pooled together for the analysis. For clarity, unclassified strains 332, LTH2070 and 21 were grouped with cluster H.
FIG. 3: Dendrogram showing estimates of genomic relationships of the strains was constructed by complete linkage hierarchical analysis. The scale represents the distance at each node. Strains were grouped on the basis of cluster branches which confidence was above 90% and with a maximum distance of 0.5 between isolates. In the matrix of gene content, a black colour indicates the presence of the gene (1) and a white colour indicates the absence of the gene (0). The genes on the left are ordered as in Table 1. More specifically, the different lines of the matrix correspond to the following gene markers: 1) LSA1641, 2) LSA1182 and LSA1183_c, 3) LSA0172, 4) LSA1731, 5) LSA0211 and LSA0212, 6) LSA1579 and LSA1580, 7) LSA0118, 8) LSA0529, 9) LSA0439, 10) LSA0572, 11) LSA0219b, 12) LSA0564_a, LSA0564_b and LSA0564_c, 13) FGP21-0001, 14) sspA and spiA, 15) FGP332-0001, 16) FGP332-0002, 17) FGP332-007, 18) FGP332-0008, 18) FGP332-0009, 19) FGP332-0010, 20) FGP332-0011, 21) FGP332-0011, 22) FGP332-0012, 23) FGP332-0013.



DETAILED DESCRIPTION OF THE INVENTION
In order to identify marker genes that could be used to characterize specifically a L. sakei strain compared with other L. sakei strains, the inventors first identified variable genomic islands in the genome of L. sakei 23K that were likely to contain variable genes displaying a high diversity according to the species (step described in Example 1). Then, they performed a classification of L. sakei strains according to the presence or absence of genes from genomic islands previously identified and compared the so obtained classification with the one obtained using several other typing techniques (step described in Example 2). Finally, they identified an optimal combination of a minimum number of genes from the genomic islands that could be used to obtain a reliable classification of L. sakei strains compared to the one obtained with a high number of variable genes (step described in Example 3).
Accordingly, according to the invention, a combination of markers enabling for characterizing a L. sakei strain comprises at least two marker genes selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29).
As used herein, the term—marker—refers to any biological, chemical or physical mean allowing identifying the presence, and possibly quantifying the expression of a target gene in a bacterial strain. Such markers are well known from one skilled in the art. Advantageously, the markers according to the invention are gene markers.
Table 1 gives the correspondence between the above referenced locus tag and sequence of marker genes and the name of the genes.







TABLE 1







Correspondence between the locus tag and the name of the gene markers.









SEQ ID NO
Locus Tag
Gene name












1
LSA1641
N-acetylmannosamine-6-phosphate 2-epimerase (N-




acetylmannosamine-6-phosphate 2 epimerase))


2
LSA1182
Putative Cytochrome P450 (authentic frameshifted gene)




C-terminal part


3
LSA1183_c
Putative Cytochrome P450 (authentic frameshifted gene)




N-terminal part


4
LSA0172
CscC-type cell-surface protein with Invasin/Mucin-like




domain and WxL domain


5
LSA1731
CscC-type cell-surface protein with Haemagluttinine-like




domain and WxL domain


6
LSA0211
CscC-type cell-surface protein with adhesion-like domain




and WxL domain (authentic frameshifted gene) N-




terminal part


7
LSA0212
CscC-type cell-surface protein with adhesion-like domain




and WxL domain (authentic frameshifted gene) C-




terminal part


8
LSA1579
Putative teichoic acid/polysaccharide export protein




complex


9
LSA1580
Putative teichoic acid/polysaccharide export protein




complex


10
LSA0118
Hypothetical protein (Putative cell-surface collagen-




binding protein)


11
LSA0529
MarR-type Transcriptional regulator, putative peroxide




stress regulator OhrR


12
LSA0439
Hypothetical extracellular lipase/esterase precursor


13
LSA0572
Threonine deaminase (Threonine ammonia-lyase)


14
LSA0219_b
Putative cyanate transport protein


15
LSA0564_a
Hypothetical small peptide


16
LSA0564_b
Hypothetical small peptide


17
LSA0564_c
Putative bacteriocin Immunity protein


18
FGP21-0001
Putative bacteriocin Immunity protein


19
sspA
Bacteriocin sakacin P precursor (Sakacin 674)


20
spiA
Sakacin P immunity protein


21
FGP332-0001
Putative 6-phospho-beta-glucosidase


22
FGP332-0002
CscC-type cell-surface protein with bacterial adhesion-




like domain and WxL domain


23
FGP332-0007
Putative autotransporter


24
FGP332-0008
Hypothetical protein


25
FGP332-0009
Hypothetical protein


26
FGP332-0010
Hypothetical protein


27
FGP332-0011
Putative transcriptional regulator, LysR family


28
FGP332-0012
Putative quinine oxidoreductase


29
FGP332-0013
Putative asparagines synthase









The presence or absence of one or a combination of some of these genes can be used to discriminate a particular L. sakei strain from another closely related strain.
As described in Example 3, the inventors have demonstrated that these 29 marker genes were optimal to characterize a L. sakei strain among other L. sakei strains. Therefore, in a preferred embodiment, said combination consists of the marker genes SEQ ID NO: 1 to SEQ ID NO: 29.
More precisely, the following correspondence between the pattern of presence or absence of the 29 marker genes defined above and 75 L. sakei strains has been obtained, as illustrated in FIG. 3:

    
    
        in strain DSM20019, the marker gene SEQ ID NO: 21 is present, and the other marker genes are absent.
        in strain CTC494, the marker gene SEQ ID NO: 21 is present, and the other marker genes are absent.
        in strain AGR46, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, and SEQ ID NO: 24 are present and the other marker genes are absent.
        in strain 156, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO:17 and SEQ ID NO: 24 are present and the other marker genes are absent.
        in strain AGR51, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO:18, SEQ ID NO: 21 and SEQ ID NO: 24 are present and the other marker genes are absent.
        in strain MF1048, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO:13, SEQ ID NO: 14, SEQ ID NO: 21, SEQ ID NO: 24 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain MF2092, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO:13, SEQ ID NO: 14, SEQ ID NO: 21, SEQ ID NO: 24 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 195, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 21, and SEQ ID NO: 25 are present and the other marker genes are absent.
        in strain CIP105422, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO; 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain CTC429, the marker genes SEQ ID NO: 1, SEQ ID NO:2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 26, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LTH673, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 25, SEQ ID NO: 26 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain 14, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 14, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 18, the marker genes SEQ ID NO:1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID N: 12, SEQ ID NO: 14, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LV52, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 14, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 26 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain LV59, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 14, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 25, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LV92, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 14, SEQ ID NO: 19, SEQ ID NO: 20, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain CRL1467, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC335, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain SF842, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain SF770, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain SF771, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC287, the marker genes SEQ ID NO: 1, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC041, the marker genes SEQ ID NO: 1, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO:15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC427, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain YME344, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain V553, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16 and SEQ ID NO: 17 are present and the other marker genes are absent.
        in strain LV34, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID N: 15, SEQ ID NO: 16 and SEQ ID NO: 17 are present and the other marker genes are absent.
        in strain 23K, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16 and SEQ ID NO: 17 are present and the other marker genes are absent.
        in strain 72, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16 and SEQ ID NO: 17 are present and the other marker genes are absent.
        in strain 504, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain TISTR911, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain YMW540, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 21, SEQ ID NO: 23 and SEQ ID NO: 24 are present and the other marker genes are absent.
        in strain 300, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21 and SEQ ID NO: 23 are present and the other marker genes are absent.
        in strain LV5, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21 and SEQ ID NO: 23 are present and the other marker genes are absent.
        in strain LTH2070, the marker genes SEQ ID NO: 1, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain JG3, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LTH1764, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain CTC6626, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 18, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain SF841, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 18, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain SF843, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC163, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain LTH5590, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC014, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain 160K, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain AGR48, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain LTH5589, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 112, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain YMH243, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 26, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 134, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain YMW557, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain L110, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LTH5588, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain MF2091, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain MF2089, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain MF2090, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 745, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 710, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 205, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain MF2088, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC494, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain 331, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID NO: 17, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain AGR53, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 11, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain 495, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain Lb706, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LV21, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 4, SEQ ID NO: 5, SEQ ID NO: 11, SEQ ID NO: 13, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain LTH677, the marker genes SEQ ID NO: 5, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID N: 17, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 332, the marker genes SEQ ID NO: 1, SEQ ID NO: 4, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 15, SEQ ID NO: 16, SEQ ID N: 17, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 25, SEQ ID NO: 26, SEQ ID NO: 27, SEQ ID NO: 28 and SEQ ID NO: 29 are present and the other marker genes are absent.
        in strain 33, the marker genes SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 13, SEQ ID NO: 18, SEQ ID NO: 22, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain 64, the marker genes SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24, SEQ ID NO: 27 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain CTC6469, the marker genes SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 12, SEQ ID NO: 13, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 24 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain LTH5728, the marker genes SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 8, SEQ ID NO: 9, SEQ ID NO: 10, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 22, SEQ ID NO: 23, SEQ ID NO: 24 and SEQ ID NO: 27 are present and the other marker genes are absent.
        in strain 21, the marker genes SEQ ID NO: 1, SEQ ID NO: 2, SEQ ID NO: 3, SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 22 and SEQ ID NO: 24 are present and the other marker genes are absent.
        in strain TISTR890, the marker genes SEQ ID NO: 5, SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 18, SEQ ID NO: 21, SEQ ID NO: 24 and SEQ ID NO: 28 are present and the other marker genes are absent.
        in strain ATCC15521, the marker genes SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 21 and SEQ ID NO: 24 are present and the other marker genes are absent.

in strain LTH675, the marker genes SEQ ID NO: 6, SEQ ID NO: 7, SEQ ID NO: 13, SEQ ID NO: 14, SEQ ID NO: 18, SEQ ID NO: 21 and SEQ ID NO: 23 are present and the other marker genes are absent.

    
    


The present invention also relates to a method for characterizing a L. sakei strain which comprises the step consisting of determining the presence or absence of at least one marker gene selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29) in said L. sakei strain.
Numerous methods allowing determining the presence or absence of a gene in a bacterial strain are well known from one skilled in the art. These methods include, without being limited, the use of an antibody specifically binding to an antigen which is constituted by the expression product of said marker gene(s), the detection of mRNA, cDNA or polypeptide from said marker gene(s), or of fragments thereof. Preferably, the presence or absence of said marker gene(s) is determined according to the invention by amplification, or by hybridization with DNA probes specific of said marker gene(s).
More preferably, the presence or absence of all of said marker genes SEQ ID NO:1 to SEQ ID NO:29 is determined.
In a particular embodiment, the combination of probes of these marker genes is assembled on a same support, preferably a standardized support. These supports are known from one skilled in the art. Their size can vary according to the apparatuses used to detect the presence or absence of said marker gene(s).
Advantageously, the combination of marker genes according to the invention is in form of a DNA matrix, comprising a support on which nucleic acids fragments likely to hybridize to target genes are deposed, preferably in a standardized way. The size of such supports varies according to the preparation and detection methods used. Such small supports are also referred to DNA array.
Another aspect of the present invention thus concerns a DNA array which comprises a combination of markers as defined above.
As used herein, the term “DNA array” refers to a set of genes, fragment of genes, oligonucleotides deposited on a support (glass slide, nylon membrane . . . ) with a high density. Numerous scientific publications about the preparation and the use of DNA arrays are available.
In another aspect of the invention, a classification of said L. sakei strain is performed by analyzing the pattern of presence or absence of the above marker gene(s) and calculating a Jaccard index or Dice coefficient or any binary distance matrix with regard to a set of reference L. sakei strains. The Jaccard index can be calculated as described in Jaccard P, 1901 (Jaccard P., Bulletin de la Societe Vaudoise des Sciences Naturelles, 37:241-272). The Dice coefficient can be calculated as described by Van Rijsbergen, 1979 (Van Rijsbergen, 1979, Information Retrieval, London: Butterworths).
As used herein, the term “classification” refers to organizing the strains in different subfamilies according to their genetic pattern. In particular, strains that display the same pattern of genes belong to the same subfamily.
As used herein, the term “reference L. sakei strains” refers to a set or collection of L. sakei strains in which the presence or absence of said marker gene(s) has already been determined according to the invention and to which said L. sakei strain to be classified is compared. Namely said reference L. sakei strains were used to constitute subfamilies or clusters of L. sakei strains in which said L. sakei strain to be classified is intended to be positioned.
Preferably, the set of reference L. sakei strains according to the invention comprises the strains disclosed in Table 2. Nevertheless, other known L. sakei strains can be used as reference strains, and are well known from one skilled in the art.
Preferably, said subfamilies of L. sakei strains constituted by said reference L. sakei strains are as defined in Example 3. However, one skilled in the art knows that clustering can evolve according to the identification of new strains. Accordingly, the use of the method according to the present invention to classify L. sakei strains in clusters which are not specifically described herein is under the scope of the present invention.
As used herein, the term “Jaccard index” is a statistic used for comparing the similarity and diversity of sample sets. The Jacquard index is defined as the size of the intersection divided by the size of the union of the sample sets, according to the following formula: J(A,B)=|A ∩ B|/|A ∪ B|.
Preferably, calculating of the Jaccard Index according to the invention is performed via the software R (R Development Core Team, 2006 R: A language and environment for statistical computing. Vienna, Austria. R Foundation for Statistical Computing) which eventually provides a clustering tree displaying the location of said L. sakei strain among the different subfamilies of reference L. sakei strains.
As used herein the term “Dice coefficient” is a similarity measure related to the Jaccard index. The Dice coefficient is similar to Jaccard index but gives twice the weight to agreements according to the formula s=2|A ∩ B|/ (|A|+|B|).
Another aspect of the present invention concerns a method of comparing at least two L. sakei strains comprising the steps consisting of:

    
    
        a) determining the presence or absence of at least one marker gene selected from the group consisting of LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29), in a first L. sakei strain, and
        b) determining the presence or absence of said at least one marker gene in a second L. sakei strain, as defined above

wherein, if the pattern of presence or absence of said at least one marker gene is different between said first and said second L. sakei strains, then the L. sakei strains are different.

    
    


Preferably, in the above method of comparing at least two L. sakei strains, the presence or absence of all of said marker genes SEQ ID NO:1 to SEQ ID NO:29 is determined.
In particular, the above method of comparing at least two L. sakei strains further comprises the step of analysing the pattern of presence or absence of said marker genes by calculating a Jaccard index, a Dice coefficient or any binary distance matrix as defined above.
The above defined methods according to the present invention are of particular interest to identify L. sakei strains in a cocktail of strains that can be used as biopreservative on meat and/or fish. Another issue of characterizing such cocktail of strains is to quantify each strain present in said cocktail.
Accordingly, another aspect of the present invention relates to a method of quantifying a specific L. sakei strain in a sample, which comprises the steps consisting of:

    
    
        a) identifying and/or characterizing the L. sakei strains present in said sample by a method as defined above,
        b) determining the marker gene(s) which is(are) differently present in said specific L. sakei strain with respect to other strains present in said sample, and
        c) amplifying said gene(s) by quantitative PCR.

Methods of quantitative PCR are well-known in the art and include real-time PCR, competitive PCR and radioactive PCR.

    
    


As used herein the term “sample” encompasses samples in which one or several L. sakei strain(s) is/are present, optionally in combination with other bacterial species.
Accordingly, the above method of quantifying a specific L. sakei strain in a sample may include a prior step of individualizing said strain from other strains present in the sample.
The following examples further describe the way the present inventors have identified said marker genes. These examples are illustrative, without being limited, of the above defined methods.
EXAMPLES
Example 1
Detection and Analysis of Genomic Islands in Lactobacillus sakei 23K Chromosome
This example describes the identification by the inventors of genomic islands in L. sakei 23K chromosome, which carry the genes used as markers in the present invention.
It is admitted that bacterial chromosomes are divided into two sets of genes: the core gene pool, that comprises genes that persisted in the species for long enough to show maintenance among strains, and the flexible gene pool, that comprises variable and auxiliary genes generally clustered in islands and often acquired by horizontal gene transfer (HGT).
The aim of the present study is therefore to identify the putative variable genomic islands in the genome of L. sakei 23K and to characterize the cellular functions, estimated to be ecologically important, that could be ascribed to HGT.
Experimental Procedures
Codon Usage Analysis
The codon usage signature of highly expressed genes was searched by using the CAI introduced by Sharp and Weng-Hsiung (Sharp and Weng-Hsiung 1987 Nucleic Acids Res. 15:1281-1295). The codon frequencies used as a reference for CAI computation were estimated on the genes encoding ribosomal proteins and aminoacyl-tRNA synthetases found in the L. sakei 23K genome. Secondly, unsupervised clustering of the coding sequences (CDSs) was performed on the basis of their usage of the synonymous codons with an algorithm designed by the inventors. This approach is based on the use of a mixture model that accounts for the choice of synonymous codons in different gene classes. Briefly, once the number of gene classes is defined by the user, the model is fitted to the sequences using an iterative algorithm and the final clustering associates each CDS to its most probable class.
Here follows a more formal description of the mixture models and of the clustering procedure: let n denote the number of codon usage class, or component, of the mixture model. Each class u is characterized by its incidence in the mixture qu and a set probabilities (fu,i,j), 1≦i≦20, 1≦j≦si where fu,i,j corresponds to the probability of using the jth synonymous codon for amino acid i and si denotes the number of synonymous codon for amino acid i. Maximum Likelihood estimates of the model parameters qu and fu,i,j are obtained using the Expectation-Maximization algorithm. Finally, the final clustering is obtained after computing the probability of each class u for each CDS k with Bayes' formula: pk(u) is proportional to qu fu,i,jc(k,i,j) where pk(u) denotes the probability that CDS k belongs to class u and c(k,i,j) counts the occurrences of the jth synonymous codon for amino-acid i in CDS k. A detailed presentation of the mixture model framework can be found in McLachlan and Peel (McLachlan and Peel 2000 In Finite mixture models. New York: Wiley-Interscience). Finally, a correspondence analysis was performed to allow the graphical display of the cloud of CDSs in a similar way as described in Médigue et al., 1991 (Médigue et al., 1991 J. Mol. Biol. 222:851-856). Short CDSs (less than 300 bps-100 aa long) and the two cysteine codons were omitted from the correspondence analysis as they do not allow accurate estimation of the associated relative codon frequencies. In addition, codon frequency estimates were smoothed by adding a pseudo-count 1 to all counts.
All analysis were performed with an ad-hoc Perl script named CODONUSAGE.PL responsible for computing A+T-content and CAI and also for calling both a C++ program that performs the new cluster analysis described above and a R script that performs the correspondence analysis with the CA function of the MULTIV package (R Development Core Team 2006 R: A language and environment for statistical computing. Vienne, Austria: R Foundation for Statistical Computing). All programs can be downloaded at URL: http://genome.jouy.inra.fr/ ̃pnicolas/codonmixture/.
Results
L. sakei Codon Usage is Shaped by Chromosome Replication and Translation Efficiency
In order to detect possible HGT events, an array of statistical measurements was used to characterize the compositional trends of the L. sakei 23K gene pool including G+C content analysis, Codon Adaptation Index (CAI) computations, correspondence analysis of synonymous codon usage and unsupervised clustering of the gene pool with a new method based on a mixture modelling of the usage of synonymous codons as described above.
Unsupervised clustering of the complete gene pool allowed identifying four groups of genes. Genes of the first and second groups account for 74.5% of the total gene pool. They are highly preferentially located on the leading strand of chromosome replication (93%) and are distinguished by their expression level as measured with the CAI (reflecting the protein expression level): group 1 genes have low or average CAI values whereas group 2 genes show high CAI values. This last group comprises typical highly expressed genes such as S-layer-like proteins, Cpl-like chaperones and highly expressed metabolically important enzymes such as those of the glycolytic pathway. Group 3 encompasses 18.5% of the total gene pool and contains mostly genes located on the lagging strand (87%).
Therefore, three groups of genes whose codon usage are strongly shaped by their orientation relative to the chromosome replication and by their putative level of translation efficiency were able to be revealed.
More importantly, the models allowed demarcating a fourth group of genes (7% of the total gene pool) with atypical A+T-rich content and low CAI value. The compositional pattern of this fourth group may have been shaped by HGT as genes with low G+C-content have been described to be related to HGT in many bacteria (Médigue al., 1991 J. Mol. Biol. 222:851-856).
Atypical A+T-Rich L. sakei 23K Genes are Clustered in Genomic Islands.
To verify the HGT-related origin of this fourth group of genes, their corresponding genetic context in the chromosome was further looked at in more details. It was presently first noticed that most of these atypical CDSs are grouped to form large islands of functionally related genes often located nearby mobile elements, a classical feature of horizontally transferred genetic clusters (Ochman et al., 2000 Nature 405:299-304).
Furthermore, genes with atypical codon usage were often encoding products homologous to proteins only found in phylogenetic distant bacterial genera.
From these data, the inventors demonstrated that the putative HGT gene pool was thus comprised of 27 genomic islands (from 1.6 to 28 kb with an average size of 7.3 kb) and of 49 single genes, covering altogether 235 kb (12.5%) of L. sakei 23K chromosome. It has also been assumed that the 27 genomic islands may pertain to the dispensable L. sakei 23K genome.
Example 2
Intra-Species Genomic Diversity and Natural Population Structure of the Meat-Borne Lactic Acid Bacterium Lactobacillus sakei 
This example describes the classification, by the inventors, of L. sakei strains using marker genes.
In this study, the inventors have performed a combination of several typing techniques including pulse-field gel electrophoresis (PFGE) genome mapping, PCR-based detection of genetic markers taken from a pool of variable genes identified by the inventors (see example 1) for hierarchical clustering of the strains, and finally, a proteomic comparison to evaluate their respective phenotypic diversity. Isolates from diverse laboratory collections corresponding to various geographical locations and to various sources of meat- or fish-related products have been specifically chosen, based on the expectation that these undomesticated strains would represent the diversity of the natural L. sakei population.
These results provide an integrated genomic-based framework for classifying the repertoire of molecular subtypes of L. sakei isolates.
Experimental Procedures
Bacterial Strains and Culture Conditions
All L. sakei and L. curvatus strains used in this study are described in Table 2. Bacterial strains were grown to midlog exponentional phase in MRS broth medium (Difco) (De man et al., 1960 J. Appl. Bacteriol. 23:133-135) at 30° C. For proteomic studies, bacterial strains were grown at 30° C. in chemically defined medium MCD (Lauret et al., 1996 Appl Environ Microbiol. 62:1922-1927) supplemented with 0.5% glucose.
Molecular Biology Techniques
Substrative suppressive hybridization (HSS) of L. sakei 332F was performed as follows: to prepare tester strain 332F, cured from its endogenous plasmid pRV500 (Alpert et al., 2003 Appl Environ Microbiol. 69:5574-5584), the parent L. sakei 332 was electroporated by the method of Berthier et al., (Berthier et al., 1996 Appl. Envir. Microbiol. 62:3037-3041) with pRV566 plasmid carrying resistance to erythromycin and derived from pRV500 replicon (Alpert et al., 2003 Appl. Envir. Microbiol. 69:5574-5584).
One Erythromycin resistant clone was further cultivated for 200 generations in MRS broth without antibiotic at 30° C. Several dilutions from the last culture were plated on MRS agar. Replica plating of 200 clones on MRS agar with or without erythromycin (5 μg/ml) allowed selecting one erythromycin sensitive clone. The loss of pR566 plasmid was verified by Southern blotting (ECLTM direct nucleic acid labelling, Amersham Biosciences) using a probe specific of the repA gene. The corresponding strain was named 332F.
HSS experiment was carried out with Clontech PCR-selectedTM bacterial genome subtraction kit according to the manufacturer recommendations and by using L. sakei 23K as driver. This technique led to the identification of 8 new genes absent from L. sakei 23K. FGP21-0001 gene from L. sakei 21 was identified after sequencing of a PCR product (LSA0565 to LSA566) giving an unexpected size and revealing a new type of bacteriocin immunity-like protein-encoding gene.







TABLE 2







L. sakei and L. curvatus strains used in this study.












Number






of


Laboratory Collection - Country
strains
Names (synonyms)
Isolated from
Reference














PART 1: L. sakei strains






INRA - Jouy en Josas - FRANCE
12
23K, 14, 18, 21, 33, 64, 72,
Various French-style
Berthier et al., 1999


Unité Flore Lactique &

112, 134, 156
fermented dry sausages
Int. J. Syst. Bacteriol.


Environnement carné



49: 997-1007




160K
Fresh horse meat
Berthier et al., 1999






Int. J. Syst. Bacteriol.






49: 997-1007




JG3
Fresh beef meat
Berthier et al., 1999






Int. J. Syst. Bacteriol.






49: 997-1007


INRA - Theix - France
11
195, 205, 300, 332
Vacuum-packed beef meat
Champomier et al., 1987


Station de Recherches sur la



Ann Inst Pasteur Microbiol.


Viande



138: 751-758




331, 495, 504, 532, 710, 741
Vacuum-packed pork meat
Montel et al., 1991






J Appl Bacteriol 70: 469-472




L110
Starter for French-style
Champomier et al., 1987





fermented dry sausage
Ann Inst Pasteur Microbiol.






138: 751-758


IFREMER - Nantes - France
5
SF770, SF771, SF841, SF842,
Smoked salmon
Joffraud et al., 2006


Laboratoire de Génie Alimentaire

SF843

Int. J. Food Microbiol.






112: 51-61


IRTA - Monells - SPAIN
10
CTC014, CTC041, CTC163,
Various Spanish-style
Hugas et al., 1995


Meat Technology Center

CTC287, CTC335, CTC427,
fermented dry sausages
J. Appl. Bacteriol.




CTC429, CTC494, CTC6469,
(including Chorizo)
79: 332-330




CTC6626

Hugas et al., 1993






Int J Food Microbiol






18: 107-113


ARC - Langford, Bristol -
6
LV5, LV21, LV92
Vacuum-packed pork meat -
Shaw et al., 1984


UNITED KINGDOM


bacon
J Appl Bacteriol 56: 25-40


Meat Research Institute

LV52, LV59
Vacuum-packed lamb meat
Shaw et al., 1984






J Appl Bacteriol 56: 25-40




LV34
Vacuum-packed beef meat
Shaw et al., 1984






J Appl Bacteriol 56: 25-40


AgResearch - Hamilton -
4
AGR46, AGR48, AGR51,
Chilled lamb meat
This study


NEW ZEALAND

AGR53


Meat Science Group


Mahidol University - Bangkok -
2
TISTR890, TISTR911
Nham (Thai-style fermented
Tanasupawat et al., 1983


THAILAND


pork sausage)
J. Gen. Appl. Microbiol.


Faculty of Science



29: 487: 506






Noonpakdee et al., 1996






Asia Pacific J. Mol. Biol.






Biotechnol. 4: 229-235


Universität Hohenheim -
9
LTH673, LTH675, LTH677,
Various German-style
Vogel et al., 1991


GERMANY

LTH5728
fermented moist-type
FEMS Microbiol Lett


Institüt für Libbensmitteltechnologie


sausages
68: 183-190




LTH1764, LTH2070
Sauerkraut
Vogel et al., 1993






J Appl Bacteriol 74: 295-300




LTH5588, LTH5589, LTH5590
Human feces
Walter et al., 2001






Appl. Envir. Microbiol.






67: 2578-2585.


Universität Berlin - GERMANY
1
CIP105422T (CCUG31331)
Raw German-style sausage
Klein et al., 1996


Institute of Meat Hygiene &



Int. J. Sys. Bacteriol.


Technology



46: 367-376


Kulmbach - GERMANY
1
Lb706
Fresh beef meat
Schillinger et al., 1987


Federal Centre for Meat Research



Food Microbiol.






4: 199-208


MATFORSK - {acute over (Å)}s - NORWAY
6
MF1048, MF2091, MF2092
Smoked salmon
This study


Norwegian Food Research Institute

MF2088, MF2089, MF209
Rakfisk (Scandinavian
This study





fermented trout)


CONICET - Tucumán -
1
CRL1467
Argentinean-style
Fontana et al., 2005


ARGENTINA


fermented dry sausage
J Microbiol Methods


Centra de Referencia para



63: 254-263


Lactobacilos (CERELA)


American Type Culture Collection
1
ATCC15521T (DSM20017)
Spoiled moto for Saké
Katagiri et al., 1934


(Original isolate from JAPAN)


manufacture
Bull. Agr. Chem. Soc. Jpn.






10: 156-157


University of Tokyo - JAPAN
5
YMN243, YME344, YMN540,
Various fresh meat
Morishita et al., 1986


National Institute of Health

YMN557, V553
products
Int. J. Food Microbiol. 3: 19-29


Part 2: L. curvatus strains


American Type Culture Collection
1
ATCC25601T (DSM20019)
Milk
Torriani et al., 1996


(Original isolate from GERMANY)


IRTA - Monells - SPAIN
1
CTC424
Spanish-style fermented dry
Hugas et al., 1993


Meat Technology Centre


sausage
Int J Food Microbiol






18: 107-113









PCR-Based Detection of Genes

The presence or absence of the flexible gene pool identified by the inventors was investigated using conventional PCR-based detection and verification of genomic islands for strain clustering as follows:
PCR template was 100 ng of chromosomal DNA extracted from the 73 L. sakei strains and from the 2 L. curvatus strains. Experiments were conducted twice to confirm the negative results. In case of weak or spurious amplifications, PCR products were sequenced to check nucleotide polymorphism between strains. If necessary, primers were redesigned.
Extraction of chromosomal DNA from L. sakei and L. curvatus was performed by the method of Anderson & McKay (Anderson and McKay, 1983 Appl Environ Microbiol. 46:549-552). For each PCR amplification, the pair of primers was designed so that the expected length of the products was less than 2 kb.
The PCR cycling conditions were 94° C. for 4 min followed by 30 cycles of 94° C. for 1 min, 55° C. for 1 min and 72° C. for 3 min. All PCR products were examined using 1% agarose gels and stained with ethidium bromide. To confirm the truncation of some genes or the products of unexpected sizes, 10 μl of the amplicons were treated with 0.1 unit of Shrimp alkaline phosphatase (USB corporation) and 1 U of exonuclase I (E. coli) (Biolabs) in 20 mM Tris-HCl pH 8.0, 10 mM MgCl2 buffer for 1 hour at 37° C., followed by 10 min inactivation at 94° C. The products were then sequenced using standard technology (http://www.the-mwg.com).
PFGE Experiments and I-CeuI Pattern Analysis.
Pulse-field gel electrophoresis and I-CeuI digestion pattern analysis were carried out as described by Dudez et al., 2002 (Dudez et al., 2002 Microbiology. 148:421-31). An average of 4 gels was performed for each strain. The distribution of the strains according to their genome size was examined using the HIST function and the probability DENSITY function of the R statistical package (R Development Core Team 2006 R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing). A Gaussian probability distribution and a smoothing bandwidth of 30 (average standard deviation of genome size estimation) were chosen for the analyses.
Clustering of Strains.
The gene contents of the strains tested were described by using a two-character matrix (genes×isolates) with 0 for absence and 1 for presence of a gene. Genes truncated by IS elements were considered as distinct genetic identities than their wild-type counterpart. Similarities between the strains were determined using the Jaccard's correlation coefficient (Jaccard P., 1901 Bulletin de la Societe Vaudoise des Sciences Naturelles, 37:241-272). The unsupervised hierarchical clustering was performed using the complete linkage on the similarity matrix. The following functions : DIST, HCLUST and DENDROGRAM of the R statistical package (R Development Core Team 2006 R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing) were used to generate the clustering dendrogram. The R package PVCLUST (Suzuki and Shimodeira, 2006 Bioinformatics. 22:1540-1542) was used for multiscale bootstrap resampling to assess the statistical stability of each node. The number of bootstrap replicates was 1,000. Approximately unbiased P-values ≧90% and Jaccard's similarity coefficient ≧50% were used to discriminate the possible strain clusters.
2D Gel Electrophoresis and Identification of Proteins by Peptide Mass Fingerprinting.
Bacterial-extract preparation, electrophoresis were performed by standard methods (Jofré et al., 2007 Res. Microbiol. 158:512-520). Gels were analyzed by Image Master software (Amersham Pharmacia Biotech). Spots were excised from Coomassie-stained gels as described by Marceau et al., (Marceau et al., 2004 Appl Environ Microbiol. 70:7260-7268) and Mass spectrometry analyses were performed as previously described by Guillot et al., (Guillot et al., 2000 Int J Food Microbiol. 55:47-51). MS-Fit (University of California San Francisco Mass Spectrometry Facility; http://prospector.ucsf.edu) and Mascot (Matrix Science Inc., Boston, Mass.; http://www.matrixscience.com/search_form_select.html), installed on a local server, were used to identify proteins from peptide mass fingerprints. All searches were performed against the L. sakei 23K database (http://www.migale.jouy.inra.fr/sakei).
Results
Selection of L. sakei Strains.
To carefully estimate the biodiversity of the natural L. sakei population, the inventors took care to analyse strains which were isolated from a variety of meat- or fish-related food products (raw or fermented) or from other sources including Sauerkraut and Human faeces. In addition, since a possible sampling bias might exist in single laboratory bacterial collections because of the isolation procedure used or the type of food materials analyzed, the L. sakei strains were selected from 14 different laboratory collections geographically scattered within Europe, Asia and New-Zealand. A total of 73 L. sakei strains were selected and analysed (Table 2) as well as two Lactobacillus curvatus strains, a close relative to L. sakei species taken as external species reference (outgroup reference).
Identification of the Main L. sakei Molecular Subtypes by PCR-Based Detection of the Flexible Gene Pool.
In example 1 describing an in silico analysis of L. sakei 23K chromosome, the inventors identified its putative flexible gene pool comprising 27 genomic islands and 49 independent genes. They decided to check the presence or absence of this pool of gene (assumed to be variable between strains) for clustering analysis of the L. sakei isolates and by using conventional PCR.
This strategy was first tested on a preliminary PCR experiments on a set of 20 strains to demonstrate the intra-species variation of these genes. Only 5 islands revealed to be highly conserved and were therefore removed from the analysis. To avoid clustering disturbance due to the highly laterally transferable mobile elements (Insertion sequences, phage and Restriction/modification systems), these genes were also discarded from the analysis. Furthermore, the inventors showed that most genes inside each genomic island were usually displaying similar patterns of variations (the whole island is usually present or absent). Therefore, to avoid a bias from the large genomic clusters (containing more genes than the smaller ones), a selection of maximum 4 genes (those eventually showing a different pattern of variation) were taken for each cluster.
In addition, 11 chromosomally-encoded genes from other L. sakei strains that were absent from L. sakei 23K chromosome, were incorporated in the analysis. These genes were partly chosen from previously published clusters and partly taken from partial genome sequencing of other L. sakei strains. This selection resulted thus in  ̃80% of genes originated from L. sakei 23K (representing 21 genomic islands and 4 independent genes) and  ̃20% of genes from other strains. The characteristics of these 60 genes are summarized in Tables 3 and 4
Based on the PCR analysis of the 60 genes, the inventors attempted to classify the L. sakei natural isolates by using unsupervised complete-linkage hierarchical clustering algorithm and by estimating P-values via multiscale bootstrap resampling to assess the uncertainty of the clustering analysis (FIG. 1). From the resulting dendrogram, at least 11 bootstrap-supported clusters of strains divided in three main groups were clearly identified by the inventors.
Cluster A comprises the reference strain 23K and clusters B to D comprise strains closely related to this cluster forming together group 1. The other groups represent clusters of strains which are hierarchically less related to group 1 and comprise the L. sakei subsp. carnosus type strain CIP 105422T (cluster G, group 2) and the L. sakei subsp. sakei type strain ATCC 15521T (cluster K, group 3), the latter being the most distantly related to group 1. Three strains (332, 21 and LTH2070) from group 3 could not be clustered with certainty between cluster H and cluster I.
Variations of Chromosome Size and Geometry Between L. sakei Genotypic Clusters.
The inventors then investigated the extent of genome size variation between L. sakei isolates by PFGE analysis of I-CeuI-digested fragments. I-CeuI-mapping of L. sakei chromosome is resulting in seven DNA fragments (Dudez et al., 2002 Microbiology. 148:421-31) of various size and is an efficient tool in resolving overall genome size and geometry between L. sakei strains.
This analysis revealed important differences in genome size between the L. sakei strains. The mean chromosome size is 2,020±30 kb for the species, but the size is ranging from 1,814±30 kb (strain CTC427) to 2,309±79 kb (strain LTH677) representing thus about  ̃25% of genome variation ( ̃500 kb). More striking is the observation that the chromosome size of L. sakei strains was not homogeneously distributed along this range (FIG. 2A) and a Gaussian probability distribution of the genome size data was suggesting a possible division of the strains into sub-populations.
To assess whether the PGFE subpopulations could be explained by the distribution of the strains in the various genotypic clusters, the inventors analysed the genome size distribution across the 11 clusters (FIG. 2B). They demonstrated that the genome size was relatively uniform among strain clusters. They also showed that L. sakei isolates from group 1 (clusters A, B, C and D) generally harbour small genome (average size of 1,915 kb±79 kb), whereas groups 2 (clusters E, F and G) and 3 (clusters H, I, J and K) comprise most of isolates with larger genome (average size of 2,055 kb±80 kb and 2,080 kb±80 kb, respectively).
This good agreement found between genotypic clusters and PGFE subpopulations prompted the inventors to analyse the correlation between the genome size of the whole chromosome and that of each of the I-CeuI-digested fragments in the 73 L. sakei isolates. A good correlation was observed, as these two parameters were showing proportional increase from the smallest L. sakei genome to the largest one.







TABLE 3







Description of genes from the flexible gene pool of L. sakei 23K strain used for clustering analysis.











Gene name or

GenBank


Genomic Island
locus-tag
Product description
Accession N°





Island 1
LSA0088
Adenine deaminase
CR936503


Island 2
LSA0118
Hypothetical protein (Putative cell-surface collagen-binding protein)
CR936503


Island 3
LSA0157
Putative hydroxyl/aromatic amino acid symporter
CR936503


Island 5
LSA0165
Putative oxidoreductase, short chain dehydrogenase/reductase family
CR936503


Island 6
LSA0172
CscC-type cell-surface protein with Invasin/Mucin-like domain and WxL
CR936503




domain



LSA0178
MarR-type Transcriptional regulator
CR936503


Island 7
LSA0211/0212
CscC-type cell-surface protein with adhesin-like domain and WxL domain
CR936503




(authentic frameshifted gene)



LSA0216
MarR-type Transcriptional regulator
CR936503


Island 8
LSA0217
Putative transcriptional regulator with a Rhodanese-like domain, ArsR family
CR936503



LSA0218
Thioredoxin, TrxA1
CR936503



LSA0219_b
Putative cyanate transport protein
CR936503


Independent gene
LSA0306
L-aspartate-beta-decarboxylase
CR936503


Independent gene
LSA0439
Hypothetical extracellular lipase/esterase precursor
CR936503


Island 11
LSA0509
2-amino-3-ketobutyrate coenzyme A ligase (Glycine acetyltransferase)
CR936503



LSA0510/511
L-threonine dehydrogenase (authentic frameshifted gene)
CR936503


Independent gene
LSA0529
MarR-type Transcriptional regulator, putative peroxide stress regulator OhrR
CR936503


Island 12
LSA0564_a to _c
Putative bacteriocin-like peptides (LSA0564_ab) and cognate immunity
CR936503




protein (LSA0564_c)
CR936503



LSA0565 to 0566
Putative bacteriocin-like peptides
CR936503



LSA0567 to 0569_b
Putative bacteriocin-like peptides (LSA0569_ab) and cognate immunity
CR936503




proteins (LSA0567 and LSA0568)


Independent gene
LSA0572
Threonine deaminase (Threonine ammonia-lyase)
CR936503


Island 14
LSA0724 to 0725
Hypothetical proteins
CR936503



LSA0727
Hypothetical cell-surface precursor
CR936503


Island 15
LSA1006
Putative zinc-containing alcohol dehydrogenase (oxidoreductase)
CR936503


Island 16
LSA1182/1183
Putative Cytochrome P450 (authentic frameshifted gene)
CR936503


Island 17
LSA1220
Triphosphoribosyl-dephospho-coenzymeA synthase
CR936503



LSA1222
Oxaloacetate decarboxylase, alpha subunit
CR936503



LSA1227
Citrate (pro-3S)-lyase ligase (citrate lyase synthetase)
CR936503



LSA1232
Citrate: Mg(2+)(H+) symporter
CR936503


Island 18
LSA1283
CscC-type cell-surface protein with WxL domain
CR936503


Island 19
LSA1509
Hypothetical protein, sigma factor-related
CR936503



LSA1510_a to _c
Putative teichoic acid/polysaccharide export protein complex
CR936503



LSA1510_d to _f
Putative Glycosyl transferases complex
CR936503



LSA1510_g
Putative priming glycosyl transferase
CR936503



LSA1512/1513
Putative polysaccharide biosynthesis protein, chain length determination
CR936503


Island 20
LSA1572
Putative teichoic acid/polysaccharide glycosyl transferase
CR936503



LSA1579/1580
Putative teichoic acid/polysaccharide export protein complex
CR936503



LSA1581
Putative teichoic acid-binding N-acetylmuramoyl L-alalanine amidase (cell
CR936503




wall hydrolase)



LSA1584/1585
Putative teichoic acid/polysaccharide glycosyl transferase
CR936503


Island 21
LSA1640
N-acetylneuraminate lyase
CR936503



LSA1641
N-acylglucosamine-6-phosphate 2-epimerase (N-acetylmannosamine-6-
CR936503




phosphate 2-epimerase)



LSA1642
Putative Solute: Na(+) symporter
CR936503


-
LSA1720
Hypothetical protein (E. coli plasmidic gene)
CR936503


Island 22
LSA1724
MarR-type Transcriptional regulator
CR936503



LSA1730
CscC-type cell-surface protein with bacterial adhesin-like domain and WxL
CR936503




domain



LSA1731
CscC-type cell-surface protein with Haemagluttinine-like domain and WxL
CR936503




domain


Island 24
LSA1806
Hypothetical protein associated with CSC-type cluster
CR936503



LSA1809
Hypothetical extracellular protein precursor associated with CSC-type cluster
CR936503


Island 27
LSA1874
MarR-type Transcriptional regulator
CR936503
















TABLE 4







Description of genes from other strains than 23K used for clustering analysis.











Gene name or

GenBank


Strain
locus-tag*
Product description
Accession N°





Lb674
sspT
Sakacin P ABC-transporter, ATP-binding and permease protein sspT
Z48542


Lb674
sspA
Bacteriocin sakacin P precursor (Sakacin 674)
Z48542


21
FGP21-0001
Putative bacteriocin immunity protein
EU391636


KG15
dsrB
Cell-surface dextransucrase precursor (sucrose 6-glycosyltransferase)
AY697434


332F
FGP332-0001
Putative 6-phospho-beta-glucosidase
EU402602


332F
FGP332-0002
CscC-type cell-surface protein with bacterial adhesin-like domain and WxL
EU402603




domain


332F
FGP332-0003
Hypothetical cell-surface protein
EU402604


332F
FGP332-0005
Putative pyridine nucleotide-disulfide oxidoreductase
EU402605


332F
FGP332-0006
Putative ferritin-like DNA-binding protein (oxidative damage protectant) (dps-
EU402605




type)


332F
lacC
Putative tagatose-6-phosphate kinase
EU402605


332F
lacG
Putative 6-phospho-beta-galactosidase
EU402605





*nomenclature used was as follow: Flexible Gene Pool-strain nam-CDS number






Global Proteomic Variability Between L. sakei Isolates from the Different Clusters.

The inventors then used two-dimensional electrophoresis to compare the proteomes of a selection of 12 strains chosen from the various genotypic clusters A to K as a final verification of the clustering. Although an average of  ̃400 spots were commonly observed in the Pi range of 4 to 7, they noticed a remarkable variation of more than 20% in this number of spots detected between the strains.
Spots representing major differences were identified by using MALDI-TOF mass spectroscopy. Most of the differences revealed to be strain specific spots of moderately or weakly expressed proteins or to proteins likely encoded by genes not present in L. sakei 23K since they could not be identified in the protein database from this reference genome. Some variations in the 2D-protein pattern were also found to be the consequence of migration differences of some strain-specific isoforms of highly expressed proteins. This information was particularly interesting because analysis of SDS-PAGE protein profiles, the technique used in other studies to define the L. sakei subspecies, is mainly based on the detection of highly expressed proteins. These data further confirmed the hierarchical clustering shown in FIG. 1.
In conclusion, the inventors provided a first insight into the possible number of molecular subtypes within the L. sakei species. From the above results this natural population can be observed as three main groups of strains, each of them subdivided into 3 to 4 clusters. There is a substantial difference in genome size between isolates of group 1 (avg. 1,915 kb) from those of groups 2 and 3 (avg.  ̃2,075 kb for both).
The inventors evaluated the extent of intra-species genomic variation of L. sakei species and generated, for the first time, a comprehensive classification of the natural isolates.
Example 3
In this example, the inventors have performed a PCR-based detection of 29 genetic markers taken from a pool of variable genes for hierarchical clustering of the strains.
All of the L. sakei and L. curvatus strains used in this study are described in Table 2.
The gene contents of the strains tested were described by using a two-character matrix (genes×isolates) with 0 for absence and 1 for presence of a gene. Similarities between the strains were determined using the Jaccard's correlation coefficient (or Jaccard index) as defined above.
A screening similar to the one described in example 2 was performed by the inventors, and resulted in the selection of only 29 marker genes summarized in Table 1.
Based on the PCR analysis of these only 29 genes, the inventors classified the L. sakei natural isolates by using unsupervised complete-linkage hierarchical clustering algorithm and by estimating P-values via multiscale bootstrap resampling to assess the uncertainty of the clustering analysis. From the resulting dendrogram, 11 bootstrap-supported clusters of strains were clearly identified, designated A to K (FIG. 3). These clusters of strains are very similar to the ones described in Example 2, obtained using about 60 marker genes.
Accordingly, the present inventors were able to identify 29 genes that were sufficient to classify correctly the L. sakei natural isolates.

## Claims

1. A method for classifying a Lactobacillus sakei strain into a subfamily which comprises the steps:
a) determining the presence or absence of the marker genes LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183 c (SEQ ID NO:3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564 a (SEQ ID NO:15), LSA0564 b (SEQ ID NO:16), LSA0564 c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29) in said Lactobacillus sakei strain;
b) analyzing the pattern of presence or absence of said marker genes of said Lactobacillus sakei strain and calculating a Jaccard index or Dice coefficient or any binary distance matrix with regard to a set of reference Lactobacillus sakei strains; and
c) classifying said Lactobacillus sakei strain into a subfamily according to the analyzed pattern of said marker genes.

2. A method of quantifying a specific Lactobacillus sakei strain in a sample containing Lactobacillus sakei strains comprising the steps:
a) genetically characterizing the Lactobacillus sakei strains present in said sample by determining the presence or absence of the marker genes LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564—a (SEQ ID NO:15), LSA0564—b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29) in said Lactobacillus sakei strains,
b) determining the marker gene(s) which is(are) differently present or absent in said specific Lactobacillus sakei strain with respect to other Lactobacillus sakei strains present in said sample, and
c) amplifying said gene(s) determined in step b) by quantitative PCR.

3. A combination of markers enabling for characterizing a Lactobacillus sakei strain which comprises the marker genes LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29).

4. A DNA array which comprises a combination of markers according to claim 3.

5. A method of comparing at least two Lactobacillus sakei strains comprising the steps:
a) determining the presence or absence of the marker genes LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29), in a first Lactobacillus sakei strain, and
b) determining the presence or absence of said marker genes in a second Lactobacillus sakei strain,
wherein, if the pattern of presence or absence of said marker genes is different between said first and said second Lactobacillus sakei strains, then the Lactobacillus sakei strains are different.

6. The method according to claim 5, which further comprises the step of analysing the pattern of presence or absence of said marker genes by calculating a Jaccard index or Dice coefficient or any binary distance matrix.

7. The method of claim 5 wherein said steps of determining the presence of absence of marker genes are performed employing a DNA array which comprises a combination of markers enabling for characterizing a Lactobacillus sakei strain which combination comprises the marker genes LSA1641 (SEQ ID NO:1), LSA1182 (SEQ ID NO:2), LSA1183_c (SEQ ID NO: 3), LSA0172 (SEQ ID NO:4), LSA1731 (SEQ ID NO:5), LSA0211 (SEQ ID NO:6), LSA0212 (SEQ ID NO:7), LSA1579 (SEQ ID NO:8), LSA1580 (SEQ ID NO:9), LSA0118 (SEQ ID NO:10), LSA0529 (SEQ ID NO:11), LSA0439 (SEQ ID NO:12), LSA0572 (SEQ ID NO:13), LSA0219b (SEQ ID NO:14), LSA0564_a (SEQ ID NO:15), LSA0564_b (SEQ ID NO:16), LSA0564_c (SEQ ID NO:17), FGP21-0001 (SEQ ID NO:18), sspA (SEQ ID NO:19), spiA (SEQ ID NO:20), FGP332-0001 (SEQ ID NO:21), FGP332-0002 (SEQ ID NO:22), FGP332-0007 (SEQ ID NO:23), FGP332-0008 (SEQ ID NO:24), FGP332-0009 (SEQ ID NO:25), FGP332-0010 (SEQ ID NO:26), FGP332-0011 (SEQ ID NO:27), FGP332-0012 (SEQ ID NO:28), and FGP332-0013 (SEQ ID NO:29).

8. The method according to claim 1, wherein the presence or absence of said marker genes is determined by amplification.

9. The method according to claim 1, wherein the presence or absence of said marker genes is determined by hybridization with probes specific for said marker genes.

10. The method according to claim 2, wherein the presence or absence of said marker genes is determined by amplification.

11. The method according to claim 2, wherein the presence or absence of said marker genes is determined by hybridization with probes specific of said marker genes.

12. The method according to claim 2 wherein a classification of said Lactobacillus sakei strain is performed by analyzing the pattern of presence or absence of said marker genes and calculating a Jaccard index or Dice Coefficient or any binary distance matrix with regard to a set of reference Lactobacillus sakei strains.

13. The method according to claim 5, wherein the presence or absence of said marker genes is determined by amplification.

14. The method according to claim 13 which further comprises the step of analysing the pattern of presence or absence of said marker genes by calculating a Jaccard index or Dice coefficient or any binary distance matrix.

15. The method according to claim 5, wherein the presence or absence of said marker genes is determined by hybridization with probes specific for said marker genes.

16. The method according to claim 15 which further comprises the step of analysing the pattern of presence or absence of said marker genes by calculating a Jaccard index or Dice coefficient or any binary distance matrix.