Harley H. McAdams,3 and
Gary L. Andersen1*
Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, California 94720,1 Frontier Research System for Extremophiles, Japan Marine Science & Technology Center, 2-15, Natushima-cho, Yokosuka 237-006, Japan,2 Department of Developmental Biology, Stanford University, Stanford, California 94305-53293
Received 24 June 2005/ Accepted 26 September 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Resistant bacteria possess a number of strategies to withstand elevated concentrations of heavy metals. Many resistance mechanisms revolve around removing the heavy metal or decreasing its toxicity (70). Alternatively, the concentration of metal entering the cytoplasm may be decreased through active (extracellular precipitation) and passive (native biosorption) processes (30, 35). Metal-chelating proteins have been reported as a means of resistance mainly in eukaryotes and also in some limited examples of prokaryotes (70). The major bacterial resistance mechanisms include (i) active efflux, (ii) transformation of the heavy metal ion to a less toxic form, for example, Cr(VI) can be reduced to Cr(III) (12), and (iii) precipitation, either intercellular or extracellular (35, 64, 66, 67).
Caulobacter spp. are extremely ubiquitous and are able to survive in low-nutrient environments (51). They have been found in freshwater, seawater, soil (51), ground water (37), wastewater (36), deep-sea sediment (38), and a deep subsurface gold mine (19) and have been noted for their ability to survive in broad environmental habitats where contamination may be present (8, 48). In addition, Caulobacter crescentus has been shown to form high-density biofilms with the potential for use in bioreactors for bioremediation (60) and has been used as a model organism to study cell cycle control (32, 40, 41). Previous knowledge of this organism, including its genome sequence (45), has provided new and extremely valuable tools to study genome-wide response to heavy metal stress. Both oligonucleotide cDNA microarrays and Affymetrix GeneChip microarrays have been used to study cell cycle regulation (31) (McGrath and McAdams, unpublished data). In this paper, we use a Caulobacter Affymetrix GeneChip microarray (Caulobacter chip) to study the transcriptional response of C. crescentus to heavy metal stress. This chip was designed by the McAdams laboratory at Stanford University in collaboration with Affymetrix. A complete description of all features of the chip will be published separately. In our work reported here, only the gene expression assay features, which are based on eight optimally selected 25-mer match/mismatch probe pairs per predicted open reading frame (ORF), were used.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Toxic metal effect on growth, survival, and morphology. Metal stock solutions were prepared by dissolving the compounds (Sigma-Aldrich) in water to 10,000 ppm, with the exception of the uranyl nitrate stock solution, which was dissolved to 100 mM (23,800 ppm). All metal stocks were sterilized by filtration through a 0.2-µm membrane. Overnight cultures were diluted in fresh M2G medium with various concentrations of cadmium sulfate (CdSO4), sodium selenite (Na2SeO3), potassium chromate (K2CrO4), potassium dichromate (K2Cr2O7), and uranyl nitrate [(UO2)(NO3)2 · 6H2O]. Growth was followed spectrophotometrically (optical density at 600 nm). One metal concentration was selected to be used for each set of microarray experiments, based on the following requirements for each metal. (i) The stressed conditions only slightly affected growth (increased doubling time by 15 to 30 min). (ii) The addition of the heavy metal compound did not result in precipitation with salts or cause other obvious changes in the medium. (iii) The metal concentration was above a level considered toxic or close to those conditions used commonly in other studies. While not all these criteria were fully satisfied, the final concentrations used in this study were a compromise of the factors described above and are reported in the results.
For morphology and motility observations, a sample was taken from mid-log-phase culture (OD600 was 0.3 to 0.4) at 0 hours after end log-phase growth (T0). Metal stock was added to the final required concentration to be used in microarray experiments. After 30 min, another sample was taken as T1. A final sample was taken 3 to 4 h poststress as T2. Ten microliters of each sample was examined under 100x phase-contrast light microscopy. Fifty microliters of each sample was used to determine CFU. Bacterial membrane integrity was assessed using a Live/Dead BacLight bacterial viability kit (Molecular Probes, OR) according to the manufacturer's recommended protocol.
RNA extraction. An overnight culture of C. crescentus CB15N was diluted into fresh M2G medium. When the culture reached exponential growth (OD600 was just over 0.3), 10 ml culture was removed (as nonstressed control). Heavy metal stock was added, and incubation was continued for another 30 min. After this period, a further 10 ml of sample was removed (as the stressed sample). Immediately after the samples were collected, they were centrifuged at 10,000 x g for 5 min and supernatant was removed. The cell pellets were frozen with liquid nitrogen and stored at 80°C. The RNA extraction protocol was described previously (18). Briefly, total RNA was extracted with Trizol reagent (Invitrogen) and any contaminating DNA was digested with DNase I. The RNA samples were further purified with acid phenol-chloroform-isoamyl alcohol (125:24:1, pH 4.5) (Ambion) extraction followed by salt-ethanol precipitation. RNA quantity was determined by OD260, and quality was determined by 2% (wt/vol) agarose gel electrophoresis and by a OD260/OD280 ratio.
Affymetrix GeneChip RNA expression analysis. Procedures for sample preparation and array processing are described fully in the Affymetrix GeneChip Expression Analysis Technical Manual and briefly described here. Transcripts of three genes (Bacillus subtilis dab, phe, and thr) were added to the total RNA as spike-in controls to monitor labeling, hybridization, and staining efficiency. To generate the spike-in control RNA, the plasmids containing B. subtilis phe, thr, and dap genes were purified from strains ATCC 87483, ATCC 87484, and ATCC 87486, respectively. Linear template DNA was generated by digesting the plasmid with restriction enzyme NotI and sense RNA produced subsequently by in vitro transcription using T3 RNA polymerase (MEGAscript T3 kit; Ambion).
Total RNA (12 µg) was primed with random primer (Invitrogen), and cDNA was synthesized with reverse transcriptase (superscript II, Invitrogen). The resulting cDNA was fragmented with DNase I (Amersham) and biotin labeled using the Enzo BioArray terminal labeling kit (Affymetrix). Biotin-labeled samples were hybridized onto the Caulobacter microarray at 50°C overnight, and chip washing and staining followed standard Affymetrix GeneChip protocols (with stringent washing at 50°C). The high-density chip was scanned using an Affymetrix Scan3000 scanner.
Microarray data analysis and identifying differentially expressed transcripts under heavy metal stress. The Caulobacter Affymetrix chip was used to assay gene expression levels for all 3,767 genes (45). For analysis of the protein-coding region, probe sets consisting of multiple (typically eight) 25-mer oligonucleotide probe pairs covering the gene were used for transcriptional interrogation. These probe sets were analyzed using the MAS5 statistical algorithm (3) for background adjustment and scaling in GCOS software (Affymetrix). Briefly, data from a minimum of three independent experiments were included as biological replicates in each comparison. Global scaling of all probe sets to a target signal intensity of 500 was applied to each chip (all microarray data are available at http://greengenes.lbl.gov/Download/Caulobacter_metal_stress_supplemental_microarray_data/). The data set was normalized using the spike-in controls mentioned above. For each comparison, the t test was performed on the data. For a probe set (gene) to be considered up-regulated under metal stress in these studies, it had to meet the following criteria. (i) The gene had to be called "present" by GCOS software for every experiment which was under metal stress. (ii) The average difference score (signal) for the gene had to be equal to or greater than 200, eliminating very low expression levels requiring more sophisticated analyses. (iii) The P value of the Student t test had to be less than 0.01, ensuring that the difference between the two conditions (nonstressed and stressed) was significant at a 99% confidence interval. (iv) The ratio of the average signals from stressed culture versus the average signals from nonstressed culture had to show at least twofold differences in expression.
The probes for predicted small protein regions were tiled at every 15 nucleotides on both strands. Since the boundary of these predicted proteins may be inaccurate, it is possible that some of the probes in the default probe set do not belong to a single transcript. It is inappropriate to use any of the existing software to obtain probe set values. At present, available microarray analysis programs do not allow dynamic definition of probe sets; thus, identification of transcripts driven by experimental data is not possible. To analyze these regions, we opted to examine their expression on a probe-by-probe basis. The background adjustment and normalization were performed using custom scripts developed by our lab. Probe-level data were examined in which individual probes meeting the following criteria were selected: (i) The average difference score (signal) for a probe was equal to or greater than 200. (ii) The P value of the Student t test was less than 0.015. (iii) The average signals from stressed culture versus the average signals from nonstressed culture were different by at least twofold to be considered for further analysis. Probes selected using these criteria were assembled in order to find continuous regions of up-regulation. The criteria defining a probe set for the small protein regions were as follows: (i) The maximum number of nonpassing probes (i.e., those not meeting the individual probe criteria listed above) within a probe set could not exceed three. (ii) The minimum length of assembled transcript was >100 bp. Probe set boundaries were also corrected according to probe behavior. For probe sets of interest, manual inspections were performed for final evaluation. This approach permitted monitoring of the transcriptional activities without bias.
Gene annotation and metabolic pathway analysis. Initial Caulobacter crescentus open reading frame annotations were taken from GenBank accession number NC_002696 (45). Clusters of orthologous gene descriptions (69) were used if they described the functions more clearly. If a clearly definitive annotation of function was not found, BLAST was performed against all bacterial genomes and Pfam analysis was run to identify any domains. The identification of pathways involved in metabolism was aided by the use of BioCyc (http://biocyc.org). Other analyses, such as operon prediction and gene ontology were improved through the use of the VIMSS (Virtual Institute for Microbial Stress and Survival) database (http://www.microbesonline.org/).
TEM observation. The samples were fixed in 4% paraformaldehyde in 0.1 M cacodylate buffer at pH 7.4 (Electron Microscopy Sciences) for 1 h at room temperature. The fixed samples were dehydrated with a graded series of ethanol and t-butyl alcohol. The samples were freeze-dried and mounted on 200-mesh copper grids for characterization with an HF2000 field-emission transmission electron microscope (Hitachi, Inc.) which was equipped with energy-dispersive X-ray spectroscopy (EDX). Transmission electron microscopy (TEM) was operated at an accelerating voltage of 200 kV. The abiotic sample was prepared by adding uranium stock solution to the sterile growth medium followed by centrifugation at 10,000 x g for 1 min.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
Differentially regulated genes common to multiple heavy metals. Transcriptional response to cadmium, chromate, dichromate, and uranium shared four up-regulated transcripts (Table 2); however, the only gene with known function was the superoxide dismutase with Mn as its cofactor (sodA). The greatest induction of this gene occurred under cadmium and chromium stress, with induction under uranium stress being lower. Superoxide dismutases are known to remove superoxide radicals which may be generated upon exposure to heavy metals (22, 23, 46, 63). While Caulobacter crescentus has three superoxide dismutase genes with different cofactors, sodA (CC1777, Mn), sodB (CC3557, Fe) and sodC (CC1579, Cu-Zn), the major gene involved in heavy metal response appears to be sodA. Under nonstressed conditions, it was at background levels but was up-regulated 19-, 14-, 9-, and 3-fold in cadmium, chromate, dichromate, and uranium stress, respectively. Transcription of sodB, on the other hand, was constitutive and was only twofold up-regulated under cadmium stress. Since sodB was induced under only cadmium stress, when sodA was already substantially up-regulated (19-fold), it is possible that sodB was up-regulated as a compensation mechanism under substantial oxidative stress. Transcription of sodC was also constitutive under nonstressed conditions but did not increase above the twofold threshold under any of the metal stress conditions. The difference in the regulation of the three enzymes may reflect their subtle differences in response to various oxidative stressors. It is known that C. crescentus sodA(Mn) and sodB(Fe) are cytosolic and sodC(Cu-Zn) is periplasmic (55). The induction of sodA and sodB but not sodC suggests that the oxidative stress imposed by the metals is intracellular rather than extracellular, which is consistent with the view that the oxidative stress originated from the reactive oxygen species generated by the interaction of metals and cellular components. The different scale of induction of sodA with various metals probably reflects the variable oxidative potentials of the metals inside the cells. Chromium and cadmium are well known for inducing oxidative stress in cells, while uranium, in our system, seemed to provoke a smaller response, although the precise intracellular quantity of metals is not directly measured in our system. As such, it is not possible to determine whether the variable oxidative response is a function of metal properties or whether it is a dose response.
|
We identified a total of 59 genes which were down-regulated at least fivefold under cadmium, chromium, or uranium stress, although we did not observe any genes commonly down-regulated more than fivefold. The most significant finding was the down-regulation of a sulfate transporter under chromium stress (see the supplemental data). Since chromate structurally resembles sulfate, and chromate typically enters the cell by the sulfate uptake system (47), we believe the response intended to reduce chromium entering into the cell. Similarly, dichromate stress provoked the down-regulation of the same transporter despite its structural differences. The reason may be rooted in chromate and dichromate equilibrium. Chromate and dichromate equilibrium are concentration and pH dependent (K = 102.2). Under our experimental conditions (pH 6 to 7 and concentration below 0.05 mM), a significant amount of chromate will be present (44) in the dichromate solution, thus eliciting a transcriptional response to the chromate stress.
Differential gene expression under cadmium stress. A total of 144 genes were up-regulated by at least twofold under cadmium stress (see the supplemental data). Several groups of annotated genes (Table 3) are discussed here.
|
Another cluster (cluster two) of efflux pumps, CC3195, CC3196, and CC3197 were also specifically up-regulated under cadmium stress; however, they were transcribed at basal but clearly above background levels under nonstressed conditions. The increase of these transcripts was moderate, about threefold under cadmium stress. It is possible that this group of proteins acted as "patrol" agents, responding to stresses if needed. On the other hand, proteins in cluster I were truly specifically induced by cadmium (expression was at background level under nonstressed conditions). It is not been experimentally demonstrated whether the two groups of transporters have subtle differences in functions. If the differential transcriptional patterns cannot be explained by functional requirement, this may be another example of supplementary induction similar to that observed with superoxide dismutase. In this scenario, it appears that the main system is specifically induced, and if it alone is not sufficient to manage the effects of intracellular stressors, a supplementary mechanism is also induced. At this point, it is not clear whether any of these efflux pumps work synergistically as essential components of one system or simultaneously to increase total output.
Cadmium is known to cause oxidative stress by depleting glutathione and protein-bound sulfhydryl groups, resulting in the production of reactive oxygen species. Consequently, it leads to enhanced lipid peroxidation, DNA damage, and altered calcium and sulfhydryl homeostasis (63). We observed many groups of genes which were up-regulated to deal with this type of stress. As described above ("Differentially regulated genes common to multiple heavy metals"), the enzyme directly involved in response to oxygen stress is the superoxide dismutase (Mn). Other proteins involved in removing toxic compounds and protecting thio groups were also up-regulated. Glutathione S-transferase has been shown to be induced by heavy metals in plants (39). Glutathione has two general functions: to remove toxic metabolites from the cell and to maintain cellular sulfhydryl groups in their reduced form. Glutathione S-transferase detoxifies xenobiotic compounds or products of oxidative stress by the covalent linking of glutathione to hydrophobic substrates (13, 17). C. crescentus glutathione S-transferases CC1316 and CC2434 were up-regulated in response to cadmium stress. Thioredoxin is a general protein disulfide reductase, believed to serve as a cellular antioxidant by reducing protein disulfide bonds produced by various oxidants. In this study, we noted three thioredoxin-coding transcripts (CC0062, CC0110, and CC3539) that were up-regulated under cadmium stress. Glutaredoxin is also known as thioltransferase; it is a small protein of approximately 100 amino acid residues which functions as an electron carrier in the glutathione-dependent synthesis of deoxyribonucleotides by the enzyme ribonucleotide reductase. Like thioredoxin, which functions in a similar way, glutaredoxin possesses an active center disulphide bond and exists in either a reduced or an oxidized form where the two cysteine residues are linked in an intramolecular disulphide bond. It is not surprising that C. crescentus glutaredoxin (CC2505) was up-regulated under cadmium stress, while the other glutaredoxin (CC0829) was not, since CC0829 has more sequence homology to E. coli glutaredoxin 3 (grx3), which has been shown to have narrower substrate specificity in vivo (54). Peptide methionine sulfoxide reductase (CC0994 and CC1039) was also up-regulated, and, together with thioredoxin, it can reverse the effects of oxidative damage on methionine residues in proteins (34, 42). Under cadmium stress, both were up-regulated slightly over twofold.
Both the riboflavin biosynthesis pathway and GTP cyclohydrolase I (folE) were up-regulated (3.5- to 4-fold) under cadmium stress (Table 3). Their up-regulation may be due to oxidative stress. Previous studies have shown that H2O2-induced expression of GTP cyclohydrolase I mRNA in vascular endothelial cells (56) is extremely sensitive to oxidative stress and is also one of the major targets of H2O2 in E. coli (33). There are several possibilities for the induction of the riboflavin biosynthesis pathway; for example, studies in E. coli have demonstrated that GTP cyclohydrolase II is induced by redox-cycling agents and is positively regulated by soxR and soxS (29), the global regulators for oxidative stress. Riboflavin is a precursor of both flavin mononucleotide (FMN) and flavin adenine dinucleotide, which are important coenzymes of several oxidoreductases (1, 52). In addition, ribonucleotide reductase was up-regulated, and, in E. coli at least, this enzyme is activated by flavins (15). It provides precursors necessary for DNA synthesis by catalyzing the reductive synthesis of deoxyribonucleotides from their corresponding ribonucleotides. This reaction also replenishes the pool of reduced thioredoxin, which also has a role in protection against oxidative stress (54). Another example of an up-regulated enzyme with FAD and FMN as cofactors is CC2129, which belongs to the NADH:flavin oxidoreductase family, members of which are capable of reducing a range of alternative electron acceptors. It is also possible that GTP cyclohydrolase II was up-regulated to protect genetic material against damage from oxygen radicals. For example E. coli GTP cyclohydrolase II can hydrolyze 8-oxo-dGTP, which is an oxidized form of dGTP; overproduction of this enzyme has been shown to reduce mutation frequency of the mutT strain (mutT protein prevents A:T to C:G transversion by hydrolyzing 8-oxo-dGTP) to almost normal levels (28).
Interestingly four genes resembling arsenic resistance system components were induced by cadmium (Table 3). These genes have homologues in Pseudomonas putida KT2440 (11), in Thiobacillus ferrooxidans (10), and in virulence plasmid pYV of Yersinia enterocolitica (43). CC1503 is most similar to arsC, an arsenic reductase that is responsible for arsenate reduction to arsenite which is then translocated out of the cell by arsB using proton motive force. However, while CC1504 would be the primary candidate for arsB according to its location and the fact that it also contains a transmembrane domain, its sequence and the domain structure are different from those of arsB. As such, its capacity to function in a similar manner to arsB or indeed its involvement in response to cadmium stress is not clear. CC1505 is most similar to arsR, which is a transcriptional repressor with a helix-turn-helix DNA-binding domain that is thought to dissociate from DNA in the presence of metal ions (74). Sequence comparison shows that the metal-binding site of arsR resembles that of cadC from Staphylococcus aureas plasmid pI258 (76) and from Bacillus firmus OF4 (7). Cadmium ion is thought to relieve repression by cadC (20, 74). Therefore, it is possible that the arsRBC operon in C. crescentus was up-regulated because cadmium binds to arsR and releases the repression. CC1506 is most similar to arsH and was highly up-regulated in C. crescentus cells under cadmium stress. Despite the fact that CC1506 contains a NADPH-dependent FMN reductase domain, its direct involvement in metal resistance is not definitive. Likewise the detailed function of other arsH genes is unclear.
Response to uranium stress. Bacteria are known to possess several mechanisms for resistance to uranium that frequently involve precipitation to reduce toxicity. For example, uranyl ions may be sequestered intracellularly by complexation with phosphate granules, as in the case of Arthrobacter spp (66) and Pseudomonas spp (64). In Citrobacter spp., inorganic phosphate is liberated from organic forms, resulting in precipitation of various uranyl phosphate crystal complexes outside the cell (35). To investigate if such mechanisms were utilized as a uranium resistance mechanism by Caulobacter crescentus, we performed TEM and EDX analysis and demonstrated that C. crescentus did not form any uranium-containing phosphate granules intracellularly. However, TEM images of whole cells of C. crescentus revealed extracellular precipitates associated with the cells (Fig. 2A). EDX spectra from cells and extracellular precipitates show that while uranium is almost absent within cells (Fig. 2B), extracellular precipitates contain high concentrations of uranium, phosphorus, and calcium (Fig. 2C), suggesting that the extracellular precipitates are composed mainly of these elements. The selected-area electron diffraction patterns of the extracellular precipitates indicated that the extracellular precipitates are amorphous, although uranyl phosphate minerals are known to readily become amorphous under the high vacuum of TEM (68). In order to examine whether the formation of the extracellular precipitates was catalyzed by C. crescentus or was simply an abiotic phenomenon, uranium-bearing precipitates formed in the same medium in the absence of C. crescentus were also characterized (Fig. 2D). EDX analysis revealed that abiotic precipitates were mainly composed of uranium and phosphorus with substantially lower calcium content than the biological precipitates (Fig. 2E), indicating either direct or indirect biological involvement in the formation of the calcium-uranium-phosphate complex.
|
Uranium stress also induced transcription of 48 transcripts which were up-regulated at least twofold (see the supplemental data); however, the response to uranium does not appear to overlap substantially with other heavy metal stresses evaluated in this study. It is believed that uranium, unlike cadmium and chromate, imposes less direct oxidative damage to cells, thus it is not surprising that most of the commonly up-regulated genes noted under cadmium and chromium stress were not up-regulated under uranium stress.
A large portion of the up-regulated genes were difficult to classify into pathways due to mainly a lack of functional annotation (Table 4), although even with clearly annotated proteins, such as two-component systems, their targets have not previously been identified. Two up-regulated two-component systems, CC1293-CC1294 and CC1304-CC1305, were identified. They are DNA-binding response regulators consisting of a CheY-like receiver domain and were up-regulated (seven- to ninefold) under only uranium stress. Both CC1293 and CC1304 have signal-sensing domains at the N-terminal, and a DNA-binding domain at the C-terminal but no sigma factor interaction domains were detected. CC1294 and CC1305 have histidine kinase domains at the C-terminal, and a HAMP membrane domain was also detected at the N-terminal of CC1294, indicating it could be a membrane protein. The E value of Pfam analysis for a membrane domain in CC1305 is too low (0.27, compared to 6.6e12 for CC1294) to believe it has a membrane component; therefore, it may be cytoplasmic. The domain structure of the proteins suggests that CC1294 and CC1305 receive their respective metabolic/environmental signals in the following manner. The histidine kinase catalyzes ATP and transfers a phosphoryl group to the response regulator (CC1293 or CC1304), resulting in activation of the DNA-binding domain that elicits the specific responseactivation or repression of the transcription of their targeted gene(s).
|
CC1293 and
CC1294 are knockout mutants in which CC1293 and CC1294 were replaced in-frame by tetracycline-resistant cassettes (59). Both mutants were generated using strain CB15N, thus the results of the mutants should be directly comparable with other data in this study. Sequencing of PCR amplification products confirmed the deletion and correct chromosomal location; furthermore, transcription of CC1293 or CC1294 was not detected in their respective mutants. Compared to that of the wild type, the growth of
CC1294 was unaffected under 200 µM uranyl nitrate stress; however, growth did slow after 120 min of uranium stress at 1 mM. Growth of
CC1293 showed no significant difference from the wild type at any uranium concentration tested. Expression of CC1293 decreased in the
CC1294 background, possibly because CC1294 is closer to the transcriptional start of the operon and replacing it with a tetracycline cassette negatively impacted the overall quantity of the CC1293 message; despite this, the regulation was not lost (Table 5). With the exception of four transcripts (Table 5) whose functions are unknown, transcripts responding to uranium maintained the same regulation patterns in the knockout mutants as in the wild type. Therefore it appears from growth and microarray data that CC1293-CC1294 is not a master regulator of uranium response despite its specific up-regulation under uranium stress.
|
|
|
In C. crescentus, it appears that the majority of up-regulated genes were in response to chromium-induced oxidative stress (Table 7). Both chromate and dichromate contain Cr(VI) and, subsequently, a large portion of the up-regulated genes (214) are common to both metal stresses; however, their chemical structures differ and this may account for the number of genes up-regulated only twofold upon exposure to one form of chromium. Following chromate exposure, 84 genes were specifically up-regulated, while 23 were specifically up-regulated following dichromate exposure (see the supplemental data).
|
DNA damage by reactive oxygen species upon chromium exposure is well documented (4, 73). Up-regulation of recA is known to be induced by DNA breakage in E. coli, and previous studies have shown that chromate-induced DNA damage strongly depends on the reactive intermediates. Frequently, chromate causes DNA single-strand breakage and 8-hydroxydeoxyguanosine formation (correlated with hydroxyl radical as the DNA-damaging species) (5, 26). In this study, DNA repair enzymes, such as CC2272 and CC2200, were up-regulated (Table 7). This suggests that chromate-induced damage in C. crescentus cells may indeed be mediated by hydroxyl radicals generated through nonspecific chromate reduction.
Studies of animal and plant cells have shown that chromium can cause membrane damage through direct or oxidative stress-mediated interactions (14, 25, 50), and our transcriptional data indicate that chromium exposure induces a membrane response. We observed induction of two OmpA family proteins [ompA mutants exhibit increased sensitivity to environmental stresses (72)], TonB-dependent outer membrane receptors, and lipopolysaccharide biosynthesis (Table 7). Our data indicate that the TonB receptor family of proteins was involved in the response to several different metals, such as chromium (Table 7), cadmium (Table 3), and uranium (Table 4), and yet, in most cases, different proteins in the family were up-regulated under specific metal stresses. However, it is not clear whether the up-regulated TonB receptor gene interacts with TonB protein, since the expression of the predicted TonB protein (CC2327) was suppressed and not activated under heavy metal stress. It is likely that the receptor merely binds to the substrates (heavy metals) and communicates extracellular environmental information.
We observed up-regulation of several genes and pathways which typically occur during stationary phase or under nutrient-limiting conditions (Table 7), including acquisition of ammonium, phosphate-starvation response, poly ß-hydroxybutyrate (PHB) biosynthesis (including genes in the fatty acid oxidation pathway which can provide precursors for synthesis of medium chain PHBs), and energy/carbon utilization (PHB depolymerase and glucan glucohydrolase).
The cells appeared to exhibit an increased demand for ammonium and glutamate, since ammonium transport, the glutamate synthase, and the histidine degradation pathway (generating ammonium and glutamate) were up-regulated (Table 7). Glutamate is one of the central amino acids that links nitrogen (glutamine synthetase-glutamine 2-oxoglutarate-aminotransferase cycle) and carbon/energy metabolism (trichloracetic acid cycle) through
-keto-glutarate. There appeared to be an increased flow of glutamate to
-keto-glutarate, as indicated by the up-regulation of the serine biosynthesis pathway. In addition, glutamate synthase contains an iron-sulfur cluster and has close structural homology with glutamine phosphoribosylpyrophosphate amidotransferase (71). This latter enzyme from B. subtilis is known to be inactivated by O2 in stationary phase (9). Although up-regulation of ntrBC and PII usually indicates nitrogen limitation response, the key enzyme (glutamine synthetase) was not up-regulated. It is possible that increased ammonium uptake was the first step towards responding to nitrogen demand, and thus it seems plausible that simply obtaining more ammonium from the environment would provide a more energy efficient response than provoking a complete nitrogen starvation response. If this hypothesis is correct, our data suggest that the three PII proteins of C. crescentus may be regulated separately by specific physiological conditions and that ntrBC may be involved in activating the transcription of not only glutamine synthetase but also other operons, such as ammonium transport. Therefore, it may be that nitrogen regulation in C. crescentus is more complicated than in E. coli, at least at the transcriptional level.
Differential gene expression under selenium stress. The response to sodium selenite was mild compared to other metal stresses. Only 12 transcripts were up-regulated and at most by four- to fivefold (see the supplemental data). All of these were also observed to be up-regulated under chromium or cadmium stresses and included membrane components, glutathione S-transferases, and transport proteins.
Conclusions. In this study, we investigated whole-genome transcriptional response of Caulobacter crescentus to the stress of several heavy metals, including chromium and uranium, which are significant environmental contaminants and a current focus of bioremediation efforts. In addition to the surprising finding that C. crescentus CB15N is tolerant to uranium, our studies combining physiology observation, transcriptional measurement, and imaging analysis clearly showed that Caulobacter formed a calcium-uranium-phosphate precipitate extracellularly, in contrast to the intracellular sequestration mechanism of other resistant bacteria, such as Arthrobacter spp. This observation was consistent with a limited response to oxidative stress, such as that seen with other metals. The stress response strategy of lowering intracellular metal concentration was also present in cadmium and chromium response. Efflux pumps were up-regulated under cadmium stress. C. crescentus does not seem to have a specific extrusion mechanism for chromium; however, the cells down-regulated a sulfate transporter, which may reduce the uptake of chromate.
In broader terms, cells exposed to cadmium share many up-regulated transcripts with those under chromium stress. Most of those up-regulated genes respond to oxidative stress, such as superoxide dismutase, glutathione S-transferase, thioredoxin, and glutaredoxin. However, on closer inspection, the individual proteins up-regulated, and the fold changes were specific to each metal (for example, different sets of glutathione S-transferase were up-regulated under cadmium and chromium stress), indicating the subtle difference of each metal stress and physiological conditions. We also observed up-regulation of TonB-dependent outer membrane receptors which may serve as sensors for environmental signals. While the detailed mechanisms of their involvement are still not known, our results suggest that they may be an important member of the response network.
We believe that not all of the observed up- or down-regulation was a direct response to the metal toxicity. This was particularly evident in the case of chromium and/or cadmium stress. While the direct involvement of an arsenic reductase operon in cadmium response is unclear, the likely mechanism is that cadmium binds to the repressor of the operon, resulting in its up-regulation. The response of C. crescentus under chromium stress was clearly different from that of other oxidative stress (for example, cadmium), yet it may be complicated by secondary responses. Future transcriptomics studies, with various concentrations of chromium, augmented with proteomic analyses may help elucidate the complex response observed to this heavy metal such as the role of cytochrome oxidases and the apparent nutrient-limitation response.
Our data have also clearly demonstrated the importance of interrogating the whole genome on both strands. We have identified antisense transcripts which are differentially regulated specific to metals, which, as either proteins or RNAs, may play an import part in the response model.
| ACKNOWLEDGMENTS |
|---|
We thank Susan Lynch and A. C. Matin at Stanford University for the generous gift of bacterial strains Escherichia coli K12 and Pseudomonas putida KT2440. We thank M. Laub at Harvard University for the generous gift of C. crescentus two-component signal transduction deletion mutants
CC1293 and
CC1294.
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
Present address: National Institute of Advanced Industrial Science and Technology (AIST), Research Center for Deep Geological Environments, AIST Tsukuba Central 7, 1-1-1 Higashi, Tsukuba, Ibaraki 305-8567, Japan. ![]()
| REFERENCES |
|---|
|
|
|---|
siki, and R. J. Reiter. 2001. Carcinogen-induced, free radical-mediated reduction in microsomal membrane fluidity: reversal by indole-3-propionic acid. J. Bioenerg. Biomembr. 33:73-79.[CrossRef][Medline]