Comparison of Fingerprinting Methods for Typing Methicillin-Resistant ...

2 downloads 0 Views 4MB Size Report
May 8, 2009 - ylococcal cassette chromosome mec (SCCmec) typing. ... restriction enzymes, was able to predict SCCmec type and MLVA type completely ...
JOURNAL OF CLINICAL MICROBIOLOGY, Oct. 2009, p. 3313–3322 0095-1137/09/$08.00⫹0 doi:10.1128/JCM.00910-09 Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Vol. 47, No. 10

Comparison of Fingerprinting Methods for Typing Methicillin-Resistant Staphylococcus aureus Sequence Type 398䌤 G. Rasschaert,1* W. Vanderhaeghen,2 I. Dewaele,1 N. Janezˇ,1 X. Huijsdens,3 P. Butaye,2,4 and M. Heyndrickx1 Institute for Agricultural and Fisheries Research (ILVO), Technology and Food Science Unit, Brusselsesteenweg 370, 9090 Melle, Belgium1; Veterinary and Agrochemical Research Centre, Groeselenberg 99, 1180 Brussels, Belgium2; National Institute of Public Health and the Environment (RIVM), P.O. Box 1, 3720 BA Bilthoven, The Netherlands3; and Ghent University, Faculty of Veterinary Medicine, Department of Bacteriology, Pathology and Poultry Diseases, Salisburylaan 133, 9820 Merelbeke, Belgium4 Received 8 May 2009/Returned for modification 23 July 2009/Accepted 17 August 2009

This study evaluates the multiple-locus variable-number tandem-repeat assay (MLVA) and pulsed-field gel electrophoresis (PFGE) when using restriction enzymes BstZI, SacII, and ApaI to fingerprint a diverse collection of methicillin (meticillin)-resistant Staphylococcus aureus (MRSA) sequence type 398 (ST398) isolates. These isolates had been characterized previously by multilocus sequence typing, spa typing, and staphylococcal cassette chromosome mec (SCCmec) typing. Typeability and discriminatory power were analyzed, and the concordance between the various methods was determined. All MRSA ST398 isolates were typeable by the MLVA and PFGE using BstZI, SacII, and ApaI. With each method, the MRSA ST398 isolates formed a separate group from the two non-ST398 MRSA strains. PFGE, performed with any of the three restriction enzymes, had the most discriminatory power, followed by MLVA, spa typing, and SCCmec typing. The MLVA showed the highest concordance with PFGE using ApaI and spa typing. As further expressed by the Wallace coefficient, the MLVA type was poorly predicted by spa typing, whereas the spa type was well predicted by MLVA. PFGE, using a combination of all three restriction enzymes, had the highest concordance with the MLVA but had a low probability of being predicted by MLVA. PFGE, using a combination of all three restriction enzymes, was able to predict SCCmec type and MLVA type completely and had a high probability of predicting spa type. Both the MLVA and PFGE could be used to discriminate among the MRSA ST398 isolates. Although the MLVA is a faster technique, PFGE had more discriminatory power than the MLVA, especially when a combination of restriction enzymes was used. methods with good discriminatory power are necessary to conduct large-scale epidemiological studies. Traditionally, human MRSA isolates have been typed by pulsed-field gel electrophoresis (PFGE), using SmaI as the restriction enzyme (19). The advantages of using PFGE are good discriminatory power and good reproducibility at the interlaboratory level when standardized protocols are used. However, AAMRSA is not typeable by this method, as the activity of SmaI is blocked due to methylation of the restriction site (2). More recently, methods based on DNA sequencing, such as MLST and spa typing, are increasingly being used to discriminate among different MRSA strains. Given their excellent interlaboratory reproducibility, online databases have been made to collate and harmonize data from various geographic regions. The drawback of MLST, which measures sequence variation at seven housekeeping loci, is its limited use with epidemiological studies due to its weak discriminatory power, time-consuming protocols, and high costs. spa typing, based on the variation in repeats present in the X-region of staphylococcal protein A, has a discriminatory power that lies between those of PFGE and MLST. Within ST398, several spa types have been distinguished, although the number of spa types seems rather limited in most countries. One promising method is the multiple-locus variable-number tandem-repeat assay (MLVA), a PCR-based method,

Infections caused by methicillin (meticillin)-resistant Staphylococcus aureus (MRSA) have been a problem in hospitals and nursing homes for many decades. These MRSA isolates are therefore called health care-associated MRSA (HA-MRSA) isolates (1). Since the early 1990s, MRSA has emerged in healthy persons without risk factors for MRSA infections. These isolates are the so-called community-associated MRSA (CA-MRSA) isolates (20). In the last few years, MRSA has been isolated from livestock animals (pigs in particular) and pig farmers (5, 6, 30). These MRSA strains are called animal-associated MRSA (AA-MRSA) strains. It seems that (livestock) animals form a new, separate reservoir. These AA-MRSA strains all appear to belong to the new clonal complex 398 (CC398), with sequence type 398 (ST398) as the basic type, as determined by multilocus sequence typing (MLST) (29). MRSA ST398 has already been isolated in Europe, Asia, and North America (32). Considering the worldwide spread of MRSA, epidemiological questions arise about its transmission within farms, among farms, and from farms to the population. Fast and inexpensive typing

* Corresponding author. Mailing address: Institute for Agricultural and Fisheries Research (ILVO), Technology and Food Science Unit, Brusselsesteenweg 370, 9090 Melle, Belgium. Phone: 32 9 272 30 26. Fax: 32 9 272 30 01. E-mail: [email protected]. 䌤 Published ahead of print on 26 August 2009. 3313

3314

RASSCHAERT ET AL.

J. CLIN. MICROBIOL.

TABLE 1. MRSA ST398 isolates used in the present study grouped according their origin and types defined by the MLVA and PFGE Strain no.a

MB4363 MB4364 MB4366 MB4367 MB4368 MB4387 MB4388 MB4389 MB4390 MB4391 MB4392 MB4394 MB4395 MB4396 MB4397 MB4398 MB4393 MB4360 MB4361 MB4362 MB4365a RIVM3 RIVM10 RIVM21 RIVM26 RIVM29 RIVM30 RIVM33 RIVM34 RIVM38 RIVM39 RIVM40 ATCC 33592 ATCC 43300

