Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, July 2006, p. 2506-2515, Vol. 50, No. 7
0066-4804/06/$08.00+0 doi:10.1128/AAC.01640-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Chemical and Biological Engineering, University of Colorado, Boulder, Colorado 80309
Received 27 December 2005/ Returned for modification 20 February 2006/ Accepted 7 April 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
We chose to study a clinical isolate, named B1, that showed high levels of resistance to five aminoglycosides tested: amikacin, gentamicin, kanamycin, streptomycin, and tobramycin. B1 was found to be of serotype O12, a serotype that has been associated with multiresistance to a number of antibiotic classes, including aminoglycosides and beta-lactams (4, 15, 46, 65, 66). B1 was chemically mutagenized and screened to identify mutants that exhibited increased levels of susceptibility to amikacin. Based on the fact that resistance can arise due to genetic mutation (2, 13, 14, 29), we expected that it would also be possible to restore sensitivity to an already resistant strain through mutation. We hypothesized that the frequency of finding mutations that increase sensitivity would be higher than the frequency of isolating mutants with increased resistance. This is based on the premise that, with only a limited number of point mutations providing a selective advantage in the presence of an antibiotic, a larger percentage of mutations may either have no impact or decrease the level of antibiotic tolerance. Although the frequency of finding sensitivity mutations is a key consideration in assessing the potential of combination therapies, the ability to identify sensitivity-restoring genes/mechanisms is the critical parameter in the design of such therapies.
We hypothesized that whole-genome transcriptional profiling with Affymetrix PAO1 GeneChips could be used for such a purpose. To test this hypothesis, we compared the expression patterns among strain B1, two of its susceptible mutants (named M5 and M31), and laboratory strain PAO1 (81). Moreover, we assessed the resistance mechanisms suggested by transcriptional profiling by MIC assays performed under a variety of conditions (in the presence of verapamil, carbonyl cyanide m-chlorophenylhydrazone [CCCP], and polymyxin B and with spheroplasts) as well as PCR-based gene identification approaches.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Disk susceptibility and MIC determination. Susceptibility to aminoglycosides (amikacin, gentamicin, kanamycin, streptomycin, tobramycin) was determined by methods similar to the disk diffusion method described previously (58) with LB agar plates. LB agar plates were inoculated by swabbing with a turbid culture. Immediately following inoculation, BBL Sensi-Discs (Fisher Scientific, Pittsburgh, PA) were placed with sterile forceps onto the plate. The plates were incubated overnight, and the diameters around the disks were measured.
MICs were determined by a standard broth dilution method in LB medium (1). The lowest concentration at which no growth was noted after 18 h was deemed the MIC. Cultures at an optical density at 600 nm of 0.6 were used for MIC studies For preparation of the medium to be used in studies of MICs determined in the presence of CCCP or verapamil, either CCCP was added to LB medium to achieve a 250 µM solution (87) or verapamil was added to LB medium to achieve a final concentration of 100 µg/ml. This medium was then used in the standard broth dilution MIC assay, with the appropriate levels of amikacin added.
For MIC assays performed with spheroplasts, spheroplasts of P. aeruginosa were formed by treating the cells with EDTA and lysozyme (77). Spheroplast formation was checked by osmotic shock (5, 59).
Random mutagenesis and sensitivity screening. Random mutagenesis with N-methyl-N'-nitro-N-nitrosoguanidine (MNNG; TCI America, Portland, OR) was carried out as described previously (53). Cells were exposed to MNNG at a concentration of 50 µg/ml for 8 min. This level of mutagenesis was sufficient to result in the survival of 50% of the exposed cells. The mutated cells were then plated onto LB agar plates. The contents of the master plates were then replica plated (44) onto LB agar plates containing increasing levels of amikacin. MIC studies were carried out with the colonies that grew on the master plate but not on the imprints to confirm sensitive phenotypes.
For studies conducted to examine the frequencies of sensitivity-restoring mutations, mutated cells were plated onto LB agar plates and incubated overnight. The cells were then patched into 384-well plates containing LB medium. These plates were then patched by using a 384-pin blot replicator into 384-well plates containing increasing amounts of amikacin. After 18 h of incubation, growth in the 384-well plates was confirmed by once again replica plating the cells from the amikacin-containing plates into plates containing LB medium with no amikacin. These plates were then incubated, and growth was determined visually after 18 h.
Detection of plasmids. Strain B1 was examined for the presence of plasmids. No plasmids were found from plasmid extractions with a QIAprep Spin Miniprep kit (QIAGEN, Valencia, CA) or the hot alkaline method of Kado and Liu (36) which has been shown to be effective for the isolation of plasmids in the size range of 2.6 to 350 MDa.
LPS gels. Lipopolysaccharide (LPS) was isolated by first lysing the cells and treating the lysate with proteinase K (26). LPS was run on 16.5% acrylamide Tris-Tricine gels (Bio-Rad, Hercules, CA) and silver stained (20).
Transcriptional profiling. Cultures were grown in medium containing amikacin at a concentration of 50% of the respective MIC and were incubated until logarithmic phase. Volumes of the culture (10 ml) were collected and immediately immersed in liquid nitrogen for 10 s. The samples were then centrifuged (5,000 x g, 10 min, 4°C) and the supernatant was removed. The cell pellets were then once again immersed in liquid nitrogen for 15 s and then stored at 80°C for later RNA extraction.
RNA was extracted from strain B1 and its mutants by using a TRIzol Max bacterial RNA isolation kit (Invitrogen Life Technologies, Carlsbad, CA), followed by further purification with an RNeasy Mini kit (QIAGEN), according to the manufacturer's specifications. RNA from strain PAO1 was extracted by using just the RNeasy Mini kit, but the samples were further purified by using an additional RNeasy mini spin column (63).
Microarray probes were prepared according to the Affymetrix (Santa Clara, CA) P. aeruginosa GeneChip expression analysis protocol, with the slight modification that 2x PCR Enhancer Solution (Invitrogen Life Technologies) was added during cDNA synthesis. Target hybridization, washing, staining, and scanning were performed by the University of Colorado DNA Microarray Facility, according to the manufacturer's specifications, by using a GeneChip hybridization oven, a GeneChip fluidics station, a GeneArray scanner, and GeneChip operating software (v1.1) (Affymetrix).
Microarray data analysis. Microarray data were analyzed with ArrayAssist (Stratagene, La Jolla, CA). The data were analyzed by first performing robust multichip averaging (6, 32, 33), followed by principal component analysis (72), hierarchical clustering (18), and t tests. Data for a third replicate of strain B1 were removed from the analysis due to the large discrepancies in the results of principal component analysis and hierarchical clustering between it and the two other replicates. Additionally, the third replicate of B1 was generated by using a microarray from a different batch of Affymetrix GeneChips. Log fold differences are of base 2.
| RESULTS |
|---|
|
|
|---|
Initially, 850 mutants were screened by using replica plating for increased susceptibility to amikacin. Of the 850 mutants, 4 were confirmed to have MICs lower than that of B1 (amikacin MIC, 25 µg/ml) (Table 1). Three mutants (M5, M50, and M52) had MICs at intermediate levels of 8, 10, and 8 µg/ml, respectively, while mutant M31 had an MIC equal to that of laboratory strain PAO1 (4 µg/ml). It is of note that we performed the MIC assays with changes in concentrations less than the standard twofold change, with significance determined by t tests.
|
To gain an idea of the number of places on the genome that could be mutated to increase susceptibility to amikacin, we compared the frequencies for which we found sensitive mutants to the frequency of finding a mutant resistant to rifampin, which is known to be due to a single point mutation. The frequency of finding B1 mutants resistant to rifampin when they were subjected to the same level of exposure to MNNG was found to be 3.4 x 106, over 4 orders of magnitude less than the frequency of finding a sensitivity-increasing mutation.
Finally, we wanted to assess in greater detail whether or not sensitivity could be further increased through recursive mutagenesis. To do so, we subjected the most sensitive mutant resulting from a single round of mutagenesis, mutant M31, to a second round of random mutagenesis with MNNG. From this additional round, a total of close to 3,000 mutants were screened for increased susceptibility to amikacin. Interestingly, none were found to have further increased sensitivity relative to that of the M31 parental strain.
Transcriptional profiling. Transcriptional profiles were obtained for strain B1, the sensitive mutants M31 and M5, and strain PAO1 grown in the presence of half of their respective amikacin MICs. Amikacin was included in the medium so that the profiles of the strains would reflect the transcription of the strains under similar levels of antibiotic stress. In all cases, cells were cultured identically and samples were obtained in early exponential phase. To determine the extent to which sensitivity mutations altered the overall transcriptional profiles, we performed hierarchical clustering and principal component analysis with the gene expression data corresponding to each strain evaluated. Hierarchical clustering results (Fig. 1a) revealed that the sensitive mutants had expression patterns more similar to that of the laboratory strain PAO1 than to that of parental isolate B1. To further assess this unexpected result, we used principal component analysis (Fig. 1b) to reduce the large amount of data generated from the microarray experiments to a few key variables, called principal components, that account for much of the variation among samples (72). Again, the principal component analysis values corresponding to replicate microarrays of the sensitive mutants and PAO1 strain clustered tightly together, while those of strain B1 were separated from the other samples along the first and second principal components (which accounted for 54% and 23% of the overall variation, respectively). This indicated that, with respect to these principal components, PAO1 and the sensitive mutants were more similar to each other than to B1. The discrepancies observed between the B1 replicates could, in part, be explained by the difficulty of working with the mRNA of this isolate. This, however, does not detract from the finding that the sensitive mutants exhibited expression patterns more similar to that of PAO1 than to those of either of the B1 replicates. This result suggests that the major variation in the overall gene expression data is the result of differences between the B1 profiles and all other profiles.
|
0.05) and showed at least a 1 log fold (base 2) difference. A twofold change in expression (equivalent to 1 log fold, base 2) has previously been found to be sufficient for the detection of nearly 95% of differentially expressed genes by the use of Affymetrix GeneChips (11). The number of genes differentially expressed between M31 and B1 was also high (350 genes). Of these sets of genes, only 101 genes were found to be differentially expressed in both mutants. It should be noted that the changes in gene expression not only might be due directly to a mutation in that gene or a regulator of that gene but also might be due to the response of the overall genetic regulatory network to such mutations. Of the genes that were differentially expressed in mutant M5 or M31 compared with their expression in B1, a vast majority were genes that were grouped into three functional categories or that had no assigned function. The largest proportion of these genes were categorized as unclassified or as having only hypothetical functions (281 genes for M5 and 134 genes for M31). The largest functional category of genes (Fig. 2) comprised genes involved in cell permeability, LPS synthesis, efflux, and the transport of small molecules. There were changes in the expression levels in 113 of these genes for M5 and 71 for M31. The second and third categories contained genes involved in transcription or translation (88 for M5 and 32 for M31) and genes involved in metabolism, catabolism, and biosynthesis (92 for M5 and 40 for M31). Of particular interest, six genes encoding for aminoglycoside-modifying enzymes (AMEs) were found in both M5 and M31 to exhibit decreased expression relative to that in B1, suggesting that both strains would be similarly altered in aminoglycoside-modifying activity. Table 2 provides a detailed summary of the selected genes belonging to the AME or permeability functional class that had significant changes in expression.
|
|
Assays and MIC studies to further characterize sensitive mutants. Several hypotheses were developed on the basis of our transcriptional profile studies. First, the observed increase in the sensitivity of M5 and M31 could be due to the decreased expression of several AME genes relative to their levels of expression in strain B1. Second, the sensitivities of M5 and M31 could be the result of mutations affecting amikacin permeability. To test these hypotheses, we assessed changes in the MICs for the mutants strains, B1, and PAO1 (i) to multiple aminoglycosides to assess AME contributions, (ii) to amikacin in the presence of compounds that effect afflux pump activity, (iii) to the polycationic antibiotic polymyxin B to assess altered LPS mechanisms, and, finally, (iv) with spheroplasts of each strain to delineate between the mechanisms present within the plasma membrane and those due to the outer membrane and LPS structures (Table 1).
To examine changes in mutants M5 and M31 that may have been due to AMEs, MIC studies with four additional aminoglycosides were carried out (49). Only the MIC of amikacin for M5 was affected. Since genes affecting the accumulation of amikacin within the cell are less likely to be amikacin specific, the increased sensitivity of M5 is likely due to the decreased activity of an AME that has the capacity to modify amikacin. In contrast, M31 had decreased resistance not only to amikacin but also to tobramycin, gentamicin, kanamycin, and streptomycin. This would indicate that M31 has mutations that confer either increased accumulation of amikacin or, possibly, decreased AME activity.
To further examine the possibility of altered amikacin accumulation, the MIC of amikacin was determined in the presence of CCCP and verapamil. CCCP is an uncoupling proton ionophore that can carry protons across the cytoplasmic membrane and that can thus inhibit efflux pumps and other transporters that are dependent upon the proton motive force as a source of energy (47, 48, 69, 87). Verapamil, a calcium channel blocker, has been noted to inhibit efflux pumps, including ATP-binding cassette (ABC) efflux pumps (12, 35, 45). CCCP (250 µM) and verapamil (100 µg/ml) had no impact upon the MIC of amikacin for PAO1, M5, M31, or B1, indicating that efflux pumps are not substantially contributing to the changes in resistance levels.
Since efflux did not appear to contribute to resistance, we next sought to determine if altered uptake was affected. Polymyxin B is a polycationic antibiotic that is thought to exhibit self-promoted uptake mechanisms for entry into the cell membrane similar to those of aminoglycosides by first interacting with LPS molecules (62). Cross-resistance to polymyxin B and aminoglycosides has been noted for a number of systems (55, 70), possibly due to overall changes in the charge of the outer membrane (83) or the presence of outer membrane proteins that may be blocking antibiotic binding to the membrane (60, 61). The MICs for B1, M5, and M31 were the same as that for PAO1, indicating that changes in sensitivity levels were likely not due to significant changes that would hinder polycationic antibiotic binding to the LPS.
Finally, to delineate between resistance conferred by the outer membrane and resistance mechanisms within the cytoplasmic membrane, spheroplasts of B1, M31, M5, and PAO1 were made; and the MIC of amikacin was then determined. Interestingly, the MICs for spheroplasts of B1, M5, and PAO1 all decreased, while the MIC for the spheroplast for M31 was the same as that for unaltered M31 (Fig. 3). It should be noted that slight changes in MICs may not be clinically relevant but may be biologically relevant. Based upon replicates of this experiment, our results have been shown to be repeatable and statistically significant. The MIC for B1 spheroplasts decreased significantly, but not to the level of those for the M5 and M31 spheroplasts. The M5 spheroplasts had the same MIC as M31 and the M31 spheroplasts. This suggests that the resistance of B1 is due in large part to altered outer membrane permeability and, to a lesser extent, internal mechanisms or mechanisms related to the cytoplasmic membrane. M5 has an increased level of amikacin resistance compared with that for M31, and this increase is likely due to a mutation in permeability. Once the permeability barrier has been removed, the MIC drops to that for M31. The unchanged MICs for M31 and the M31 spheroplasts indicate that the susceptibility of M31 is due to mutations that affect cell permeability. The fact that the MICs for M5 and M31 spheroplasts did not drop to the levels for the PAO1 spheroplasts indicates that there are still mechanisms aside from the permeability of the outer membrane that are contributing to resistance. Furthermore, the internal resistance mechanisms of M5 and M31 are attenuated compared with those of B1.
|
| DISCUSSION |
|---|
|
|
|---|
A notable detriment to the use of chemical mutagenesis is the lack of an efficient way to determine the mutations that occur. One of the promises of transcriptional profiling is an improved ability to decipher the genetic basis of relevant phenotypes through comparisons of gene expression profiles. Thus, we next sought to determine if transcriptional profiling could be used to identify relevant genetic alterations in two sensitive mutants of the B1 isolate. We reasoned that if only minor changes were required for the restoration of sensitivity, then the transcriptional profiles of sensitive mutants should strongly correlate with the profile of the parental resistant strain and not with the profile of the laboratory PAO1 strain. We hypothesized that the profiles of sensitive mutants would differ in the expression of a limited number of genes, which might allow prediction of the sensitivity-restoring mechanism or mechanisms.
From hierarchical clustering and principal component analysis of gene expression data obtained from strain B1, mutants M5 and M31, and strain PAO1, we found that the sensitive mutants had significant changes in gene expression and that their gene expression patterns were more similar to that of PAO1 than to that of their parent isolate, B1. Upon examination of the genes that were differentially expressed in the sensitive mutants compared with their expression in B1, we found that (i) a large number of genes had altered expression in both mutants and (ii) the overall changes in gene expression were similarly distributed among different functional categories for both M5 and M31.
According to GeneChip analysis, several AMEs exhibited decreased expression in both M5 and M31. Among them were the genes for aminoglycoside 3'-N-acetyltransferases [aac(3)-Ib, aac(3)-IIIb, and aac(3)-IIIc], aminoglycoside adenyltransferases from integron cassettes (aadA6 and aadB), as well as an aminoglycoside-4'-adenyltransferase (ant(4')-IIa). Amikacin is believed to be protected by steric hindrance or folding from modification by aminoglycoside 3'-N-acetyltransferases (54). The aadA6 gene has been shown to render resistance to streptomycin and spectinomycin (57, 79), while aadB has been linked to gentamicin, kanamycin, and tobramycin resistance (56, 79). The AME gene with the largest difference in expression between B1 and the sensitive mutants was ant(4')-IIa, which encodes for an AME that has been shown to confer resistance to both tobramycin and amikacin (79). However, the level of amikacin resistance displayed by B1 is an order of magnitude lower than the MICs previously reported for other P. aeruginosa strains expressing this gene (34, 75, 78). Aminoglycoside-resistant isolates of serotype O12 have been reported to react with probes for the genes ant(3') and aac(6')-I and -II (46).
When attempting to confirm the presence of genes encoding for AMEs by PCR with AME-specific primers, we were able to amplify by PCR an aac(6)-Ib gene from strain B1 and mutants M5 and M31. This gene has been associated with amikacin resistance, with the reported levels of amikacin resistance of strains carrying this gene being similar to the levels found in strain B1 and the sensitive mutants (22). It is feasible that a decrease in expression due to a mutation in a regulator of this gene could cause a decrease in tolerance to amikacin or that a mutation in this gene could also affect its specificity for amikacin (41). Probes for this gene are not available on the Affymetrix P. aeruginosa GeneChip, which explains its absence from our transcriptional profiling studies. Based on this analysis and the high level of resistance to several aminoglycosides displayed by B1, we suspect that B1 does contain a set of AMEs that is partially responsible for its increased overall aminoglycoside tolerance but that these enzymes are not effective at modifying amikacin or are poorly expressed.
Transcriptional profiling also suggested that genes related to the cellular membrane or transport of small molecules were of interest. Of particular note was the change in expression of genes involved in O-antigen synthesis and assembly. The Affymetrix PAO1 GeneChip has probes for O-antigen genes from serotypes O6 (3) and O11 (17), in addition to PAO1 serotype O5 (8). Although strain B1 and its sensitive mutants are of serotype O12, genes involved in O-antigen biosynthesis (wbp genes) and assembly (wzx and wzz) from the O5, O6, and O11 serogroups showed 2- to 4-log-fold decreases in expression levels in sensitive mutants M5 and M31 compared with that in B1. The changes noted in genes from several serotypes may be explained by the fact that hybridization occurs for sequences with greater than approximately 70% identity. Notably, genes related to the B-band O antigen had significant changes in expression, suggesting that altered B-band synthesis could be responsible for the increased sensitivities of the M5 and M31 mutants. There are mixed reports of the impact that the loss of O antigens has upon resistance. The loss of the B-band O antigen has been shown to result in an increase in resistance levels to aminoglycosides in P. aeruginosa (7, 38), possibly because the B band is highly anionic and may thus attract and attach to the highly cationic aminoglycosides. On the other hand, the loss of O-specific antigen, and possibly other parts of the core region, has also been shown to decrease resistance to gentamicin, with the explanation that negatively charged sites of the lipid A of LPS may be involved in aminoglycoside uptake (76). Despite the changes in gene expression noted, sodium dodecyl sulfate-polyacrylamide gel electrophoresis analysis of LPS showed that all strains had smooth LPS, with no detectable differences in the LPS of B1 and the sensitive mutants.
Other genes that were identified as having significant changes in expression included oprF, oprG, and oprI, which encode outer membrane protein or lipoprotein precursors. These genes were all shown to have higher levels of expression in mutants M5 and M31 than in strain B1. OprF is a major porin of P. aeruginosa (25), and a decrease in expression of this porin has been observed in antibiotic-resistant strains of P. aeruginosa (37, 68) (the strains were not tested for aminoglycoside resistance). Additionally, M5 and M31 also have decreased expression of probable ABC-type transporters. Certain ABC transporters have been found to serve as efflux pumps for antibiotics and have been shown to increase the levels of resistance to a number of antibiotics (19, 31, 64, 67). Recently, inactivation of the glmR gene, which is believed to be involved in amino sugar metabolism, has been shown to increase sensitivity to aminoglycosides (71). This gene was found to have no significant change in expression between B1, M5, and M31.
Transcriptional profiling results suggested that the increased sensitivity of M5 and M31 could be due in part to modified AME activity or changes in amikacin permeation. We found that M5 had a change in its MIC for amikacin but not for other aminoglycosides, while M31 had a decreased tolerance for other aminoglycosides, in addition to amikacin. Additionally, spheroplasts of M5 had decreased tolerance to amikacin, while spheroplasts of M31 had no changes in amikacin resistance. We rationalized that mutations affecting AMEs would likely be more specific for different aminoglycosides, while mutations affecting permeability would likely be less specific and would possibly be indicated by changes in MICs between untreated cells and spheroplasts of cells. Thus, our results imply that M5 has mutations that affected the ability of AMEs to modify amikacin and that M31 has mutations that affect amikacin permeability and possibly mutations that affect AME activity.
Conclusion. This work describes an assessment of a reverse engineering approach for the identification and characterization of sensitive mutants of a resistant P. aeruginosa isolate. Our results indicate that the frequency of identification of sensitive mutants was unexpectedly high and was several orders of magnitude higher than the reported frequencies of finding resistant mutants. From the frequencies of detection of mutations, along with the various levels of amikacin resistance among the mutants that were identified, it was shown that there are multiple ways in which sensitivity could be increased. This result was bolstered by the dramatic and different changes in gene expression observed for each sensitive mutant, as well as the findings of detailed MIC assays performed with the two mutants examined in greater detail.
A key issue in our evaluation was whether or not whole-genome transcriptional profiling could be used to gain insight into sensitivity-restoring genes/pathways, which is the critical information required for the development of new therapies. Unexpectedly, we found that the changes leading to sensitivity resulted in dramatic but only partially coordinated alterations in gene expression between the two sensitive mutants. Furthermore, the transcriptional profiles of the sensitive mutants more closely resembled that of sensitive laboratory strain PAO1 than that of the resistant parental isolate. Analysis of such transcriptional profiles indicated that changes in amikacin permeability and/or modification by AMEs was the most likely source of the increased sensitivity in mutants M5 and M31, which was later confirmed by several additional assays. Although the transcriptional profiling results provided interesting general insights into sensitivity-restoring mechanisms, the surprisingly large number of genes displaying significant changes in expression prohibited the identification of any single target gene or mutation. Furthermore, this study was hindered by the lack of an array containing probes for all of the genes found within the resistant isolate but not found within the PAO1 genome. Future applications might benefit from the use of a combination of this approach with the use of conventional genetic strategies involving screening of a knockout library to study the impact that gene disruption may have upon sensitivity or screening of an overexpression library in which the effects of increased copy numbers of genes could be examined.
| ACKNOWLEDGMENTS |
|---|
We thank Mike Vasil (Department of Microbiology, University of Colorado Health Sciences Center) for supplying strains evaluated in this study. We also thank Shelley Copley (Department of Molecular, Cellular, and Developmental Biology, CUBoulder) and Mike Lynch for thoughtful suggestions throughout this work.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Clin. Vaccine Immunol. | Clin. Microbiol. Rev. |
|---|---|
| J. Clin. Microbiol. |