Skip to main content

Characterization of Salmonella Typhimurium isolates from domestically acquired infections in Finland by phage typing, antimicrobial susceptibility testing, PFGE and MLVA



Salmonella enterica spp. enterica serotype Typhimurium (STM) is the most common agent of domestically acquired salmonellosis in Finland. Subtyping methods which allow the characterization of STM are essential for effective laboratory-based STM surveillance and for recognition of outbreaks. This study describes the diversity of Finnish STM isolates using phage typing, antimicrobial susceptible testing, pulsed-field gel electrophoresis (PFGE) and multilocus variable-number tandem repeat analysis (MLVA), and compares the discriminatory power and the concordance of these methods.


A total of 375 sporadic STM isolates were analysed. The isolates were divided into 31 definite phage (DT) types, dominated by DT1 (47 % of the isolates), U277 (9 % of the isolates) and DT104 (8 % of the isolates). Of all the isolates, 62 % were susceptible to all the 12 antimicrobials tested and 11 % were multidrug resistant. Subtyping resulted in 83 different XbaI-PFGE profiles and 111 MLVA types. The three most common XbaI-PFGE profiles (STYM1, STYM7 and STYM8) and one MLVA profile with three single locus variants accounted for 56 % and 49 % of the STM isolates, respectively. The studied isolates showed a genetic similarity of more than 70 % by XbaI-PFGE. In MLVA, 71 % of the isolates lacked STTR6 and 77 % missed STTR10p loci. Nevertheless, the calculated Simpson’s diversity index for XbaI-PFGE was 0.829 (95 % CI 0.792−0.865) and for MLVA 0.867 (95 % CI 0.835−0.898). However, the discriminatory power of the 5-loci MLVA varied among the phage types. The highest concordance of the results was found between XbaI-PFGE and phage typing (adjusted Wallace coefficient was 0.833 and adjusted Rand coefficient was 0.627).


In general, the calculated discriminatory power was higher for genotyping methods (MLVA and XbaI-PFGE) than for phenotyping methods (phage typing). Overall, comparable diversity indices were calculated for PFGE and MLVA (both DI > 0.8). However, MLVA was phage type dependent providing better discrimination of the most common phage types. Furthermore, 5-loci MLVA was a less laborious method and easier to interpret than XbaI-PFGE. Thus, the laboratory-based surveillance of the Finnish human STM infections has been conducted with a combination of phage typing, antimicrobial susceptibility testing and 5-loci MLVA since January 2014.


By a rough estimation, 5.1 million foodborne non-typhoidal salmonella infections and consequently 8400 deaths occur in Europe annually [1]. Salmonella enterica ssp. enterica serotype Typhimurium (STM) is, after the serotype Enteritidis, the second most commonly isolated serotype from human salmonelloses in Europe [2]. In Finland, the total number of reported salmonella infections ranged from approximately 2200 to 3100 cases per year in 2007–2012 [3]. Of these infections, 10–30 % were considered domestically acquired and serotype Typhimurium dominated among the domestic isolates [3]. Most salmonelloses are characterized by mild-to-moderate self-limited diarrhoea but also serious disease resulting in death has been reported especially in connection with STM infection in children less than 5 years old [4, 5]. Salmonella infections, including those caused by STM, are considered sporadic and their sources remain unknown. However, national outbreaks [6, 7] and large outbreaks across the national borders associated with STM have also been reported [8, 9]. The STM outbreaks are frequently associated with consumption of contaminated raw vegetables, fruits, poorly cooked meat or meat products and eggs [1013].

Molecular subtyping of bacterial isolates is crucial when monitoring the circulation of bacterial strains across the geographic regions. Subtyping methods must fulfil many requirements such as high discriminatory power, easy performance, clear interpretation and possibility for international standardization. Currently, pulsed-field gel electrophoresis (PFGE), which has been internationally standardized, is considered as the golden standard for genotyping of Salmonella and it is the only generic molecular method suitable for all Salmonella serotypes [1416]. However, not all the isolates unrelated to an outbreak are different by PFGE from the outbreak-related isolates and PFGE has been shown to suffer from a high variation between the different laboratories. Thus, new typing methods have been developed around the world. Among these, multilocus variable-number tandem repeat analysis (MLVA) has been used as a typing tool for tracing back the sources of foodborne disease outbreaks, in bacterial surveillance and in research [17]. MLVA can be used to discriminate between genetically closely related strains [18]. A database of tandem repeats for several completed Salmonella genome sequences is publicly available at [19]. These variable-number tandem repeat (VNTR) markers have been applied for typing and clustering of several Salmonella serotypes, for example, Dublin [20], Enteritidis [21], Gallinarum [22], Montevideo [23], Newport [24], Paratyphi A [25], Typhi [26] and Typhimurium [27]. MLVA typing has also been used successfully in several STM outbreak investigations [2830]. Reduced typing time, easy handling and clearer interpretation of MLVA results, compared to those of PFGE, are beneficial in an outbreak situation. The major disadvantage of MLVA is serotype-specificity. There is, for example, no MLVA protocol which is suitable for all the over 2600 Salmonella serotypes and, currently, internationally standardized MLVA protocols are available for only a few serotypes. In Europe, the 5-loci MLVA scheme developed by Lindstedt et al. 2004, has been widely adopted for STM by several national public health laboratories [17]. The MLVA loci nomenclature devised by Larson et al., 2009 allows the normalisation of the raw data and enables a direct comparison of the MLVA results between laboratories [31].