Origin

Pig, farm A, Belgium Pig, farm A, Belgium Pig, farm B, Belgium Pig, farm B, Belgium Pig, farm C, Belgium Nose of a pig farmer, farm D, Belgium Nose of a pig farmer, farm E, Belgium Nose of a pig farmer, farm F, Belgium Nose of a staff member, pig farm G, Belgium Environment of pig farm G, Belgium Environment of pig farm H, Belgium Pig, farm I, Belgium Pig, farm I, Belgium Nose of a pig farmer, farm J, Belgium Pig, farm J, Belgium Pig, farm J, Belgium Nose of a veterinarian, Belgium Cow with subclinical mastitis, farm K, Belgium Cow with mastitis, farm L, Belgium Cow with mastitis, farm M, Belgium Chicken, farm N, Belgium Nose of a pig farmer, farm O, The Netherlands Wound/abscess of a pig farmer, farm P, The Netherlands Pig, farm Q, The Netherlands Pig, farm R, The Netherlands Pig, farm S, The Netherlands Pig, farm T, The Netherlands Pig, farm U, The Netherlands Pig, farm V, The Netherlands Pig, farm W, The Netherlands Pig, farm W, The Netherlands Pig, farm X, The Netherlands Human MRSA reference strain Human MRSA reference strain

PFGE type using:

Spa type

SCCmec type

MLVA type

BstZI

SacII

ApaI

t034 t034 t011 t011 t011 t034 t011 t011 t011 t011 t011 t011 t011 t2370 t011 t011 t567 t011 t011 t011 t1456 t108 t567

IVa IVa IVa IVa IVa V IVa V V V IVa V V V IVa IVa V V IVa IVa IVa V V

A A B B B C B C C C B C C D B B E F B B G H I

BstZI-1 BstZI-1 BstZI-2 BstZI-2 BstZI-3 BstZI-4 BstZI-2 BstZI-5 BstZI-5 BstZI-6 BstZI-1 BstZI-5 BstZI-5 BstZI-6 BstZI-7 BstZI-7 BstZI-8 BstZI-9 BstZI-3 BstZI-2 BstZI-10 BstZI-11 BstZI-12

Sac-1 Sac-1 Sac-2 Sac-2 Sac-3 Sac-4 Sac-2 Sac-5 Sac-5 Sac-6 Sac-1 Sac-5 Sac-5 Sac-7 Sac-8 Sac-8 Sac-9 Sac-10 Sac-3 Sac-2 Sac-11 Sac-9 Sac-12

Apa-1 Apa-1 Apa-2 Apa-2 Apa-3 Apa-4 Apa-5 Apa-6 Apa-7 Apa-8 Apa-2 Apa-6 Apa-6 Apa-6 Apa-2 Apa-2 Apa-9 Apa-10 Apa-11 Apa-2 Apa-2 Apa-7 Apa-12

t108 t108 t1254 t1254 t108 t108 t011 t011 t567 t037 t007

V V IVa IVa V V IVa IVa V

H H B B H H B B G J K

BstZI-8 BstZI-13 BstZI-1 BstZI-1 BstZI-8 BstZI-8 BstZI-1 BstZI-1 BstZI-8 BstZI-14 BstZI-15

Sac-9 Sac-13 Sac-1 Sac-1 Sac-9 Sac-9 Sac-1 Sac-1 Sac-9 Sac-14 Sac-15

Apa-7 Apa-13 Apa-2 Apa-2 Apa-7 Apa-7 Apa-2 Apa-2 Apa-12 Apa-14 Apa-15

a MB, Molecular Bacteriology strain collection at ILVO-T&V laboratories (MB4365 originates from the work of Persoons et al. 关22兴; MB4387, MB4388, and MB4389 were provided by the Department of Microbiology, Université Libre de Bruxelles, Ho ˆpital Erasme, Belgium); RIVM, National Institute for Public Health and the Environment; ATCC, American Type Culture Collection.

based on the analysis of the number of repeats in the variablenumber tandem-repeat regions of various individual genes. This method has proven to be useful for typing both Staphylococcus aureus and clinical MRSA isolates with good reproducibility and good discriminatory power. Because the MLVA is also simple, inexpensive, and easy to interpret, it is useful as a typing method for large-scale epidemiological studies (10, 11, 12, 15, 16, 17, 23, 24, 27). This study aimed to investigate various methods for typing MRSA ST398 isolates. An MLVA, consisting of a selection of primers from three existing MLVA systems, was tested with a collection of MRSA ST398 isolates. In addition, PFGE with restriction enzymes other than SmaI was performed with this set of isolates. These isolates had been previously characterized by MLST, spa typing, and staphylococcal cassette chromosome mec (SCCmec) typing. Typeability and discriminatory power were analyzed for all methods, and the concordance among the different methods was determined. MATERIALS AND METHODS Selection of the strains. This study used 34 MRSA strains (Table 1). MLST typing revealed that 32 isolates belong to MRSA ST398 (8; http://saureus.mlst

.net). The isolates had various spa types (www.ridom.de/staphtype), SCCmec types, and origins (18, 21, 33). The strains were isolated from humans who were in contact with pigs, chickens, cows, and the farm environment. Twenty-one strains were isolated in Belgium during studies performed in 2007 and 2008 (7, 22, 31). Another 11 strains originated in The Netherlands (5, 28). The strains ATCC 33592 and ATCC 43300 were included as non-ST398 MRSA reference strains. MLVA. Chromosomal DNA was prepared by suspending a few colonies from a tryptic soy agar plate in water containing 10% lysostaphin (1 mg/ml). After incubation for 10 min at 37°C and centrifugation for 2 min at 14,000 ⫻ g, DNA was prepared further according the protocol of Flamm et al. (9). First, the primer pairs listed in Table 2 were used separately in order to evaluate the typeability and discriminating power of each primer pair on the collection of 34 MRSA isolates. Each PCR mixture contained 1⫻ PCR buffer, 1.3 mM MgCl2, 2.5 U AmpliTaq DNA polymerase (Applied Biosystems, Foster City, CA), 200 ␮M of each deoxynucleoside triphosphate, 1.5 ␮M of a primer pair, and 2 ␮l of template DNA (25 ng/␮l) in a final volume of 25 ␮l. PCR was performed with predenaturation at 94°C for 5 min, followed by 25 cycles of 30 s at 94°C, 30 s at 58°C, and 60 s at 72°C, with a final extension at 72°C for 5 min on a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA). Ten microliters of PCR product was used for electrophoresis in 2% SeaKem LE agarose with 1⫻ Tris-borate-EDTA (TBE; pH 8) for 60 min at 100 V. Only primers that discriminated among the 32 MRSA ST398 isolates were included in the final MRSA ST398 MLVA. This multiplex PCR was performed under the same conditions described above. Eighteen microliters of PCR product was electrophoresed in 1.5% SeaKem LE agarose with 1⫻ TBE (pH 8) for

