RecA gene sequence and Multilocus Sequence Typing for species-level resolution of Burkholderia cepacia complex isolates

Aim:  To identify, by means of recA sequencing and multilocus sequence typing (MLST), Burkholderia cepacia complex (BCC) isolates of environmental and clinical origin, which failed to be identified by recA RFLP and species-specific PCR.
Methods and Results:  By using recA sequence-based identification, 17 out of 26 BCC isolates were resolved at the level of species and lineage (ten Burkholderia cenocepacia IIIB, two Burkholderia arboris and five Burkholderia lata ). By using MLST method, 24 BCC isolates were identified. MLST confirmed recA sequence results, and, furthermore, enabled to identify isolates of the BCC5 group, and showed relatedness with Burkholderia contaminans for one of the two isolates not identified.
Conclusions:  recA sequence-based identification allowed to resolve, at the level of species and lineage, 65·4%, of the BCC isolates examined, whilst MLST increased this percentage to 88·5%.
Significance and Impact of the Study:  BCC isolates previously not resolved by recA RFLP and species-specific PCR were successfully identified by means of recA sequencing and MLST, which represent the most appropriate methods to identify difficult strains for epidemiological purposes and cystic fibrosis patients management.     

Detection of cultured and uncultured Burkholderia cepacia complex bacteria naturally occurring in the maize rhizosphere

The species composition of a Burkholderia cepacia complex population naturally occurring in the maize rhizosphere was investigated by using both culture-dependent and culture-independent methods. B. cepacia complex isolates were recovered from maize root slurry on the two selective media Pseudomonas cepacia azelaic acid tryptamine (PCAT) and trypan blue tetracycline (TB-T) and subjected to identification by a combination of restriction fragment length polymorphism (RFLP) analysis and species-specific polymerase chain reaction (PCR) tests of the recA gene. DNA extracted directly from root slurry was examined by means of nested PCR to amplify recA gene with species-specific B. cepacia complex primers and to obtain a library of PCR amplified recA genes. Using the culture-dependent method the species Burkholderia cepacia, Burkholderia cenocepacia, Burkholderia ambifaria and Burkholderia pyrrocinia were identified, whereas using the culture-independent method also the species Burkholderia vietnamiensis was detected. The latter method also allowed us to highlight a higher diversity within the B. cenocepacia species. In fact, by using the culture-independent method the species B. cenocepacia recA lineages IIIA and IIID besides B. cenocepacia recA lineage IIIB were detected. Moreover, higher heterogeneity of recA RFLP patterns was observed among clones assigned to the species B. cenocepacia than among B. cenocepacia isolates from selective media.     

Identification and onion pathogenicity of Burkholderia cepacia complex isolates from the onion rhizosphere and onion field soil

Burkholderia cepacia complex strains are genetically related but phenotypically diverse organisms that are important opportunistic pathogens in patients with cystic fibrosis (CF,) as well as pathogens of onion and banana, colonizers of the rhizospheres of many plant species, and common inhabitants of bulk soil. Genotypic identification and pathogenicity characterization were performed on B. cepacia complex isolates from the rhizosphere of onion and organic soils in Michigan. A total of 3,798 putative B. cepacia complex isolates were recovered on Pseudomonas cepacia azelaic acid tryptamine and trypan blue tetracycline semiselective media during the 2004 growing season from six commercial onion fields located in two counties in Michigan. Putative B. cepacia complex isolates were identified by hybridization to a 16S rRNA gene probe, followed by duplex PCR using primers targeted to the 16S rRNA gene and recA sequences and restriction fragment length polymorphism analysis of the recA sequence. A total of 1,290 isolates, 980 rhizosphere and 310 soil isolates, were assigned to the species B. cepacia (160), B. cenocepacia (480), B. ambifaria (623), and B. pyrrocinia (27). The majority of isolates identified as B. cepacia (85%), B. cenocepacia (90%), and B. ambifaria (76%) were pathogenic in a detached onion bulb scale assay and caused symptoms of water soaking, maceration, and/or necrosis. A phylogenetic analysis of recA sequences from representative B. cepacia complex type and panel strains, along with isolates collected in this study, revealed that the B. cenocepacia isolates associated with onion grouped within the III-B lineage and that some strains were closely related to strain AU1054, which was isolated from a CF patient. This study revealed that multiple B. cepacia complex species colonize the onion rhizosphere and have the potential to cause sour skin rot disease of onion. In addition, the onion rhizosphere is a natural habitat and a potential environmental source of B. cenocepacia.     