The aim of this study was to characterize Finnish STM isolates using several pheno- and genotyping methods. For that, 375 sporadic STM isolates from Finns who had not been abroad recently were characterized by phage typing, antimicrobial susceptibility testing, XbaI-PFGE and 5-loci MLVA during the 5-year period between November 2007 and October 2012. The specific aim was to determinate (i) the genetic diversity among the Finnish STM isolates, especially those of endemic definitive phage type DT1, and to detect clonal clusters, (ii) to compare the discriminatory power of XbaI-PFGE and 5-loci MLVA and (iii) to investigate whether the 5-loci MLVA results are in concordance with the phage typing or XbaI-PFGE results. In addition, the potential and possible added value of 5-loci MLVA as a typing tool in a laboratory for the STM surveillance was evaluated.


Overall diversity of the isolates

The 375 STM isolates from Finns, who had not travelled abroad recently were typed by phage typing, antimicrobial susceptibility testing, PFGE and MLVA (5 primer pairs in Table 1). Among the isolates, 31 distinct phage types were detected (Table 2). Of the isolates, 65 % (245 isolates) were assigned to the three most common DTs: DT1 (48 % of the isolates), U277 (9 % of the isolates) and DT104 (8 % of the isolates). Each of the remaining phage types contained less than 4 % of the isolates (Table 2). The isolates associated with the phage type called “reacts but does not conform” (RDNC) accounted for 10 % (38 isolates). Based on similar phage reactions, these isolates were divided into five groups RDNC1-5 (22 isolates). In addition, 16 isolates showing each a unique RDNC pattern were found. Of all the isolates, 62 % (231 isolates) were susceptible to all antimicrobial agents tested. Within the most common phage types, 74 % (133 of 180isolates) of DT1, 87 % (29 of 33) of U277 but only 3 % (1 of 32) of DT104 isolates was susceptible to all antimicrobials tested. On the whole, 11 % (43 isolates) were multidrug resistant (MDR) (Table 2). The three most common MDR patterns were ACSSuTNx (14 isolates), ACSSuT (10 isolates) and ASSuT (7 isolates), representing 8 % of the isolates. Of the DT104 isolates with any resistance marker, 45 % (14 isolates) displayed resistance to six antimicrobials (ACSSuTNx). All of the 14 Nx resistant isolates had a decreased sensitivity to ciprofloxacin (MIC values were ranging from 0.125 to 0.19 mg/L). In addition to the 14 DT104isolates, a decreased sensitivity to ciprofloxacin was observed in two DT120 (0.25–0.38 mg/L), two U311 (0.19–0.38 mg/L) and one DT104B (3 mg/L) isolates. Furthermore, one isolate of phage type U311 showed combined resistance to more than six antimicrobials including an intermediate resistance to cefotaxime. None of the isolates were resistant to imipenem.

Table 1 Characteristics of each MLVA loci for S. Typhimurium
Table 2 Distribution of phage types, multidrug resistance, PFGE profiles and MLVA types among the 375 STM isolates from the domestic infections

In the molecular analyses, 83 different XbaI-PFGE and 111 different MLVA profiles were generated (Table 2). Among DT1 isolates, the most common XbaI-PFGE profile STYM1 contained 10 MLVA profiles, of which the MLVA profile 3-16-NA-NA-0311 was the most prevalent. The U277 isolates were divided into six XbaI-PFGE and seven MLVA profiles (Table 2). The most common XbaI-PFGE profile STYM8 showed five different MLVA profiles, of which again the 3-16-NA-NA-0311 was the most common. The DT104 isolates were divided into six XbaI-PFGE and 17 MLVA profiles (Table 2). Here, the most common XbaI-PFGE profile STYM7 contained 13 different MLVA profiles, of which the profile 3-13-5-12-0311 was the most common.

Although higher number of MLVA types (n = 111) than XbaI-PFGE profiles (n = 83) were generated (p ≤ 0.05), the calculated diversity indices for MLVA (Simpson’s DI = 0.867 with 95 % CI 0.835−0.898; Shannon’s DI, H’ = 4.697) and for XbaI-PFGE (Simpson’s DI = 0.829 with 95 % CI 0.792−0.865; Shannon’s DI, H’ = 4.207) were comparable and the confidence intervals partly overlapped (Table 3). The discriminatory power of MLVA was better for the three most common phage types compared to XbaI-PFGE. In general, the MLVA data obtained in this study suggest that MLVA better performance compared to XbaI-PFGE varies with phage type (Table 4). In general, the MLVA data that we have obtained suggest that MLVA is phage type dependent. The discriminatory power of MLVA was better for the five most common phage types compared to XbaI-PFGE (Table 4). In addition, the Simpson’s diversity indices for each of the five MLVA loci varied being highest for locus STTR5 (DI = 0.789; 95 % CI 0.780−0.811) and lowest for locus STTR9 (DI = 0.328; 95 % CI 0.298-0.358) (Table 5). The highest variability in repeat numbers at one locus was found in locus STTR6 (21 different repeats) and the lowest in locus STTR9 (7 different repeats) (Table 5).