VOL. 47, 2009

MRSA ST398 FINGERPRINTING METHODS

3315

TABLE 2. Primers tested for discriminatory power with the MLVA Locus

clfA clfB spa ssp sdr sdr fnbAB cna SIRU01 SIRU05 SIRU07 SIRU13 SIRU15 SIRU16 SIRU21 a

Primer

Primer sequence (5⬘–3⬘)

Reference

No. of types

Discriminatory power (DI)

clfA-F clfA-R clfB-F clfB-R spa-F spa-R sspA-F sspA-R sdrCDE-F sdrCDE-R sdrCDE2-F sdrCDE2-R fnBP-F fnBP-R cna-F cna-R SIRU01-F SIRU01-R SIRU05-F SIRU05-R SIRU07-F SIRU07-R SIRU13-F SIRU13-R SIRU15-F SIRU15-R SIRU16-F SIRU16-R SIRU21-F SIRU21-R

GATTCTGACCCAGGTTCAGA CTGTATCTGGTAATGGTTCTTT ATGGTGATTCAGCAGTAAATCC CATTATTTGGTGGTGTAACTCTT AGCACCAAAAGAGGAAGACAA GTTTAACGACATGTACTCCGT ATCMATTTYGCMAAYGATGACCA TTGTCTGAATTATTGTTATCGCC GTAACAATTACGGATCATGATG TACCTGTTTCTGGTAATGCTTT TTACGGATCATGATGATTTCA CAYTACCTGTTTCTGGTAATGCTT GGTCAAGCRCAAGGAACCART AATAATCCGCCGAACAACAT AAAATGACAAAATGGCAAG CAGGTTTAGTTGGTGGTGTT CATGAGCAGTGCCTCCTTTA CGCCTTGTATCTTAAATTTGTTG AAAGCGTAAACACTCCGAACA CATCGGGTTCACTTCTTTCTT CGCAAAATTTATTTCGTCGT TTGGCATCATTTAGCAGACG GGGGGAAATTCTAAGCAACC CGAAATTTTCCACGTCGATT AATTATAGCATATTAGAGCCCCTTA ACGTAAAGGTCGCGACAAAA TCAGCTTGTGACTCTTCTTTCG AAGAACGTGTTGAAGGATTTTT AGCAGTAGTGCCGTTTGCTT GCTCAAGCACCAAAAGAGGA

23

3

0.65

23

5

0.49

23

6

0.59

23

1

0

23

—a



10

5

0.50

10

1

0

10

1

0

11

1

0

11

0

0

11

1

0

11

1

0

11

1

0

11

0

0

11

6

0.59

—, for this primer pair, nonspecific bands were obtained.

240 min at 140 V. The DNA size standard was a 100-bp DNA ladder (Invitrogen, Carlsbad, CA). The gels were stained in ethidium bromide, visualized on an UV transilluminator, and digitally photographed. Gel patterns were analyzed with BioNumerics version 5.10 (Applied Maths, Sint-Martens-Latem, Belgium) and clustered using Dice’s coefficient (tolerance, 1%; optimization, 0.1% to 1%) and the unweighted-pair group method using average linkages (UPGMA). PFGE. The plugs were prepared according the protocol of Struelens et al. (26), with modifications. Briefly, an overnight-grown culture was suspended in EET buffer (100 mM EDTA, 10 mM EGTA, 10 mM Tris-HCl [pH 8]) and adjusted to an optical density at 600 nm value of 0.9. The suspension was mixed with equal volumes of a 2% solution of low-melting-temperature agarose (InCert agarose; Cambrex Bio Science Rockland, Rockland, MA) in EET buffer. After solidification, the agarose plugs were incubated for 4 h at 37°C in EET buffer containing 1 mg lysozyme and 50 ␮g lysostaphin per ml. The plugs were transferred into EET buffer containing 1% sodium dodecyl sulfate and 1 mg proteinase K per ml of buffer and incubated overnight at 50°C. Afterwards, the plugs were washed thoroughly with TE buffer (10 mM Tris-HCl, 1 mM EDTA [pH 8]). The restriction enzyme chosen was isoschizomer XmaI, because it is a methylation-insensitive version of SmaI. XmaI cuts the same recognition sequence as SmaI, but at a different position. The other restriction enzymes were chosen based on the restriction digest and PFGE program provided at http://insilico.ehu .es/ (3). DNA was digested overnight with 40 U SmaI (Invitrogen, Carlsbad, CA), SacII, XmaI (New England Biolabs, Ipswich, MA), ApaI, or BstZI (Promega, Madison, WI) at the recommended assay temperatures. DNA fragments were separated on a CHEF Mapper system (Bio-Rad Laboratories, Hemel Hempstead, United Kingdom) in a 1% SeaKem gold agarose gel (Lonza, Rockland, MA). The running conditions were 6 V/cm in 0.5⫻ TBE (45 mM Tris, 45 mM boric acid, 1 mM EDTA [pH 8]) at 14°C. For separation of SmaI and XmaI DNA fragments, runs lasted 18 h, with switch times from 2.63 s to 63.8 s. ApaI fragments were separated during runs of 16 h, with switch times of 0.22 s to 17.33 s. DNA fragments digested with BstZI and SacII were separated during a run of 20 h, with an initial switch time of 0.46 s and a final switch time of 35.38 s for BstZI and a switch time of 44.69 s for SacII. Gels were stained with ethidium bromide, destained in water, and digitally captured under UV light. Gel images

