Previous Article | Next Article ![]()
Journal of Clinical Microbiology, January 2008, p. 118-129, Vol. 46, No. 1
0095-1137/08/$08.00+0 doi:10.1128/JCM.01685-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Laboratory of Molecular Genetics, Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Rua da Quinta Grande, nr. 6, 2780-156 Oeiras, Portugal,1 Knowledge Discovery in Bioinformatics Group, Instituto de Engenharia de Sistemas e Computadores: Investigação e Desenvolvimento, Rua Alves Redol nr. 9, Apartado 13069, 1000-029 Lisboa, Portugal,2 Department of Infectious Disease Epidemiology, Faculty of Medicine, Imperial College London, Old Medical School Building, St Mary's Hospital, Norfolk Place, London W2 1PG, United Kingdom,3 UEI Micobactérias, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Rua da Junqueira, nr. 96, 1349-008 Lisboa, Portugal,4 Centro de Recursos Microbiológicos, Faculdade de Ciências e Tecnologia, Quinta da Torre, 2829-516 Monte de Caparica, Portugal,5 Laboratory of Microbiology, The Rockefeller University, New York, New York 100216
Received 23 August 2007/ Returned for modification 13 October 2007/ Accepted 25 October 2007
|
|
|---|
|
|
|---|
The basic principle of epidemiological typing is that isolates of an infectious agent that belong to the same chain of transmission are clonally related; that is, they descend from a common ancestor. During their evolutionary history, isolates within a given species diversify through point mutations, recombination, or the acquisition/deletion of mobile genetic elements, giving rise to extensive genomic and phenotypic diversity. Since the rates for mutation, recombination, and transfer of mobile elements as well as environmental pressures vary from species to species, the typing method selected and the threshold of marker similarity used to define a clone need to be adjusted to the species under investigation. Furthermore, space and time also need to be considered when the optimal epidemiological markers are selected. Ideally, typing methods should be reproducible and stable; should have a high discriminatory power and epidemiological concordance; should be versatile, easy to perform, easy to interpret; and should be cost-effective and time effective.
Staphylococcus epidermidis is one of the most important nosocomial pathogens associated with catheter-related and other indwelling medical device-related infections. Moreover, approximately 70% of the S. epidermidis strains isolated in the hospital environment have acquired resistance to methicillin, and the majority of them are also resistant to almost all classes of antimicrobial agents (10, 20, 26). The most widely used method for characterizing S. epidermidis isolates is pulsed-field gel electrophoresis (PFGE), which has mainly been deployed to address short-term epidemiology issues. The molecular characterization of nosocomial S. epidermidis isolates by PFGE has revealed considerable genetic diversity within the population (1, 6, 11, 25, 31). Despite the diversity that has been observed, PFGE has been used to track the dissemination of particular S. epidermidis strains between different patients (27), wards (47), and hospitals (35), thus providing a valuable tool for the control of S. epidermidis infections in clinical settings.
It was only recently that some key features of the evolutionary history, population structure, and long-term global epidemiology of S. epidermidis were revealed when an improved multilocus sequence typing (MLST) scheme was developed for S. epidermidis (41). Despite the distinctive genetic diversity of this species, the molecular characterization of a geographically and temporally diverse collection of S. epidermidis isolates by MLST revealed that the nosocomial S. epidermidis population was composed of nine clonal lineages that were found to be spread worldwide (30). Furthermore, in contrast to Staphylococcus aureus, in which the early stages of clonal diversification were found to occur mainly by point mutations (8), analysis of S. epidermidis MLST data indicates that in S. epidermidis this occurs primarily by recombination (30).
Another method that has recently helped to clarify the molecular epidemiology and evolution of methicillin-resistant S. epidermidis (MRSE) strains is the typing of the staphylococcal chromosomal cassette mec (SCCmec), which is the mobile element carrying the determinant of methicillin resistance (mecA) (21). The molecular characterization of this element in diverse collections of S. epidermidis isolates identified five SCCmec types previously identified in S. aureus (18, 19, 32), as well as a high number of nontypeable and new SCCmec types (15, 16, 28, 30), indicating a high degree of genetic diversity within the SCCmec elements carried by S. epidermidis. Moreover, several findings, such as the high degree of variability observed in the specific region of integration of SCCmec (orfX) (28) and the high degree of diversity found in the type of SCCmec associated with a particular PFGE type or sequence type (ST) (28, 30), support the frequent acquisition and/or loss of SCCmec by S. epidermidis.
Recent discoveries regarding the population structure of S. epidermidis, namely, that the genetic diversity observed between S. epidermidis strains may derive from the rapid evolution of the chromosome through recombination and the frequent transfer of SCCmec, have a pronounced effect on the interpretation and understanding of the data generated by typing methods. Although a large volume of data has recently been collected through the characterization of S. epidermidis strains by PFGE, MLST, and SCCmec typing, no study in which the three methods were compared has been performed. Moreover, no guidelines exist on the optimal epidemiological markers or the threshold of marker similarity to be used for the definition of an S. epidermidis clone.
In the present study we assess quantitatively the strength of the correlation between the results of PFGE, MLST, and SCCmec typing when characterizing S. epidermidis strains through the calculation of measures of concordance: the adjusted Rand coefficient (AR) (17, 36) and the Wallace coefficient (W) (45). We anticipate that the finding of any congruence between the type assignments determined by these methods, which target different regions of the chromosome that evolve at different rates, will give an indication of the most appropriate typing method or combination of methods to be used for the definition of a clone in S. epidermidis.
|
|
|---|
S. aureus strains COL (33), N315 (24), HU25 (33), and WIS (19) were included as controls for SCCmec typing.
DNA preparation. Agarose disks containing chromosomal DNA and genomic DNA for PCR were prepared as described previously (5, 28, 29).
PFGE. SmaI PFGE patterns were determined by conventional protocols (5, 29). Visual assignment of the PFGE types and subtypes was performed according to established criteria (40). PFGE types were identified by letters, and subtypes were identified by letters followed by a numeric subscript.
A database containing all the SmaI PFGE patterns was created by using Bionumerics software (version 4.5) from Applied Maths (Sint-Martens-Latem, Belgium), where band patterns over the multiple gels were normalized and compared. Clustering was performed by using the Dice similarity coefficient and the unweighted pair group method with arithmetic means (UPGMA), with 1.3% of tolerance and 0.8% optimization. PFGE types were automatically assigned (the PFGE type AA method) by using cutoff similarity values that varied from 60 to 100%. The types obtained by the PFGE type AA method are represented by numbers.
MLST. MLST was carried out by using the new MLST scheme described by Thomas et al. (41). The internal fragments of seven housekeeping genes were amplified by PCR with primers whose sequences match highly conserved regions, and both strands of all amplicons were sequenced with an ABI Prism 3700 DNA sequencer by using BigDye (version 3) fluorescent terminators. Numbers for alleles and STs were assigned according to the S. epidermidis MLST database (http://sepidermidis.mlst.net/).
eBURST algorithm. The most likely patterns of evolutionary descent in our collection were assessed by using the eBURST algorithm (http://eburst.mlst.net) (9) with previously defined settings (30). Clonal complexes are represented by the abbreviation CC, and singletons are represented by the abbreviation S. CC2 was subdivided into clusters I and II, as proposed previously (30), and cluster II was further separated into subclusters, represented by the cluster number followed by the subgroup founder number (II-5, II-6, II-85, and II-89).
Analysis of SCCmec structure. The structures of the ccr and the mec complexes were determined by conventional PCRs, as described previously (19, 32). SCCmec types were defined by the combination of the type of ccr complex and the class of mec complex. SCCmec was considered nontypeable when the ccr or the mec complex, or both, were nontypeable, according to previously defined criteria (30).
Genotypic diversity. Genotypic diversity was calculated by using Simpson's index of diversity (SID) (39). Confidence intervals were calculated as described by Grundmann et al. (13).
Clustering concordance coefficients: AR and W. For comparison of the congruence between the type assignments determined by the different typing methods, AR (17, 36) and W (45) were calculated and applied as described previously (3). AR values close to 0 indicate a lack of congruence between the methods, while values close to 1 indicate a high level of congruence. W values close to 1 indicate that the results obtained by a given typing method enable one to make a good prediction of the clustering determined by the other method.
|
|
|---|
![]() View larger version (73K): [in a new window] |
FIG. 1. Dendrogram, constructed by the UPGMA method, representing the similarity of the S. epidermidis SmaI PFGE DNA restriction profiles identified by the PFGE type AA method. Representatives of all the 77 PFGE types identified by the PFGE type AA method among the collection of 216 isolates when a cutoff of 79% was used are shown. Information on the PFGE type obtained from the PFGE VA type, SCCmec type, ST, and CC for each strain is included.
|
|
View this table: [in a new window] |
TABLE 1. PFGE, MLST, and SCCmec results for strains belonging to CC2, minor CCs, and singletons identified in the strain collection
|
The highest overall congruence between the PFGE type VA method and the comparison performed with Bionumerics software was obtained by using a cutoff similarity value of 79% (ARcutoff 79% = 0.3576). It was also at the cutoff value of approximately 79% that the highest value of W in both reading directions was obtained, meaning that this is the cutoff value at which both the PFGE type AA and the PFGE type VA method assignments provide the best prediction possible of the assignment previewed by the other method (WPFGE type VA
PFGE type AA = 0.3970; WPFGE type AA
PFGE type VA = 0.3662).
By considering a cutoff similarity value of 79% in the UPGMA dendrogram, 77 PFGE types (PFGE types 0001 to 0077) were defined (Fig. 1). The most common PFGE types in the collection analyzed were PFGE types 0016 (16%) and 0010 (7%), followed by PFGE types 0001 (6%) and 0024 (5%). The remaining PFGE types were represented by less than 5% of the population. The majority of the strains assigned to the most-represented PFGE types by automatic classification (PFGE types 0001, 0010, and 0016) were also classified by the visual assignment as belonging to the most common PFGE types (PFGE types A, D, C, and FF) (see the data delimited by the gray square in Table SA1 in the supplemental material). For the subsequent analysis we maintained the visual assignment of SmaI PFGE DNA band patterns.
Molecular characterization by MLST. Seventy-four different STs were identified among the 216 S. epidermidis isolates. The most represented ST was ST2, comprising 31% of all isolates (66 of 216 isolates), followed by ST59 (6%) and ST23 (5%). Each of the remaining STs comprised less than 5% of the isolates (Tables 1).
The eBURST algorithm separated the 74 STs into one major CC (CC2), 8 minor CCs (CCs 1, 11, 21, 23, 33, 42, 49, and 66), and 13 singletons.
We have previously proposed that CC2 be subdivided into two different clusters (clusters I and II), based on the finding that two groups of STs within this CC had contrasting recombination rates, distinct phylogenetic congruence, and different levels of complexity on descendence patterns (30). Turner and colleagues (42) recently suggested that in species for which a large CC was defined by eBURST analysis due to high recombination rates, like S. epidermidis, radial subgroups within the CC corresponded reasonably well with the true ancestry, whereas the interlinking of those subgroups was probably incorrect. In order to illustrate better the true ancestry, we propose that cluster II of the large CC2 be further subdivided into four subclusters that correspond to the most radial subgroups, each containing as its ancestor a subgroup founder: subgroups II-5, II-6, II-85, and II-89 (Fig. 2).
![]() View larger version (16K): [in a new window] |
FIG. 2. Proposal for subdivision of S. epidermidis CC2 identified by eBURST analysis. Each ST is represented by a black dot. White and gray dots correspond to group and subgroup founders, respectively. Single-locus variants are linked by lines, and the CC or eBURST group corresponds to the group of connected STs. It has been proposed previously (30) that CC2 be subdivided into two clusters (clusters I and II) for phylogenetic reasons. In order to illustrate better the true ancestry, we propose that cluster II of the large CC2 be further subdivided into four subclusters that corresponded to the most radial subgroups, each containing as its ancestor a subgroup founder: II-5, II-6, II-85, and II-89 (encircled by gray lines).
|
Nearly half (41%) of the 138 S. epidermidis isolates harbored SCCmec type IV, whereas 27% carried SCCmec type III. On the other hand, a limited number of isolates carried SCCmec type V, I, or II (6%, 4%, and 4%, respectively). In addition, as many as 18% of the isolates analyzed carried nontypeable SCCmec structures (9%) or SCCmec structures with new associations between the ccr complex and the mec complex (9%).
Genotypic diversity. The most discriminatory method used was PFGE (SID = 96.43% for the PFGE type and SID = 99.88% for the PFGE subtype), followed by MLST (SID = 89.64%), eBURST analysis (SID = 77.70%), and SCCmec typing (SID = 75.55%) (Table 2). When the information provided by the PFGE type or subtype was complemented with the information provided by SCCmec typing for 138 isolates, no significant increase in the discriminatory power was obtained (Table 2).
|
View this table: [in a new window] |
TABLE 2. Number of types and genotypic diversity of each typing method
|
With respect to the SCCmec type distribution among the MRSE isolates in this study, we found that strains with the same CC, ST, PFGE type, and even PFGE subtype could harbor different SCCmec types. Strains of ST2 and ST22 had the highest number of different SCCmec types (five different types each), whereas strains belonging to the remaining 72 STs carried between one and three different SCCmec types.
SCCmec type IV was carried by strains belonging to as many as seven different CCs (CCs 1, 2, 11, 21, 23, 49, and 66) and three singletons (S65, S72, and S82), including 37 STs and 31 PFGE types. Conversely, SCCmec types I, II, III, and V were present in a limited number of different genetic backgrounds: SCCmec type I was present in strains within CC42 (PFGE type C), SCCmec type II was detected in strains belonging to CC2 cluster I and cluster II-5 (PFGE types FF, H, and V) and CC11 (PFGE type X), SCCmec type III was found to be exclusively associated with strains belonging to CC2 (PFGE types A, C, D, FF, G, KK, N, PP, QQ, and V), and SCCmec type V was carried by strains within CC2 clusters I and II-5 (PFGE types Q and C) and S56 (PFGE type T).
Quantitative measure of clustering concordance between the different typing methods.
To measure the congruence between the type assignments of the different typing methods, AR and W were calculated for a subset of 138 isolates for which the results of all typing methods were available. The values of AR obtained indicated that among the methods used in this study to type the S. epidermidis strains, the highest overall congruence obtained was between MLST and PFGE (ARST
PFGE type VA = 0.3119; ARST
PFGE type AA = 0.5152) (Table 3). The overall congruence between SCCmec typing and MLST was comparably low (ARSCCmec
ST = 0.1879), and was even lower between SCCmec typing and PFGE (ARSCCmec
PFGE type VA = 0.0833; ARSCCmec
PFGE type AA = 0.1350).
|
View this table: [in a new window] |
TABLE 3. AR for the methods used to characterize the 138 S. epidermidis isolates
|
ST < 0.65); however, a high value of W was obtained for the clustering performed by PFGE and eBURST analysis (WPFGE type VA
CC = 0.7914; WPFGE type AA
CC = 0.8529), indicating that two strains with the same PFGE type have a high probability of belonging to the same CC. When the information provided by PFGE type was complemented with the information from SCCmec typing, an even better correlation with the eBURST analysis clustering was obtained (WPFGE type VA + SCCmec
CC = 0.9212; WPFGE type AA + SCCmec
CC = 0.9379). |
View this table: [in a new window] |
TABLE 4. W for the methods used to characterize 138 S. epidermidis isolates
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Representation of the strength of concordance between the typing methods used to type a subset of 138 S. epidermidis isolates, as measured by the use of W. Circles include the type assignments performed by the typing methods used in this study: PFGE type VA, SCCmec type, ST, or CC or combinations of these typing methods. Arrows of different widths correspond to different values of W (see legend).
|
|
|
|---|
In an attempt to detect an accurate phylogenetic signal that would give confidence in the interpretation of the data and provide valuable information on the most appropriate method to be used for the definition of a clone, we performed a quantitative assessment of the congruence between three distinct typing methods when they were used to characterize the same S. epidermidis isolate collection.
SCCmec acquisition and recombination as sources of genetic diversity. The molecular characterization of S. epidermidis isolates by PFGE, MLST, and SCCmec typing revealed a high degree of genetic diversity (0.75 < SID < 0.99), as was previously described by other authors using different typing methods (7, 11, 15, 16, 28, 30, 31, 41, 43, 46, 49). Nevertheless, the level of genetic diversity obtained in this particular study may be biased toward a higher level of diversity in relation to what would usually be found in the clinical setting, since the sample analyzed was initially selected to include a very diverse collection of isolates (30).
To date, the origins of S. epidermidis genetic diversity have been explored only partially. However, according to the most recent findings on the epidemiology of S. epidermidis (28, 30), the high rates of recombination and the frequent acquisition of mobile genetic elements appear to be the most likely driving forces generating genomic polymorphisms in this species.
PFGE demonstrated the most diversity (SID = 99.88), reflecting the sensitivity of the method to detecting changes throughout the genome. The extraordinarily fast pace at which S. epidermidis macrorestriction patterns evolve is well illustrated by the finding of chromosomal rearrangements in S. epidermidis strains isolated from single patients during infection episodes (11, 49). The frequent insertion/excision of IS256 (49) and SCCmec (28, 30) may be contributing to the extensive polymorphism in the number and location of SmaI restriction sequences in the S. epidermidis chromosome (28, 49).
The genotypic diversity of MLST observed (SID = 89.64%) may be a result of the shuffling of alleles among different S. epidermidis strains, promoted by recombination, which, in conjunction with the existence of a diverse pool of alleles for each locus, may have led to the emergence of a high number of STs (30).
The origin of the high degree of genetic diversity of SCCmec carried by S. epidermidis observed in this and other studies (15, 16, 28, 30) is still unclear. The existence of regions of homology between different SCC elements (4, 22) suggests that the genetic diversity in SCCmec may mainly be generated by recombination events. It is possible that the high recombination rates and the high number of SCCmec acquisitions found in S. epidermidis (30), together with a high frequency of SCCmec carriage (10, 20, 26), may increase the potential for the generation of new SCCmec structures in this species.
Despite the genetic diversity in the structure of SCCmec carried by S. epidermidis observed in this study, SCCmec type IV was carried by nearly half of the isolates, comprising as many as seven different S. epidermidis clonal lineages. The predominance and widespread distribution of SCCmec type IV in S. epidermidis suggests that this species may be a potential reservoir of SCCmec type IV for other staphylococcal species.
Usefulness of automated assignment of PFGE types.
In the present study, we attempted to give some indication of the most appropriate parameters for automated PFGE type assignment of S. epidermidis SmaI restriction band patterns using Bionumerics software. Our results showed that the parameters for which the type obtained by the PFGE type AA method best matched that obtained by the PFGE type VA method were a tolerance of 1.3%, an optimization of 0.8%, and a cutoff value of 79%, which agree with the settings defined for other gram-positive pathogens, such as Streptococcus pneumoniae (2) and Streptococcus pyogenes (3). Although it can be argued that the values of the concordance coefficients between the types obtained by PFGE type VA and AA obtained by using the parameters mentioned were low (AR = 0.3576; WPFGE type VA
AA = 0.3970; WPFGE type AA
VA = 0.3662), we should call the attention to the fact that we were analyzing a collection of isolates biased toward a higher degree of genetic diversity. The existence of such diversity in the SmaI PFGE restriction band patterns makes automatic band matching between different strains particularly difficult to obtain. Moreover, the correctness of the automatic band matching performed may have been influenced by differences in the PFGE gel runs that may have generated systematic band shifts, which might have been interpreted as differences. We anticipate that higher values of concordance would be obtained if a more homogeneous collection of isolates was analyzed or if the PFGE gels were rerun under more standardized conditions.
The conditions for the automation-based PFGE type definition proposed here are specific to this study collection; however, we predict that the same or similar settings should perform equally well for data generated by the same PFGE protocol, which will generate similarly resolved band patterns. The use of the PFGE type AA method, specifically in a species holding such a high degree of genetic diversity at the level of macrorestriction band patterns, such as S. epidermidis, is particularly welcome, since the PFGE type VA method can become particularly difficult and subject to error, especially when large numbers of PFGE macrorestriction patterns are compared. In this study, the usefulness of the automated system is well illustrated by the higher values of the concordance coefficients obtained between PFGE and the other methods when the PFGE type AA method was performed compared to the values obtained by the PFGE type VA method, which indicates that some errors may have been introduced during the visual inspection of the band patterns.
Ability to predict CC from the combination of PFGE and SCCmec typing results. The more significant and interesting results from our study were obtained when AR and W were calculated for the quantitative measurement of the correlation between the different methods or combinations of methods. No good correlation between PFGE type, ST, and SCCmec type was detected (AR and W < 0.75). The lack of correlation between the three typing methods is apparent in the dendrogram constructed from the types obtained by PFGE type AA, where strains with the same SCCmec, ST, or CC are clustered apart (Fig. 1). The high levels of recombination that occur in S. epidermidis species (30) and the frequent acquisition of SCCmec by strains harboring the same ST and even the same PFGE type (28, 30) may be the explanation for this lack of an association.
Although no good concordance between the three typing methods tested was established, an unexpectedly high correlation was found between the type obtained by PFGE type VA or AA and the CC, and the strength of this correlation was augmented when the data generated by PFGE type assignment were complemented with those generated by SCCmec typing. This result suggests that the same strong clonal signal is being detected by both the PFGE type-SCCmec combination and eBURST analysis of the MLST data.
The explanation for such a correlation is not obvious; however, we can speculate that despite the frequent insertion/excision of SCCmec from the S. epidermidis chromosome, the rate of SCCmec acquisition by S. epidermidis must be lower than the rate of occurrence of events necessary for the emergence of new PFGE subtypes, allowing the detection of a strong phylogenetic signal when information from PFGE and SCCmec typing are combined. This indicates that, at least transiently, a specific SCCmec type is associated with a strain with a characteristic SmaI PFGE restriction profile. The ability of PFGE to successfully identify genetic lineages was already demonstrated for clonal species like S. aureus (34). However, it was surprising to verify that such an ability was also observed in a species like S. epidermidis, which holds such a high degree of genetic diversity in its SmaI DNA restriction band patterns. The fact that the strength of this correlation was amplified when the SCCmec typing results were added may be related to the existence of exclusive restriction-modification systems. The constraint on the exchange of genetic material between strains of different lineages provided by restriction-modification systems may have promoted the existence of specific sets of SCCmec types inside each CC, as suggested previously (23, 30, 44).
A proposal for clone definition. Taking into consideration these new observations, we propose that clones within the S. epidermidis species be defined by the combination of the PFGE type followed by the SCCmec type, since among all the methods and combinations tested, the results of these methods were the most concordant with those of all other typing methods, providing reliable information on the short-term epidemiology and the ability to predict with consistency the long-term evolutionary history. Besides providing the more meaningful epidemiological information, both methods have high discriminatory powers, are easy to perform in the clinical laboratory, and are cost-effective. The major weakness in our proposal of clone definition probably resides in the difficulty of achieving the interlaboratory reproducibility of PFGE band patterns, which can hinder the interchange of data between different laboratories. For S. aureus, it was demonstrated that the laboratory-to-laboratory reproducibility of PFGE is possible if some minor technical problems are overcome (5). In order to increase the interlaboratory reproducibility and standardize the analysis of the S. epidermidis SmaI PFGE band patterns, we recommend that identical PFGE protocols be applied, that prototype strains of the most important S. epidermidis PFGE types be included in each PFGE run, and that automated analysis of the PFGE band patterns with the same or similar settings be performed.
For the sake of future comparisons we propose that S. epidermidis clones be identified by the type obtained by the PFGE type AA method, followed by the SCCmec type, and the CC number (with reference to cluster/subclusters in case of CC2) (e.g., 28-V-CCII5). However, the fact that S. epidermidis has a highly dynamic chromosome, with the consequent frequent emergence of new PFGE types and the finding of new SCCmec structures, should be kept in mind. Therefore, the determination of STs for the new PFGE type-SCCmec type combinations should be performed.
The work presented here illustrates well how the integration of data from a correctly curated database of PFGE band patterns with data obtained by other typing methods and the knowledge of congruence between them can guide the choice of the most appropriate and effective approach for inferring clonal relationships between S. epidermidis isolates. We believe that the application of this same methodology to other important clinically relevant bacteria will allow a broader view of how the results of different typing methods are related and how they can be used to effectively differentiate clonally related strains from unrelated strains.
We thank K. Hiramatsu for providing strain N315 and T. Ito and W. Grubb for providing strain WIS, which were included in this study.
Published ahead of print on 7 November 2007. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
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»