Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2005, p. 3095-3100, Vol. 43, No. 7
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.7.3095-3100.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Microbiology, Faculty of Biotechnology, Jagiellonian University, 30-387 Cracow, Poland,1 National Institute of Public Health, 00-725 Warsaw, Poland,2 Department of Pathology, Hershey Medical Center, P.O. Box 850, Hershey, Pennsylvania 170333
Received 10 November 2004/ Returned for modification 7 January 2005/ Accepted 15 March 2005
|
|
|---|
|
|
|---|
Pulsed-field gel electrophoresis (PFGE) has been considered to be the "gold standard" in typing of a variety of bacteria, including S. aureus. Being highly discriminatory, it is an excellent tool to analyze outbreaks and center-to-center strain transmission events (3, 21, 22, 26), but it has also been used successfully in large-scale S. aureus (mostly MRSA) epidemiological investigations (6, 15, 16). Such studies have been further developed by application of multilocus sequence typing (MLST), which measures sequence variation at seven housekeeping loci (8). MLST characterized lower resolution and, as a DNA sequencing-based method, gives full reproducibility of results between laboratories. The recently introduced spa typing method analyzes sequence polymorphism at a single locus, and with its discriminatory power, between those of PFGE and MLST (2, 6, 11, 19, 29), it seems to be useful not only in macroevolution but also in smaller-scale studies (14). Unfortunately, all of these methods are technically demanding, expensive, and time-consuming and, in the case of PFGE, interlaboratory comparisons of results are difficult (16, 17, 18).
Rapid, easy, and relatively inexpensive typing techniques are those that are based on the PCR method. Several of these have been used in S. aureus outbreak analyses; however, all have important limitations, such as low reproducibility or limited discriminatory power (4, 7, 11, 31). It is possible that these shortfalls will be eliminated with new S. aureus (or MRSA)-specific methods, such as triplex PCR for spa and coa genes, the hypervariable region adjacent to mecA (28, 32), or the restriction profile analysis of the repetitive element STAR (20). Recently, multiplex PCR referred to as MLVA (25) has been developed in our laboratory. It analyzes the variation in number of repeats in seven individual genes (sspA, spa, sdrC, sdrD, sdrE, clfA, and clfB) and in the initial study was found to be comparable in discriminatory power with PFGE. The aim of the present study was to determine the congruence between the isolate or clonal groupings recognized by MLVA, and PFGE, spa typing, and MLST and to propose possible MLVA clustering criteria based on comparison of MLVA and PFGE DNA banding patterns.
|
|
|---|
PFGE. Preparation of genomic DNA of the isolates, followed by SmaI (MBI Fermentas, Vilnius, Lithuania) digestion, and separation in a CHEF DR II apparatus (Bio-Rad, Hercules, Calif.) was performed as described previously (5). SmaI PFGE patterns were saved in the TIF format, exported to the Molecular Analyst database (Bio-Rad), and analyzed by using the Dice coefficient-UPGMA (unweighted pair-group method with arithmetic averages). A dendrogram was generated to examine relatedness of PFGE profiles for all study isolates, and cutoff levels of 75 and 92% were applied to this dendrogram (Fig. 1 and 2, respectively). With the 75% cutoff, isolates differing by up to six DNA fragments were clustered together, whereas the 92% cutoff corresponded to the difference of up to three bands within a cluster.
![]() View larger version (38K): [in a new window] |
FIG. 1. PFGE (left) and MLVA (right) dendrograms of the study isolates generated by the UPGMA algorithm. Isolate clusters were delineated with a 75 and 70% similarity cutoff values for PFGE and MLVA, respectively, as indicated by vertical lines. Double arrows connect the corresponding clusters discerned by both methods (different arrow forms are used only for clarity of visualization).
|
![]() View larger version (33K): [in a new window] |
FIG. 2. PFGE (left) and MLVA (right) dendrograms of the study isolates. Isolate clusters were delineated with a 92 and 77% similarity cutoff values for PFGE and MLVA, respectively, as indicated by vertical lines. Double arrows connect the corresponding clusters discerned by both methods (different arrow forms are used only for clarity of visualization).
|
MLVA typing. MLVA typing was performed as described previously (25). Gel images were exported as TIF files for further analysis by using the Molecular Analyst database. Construction of a dendrogram of banding patterns by using UPGMA was performed by using the Dice coefficient.
spa typing. Amplification of the spa gene X region was performed as described previously (27), and amplicons were then sequenced by using an ABI 377 sequencer (Applied Biosystems, Foster City, Calif.). The spa types were determined with the Ridom SpaServer (12). The spa types with identical or similar repeat profiles were grouped into clusters according to the method of Koreen et al. (14). Any two spa types that differed in the number of repeats but contained many identical repeats in common or showed a single deletion of the internal spa sequence were classified into the same cluster.
MLST. MLST was performed according to protocol described by Enright et al. (8). Sequences of each locus were submitted to the Internet database (www.mlst.net), and resulting allelic profiles were assigned to particular sequence types (ST) for each isolate. START software (13) was used to classify different STs into clusters or clonal complexes (CCs) of phylogenetic relationships. Such clusters were composed of two or more isolates of the same STs or STs which differed at a single locus (single-locus variants) or two loci (double-locus variants) (10).
Calculation of concordance. Intermethod concordance was calculated as the maximum percentage proportion of isolates grouped together into unique patterns/profiles or clusters by two methods compared (24).
|
|
|---|
When the cutoff of 92% was applied, nine PFGE clusters comprising 28 isolates were distinguished and designated from P1' to P9' (Fig. 2). Of the clusters identified with the 75% cutoff, only two clusters, P3 and P8, remained unchanged as clusters P2' and P5', respectively. Clusters P2, P5, and P7 were reduced by three, four, and one isolate, respectively, into clusters P1', P3', and P4'. Cluster P9 was split into four smaller onesP6', P7', P8', and P9'and two unique isolates. Finally, all isolates of clusters P1, P4, and P6 were classified as nonrelated by PFGE with the raised cutoff value (Fig. 2).
spa typing. Our analysis yielded 20 spa types among 59 S. aureus isolates, including six new types (t386 to t391) (Fig. 3). With the use of the criteria by Koreen et al. (14), isolates with similar spa repeat profiles were grouped into four clusters designated S1 to S4. Only a single isolate, isolate BN4 (t159), was classified outside any cluster.
![]() View larger version (61K): [in a new window] |
FIG. 3. Comparison of PFGE and MLVA results (cutoffs 75 and 70%, respectively) with those obtained by MLST and spa typing. The order in which isolates are listed from top to bottom is based on the MLST dendrogram (left side of the figure). Gaps within spa profiles indicate possible deletion/insertion events in the spa locus and were introduced in order to optimize the alignment of the profiles. In PFGE and MLST columns only cluster designations were used; the unique PFGE and MLVA types were not included for clarity.
|
MLVA results and their comparison with the PFGE data. (i) Criteria for defining MLVA clusters. MLVA produced 40 different DNA banding patterns among 59 S. aureus isolates. In order to establish criteria for clustering the patterns into similarity groups, the MLVA typing results were compared to those obtained by PFGE. With the PFGE similarity cutoff of 75%, a comparable set of MLVA clusters was observed when the cutoff value between two MLVA patterns was set up at the level of 70% (Fig. 1). Isolates belonging to the same cluster differed by up to six bands, and isolates classified into different clusters differed by more than six bands. A total of eight MLVA clusters that comprised 40 isolates were distinguished. Almost all clusters delineated in the PFGE dendrogram had corresponding clusters in the MLVA dendrogram. PFGE clusters P2, P8, and P9 were identical to the MLVA clusters M4, M5, and M8, respectively (Fig. 1). A very good correlation was also observed between PFGE clusters P3, P5, P6, and P7 and MLVA clusters M6, M1, M7, and M3, respectively. The only difference was that MLVA groups were smaller by one or two isolates than their PFGE counterparts; however, none of the "missing" isolates was clustered by MLVA together with any other isolates (isolate A005a, PFGE cluster P3; 2689/98, PFGE cluster P5; 2690 and MR47, PFGE cluster P6; and H390, PFGE cluster P7). Isolates in PFGE clusters P1 (isolates 3301 and 3497) and P4 (isolates EMRSA-16, 2684/98, and 2956/01) were all classified as unique isolates by MLVA analysis. The only isolates defined as similar by MLVA and nonrelated by PFGE were two isolates of the MLVA cluster M2 (isolates 3521 and 3498) (Fig. 1).
For the nine PFGE clusters defined by the 92% cutoff, a corresponding set of nine MLVA clusters, M1' to M9', was distinguished when the MLVA cutoff level was raised to 77% (Fig. 2). These clusters comprised 39 isolates, and isolates within a single cluster differed from each other by up to four bands. Clusters M1', M2', M3', M4', M6', and M7' were identical to those from the MLVA analysis with the 70% cutoff. The original cluster M5 was reduced by one isolate into M5', and cluster M8 was divided into two smaller ones, M8' and M9'. All but one PFGE clusters correlated with particular MLVA clusters; however, none of these were identical to its MLVA counterpart. PFGE clusters P1', P2', P3', P4', P5', and P9' corresponded to MLVA clusters M4', M6', M1', M3', M5', and M9', respectively, and were usually smaller by one to three isolates than the MLVA ones. Isolates from clusters P7' and P8' were classified together by MLVA into cluster M8'. The only exception was PFGE cluster P6', which had no a corresponding MLVA cluster, and its two isolates were split between MLVA clusters M8' and M9' (Fig. 2). On the other hand, isolates of the MLVA clusters M2' and M7' were defined as nonrelated by PFGE with the 92% cutoff.
We also searched for a level of variability in PFGE patterns that would always allow classifying any two isolates to corresponding clusters by using MLVA and PFGE methods. Any two isolates with either indistinguishable PFGE patterns (the cutoff 92% clusters P1', P7', and P9', and isolates MR24 and N39 of cluster P8') or differing by a single band (cluster P3') were always assigned to a single, separate MLVA cluster. These differed by no more than four bands within the MLVA clusters.
(ii) Comparison of MLVA with spa and MLST clusters. In all cases, S. aureus isolates of a given MLVA cluster, defined by the 70% cutoff value, were grouped together within the same spa cluster and MLST clonal complex or cluster (Fig. 3). Moreover, such isolates were usually indistinguishable from each other by spa typing and MLST. Exceptions included isolates 3121 (cluster M3) and N39 (cluster M8) that differed in spa types from other isolates in their clusters, and isolates MR1003 (cluster M1), and 2260/98, 3028 and 303/00 (cluster M4) that varied with respect to their STs.
Of the four clusters revealed by spa typing, almost each one was split by MLVA into several clusters and singular isolates. The one exception was the spa cluster S4 with 13 isolates, 12 of which corresponded to the MLVA cluster M8 and one isolate (isolate 3483) was sorted as unique. Similar observations were made when the six MLST clusters were compared to MLVA clusters. The CC5 complex, corresponding to the spa cluster S4, matched precisely the MLVA cluster M8 and the single isolate 3483. Isolate 3483 slightly differed from remaining isolates of the group in spa typing and MLST. The only other MLST clusters that overlapped with MLVA (and spa) clusters were those occurring less often among isolates analyzed: cluster ST80 (corresponding to M5 and S3) and CC45 (M6 and the spa subcluster S2'). CC30 was split completely into unique isolates by MLVA, as was the case of CC1 and all other typing methods used in the present study.
Concordance between methods. The concordance values between the typing methods compared are listed in Table 1. The highest levels of correlation were found between the results produced on one hand by MLVA and PFGE (72.9%), and, on the other, by spa typing and MLST (72.9%). Of the remaining combinations of method pairs, a higher correlation level was observed only in the case of MLVA and spa typing (62.7%).
|
View this table: [in a new window] |
TABLE 1. Correlation between four typing methods for S. aureus
|
|
|
|---|
In general, typing techniques with high discriminatory power have a better level of concordance between themselves than with those of lower resolution potential. By analogy, methods with low discriminatory power correlate better with each other than with highly discriminatory techniques. According to this rule we have observed a high level of congruence between PFGE and MLVA on one hand and spa typing and MLST on the other. Comparability of the resolution potential between MLVA and PFGE has already been demonstrated in our earlier study (25). What is interesting, however, is that the level of concordance between typing techniques can vary and depends on the collection of isolates used for investigation. With a highly diverse collection of isolates, one may observe a good correlation even between typing techniques which differ greatly in discriminatory power. In this way we can explain the difference in the concordance level between PFGE and MLST obtained in the present study (32.2%) and in the analysis performed by Grundmann et al. (67%) (11). In these two studies two different types of S. aureus collections were used. Whereas our collection included only nosocomial and mostly related MRSA isolates, Grundmann et al. (11) studied a collection of mainly MSSA isolates circulating in the community.
The concordance level between MLVA and PFGE was remarkably high. The comparison of the results obtained by both methods was performed at two similarity cutoff levels in the PFGE analysis: 75 and 92%. With the PFGE cutoff of 75%, which referred to the interpretation criteria by Tenover et al. (30), isolates that differed by no more than six bands in MLVA (MLVA cutoff of 70%) could be regarded as related in PFGE and therefore classified into the same MLVA cluster. MLVA differentiated all isolates of different PFGE clusters, and almost all clusters distinguished with one of the methods had their counterparts in the second analysis. Most differences observed in PFGE and MLVA clustering were not accidental but could be explained by other observations. Sporadic isolates of PFGE clusters P1, P3, P4, P5, P6, and P7 that were segregated as unique by MLVA differed slightly from each other or from remaining isolates of their clusters also in spa types (isolates 2689/98, 2690, and A005a), STs (isolates 3301 and 3497), or both (isolates 2684/98 and EMRSA-16). On the other hand, two isolates of the MLVA cluster M2 (isolates 3521 and 3498), which were separated from each other by PFGE, turned out to be the same spa type and ST. The PFGE interpretation criteria by Tenover et al. (30) allow for identifying broader isolate clusters that have been useful in larger-scale epidemiological or phylogenetic studies (6, 15, 16). However, short-term epidemiological studies aimed mostly at outbreak investigations usually require more stringent criteria, especially in the case of highly clonal organisms, such as S. aureus.
Using the PFGE cutoff of 92%, which reflected up to a three-band difference between two PFGE patterns, it was possible to raise the MLVA cutoff level to 77% in order to obtain the most stringent criteria in MLVA interpretation. Such MLVA clusters grouped DNA patterns that differed by up to four bands from each other. Under these conditions, almost each of the PFGE clusters had a corresponding MLVA cluster, with the only exception of the PFGE cluster P6', the two isolates of which were split into MLVA clusters M8' and M9'. With the MLVA cutoff of 70%, these two clusters were combined into a single, broad M8, and it is possible that with more diverse isolates of this phylogenetic lineage included in the analysis, clusters M8' and M9' would form a single one also under the stringent conditions. The even better results were obtained when the MLVA cutoff 77% clusters were assigned to PFGE clusters that grouped isolates of up to one band difference only. In this case, all PFGE clusters had their MLVA counterparts. All of the observations discussed above suggest that MLVA, with its resolution power and clustering capability comparable to PFGE, may be a good tool for the same applications, especially in outbreak analyses and monitoring shorter-time interhospital spread of S. aureus strains.
Lower levels of correlation between MLVA (cutoff, 70%) and spa typing and MLST at the cluster level could be attributed to obvious differences in discriminatory power between these methods; however, isolates of a specific MLVA type were usually of the same spa type and ST. Exceptional cases of variation (isolates 3121 and N39 in spa and isolates MR1003, 2260/98, 3028, and 303/00 in MLST) were still confined to the same spa cluster (S2 or S4) or MLST clonal complex or cluster (CC8) and were either due to a single deletion/insertion event within the spa locus or due to changes in only one or two of the loci analyzed by MLST. The correlation between MLVA and spa typing at the DNA banding pattern level was especially noteworthy. In order to show dynamics of clone circulation, MLVA and spa typing may be combined, which would decrease the time and cost of typing procedures and allow for interlaboratory comparisons of results. Since all spa types corresponded to specific MLVA patterns, it may be proposed that representatives of MLVA types by spa gene may be routinely typed by using partial sequencing. On the other hand, because spa typing results correlated very well with the MLST data, the representatives of spa types should be good candidates for MLST analysis.
The good agreement observed between MLVA and PFGE on one hand and spa typing and MLST on the other indicates a low level of genetic recombination in S. aureus. Its clonal population structure arises mostly from point mutations (10), which explains the presence of a phylogenetic signal within the S. aureus population (14). However, long-term evolution may be also influenced by chromosomal replacements that complicate the view of the relationships between distantly related isolates. Clusters revealed by spa typing and MLST have not always corresponded in our study. The largest MLST clonal complex CC8 was also the most diverse with six different STs, including ST8 and ST239. This cluster was split into two spa clusters, S1 and S2 (with ST8 and ST239 isolates, respectively), of which cluster S2 was also identified in another MLST group, CC30 (with ST30). Isolates of ST8 and ST239 were considered closely related because they were single-locus variants of each other, differing only in the arcC locus. However, analysis of Robinson and Enright, which included spa and seven surface protein-encoding genes sas, revealed that a large (
557-kb) chromosomal fragment containing all of these genes and arcC was in ST239 identical to that of ST30 and not that of ST8. Therefore, a chromosomal replacement between the ST30 parent and the ST8 parent which resulted in the ST239 mosaic chromosome (23) was most likely responsible for the observed incongruence. The second case was found for the spa cluster S3 and the MLST complex CC1 and cluster ST80. The CC1 isolate 3502 was indistinguishable from isolate BN4 by MLST but different by spa type. On the other hand, its spa type t386 was closely related to the type t044, which is characteristic for three ST80 isolates (cluster S3). There are no detailed data on genome structure of these clones; therefore, we suggest that a chromosomal replacement encompassing the spa locus, followed by its mutational changes, could have also occurred in the evolution of S. aureus CC1 (ST1) clone.
In summary, the main objective of the present study was to compare MLVA with other approaches of typing of S. aureus clinical isolates. Its discriminatory power was found to be lower but comparable with PFGE; moreover, both of these methods were highly concordant in terms of discerning clusters of related and probably related isolates. Therefore, we conclude that MLVA may appear to be useful and reliable in shorter-term epidemiological investigations of S. aureus. The obvious advantage of this method is that it is PCR based, which makes it available for many laboratories, including hospital-based ones. Its good correlation with methods of lower resolution, such as spa typing, suggests that MLVA may also be useful in indicating clonal representatives for larger-scale analyses performed at the Reference Center level. However, this method needs further validation against other typing methods, especially PFGE, with use of other collections of S. aureus isolates. Studies focused on MSSA isolates would also be of interest.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»