were visually analyzed with BioNumerics version 5.10 (Applied Maths, SintMartens-Latem, Belgium) using Salmonella Braenderup H9812 digested with XbaI as a normalization reference. The similarity between the fingerprints was calculated using the Dice coefficient (with an optimization of 1% and a position tolerance of 0.5% to 1%), and they were grouped together according to their similarities using UPGMA. Delineation of MLVA types and pulsotypes. An MLVA type was assigned on the basis of a difference of at least one band. As this method relies on the number of the tandem repeats in the variable-number tandem-repeat regions of different individual genes, even a small shift in the position of one band was considered to show another MLVA type. MLVA types are indicated by capital letters, as shown below. For each restriction enzyme, a corresponding pulsotype was assigned, which was based on the difference in the presence or absence of at least one band. This was indicated by the name of the restriction enzyme followed by a number (e.g., BstZI-1). Calculation of discriminatory power and concordance. The discriminatory power (D) was calculated according to the formula described by Hunter and Gaston (13, 14). The concordance of different typing techniques was calculated with the software described by Carric¸o et al. (4). Rand’s index is an index that shows the proportion of agreement, whereas the adjusted Rand’s index (R) shows the proportion of agreement corrected for the presence of chance agreement. Wallace’s coefficient (W) indicates the probability that two strains classified as the same type by one method will also be classified as the same one when using the other method. For these calculations, isolates originating from the same farm and having the same spa type, SCCmec type, MLVA type, and pulsotype were considered to be the same strain and were excluded from further analysis.

RESULTS Optimizing the MLVA. All primer pairs generated one band for at least one of the two reference strains, except the primers

3316

RASSCHAERT ET AL.

J. CLIN. MICROBIOL.

FIG. 1. Dendrogram and fingerprints for the MLVA with 32 MRSA ST398 isolates and two reference MRSA strains. The similarities between the fingerprints were calculated using the Dice coefficient (optimization, 1%; position tolerance, 0.1% to 1%), and the fingerprints were grouped according to their similarities by use of the UPGMA algorithm.

for sdr, which generated multiple bands. Data for the number of types for the 32 MRSA ST398 isolates and the discriminatory power of each primer pair are given in Table 2. The primer pairs SIRU5 and SIRU16 generated no amplicons for the MRSA ST398 isolates. Primer pairs sspA, fnbAB, cna, SIRU1, SIRU7, SIRU13, and SIRU15 did not discriminate among the MRSA ST398 isolates. The ClfA and ClfB primer pairs generated three and five different amplicons, respectively, although the difference in size for clfA was minimal. Both primer pairs spa and SIRU21 had good discriminating power: both found six types among the tested MRSA ST398 isolates. For each isolate, primer pair SIRU21 consistently generated a band that was approximately 50 bp bigger than the bands generated by the spa primers. The locus of SIRU21 has also been identified as spa (12). Last to be evaluated were the two primer pairs sdrCDE and sdrCDE2, used to detect the poly-

morphism of three individual genes, namely sdrC, sdrD, and sdrE. Staphylococcus aureus strains do not always possess all three of these genes (25). The 32 MRSA ST398 strains possess only the sdrC and sdrE genes, as tested in this study according to the protocol of Sabat et al. (25) (data not shown). Given this information, only two amplicons for each isolate can be expected when using these two primer pairs. Only the sdrCDE2 primer pair generated two bands for each MRSA ST398 isolate, and five different types were obtained. The following primer pairs were chosen for their discriminatory power and included in the ST398 MLVA: SIRU21, ClfA, ClfB, and sdrCDE2. The best results were obtained by running the sdrCDE2 primer set separately from the other three primer sets and mixing 10 ␮l of the sdrCDE2 PCR product with 8 ␮l of the multiplex PCR product before electrophoresis.

VOL. 47, 2009

MRSA ST398 FINGERPRINTING METHODS

3317

FIG. 2. Dendrogram and fingerprints for PFGE using restriction enzyme BstZI with 32 MRSA ST398 isolates and two reference MRSA strains. The similarities among the fingerprints were calculated using the Dice coefficient (optimization, 1%; position tolerance, 0.5% to 1%), and the fingerprints were grouped according to their similarities by use of the UPGMA algorithm. The vertical gray line shows the delineation level of 97%.

Analysis of the MRSA ST398 strains using the MLVA. All 34 isolates were typeable using the MLVA described above. The MRSA ST398 isolates clustered together with a similarity of 33%. As some differences were so small, it was not possible to determine a cutoff value in order to automatically categorize the types with the BioNumerics software. For example, MB4393 and RIVM10 were 100% related according to the BioNumerics software, but it was possible to visually distinguish a small shift due to a smaller clfB band. Using visual categorization, 10 different MLVA types were obtained within the ST398 isolates (Fig. 1). Six MLVA types consisted of only one isolate. Two isolates belonged to MLVA type A, but these two isolates originated from the same farm. MLVA type H consisted of five isolates, all isolated in The Netherlands, whereas MLVA types B and C consisted of 13 and 6 isolates, respectively. These two last types were very similar; the only difference was a small shift of the clfA band. Analysis of the MRSA ST398 isolates using PFGE. None of the AA-MRSA isolates was typeable by SmaI digestion, which confirmed this characteristic as being typical for MRSA CC398. However, they were all typeable using the restriction

enzymes BstZI, SacII, and ApaI. XmaI, an alternative to SmaI, as mentioned above, was also evaluated. However, the fingerprints obtained yielded several weak bands, even after optimization of the restriction conditions. Because of this, XmaI was excluded from the study. A cutoff value of 97% for delineation of the different pulsotypes was determined for each restriction enzyme, according to the criteria of the delineation of pulsotypes as described above in Materials and Methods. The fingerprints of the MRSA ST398 isolates obtained from PFGE using restriction enzyme BstZI clustered together with a similarity of 79%, whereas the two reference non-ST398 MRSA isolates had a similarity with the MRSA ST398 isolates of less than 57%. Based on the cutoff value of 97%, 13 BstZI pulsotypes were determined within the MRSA ST398 isolates (Fig. 2). Using SacII, the AA-MRSA isolates clustered together with a similarity of 76%, whereas the two reference strains were an outgroup, with a similarity of 61% with the MRSA ST398 isolates (Fig. 3). Thirteen SacII pulsotypes were distinguished among the AA-MRSA isolates, with a cutoff value of 97%. The