Table 3 Discriminatory power of three typing methods for 375 sporadic S. Typhimurium isolates
Table 4 Simpson’s diversity indices of Xbal-PFGE and 5-loci MLVA among the five most common phage types
Table 5 Simpson’s diversity index of each MLVA loci

Relatedness of the sporadic isolates by PFGE and MLVA

UPGMA clustering was done by construction of dendrograms for both PFGE and MLVA profiles (Figs. 1 and 2). The cluster analysis of XbaI-PFGE results demonstrated a genetic similarity of more than 70 % among the studied isolates (Fig 1). In 11 cases the PFGE banding patterns differed from each other only by the thickness of one band. The most common XbaI-PFGE profile STYM1 (n = 147) was evenly distributed during the study period. Nonetheless, 50 isolates showed a unique XbaI-PFGE profile. The cluster analysis of MLVA types revealed that no particular clustering emerged between the STM isolates. The majority of the isolates had MLVA type 3-16-NA-NA-0311 but also 89 unique MLVA types were detected (Fig. 2). The most common feature among the Finnish STM isolates was the lack of loci STTR6 and STTR10p, both loci were absent in 70 % of the isolates (260 isolates).

Fig. 1
figure 1

A UPGMA dendrogram of XbaI-PFGE pattern of 375 STM isolates. Each XbaI-PFGE profile (n = 83) is presented only once and the number of isolates column provides information on frequency of each XbaI-PFGE profile. STYM1 (n = 147) was the most common profile and 50 isolates had a unique profile

Fig. 2
figure 2

A UPGMA dendrogram of the MLVA types of 375 STM isolates. Each MLVA type (n = 111) is presented only once and the number of isolates column provides information on frequency of each MLVA type. MLVA type 3-16-NA-NA-0311 (n = 129) was the most common and 89 isolates had a unique MLVA type

Concordance of MLVA with phage typing and PFGE

To assess the concordance between the three subtyping methods, adjusted Wallace (AW) and adjusted Rand (AR) coefficients were calculated (Table 6). A poor directional correlation between 5-loci MLVA and XbaI-PFGE and phage typing was found. The probability of two isolates that have the same MLVA type to share the same XbaI-PFGE profile was only 30 % meaning that neither phage type nor PFGE profile could be predicted directly from the MLVA type. According to the AW and AR coefficients, the highest concordance was between XbaI-PFGE and phage typing (AW = 0.833; AR = 0.627).

Table 6 Concordance between the typing methods presented by adjusted Wallace (AW) and adjusted Rand (AR) coefficient


STM is the most frequently isolated Salmonella serotype from the Finns who suffer from salmonellosis and have not travelled abroad prior to their symptoms [3]. Thus, accurate surveillance of STYM is of public health importance. Since the 1960s, STM has caused the majority of the domestic salmonellosis in Finland (personal communication, Anja Siitonen, THL). Salmonella isolates, however, are rare in Finnish production animals (cattle, poultry and swine) and only about 1 % of the animals tested are salmonella positive [32]. In addition, the phage types circulating among the animals are well known [32, 33]. For example, STM DT1 which was the most common cause of domestic human salmonella infections in this study, is endemic in Finland and isolates with this phage type have been isolated from production animals (poultry and cattle) and wild animals since the 1980s [3, 33, 34]. Phage types DT40 and U277 have been isolated from domestic turkey, wild birds, hedgehogs and regularly from human infections as well [32].

More than 60 % of all the clinical isolates and almost 75 % of all the DT1 isolates were susceptible to all 12 antimicrobials tested in this study. Similar susceptibility results were observed during previous years as well [35]. This could result from moderate usage of antibiotics in Finnish animal husbandry [36]. Nonetheless, 11 % of our isolates were found to be MDR including nine different phage types (DT104, DT104B and DT193 being the most common). Several studies have shown MDR to be more common among certain STM phage types e.g. in DT104 [37, 38]. Because MDR Salmonella is rare among domestic production animals, the occurrence of MDR isolates among Finnish patients could indicate that at least some of the domestically acquired infections might actually be due to imported foods sold in supermarkets and restaurants [33]. According to a study in Finland, the risk of receiving Salmonella from imported beef or beef products was close to 1 % [39]. On the other hand, in a recent source attribution study in Sweden it was estimated that almost 6.5 % of the domestic salmonella infections are due to imported foods [40]. Alternatively, some of the infections could be secondary infections from patients, who have contracted their salmonella infection abroad. Furthermore, about 5 % of the domestic isolates had decreased sensitivity to ciprofloxacin, defined as MIC ≥ 0.125mg/L, which has in previous studies been associated with travelling abroad [41].

