Department of Molecular Biology and Biotechnology, Krebs Institute for Biomolecular Research, The University of Sheffield, Sheffield, United Kingdom1
Received 22 June 2004/ Accepted 27 September 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Although Zn is classified as essential, excess Zn is toxic. The diverse toxic effects have been attributed to the interactions with sulfydryl groups in a wide range of proteins, especially those proteins involved with respiration in electron transport systems (22). In Escherichia coli, membrane uptake systems, including the primary importer ZnuABC, a high-affinity periplasmic binding protein-dependent transport system for Zn (32), and ZupT (zinc uptake transporter) (16), a member of the ZIP (ZRT, IRT-like protein) family, provide adequate levels of Zn, even in the face of variable extracellular concentrations. Intracellular pools must be regulated, in part by sequestration. However, information on the E. coli proteins ZraP (zinc resistance-associated protein) (31) and YdaE (5) that may process or chaperone Zn for intracellular functions is rudimentary. In other bacteria, metallothioneins play a role in detoxification of Zn (34). The levels of Zn are also maintained by the export of excess ions from the cell by both ZntA, a P-type Zn-exporting ATPase (3), and members of the CDF (cation diffusion facilitator) family of proteins, including ZitB (formerly YbgR), a metal ion transporter (15).
Genome-wide transcriptional analysis with microarray technology is a powerful tool for determining how gene expression within the cell is affected by the external growth conditions. To date, genome-wide microarray studies on the effect of zinc on E. coli have been carried out only with batch cultures (8) by using zinc-adapted cells. The problem with these cultures is that the growth conditions cannot be strictly monitored and maintained at constant levels as the pH and the dissolved oxygen and nutrient concentrations change during growth. The regulation of gene expression is therefore affected not only by the difference in growth conditions between a control and an experimental culture but also by culture variables, including, perhaps most importantly, the specific growth rate. This variation in gene regulation can obscure the interpretation of gene expression in response to a supposedly single difference in culture growth conditions.
In this study, E. coli cultures were grown in chemostats at a fixed specific growth rate, temperature, and pH so that the effects on the transcriptional profile in the presence and absence of additional ZnSO4 in the medium could be compared by microarray analysis. We found significant changes in gene transcription in response to elevated levels of Zn. The transcriptional response to zinc stress involves not only established zinc transport systems but also membrane transporters and sensors previously not linked to zinc metabolism.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Growth in liquid batch cultures was measured in sidearm flasks by using a Klett colorimeter and a red filter (number 66; Manostat Corporation). The colorimeter was zeroed by using minimal medium.
ZnSO4 chemostat experiments and RNA isolation. Cells were grown in two parallel chemostats under batch conditions overnight. ZnSO4 (final concentration, 0.2 mM) was then added to one medium reservoir before continuous culturing began. Cells in both chemostats were grown under the conditions described above to the steady state, as verified by collecting at least 5 culture volumes and measuring the optical density at 600 nm with a Jenway 6100 spectrophotometer (path length, 1 cm) before cells were collected by centrifugation. Samples were taken for subsequent microarray analysis by harvesting 10 ml of culture directly into RNA Protect (QIAGEN) to stabilize the RNA. Total RNA was purified by using a QIAGEN RNeasy Mini kit as recommended by the supplier. The RNA concentration and purity were determined by using a Beckman DU 650 spectrophotometer (37).
Preparation of labeled cDNA and hybridization. Equal quantities of RNA from control and ZnSO4-supplemented cultures were labeled by using nucleotide analogues of dCTP containing either the Cy3 or Cy5 fluorescent dye. For each microarray slide, one sample was labeled with Cy3-dCTP, while another sample was labeled with Cy5-dCTP. Dye swap experiments were performed for each pair of samples to compensate for the different efficiencies of incorporation of the labeled nucleotides (41). The slides used were E. coli K-12 Pan arrays purchased from MWG Biotech. These slides contained 4,288 gene-specific oligonucleotide probes representing the complete E. coli K-12 genome.
cDNA synthesis (with MWG Biotech's protocol) was carried out by using 12 µg of starting material, which was primed with 9 µg of pd(N)6 random hexamers (Amersham Biosciences). Reaction mixtures (20 µl) containing 0.5 mM dATP, 0.5 mM dTTP, 0.5 mM dGTP, 0.2 mM dCTP, and 1 mM Cy3- or 1 mM Cy5-dCTP were incubated for 2 h at 42°C with 200 U of Superscript II RNase H reverse transcriptase (Invitrogen). Following synthesis, cDNA was purified by using a PCR purification kit (QIAGEN) to remove the unincorporated deoxynucleoside triphosphates, fluorescent dye, and primers. Equal volumes of cDNA were combined and evaporated for approximately 45 min with an SPD121P Speed Vac (Thermo Savant). For hybridization to the microarray slides, cDNA was resuspended in an appropriate volume of salt-based hybridization buffer (provided by MWG Biotech). Prior to addition to the slides, cDNA samples were heated to 95°C for 3 min. The slides were placed in MWG Biotech hybridization chambers and incubated for 16 to 24 h in a water bath at 42°C.
Washing and scanning of slides. Following incubation, the slides were washed with decreasing salt concentrations, as recommended by the supplier, ranging from 2x SSC-1% sodium dodecyl sulfate to 0.5x SSC at 37°C with gentle agitation (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate). The slides were dried by centrifugation at 1,200 rpm (254 x g) for 5 min and subsequently scanned with an Affymetrix 428 scanner.
Data analysis.
The average signal intensity and local background correction were obtained by using commercially available software from Biodiscovery, Inc. (Imagene, version 4.0, and GeneSight, version 3.5). The mean values from each channel were log2 transformed and normalized by using the subtract-by-mean method to remove intensity-dependent effects in the log2 (ratio) values. The Cy3/Cy5 fluorescence ratios were calculated from the normalized values. Biological experiments (i.e., chemostat growth with and without 0.2 mM ZnSO4) were carried out three times, and a dye swap analysis was performed for each experiment, which provided three technical repeats, one for each of the three biological repeats. Data from the independent experiments were combined. Genes that were differentially regulated
twofold and displayed a P value of
0.05 (as determined by a t test) were defined as being statistically significantly differentially transcribed. This experiment (raw data only) has been loaded into ArrayExpress with accession number E-MEXP-222 (http://www.ebi.ac.uk).
Real-time PCR. RNA was extracted as described above. cDNA synthesis was carried out by using 4 µg of starting material primed with 9 µg of pd(N)6 random hexamers (Amersham Biosciences). Reaction mixtures (20 µl) containing 0.5 mM dATP, 0.5 mM dTTP, 0.5 mM dGTP, and 0.5 mM dCTP were incubated for 2 h at 42°C with 200 U of Superscript II RNase H reverse transcriptase (Invitrogen). Following synthesis, cDNA was purified by using a PCR purification kit (QIAGEN) to remove unincorporated deoxynucleoside triphosphates and primers.
Gene-specific primers were designed to amplify 50- to 150-nucleotide fragments of target genes by using the Primer 3 software (36). Each reaction was carried out in a 25-µl (total volume) mixture in a 96-well optical reaction plate (Applied Biosystems). Each well contained 50% SYBR green PCR master mixture (AB Applied Biosystems), 12.5 pmol of each of the two primers, and 5 µl of cDNA sample. PCR amplification was carried out by using an ABI 7700 thermocycler (PE Applied Biosystems) and the following thermal cycling conditions: 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. The data were analyzed by using the SEQUENCE DETECTOR SYSTEM software (PE Applied Biosystems) and were further processed with Microsoft EXCEL. A standard curve was established by using genomic DNA for each gene studied to confirm that the primers amplified at the same rate and to validate the experiment. The relative levels of expression of genes of interest in zinc-grown E. coli compared with untreated controls were calculated by using the protocol for the standard curve method in User Bulletin #2 (ABI Prism 7700 sequence detection system) supplied by Applied Biosystems. No-template reactions were included as negative controls.
| RESULTS |
|---|
|
|
|---|
|
Genes up-regulated by ZnSO4. Table 2 lists all 64 genes that were significantly up-regulated along with the functions of their products and their levels of regulation. Encouragingly, several genes that we expected to see up-regulated because they are known from other studies to be induced by zinc were found to be highly responsive. These included zntA, whose product is the primary zinc ATPase transport pump (3), zraP, which encodes a possible zinc chaperone as it is known to bind to zinc in the periplasm (25), and hydG, which is one component of a two-component response regulator for zraP (25). More interestingly, there were many genes involved in transport, sensing, or regulation that have not previously been implicated in zinc tolerance but were highly up-regulated in response to zinc in this study.
|
|
acrAB mutant (2, 30). It is interesting that MdtAC alone confers bile salt resistance but not novobiocin resistance, suggesting that MdtB extends the substrate specificity of the transporter (30). mdtD was also up-regulated; although this gene is a member of the same operon, it is not involved in novobiocin and deoxycholate resistance, but it does encode a transporter of the major facilitator type. Other up-regulated transporters include acrD (6.2-fold increase), encoding an aminoglycoside efflux pump. It has been suggested that AcrD can become associated with a member of the membrane fusion protein family and an outer membrane channel, encoded by genes elsewhere on the chromosome, to function as part of a multisubunit transporter complex (35). Also, amtB was induced (2.5-fold increase); this gene encodes an ammonium-methylammonium transporter and is linked with another up-regulated gene, glnK (4.1-fold increase), encoding a small cytoplasmic signal transduction protein (39).
Three genes encoding proteins involved with flagella were significantly up-regulated. These genes were fliM, which is involved in flagellar synthesis, motor switching, and energizing (4.1-fold increase), flgB, which encodes the basal body rod subunit (3.2-fold increase), and motB, whose product is involved in flagellar rotation (2.3-fold increase). These genes were not known to be regulated by zinc, but fliC has previously been shown to be up-regulated by Al (17).
The requirement in zinc tolerance for BasR-BasS, a two-component sensor regulation system, is unclear. Both genes were significantly up-regulated in this study (Fig. 2C), but to date little is known about this system in E. coli. The equivalent system in Salmonella, PmrA-PmrB, has been better studied and is required for resistance to polymyxin B and other microbial compounds. This two-component system controls transcription of several loci (Fig. 3), including pmrF (pbgPE operon) and pmrG, whose products modify lipopolysaccharide (LPS) and therefore are essential for polymyxin B resistance. Transcription of pmrCAB is also controlled, and the system is autoregulated (24). The up-regulation of basR (3.4-fold increase) and basS (2.4-fold increase) was confirmed by the up-regulation of numerous genes, all of which are E. coli equivalents of the Salmonella pmrF, pmrE, and pmrCAB genes. Unlike the Salmonella pmrF operon, which contains seven genes, the gene set in E. coli consists of six genes, all of which were shown to be up-regulated in this study (Fig. 2D) and some of which are known to be involved in lipopolysaccharide modification in E. coli. The first gene, arnB (6.2-fold increase), encodes an aminotransferase that is involved in aminoarabinose modification of lipid A (6). arnA (2.2-fold increase) encodes a UDP-glucuronate dehydrogenase that converts UDP-glucuronic acid to UDP-Ara4O (uridine 5'-ß-L-threo-pentapyranosyl-4"-ulose diphosphate) (7), and arnT (1.9-fold up-regulated) (data not shown) encodes a 4-amino-4-deoxy-L-arabinose:lipid A transferase that utilizes the novel glycolipid donor undecaprenyl phosphate-L-4-amino-4-deoxy-L-arabinose (40). The other three genes in the predicted operon, arnC (11-fold increase), yfbH (4.5-fold increase), and yfbJ (2.0-fold increase), were all up-regulated, but at present their functions are unknown. Like Salmonella pmrC, the homologue in E. coli, eptA was also up-regulated (3.5-fold increase) and encodes a lipid A phosphoethanolamine transferase. The homolog of pmrE in E. coli, ugd, was up-regulated (5.0-fold increase) and is believed to encode a UDP-glucose 6-dehydrogenase.
|
Interestingly, up-regulation in most of the operons or groups of adjacent genes identified resulted in similar patterns of expression. For the majority of genes in operons, all genes showed up-regulation, although the induction ratios decreased for the downstream genes (Fig. 2).
In addition to the various genes mentioned above, there are many "y" genes of unknown function that either are unclassified or have unknown functions.
Genes down-regulated by ZnSO4. Table 3 lists all 58 significantly down-regulated genes along with the functions of their products and their levels of regulation. The percentage of genes whose functions are unclassified or unknown that were down-regulated was larger than the percentage of genes whose functions are unclassified or unknown that were up-regulated, and the most down-regulated gene was yrbL (7.7-fold decrease). Although cells in the chemostat were constantly growing at a steady state and growth was limited by a single nutrient (glycerol in this case), some down-regulated genes are known to be induced under stationary-phase conditions; these genes include csiE (2.7-fold decrease), osmE (2.9-fold decrease), osmC (2.6-fold decrease; involved in the detoxification of organic hydroperoxides [1]), uspB (2.1-fold decrease; required for stationary-phase resistance to ethanol [9]), and gadE (formerly yhiE) (2.7-fold decrease; encodes a protein that activates glutamate decarboxylase-dependent acid resistance [26]).
|
-aminobutyrate, and gadC (5.9-fold decrease) is predicted to encode a putative glutamate-
-aminobutyrate antiporter (28). It has been proposed that these two proteins function together to help maintain the intracellular pH when cells are exposed to extremely acidic conditions (38). The expression of gadBC, however, requires gadE, which is also down-regulated (26). hdeA (5.4-fold decrease) encodes a periplasmic chaperone of acid-denatured proteins, and hdeB encodes a structural homologue of HdeA, which may form heterodimers with HdeA in order to prevent aggregation of the periplasmic proteins that become denatured in acidic conditions (13). hdeD (2.8-fold decrease) was also repressed; this gene encodes a putative membrane transporter that is suggested to be involved in acid stress resistance (27). Other transport systems shown to be down-regulated include csgF (6.2-fold decrease), csgG (2.7-fold decrease), ompW (2.4-fold decrease) encoding an outer membrane protein, and lsrB (2.5-fold decrease), which is part of a system for autoinducer 2 uptake in Salmonella. Also, ugpA (2.2-fold decrease) and ugpE (2.4-fold decrease) were down-regulated; these genes encode hydrophobic integral inner membrane proteins that form a substrate-translocating complex for the uptake of sn-glycerol-3-phosphate (20).
Validation of the microarray gene profiling by quantitative real-time PCR (RT-PCR). The fact that our analysis identified a number of genes (e.g., zntA, zraP, hydG, and spy) previously known to be induced (directly or indirectly) by Zn(II) in E. coli validates the ability of the microarray experiments to identify candidate genes important in resistance to toxic concentrations of extracellular Zn(II).
To verify further the results obtained by the microarray analysis, several genes that exhibited induction with 0.2 mM ZnSO4 were examined by quantitative real-time PCR to independently determine relative mRNA levels. The levels of up-regulation determined by RT-PCR (mean ± standard deviation) were as follows; cusF, 89 ± 0.7; basR, 4.6 ± 0.2; ais, 19 ± 0.02; asr, 38 ± 0.03; and cpxP, 14 ± 0.3. These values correspond to increases in the microarray analysis of 48-, 3-, 16-, 25-, and 11-fold, respectively. The mRNA levels from asnS (internal control) were unchanged as determined by both RT-PCR and the array analysis. The values obtained by RT-PCR strongly support the fidelity of the cDNA microarray analysis, taking into account any variation expected due to the differences between the two methodologies.
Hypersensitivity to ZnSO4 of strains deficient in genes up-regulated in the microarray experiments. In this study we identified many genes that were up-regulated by zinc, and some of these genes (notably zntA) are known to be important for conferring resistance to high zinc concentrations. To assess the importance of other up-regulated genes, a selection of knockout mutants with mutations in genes shown to be up-regulated in the microarray work were obtained. These mutants included strains lacking basR, basS, fliM, ycdM, and yibD. Triplicate liquid culture growth curves (Fig. 4) revealed that all five mutants were hypersensitive to ZnSO4 when their growth was compared to wild-type growth. In GGM containing no additional ZnSO4 (Fig. 4A and B) there was no difference in the growth between any of the mutant strains and the wild type. When 0.15 mM (final concentration) ZnSO4 was added (data not shown), slight hypersensitivity was revealed by a small reduction in growth for the basR and basS mutants compared to the other strains. Addition of 0.175 mM (final concentration) ZnSO4 (Fig. 4C and D) to the medium resulted in a large reduction in growth (particularly for the basS mutant strain) and a much longer lag phase for the basR and basS mutants. Mutants with mutations in ycdM and yibD also showed hypersensitivity to ZnSO4, and there was an extended lag phase of growth similar to that observed for the basR and basS mutants. Finally, when ZnSO4 was added to a final concentration of 0.2 mM (Fig. 4E and F), all mutants showed hypersensitivity. All mutants except the fliM mutant exhibited extended lag phases during growth and reductions in the amount of growth. The fliM mutant had an extended lag phase that was shorter than those of the other mutants but showed no reduction in the amount of growth. To eliminate the possibility that a Tn5 insertion anywhere in the chromosome does not result in Zn hypersensitivity, growth curves for metA::Tn5 and qor::Tn5 mutant strains in liquid culture were determined. Both strains showed no heightened sensitivity (compared to wild-type strain MG1655) to a range of ZnSO4 concentrations up to 0.2 mM.
|
|
| DISCUSSION |
|---|
|
|
|---|
Previous work on metal ion tolerance in E. coli by Brocklehurst and Morby (8) involved microarray analysis of preadapted wild-type cells grown (under batch conditions) to the stationary phase. The work was carried out in LB medium in which the exact amounts of bioavailable metal ions were unknown. The use of a minimal medium, particularly a medium lacking inorganic phosphate that might form insoluble products with metal ions, allows maximum availability of added metal ions (3, 21). We did not observe any of the genes that were shown to be zinc responsive in the study of Brocklehurst and Morby (8), nor did the experiments of these authors reveal any genes that are known through other studies to be responsive to zinc. However, Brocklehurst and Morby (8) did find insA and insB genes that were up-regulated, which was observed with many different stresses, not just with stress due to excess zinc. In this work, the dilution rate, nutrient limitation, and pH were all controlled and/or monitored and were shown to be equal for both control and zinc-treated cultures. In agreement with Hayes et al. (19), we consider the use of chemostat cultures to be highly desirable in microarray analysis for identifying a gene expression change as a result of a single change in the growth conditions.
The microarray analysis demonstrated that the mdtABC operon was up-regulated in response to stress caused by excess zinc. This RND-type efflux system has been implicated in conferring resistance to certain antibiotics, including novobiocin and the bile salt component deoxycholate (2). Our finding that mdt is up-regulated in response to zinc therefore has potentially important implications for the influence of metal stresses on bacterial resistance to antibiotics. The presumed transport function of Mdt might indicate that it transports zinc directly, resulting in net export of the metal. An alternative possibility (of many) is that the transported solute moderates zinc toxicity, perhaps by binding the metal. The regulation of mdt is thought to be achieved by the two-component BaeRS system, encoded by genes immediately downstream of mdtABC. The signals sensed by BaeRS are not known. However, baeRS and ultimately mdt are clearly up-regulated by zinc, but it is possible that mdt transcription is not dependent solely on BaeRS.
Zinc supplementation in chemostat cultures was sufficient to increase expression of basR and basS (pmrA and pmrB, respectively, in Salmonella). We propose that high extracellular concentrations of zinc ions may be sensed by BasS (Fig. 3B) in a way that is similar to way that high iron concentrations are sensed by PmrB in Salmonella (Fig. 3A). This sensing of high zinc concentrations by BasS could in turn (as it does in Salmonella) cause phosphorylation of BasR and ultimately up-regulation of numerous genes (Fig. 3B). The up-regulation of these genes is an important finding, since this two-component system regulates, for example, resistance to polymyxin B and cationic peptides in Pseudomonas aeruginosa and Salmonella and also controls virulence in Erwinia carotovora (23). Interestingly, the importance of zinc in extracellular polysaccharide synthesis has recently been independently described by Hagiwara et al. (18). The Rcs signaling system in E. coli that is involved in capsular polysaccharide synthesis has recently been shown to be zinc responsive. Propagation of the zinc-responsive Rcs system was also largely dependent on PhoPQ, while in Salmonella the same two-component regulator system is implicated in the PmrAB mechanism. Thus, there is a recent confluence of data that strongly point to links among the sensing of and tolerance to divalent metal cations, including Zn and Mg(II), extracellular polysaccharide synthesis, the sensitivity of LPS to metal ions, and ultimately the sensitivity to antibiotics. Teleologically, a reason for basRS up-regulation by zinc ions might be related to the requirement for this system in modifying LPS, which is stabilized by divalent cations. Indeed, for yfb and arn (pbg in Salmonella), whose products are required for the modification of lipid A, there was a dramatic increase in expression, which was induced by zinc in these microarray experiments. Millimolar concentrations of magnesium stabilize the outer membrane (as judged by decreased susceptibility of periplasmic proteins to loss), while very high concentrations of calcium (20 mM) and magnesium (100 mM) increase the outer membrane permeability (and are commonly used, of course, in preparing cells that are competent for transformation by exogenous DNA). We speculate that the role of increased basRS transcription is to increase synthesis or effect modification of the LPS either to stabilize it in the face of high zinc concentrations or to act as a biosorbent of free metal ions.
The results of the hypersensitivity experiments with E. coli strains lacking some of the up-regulated genes identified, including basR, basS, yibD, and ycdM, support the hypothesis that cells grown with excess zinc need these up-regulated genes to respond to the stress. These mutant strains all showed hypersensitivity to zinc compared to the wild-type strain, with extended lag phases and lower growth yields, and it can therefore be suggested that these genes are important for zinc homeostasis within the cell. On the basis of present knowledge, it is not clear why ycdM and yibD (or fliM [Fig. 4]) appear to have roles in Zn tolerance. However, microarray analyses provide investigators with clues to the roles of genes with unknown functions or new roles for previously studied genes.
In summary, our results endorse the power of the microarray approach for identifying new genes required for zinc tolerance. Many genes have been identified as genes that are up-regulated in response to zinc, some of which were expected; however, the vast majority of the genes are not currently known to be involved in zinc homeostasis in E. coli. These findings may have important implications for understanding not only zinc homeostasis within the cell but also how bacterial resistance to antibiotics is modulated by metal ions.
| ACKNOWLEDGMENTS |
|---|
We thank Julie Scholes (Department of Plant and Animal Sciences, University of Sheffield) for her help and for use of resources for the RT-PCR experiments. We also thank F. Blattner and the E. coli Genome Project, University of Wisconsin-Madison, for kindly providing the MG1655 mutant strains and Tony Gordon and Pam Trickett for assistance with microarray profiling. We are grateful to Martin Hughes (King's College, London, United Kingdom) for his critical reading of the manuscript.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
s-regulated gene in Escherichia coli which is required for stationary-phase resistance to ethanol. J. Bacteriol. 180:6140-6147.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Appl. Environ. Microbiol. | Infect. Immun. | Eukaryot. Cell |
|---|---|---|
| Mol. Cell. Biol. | J. Virol. | Microbiol. Mol. Biol. Rev. |
| ALL ASM JOURNALS |