3318

RASSCHAERT ET AL.

J. CLIN. MICROBIOL.

FIG. 3. Dendrogram and fingerprints for PFGE using restriction enzyme SacII with 32 MRSA ST398 isolates and two reference MRSA strains. The similarities among the fingerprints were calculated using the Dice coefficient (optimization, 1%; position tolerance, 0.5% to 1%), and the fingerprints were grouped according to their similarities by use of the UPGMA algorithm. The vertical gray line shows the delineation level of 97%.

ST398 cluster could be further subdivided into two major groups (cutoff value of 80%), one group of isolates with SCCmec type V, and another group of isolates with SCCmec type IVa. The same isolates grouped together with the restriction enzyme SacII as those that grouped with restriction enzyme BstZI, with the exception of three pulsotypes. Pulsotype Sac-9 contained pulsotypes BstZI-8 and BstZI-11; thus, in this case, BstZI was more discriminating. In contrast, restriction enzyme SacII was more discriminating for pulsotypes Sac-6 and Sac-7, which shared the same fingerprints when using restriction enzyme BstZI. The AA-MRSA isolates clustered together with a similarity of 71% when using restriction enzyme ApaI (Fig. 4). Again, the two reference strains had a similarity of 62% or less with the MRSA ST398 cluster. As with restriction enzyme SacII, the ST398 cluster could be further subdivided into two groups (cutoff value of 80%), separating isolates carrying SCCmec type IVa and those carrying SCCmec type V. With the cutoff

value of 97%, 13 different Apa pulsotypes were distinguished. However, there was only little correlation between the pulsotypes obtained by BstZI or SacII and the pulsotypes obtained by ApaI. For example, pulsotype Apa-2 contained four different BstZI/SacII pulsotypes. On the other hand, isolates belonging to pulsotype BstZI-8 or Sac-9 were further divided into three pulsotypes using ApaI. Discriminatory power and concordance. The discriminatory indices (DI) of the methods, and those of the combination of the methods, are shown in Table 3. The discriminatory power of SCCmec typing, spa typing, and the MLVA for ST398 were quite low, as follows: 0.51, 0.74, and 0.81, respectively. The combination of spa typing and the MLVA resulted in more discriminatory power (D ⫽ 0.87) than that of the MLVA used alone but was still lower than that of PFGE when using a single restriction enzyme (D ⫽ 0.88 to 0.92). For this data set, the DI of PFGE with the restriction enzymes BstZI and ApaI was as high as the DI of PFGE using all three restriction enzymes

VOL. 47, 2009

MRSA ST398 FINGERPRINTING METHODS

3319

FIG. 4. Dendrogram and fingerprints for PFGE using restriction enzyme ApaI with 32 MRSA ST398 isolates and two reference MRSA strains. The similarities among the fingerprints were calculated using the Dice coefficient (optimization, 1%; position tolerance, 0.5% to 1%), and the fingerprints were grouped according to their similarities by use of the UPGMA algorithm. The vertical gray line shows the delineation level of 97%.

together or using those enzymes in combination with the MLVA (D ⫽ 0.97). The combination of PFGE (when using all restriction enzymes or BstZI and ApaI) and spa typing had the most discriminatory power (D ⫽ 0.98). The following isolates were indistinguishable by all means applied: MB4392 and RIVM38; MB4362 and MB4366; RIVM29 and RIVM30; and RIVM21, RIVM33, and RIVM34. The concordance among different typing methods is shown in Tables 4 and 5. For the MLVA, the highest concordance was with PFGE using ApaI and with spa typing. As further expressed by the Wallace coefficient, the MLVA type was completely predicted by PFGE using a combination of all the restriction enzymes (W ⫽ 1). The MLVA type was poorly predicted only by spa typing (W ⫽ 0.48), whereas the spa type was well predicted by MLVA (W ⫽ 0.68). As already indicated above, the concordance between PFGE using restriction enzyme BstZI and PFGE using SacII was very high (Rand’s index ⫽ 0.98; R ⫽ 0.89). PFGE, using

a combination of all three restriction enzymes, was able to predict the SCCmec type and MLVA type completely (W ⫽ 1.00) and had a high probability of predicting the spa type (W ⫽ 0.64). SCCmec type was a poor predictor for all typing methods (W ⫽ 0.06 to 0.38), although it was well predicted by MLVA (W ⫽ 0.98) and completely predicted using all PFGE methods (W ⫽ 1.00). DISCUSSION AA-MRSA strains of CC398 are not typeable using PFGE with SmaI digestion, the gold standard for MRSA typing. For this reason, until now, sequence-based typing methods, such as MLST and spa typing, have been used for typing AA-MRSA in epidemiological studies. Although these methods have the advantage of reliability and excellent interlaboratory reproducibility, their discriminatory power is weak and, consequently, less useful for short-term epidemiological studies such as an

3320

RASSCHAERT ET AL.

J. CLIN. MICROBIOL.