Discriminatory power of a method, defined as the ability to distinguish between unrelated isolates, is commonly obtained by calculation of Simpson’s diversity index. In order to get a more objective measure of the subtyping methods used, two different indices, Simpson’s and Shannon’s, were calculated. Both diversity indices provided similar degree of diversity and Shannon’s diversity index (J) also indicated that the different PTs, XbaI-PFGE profiles and MLVA types were equally distributed in the STM population during the study period. In addition, comparable discriminatory power was evidenced for both methods (both DI > 0.8). Sporadic domestic STM isolates showed only limited amount of diversity. The studied isolates showed a genetic similarity of more than 70 % by XbaI-PFGE. In MLVA, absence of two loci was characteristic for domestic STM isolates as 71 % of the isolates lacked STTR6 and 77 % were missing STTR10p loci. Despite of missing loci, 111 distinct MLVA types were observed. According to other studies, MLVA was highly discriminatory (DI > 0.9) for STM [30, 42]. Low diversity observed, especially among the domestic DT1, might result from factors such as geographical isolation, endemic nature of the isolates and lacking of two loci (STTR6 and STTR10pl) in Finnish STM, rather than from the low discriminatory ability of MLVA method as such. This is supported by the fact that the PFGE typing showed low diversity as well. The DT101 isolates in New Zealand [30] and the STM isolates in Malaysia [43] have also been shown to lack STTR6 and STTR10p loci.

Based on the results of this study, The National Salmonella Reference Laboratory at THL has since January 2014 used the 5-loci MLVA as the primary genotyping method for STYM isolates and XbaI-PFGE as the secondary method for epidemiological investigations. Despite the comparable degree of discriminatory power for XbaI-PFGE and MLVA, the latter was more cost-effective, faster to perform and above all, easier to interpret and apply in international comparison. However, in order to interpret 5-loci MLVA results appropriately in outbreak situation, more information on VNTR mutation rate and their stability is needed. Some studies have been conducted to determine the VNTR stability. For example, the influence of multiple freeze-thaw cycles on single-colony isolates and the effects of external stress to bacterial genotype have been studied. A study performed with Escherichia coli O157:H7 revealed that certain stress factors (temperature, UV-light and starvation) can affect the mutation rate of multiple tandem repeats in one locus more [44]. Another stability study revealed that changes were observed mainly in MLVA loci STTR6, STTR10 and STTR5, but a single change was also detected in STTR9 [38]. In conclusion, the majority of the stability studies agree that the order of instability among STM loci is following STTR6 > STTR10 > STTR5 > STTR9 > STTR3. Nevertheless, as with any typing method, the results obtained with MLVA must be viewed in conjunction with the epidemiological data. Importantly, international criteria for cluster and outbreak investigations should be agreed on.


In Finland, domestically acquired infections caused by salmonella serotype Typhimurium account for approximately 30 % of all domestic salmonella infections. In this study, the three most common phage types, XbaI-PFGE and 5-loci MLVA profiles comprised 65, 56 and 49 % of the 375 sporadic STM isolates, respectively. More than 60 % of the isolates were susceptible to all 12 antimicrobials tested. For the most common phage types, the 5-loci MLVA had better discriminatory power than XbaI-PFGE. Furthermore, in a daily use, MLVA is simple, fast, robust and cost-effective technique for subtyping STM isolates. The MLVA results are easy to interpret and reproduce, and importantly, the protocol is internationally harmonised allowing the comparison of the results between laboratories. In Finland, MLVA is used as routine subtyping method for surveillance and cluster detection of domestic STM isolates. Nonetheless, phenotyping methods (phage typing and antimicrobial susceptibility testing) are useful for outbreak detection, especially when analyzing isolates with a common MLVA type.


Bacterial isolates

Salmonella Typhimurium isolates were received from the Finnish routine clinical microbiology laboratories that are obliged to send all their Salmonella isolates to THL based on the Finnish Communicable Disease Act. During the 5-year period (November 2007 to October 2012), all of the 455 STM isolates that had antigen structure 4,5,12:i:1,2 and that were associated with domestically acquired salmonelloses were subtyped at the Bacterial Infections Unit, National Institute for Health and Welfare (THL), Finland. In this study, a Salmonella isolate was called domestic if it was isolated from a person who had no travel history in one week prior to his or her symptoms. In case of multiple STM findings in the person, the first finding was chosen for this study. Eighteen percent (80/455) of the original isolates were associated with a known STM infection epidemic or a cluster or were family-related (based on the same family name) and were therefore excluded from the study, leaving 375 sporadic isolates for the analyses.

Phenotypic characterization

