Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparison of spa Types, SCCmec Types and Antimicrobial Resistance Profiles of MRSA Isolated from Turkeys at Farm, Slaughter and from Retail Meat Indicates Transmission along the Production Chain

Abstract

The prevalence of MRSA in the turkey meat production chain in Germany was estimated within the national monitoring for zoonotic agents in 2010. In total 22/112 (19.6%) dust samples from turkey farms, 235/359 (65.5%) swabs from turkey carcasses after slaughter and 147/460 (32.0%) turkey meat samples at retail were tested positive for MRSA. The specific distributions of spa types, SCCmec types and antimicrobial resistance profiles of MRSA isolated from these three different origins were compared using chi square statistics and the proportional similarity index (Czekanowski index). No significant differences between spa types, SCCmec types and antimicrobial resistance profiles of MRSA from different steps of the German turkey meat production chain were observed using Chi-Square test statistics. The Czekanowski index which can obtain values between 0 (no similarity) and 1 (perfect agreement) was consistently high (0.79–0.86) for the distribution of spa types and SCCmec types between the different processing stages indicating high degrees of similarity. The comparison of antimicrobial resistance profiles between the different process steps revealed the lowest Czekanowski index values (0.42–0.56). However, the Czekanowski index values were substantially higher than the index when isolates from the turkey meat production chain were compared to isolates from wild boar meat (0.13–0.19), an example of a separated population of MRSA used as control group. This result indicates that the proposed statistical method is valid to detect existing differences in the distribution of the tested characteristics of MRSA. The degree of similarity in the distribution of spa types, SCCmec types and antimicrobial resistance profiles between MRSA isolates from different process stages of turkey meat production may reflect MRSA transmission along the chain.

Introduction

Staphylococcus (S.) aureus is a common cause of food poisoning due to the production of various enterotoxins. S. aureus is a frequent colonizer of the skin and mucous membranes and therefore, personnel and food-producing animals are the main sources of S. aureus in food [1]. The control of S. aureus is routinely considered in the food producing industry if standard food safety management systems are operated. In recent years, methicillin-resistant Staphylococcus aureus (MRSA), previously known as a multidrug resistant pathogen causing severe healthcare associated and community acquired infections, [2] has been observed worldwide in livestock husbandry as well as in food of different animal origins raising concerns about a possible farm to fork transmission.