restriction enzymes had higher discrimination than the use of only one restriction enzyme. The discriminatory power of PFGE using only BstZI and ApaI was equally as good as using all three restriction enzymes. The discriminatory power of the MLVA, the second best typing method, was only a bit higher than that of spa typing, although spa typing was also able to discriminate among identical MLVA types in a few cases and within a pulsotype on one occasion (using all three restriction enzymes or BstZI and ApaI). In this set of AA-MRSA isolates, the MLVA was not able to distinguish among isolates with the same PFGE profile. In agreement with our results, previously published reports of MLVA typing of clinical MRSA isolates had more discriminatory power than those of spa typing and MLST (12, 15, 16, 17). In contrast, most studies comparing the MLVA with PFGE have shown that the discriminatory power of PFGE was comparable to that of the MLVA or only slightly better. Luczak-Kadlubowska et al. (16) and Sabat et al. (23) demonstrated that their MLVA system had a comparable discriminatory power to that of PFGE. The same group had reported previously that PFGE was superior to the MLVA. These findings were based on the setup of the study and the selection of the MRSA isolates (17). The results of the MLVA based on staphylococcal interspersed repeat units (11) depended on the set of isolates: for some isolates, the MLVA was able to distinguish within the same PFGE profile, but for others, the MLVA was less discriminatory than PFGE. It should be kept in mind that in the publications mentioned above, the MLVA consisted of a different combination of primer sets than the MLVA described here. Furthermore, in those publications, PFGE was performed using only the restriction enzyme SmaI. The MLVA and PFGE were also used on a wide range of MLST and spa types belonging to HA- and CA-MRSA but not to MRSA ST398. Those studies aimed to investigate whether the faster MLVA technique was able to predict PFGE typing, MLST, or spa typing. In order to do so, relatively low similarity cutoff values between 70% and 80% were used for the MLVA and PFGE (17, 27). This contrasts with the present study, during which the MLVA and PFGE techniques were evaluated for their power to distinguish among MRSA ST398 isolates. To this end, isolates that differed even slightly were considered to belong to different types. For the MLVA, some differences were small. It is thus advisable to use capillary electrophoresis for better resolution between bands.

TABLE 3. Discriminatory power of the methods, apart and in combination, evaluated with 27 MRSA ST398 strains No. of No. of No. of Cluster unique clustered clusters size isolates isolates

Methoda

No. of types

SCCmec typing spa typing MLVA spa typing ⫹ MLVA PFGE (BstZI) PFGE (SacII) PFGE (ApaI) PFGE (BstZI ⫹ SacII) PFGE (BstZI ⫹ ApaI) PFGE (RE) PFGE (RE) ⫹ MLVA PFGE (RE) ⫹ spa typing

2 7 9 11 13 13 13 14

0 2 5 7 7 8 9 9

27 25 22 20 20 19 18 18

2 5 4 4 6 5 4 5

12–15 2–13 2–10 2–8 2–5 2–6 2–8 2–5

0.51 0.74 0.81 0.87 0.92 0.91 0.88 0.92

20

16

11

4

2–4

0.97

20 20

16 16

11 11

4 4

2–4 2–4

0.97 0.97

21

16

11

5

2–3

0.98

a

DI

RE, combination of all three restriction enzymes.

examination of contamination routes on farms. To our knowledge, this is the first study exploring alternative fingerprinting techniques for subtyping MRSA ST398. In this study, PFGE using the restriction enzymes XmaI (an isoschizomer of SmaI), BstZI, SacII, and ApaI was tested to obtain the enzymes’ discriminatory power with a set of 32 MRSA ST398 isolates. In addition, three recently developed MLVA techniques, successfully applied to Staphylococcus aureus and more specifically to HA-MRSA and CA-MRSA, were evaluated and optimized for use with ST398 strains. Both the initial Polish system (15, 16, 17, 23, 24) and the Swiss system (10) are based on tandem repeats in genes, whereas the British system is based on staphylococcal interspersed repeat units, called SIRUs (11, 12). In this study, only primers with discriminatory power were included in this optimized MLVA specific for AA-MRSA. All isolates were typeable using both the MLVA and PFGE. However, fingerprints obtained by PFGE using XmaI had weak patterns and were thus excluded from the study. The discriminatory power of PFGE, when using any of the remaining three restriction enzymes, was higher than that of SCCmec typing, spa typing, the MLVA, or spa typing in combination with the MLVA. Within PFGE, the combination of all three

TABLE 4. Concordance among typing methods according to Rand’s and adjusted Rand’s indices Rand’s and adjusted Rand’s index values obtained using the indicated typing methods spa Typing method

SCCmec MLVA PFGE using PFGE using PFGE using PFGE using a

BstZI SacII ApaI REa

SCCmec

Rand’s index

Adjusted Rand’s index

0.54 0.80 0.73 0.74 0.72 0.74

0.07 0.44 0.11 0.14 0.14 0.08

Rand’s index

Adjusted Rand’s index

0.70 0.59 0.60 0.63 0.54

0.38 0.17 0.19 0.25 0.07

RE, combination of all three restriction enzymes.

MLVA Rand’s index

Adjusted Rand’s index

0.82 0.83 0.85 0.84

0.26 0.30 0.44 0.25

PFGE using BstZI Rand’s index

Adjusted Rand’s index

0.98 0.86 0.95

0.89 0.24 0.54

PFGE using SacII Rand’s index

Adjusted Rand’s index

0.87 0.94

0.31 0.49

PFGE using ApaI Rand’s index

Adjusted Rand’s index

0.91

0.38

VOL. 47, 2009

MRSA ST398 FINGERPRINTING METHODS

3321

TABLE 5. Wallace’s coefficient values among typing methods Wallace’s coefficient value obtained using the indicated typing method Typing method spa

spa SCCmec MLVA PFGE using PFGE using PFGE using PFGE using a

BstZI SacII ApaI REa

SCCmec

MLVA

0.56

0.48 0.38

0.30 0.68 0.46 0.50 0.45 0.64

0.98 1.00 1.00 1.00 1.00

0.57 0.59 0.67 1.00

PFGE using BstZI

PFGE using SacII

PFGE using ApaI

PFGE using REa

0.14 0.16 0.24

0.17 0.19 0.29 0.96

0.2 0.25 0.42 0.39 0.44

0.08 0.06 0.17 0.39 0.34 0.26

0.84 0.26 1.00

0.33 1.00

1.00

RE, combination of all three restriction enzymes.