All the isolates were serotyped by slide agglutination according to the White-Kauffmann-Le Minor scheme [45]. The definite phage typing of STM isolates was performed using a standard set of 38 typing phages received from the Public health England [46, 47]. The isolates showing a pattern that did not conform to any recognized phage patterns were designated as “reacts but did not conform” (RDNC). Furthermore, the isolates belonging to RDNC were additionally grouped based on their phage reactions if more than one isolate with similar phage reactions (difference of 1−3 reactions) were detected (RDNC1-5). Isolates that did not react with any of the typing phages were designated as “not typeable” (NT). The antimicrobial susceptibility to 12 antimicrobials was determined by the agar diffusion method on Müller-Hinton II agar for ampicillin (A) (10 μg), chloramphenicol (C) (30 μg), streptomycin (S) (10 μg), sulphonamide (Su) (300 μg), tetracycline (T) (30 μg), ciprofloxacin (Cp) (5 μg), trimethoprim (Tm) (5 μg), gentamicin (G) (10 μg), nalidixic acid (Nx) (30 μg), cefotaxime (Ct) (5 μg), mecillinam (M) (10 μg) and imipenem (I) (10 μg). During 2007–2010, the protocols and clinical breakpoints of the Clinical and Laboratory Standards Institute (CLSI) [48] and during 2011–2012 those of the European Committee on Antimicrobial Susceptibility Testing (EUCAST) ( were applied. Minimal inhibitory concentration (MIC) for ciprofloxacin (from 0.002 to 32 mg/L) was detected by E-Test (Biomérieux, Marcy l’Étoile, France) for the isolates that were resistant (R) or intermediate resistant (I) to Nx. MIC breakpoint ≤1 mg/L was interpreted as susceptible [48]. Multidrug resistance (MDR) was defined as resistance to four or more antimicrobials.

Genotypic characterization of the isolates


For PFGE typing, the isolates were cultivated overnight on trypticase soy agar (TSA) plates, and PFGE was performed according to the internationally standardized PulseNet protocol [49] using 15 U XbaI restriction enzyme (Roche, Basel, Switzerland) for each plug. XbaI-digested S. Braenderup (H9812) served as a DNA size marker. Banding patterns were analysed using BioNumerics v.6.6 (AppliedMaths, Kortrijk, Belgium). The bands within a size range of 33 kb and 1135 kb were included in the analysis, and isolates differing even in one banding position or in thickness of the band were assigned as different PFGE types. An unweighted pair group method with arithmetic mean (UPGMA) dendrogram was constructed using a Dice coefficient.


The MLVA was performed as described in [27, 31, 50] with some modifications (Table 1). Isolates were grown overnight on nutrient agar plates, and 1-2 colonies were suspended into 500 μl sterile water and boiled for 5 min. After a quick centrifugation, 1 μl of the supernatant was used as template in each PCR reaction. For ABI3730xl (G5 filter) suitable fluorescence-labelled forward primers (STTR3-F-NED, STTR5-F-NED, STTR6-F-FAM, STTR9-F-FAM and STTR10F-PET) were ordered from Applied Biosystems (CA, USA) and the unlabelled reverse primers from Oligomer (Espoo, Finland). A 5-plex PCR reaction was performed with a Qiagen multiplex kit (Hilden, Germany) in a total volume of 25 μl and included 2.50 pmol of primers STTR3-F, STTR3-R, STTR6-F, and STTR6-R and 1.25 pmol of primers STTR5-F, STTR5-R, STTR9-F, STTR9-R, STTR10pl-F, and STTR10pl-R. Amplification was performed with a Dyad Peltier thermal cycler (Bio-Rad, Hercules, USA), starting with 15 min at 95 °C, followed by 30 cycles of 30 s at 94 °C, 90 s at 60 °C, and 90 s at 72 °C and ending with an extension step for 10 min at 72 °C. The PCR products were diluted in a ratio of 1:85 in sterile ddH2O. Total of 10 μl of formamide/size standard mixture containing 2 μl of internal size standard GeneScan™ 600 LIZ ® (Applied Biosystem, Forter City, CA, USA) and 1 ml Hi-Di™ formamide (Applied Biosystems, Foster City, CA, USA) was applied to 96-well plates. Then, 1 μl of diluted PCR product (1:85) was added into each well. The PCR products were separated with an ABI3730xl automated DNA analyzer (Applied Biosystems, CA, USA). The size and dye label associated with each amplicon were determined using Peak Scanner v.1.0 software (Applied Biosystems, Foster City, CA, USA). Initially, a set of 33 calibration STM isolates and a correction table provided by Dr. Eva Møller Nielsen, Statens Serum Institute, Denmark, were used to calibrate the method and to normalise the raw data for the correct determination of the number of repeat units in each locus. Based on the fragment length, repeat numbers were assigned for each isolate by using arbitrary numbers [31]. Unique allelic combinations were assigned as a separate MLVA type, and all MLVA types were reported as repeat number values in the following order: STTR9-STTR5-STTR6-STTR10-STTR3. A null amplification product was considered a distinct allele after confirmation by a repeated single PCR assay and markedas -2 in the Bionumerics. MLVA results were stored at Bionumerics v.6.6 (Applied Maths, Kortrijk, Belgium). A UPGMA dendrogram was constructed using the categorical values coefficient.


Discriminatory power of phage typing, XbaI-PFGE and 5-loci MLVA were assessed using Simpson’s index of diversity (DI) and the Shannon’s diversity index (H’), which are indicator of species richness (i.e. number of subtypes), and equitability (J), which is a measure of the evenness of subtype distribution, via the online tool ( Additionally, the diversity index was assigned for each MLVA loci ( Simpson’s DIs ranges from zero (no diversity) to one (high diversity). Confidence intervals were calculated for Simpson’s DIs as described earlier [51]. The χ 2 test was used to determine the significance of the difference between typing methods, and a p-value of <0.05 was considered indicative of a significant difference. The concordance of different typing methods was determined using adjusted Rand (AR) and adjusted Wallace (AW) coefficients ( AR is a coefficient suitable for quantitative evaluation of the concordance between different microbial typing methods. AR values range from zero to one (global congruence of typing method is high). Wallace’s coefficient indicates the probability at which two isolates of the same type are also classified as the same by another method.


  1. Majowicz SE, Musto J, Scallan E, Angulo FJ, Kirk M, O’Brien SJ, et al. The global burden of nontyphoidal Salmonella gastroenteritis. Clin Infect Dis. 2010;50(6):882–9.

    Article  PubMed  Google Scholar 

  2. EFSA (European Food Safety Authority), ECDC (European Centre for Disease Prevention and Control). Trends and Sources of Zoonoses, Zoonotic Agents and Food-borne Outbreaks in 2012. EFSA Journal. 2014;12(2):3547.

    Google Scholar 

  3. Rimhanen-Finne R, Salmenlinna S, Siitonen A, Nohynek H, et al. Salmonella. In: Jaakola S, Lyytikäinen O, Rimhanen-Finne R, editors. Infection diseases in Finland 2013. Helsinki, Finland: National Institute for Health and Welfare of Finland (THL); 2014. p. 13–5. Available in:

    Google Scholar 

  4. Bern C, Martines J, de Zoysa I, Glass RI. The magnitude of the global problem of diarrhoeal disease: a ten-year update. Bull World Health Organ. 1992;70(6):705–14.

    CAS  PubMed Central  PubMed  Google Scholar 

  5. Barton Behravesh C, Jones TF, Vugia DJ, Long C, Marcus R, Smith K, et al. Deaths associated with bacterial pathogens transmitted commonly through food: foodborne diseases active surveillance network (FoodNet), 1996-2005. J Infect Dis. 2011;204(2):263–7.

    Article  PubMed  Google Scholar 

  6. Ethelberg S, Sorensen G, Kristensen B, Christensen K, Krusell L, Hempel-Jorgensen A, et al. Outbreak with multi-resistant Salmonella Typhimurium DT104 linked to carpaccio, Denmark, 2005. Epidemiol Infect. 2007;135(6):900–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Cavallaro E, Date K, Medus C, Meyer S, Miller B, Kim C, et al. Salmonella typhimurium infections associated with peanut products. N Engl J Med. 2011;365(7):601–10.

    Article  CAS  PubMed  Google Scholar 

  8. Bruun T, Sorensen G, Forshell LP, Jensen T, Nygard K, Kapperud G, et al. An outbreak of Salmonella Typhimurium infections in Denmark, Norway and Sweden, 2008. Euro Surveill. 2009;14(10):19147.

    PubMed  Google Scholar 

  9. Nygard K, Lindstedt BA, Wahl W, Jensvoll L, Kjelso C, Molbak K, et al. Outbreak of Salmonella Typhimurium infection traced to imported cured sausage using MLVA-subtyping. Euro Surveill. 2007;12(3), E070315.5.

    CAS  PubMed  Google Scholar 

  10. Takkinen J, Nakari UM, Johansson T, Niskanen T, Siitonen A, Kuusi M. A nationwide outbreak of multiresistant Salmonella Typhimurium in Finland due to contaminated lettuce from Spain, May 2005. Euro Surveill. 2005;10(6), E050630.1.

    PubMed  Google Scholar 

  11. Torres MI, Lewis P, Cook L, Cook P, Kardamanidis K, Shadbolt C, et al. An outbreak of Salmonella Typhimurium linked to a kebab takeaway shop. Commun Dis Intell Q Rep. 2012;36(1):101–6.

    PubMed  Google Scholar 

  12. Moffatt CR, Appuhamy R, Kaye A, Carswell A, Denehy D. An outbreak of Salmonella Typhimurium phage type 135a gastroenteritis linked to eggs served at an Australian Capital Territory cafe. Commun Dis Intell Q Rep. 2012;36(3):E281–5.

    PubMed  Google Scholar 

  13. Ziehm D, Dreesman J, Rabsch W, Fruth A, Pulz M, Kreienbrock L, et al. Subtype specific risk factor analyses for sporadic human salmonellosis: a case-case comparison in Lower Saxony, Germany. Int J Hyg Environ Health. 2013;216(4):428–34.

    Article  CAS  PubMed  Google Scholar 

  14. Peters TM. Pulsed-field gel electrophoresis for molecular epidemiology of food pathogens. Methods Mol Biol. 2009;551:59–70.

    Article  CAS  PubMed  Google Scholar 

  15. Wattiau P, Boland C, Bertrand S. Methodologies for Salmonella enterica subsp. enterica subtyping: gold standards and alternatives. Appl Environ Microbiol. 2011;77(22):7877–85.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Camarda A, Circella E, Pupillo A, Legretto M, Marino M, Pugliese N. Pulsed-Field Gel Electrophoresis of Salmonella enterica. Methods Mol Biol. 2015;1301:191–210.

    Article  PubMed  Google Scholar 

  17. Lindstedt BA, Torpdahl M, Vergnaud G, Le Hello S, Weill FX, Tietze E, et al. Use of multilocus variable-number tandem repeat analysis (MLVA) in eight European countries, 2012. Euro Surveill. 2013;18(4):20385.

    CAS  PubMed  Google Scholar 

  18. van Belkum A, Scherer S, van Alphen L, Verbrugh H. Short-sequence DNA repeats in prokaryotic genomes. Microbiol Mol Biol Rev. 1998;62(2):275–93.

    PubMed Central  PubMed  Google Scholar 

  19. Denoeud F, Vergnaud G. Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource. BMC Bioinformatics. 2004;5:4.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Kjeldsen MK, Torpdahl M, Campos J, Pedersen K, Nielsen EM. Multiple-locus variable-number tandem repeat analysis of Salmonella enterica subsp. enterica serovar Dublin. J Appl Microbiol. 2014;116(4):1044–54.

    Article  CAS  PubMed  Google Scholar 

  21. Boxrud D, Pederson-Gulrud K, Wotton J, Medus C, Lyszkowicz E, Besser J, et al. Comparison of multiple-locus variable-number tandem repeat analysis, pulsed-field gel electrophoresis, and phage typing for subtype analysis of Salmonella enterica serotype Enteritidis. J Clin Microbiol. 2007;45(2):536–43.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Bergamini F, Iori A, Massi P, Pongolini S. Multilocus variable-number of tandem-repeats analysis of Salmonella enterica serotype Gallinarum and comparison with pulsed-field gel electrophoresis genotyping. Vet Microbiol. 2011;149(3-4):430–6.

    Article  CAS  PubMed  Google Scholar 

  23. Harada T, Sakata J, Kanki M, Seto K, Taguchi M, Kumeda Y. Molecular epidemiological investigation of a diffuse outbreak caused by Salmonella enterica serotype Montevideo isolates in Osaka Prefecture, Japan. Foodborne Pathog Dis. 2011;8(10):1083–8.

    Article  CAS  PubMed  Google Scholar 

  24. Davis MA, Baker KN, Call DR, Warnick LD, Soyer Y, Wiedmann M, et al. Multilocus variable-number tandem-repeat method for typing Salmonella enterica serovar Newport. J Clin Microbiol. 2009;47(6):1934–8.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Tien YY, Wang YW, Tung SK, Liang SY, Chiou CS. Comparison of multilocus variable-number tandem repeat analysis and pulsed-field gel electrophoresis in molecular subtyping of Salmonella enterica serovars Paratyphi A. Diagn Microbiol Infect Dis. 2011;69(1):1–6.

    Article  CAS  PubMed  Google Scholar 

  26. Octavia S, Lan R. Multiple-locus variable-number tandem-repeat analysis of Salmonella enterica serovar Typhi. J Clin Microbiol. 2009;47(8):2369–76.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Lindstedt BA, Vardund T, Aas L, Kapperud G. Multiple-locus variable-number tandem-repeats analysis of Salmonella enterica subsp. enterica serovar Typhimurium using PCR multiplexing and multicolor capillary electrophoresis. J Microbiol Methods. 2004;59(2):163–72.

    Article  CAS  PubMed  Google Scholar 

  28. Ross IL, Davos DE, Mwanri L, Raupach J, Heuzenroeder MW. MLVA and phage typing as complementary tools in the epidemiological investigation of Salmonella enterica serovar Typhimurium clusters. Curr Microbiol. 2011;62(3):1034–8.

    Article  CAS  PubMed  Google Scholar 

  29. Best EL, Lindstedt BA, Cook A, Clifton Hadley FA, Threlfall EJ, Liebana E. Multiple-locus variable-number tandem repeat analysis of Salmonella enterica subsp. enterica serovar Typhimurium: comparison of isolates from pigs, poultry and cases of human gastroenteritis. J Appl Microbiol. 2007;103(3):565–72.

    Article  CAS  PubMed  Google Scholar 

  30. Dyet KH, Turbitt E, Carter PE. Multiple-locus variable-number tandem-repeat analysis for discriminating within Salmonella enterica serovar Typhimurium definitive types and investigation of outbreaks. Epidemiol Infect. 2011;139(7):1050–9.

    Article  CAS  PubMed  Google Scholar 

  31. Larsson JT, Torpdahl M, Petersen RF, Sorensen G, Lindstedt BA, Nielsen EM. Development of a new nomenclature for Salmonella typhimurium multilocus variable number of tandem repeats analysis (MLVA). Euro Surveill. 2009;14(15):19174.

    PubMed  Google Scholar 

  32. Siitonen A, Kuronen H, Pelkonen S. Kotimainen salmonellatilanne – katsaus, Salmonella in Finland – Overview. Suomen Eläinlääkärilehti. 2008;114(3):139–44.

    Google Scholar 

  33. Hallanvuo, S., Johansson, T. Elintarvikkeiden mikrobiologiset vaarat, Microbiological hazards in food, Evira publications, 1/2010; 67-78. Available in:

  34. Lindqvist N, Heinikainen S, Siitonen A, Pelkonen S. Molecular characterization of Salmonella enterica subsp. enterica serovar Typhimurium DT1 isolates. Epidemiol Infect. 2004;132(2):263–72.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. EFSA (European Food Safety Authority) and ECDC (European Centre for Disease Prevention and Control). EU Summary Report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2013. EFSA Journal, 2015;13(2):4036, 178 pp. Available in:

  36. European Medicines Agency: European Surveillance of Veterinary Antimicrobial Consumption (2013) Sales of veterinary antimicrobial agents in 25 EU/EEA countries in 2011. EMA/236501/2013 2013. Available in:

  37. Threlfall EJ. Epidemic salmonella typhimurium DT 104--a truly international multiresistant clone. J Antimicrob Chemother. 2000;46(1):7–10.

    Article  CAS  PubMed  Google Scholar 

  38. Wuyts V, Mattheus W, De Laminne de Bex G, Wildemauwe C, Roosens NH, Marchal K, et al. MLVA as a tool for public health surveillance of human Salmonella Typhimurium: prospective study in Belgium and evaluation of MLVA loci stability. PLoS One. 2013;8(12), e84055.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Tuominen P, Ranta J, Maijala R. Salmonella risk in imported fresh beef, beef preparations, and beef products. J Food Prot. 2006;69(8):1814–22.

    CAS  PubMed  Google Scholar 

  40. Wahlstrom H, Andersson Y, Plym−Forshell L, Pires SM. Source attribution of human Salmonella cases in Sweden. Epidemiol Infect. 2011;139(8):1246–53.

    Article  CAS  PubMed  Google Scholar 

  41. Hakanen A, Kotilainen P, Huovinen P, Helenius H, Siitonen A. Reduced fluoroquinolone susceptibility in Salmonella enterica serotypes in travelers returning from Southeast Asia. Emerg Infect Dis. 2001;7(6):996–1003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Chiou CS, Hung CS, Torpdahl M, Watanabe H, Tung SK, Terajima J, et al. Development and evaluation of multilocus variable number tandem repeat analysis for fine typing and phylogenetic analysis of Salmonella enterica serovar Typhimurium. Int J Food Microbiol. 2010;142(1-2):67–73.

    Article  CAS  PubMed  Google Scholar 

  43. Ngoi ST, Lindstedt BA, Watanabe H, Thong KL. Molecular characterization of Salmonella enterica serovar Typhimurium isolated from human, food, and animal sources in Malaysia. Jpn J Infect Dis. 2013;66(3):180–8.

    Article  PubMed  Google Scholar 

  44. Cooley MB, Carychao D, Nguyen K, Whitehand L, Mandrell R. Effects of environmental stress on stability of tandem repeats in Escherichia coli O157:H7. Appl Environ Microbiol. 2010;76(10):3398–400.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Issenhuth-Jeanjean S, Roggentin P, Mikoleit M, Guibourdenche M, de Pinna E, Nair S, et al. Supplement 2008–2010 (no. 48) to the White-Kauffmann-Le Minor scheme. Res Microbiol. 2014;165(7):526–30.

    Article  PubMed  Google Scholar 

  46. Anderson ES, Ward LR, Saxe MJ, de Sa JD. Bacteriophage-typing designations of Salmonella typhimurium. J Hyg (Lond). 1977;78(2):297–300.

    Article  CAS  Google Scholar 

  47. Rabsch W. Salmonella typhimurium phage typing for pathogens. Methods Mol Biol. 2007;394:177–211.

    Article  CAS  PubMed  Google Scholar 

  48. Clinical and Laboratory Standards Institute (CLSI): Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria that Grow Aerobically. Available in: 2009.

  49. PulseNet: One-day (24–48 h) standardized laboratory protocol for molecular subtyping of Escherichia coli O157:H7, non-typhoidal Salmonella serotypes, and Shigella sonnei by pulsed field gel electrophoresis (PFGE). Available in: 2002.

  50. European Centre of Disease Prevention and Control (ECDC): Laboratory standard operating procedure for MLVA of Salmonella enterica serotype Typhimurium. Available in: 2011.

  51. Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26(11):2465–6.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references


Funding for this research was provided by The Academy of Finland, ELVIRA project (project number 117897). The authors thank the personnel of Statens Serum Institute, Denmark. Especially Dr. Eva Møller-Nielsen and Jonas Larsson are thanked for providing the set of the Salmonella Typhimurium reference strains for harmonization of the MLVA method and helping to calculate the real numbers of repeats. Pekka Ellonen and the other personnel at FIMM are acknowledged for their help in fragment analysis. The skilful technical assistance of the staff at the Bacteriology Unit, National Institute for Health and Welfare (THL), Helsinki, Finland is gratefully acknowledged.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Taru Lienemann.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TL participated in the design of the study, did or supervised the MLVA and PFGE, carried out the data analysis and drafted the manuscript. AK performed or supervised the phage typing and the antimicrobial susceptible testing. JH did MLVA for the isolates in 2012 and helped with the MLVA data analysis. KH contributed to the study design and drafting of the manuscript. AS obtained the funding, planned the study, and edited the manuscript. All authors have read and approved the final manuscript.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lienemann, T., Kyyhkynen, A., Halkilahti, J. et al. Characterization of Salmonella Typhimurium isolates from domestically acquired infections in Finland by phage typing, antimicrobial susceptibility testing, PFGE and MLVA. BMC Microbiol 15, 131 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: