Assessing the reliability of eBURST using simulated populations with known ancestry

Background  

The program eBURST uses multilocus sequence typing data to divide bacterial populations into groups of closely related strains (clonal complexes), predicts the founding genotype of each group, and displays the patterns of recent evolutionary descent of all other strains in the group from the founder. The reliability of eBURST was evaluated using populations simulated with different levels of recombination in which the ancestry of all strains was known.     

eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.     

Sequence Analysis of 96 Genomic Regions Identifies Distinct Evolutionary Lineages within CC156, the Largest Streptococcus pneumoniae Clonal Complex in the MLST Database

Multi-Locus Sequence Typing (MLST) of Streptococcus pneumoniae is based on the sequence of seven housekeeping gene fragments. The analysis of MLST allelic profiles by eBURST allows the grouping of genetically related strains into Clonal Complexes (CCs) including those genotypes with a common descent from a predicted ancestor. However, the increasing use of MLST to characterize S. pneumoniae strains has led to the identification of a large number of new Sequence Types (STs) causing the merger of formerly distinct lineages into larger CCs. An example of this is the CC156, displaying a high level of complexity and including strains with allelic profiles differing in all seven of the MLST loci, capsular type and the presence of the Pilus Islet-1 (PI-1). Detailed analysis of the CC156 indicates that the identification of new STs, such as ST4945, induced the merging of formerly distinct clonal complexes. In order to discriminate the strain diversity within CC156, a recently developed typing schema, 96-MLST, was used to analyse 66 strains representative of 41 different STs. Analysis of allelic profiles by hierarchical clustering and a minimum spanning tree identified ten genetically distinct evolutionary lineages. Similar results were obtained by phylogenetic analysis on the concatenated sequences with different methods. The identified lineages are homogenous in capsular type and PI-1 presence. ST4945 strains were unequivocally assigned to one of the lineages. In conclusion, the identification of new STs through an exhaustive analysis of pneumococcal strains from various laboratories has highlighted that potentially unrelated subgroups can be grouped into a single CC by eBURST. The analysis of additional loci, such as those included in the 96-MLST schema, will be necessary to accurately discriminate the clonal evolution of the pneumococcal population.     

Multilocus sequence type system for the plant pathogen Xylella fastidiosa and relative contributions of recombination and point mutation to clonal diversity

Multilocus sequence typing (MLST) identifies and groups bacterial strains based on DNA sequence data from (typically) seven housekeeping genes. MLST has also been employed to estimate the relative contributions of recombination and point mutation to clonal divergence. We applied MLST to the plant pathogen Xylella fastidiosa using an initial set of sequences for 10 loci (9.3 kb) of 25 strains from five different host plants, grapevine (PD strains), oleander (OLS strains), oak (OAK strains), almond (ALS strains), and peach (PP strains). An eBURST analysis identified six clonal complexes using the grouping criterion that each member must be identical to at least one other member at 7 or more of the 10 loci. These clonal complexes corresponded to previously identified phylogenetic clades; clonal complex 1 (CC1) (all PD strains plus two ALS strains) and CC2 (OLS strains) defined the X. fastidiosa subsp. fastidiosa and X. fastidiosa subsp. sandyi clades, while CC3 (ALS strains), CC4 (OAK strains), and CC5 (PP strains) were subclades of X. fastidiosa subsp. multiplex. CC6 (ALS strains) identified an X. fastidiosa subsp. multiplex-like group characterized by a high frequency of intersubspecific recombination. Compared to the recombination rate in other bacterial species, the recombination rate in X. fastidiosa is relatively low. Recombination between different alleles was estimated to give rise to 76% of the nucleotide changes and 31% of the allelic changes observed. The housekeeping loci holC, nuoL, leuA, gltT, cysG, petC, and lacF were chosen to form the basis of a public database for typing X. fastidiosa (www.mlst.net). These loci identified the same six clonal complexes using the strain grouping criterion of identity at five or more loci with at least one other member.     

Displaying the relatedness among isolates of bacterial species -- the eBURST approach

Determining the most appropriate way to represent the relationships between bacterial isolates is complicated by the differing rates of recombination within species. In many cases, a bifurcating tree can be positively misleading. The recently described program eBURST can be used with multilocus data to define groups or clonal complexes of related isolates derived from a common ancestor, the patterns of descent linking them together, and the ancestral genotype. eBURST has recently been extensively updated to include additional tools for exploring the relationships between isolates. We discuss the advantages of this approach and describe its use to explore patterns of descent within clonal complexes identified using multilocus sequence typing.     

Population Structure of Clinical Vibrio parahaemolyticus from 17 Coastal Countries,Determined through Multilocus Sequence Analysis

Vibrio parahaemolyticus is a leading cause of food-borne gastroenteritis worldwide. Although this bacterium has been the subject of much research, the population structure of clinical strains from worldwide collections remains largely undescribed, and the recorded outbreaks of V. parahaemolyticus gastroenteritis highlight the need for the subtyping of this species. We present a broad phylogenetic analysis of 490 clinical V. parahaemolyticus isolates from 17 coastal countries through multilocus sequence analysis (MLST). The 490 tested isolates fell into 161 sequence types (STs). The eBURST algorithm revealed that the 161 clinically relevant STs belonged to 8 clonal complexes, 11 doublets, and 94 singletons, showing a high level of genetic diversity. CC3 was found to be a global epidemic clone of V. parahaemolyticus, and ST-3 was the only ST with an international distribution. recA was observed to be evolving more rapidly, exhibiting the highest degree of nucleotide diversity (0.028) and the largest number of polymorphic nucleotide sites (177). We also found that the high variability of recA was an important cause of differences between the results of the eBURST and ME tree analyses, suggesting that recA has a much greater influence on the apparent evolutionary classification of V. parahaemolyticus based on the current MLST scheme. In conclusion, it is evident that a high degree of genetic diversity within the V. parahaemolyticus population and multiple sequence types are contributing to the burden of disease around the world. MLST, with a fully extractable database, is a powerful system for analysis of the clonal relationships of strains at a global scale. With the addition of more strains, the pubMLST database will provide more detailed and accurate information, which will be conducive to our future research on the population structure of V. parahaemolyticus.     

The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis.

Multilocus sequence typing (MLST) is a recently developed nucleotide sequence-based method for the definitive assignment of isolates within bacterial populations to specific clones. MLST uses the same principles as multilocus enzyme electrophoresis and provides data that can be used to investigate aspects of the population genetics and evolution of bacterial species. We used an MLST data set consisting of the sequences of approximately 450-bp fragments from seven housekeeping loci from a large strain collection of Neisseria meningitidis to estimate the relative impact of recombination compared with point mutation in the diversification of N. meningitidis clonal complexes. 126 meningococcal isolates were assigned to 10 clonal complexes, 9 of which contained minor clonal variants. The allelic variation within each complex was classified as a recombinational exchange or a putative point mutation through a comparison of the sequences of each variant allele with that of the allele typically found in the clonal complex. The nine clonal complexes contained a total of 23 allelic variants, and analysis of the sequences of these variant alleles revealed that a single nucleotide site in a meningococcal housekeeping gene is at least 80-fold more likely to change as a result of recombination than as a result of mutation. This value is estimated to be 10-50-fold for Escherichia coli and approximately 50-fold for Streptococcus pneumoniae.     

Inferring a population structure for Staphylococcus epidermidis from multilocus sequence typing data

Despite its importance as a human pathogen, information on population structure and global epidemiology of Staphylococcus epidermidis is scarce and the relative importance of the mechanisms contributing to clonal diversification is unknown. In this study, we addressed these issues by analyzing a representative collection of S. epidermidis isolates from diverse geographic and clinical origins using multilocus sequence typing (MLST). Additionally, we characterized the mobile element (SCCmec) carrying the genetic determinant of methicillin resistance. The 217 S. epidermidis isolates from our collection were split by MLST into 74 types, suggesting a high level of genetic diversity. Analysis of MLST data using the eBURST algorithm revealed the existence of nine epidemic clonal lineages that were disseminated worldwide. One single clonal lineage (clonal complex 2) comprised 74% of the isolates, whereas the remaining isolates were clustered into 8 minor clonal lineages and 13 singletons. According to our evolutionary model, SCCmec was acquired at least 56 times by S. epidermidis. Although geographic dissemination of S. epidermidis strains and the value of the index of association between the alleles, 0.2898 (P < 0.05), support the clonality of S. epidermidis species, examination of the sequence changes at MLST loci during clonal diversification showed that recombination gives rise to new alleles approximately twice as frequently as point mutations. We suggest that S. epidermidis has a population with an epidemic structure, in which nine clones have emerged upon a recombining background and evolved quickly through frequent transfer of genetic mobile elements, including SCCmec.     

Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing

Vibrio parahaemolyticus is an important human pathogen whose transmission is associated with the consumption of contaminated seafood. There is a growing public health concern due to the emergence of a pandemic strain causing severe outbreaks worldwide. Many questions remain unanswered regarding the evolution and population structure of V. parahaemolyticus. In this work, we describe a multilocus sequence typing (MLST) scheme for V. parahaemolyticus based on the internal fragment sequences of seven housekeeping genes. This MLST scheme was applied to 100 V. parahaemolyticus strains isolated from geographically diverse clinical (n = 37) and environmental (n = 63) sources. The sequences obtained from this work were deposited and are available in a public database (http://pubmlst.org/vparahaemolyticus). Sixty-two unique sequence types were identified, and most (50) were represented by a single isolate, suggesting a high level of genetic diversity. Three major clonal complexes were identified by eBURST analysis. Separate clonal complexes were observed for V. parahaemolyticus isolates originating from the Pacific and Gulf coasts of the United States, while a third clonal complex consisted of strains belonging to the pandemic clonal complex with worldwide distribution. The data reported in this study indicate that V. parahaemolyticus is genetically diverse with a semiclonal population structure and an epidemic structure similar to that of Vibrio cholerae. Genetic diversity in V. parahaemolyticus appears to be driven primarily by frequent recombination rather than mutation, with recombination ratios estimated at 2.5:1 and 8.8:1 by allele and site, respectively. Application of this MLST scheme to more V. parahaemolyticus strains and by different laboratories will facilitate production of a global picture of the epidemiology and evolution of this pathogen.     

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.     

The Bacillus cereus group: novel aspects of population structure and genome dynamics

AIMS: To provide new insights into the population and genomic structure of the Bacillus cereus group of bacteria. METHODS AND RESULTS: The genetic relatedness among B. cereus group strains was assessed by multilocus sequence typing (MLST) using an optimized scheme based on seven chromosomal housekeeping genes. A set of 48 strains from different clinical sources was included, and six clonal complexes containing several genetically similar isolates from unrelated patients were identified. Interestingly, several clonal groups contained strains that were isolated from similar human sources. Furthermore, comparative whole genome sequence analysis of 16 strains led to the discovery of novel ubiquitous genome features of the B. cereus group, such as atypical group II introns, IStrons, and hitherto uncharacterized repeated elements. CONCLUSIONS: The B. cereus group constitutes a coherent population unified by the presence of ubiquitous and specific genetic elements which do not show any pattern, either in their sequences or genomic locations, which allows to differentiate between the member species of the group. Nevertheless, the population is very dynamic, as particular lineages of clinical origin can evolve to form clonal complexes. At the genome level, the dynamic behaviour is indicated by the presence of numerous mobile and repeated elements. SIGNIFICANCE AND IMPACT OF THE STUDY: The B. cereus group of bacteria comprises species that are of medical and economic importance. The MLST data, along with the primers and protocols used, will be available in a public, web-accessible database (http://mlstoslo.uio.no).     

Multilocus sequence typing as an approach for population analysis of Medicago-nodulating rhizobia

Multilocus sequence typing (MLST), a sequence-based method to characterize bacterial genomes, was used to examine the genetic structure in a large collection of Medicago-nodulating rhizobial strains. This is the first study where MLST has been applied in conjunction with eBURST analysis to determine the population genetic structure of nonpathogenic bacteria recovered from the soil environment. Sequence variation was determined in 10 chromosomal loci of 231 strains that predominantly originated from southwest Asia. Genetic diversity for each locus ranged from 0.351 to 0.819, and the strains examined were allocated to 91 different allelic profiles or sequence types (STs). The genus Medicago is nodulated by at least two groups of rhizobia with divergent chromosomes that have been classified as Sinorhizobium meliloti and Sinorhizobium medicae. Evidence was obtained that the degree of genetic exchange among the chromosomes across these groups is limited. The symbiosis with Medicago polymorpha of nine strains placed in one of these groups, previously identified as S. medicae, ranged from ineffective to fully effective, indicating that there was no strong relationship between symbiotic phenotype and chromosomal genotype.     

Population structure of Pseudomonas aeruginosa from five Mediterranean countries: evidence for frequent recombination and epidemic occurrence of CC235

Several studies in recent years have provided evidence that Pseudomonas aeruginosa has a non-clonal population structure punctuated by highly successful epidemic clones or clonal complexes. The role of recombination in the diversification of P. aeruginosa clones has been suggested, but not yet demonstrated using multi-locus sequence typing (MLST). Isolates of P. aeruginosa from five Mediterranean countries (n = 141) were subjected to pulsed-field gel electrophoresis (PFGE), serotyping and PCR targeting the virulence genes exoS and exoU. The occurrence of multi-resistance (≥3 antipseudomonal drugs) was analyzed with disk diffusion according to EUCAST. MLST was performed on a subset of strains (n = 110) most of them had a distinct PFGE variant. MLST data were analyzed with Bionumerics 6.0, using minimal spanning tree (MST) as well as eBURST. Measurement of clonality was assessed by the standardized index of association (I_A ^S). Evidence of recombination was estimated by ClonalFrame as well as SplitsTree4.0. The MST analysis connected 70 sequence types, among which ST235 was by far the most common. ST235 was very frequently associated with the O11 serotype, and frequently displayed multi-resistance and the virulence genotype exoS ⁻/exoU ⁺. ClonalFrame linked several groups previously identified by eBURST and MST, and provided insight to the evolutionary events occurring in the population; the recombination/mutation ratio was found to be 8.4. A Neighbor-Net analysis based on the concatenated sequences revealed a complex network, providing evidence of frequent recombination. The index of association when all the strains were considered indicated a freely recombining population. P. aeruginosa isolates from the Mediterranean countries display an epidemic population structure, particularly dominated by ST235-O11, which has earlier also been coupled to the spread of ß-lactamases in many countries.     

Multilocus Sequence Typing of Leuconostoc gelidum subsp. gasicomitatum,a Psychrotrophic Lactic Acid Bacterium Causing Spoilage of Packaged Perishable Foods

Leuconostoc gelidum subsp. gasicomitatum is a psychrotrophic lactic acid bacterium (LAB) that causes spoilage of a variety of modified-atmosphere-packaged (MAP) cold-stored food products. During the past 10 years, this spoilage organism has been increasingly reported in MAP meat and vegetable products in northern Europe. In the present study, the population structure within 252 L. gelidum subsp. gasicomitatum strains was determined based on a novel multilocus sequence-typing (MLST) scheme employing seven housekeeping genes. These strains had been isolated from meat and vegetable sources over a time span of 15 years, and all 68 previously detected pulsed-field gel electrophoresis (PFGE) genotypes were represented. A total of 46 sequence types (STs) were identified, with a majority of the strains (>60%) belonging to three major STs, which were grouped into three clonal complexes (CCs) and 17 singletons by Global Optimal eBURST (goeBURST). The results by Bayesian analysis of population structure (BAPS) mostly correlated with the grouping by goeBURST. Admixture analysis by BAPS indicated a very low level of exchange of genetic material between the subpopulations. Niche specificity was observed within the subpopulations: CC1 and BAPS cluster 1 consisted mostly of strains from a variety of MAP meats, whereas vegetable strains grouped together with strains from MAP poultry within CC2 and BAPS cluster 2. The MLST scheme presented in this study provides a shareable and continuously growing sequence database enabling global comparison of strains associated with spoilage cases. This will further advance our understanding of the microbial ecology of this industrially important LAB.     

Genetic diversity of Campylobacter jejuni isolates from farm animals and the farm environment

The genetic diversity of Campylobacter jejuni isolates from farm animals and their environment was investigated by multilocus sequence typing (MLST). A total of 30 genotypes, defined by allelic profiles (assigned to sequence types [STs]), were found in 112 C. jejuni isolates originating in poultry, cattle, sheep, starlings, and slurry. All but two of these genotypes belonged to one of nine C. jejuni clonal complexes previously identified in isolates from human disease and retail food samples and one clonal complex previously associated with an environmental source. There was some evidence for the association of certain clonal complexes with particular farm animals: isolates belonging to the ST-45 complex predominated among poultry isolates but were absent among sheep isolates, while isolates belonging to the ST-61 and ST-42 complexes were predominant among sheep isolates but were absent from the poultry isolates. In contrast, ST-21 complex isolates were distributed among the different isolation sources. Comparison with MLST data from 91 human disease isolates showed small but significant genetic differentiation between the farm and human isolates; however, representatives of six clonal complexes were found in both samples. These data demonstrate that MLST and the clonal complex model can be used to identify and compare the genotypes of C. jejuni isolates from farm animals and the environment with those from retail food and human disease.     

Multiple-locus sequence typing analysis of Bacillus cereus and Bacillus thuringiensis reveals separate clustering and a distinct population structure of psychrotrophic strains

We used multilocus sequence typing (MLST) to characterize phylogenetic relationships for a collection of Bacillus cereus group strains isolated from forest soil in the Paris area during a mild winter. This collection contains multiple strains isolated from the same soil sample and strains isolated from samples from different sites. We characterized 115 strains of this collection and 19 other strains based on the sequences of the clpC, dinB, gdpD, panC, purF, and yhfL loci. The number of alleles ranged from 36 to 53, and a total of 93 allelic profiles or sequence types were distinguished. We identified three major strain clusters-C, T, and W-based on the comparison of individual gene sequences or concatenated sequences. Some less representative clusters and subclusters were also distinguished. Analysis of the MLST data using the concept of clonal complexes led to the identification of two, five, and three such groups in clusters C, T, and W, respectively. Some of the forest isolates were closely related to independently isolated psychrotrophic strains. Systematic testing of the strains of this collection showed that almost all the strains that were able to grow at a low temperature (6 degrees C) belonged to cluster W. Most of these strains, including three independently isolated strains, belong to two clonal complexes and are therefore very closely related genetically. These clonal complexes represent strains corresponding to the previously identified species Bacillus weihenstephanensis. Most of the other strains of our collection, including some from the W cluster, are not psychrotrophic. B. weihenstephanensis (cluster W) strains appear to comprise an effectively sexual population, whereas Bacillus thuringiensis (cluster T) and B. cereus (cluster C) have clonal population structures.     

Genetic Diversity of Campylobacter jejuni Isolates from Farm Animals and the Farm Environment

The genetic diversity of Campylobacter jejuni isolates from farm animals and their environment was investigated by multilocus sequence typing (MLST). A total of 30 genotypes, defined by allelic profiles (assigned to sequence types [STs]), were found in 112 C. jejuni isolates originating in poultry, cattle, sheep, starlings, and slurry. All but two of these genotypes belonged to one of nine C. jejuni clonal complexes previously identified in isolates from human disease and retail food samples and one clonal complex previously associated with an environmental source. There was some evidence for the association of certain clonal complexes with particular farm animals: isolates belonging to the ST-45 complex predominated among poultry isolates but were absent among sheep isolates, while isolates belonging to the ST-61 and ST-42 complexes were predominant among sheep isolates but were absent from the poultry isolates. In contrast, ST-21 complex isolates were distributed among the different isolation sources. Comparison with MLST data from 91 human disease isolates showed small but significant genetic differentiation between the farm and human isolates; however, representatives of six clonal complexes were found in both samples. These data demonstrate that MLST and the clonal complex model can be used to identify and compare the genotypes of C. jejuni isolates from farm animals and the environment with those from retail food and human disease.     

Genetic diversity of clinical isolates of <Emphasis Type="Italic">Bacillus cereus</Emphasis>using multilocus sequence typing

Background  

Bacillus cereus is most commonly associated with foodborne illness (diarrheal and emetic) but is also an opportunistic pathogen that can cause severe and fatal infections. Several multilocus sequence typing (MLST) schemes have recently been developed to genotype B. cereus and analysis has suggested a clonal or weakly clonal population structure for B. cereus and its close relatives B. anthracis and B. thuringiensis. In this study we used MLST to determine if B. cereus isolates associated with illnesses of varying severity (e.g., severe, systemic vs. gastrointestinal (GI) illness) were clonal or formed clonal complexes.     

Multiple-Locus Sequence Typing Analysis of Bacillus cereus and Bacillus thuringiensis Reveals Separate Clustering and a Distinct Population Structure of Psychrotrophic Strains

We used multilocus sequence typing (MLST) to characterize phylogenetic relationships for a collection of Bacillus cereus group strains isolated from forest soil in the Paris area during a mild winter. This collection contains multiple strains isolated from the same soil sample and strains isolated from samples from different sites. We characterized 115 strains of this collection and 19 other strains based on the sequences of the clpC, dinB, gdpD, panC, purF, and yhfL loci. The number of alleles ranged from 36 to 53, and a total of 93 allelic profiles or sequence types were distinguished. We identified three major strain clusters—C, T, and W—based on the comparison of individual gene sequences or concatenated sequences. Some less representative clusters and subclusters were also distinguished. Analysis of the MLST data using the concept of clonal complexes led to the identification of two, five, and three such groups in clusters C, T, and W, respectively. Some of the forest isolates were closely related to independently isolated psychrotrophic strains. Systematic testing of the strains of this collection showed that almost all the strains that were able to grow at a low temperature (6°C) belonged to cluster W. Most of these strains, including three independently isolated strains, belong to two clonal complexes and are therefore very closely related genetically. These clonal complexes represent strains corresponding to the previously identified species Bacillus weihenstephanensis. Most of the other strains of our collection, including some from the W cluster, are not psychrotrophic. B. weihenstephanensis (cluster W) strains appear to comprise an effectively sexual population, whereas Bacillus thuringiensis (cluster T) and B. cereus (cluster C) have clonal population structures.     

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.