Even when combining all methods, some isolates from different origins were indistinguishable from each other: pig isolates MB4392 and RIVM38 (one isolated in Belgium, the other in The Netherlands), MB4362 and MB4366 (the first isolated from a cow, the second from a pig), RIVM29 and RIVM30 (both isolated from Dutch pigs), and finally, RIVM21, RIVM33, and RIVM34 (also found in Dutch pig herds). It may be possible to distinguish between these isolates using other techniques. Another possibility, given the large number of animals transported within and between Belgium and The Netherlands, is that those isolates are the same strain circulating in different countries and in different animal species. The concordance among the different methods was calculated according to Rand’s and adjusted Rand’s indices and Wallace coefficients. As this is calculated for only a limited set of data, it should be interpreted carefully. According to the obtained results, the spa type was best predicted by MLVA, whereas the MLVA type was poorly predicted only by spa typing. The MLVA type was completely predicted by PFGE using a combination of all restriction enzymes (W ⫽ 1) but not vice versa. Tenover et al. (27) also concluded that the MLVA cannot be used to predict PFGE types. In conclusion, both the MLVA and PFGE can be used to discriminate among isolates belonging to the newly emerging MRSA ST398 type. All isolates were typeable using both of the fingerprinting techniques and had better discriminatory power than the commonly used techniques. Although the MLVA is a faster technique, PFGE had more discriminatory power than the MLVA, especially when two restriction enzymes (i.e., BstZI and ApaI) were used. Both techniques can be used for epidemiological tracing of the MRSA ST398 isolates. ACKNOWLEDGMENTS We are grateful to M. Struelens from the Laboratoire de Re´fe´rence MRSA—Staphylocoques, Department of Microbiology, Universite´ Libre de Bruxelles, Ho ˆpital Erasme, Brussels, Belgium, for providing strains MB4387, MB4388, and MB4389. We thank Miriam Levenson for critical reading of the manuscript. REFERENCES 1. Barrett, F. F., R. F. McGehee, Jr., and M. Finland. 1968. Methicillin-resistant Staphylococcus aureus at Boston City Hospital. Bacteriologic and epidemiologic observations. N. Engl. J. Med. 279:441–448. 2. Bens, C. P. M., A. Voss, and C. H. W. Klaassen. 2006. Presence of a novel DNA methylation enzyme in methicillin-resistant Staphylococcus aureus isolates associated with pig farming leads to uninterpretable results in standard pulsed-field gel electrophoresis analysis. J. Clin. Microbiol. 44:1875–1876. 3. Bikandi, J., R. S. Millan, A. Rementeria, and J. Garaizar. 2004. In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction. Bioinformatics 20:798–799.

4. Carric¸o, J. A., C. Silva-Costa, J. Melo-Cristino, F. R. Pinto, H. de Lencastre, J. S. Almeida, and M. Ramirez. 2006. Illustration of a common framework for relating multiple typing methods by application to macrolide-resistant Streptococcus pyogenes. J. Clin. Microbiol. 44:2524–2532. 5. de Neeling, A. J., M. J. M. van den Broek, E. C. Spalburg, M. G. van Santen-Verheuvel, W. D. C. Dam-Deisz, H. C. Boshuizen, A. W. van de Giessen, E. van Duijkeren, and X. W. Huijsdens. 2007. High prevalence of methicillin resistant Staphylococcus aureus in pigs. Vet. Microbiol. 122:366– 372. 6. Denis, O., C. Suetens, M. Hallin, B. Catry, I. Ramboer, M. Dispas, G. Willems, B. Gordts, P. Butaye, and M. J. Struelens. 2009. Methicillin-resistant Staphylococcus aureus ST398 in swine farm personnel, Belgium. Emerg. Infect. Dis. 15:1098–1101. 7. Dewaele, I., I. De Man, A. Stael, P. Delputte, P. Butaye, G. Vlaemynck, L. Herman, M. Heyndrickx, and G. Rasschaert. 2008. Methicillin-resistant Staphylococcus aureus (MRSA) on Belgian pig farms, abstr. S:8:5, p. 23. Abstr. ASM Conf. Antimicrob. Resist. Zoonotic Bacteria Foodborne Pathog. 8. Enright, M. C., N. P. J. Day, C. E. Davies, S. J. Peacock, and B. G. Spratt. 2000. Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J. Clin. Microbiol. 38:1008–1015. 9. Flamm, R. K., D. J. Hinrichs, and M. F. Thomashow. 1984. Introduction of pAM␤1 into Listeria monocytogenes by conjugation and homology between native L. monocytogenes plasmids. Infect. Immun. 44:157–161. 10. Francois, P., A. Huyghe, Y. Charbonnier, M. Bento, S. Herzig, I. Topopski, B. Fleury, D. Lew, P. Vaudaux, S. Harbarth, W. van Leeuwen, A. van Belkum, D. S. Blanc, D. Pittet, and J. Schrenzel. 2005. Use of an automated multiple-locus, variable-number tandem repeat-based method for rapid and high-throughput genotyping of Staphylococcus aureus isolates. J. Clin. Microbiol. 43:3346–3355. 11. Hardy, K. J., B. A. Oppenheim, S. Gossain, F. Gao, and P. Hawkey. 2006. Use of variations in staphylococcal interspersed repeat units for molecular typing of methicillin-resistant Staphylococcus aureus strains. J. Clin. Microbiol. 44:271–273. 12. Hardy, K. J., D. W. Ussery, B. A. Oppenheim, and P. M. Hawkey. 2004. Distribution and characterization of staphylococcal interspersed repeat units (SIRUs) and potential use for strain differentiation. Microbiology 150:4045– 4052. 13. Hunter, P. 1990. Reproducibility and indices of discriminatory power of microbial typing methods. J. Clin. Microbiol. 28:1903–1905. 14. Hunter, P. R., and M. A. Gaston. 1988. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J. Clin. Microbiol. 26:2465–2466. 15. Karynski, M., A. J. Sabat, J. Empel, and W. Hryniewicz. 2008. Molecular surveillance of methicillin-resistant Staphylococcus aureus by multiple-locus variable number tandem repeat fingerprinting (formerly multiple-locus variable number tandem repeat analysis) and spa-typing in hierarchic approach. Diagn. Microbiol. Infect. Dis. 62:255–262. 16. Luczak-Kadlubowska, A., A. Sabat, A. Tambic-Andrasevic, M. Payerl-Pal, J. Krzyszton-Russjan, and W. Hryniewicz. 2008. Usefulness of multiple-locus VNTR fingerprinting in detection of clonality of community- and hospitalacquired Staphylococcus aureus isolates. Antonie van Leeuwenhoek 94:543– 553. 17. Malachowa, N., A. Sabat, M. Gniadkowski, J. Krzyszton-Russjan, J. Empel, J. Miedzobrodzki, K. Kosowska-Shick, P. C. Appelbaum, and W. Hryniewicz. 2005. Comparison of multiple-locus variable-number tandem-repeat analysis with pulsed-field gel electrophoresis, spa typing, and multilocus sequence typing for clonal characterization of Staphylococcus aureus isolates. J. Clin. Microbiol. 43:3095–3100. 18. Milheiric¸o, C., D. C. Oliveira, and H. Lencastre. 2007. Multiplex PCR strategy for subtyping the staphylococcal cassette chromosome mec type IV