First reported from pigs in the Netherlands [3] and France [4] a distinct MRSA lineage, Clonal Complex (CC) 398, has emerged in food producing animals in Europe especially in herds of pigs [5][8], veal calves [9] broiler flocks [10], [11] and turkeys [12]. Therefore, the term “livestock-associated MRSA” (LA-MRSA) was introduced considering livestock to form a new and separate reservoir for MRSA [13]. In Asian countries, however, sequence type ST9, a separate genetic linage, is predominating among MRSA isolates from livestock animals [14], [15]. Different DNA sequencing methods are used for typing MRSA strains. In order to define MRSA clones, Multilocus sequence typing (MLST), a method of classifying MRSA strains by the allelic profile of seven housekeeping genes, is used in conjunction with PCR analysis of the staphylococcal chromosomal cassette mec (SCCmec), a mobile genetic element that contains the mec A gene encoding for resistance to methicillin [16]. 11 different SCCmec types have been described, so far. The class of mec gene complex and the type of ccr gene complex carrying a set of recombinase genes responsible for integration and excision of the cassette characterize the different types of SCCmec elements [17]. Whereas SCCmec I-X harbor mecA SCCmec XI carries a divergent mecA homologue (mecALGA251) [18]. Spa typing differentiates MRSA strains by the number of tandem repeats and the sequence variation in region X of the protein A gene (spa) and can be used for reliable and discriminatory typing of MRSA [19]. As particular MLST have shown to be associated with specific repeats and repeat successions it is, with few exceptions, possible to infer an MLST type from the spa type. (http://www.spaserver.ridom.de). The frequent use of antimicrobials in animal production is suspected to facilitate the emergence and spread of MRSA due to antimicrobial selection pressure [20][22]. High stocking density in intensive food animal production holdings and intensive animal trading promote the rapid spread of MRSA between livestock populations [23], [24]. LA-MRSA strains have also been detected in raw meat at retail including beef, veal, pork and poultry [25][33] indicating potential transmission along the chain due to cross contamination during slaughter and processing. However, the extent of this transmission is so far poorly understood.

In Germany, the national monitoring for zoonotic agents aims at characterizing the prevalence of potential zoonotic pathogens at different stages of various food chains. The monitoring is part of the official control of foodstuffs and fulfills the requirements of EU Directive 2003/99/EC [34]. In 2010, the turkey meat production chain was addressed in this monitoring scheme.

The objective of the present study was to use data from the national monitoring of zoonotic agents in the food chain to obtain a comprehensive insight into the presence and transmission of MRSA in the German turkey meat production chain. A new approach is proposed for analyzing a cross sectional MRSA data set from different stages of the food chain in order to draw conclusions on potential farm to fork transmission. For this purpose, the prevalence of MRSA and the distribution of spa types, SCCmec types and antimicrobial resistance profiles among MRSA isolated from different steps of the turkey meat production chain were compared. It is proposed that the degree of similarity in the distribution of spa types, SCCmec types and antimicrobial resistance profiles between the samples from the three process steps may be interpreted as reflecting MRSA transmission along the chain.

Materials and Methods

1. Study Design

Sampling was conducted in 2010 by the competent authorities of the federal states according to a pre-defined protocol in the framework of the national monitoring for zoonotic agents. All participating competent authorities are listed in table S1. Dust samples from 112 German turkey flocks were collected in order to quantify the presence of MRSA in primary production and to assess the introduction of MRSA into the slaughterhouses. Samples at slaughterhouses (n = 359) were analyzed to estimate the transfer to carcasses during slaughter and to determine the transmission of MRSA from carcasses to fresh turkey meat during further processing. Finally, 460 turkey meat portions were sampled to evaluate the MRSA exposure of consumers via contaminated turkey meat.

Turkey pens were sampled by pooling 5 dust swab samples, collected from different sections representing an area of 500 cm2, each. At the slaughterhouse, at least 30 g neck skin was sampled from turkey carcasses after slaughter and chilling, but prior to further processing. Samples of 25 g of fresh turkey meat (with or without skin) were collected at retail. In order to ensure a high level of representativity, the distribution of the samples in primary production and at slaughter across Germany was proportional to the number of turkey flocks and the slaughter capacity of the respective federal state. Meat samples at retail were distributed according to the human population size of the executive federal state. A more detailed description of the principles of the national monitoring for zoonotic agents has been published before [35].

2. MRSA Isolation

MRSA were isolated by the regional laboratories according to the recommended method of the National Reference Laboratory (NRL) for staphylococci including S. aureus at the Federal Institute for Risk Assessment (BfR). The dust samples were pooled per turkey house in 100 ml Mueller Hinton broth supplemented with 6.5% NaCl for pre-enrichment. Neck skin samples (at least 30 g), fresh meat (25 g) and meat preparations (25 g) were pre-enriched in 225 ml Mueller Hinton broth supplemented with 6.5% NaCl. After incubation for 16–20 h at 37°C, 1 ml pre-enrichment broth was transferred into 9 ml of tryptic soy broth supplemented with 50 mg/l aztreonam and 3.5 mg/l cefoxitin. After incubation of this selective-enrichment broth for a further 16–20 h at 37°C one loopful was plated onto sheep blood agar and chromogenic MRSA screening agar respectively, and incubated for 24–48 h at 37°C. Presumptive MRSA isolates were sent to the National Reference Laboratory (NRL) for staphylococci including S. aureus at the Federal Institute for Risk Assessment (BfR) for MRSA confirmation and characterization. The number of MRSA isolates included in further analyses is not exactly congruent to the amount of positive samples obtained within the national monitoring for zoonotic agents because first, the NRL did not always receive the corresponding isolate from the competent authorities of the federal states or second, isolates which did not exactly correspond to the monitoring sampling plan in terms of completeness of data reporting to the national level but were obtained from the correct matrix were excluded from prevalence estimations but included in further typing and strain comparisons.

Twenty one MRSA isolates from wild boar meat within the national monitoring for zoonotic agents of 2011 were used in the analyses as a control group (data not shown in detail). The control group was selected to ensure wide differences with the population under study concerning the distribution of MRSA strains in order to evaluate if the used analytical approach is appropriate to differentiate between the matrices.

3. Molecular Typing

Presumptive MRSA isolates were confirmed by an in-house multiplex PCR simultaneously targeting the 23S rDNA specific for Staphylococcus species [36], the nuclease gene nuc which is specific for S. aureus, and the resistance gene mecA [37]. Template DNA was extracted using the “RTP Bacteria DNA Mini Kit” (Invitek, Berlin, Germany). All MRSA isolates were further characterized using spa typing [38] and SCCmec-typing [39]. The method applied for typing of the SCCmec differentiates SCCmec types I to V and their subtypes. However, isolates of the CC398 characterized as type III by the method have been shown to rather be a variant of type V [40]. The software Ridom Staphytype (Ridom GmbH, Würzburg, Germany) was used to assign spa types. Spa types which have not been identified and assigned to a clonal complex (CC) by the NRL before were additionally subjected to multilocus sequence typing (MLST) [41].

4. Antimicrobial Susceptibility Testing

All isolates were tested for the susceptibility to antimicrobials using broth microdilution in accordance with Clinical and Laboratory Standards Institute guidelines [42]. Commercial microtitre plates were used (TREK Diagnostic Systems, Magellan Biosciences, West Sussex, England). Minimum inhibitory concentrations (MIC) were evaluated according to epidemiological cut-off values (ECOFFs) published for MRSA and S. aureus by the European committee for antimicrobial susceptibility testing (www.eucast.org). MIC values above the ECOFFs indicated microbiological resistance. MIC lower or equal to the ECOFFs characterised susceptible strains. S. aureus strain ATCC 25923 was used for quality assurance Resistance testing included gentamicin, kanamycin, streptomycin, chloramphenicol, ciprofloxacin, tetracycline, clindamycin, erythromycin, mupirocin, linezolid, vancomycin, quinupristin/dalfopristin, penicillin, fusidic acid, cefoxitin, trimethoprim, sulfamethoxazole, rifampicin and tiamulin.

5. Statistical Analysis

The chi square test of homogeneity was used to analyze differences in the distribution of spa types and antibiotic resistance profiles between MRSA strains from the turkey flocks, carcasses at slaughter and meat. Isolates were grouped according to their spa types and antibiotic resistance profiles to assure appropriate numbers of isolates in all categories. All spa types were aggregated in accordance to their frequency of occurrence. The phenotypic antimicrobial resistance profiles were grouped by hierarchical cluster analysis using Ward’s minimum variance and squared Euclidean distance. The MIC values for each isolate were categorized into resistant or susceptible according to the ECOFFs to generate a binary data set. The final amount of clusters was determined using the Pseudo-F [43] and Pseudo-T [44] statistics. Both tests indicate possible breakpoints for splitting the data into the appropriate amount of clusters. The distribution of SCCmec types in the different matrices were compared using Fisher’s exact test as 33.3% of the cells of the contingency table had an expected value below 5. P-values of <0.05 were considered statistically significant. Chi square test, Fisher’s exact test and cluster analysis were calculated using the statistical software package SPSS 18.0 (SPSS Inc. Munich, Germany). Pseudo-F and Pseudo-T statistics were performed using SAS/STAT software 9.2 (SAS Institute Inc., Cary, NC, USA).

The degree of similarity between the frequency distributions of spa types, SCCmec types and resistance profiles of MRSA among the sample sets from the turkey primary production, carcasses at slaughterhouse and turkey meat at retail was estimated using the Czekanowski index or proportional similarity index (PSI) [45]. It is calculated by:where pi and qi represent the proportion of strains out of all strains among the data sets P and Q which agree in the realization i of the variable of interest. The values for PS range from 1 for identical frequency distributions of the variable of interest to zero for no similarities between the data sets. Since the size of the samples is rather small, a realization of the PSI index may deviate largely from its true value. Thus, the PSI was bootstrapped obtaining a probability density distribution from which we derived the 95% confidence interval for the PSI. The statistic open source software R (available at: http://www.R-project.org) was used to calculate the approximate confidence interval of the Czekanowski index using the non-parametric boostrap BCα method utilizing 2000 iterations [46].

Results

Twenty two (19.6%) of 112 dust samples from the turkey primary production, 235 (65.5%) of 359 turkey carcasses after slaughter and 147 (32.0%) of 460 turkey meat samples at retail were tested positive for MRSA [47]. A set of 32 isolates from dust samples, 248 isolates from turkey carcasses and 241 isolates from turkey meat was used for further laboratory analyses (Table 1).

thumbnail
Table 1. MRSA prevalence and distribution of spa types, SCCmec types and antimicrobial resistance clusters at different steps of the German turkey meat production chain in 2010.

https://doi.org/10.1371/journal.pone.0096308.t001

A total of 16 different spa types were identified. The number of different spa types increased during processing from 5 different types in dust samples over 8 in carcasses to 15 different types in meat samples. The proportion of strains assigned to CC398 ranged between 85.9 and 90.6%. Among CC398, t011 (43.8–46.9%) and t034 (32.0–43.8%) were the predominating spa-types on every process step. Spa types t1430 (4.0–6.3%) and t002 (3.1–9.1%) were dominating within the group of non CC398 strains.

Most of the strains carried SCCmec-type V (58.1–71.9%) followed by type IVa (19–27.0%). Type III (0–1.2%) was identified sporadically (Table 1). However, there is evidence in former literature that CC398 strains which were identified as SCCmec type III by the typing scheme of Zhang et al. [39] are rather assigned to a separate variant of SCCmec type V [40], [48], [49]. In 5.7–17.6% of the strains the SCCmec type could not be identified by the method used.

Susceptibility to 19 different antimicrobial agents was determined (Figure 1). Throughout the turkey production chain, the vast majority of isolates was resistant to tetracycline (98.8%–100%). High resistance rates were obtained to clindamycin (79.4–93.8%), erythromycin (73.8–87.5%), trimethoprim (65.7–78.1%), quinupristin/dalfopristin (62.2–66.1%) and tiamulin (52.3–65.6%). Resistances to mupirocin, linezolid, sulfamethoxazole and rifampicin were observed sporadically in individual isolates from all steps of the process chain. All isolates were susceptible to vancomycin. Resistance to tiamulin (62.2 versus 8.2%), gentamicin (25.2 versus 6.6%) and trimethoprim (72.0 versus 36.1%) was considerably more frequent among CC398 than among non-CC398 strains. Resistance to ciprofloxacin was common among non-CC398 strains (98.4 versus 26.1% in CC398 strains) (Figure 1).

thumbnail
Figure 1. Antimicrobial resistance of MRSA in the German turkey meat production chain.

Distribution of antimicrobial resistance of MRSA strains separated into CC398 and non CC 398 strains as well as different steps of the turkey meat production chain isolated from dust samples at turkey primary production (n = 32), carcasses at slaughter (n = 248) and meat at retail (n = 241). The MRSA strains were isolated in the course of the national monitoring for zoonotic agents in Germany in 2010.

https://doi.org/10.1371/journal.pone.0096308.g001

All 521 MRSA strains were included in further similarity estimations. In accordance to the frequency of their occurrence all spa types were aggregated in 4 different categories for further statistical analysis. The most prevalent spa types t011 and t034 built their own group whereas rare spa types of CC398 and all non CC398 strains were summarized in separate groups. The chi square distribution of the spa type groups did not significantly differ between primary production, carcasses at slaughter and meat at retail (p = 0.06). Likewise, no significant difference was identified in the distribution of SCCmec types between the origins using fisher’s exact test (p = 0.095). A total of 101 different resistance profiles were identified among the MRSA isolates including resistance to 2 to 12 different antimicrobial substances. The hierarchical cluster algorithm of Wards minimum variance combined with squared Euclidean distance separated the antimicrobial resistance profiles into homogenous clusters. Identical resistance phenotypes did not appear in more than one cluster. Based on the Pseudo-F and Pseudo-T statistics the 3 cluster solution containing 33, 44 and 24 different phenotypic resistance profiles, respectively, was identified to best describe the binary data set. Detailed characteristics of the cluster composition, concerning antimicrobial resistance and the distribution of groups of spa types and SCCmec types, is summarized in table S2. The antimicrobial resistance clusters did not significantly differ in their chi square distribution between the MRSA samples from the three origins (p = 0.295).

The distribution of spa types, SCCmec types and antimicrobial resistance profiles within the sample collections from the three process steps and the control group were compared pair wise using the Czekanowski index (Table 2). High index values were obtained for the distribution of spa types (PSI 0.79–0.86) among MRSA from the turkey meat chain. The comparison of the distribution of antimicrobial resistance profiles resulted in the lowest index values (PSI 0.42–0.56). The distribution of spa types and antimicrobial resistance profiles showed remarkably higher similarity between the different production steps of the turkey meat chain as to samples from the control group (PSI 0.55–0.56 and 0.13–0.19 resp.). High similarity in the distributions of SCCmec types was calculated between all process steps of the turkey meat production chain (PSI 0.85–0.91). However, a strong association was also received with SCCmec types of the control group (PSI 0.83–0.85).

thumbnail
Table 2. Similarity matrix of spa types, SCCmec types and resistance profiles of MRSA isolated from the German turkey meat production chain in the course of the national monitoring for zoonotic agents in 2010 (95% confidence intervals).

https://doi.org/10.1371/journal.pone.0096308.t002

Discussion

In the present study, a new approach is proposed for analyzing a cross sectional set of MRSA isolates originating from three consecutive stages of the turkey meat production chain in order to draw conclusions on a potential farm to fork transmission. In the course of the German national monitoring for zoonotic agents in 2010 MRSA was isolated at all stages of the turkey meat production chain with prevalences ranging from 19.6% to 65.5%. To our knowledge, this is the first representative national MRSA prevalence study in the turkey production chain. In a regional prevalence study among fattening turkeys in southern Germany in 2009, a considerably higher prevalence of 90% MRSA positive flocks was observed using the same sampling procedure [12]. The difference might be explained by the regional restriction of sampling and the small sample size in that study. The proportion of positive meat samples is in line with results from the Netherlands [25]. Outside of Europe, low MRSA contamination rates of 3.85% [32] and 1.7% [50] were reported among US turkey meat.

The high MRSA prevalence in turkey carcasses after slaughter in comparison to the flock prevalence is in contrast to the situation in pigs [12], [51] and indicates that the turkey slaughter process may play an important role in the transmission of MRSA. Turkeys are slaughtered highly automated at a speed of line up to 3,600 turkey hens and up to 2,700 turkey toms per hour which leads to a permanent introduction of MRSA into the poultry processing plants [52]. During the process, MRSA on animal surfaces can get transmitted via direct contact or indirect via surface processing machinery, scalding water or the hands of staff. Scalding takes place at a constant water temperature between 50 and 65°C for 60 to 210 sec [52]. Although the surface of the carcasses is exposed to a heat treatment during scalding, the temperature and duration of the process might be insufficient to substantially reduce superficial MRSA counts. The selective growth of S. aureus after the elimination of less heat resistant microbial flora in the scalding water has been discussed [53]. As bacterial counts increase in the tanks throughout the slaughter day scalding can contribute to cross contamination [54]. After scalding, the birds go through the plucking machines consisting of revolving drums with rubber beaters or discs with plucking fingers. The birds are flailed and scraped for 30–90 sec while being sprayed with warm or cold water [52]. Plucking equipment is difficult to clean and a persisting microbiological flora can get established [55]. Cross contamination during slaughter and meat processing might lead to an extensive distribution of spa types between different animals and slaughter flocks. In addition, the increase in manual handling during processing facilitates the entry of human MRSA strains into the production units. This can explain the increase in the variability of spa types along the chain and is in line with the increase in the proportion of non CC398 strains in meat samples compared to dust or carcasses. Spa types t002 and t1430 were also present in primary production and therefore probably have been transmitted along the food chain. In contrast, spa types t010, t015 were first observed in meat samples.

The majority of MRSA from the German turkey production chain was assigned to the livestock associated CC398 with the predominant spa types t011 and t034. This is in line with results from other livestock like veal calves [9], dairy cattle [56], [57] and pigs [6] as well as in food [25]. In the present study, 37 of the 521 MRSA strains (7.1%) were identified as t002. This spa type t002 is assigned to CC5. In Germany, CC5 is one of the epidemic MRSA strains among humans [58]. Finding t002 in turkey flocks and in turkey meat is in line with other studies from central Europe [12], [25], [59]. So far, it is not known, whether this strain originates from the “human” strain and is introduced into the food chain on different levels or whether it got established in the turkey population and is transmitted along the chain. Detailed molecular-epidemiological investigations are needed to compare strains both from human and farm to fork origin. In the present study, 4.2% of the MRSA isolates were characterized as spa type t1430, a MRSA strain which was also frequently isolated from chicken meat [25] and broilers at slaughter [60] in the Netherlands. However, it was has also been detected in turkey flocks at farm level [12]. The strain is assigned to ST9, a lineage genetically unrelated to ST398. ST9 is the predominating sequence type among MRSA from pigs in Asian countries [14], [15], [61][64]. Outside of Europe, MRSA contamination was reported among US turkey meat [32], [50]. In both surveys, all isolates belonged to USA 300 (ST8), the most common community associated MRSA strain in the USA, suggesting human contamination during processing.

The frequent use of antimicrobials at farm is discussed as a risk factor for the wide dissemination of MRSA in livestock production chains [65]. In recent studies antimicrobials were identified to be used in more than 90% of the investigated turkey flocks and animals received on average 33 daily doses of antimicrobials during raising and fattening [66]. With a share of 21% β-lactams were most often used followed by polypeptides (15.2%), macrolides (13.4%), tetracyclines and aminoglycosides (12.4% both). Fluoroquinolones were used in 6.5% of the investigated flocks. The common application of antimicrobials via drinking water bears the risk of under dosing of individual animals and contamination of the barn environment with antimicrobials which also facilitates the selection of resistance [67].

Cluster analysis was used to better describe the multidimensional data set of antibiotic resistance profiles grouping all MRSA strains within 3 different clusters. As the ordinal MIC values generated by two-fold dilutions in substance concentration are difficult to describe by cluster analysis a binary interpretation of the data set was used. Ward’s minimum variance with squared Euclidian distance was proven to be the best method to produce well separated cluster in binary antimicrobial resistance data sets [68], [69]. No resistance phenotype simultaneously appeared in several clusters. The distribution of spa types, SCCmec types and the three clusters of antimicrobial resistance types did not significantly differ in the MRSA samples from the three origins. The chi square value was approaching significance with respect to the spa-types, which was presumably due to the slightly higher proportion of other CC398 and non CC398. However, considering all three features it cannot be rejected on the basis of the included data that the MRSA isolates from different steps of the turkey meat production chain originate from the same population of strains. This result might rather indicate farm to fork transmission of MRSA of the same pool of strains than development of separate MRSA populations at each step of the chain. The calculation of the Czekanowski index for spa type and SCCmec type data results in consistently high similarity values between the matrices whereas the comparison of antimicrobial resistance phenotypes observed medium index values. Higher values of similarity were obtained between the adjacent process steps primary production/slaughter and slaughter/meat than between samples from primary production and meat. This result was expected as an increase in the variability of the MRSA isolates might be conceivable at each process stage due to external introduction of new strains via human or environmental contamination or due to spontaneous mutations in the strains.

The lower values of similarity between the distribution of spa types and antimicrobial resistance profiles of samples from the turkey meat production chain and the control group indicate that that the proposed statistical method is valid to detect existing differences in the distribution of these characteristics of MRSA.

Concerning SCCmec types, high index values were also observed in comparison to the control group which might be explained by the insufficient discriminatory power of SCCmec typing. In addition, MRSA isolates with not typeable SCCmec cassettes were considered as equal that might lead to an overestimation of similarity.

It can be concluded that MRSA is present at every step of the turkey meat production chain in Germany. Using the Czekanowski index it is possible to quantify the similarity of the distribution of spa types, SCCmec types and antimicrobial resistance phenotypes between MRSA data sets from different stages of turkey meat production chain. Combined with chi square statistics, the high level of similarity suggests MRSA transmission along the chain.

Supporting Information

Table S1.

List of competent authorities of the German federal states which were responsible for collecting the samples.

https://doi.org/10.1371/journal.pone.0096308.s001

(PDF)

Table S2.

Distribution of resistance against 19 different antimicrobials, grouped spa types and SCCmec types within the binary phenotypic resistance clusters of 521 MRSA isolates sampled at different steps of the German turkey meat production chain in 2010.

https://doi.org/10.1371/journal.pone.0096308.s002

(PDF)

Author Contributions

Conceived and designed the experiments: AK BAT KA. Performed the experiments: BK. Analyzed the data: BV. Contributed reagents/materials/analysis tools: JB. Wrote the paper: BV. Designed the software used for bootstrap method: JB. Administration and organisation of the project: AF.

References

  1. 1. Hennekinne JA, De Buyser ML, Dragacci S (2012) Staphylococcus aureus and its food poisoning toxins: Characterization and outbreak investigation. FEMS Microbiology Reviews 36: 815–836.
  2. 2. Köck R, Becker K, Cookson B, van Gemert-Pijnen JE, Harbarth S, et al. (2010) Methicillin-resistant Staphylococcus aureus (MRSA): burden of disease and control challenges in Europe. Euro surveillance: bulletin européen sur les maladies transmissibles  =  European communicable disease bulletin 15: 19688.
  3. 3. Voss A, Loeffen F, Bakker J, Klaassen C, Wulf M (2005) Methicillin-resistant Staphylococcus aureus in Pig Farming. Emerging Infectious Diseases 11: 1965–1966.
  4. 4. Armand-Lefevre L, Ruimy R, Andremont A (2005) Clonal comparison of Staphylococcus from healthy pig farmers, human controls, and pigs. Emerging Infectious Diseases 11: 711–714.
  5. 5. Crombé F, Willems G, Dispas M, Hallin M, Denis O, et al. (2012) Prevalence and antimicrobial susceptibility of methicillin-resistant Staphylococcus aureus among pigs in Belgium. Microbial Drug Resistance 18: 125–131.
  6. 6. EFSA (2009) Analysis of the baseline survey on the prevalence of methicillin-resistant Staphylococcus aureus (MRSA) in holdings with breeding pigs, in the EU, 2008, Part A: MRSA prevalence estimates; on request from the European Commission. The EFSA Jounal 11: 1376.
  7. 7. Khanna T, Friendship R, Dewey C, Weese JS (2008) Methicillin resistant Staphylococcus aureus colonization in pigs and pig farmers. Vet Microbiol 128: 298–303.
  8. 8. Smith TC, Male MJ, Harper AL, Kroeger JS, Tinkler GP, et al.. (2009) Methicillin-resistant Staphylococcus aureus (MRSA) strain ST398 is present in midwestern U.S. swine and swine workers. PloS one 4.
  9. 9. Graveland H, Wagenaar JA, Heesterbeek H, Mevius D, van Duijkeren E, et al.. (2010) Methicillin resistant Staphylococcus aureus ST398 in veal calf farming: human MRSA carriage related with animal antimicrobial usage and farm hygiene. PloS one 5.
  10. 10. Nemati M, Hermans K, Lipinska U, Denis O, Deplano A, et al. (2008) Antimicrobial resistance of old and recent Staphylococcus aureus isolates from poultry: First detection of livestock-associated methicillin-resistant strain ST398. Antimicrobial Agents and Chemotherapy 52: 3817–3819.
  11. 11. Persoons D, Van Hoorebeke S, Hermans K, Butaye P, De Kruif A, et al. (2009) Methicillin-resistant Staphylococcus aureus in poultry. Emerging Infectious Diseases 15: 452–453.
  12. 12. Richter A, Sting R, Popp C, Rau J, Tenhagen BA, et al. (2012) Prevalence of types of methicillin-resistant Staphylococcus aureus in turkey flocks and personnel attending the animals. Epidemiology and Infection 140: 2223–2232.
  13. 13. Reischl U, Frick J, Hoermansdorfer S, Melzl H, Bollwein M, et al.. (2009) Single-nucleotide polymorphism in the SCCmec-orfX junction distinguishes between livestock-associated MRSA CC398 and human epidemic MRSA strains. Euro surveillance: bulletin european sur les maladies transmissibles  =  European communicable disease bulletin 14.
  14. 14. Anukool U, O’Neill CE, Butr-Indr B, Hawkey PM, Gaze WH, et al. (2011) Meticillin-resistant Staphylococcus aureus in pigs from Thailand. International Journal of Antimicrobial Agents 38: 86–87.
  15. 15. Wagenaar JA, Yue H, Pritchard J, Broekhuizen-Stins M, Huijsdens X, et al. (2009) Unexpected sequence types in livestock associated methicillin-resistant Staphylococcus aureus (MRSA): MRSA ST9 and a single locus variant of ST9 in pig farming in China. Veterinary Microbiology 139: 405–409.
  16. 16. Cookson BD, Robinson DA, Monk AB, Murchan S, Deplano A, et al. (2007) Evaluation of molecular typing methods in characterizing a European collection of epidemic methicillin-resistant Staphylococcus aureus strains: the HARMONY collection. J Clin Microbiol 45: 1830–1837.
  17. 17. (IWG-SCC) IWGotCoSCCE (2009) Classification of staphylococcal cassette chromosome mec (SCCmec): guidelines for reporting novel SCCmec elements. Antimicrobial Agents and Chemotherapy 53: 4961–4967.
  18. 18. Garcia-Alvarez L, Holden MT, Lindsay H, Webb CR, Brown DF, et al. (2011) Meticillin-resistant Staphylococcus aureus with a novel mecA homologue in human and bovine populations in the UK and Denmark: a descriptive study. Lancet Infect Dis 11: 595–603.
  19. 19. Frenay HM, Bunschoten AE, Schouls LM, van Leeuwen WJ, Vandenbroucke-Grauls CM, et al. (1996) Molecular typing of methicillin-resistant Staphylococcus aureus on the basis of protein A gene polymorphism. Eur J Clin Microbiol Infect Dis 15: 60–64.
  20. 20. de Neeling AJ, van den Broek MJ, Spalburg EC, van Santen-Verheuvel MG, Dam-Deisz WD, et al. (2007) High prevalence of methicillin resistant Staphylococcus aureus in pigs. Vet Microbiol 120: 366–372.
  21. 21. Pires SM, Evers EG, van PW, Ayers T, Scallan E, et al. (2009) Attributing the human disease burden of foodborne infections to specific sources. Foodborne PathogDis 6: 417–424.
  22. 22. WHO (2001) World Health Organisation Global Principles for the Containment of Antimicrobial Resistance in Animals intended for Food. http://wwwwhoint/drugresistance/WHO_Global_Strategyhtm/en/.
  23. 23. Alt K, Fetsch A, Schroeter A, Guerra B, Hammerl J, et al. (2011) Factors associated with the occurrence of MRSA CC398 in herds of fattening pigs in Germany. BMC Veterinary Research 7: 69.
  24. 24. Broens EM, Graat EAM, Van der Wolf PJ, van de Giessen AW, de Jong MCM (2011) Prevalence and risk factor analysis of livestock associated MRSA-positive pig herds in The Netherlands. Preventive Veterinary Medicine 102: 41–49.
  25. 25. de Boer E, Zwartkruis-Nahuis JTM, Wit B, Huijsdens XW, de Neeling AJ, et al. (2009) Prevalence of methicillin-resistant Staphylococcus aureus in meat. International Journal of Food Microbiology 134: 52–56.
  26. 26. Hanson BM, Dressler AE, Harper AL, Scheibel RP, Wardyn SE, et al. (2011) Prevalence of Staphylococcus aureus and methicillin-resistant Staphylococcus aureus (MRSA) on retail meat in Iowa. Journal of Infection and Public Health 4: 169–174.
  27. 27. Lim SK, Nam HM, Park HJ, Lee HS, Choi MJ, et al. (2010) Prevalence and characterization of methicillin-resistant Staphylococcus aureus in raw meat in Korea. Journal of Microbiology and Biotechnology 20: 775–778.
  28. 28. Lozano C, López M, Gómez-Sanz E, Ruiz-Larrea F, Torres C, et al. (2009) Detection of methicillin-resistant Staphylococcus aureus ST398 in food samples of animal origin in Spain. Journal of Antimicrobial Chemotherapy 64: 1325–1326.
  29. 29. O’Donoghue M, Chan M, Ho J, Moodley A, Boost M. Prevalence of Mehicillin-Resistant Staphylococcus aureus in Meat from Hong Kong Shops and Markets. 27.
  30. 30. Pu S, Han F, Ge B (2009) Isolation and characterization of methicillin-resistant Staphylococcus aureus strains from louisiana retail meats. Applied and Environmental Microbiology 75: 265–267.
  31. 31. Van Loo IHM, Diederen BMW, Savelkoul PHM, Woudenberg JHC, Roosendaal R, et al. (2007) Methicillin-resistant Staphylococcus aureus in meat products, the Netherlands. Emerging Infectious Diseases 13: 1753–1755.
  32. 32. Waters AE, Contente-Cuomo T, Buchhagen J, Liu CM, Watson L, et al. (2011) Multidrug-resistant staphylococcus aureus in US meat and poultry. Clinical Infectious Diseases 52: 1227–1230.
  33. 33. Weese JS, Reid-Smith R, Rousseau J, Avery B (2010) Methicillin-resistant Staphylococcus aureus (MRSA) contamination of retail pork. Canadian Veterinary Journal-Revue Veterinaire Canadienne 51: 749–752.
  34. 34. EC (2003) Directive 2003/99/EC of the European Parliament and of the Council of 17 November 2003 on the monitoring of zoonoses and zoonotic agents. Official Journal of the European Union, L325/31. http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2003:325:0031:0040:EN:PDF.
  35. 35. Käsbohrer A, Wegeler C, Tenhagen BA (2009) EU-weite und nationale Monitoringprogramme zu Zoonoseerregern in Deutschland. Journal für Verbraucherschutz und Lebensmittelsicherheit. EU-weite und nationale Monitoringprogramme zu Zoonoseerregern in Deutschland. 41–45.
  36. 36. Straub JA, Hertel C, Hammes WP (1999) A 23S rDNA-targeted polymerase chain reaction-based system for detection of Staphylococcus aureus in meat starter cultures and dairy products. Journal of Food Protection 62: 1150–1156.
  37. 37. Poulsen AB, Skov R, Pallesen LV (2003) Detection of methicillin resistance in coagulase-negative staphylococci and in staphylococci directly from simulated blood cultures using the EVIGENE MRSA Detection Kit. Journal of Antimicrobial Chemotherapy 51: 419–421.
  38. 38. Shopsin B, Gomez M, Montgomery SO, Smith DH, Waddington M, et al. (1999) Evaluation of protein A gene polymorphic region DNA sequencing for typing of Staphylococcus aureus strains. Journal of Clinical Microbiology 37: 3556–3563.
  39. 39. Zhang K, McClure JA, Elsayed S, Louie T, Conly JM (2005) Novel multiplex PCR assay for characterization and concomitant subtyping of staphylococcal cassette chromosome mec types I to V in methicillin-resistant Staphylococcus aureus. Journal of Clinical Microbiology 43: 5026–5033.
  40. 40. Jansen MD, Box ATA, Fluit AC (2009) SCCmec typing in methicillin-resistant Staphylococcus aureus strains of animal origin. Emerging Infectious Diseases 15: 136.
  41. 41. Enright MC, Day NPJ, Davies CE, Peacock SJ, Spratt BG (2000) Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. Journal of Clinical Microbiology 38: 1008–1015.
  42. 42. CLSI (2006) Performance Standards for Antimicrobial Disk Susceptibility Tests; Approved Standard.
  43. 43. Calinski T, Harabasz J (1974) A dendrite method for cluster analysis. Communications in Statistics: Taylor & Francis. 1–27.
  44. 44. Duda RO, Hart PE (1973) Pattern Classification and Scene Analysis. New York: Wiley.
  45. 45. Rosef O, Kapperud G, Lauwers S, Gondrosen B (1985) Serotyping of Campylobacter jejuni, Campylobacter coli, and Campylobacter laridis from domestic and wild animals. Applied and Environmental Microbiology 49: 1507–1510.
  46. 46. Efron B, Tibshirani R (1986) Bootstrap methods for standard errors, confidence invervals and other measures of statistical accuracy. Statistical Science 1: 54–75.
  47. 47. Bundesamt für Verbraucherschutz und Lebensmittelsicherheit (2012) Berichte zur Lebensmittelsicherheit 2010. BVL Reporte http://www.bvl.bund.de/SharedDocs/Downloads/01_Lebensmittel/04_Zoonosen_Monitoring/Zoonosen_Monitoring_Bericht_2010.pdf?__blob=publicationFile&v=6.
  48. 48. Kreausukon K, Fetsch A, Kraushaar B, Alt K, Müller K, et al. (2012) Prevalence, antimicrobial resistance, and molecular characterization of methicillin-resistant Staphylococcus aureus from bulk tank milk of dairy herds. Journal of Dairy Science 95: 4382–4388.
  49. 49. Argudin MA, Rodicio MR, Guerra B (2010) The emerging methicillin-resistant Staphylococcus aureus ST398 clone can easily be typed using the Cfr9I SmaI-neoschizomer. Lett Appl Microbiol 50: 127–130.
  50. 50. Bhargava K, Wang X, Donabedian S, Zervos M, da Rocha L, et al. (2011) Methicillin-resistant staphylococcus aureus in retail meat, Detroit, Michigan, USA. Emerging Infectious Diseases 17: 1135–1137.
  51. 51. Lassok B, Tenhagen BA (2013) From pig to pork: Methicillin-resistant staphylococcus aureus in the pork production chain. Journal of Food Protection 76: 1095–1108.
  52. 52. Löhren U (2012) Overview on current practices of poultry slaughtering and poultry meat inspection Supporting Publications. 2012: EN–298 http://wwwefsaeuropaeu/en/supporting/doc/298epdf.
  53. 53. Hentschel S, Kusch D, Sinell HJ (1979) Staphylococcus aureus in Poultry-biochemical characteristics, antibiotic resistance and phage pattern. Zentralblatt für Bakteriologie, Parasitenkunde, Infektionskrankheiten und Hygiene. Erste Abteilung Originale. Reihe B: Hygiene, Betriebshygiene, präventive Medizin 07, 168 (5–6): 546–561.
  54. 54. Groβklaus D, Lessing G (1972) Hygieneprobleme beim Schlachtgeflügel. Fleischwirtschaft 52: 1011–1013.
  55. 55. Berrang ME, Buhr RJ, Cason JA, Dickens JA (2001) Broiler carcass contamination with Campylobacter from feces during defeathering. Journal of Food Protection 64: 2063–2066.
  56. 56. Spohr M, Rau J, Friedrich A, Klittich G, Fetsch A, et al. (2011) Methicillin-resistant Staphylococcus aureus (MRSA) in three dairy herds in southwest Germany. Zoonoses and Public Health 58: 252–261.
  57. 57. Vanderhaeghen W, Hermans K, Haesebrouck F, Butaye P (2010) Methicillin-resistant Staphylococcus aureus (MRSA) in food production animals. Epidemiol Infect 138: 606–625.
  58. 58. Robert Koch-Institut (2011) Auftreten und Verbreitung von MRSA in Deutschland 2010. Epidemiologisches Bulletin 26: 233–244.
  59. 59. Feβler AT, Kadlec K, Hassel M, Hauschild T, Eidam C, et al. (2011) Characterization of methicillin-resistant Staphylococcus aureus isolates from food and food products of poultry origin in Germany. Applied and Environmental Microbiology 77: 7151–7157.
  60. 60. Mulders MN, Haenen APJ, Geenen PL, Vesseur PC, Poldervaart ES, et al. (2010) Prevalence of livestock-associated MRSA in broiler flocks and risk factors for slaughterhouse personnel in the Netherlands. Epidemiology and Infection 138: 743–755.
  61. 61. Cui S, Li J, Hu C, Jin S, Li F, et al. (2009) Isolation and characterization of methicillin-resistant Staphylococcus aureus from swine and workers in China. Journal of Antimicrobial Chemotherapy 64: 680–683.
  62. 62. Larsen J, Imanishi M, Hinjoy S, Tharavichitkul P, Duangsong K, et al.. (2012) Methicillin-resistant Staphylococcus aureus ST9 in pigs in Thailand. PloS one 7.
  63. 63. Neela V, Zafrul AM, Mariana NS, Van Belkum A, Liew YK, et al. (2009) Prevalence of ST9 methicillin-resistant Staphylococcus aureus among pigs and pig handlers in Malaysia. Journal of Clinical Microbiology 47: 4138–4140.
  64. 64. Tsai HY, Liao CH, Cheng A, Liu CY, Huang YT, et al. (2012) Isolation of meticillin-resistant Staphylococcus aureus sequence type 9 in pigs in Taiwan. International Journal of Antimicrobial Agents 39: 449–451.
  65. 65. Schwarz S, Chaslus-Dancla E (2001) Use of antimicrobials in veterinary medicine and mechanisms of resistance. Veterinary Research 32: 201–225.
  66. 66. State Office for Consumer Protection and Food Safety in Lower Saxony G (2011) Bericht über den Antibiotikaeinsatz in der landwitschaftlichen Nutztierhaltung in Niedersachsen November 2011. http://wwwmlniedersachsende/portal/livephp?navigation_id=27751&article_id=102202&_psmand=7.
  67. 67. Richter A, Hafez HM, Böttner A, Gangl A, Hartmann K, et al.. (2009) Applications of antibiotics in poultry. Tierarztliche Praxis Ausgabe G: Grosstiere – Nutztiere Verabreichung von Antibiotika in Geflügelbeständen. 321–329.
  68. 68. Berge ACB, Atwill ER, Sischo WM (2003) Assessing antibiotic resistance in fecal Escherichia coli in young calves using cluster analysis techniques. Preventive Veterinary Medicine 61: 91–102.
  69. 69. Milligan GW (1981) A monte carlo study of thirty internal criterion measures for cluster analysis. Psychometrika 46: 187–199.