Multilocus sequence typing reveals high genetic diversity and epidemic population structure for the fish pathogen Yersinia ruckeri

Yersinia ruckeri is the causative agent of enteric redmouth in fish and one of the major bacterial pathogens causing losses in salmonid aquaculture. Previously typing methods, including restriction enzyme analysis, pulsed-field gel electrophoresis and multilocus enzyme electrophoresis (MLEE) have indicated a clonal population structure. In this work, we describe a multilocus sequence typing (MLST) scheme for Y.ruckeri based on the internal fragment sequence of six housekeeping genes. This MLST scheme was applied to 103 Y.ruckeri strains from diverse geographic areas and hosts as well as environmental sources. Sequences obtained from this work were deposited and are available in a public database (http://publmst.org/yruckeri/). Thirty different sequence types (ST) were identified, 21 of which were represented by a single isolate, evidencing high genetic diversity. ST2 comprised more than one-third of the isolates and was most frequently observed among isolates from trout. Two major clonal complexes (CC) were identified by eBURST analysis showing a common evolutionary origin for 94 isolates forming 21 STs into CC1 and for 6 isolates of 6 STs in the CC2. It was also possible to associate some unique ST with isolates from recent outbreaks in vaccinated salmonid fish.     

Diversity and distribution of Burkholderia cepacia complex in the rhizosphere of rice and maize

A survey of Burkholderia cepacia complex (Bcc) species was conducted in agricultural fields within Hangzhou, China. Out of the 251 bacterial isolates recovered on the selective media from the rhizosphere of rice and maize, 112 of them were assigned to Bcc by PCR assays. The species composition of the Bcc isolates was analyzed by a combination of recA-restriction fragment length polymorphism assays, species-specific PCR tests and recA gene sequencing. The results revealed that the majority belong to B. cepacia, Burkholderia cenocepacia recA lineage IIIB, Burkholderia vietnamiensis and Burkholderia pyrrocinia. Burkholderia cenocepacia and B. vietnamiensis dominated the rhizosphere of maize and rice, respectively, indicating that species composition and abundance of Bcc may vary dramatically in different crop rhizospheres. In addition, one isolate (R456) formed a single discrete cluster within the phylogenetic analysis of the Bcc recA gene, and it may belong to a new genomovar.     

Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones

Bartonella henselae is a zoonotic pathogen and the causative agent of cat scratch disease and a variety of other disease manifestations in humans. Previous investigations have suggested that a limited subset of B. henselae isolates may be associated with human disease. In the present study, 182 human and feline B. henselae isolates from Europe, North America and Australia were analysed by multi-locus sequence typing (MLST) to detect any associations between sequence type (ST), host species and geographical distribution of the isolates. A total of 14 sequence types were detected, but over 66% (16/24) of the isolates recovered from human disease corresponded to a single genotype, ST1, and this type was detected in all three continents. In contrast, 27.2% (43/158) of the feline isolates corresponded to ST7, but this ST was not recovered from humans and was restricted to Europe. The difference in host association of STs 1 (human) and 7 (feline) was statistically significant (P< or =0.001). eBURST analysis assigned the 14 STs to three clonal lineages, which contained two or more STs, and a singleton comprising ST7. These groups were broadly consistent with a neighbour-joining tree, although splits decomposition analysis was indicative of a history of recombination. These data indicate that B. henselae lineages differ in their virulence properties for humans and contribute to a better understanding of the population structure of B. henselae.     

Genotypes and toxin gene profiles of Staphylococcus aureus clinical isolates from China

A total of 108 S. aureus isolates from 16 major hospitals located in 14 different provinces in China were characterized for the profiles of 18 staphylococcal enterotoxin (SE) genes, 3 exfoliatin genes (eta, etb and etd), and the toxic shock syndrome toxin gene (tsst) by PCR. The genomic diversity of each isolate was also evaluated by pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and accessory gene regulator (agr) typing. Of these strains, 90.7% (98/108) harbored toxin genes, in which tsst was the most prevalent toxin gene (48.1%), followed by sea (44.4%), sek (42.6%) and seq (40.7%). The see and etb genes were not found in any of the isolates tested. Because of high-frequency transfer of toxin gene-containing mobile genetic elements between S. aureus strains, a total of 47 different toxin gene combinations were detected, including a complete egc cluster in 19 isolates, co-occurrence of sea, sek and seq in 38 strains, and sec and sel together in 11 strains. Genetic typing by PFGE grouped all the strains into 25 clusters based on 80% similarity. MLST revealed 25 sequence types (ST) which were assigned into 16 clonal complexes (CCs) including 2 new singletons. Among these, 11 new and 6 known STs were first reported in the S. aureus strains from China. Overall, the genotyping results showed high genetic diversity of the strains regardless of their geographical distributions, and no strong correlation between genetic background and toxin genotypes of the strains. For genotyping S. aureus, PFGE appears to be more discriminatory than MLST. However, toxin gene typing combined with PFGE or MLST could increase the discriminatory power of genotyping S. aureus strains.     

Multi-locus sequence typing (MLST) of non-fermentative Gram-negative bacilli isolated from bloodstream infections in southern Poland

Non-fermentative Gram-negative bacilli are now one of the most important causes of severe infections in Polish hospitals. Acinetobacter species are serious concern because of the high prevalence of multi-drug resistance among strains. Resistance profiles for 53 Gram-negative non-fermentative blood isolates were done. MLST was carried out using 44 strains representing the most commonly isolated species: A. baumannii, P. aeruginosa, and S. maltophilia. MLST revealed that all 22 A. baumannii belonged to sequence type (ST) 2. The P. aeruginosa isolates belonged to 10 different STs. Four S. maltophilia isolates matched STs present in the database (ST4, ST15, ST116, ST142), seven isolates showing novel sequence types. Among P. aeruginosa and S. maltophilia PFGE confirmed the genetical variety of strains.     

Differential interaction of bacterial species from the Burkholderia cepacia complex with human airway epithelial cells

To increase knowledge of the pathogenic potential of the Burkholderia cepacia complex (BCC), we investigated the effects of reference strains of the nine BCC species on human bronchial epithelial cells in vitro. B. multivorans exhibited the highest rates of adherence to and internalization by host cells. Two out of three clinical isolates recovered from cystic fibrosis patients confirmed the B. multivorans high adhesiveness. All four B. multivorans isolates exhibited an aggregated pattern of adherence but any of them expressed cable pili. When bacteria were centrifuged onto cell cultures to circumvent their poor adhesiveness, B. pyrrocinia exhibited the highest internalization rate, followed by B. multivorans. The percentages of apoptotic cells in cultures infected with B. cepacia, B. multivorans, B. cenocepacia (subgroups IIIA and IIIB), B. stabilis and B. vietnamiensis were significantly higher than in control non-infected cultures. All nine BCC species triggered a similar release of the inflammatory cytokine IL-8, that was not reduced by cell treatment with cytochalasin D. Hence, our data demonstrate, for the first time, that all BCC species exhibit a similar ability to induce the expression of host immune mediators whereas they differ on their ability to adhere to, invade and kill airway epithelial cells.     

Molecular characterization of Enterococcus faecium isolated from hospitalized patients in Korea

AIMS: Multilocus sequence typing (MLST) was performed for vancomycin-resistant Enterococcus faecium (VREF) from diverse geographical areas in Korea to obtain insights into the genetic relationships with other molecular profiles. To understand the diversity of lineages, vancomycin-susceptible E. faecium (VSEF) were included. METHODS AND RESULTS: A total of 60 E. faecium isolates were analysed by MLST and esp profile. Molecular typing of Tn1546 of 30 VREF strains was evaluated by overlapping PCR of Tn1546 and DNA sequencing. Seven sequence types (ST) were found among 30 VSEF isolates, and four STs were found among 30 VREF isolates. The types most frequently encountered were ST 78 (26 isolates) and ST 203 (16 isolates). Of the 60 E. faecium isolates, 35 isolates were positive for the esp gene. On molecular typing of Tn1546, all VREF isolates were divided into four main types. Strains with the same ST showed divergence in Tn1546 types and strains with the same Tn1546 type represented different STs. CONCLUSIONS: An association between Tn1546 typing and MLST was not found. SIGNIFICANCE AND IMPACT OF THE STUDY: These results suggest that the horizontal spread of Tn1546 between strains plays a major role in the dissemination of vancomycin resistance in Korea.     

Clonal Analysis of Meningococci during a 26 Year Period Prior to the Introduction of Meningococcal Serogroup C Vaccines

Meningococcal disease remains a public health burden in the UK and elsewhere. Invasive Neisseria meningitidis, isolated in Scotland between 1972 and 1998, were characterised retrospectively to examine the serogroup and clonal structure of the circulating population. 2607 isolates causing invasive disease were available for serogroup and MLST analysis whilst 2517 were available for multilocus sequence typing (MLST) analysis only. Serogroup distribution changed from year to year but serogroups B and C were dominant throughout. Serogroup B was dominant throughout the 1970s and early 1980s until serogroup C became dominant during the mid-1980s. The increase in serogroup C was not associated with one particular sequence type (ST) but was associated with a number of STs, including ST-8, ST-11, ST-206 and ST-334. This is in contrast to the increase in serogroup C disease seen in the 1990s that was due to expansion of the ST-11 clonal complex. While there was considerable diversity among the isolates (309 different STs among the 2607 isolates), a large proportion of isolates (59.9%) were associated with only 10 STs. These data highlight meningococcal diversity over time and the need for ongoing surveillance during the introduction of new meningococcal vaccines.     

Population structure and evolution of pathogenicity of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is an enteric human pathogen but is widespread in the environment. Pathogenicity is determined by a number of virulence factors, including the virulence plasmid pYV, the high-pathogenicity island (HPI), and the Y. pseudotuberculosis-derived mitogen (YPM), a superantigen. The presence of the 3 virulence factors varies among Y. pseudotuberculosis isolates. We developed a multilocus sequence typing (MLST) scheme to address the population structure of Y. pseudotuberculosis and the evolution of its pathogenicity. The seven housekeeping genes selected for MLST were mdh, recA, sucA, fumC, aroC, pgi, and gyrB. An MLST analysis of 83 isolates of Y. pseudotuberculosis, representing 19 different serotypes and six different genetic groups, identified 61 sequence types (STs) and 12 clonal complexes. Out of 26 allelic changes that occurred in the 12 clonal complexes, 13 were mutational events while 13 were recombinational events, indicating that recombination and mutation contributed equally to the diversification of the clonal complexes. The isolates were separated into 2 distinctive clusters, A and B. Cluster A is the major cluster, with 53 STs (including Y. pestis strains), and is distributed worldwide, while cluster B is restricted to the Far East. The YPM gene is widely distributed on the phylogenetic tree, with ypmA in cluster A and ypmB in cluster B. pYV is present in cluster A only but is sporadically absent in some cluster A isolates. In contrast, an HPI is present only in a limited number of lineages and must be gained by lateral transfer. Three STs carry all 3 virulence factors and can be regarded as high-pathogenicity clones. Isolates from the same ST may not carry all 3 virulence factors, indicating frequent gain or loss of these factors. The differences in pathogenicity among Y. pseudotuberculosis strains are likely due to the variable presence and instability of the virulence factors.     

DNA identification and characterization of Campylobacter jejuni and Campylobacter coli isolated from caecal samples of chickens in Grenada

Aims: To speciate Campylobacter strains from the caeca of chickens in Grenada using PCR and to evaluate DNA‐based typing methods for the characterization of these isolates. Methods and Results: Isolates were speciated with two multiplex PCR assays and were typed with flaA‐RFLP, pulsed‐field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). Results confirmed that Campylobacter coli strains were more predominant than Campylobacter jejuni strains. From 56 isolates, 18 were misidentified using biochemical tests. PFGE typing gave the highest discriminatory power among the methods used (Simpson’s index of diversity, D = 0·9061). However, the combination of flaA‐RFLP, PFGE and MLST results gave the highest discrimination for subtyping of these isolates (D = 0·9857). A band position tolerance of 4% in Bio Numerics was the most appropriate for the analysis of this database. MLST profiles were generally concordant with PFGE and/or flaA‐RFLP types. Several isolates exhibited new MLST sequence types (STs), and 43 of the 49 Camp. coli strains belonged to the ST‐828 clonal complex. Conclusions: Campylobacter coli was the most prevalent species isolated from broilers and layers in Grenada, and a combination of restriction and sequence methods was most appropriate for the typing of Camp. coli isolates. Campylobacter coli STs clustered with described poultry‐associated Camp. coli STs by phylogenetic analysis. Significance and Impact of the Study: Further studies to understand the predominance of Camp. coli within Campylobacter spp. from chickens in Grenada may help elucidate the epidemiology of these pathogens in chickens.     

Multilocus sequence typing of Campylobacter jejuni, and the correlation between clonal complex and pulsed-field gel electrophoresis macrorestriction profile

The characterization of Campylobacter jejuni has been significantly improved by the use of multilocus sequence typing (MLST), which allows the relationship between isolates to be determined. The sequence types (STs) of 261 isolates of C. jejuni from New Zealand were determined. Isolates were obtained from a range of sources including chicken meat, cattle, pigs, duck, sheep, water and human infections. Thirty-two new alleles and 44 new STs were identified. Comparison of the MLST data and pulsed-field gel electrophoresis macrorestriction profiles showed that the macrorestriction profiles were good predictors of the clonal complex (CC) but not ST. All the major CCs identified elsewhere in the world were found in New Zealand as well as the association of certain CCs with particular animal niches. The majority of new STs identified were from river water isolates.     

Multilocus sequence typing scheme for bacteria of the Bacillus cereus group

In this study we developed a multilocus sequence typing (MLST) scheme for bacteria of the Bacillus cereus group. This group, which includes the species B. cereus, B. thuringiensis, B. weihenstephanensis, and B. anthracis, is known to be genetically very diverse. It is also very important because it comprises pathogenic organisms as well as bacteria with industrial applications. The MLST system was established by using 77 strains having various origins, including humans, animals, food, and soil. A total of 67 of these strains had been analyzed previously by multilocus enzyme electrophoresis, and they were selected to represent the genetic diversity of this group of bacteria. Primers were designed for conserved regions of housekeeping genes, and 330- to 504-bp internal fragments of seven such genes, adk, ccpA, ftsA, glpT, pyrE, recF, and sucC, were sequenced for all strains. The number of alleles at individual loci ranged from 25 to 40, and a total of 53 allelic profiles or sequence types (STs) were distinguished. Analysis of the sequence data showed that the population structure of the B. cereus group is weakly clonal. In particular, all five B. anthracis isolates analyzed had the same ST. The MLST scheme which we developed has a high level of resolution and should be an excellent tool for studying the population structure and epidemiology of the B. cereus group.     

MLST clustering of Campylobacter jejuni isolates from patients with gastroenteritis,reactive arthritis and Guillain–Barré syndrome

Aims: To determine the diversity and population structure of Campylobacter jejuni (C. jejuni) isolates from Danish patients and to examine the association between multilocus sequence typing types and different clinical symptoms including gastroenteritis (GI), Guillain–Barré syndrome (GBS) and reactive arthritis (RA). Methods and Results: Multilocus sequence typing (MLST) was used to characterize 122 isolates, including 18 from patients with RA and 8 from patients with GBS. The GI and RA isolates were collected in Denmark during 2002–2003 and the GBS isolates were obtained from other countries. In overall, 51 sequence types (STs) were identified within 18 clonal complexes (CCs). Of these three CCs, ST‐21, ST‐45 and ST‐22 clonal complexes accounted for 64 percent of all isolates. The GBS isolates in this study significantly grouped into the ST‐22 clonal complex, consistent with the PubMLST database isolates. There was no significant clustering of the RA isolates. Conclusions: Isolates from Denmark were found to be highly genetically diverse. GBS isolates grouped significantly with clonal complex ST‐22, but the absence of clustering of RA isolates indicated that the phylogenetic background for this sequela could not be reconstructed using variation in MLST loci. Possibly, putative RA‐associated genes may vary, by recombination or expression differences, independent of MLST loci. Significance and Impact of the Study: MLST typing of C. jejuni isolates from Danish patients with gastroenteritis confirmed that the diversity of clones in Denmark is comparable to that in other European countries. Furthermore, a verification of the grouping of GBS isolates compared to RA isolates provides information about evolution of the bacterial population resulting in this important sequela.     

Diversification and Distribution of Ruminant Chlamydia abortus Clones Assessed by MLST and MLVA

Chlamydia abortus, an obligate intracellular bacterium, is the most common infectious cause of abortion in small ruminants worldwide and has zoonotic potential. We applied multilocus sequence typing (MLST) together with multiple-locus variable-number tandem repeat analysis (MLVA) to genotype 94 ruminant C. abortus strains, field isolates and samples collected from 1950 to 2011 in diverse geographic locations, with the aim of delineating C. abortus lineages and clones. MLST revealed the previously identified sequence types (STs) ST19, ST25, ST29 and ST30, plus ST86, a recently-assigned type on the Chlamydiales MLST website and ST87, a novel type harbouring the hemN_21 allele, whereas MLVA recognized seven types (MT1 to MT7). Minimum-spanning-tree analysis suggested that all STs but one (ST30) belonged to a single clonal complex, possibly reflecting the short evolutionary timescale over which the predicted ancestor (ST19) has diversified into three single-locus variants (ST86, ST87 and ST29) and further, through ST86 diversification, into one double-locus variant (ST25). ST descendants have probably arisen through a point mutation evolution mode. Interestingly, MLVA showed that in the ST19 population there was a greater genetic diversity than in other STs, most of which exhibited the same MT over time and geographical distribution. However, the evolutionary pathways of C. abortus STs seem to be diverse across geographic distances with individual STs restricted to particular geographic locations. The ST30 singleton clone displaying geographic specificity and represented by the Greek strains LLG and POS was effectively distinguished from the clonal complex lineage, supporting the notion that possibly two separate host adaptations and hence independent bottlenecks of C. abortus have occurred through time. The combination of MLST and MLVA assays provides an additional level of C. abortus discrimination and may prove useful for the investigation and surveillance of emergent C. abortus clonal populations.     

Multilocus sequence-typing for characterization of Moscow strains of Haemophilus influenzae type b

Haemophilius influenzae, type b (Hib) bacteria, were genotyped by multilocus sequence typing (MLST) using 5 loci (adk, fucK, mdh, pgi, recA). 42 Moscow Hib strains (including 38 isolates form cerebrospinal fluid of children, who had purulent meningitis in 1999-2001, and 4 strains isolated from healthy carriers of Hib), as well as 2 strains from Yekaterinburg were studied. In MLST a strain is characterized, by alleles and their combinations (an allele profile) referred to also as sequence-type (ST). 9 Sts were identified within the Russian Hib bacteria: ST-1 was found in 25 strains (57%), ST-12 was found in 8 strains (18%), ST-11 was found in 4 strains (9%) and ST-15 was found in 2 strains (4.5%); all other STs strains (13, 14, 16, 17, 51) were found in isolated cases (2.3%). A comparison of allelic profiles and of nucleotide sequences showed that 93% of Russian isolates, i.e. strain with ST-1, 11, 12, 13, 15 and 17, belong to one and the same clonal complex. 2 isolates from Norway and Sweden from among 7 foreign Hib strains studied up to now can be described as belonging to the same clonal complex; 5 Hib strains were different from the Russian ones.     

Genetic diversity of Staphylococcus equorum isolates from Saeu-jeotgal evaluated by multilocus sequence typing

Staphylococcus equorum, the predominant bacterial species detected in Saeu-jeotgal, a Korean high-salt fermented seafood, is a candidate starter bacterium for Saeu-jeotgal fermentation. A multilocus sequence typing (MLST) scheme was developed to evaluate the genetic diversity and background of S. equorum strains isolated from Saeu-jeotgal. A total of 135 strains, including 117 isolates from Saeu-jeotgal, and others from Myeolchi-jeotgal, sausage, cheese and horse skin, were subjected to MLST, and the internal fragments of seven housekeeping genes, aroE, dnaJ, glpF, gmk, hsp60, mutS, and pta, were compared. This MLST scheme produced 45 sequence types (STs) and the eBURST algorithm clustered the STs into nine clonal groups and seven singletons. Clonal group 1, the major group, consisted of 30 isolates from cheese, Saeu-jeotgal and sausages, which were classified into 12 STs. The predominant ST, ST26, comprised 25 isolates and presented as a singleton. Most of the isolates from Myeolchi-jeotgal and sausages clustered on two different branches of a phylogenetic tree generated with a cluster analysis using the maximum likelihood algorithm. This MLST scheme established the genetic backgrounds of S. equorum strains isolated from different types of food. Among the housekeeping genes used for MLST, gmk had the fewest allele types and fairly low sequence identities (74.0–90.0 %) within the Staphylococcus species. Therefore, sequence analyses of the gmk gene and 16S rRNA gene can be used for the accurate and rapid identification of S. equorum.