3322

19.

20.

21.

22. 23.

24.

25.

RASSCHAERT ET AL.

in methicillin-resistant Staphylococcus aureus: ‘SCCmec IV multiplex.’ J. Antimicrob. Chemother. 60:42–48. Murchan, S., M. E. Kaufmann, A. Deplano, R. de Ryck, M. Struelens, C. E. Zinn, V. Fussing, S. Salmenlinna, J. Vuopio-Varkila, N. El Solh, C. Cuny, W. Witte, P. T. Tassios, N. Legakis, W. van Leeuwen, A. van Belkum, A. Vindel, I. Laconcha, J. Garaizar, S. Haeggman, B. Olsson-Liljequist, U. Ransjo, G. Coombes, and B. Cookson. 2003. Harmonization of pulsed-field gel electrophoresis protocols for epidemiological typing of strains of methicillin-resistant Staphylococcus aureus: a single approach developed by consensus in 10 European laboratories and its application for tracing the spread of related strains. J. Clin. Microbiol. 41:1574–1585. Naimi, T. S., K. H. LeDell, K. Como-Sabetti, S. M. Borchardt, D. J. Boxrud, J. Etienne, S. K. Johnson, F. Vandenesch, S. Fridkin, C. O’Boyle, R. N. Danila, and R. Lynfield. 2003. Comparison of community- and health careassociated methicillin-resistant Staphylococcus aureus infection. JAMA 290: 2976–2984. Oliveira, D. C., and H. Lencastre. 2002. Multiplex PCR strategy for rapid identification of structural types and variants of the mec element in methicillin-resistant Staphylococcus aureus. Antimicrob. Agents Chemother. 46: 2155–2161. Persoons, D., S. Van Hoorebeke, K. Hermans, P. Butaye, A. de Kruif, F. Haesebrouck, and J. Dewulf. 2009. Emerg. Infect. Dis. 15:452–453. Sabat, A., J. Krzyszton-Russjan, W. Strzalka, R. Filipek., K. Kosowska, W. Hryniewicz, J. Travis, and J. Potempa. 2003. New method for typing Staphylococcus aureus strains: multiple-locus variable-number tandem repeat analysis of polymorphism and genetic relationships of clinical isolates. J. Clin. Microbiol. 41:1801–1804. Sabat, A., N. Malachowa, J. Miedzobrodzki, and W. Hryniewicz. 2006. Comparison of PCR-based methods for typing Staphylococcus aureus isolates. J. Clin. Microbiol. 44:3804–3807. Sabat, A., D. C. Melles, G. Martirosian, H. Grundmann, A. van Belkum, and W. Hryniewicz. 2006. Distribution of the serine-aspartate repeat proteinencoding sdr genes among nasal carriage and invasive Staphylococcus aureus strains. J. Clin. Microbiol. 44:1135–1138.

J. CLIN. MICROBIOL. 26. Struelens, M. J., A. Deplano, C. Godard, N. Maes, and E. Serruys. 1992. Epidemiologic typing and delineation of genetic relatedness of methicillinresistant Staphylococcus aureus by macrorestriction analysis of genomic DNA by using pulsed-field gel electrophoresis. J. Clin. Microbiol. 30:2599– 2605. 27. Tenover, F. C., R. R. Vaughn, L. K. McDougal, G. E. Fosheim, and J. E. McGowan. 2007. Multiple-locus variable-number tandem-repeat assay analysis of methicillin-resistant Staphylococcus aureus strains. J. Clin. Microbiol. 45:2215–2219. 28. van Belkum, A., D. C. Melles, J. K. Peeters, W. B. van Leeuwen, E. van Duijkeren, X. W. Huijsdens, E. Spalburg, A. J. de Neeling, and H. A. Verbrugh on behalf of the Dutch Working Party on Surveillance and Research of MRSA-SOM. 2008. Methicillin-resistant and -susceptible Staphylococcus aureus sequence type 398 in pigs and humans. Emerg. Infect. Dis. 14:479– 483. 29. van Loo, I., X. Huijsdens, E. Tiemersma, A. de Neeling, N. van de SandeBruinsma, D. Beaujean, A. Voss, and J. Kluytmans. 2007. Emergence of methicillin-resistant Staphylococcus aureus of animal origin in humans. Emerg. Infect. Dis. 13:1834–1839. 30. Voss, A., F. Loeffen, J. Bakker, C. Klaassen, and M. Wulf. 2005. Methicillinresistant Staphylococcus aureus in pig farming. Emerg. Infect. Dis. 11:1965– 1966. 31. Willems, G., M. Dispas, O. Denis, M. Hallin, C. Suetens, B. Gordts, M. Struelens, and P. Butaye. 2007. Characterization of MRSA from pigs in Belgium, p. 18. Abstr. 2nd Symp. Antimicrob. Resist. Anim. Environ. (ARAE). 32. Wulf, M. W. H., and A. Voss. 2008. MRSA in livestock animals—an epidemic waiting to happen? Clin. Microbiol. Infect. 14:519–521. 33. Zhang, K., J. A. McClure, S. Elsayed, T. Louie, and J. M. Conly. 2005. Novel multiplex PCR assay for characterization and concomitant subtyping of staphylococcal cassette chromosome mec types I to V in methicillin-resistant Staphylococcus aureus. J. Clin. Microbiol. 43:5026–5033.