Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2004, p. 2054-2064, Vol. 42, No. 5
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.5.2054-2064.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Département "Ecosystèmes et Epidémiologie des Maladies Infectieuses,",1 Génopole, Intégration et Analyse Génomique,2 Laboratoire de Génomique des Micro-Organismes Pathogènes, Institut Pasteur, Paris, France,3 Laboratoire de Microbiologie, Hôpital d'Enfants de Tunis, Tunis, Tunisia4
Received 30 September 2003/ Returned for modification 18 November 2003/ Accepted 31 December 2003
|
|
|---|
|
|
|---|
S. aureus typing is a useful adjunct in several clinical settings, in addition to its use during dramatic acute outbreaks. Despite the use of several phenotypic and genotypic methods (antibiotyping, phage typing, multilocus enzyme electrophoresis, restriction analysis of cellular DNA, analysis of PCR products, and multilocus sequence typing) (10, 13, 22, 24, 31, 32, 35, 36), indistinguishable or closely related isolates have been detected not only among those responsible for outbreaks, but also among those isolated in different countries, at time intervals of several years, and without any obvious epidemiological links. Indeed, Oliveira et al. (27) identified five major pandemic MRSA clones that accounted for almost 70% of the 3,000 isolates analyzed.
The whole genome sequencing of seven S. aureus strains (N315 [19], Mu50 [19], COL [http://www.tigr.org/tdb/], MW2 [2], NCTC8325 [http://www.genome.ou.edu/staph.html], methicillin-susceptible S. aureus strain 476 [http://www.sanger.ac.uk/Projects/S_aureus/], and epidemic MRSA (EMRSA) 16 strain 252 [http://www.sanger.ac.uk/Projects/S_aureus/]) revealed the presence of large amounts of well-conserved DNA regions in the chromosomes. Fitzgerald et al. (11) demonstrated that 2,198 (78%) of the 2,817 COL chromosomal open reading frames (ORFs) represented on a DNA microarray were shared by the 36 analyzed S. aureus isolates from various sources, which belonged to 14 multilocus enzyme electrophoretic types. Ten of the 18 large regions of difference carry genes that encode putative virulence factors and proteins that mediate antibiotic resistance.
The aim of the present study was to design a DNA macroarray with several intragenic PCR amplicons to identify S. aureus at the species level and to type S. aureus isolates. To evaluate the DNA macroarray's usefulness for typing and for the investigation of a putative pathogenicity index correlated with bone infections (BIs), we probed it with cellular DNAs from 80 clinical isolates that were previously typed by the determination of their antibiograms and SmaI restriction patterns. These included unrelated isolates responsible for BIs and isolates from nasal samples of uninfected carriers to check whether these two categories of isolates could be distinguished.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Relevant characteristics of S. aureus clinical isolates
|
Comparative genome analysis, primer design, and PCR amplification. For the annotation and comparative analysis of the available genome sequences from the seven S. aureus isolates cited above, the program CAAT-Box (12) was used. Genes whose nucleotide sequences exhibited <80% similarity were considered distinct. CAAT-Box uses the BLAST program, which presents the area of least similarity with the rest of the genome. The Primer3 program (http://www.broad.mit.edu/cgi-bin/primer/primer3-www.cg) identifies primer pairs in this specific area which are unlikely to produce nonspecific amplifications with regard to the seven sequenced S. aureus genomes. The criteria used by CAAT-Box and Primer3 were as follows: match threshold, 21; maximum length of nonspecific PCR products, 3,000 bases; minimum PCR product length, 250 bases; optimum PCR product length, 400 to 500 bases; primer size, 18, 20, or 25 bases (minimum, optimum, and maximum sizes); primer melting temperature (Tm), 51, 55, or 60°C; % G+C, 25, 50, or 80%; maximum difference in Tm for a primer pair, 5°C.
Each of the 478 selected genes encoded at least 150 amino acids. Primers were designed to amplify a fragment of 400 to 500 bp specific for each gene. Each PCR was performed in a 100-µl reaction volume containing 10 to 20 ng of DNA and a 1 µM concentration of each primer (Eurogentec, Liege, Belgium). The conditions used were an initial cycle of 5 min at 94°C, followed by 35 cycles of 1 min at 94°C, 1 min at 50°C, and 1 min at 72°C, with a final extension step of 7 min at 72°C. The concentration and size of each PCR product were verified by electrophoresis using agarose gels.
Array construction. For array preparation, high-density nylon Performa membranes (Genetix, New Milton, United Kingdom) were soaked in TE solution (10 mM Tris [pH 7.6], 1 mM EDTA). Double spot blots of each PCR product were printed (50 ng of DNA in PCR buffer per spot) by a Qpix robot (Genetix). After spot deposition, DNAs were denatured and fixed on the membranes by incubation for 15 min in 0.5 M NaOH-1.5 M NaCl. The membranes were then washed briefly in distilled water and stored wet at 20°C until use.
Hybridization.
The cellular DNAs of the S. aureus strains (50 ng) were labeled by use of a random priming DNA labeling kit (Roche Diagnostics GmbH, Penzberg, Germany) and 50 µCi of 5'-[
-33P]dCTP (Amersham, Piscataway, N.J.). Labeled probes were purified by use of a QIAquick nucleotide removal kit (Qiagen). The membranes were moistened in 2x SSC (0.3 M NaCl, 0.03 M sodium citrate) and prehybridized for 1 h in 10 ml of 5x SSPE (0.9 M NaCl, 6 mM NaH2PO4, 7.5 mM EDTA, pH 8), 4% sodium dodecyl sulfate, 1x Denhardt solution (0.02% Ficoll, 0.02% polyvinylpyrrolidone, 0.02% bovine serum albumin), and 1 mg of denatured salmon sperm DNA. Hybridization was performed overnight at 65°C. Membranes were washed twice at room temperature and twice at 65°C in 0.5x SSPE-0.2% sodium dodecyl sulfate. Arrays were then sealed in polypropylene bags and exposed to a PhosphorImager screen for 24 h.
Verification of specificity of DNA macroarray. Of the 478 DNA fragments amplified, 106 were randomly chosen and sequenced. Sequencing of the PCR products was done with an ABI3700 capillary sequencer. For a test of correct spotting, the membranes loaded with the amplicons were hybridized with the cellular DNAs of the S. aureus strains used as substrates in PCR amplifications. For 465 amplicons, the results were as expected, i.e., specific. Thirteen of the 478 genes selected were eliminated, either because two nonspecific DNA bands were amplified (1 gene) or because hybridization experiments revealed false-positive or -negative results (10 and 2 genes, respectively). The characteristics of the amplicons and the strains used as substrates, as well as the sequences of the primers and their positions on the genome, are shown in Table S1 (http://genopole.pasteur.fr/staph/).
Data analysis.
For scanning, a Typhoon 9400 PhosphorImager (Molecular Dynamics) was used. Array Vision software (Imaging Research) was used for the quantification of the hybridization intensities and for normalization. For each spot, the hybridization intensity value was normalized by dividing it by the average of all significant intensity values on each membrane. For gene content analysis, a reference array was built by combining the average normalized data of two replicate hybridization experiments with the cellular DNAs of the strains used as substrates for PCR amplification. When a gene was known to be present either as a single copy or as multiple copies, the lowest significant intensity value corresponding to a single-copy gene was chosen. When a gene was known to be present in the tested strain used as a substrate, such as in the five strains whose genomes have been sequenced (N315 [19], Mu50 [19], COL [http://www.tigr.org/tdb/], MW2 [2], and NCTC8325 [http://www.genome.ou.edu/staph.html]), the ratio between the normalized signal intensity of the gene hybridized with the tested strain and that of the reference array was always higher than 0.3. Thus, the threshold for the presence of a gene or a variant related by at least 80% similarity was defined as 0.3. The data were then converted into a binary score as follows: at
0.3, a gene was scored as present (score = 1), and at <0.3, a gene was scored as absent (score = 0).
The binary data were used to cluster the isolates hierarchically, using the program J-Express (9). The threshold adopted to distribute the isolates into clusters was that which enabled each of the outbreak-related isolates belonging to SmaI genotypes 100 or 101 (Table 1) to be grouped and distinguished from any of the other isolates.
Comparative analysis of the gene contents for different categories of isolates.
When categories of n and m isolates are compared, the probability that a given gene is present by chance in n1 isolates of the first category of isolates and n2 isolates of the second category is given by the following binomial formula:
![]() |
|
|
|---|
Furthermore, 67 genes that were not detected in these seven S. aureus genomes were also spotted on the array because they encoded specific groups of proteins. (i) Genes encoding staphylococcal and enterococcal proteins mediating drug resistance were included. Thirteen antibiotic resistance genes were identified in gram-positive species other than S. aureus, as follows: Staphylococcus hyicus, tetL; Staphylococcus cohnii, vatC and vgbB; Staphylococcus epidermidis, fos and lnuA; Enterococcus faecium, vatD, vatE, msrA, lnuB, and vanA; Enterococcus faecalis, vanB and lsa; and Enterococcus gallinarum, vanC. These genes were chosen because of their possible transfer to S. aureus. (ii) Genes encoding factors known to be involved in S. aureus pathogenicity and structurally related proteins (e.g., toxins, adhesins, and enzymes involved in the biosynthesis of capsule or slime) were also included. (iii) Finally, genes encoding proteins produced by mobile elements (transposons, insertion sequences, and plasmids) were spotted on the array. The negative control consisted of an amplicon corresponding to the Staphylococcus intermedius-specific nucI gene (6).
Thus, a total of 465 amplicons were spotted on the membranes. S. aureus strains N315, Mu50, COL, MW2, and NCTC8325 were used to amplify 385 intragenic fragments. The 80 other genes were previously amplified from 39 other strains (see Table S1 in the supplemental material [http://genopole.pasteur.fr/staph/]).
Distribution of the 465 genes among the 80 S. aureus clinical isolates analyzed. The gene content of each of the 80 isolates is given in Table S2 in the supplemental material (http://genopole.pasteur.fr/staph/). Of the 92 genes shared by the seven sequenced genomes and used in the macroarray, 76, including S. aureus nuc and sodM, were detected in all isolates analyzed. Therefore, a total of 388 genes of this set were useful for typing.
Antibiotic resistance genes and phenotypes. An analysis of the data reported in Table 1 enabled us to check whether the genes detected by hybridization were correlated with their phenotypic expression in the isolates. As shown in Table 2, for 79 of the 80 isolates, each antibiotic resistance gene detected was associated with the corresponding phenotype. A single mecA+ isolate was susceptible to ß-lactams.
|
View this table: [in a new window] |
TABLE 2. Antibiotic resistance genes and their corresponding phenotypes in each of the 80 isolates
|
The S. aureus fosB gene, included in the arrays because of its similarity to fos, was found in 69 of the isolates, independent of their phenotypes of resistance to fosfomycin.
The combinations of genes carried by the transposons Tn554 (spe, ermA, tnpA, and tnpB), Tn5406 (vgaAv, tnpA, and tnpB), and Tn4001 (aacA-aphD and IS256 tnp) were found in the isolates exhibiting the antibiotic resistance phenotypes mediated by these transposons. The genes blaZ and tnp480, which are cocarried by Tn552, were associated with only 28 of the 70 isolates containing blaZ. As was stated previously (8), the genes aadE, sat4, and aphA-3, initially found in Tn5405, were always combined, and they were found in seven isolates in this study. This last combination was occasionally associated with other Tn5405 genes, i.e., orfX (two isolates), orfX and IS1182 tnp (four isolates), or orfX, IS1182 tnp, and IS1181 tnp (one isolate).
Distribution of genes in mecA+ isolates and isolates lacking mecA. As shown in Table 1, 36 of the 80 tested isolates were mecA+ and 44 lacked mecA. Several genes, including those coding for antibiotic resistance and putative virulence factors, had a distribution which was significantly different (pg < 0.1) for the two categories of isolates. The distribution of genes encoding putative toxins or adhesins is reported in Table 3. Interestingly, the enterotoxin-encoding genes seg, sei, sem, sen, and seo, codetected in the same pathogenicity island of the S. aureus N315 and Mu50 strains (19), were always associated with each other in our isolates and were significantly predominant in the mecA-negative isolates (30 of 44 isolates) compare to the mecA+ isolates (1 of 36 isolates).
|
View this table: [in a new window] |
TABLE 3. Comparative analysis of the mecA+ and mecA-negative isolates included in this study
|
Clustering of the 80 S. aureus clinical isolates on the basis of their gene contents, as investigated with the DNA macroarray designed for this study. The hierarchical clustering of the isolates by neighbor joining is represented in the dendrogram shown in Fig. 1. First we checked whether the outbreak-related isolates (shown in gray boxes in the figure) were more closely linked to each other than to any of the other isolates.
![]() View larger version (44K): [in a new window] |
FIG. 1. Hierarchical clustering of the 80 S. aureus isolates investigated according to their gene contents by the J-Express program (9). The threshold was chosen to distinguish each of the outbreak-related isolates belonging to SmaI genotype 100 or 101 (clusters 9 and 28, respectively) from any of the other isolates.
|
The other outbreak-related isolates belonged to SmaI subtype 39aA (IPF 555, IPF 557, and IPF 562) (14) or 45aA (BM 9586, BM 12184, nad BM 12188) (22) and were isolated in hospital C (Paris) in 1999 and 1987, respectively (Table 1). Isolates belonging to SmaI genotype 39A were phenotypically recognizable because of their decreased susceptibility to glycopeptides. Those belonging to SmaI genotype 45A and phage type 77 were initially discovered in 1987, during the emergence of resistance to fluoroquinolones in French hospitals. Such isolates preexisted in European hospitals before these outbreaks, but at very low frequencies. In this study, we analyzed 24 mecA+ isolates belonging to SmaI genotypes 39A and 45A that were isolated in several European countries and at time intervals of several years. These endemic isolates, which are considered putatively related according to their SmaI genotypes, were more linked to each other than to any of the 56 other isolates (Fig. 1). Note that some of them are clearly divergent in the dendrogram and that the mode of their linkage is not correlated to their SmaI genotype, but those considered to be outbreak related are closely linked.
Clustering of the 80 clinical isolates after choice of threshold for hierarchical clustering dendrogram. For the distribution of the isolates into clusters, it was necessary to choose a threshold for the dendrogram. For this purpose, the threshold adopted was that which enabled each of the outbreak-related isolates belonging to SmaI genotype 100 or 101 to be distinguished from any of the other isolates. These isolates were taken into consideration because they were not detected before the outbreaks, in contrast to the SmaI subtype 39aA or 45aA outbreak-related isolates. The choice of this threshold enabled the discrimination of 52 clusters belonging to 45 SmaI genotypes among the 80 isolates (Fig. 1). In Table S2 in the supplemental material (http://genopole.pasteur.fr/staph/), the genes are listed according to the clusters to which they belong.
With the selected threshold, a total of five clusters were found among the 10 SmaI type 39A isolates and eight clusters were found among the 14 SmaI type 45A isolates (Fig. 1). Among these isolates, which are endemic in European cities, those collected in the same hospital or city were not necessarily the most closely linked. The three outbreak-related SmaI subtype 39aA isolates collected in hospital C (Paris) in 1999 (IPF 555, IPF 557, and IPF 562) are linked in cluster 41, which includes another SmaI subtype 39aA isolate (BM 12612) collected at Villiers St. Denis in 1998. Moreover, four of five isolates belonging to two SmaI subtypes, 45aA and 45dA, and collected in three French hospitals in 1987 are within cluster 46 (BM 9290, BM 9343, BM 9586, and BM 12184). The fifth isolate, BM 12188, located in the separate but close cluster 47, was distinguishable by the lack of five drug resistance genes, namely blaZ, qacA, qacC, CZ040, and CZ041, encoding ß-lactamase, resistance to antiseptics, organomercurial lyase, and mercuric reductase, respectively. Figure 2 shows the images resulting from scanning of the two DNA macroarrays hybridized with the total cellular DNAs from the BM9290 and BM12188 isolates (Table 1).
![]() View larger version (90K): [in a new window] |
FIG. 2. Images resulting from scanning of the DNA macroarrays hybridized with the total cellular DNAs from two isolates. (A) Isolate BM9290 (cluster 46). (B) Isolate BM12188 (cluster 47). Even though they belonged to the same SmaI subtype (45aA), the two isolates were found in two close but separate clusters (Fig. 1) due to the lack, in BM12188, of the following five drug resistance genes: blaZ, qacA, qacC, CZ040, and CZ041, encoding ß-lactamase, resistance to antiseptics, organomercurial lyase, and mercuric reductase, respectively.
|
|
|
|---|
The detection of antibiotic resistance genes is particularly interesting when these genes mediate low antibiotic resistance levels that are not reproducibly detectable by antibiograms. This level of detection also contributes to the selection of isolates that carry genes that have not yet been described. By hybridization with 400- to 500-bp amplicons, mutations in preexisting genes associated with antibiotic resistance cannot be visualized and would necessitate hybridization with oligonucleotides. For 79 of the 80 clinical isolates tested, the resistance phenotype conferred by each of the detected resistance genes was expressed, whereas one mecA+ isolate was susceptible to ß-lactams. This high correlation demonstrated an extensive and satisfactory choice of antibiotic resistance genes spotted on the membranes. For the two related streptogramin A-resistant isolates, the lack of any known staphylococcal or enterococcal gene conferring resistance to this antibiotic is probably due to the presence of a gene(s) that has not yet been described.
The assessment of the presence of all known S. aureus genes encoding putative virulence factors may contribute to the determination of the pathogenic potential correlated with particular types of infection and to the identification of emerging pathotypes. In this study, we checked whether some genes were more prevalent in isolates responsible for BIs than in isolates from uninfected NCs. For this purpose, only unrelated isolates from our collection were included. This constraint explains why the numbers of isolates analyzed were 16 BI isolates and 12 NC isolates. Despite the fact that BIs were contracted by children outside the hospital, several patients were infected by S. aureus isolates that were considered monoclonal on the basis of their SmaI patterns. Although a few genes, including sdrD, encoding a putative SD adhesin, appeared predominant in one of the two categories of isolates, the differences were not significant when the 388 genes used for typing were taken into account for the calculation of the probability that a given gene is present by chance. Thus, a larger number of unrelated isolates from various sources merits further analysis. However, when only the 11 genes encoding putative adhesins were taken into account, the higher prevalence of sdrD in BI isolates than in NC isolates became significant. Some SD proteins were shown to bind fibrinogen (ClfA [21], ClfB [25], and SdrG [16]) or bone sialoprotein (Bbp) (33), but the ability of SdrD to bind a matrix protein(s) has not been investigated. The impact of sdrD inactivation merits evaluation in an animal model of BIs.
The significantly distinct distribution of some genes encoding enterotoxins or adhesins among the mecA+ and mecA- negative isolates in this study (Table 3) may not be the case among isolates from various sources. Indeed, most of the 80 isolates tested were collected in France and Tunisia, and the mecA+ isolates belonged to a limited number of SmaI genotypes. Nevertheless, the low frequency of cna detection in mecA+ isolates has been reported already by Booth et al. (5).
Due to the use of a large number of genes for typing (388), all 80 isolates tested were typeable. A method based on the analysis of a large number of genes was expected to yield more discrimination between the isolates than the typing methods based on sequencing of a limited number of genes or on the analysis of SmaI patterns, which depends on the number and locations of SmaI sites in the genome. This was confirmed by this study, for the mecA+ isolates were endemic to several European cities and were collected at large time intervals (SmaI genotypes 39A and 45A). Among the latter isolates, those considered to be outbreak related in the same hospital were found in the same or in a close cluster(s): cluster 41 or 46-47. In such a context, the typing method proposed in this study provides more discrimination of the isolates responsible for acute outbreaks than the determination of SmaI patterns. For the other isolates, if we excluded the three pairs which were linked in the same cluster despite belonging to unrelated SmaI types, our results revealed a correlation between the modes of isolate clustering based on the two typing methods, i.e., the analysis of gene contents and the SmaI patterns. Indeed, the isolates belonging to the same or related SmaI types appeared to be more linked to each other than to those belonging to unrelated SmaI types.
In conclusion, the typing method proposed here performed better than that based on the analysis of SmaI patterns, in particular for distinguishing outbreak-related isolates from those that are endemic to a particular area. It also has the advantages of being faster and providing additional information concerning the gene contents of interest. This macroarray should be updated when additional genes are described and also needs to be validated for the analysis of the transcription of genes in order to evaluate the levels of gene expression which may be correlated with particular types of infections. The method described here can also be performed with glass slides and fluorescent labeling in order to be more amenable to automation for routine analyses.
We thank the biologists and NARSA who provided several of the strains used in this study and Iain Old for reviewing the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»