|Home | About | Journals | Submit | Contact Us | Français|
The biological impact of 24-h (“chronic”) chromium(VI) [Cr(VI) or chromate] exposure on Shewanella oneidensis MR-1 was assessed by analyzing cellular morphology as well as genome-wide differential gene and protein expression profiles. Cells challenged aerobically with an initial chromate concentration of 0.3 mM in complex growth medium were compared to untreated control cells grown in the absence of chromate. At the 24-h time point at which cells were harvested for transcriptome and proteome analyses, no residual Cr(VI) was detected in the culture supernatant, thus suggesting the complete uptake and/or reduction of this metal by cells. In contrast to the untreated control cells, Cr(VI)-exposed cells formed apparently aseptate, nonmotile filaments that tended to aggregate. Transcriptome profiling and mass spectrometry-based proteomic characterization revealed that the principal molecular response to 24-h Cr(VI) exposure was the induction of prophage-related genes and their encoded products as well as a number of functionally undefined hypothetical genes that were located within the integrated phage regions of the MR-1 genome. In addition, genes with annotated functions in DNA metabolism, cell division, biosynthesis and degradation of the murein (peptidoglycan) sacculus, membrane response, and general environmental stress protection were upregulated, while genes encoding chemotaxis, motility, and transport/binding proteins were largely repressed under conditions of 24-h chromate treatment.
The metal oxyanion chromate (CrO42−) is a widespread environmental contaminant due to its prevalent use in industrial and defense applications such as tanning, electroplating, paint pigment manufacturing, stainless steel welding, and nuclear weapons production (25, 26). The hexavalent form of chromium, Cr(VI), is highly soluble and toxic, with chronic exposure leading to mutagenesis and carcinogenesis. Cr(VI)-induced apoptosis, for example, was demonstrated in p53 human bronchoalveolar cells (46), and Cr(VI) exposure results in a spectrum of genomic damage in cultured cells including DNA single-strand and double-strand breaks, binding of amino acids and proteins to DNA, DNA interstrand cross-links, and Cr-DNA adducts (11, 27, 28, 43, 47, 53, 59, 64, 65). Cr toxicity is also associated with the generation of reactive oxygen intermediates during the intracellular partial reduction of Cr(VI) to the unstable intermediate Cr(V) by various in vivo nonspecific reductants (e.g., glutathione, NADH, NADPH, and cysteine) or cellular one-electron reductases (16, 27, 50). The other most stable, common form of chromium, trivalent Cr(III), is considered less toxic than Cr(VI) because of its tendency to form insoluble hydrated Cr3+ complexes, which cannot cross cell membranes. However, Cr(III) was shown to cause DNA damage and inhibit topoisomerase DNA relaxation activity in bacteria (40).
The adverse biological impact of Cr(VI) is attributable to the cellular uptake process. Chromate is transported across eukaryotic and prokaryotic cellular membranes via surface anion transport systems, namely, the sulfate transport system (16, 36, 38). Microorganisms have evolved diverse resistance mechanisms to cope with chromate toxicity. These detoxification strategies include biosorption, diminished intracellular accumulation through either direct obstruction of the ion uptake system or active chromate efflux, precipitation, and reduction of Cr(VI) to the less toxic, less mobile Cr(III) (reviewed in reference 16). Plasmid-determined resistance to chromate, for example, has been shown to occur in bacteria, including species of Pseudomonas (9, 17, 49) and Alcaligenes (36). A hydrophobic protein with 12 proposed transmembrane-spanning domains, designated ChrA, was found to be responsible for the plasmid-specified resistance phenotype in these organisms (16, 18, 37) and appears to function as a secondary transport system for the extrusion of chromium ions (3).
The in situ microbial catalysis of Cr(VI) reduction to sparingly soluble, less bioavailable Cr(III) has been proposed as a potential remediation strategy for Cr(VI)-contaminated subsurface environments. Shewanella oneidensis MR-1, a facultatively anaerobic γ-proteobacterium, possesses diverse metal-reducing capabilities, including the ability to transform Cr(VI) to Cr(III) (34, 56). As a result, its potential utility in the bioremediation of dissolved metal ions prompted the complete sequencing of the S. oneidensis MR-1 genome (23). Predicting the utility of S. oneidensis MR-1 for remediating metal-contaminated sites requires an understanding of the gene/protein components and cellular pathways enabling heavy metal resistance and biotransformation. In addition, the efficacy of in situ chromate bioremediation will depend on the capacity of remediating bacteria to cope with and perhaps minimize the cellular effects of heavy metal toxicity. Knowledge of the molecular and physiological response of MR-1 to Cr(VI) toxicity, in particular prolonged or chronic exposures, remains limited, however.
The primary goal of this work was to gain insight into global changes in the mRNA and protein expression patterns that occur in S. oneidensis MR-1 cells at a time point (in this case, 24 h) marked by the complete removal of chromate from the culture medium. We refer to this prolonged exposure to Cr(VI) and its derivatives as a “chronic” challenge to distinguish it from acute exposures before the onset of Cr(VI) reduction. Both transcriptome profiling and whole-cell proteomic analysis revealed that the predominant molecular response to chronic chromate exposure in complex medium was the induction of prophage-related genes and their encoded products as well as a number of functionally undefined hypothetical genes clustered within the lambda- and mu-like integrated phage regions of the MR-1 genome. In addition, a number of genes and/or their corresponding proteins with annotated functions in DNA metabolism (topB, hsdM-2, uvrD, and recO), electron transport (NADH dehydrogenase gene cluster, SO3056, and SO4360), cell division (era, ftsAZ, and ftsL), biosynthesis of murein (peptidoglycan) sacculus (rodA, mreD, murE, murF, murG, and mraY), and protein fate and general environmental stress responses (htpG, SO3391, hslVU, ibpA, groES, and groEL) were induced, while genes encoding regulatory, chemotaxis, motility, and transport/binding proteins were generally repressed. The accompanying physiological response to chronic Cr(VI) exposure included marked changes in cellular morphology as revealed by scanning confocal microscopy and atomic force microscopy (AFM). This study combines cell imaging, transcriptional measurement, and proteomic characterization to provide insight into the response and susceptibility of S. oneidensis MR-1 to 24-h chromate exposure.
The S. oneidensis wild-type strain MR-1 (35) was cultivated aerobically at 30°C in Luria-Bertani (LB) medium (pH 7.2) or in LB medium amended with 0.3 mM potassium chromate (K2CrO4; Sigma-Aldrich Co., St. Louis, MO) with shaking at 250 rpm in a New Brunswick Scientific (Edison, NJ) incubator. A starter culture was prepared by inoculating LB medium amended with 0.3 mM chromate with a single colony of S. oneidensis MR-1 from a freshly streaked LB plate and allowing the cells to grow in the presence of chromate for 24 h. An aliquot of this starter culture, giving an initial optical density at 600 nm (OD600) of 0.08, was used to inoculate fresh LB medium containing either 0 or 0.3 mM potassium chromate. OD600 was monitored in triplicate using either a Spectronic 20D+ spectrophotometer (Thermo Electron Cooperation, Waltham, MA) or the Bioscreen C microbiological culture system (Growth Curves USA, Piscataway, NJ) as described elsewhere (13). SigmaPlot version 8.0 (SPSS Inc., Chicago, IL) was used to plot growth curves.
Extracellular Cr(VI) was quantified spectrophotometrically at a wavelength of 540 nm using the 1,5-diphenylcarbazide (DPC) method as described elsewhere (39). Cultures were assayed at different time points (0, 3, 6, 9, 12, 15, 18, 21, and 24 h) during growth to determine the amount of residual Cr(VI) remaining in the medium by measuring absorbance at 540 nm using a Varian (Cary-1E) UV-visible spectrophotometer (Hewlett-Packard, Wilmington, DE). Uninoculated LB medium containing 0.3 mM chromate served as the abiotic control.
S. oneidensis MR-1 cells were grown aerobically for 24 h in LB medium in the presence or absence of 0.3 mM chromate. Cells were then harvested and stained using the LIVE/DEAD BacLight bacterial viability kit (Molecular Probes, Eugene, OR) in accordance with the manufacturer's instructions. Following a wash in distilled water, cells were filtered onto Isopore membrane filters (Millipore Corporation, Bedford, MA). The stained cells were viewed by confocal laser scanning microscopy using a TCS SP2 microscope (Leica Microsystems, Inc., Exton, PA). Z-series optical sections of cells were taken at 0.5 μm spacing and processed using Leica confocal software.
Following 24-h exposure in LB medium alone (control sample) or LB medium containing 0.3 mM chromate (experimental sample), MR-1 cells were fixed in 4% paraformaldehyde for 16 to 18 h at ambient temperature. After fixation, cells were pelleted by centrifugation at 3,000 rpm for 3 min and then resuspended in distilled water. A drop of the resuspended cells was placed on a polylysine-coated microscope slide and left undisturbed for 15 min. Slides were washed briefly with distilled water, dried, and mounted for imaging. AFM images were recorded using a PicoPlus atomic force microscope (Molecular Imaging Inc., Tempe, AZ) with a 100-μm scanning head. The instrument was operated in MacMode using silicon cantilevers (type II Maclevers; Molecular Imaging, Tempe, AZ) with a spring constant of 2.8 N/m at a scan speed of 0.8 lines/second at 512 pixels per line scan. All of the images presented are first order flattened.
For transcriptome characterization, a starter culture, which was preexposed to 0.3 mM chromate for 24 h, was used to inoculate 100 ml of LB medium only (control samples) or LB medium containing 0.3 mM chromate (experimental samples) in 250-ml sidearm Pyrex flasks. A total of three separate control cultures and three separate experimental cultures were permitted to grow aerobically in parallel at 30°C for 24 h, at which time cells were harvested for total RNA extraction by centrifugation at maximum speed in a 5415R centrifuge (Eppendorf, Westbury, N.Y.) for 30 s at 4°C and washed once in ice-cold LB medium to remove residual chromate. RNA isolation, fluorescein-labeled cDNA preparation, and probe purification were carried out essentially as described previously (13).
An S. oneidensis MR-1 whole-genome open reading frame array containing approximately 95% of the total predicted MR-1 gene content was used to examine the global transcriptional response to 24-h chromate exposure. Construction of this array has been described previously (13, 20). Global gene expression profiling was performed using six independent microarray experiments (three biological replicates and two dye reversal reactions) as described in our previous study (14). Microarray prehybridization, hybridization, and posthybridization washings were carried out as described previously (13). The arrays were scanned as described elsewhere (52). Image quantification, data normalization, and analysis of gene expression data for statistical significance have been described previously (13). Transcripts exhibiting a statistically significant change in expression (P < 0.05) and a twofold or greater change in magnitude were considered for further analysis (48). The complete microarray data set for this study is available as supplemental Table S1 and can be accessed online(http://compbio.ornl.gov/shewanella_metal_stress/chronic).
Reverse transcriptase, quantitative real-time PCR (RT-qPCR) was used to provide an independent assessment of gene expression for selected MR-1 genes. Six selected genes (SO0404, SO2426, SO2823, SO2945, SO3056, and SO3585) found to be differentially expressed by microarray hybridization were examined by comparative RT-qPCR as described elsewhere (13, 61). These genes are predicted to encode two hypothetical proteins (SO0404 and SO2945), two signal transduction response regulators (SO2426 and SO2823), a tetraheme cytochrome c (SO3056), and a putative azoreductase (SO3585). Relative expression patterns for each selected gene were independently confirmed using the following primer pairs, which were designed using the program Primer3 (http://www-genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi): SO0404, 5′-AGTATAACCAAGCGCCAGTA and 3′-GCATCGGTATTAACTTGCTC; SO2426, 5′-GCAGAAGGATTTAGGTCGAT and 3′-GCCACAACTAATTTCATGCG; SO2823, 5′-CGACACCTTACCTAAAATCG and 3′-GGCATCTATGGGTTTGAGTA; SO2945, 5′-CTGACTTGCGGGATAAATAC and 3′-GGTAAGTGACGAACATCCAT; SO3056, 5′-GATGCAGTAGTTGTCATGGA and 3′-CAGACTATCGTGGGGATTAG; SO3585, 5′-CGAGGCTATCCATCACTTAG and 3′-TGGAAAACACGATAAAGACC.
All chemical reagents were obtained from Sigma Chemical Co. (St. Louis, MO) unless stated otherwise. Modified sequencing grade trypsin (Promega, Madison, WI) was used for all digestions. High-pressure LC (HPLC) grade water and acetonitrile (ACN) were acquired from Burdick & Jackson (Muskegon, MI), and 99% formic acid was purchased from EM Science (Darmstadt, Germany).
S. oneidensis MR-1 cells analyzed by mass spectrometry were cultivated as described above for the transcriptome characterization. Following 24-h growth in the absence or presence of 0.3 mM chromate, S. oneidensis cells were lysed using the method described by Brown et al. (14). Samples were separated into crude and membrane protein fractions using high-speed ultracentrifugation (100,000 × g for 60 min at 4°C) and were quantified using bicinchoninic acidanalysis (Pierce, Rockford, IL). Approximately 2 mg of protein from each proteome fraction for the two different growth conditions was digested with sequencing grade trypsin at 1:100 (wt/wt) as described previously (14). Samples were immediately desalted using Sep-Pak Plus C18 solid-phase extraction (Waters, Milford, MA). Using a Savant SpeedVac (Thermo Electron Corporation, Waltham, MA), samples were concentrated and solvent exchanged into 0.1% formic acid in water to ~10 μg/μl starting material and then filtered, aliquoted, and frozen at −80°C until ready for LC-MS/MS analysis.
The experimental (treated with chromate for 24 h) and control (untreated) proteome samples were analyzed in duplicate via two-dimensional (2-D) LC-MS/MS analysis with an Ultimate HPLC (LC Packings, a division of Dionex, San Francisco, CA) connected to an LTQ linear trapping quadrupole (Thermo Finnigan, San Jose, CA). The HPLC pump was set at a flow rate of ~100 μl/min, which was split precolumn to achieve a final flow rate of ~200 nl/min at the nanospray tip. A split-phase column was packed as described previously (14) with approximately 500 μg of each sample loaded onto the column via a pressure cell. This loaded column was then placed behind a PicoFrit tip (inner diameter [I.D.], 100 μm; 15-μm I.D. at the tip; New Objective, Woburn, MA) packed in-house via a pressure cell with ~15 cm of the reverse phase (Aqua C18 5-μm 200A Phenomenex). The columns were positioned in front of the LTQ on a nanospray source (Thermo Finnigan).
All control and experimental samples were subjected to a 24-h 12-step 2-D analysis with an increasing amount (0 to 500 mM) of ammonium acetate via salt pulses accompanied by a 2-h 100% aqueous solvent (95% H2O-5% ACN-0.1% formic acid)-to-50% organic solvent (30% H2O-70% ACN-0.1% formic acid) gradient (method adapted from reference 29). The LTQ was run in the data-dependent MS/MS mode during the entire chromatographic procedure. Chromatographic methods and HPLC columns were identical for all sample analyses. The LC-MS/MS system was fully automated and under direct control of the Xcalibur software system (Thermo Finnigan).
The S. oneidensis MR-1 protein database used in MS/MS spectrum searches consisted of the latest genome annotation (http://www.tigr.org/) along with 36 common contaminating proteins (e.g., trypsin and keratin from the sample preparation procedure); the database is available for download at the project website http://compbio.ornl.gov/shewanella_metal_stress/chronic/databases/. MS/MS spectra were searched against the MR-1 database using SEQUEST (Thermo Electron) with cross-correlation scores of at least 1.8 (+1), 2.5 (+2), and 3.5 (+3), as described previously (14). Output data files were filtered at the one- and two-peptide levels, sorted, and compared as detailed by Brown et al. (14). Supplemental Table S2 (available online at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental) contains a list of all proteins identified at the two-peptide level organized according to the percentage of the protein sequence identified (percent sequence coverage), the number of peptides identified for the protein (peptide count), and the number of spectra that were confidently identified for the proteins (spectral count). Differentially expressed proteins under chromate conditions were identified using a semiquantitation method based on the following criteria: greater-than-40% sequence coverage, five or more unique peptides, and/or a twofold difference in the number of mass spectra identified under the control conditions versus the experimental conditions (14, 55, 61).
Our previous work investigating the molecular response of S. oneidensis MR-1 to an acute chromate challenge (a maximum exposure time of 90 min) indicated a MIC of chromate (K2CrO4) of approximately 2 mM under aerobic conditions in LB medium (14). Growth studies were conducted to determine whether preexposure to a sublethal dose of chromate, in this case 0.3 mM, increased cellular tolerance to this heavy metal. Preexposed cells attained a slightly higher OD600 at the 24-h time point compared to cells that had not been preexposed to Cr(VI) (Fig. (Fig.1).1). Similar to our previous observations (14), cells exhibited a biphasic OD pattern in the presence of 0.3 mM Cr(VI), regardless of whether cells had been preexposed (Fig. (Fig.1).1). Cr(VI) disappearance continued to occur during periods where the OD of chromate-treated cultures either changed very little (resembling a stationary-phase plateau) or actually declined, while the second phase of OD increase (at ~18 h) coincided with the appearance of cells displaying extremely elongated or filamentous morphology (data not shown for the 18-h time point). Hence, the increases in OD600 observed after about 18 h may not be an accurate reflection of cell growth but may be due largely to the drastic changes in cell size. This was shown to be the case for Escherichia coli K-12 cells, which also exhibited biphasic growth kinetics in chromate-amended LB medium (2).
The impact of Cr(VI) exposure on the cell morphology and surface topology of S. oneidensis cells was examined by confocal laser scanning microscopy and AFM. The results indicated that cells from 24-h cultures with no detectable extracellular Cr(VI) formed long, apparently aseptate filaments in contrast to the normal rod-shaped morphology exhibited by untreated cells (Fig. (Fig.2).2). As shown in Fig. Fig.2C,2C, the majority (approximately 93%) of examined control cells (grown in LB medium only) appeared normal in size (~1 to 4 μm in length) and were motile, while only a small percentage (7%) fell within the size range of 5 to 7 μm. However, only about 58% of the cells subjected to 24-h chromate exposure had sizes within the normal range (i.e., ~1 to 4 μm in length), while 35% were ~8 μm or longer (in some cases ~25 μm) (Fig. (Fig.2C).2C). Cells that were longer than normal were also nonmotile, and some appeared as aggregated filaments. These filamentous cells, in which septation was apparently impeded, were observed under conditions of complete Cr(VI) disappearance in the growth medium as indicated by the colorimetric DPC assay and not under acute chromate stress conditions characterized by short exposure times (e.g., 90 min) and no detectable Cr(VI) transformation (results not shown). This marked effect of chromate on cell morphology has been observed for other bacteria. Nonadapted E. coli K-12 cells exposed to 250 μM potassium chromate (K2CrO4) exhibited extreme filamentous morphology within 3 h of Cr(VI) challenge (2). An enlarged or elongated cell size also has been observed following exposure to other stress conditions, such as high salt (57), low-temperature growth at 3°C (1), and UV radiation (42), and likely constitutes a general cellular response to environmental stress. The increased cell size of Staphylococcus aureus under high NaCl was correlated with shorter interpeptide bridges, less cross-linked peptidoglycan, and retarded cell separation (57).
Consistent with the bright-field imaging results, S. oneidensis filamentation (Fig. (Fig.3B)3B) was observed using AFM following a 24-h exposure to 0.3 mM chromate, which contrasted with the normal rod-shaped morphologies of untreated control cells. In some cases, damage to the cell surface was severe enough to lead to the apparent disintegration of the outer membrane and hence lysis of the bacterium (Fig. (Fig.3B),3B), revealing the underlying cytoplasmic contents. AFM indicated a similar breach in the cell membrane integrity of Escherichia coli cells subjected to different concentrations of the β-lactam antibiotic cefodizime (10).
Microarray profiling and 2-D LC-MS/MS were performed to analyze the transcriptome and proteome, respectively, of S. oneidensis MR-1 at the 24-h time point following initial exposure to 0.3 mM chromate. Gene and protein expression profiles of cells exposed to chromate were compared to those of untreated control cells grown in parallel (see Materials and Methods for experimental details). Approximately 26% (n = 1,215) of the total predicted S. oneidensis MR-1 genes (n = 4,648) with representative probes on the array displayed at least a twofold statistically significant (P < 0.05) change in transcript expression at the 24-h time point following chromate exposure. RT-qPCR was used to provide an independent assessment of mRNA expression for targeted genes analyzed in the transcriptome experiment. The six genes (changes[n-fold] given in parentheses) selected for RT-qPCR analysis exhibited either upregulated or downregulated expression patterns as identified by microarray hybridization: SO0404 (14.4), SO2426 (0.45), SO2823 (13.4), SO2945 (271), SO3056 (8.5), and SO3585 (0.7). Comparison of gene expression measurements determined by RT-qPCR and microarray analysis indicated that the two different data sets were highly correlated, with a Pearson correlation coefficient, r, of 0.95.
At the protein level, a total of 2,313 gene products, representing 47% of the predicted MR-1 proteome, were identified at the two-peptide level in duplicate analyses of control samples and samples exposed to chromate for 24 h (Table (Table1).1). A total of 3,051 proteins were identified using the less stringent one-peptide filter level; however, a one-peptide filter level for the identification of proteins results in a dramatically higher false-positive discovery rate (44), so a thorough analysis of only the two-peptide data is presented here. The levels of reproducibility between replicate analyses of the control proteome and the experimental (chromate-treated) proteome on the LTQ instrument were 75.6% and 77.2%, respectively. For the observed proteome identified in this study, 109 protein species were found to be differentially expressed under prolonged Cr(VI) exposure, with 56 proteins displaying increased abundance and 53 showing decreased abundance (Table (Table1).1). Supplementary transcriptome and proteome data (i.e., a list of all significantly expressed genes, the complete raw and filtered proteome data, and a list of differentially expressed proteins) can be accessed online at http://compbio.ornl.gov/shewanella_metal_stress/chronic/.
Proteins identified at the two-peptide level under the two different growth conditions were organized in Table Table11 according to the functional categories assigned by The Institute for Genomic Research (see www.tigr.org; Comprehensive Microbial Resource). Using multidimensional HPLC-MS/MS, we identified more than 75% of the predicted MR-1 proteins for the following five functional categories: amino acid biosynthesis; biosynthesis of cofactors, prosthetic groups, and carriers; protein synthesis; purines, pyrimidines, nucleosides, and nucleotides; and transcription. Under 24-h Cr(VI) exposure, more proteins were identified in the functional categories of amino acid biosynthesis; biosynthesis of cofactors, prosthetic groups, and carriers; purines, pyrimidines, nucleosides, and nucleotides; and protein synthesis than were found under the corresponding control condition (Table (Table11).
Sequence analysis of the MR-1 genome revealed the presence of an integrated lambda-like phage (LambdaSo; 51,857 bp) and two phylogenetically distinct phages related to the E. coli mu (MuSo1 [34,551 bp] and MuSo2 [35,666 bp]) (23). The lambda-like phage genome is also present in MR-1 in a nonintegrated form (23). There are 75, 42, and 53 open reading frames (ORFs) annotated as LambdaSo, MuSo1, and MuSo2 genes, respectively (23).
A previous study focusing on S. oneidensis MR-1 demonstrated the induction of a large number of prophage-related genes in response to UV radiation, particularly those genes from the integrated LambdaSo genome, and the presence of phage particles in UV-irradiated MR-1 cultures (42). Based on transcriptome analysis, the genomic response of MR-1 to ionizing radiation (40 Gy) was found to be very similar to its response to UV radiation (41). Similarly, we observed the strong induction of numerous prophage-related genes in MR-1 cells exposed to chromate for 24 h (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental; Table Table2),2), suggesting that prolonged Cr(VI) exposure and/or the accumulation of intracellular chromium may induce the lytic cycle of lysogenic bacteriophage in MR-1. Overall, 16 (21%), 2 (5%), and 10 (19%) ORFs annotated as LambdaSo, MuSo1, and MuSo2 genes, respectively, were significantly induced (more than twofold; P < 0.05) in response to prolonged Cr(VI) exposure. This molecular response was in striking contrast to the differentially expressed genes/proteins characterizing the cellular response to a 90-min acute 1 mM chromate challenge, during which a very small subset of predicted prophage genes (i.e., six) displayed a moderate two- to fourfold induction (14).
Gene products for 14 ORFs with annotations corresponding to mobile and extrachromosomal element functions were confidently identified as being upregulated in response to prolonged Cr(VI) treatment based on both microarray analysis and multidimensional HPLC-MS/MS (Table (Table2).2). Ten of these 14 proteins were encoded in the LambdaSo genome, whereas 1 and 2 gene products were annotated as prophage MuSo1 and MuSo2 proteins, respectively. The majority of these genes encoded such prophage structural proteins as minor and major tail proteins, tail assembly components, and the major head subunit (Table (Table2).2). Six additional ORFs (SO2941, SO2969, SO2973, SO2978, SO2985, and SO3006) with prophage LambdaSo-related functions displayed significant increases (more than twofold; P < 0.05) in mRNA expression (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). The corresponding proteins for those genes were detected only under Cr(VI) conditions by HPLC-MS/MS analysis but failed to meet our filtering criteria for determining differential expression (see Materials and Methods and supplemental Table S2): SO2941 (~20% sequence coverage), SO2969 (protein not detected), SO2973 (~8.2%), SO2978 (~20%), SO2985 (~15%), and SO3006 (~15%). The genes encoded a putative LambdaSo-associated lysozyme (SO2973, 1,053.6-fold), tail assembly protein I (SO2941; 366.5-fold), a putative holin (SO2969, 49.2-fold), a site-specific recombinase (SO2978; 23.8-fold), replication protein O (SO2985, 4.4-fold), and a type II DNA modification methyltransferase (SO3006, 3.8-fold). Other upregulated prophage-related genes had predicted functions in virion morphogenesis (SO2690, 65.6-fold), DNA transposition (SO0644, 11.9-fold; SO2655, 3.6-fold) and circulation (SO2698, 4.1-fold), positive regulation of late transcription (SO2668, 16.7-fold), baseplate (SO2700, 4.3-fold) and tail assembly (SO2699, 10.3-fold; SO2704, 17.2-fold), and assembly of the major head subunit (SO0675, 5.9-fold; SO2685, 24.5-fold), as well as assembly of other structural components (SO2681, 6.2-fold; SO2684, 6.8-fold) (see supplemental Table S1).
Of the differentially expressed proteins determined as having increased abundance under Cr(VI) conditions, 36% corresponded to hypothetical or conserved hypothetical proteins (Table (Table2).2). Five of these proteins (SO2660, SO2663, SO2667, SO2673, and SO2688) are encoded by genes from the MuSo2 genome, which implied their potential function in prophage activation and synthesis. The majority of the upregulated hypothetical proteins (i.e., SO2942, SO2944 to -46, SO2950, SO2951, SO2955, SO2979, SO2980, SO2982, SO2988, SO3001, and SO3008) were derived from genes located in the LambdaSo genome, while no potentially MuSo1-related hypothetical or conserved hypothetical proteins were measured as being differentially expressed under our experimental conditions (Table (Table2).2). Transcriptome analysis revealed a considerably greater number of upregulated hypothetical and conserved hypothetical genes clustered within the three different prophage regions (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). Of the 258 induced hypothetical and conserved hypothetical genes observed, 86 (33%) were distributed among the prophage genomes of MuSo1 (12 genes), MuSo2 (30 genes), and LambdaSo (44 genes).
The gene and protein expression data strongly suggest that, similar to UV irradiation (42), prolonged exposure to chromate and its derivatives may activate the lytic cycles of some or all three of the MR-1 prophages, leading to prophage-mediated cell lysis. At this point, it is not clear whether the induction of MR-1 prophage-related genes/structural proteins is a response to an extended Cr(VI) exposure per se or to the possible intracellular accumulation of chromium, particularly reduced Cr(III). S. oneidensis MR-1 cells exposed to Cr(VI) have been shown to precipitate reduced chromium both extracellularly on the cell surface and as electron-dense globules inside cells (30). Qiu et al. (42) suggested that prophage activation was the major lethal factor in S. oneidensis MR-1 following UV C or UV B irradiation. Our results certainly point to prophage activation as a major contributor to the toxic effects of Cr(VI) under conditions of prolonged exposure and reduction.
Bacterial cell growth occurs by the coordinated periodic alternation of morphogenesis (elongation) and cell division events (19). Central to bacterial morphogenetic processes is biosynthesis and turnover of the peptidoglycan (murein) sacculus, which functions as the principal stress-bearing and shape maintenance component of the cell wall (reviewed in reference 60). The Mur ligases and MurG (a glycosyltransferase) are key enzymes in peptidoglycan biosynthesis and are essential for bacterial cell growth, while formation of the E. coli cell septum involves proteins encoded by fts (temperature-sensitive filamentation) genes (for reviews, see references 15 and 63). Hierarchical assembly of the divisome is initiated by localization of the FtsZ protein near the site of cell division and formation of a ring-like polymeric structure (8), followed by the sequential recruitment of the remaining cell division proteins to the septum site. Cells defective in the expression of FtsI, FtsL, and FtsQ give rise to abnormal morphologies characterized by long, aseptate filaments that eventually lyse (22).
An observation to emerge from the transcriptome study was the upregulation of seven ORFs within the cell wall biosynthesis and cell division gene cluster corresponding to the SO4214-SO4228 region on the MR-1 chromosome. This chromosomal region comprises 15 ORFs with the following gene order: SO4228, SO4227, ftsL (SO4226), ftsI (SO4225), murE (SO4224), murF (SO4223), mraY (SO4222), murD (SO4221), ftsW (SO4220), murG (SO4219), murC (SO4218), ftsQ (SO4217), ftsA (SO4216), ftsZ (SO4215), andlpxC (SO4214) (23) (www.tigr.org). These 15 S. oneidensis MR-1 genes are transcribed from the same DNA strand with very small, if any, intergenic regions, suggesting that the cluster may constitute a single operon. Seven of the 15 genes exhibited statistically significant (P < 0.05) but moderate increases in mRNA expression at the 24-h time point of Cr(VI) exposure: ftsL (2.9-fold), murE (3.4-fold), murF (2.0-fold), mraY (5.4-fold), murG (3.4-fold), ftsA (2.1-fold), and ftsZ (2.1-fold). The corresponding proteins for these genes were not identified as being upregulated by our shotgun proteomics methodology (see supplemental Table S2 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental).
In addition to certain fts and mur genes, S. oneidensis ORFs coding for homologs of E. coli proteins implicated or known to function in determining the rod-like morphology of the cell (rodA and mreD), cell cycle regulation (era), and intracellular septation (ispZ) were induced two- to threefold at the transcript level under our conditions of Cr(VI) exposure. Disruption of the rodA locus in Streptococcus thermophilus resulted in abnormal cell morphologies as well as an increased susceptibility of the mutant to oxidizing agents (e.g., hydrogen peroxide), indicating a possible role for RodA in oxidative stress defense (51). In addition, depleting the cellular concentration of Era or impairing its GTPase activity results in bacterial cell cycle arrest and causes elongation of E. coli cells (12, 21). Taken together, the abnormally elongated morphologies of S. oneidensis cells exposed to prolonged Cr(VI) treatment and upregulated transcript levels for a number of genes putatively involved in cell wall biosynthesis, septation, and morphogenesis might indicate interference with normal cell division by the accumulation of intracellular chromium species. However, more research is needed to substantiate this possible explanation.
The magnitude of the stress imposed by heavy metals and an organism's resistance ability will necessarily intersect metal reduction processes. A recent transcriptome study of S. oneidensis MR-1 exposed to a diverse array of metal (Cr was not included) and nonmetal terminal electron acceptors pointed to the induction of detoxification and stress response genes as playing an important role in the adaptation of MR-1 under anoxic metal-reducing conditions (6). Other transcriptional-profiling investigations of S. oneidensis and Caulobacter crescentus also demonstrated the induction of oxidative stress protection genes in response to chromate treatment (7, 14, 24).
The sensitivity of S. oneidensis MR-1 correlated with the induction of stress-associated genes in response to 24-h Cr(VI) exposure. A number of molecular chaperones or environmental stress response genes were upregulated under our experimental conditions: groES (5.8-fold), groEL (5.8-fold), htpG (4.4-fold), ibpA (4.8-fold), and hslU (4.7-fold) (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). The proteins encoded by these genes mediate the correct folding of a fraction of de novo polypeptides and function in the reactivation or degradation of proteins damaged by different types of stressors. In addition, the gene encoding iron-dependent superoxide dismutase, sodB, was induced 2.3-fold under our experimental conditions. Other transcriptional changes included upregulation of genes coding for ATP-dependent proteases (SO3391, hslV) and a number of ribosomal proteins (SO0226, SO0227, SO0231 to SO0240, SO0243, SO0244, and SO0246 to SO0249).
The membrane response was characterized by changes in the expression of genes encoding outer membrane structural components and polysaccharide biosynthesis proteins. Induced genes included those encoding putative outer membrane porins (encoded by SO0312 and SO1821), OmpW (encoded by SO1673), and three proteins with functions related to polysaccharide biosynthesis (encoded by SO3158, SO3181, and SO3185). Proteomic analysis indicated increases in the synthesis of SO3185 and SO3189, both annotated as polysaccharide biosynthesis proteins, while four predicted lipoproteins (SO2570, NlpD, SOA0110, and SOA0112) and an OmpA family protein (SO3969) belonging to the functional category of cell envelope proteins showed decreased abundance under Cr(VI) conditions (Tables (Tables22 and and3).3). Located immediately upstream of nlpD (SO3433) in the MR-1 chromosome is the gene pcm, which is predicted to encode protein-l-isoaspartate O-methyltransferase, an enzyme involved in protein modification and repair. The ORF coding for the RNA polymerase sigma factor RpoS, which controls the expression of many stationary-phase-induced genes, is positioned just downstream of nlpD and is transcribed in the same direction as nlpD and pcm. Proteomic analysis revealed that, in addition to the lipoprotein NlpD, expression of protein-l-isoaspartate O-methyltransferase was downregulated under conditions of 24-h Cr(VI) exposure (Table (Table3).3). This is of interest because, with age or under stress conditions, proteins are susceptible to various spontaneous or deleterious covalent modifications such as deamidation, the conversion of asparagines into aspartyls and isoaspartyls, which can result in loss of protein function. Pcm functions in repairing damaged proteins by selectively methylating atypical l-isoaspartyl sites and converting them back to l-aspartyls (5, 33). The enzyme was shown to enhance the survival of stationary-phase E. coli subjected to a secondary environmental stress (58). The physiological significance of decreased synthesis of S. oneidensis Pcm under the Cr(VI) conditions used here is unclear.
Chromate and its derivatives have been shown to cause Cr-DNA adducts and DNA-DNA cross-links (11, 27, 28, 43, 47, 53, 59, 64, 65), which can disrupt regulation of gene expression or change DNA topology, leading to altered cell function and perhaps even cell death. Cr(III) can cause DNA damage and negatively affect DNA topology by directly inhibiting topoisomerase DNA relaxation activity (40). Five of the 56 upregulated proteins detected by HPLC-MS/MS under conditions of 24-h Cr(VI) exposure had predicted functions related to DNA metabolism: site-specific recombinase (SO3013), DNA topoisomerase III (TopB, SO3061), type I restriction-modification system subunit (HsdM-2, SO4265), type II restriction endonuclease (SOA0003), and type II DNA modification methyltransferase (SOA0004) (Table (Table2).2). At the transcription level, we observed the low-level induction of three genes annotated as helicase genes in response to Cr(VI) conditions: SO0368 (encoding helicase; 1.7-fold), uvrD (encoding DNA helicase II; 2.1-fold), and hrpA (encoding ATP-dependent helicase; 2.1-fold) (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). DNA helicases catalyze the unwinding of energetically stable duplex DNA and, as a result, play important roles in nearly all aspects of DNA metabolism, including replication, recombination, repair, and transcription (reviewed in reference 54). Additionally, recO, which encodes DNA repair protein RecO, was upregulated twofold under conditions of prolonged Cr(VI) challenge. DNA helicases in Pseudomonas aeruginosa have been implicated in chromate resistance via their involvement in repairing DNA damage (31) and perhaps play a similar role in Cr(VI)-stressed S. oneidensis.
Many facultative anaerobic bacteria synthesize specific primary dehydrogenases and terminal reductases or oxidases as a mechanism for appropriately adapting their respiratory electron-transport chains to environmental changes (45). Sequence analysis of the MR-1 genome revealed extensive duplication of genes involved in electron transport (23). These genes no doubt are important in the ability of S. oneidensis to function as a “respiratory generalist” (23) by allowing MR-1 to adapt to specific growth conditions.
Under our chromate treatment conditions, 9 of the 13 nuo genes (SO1009 to SO1021) encoding subunits of the heteromultimeric NADH dehydrogenase I were upregulated, as identified by microarray analysis: nuoCD, 2.5-fold; nuoE, 2.8-fold; nuoF, 4.4-fold; nuoG, 3.1-fold; nuoI, 4.1-fold; nuoJ, 6.0-fold; nuoL, 3.7-fold; nuoM, 6.1-fold; nuoN, 3.5-fold (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). Transcripts for the other four nuo genes (i.e., nuoA, nuoB, nuoH, and nuoK) either showed decreased abundance under Cr(VI) conditions or failed to meet our criterion of at least a twofold change in expression. None of the NADH dehydrogenase subunits were found to be differentially expressed at the protein level. Additionally, we observed a two- to fivefold induction of genes encoding ubiquinol-cytochrome c reductase (petC, SO0610), tetraheme cytochrome c (SO3056), and cytochrome d ubiquinol oxidase (cydA, SO3286) under conditions of 24-h Cr(VI) exposure (see supplemental Table S1). The S. oneidensis oxidase-encoding gene cydA also was shown to be induced during Cr(VI) reduction in another study (7). Preliminary growth and reduction studies of an MR-1 strain harboring an in-frame deletion of petC indicated that the mutant showed severe growth deficiency compared to wild-type MR-1 under aerobic conditions in the presence of 0.3 mM chromate (S. Barua and K. Chourey, unpublished data).
Five proteins with annotated functions in energy metabolism were identified as being downregulated under chronic Cr(VI) conditions: NADH:ubiquinone oxidoreductase, alpha subunit (NqrA-1); alanine dehydrogenase (Ald); electron transfer flavoprotein, beta subunit (EtfB); coniferyl aldehyde dehydrogenase (SO3683); and cytochrome c552 nitrite reductase (SO3980) (Table (Table3).3). Two of these downregulated proteins (NqrA-1 and EtfB) are involved in electron transport processes. The transcriptomic and proteomic data suggest that S. oneidensis modulates expression of electron transport chain components in response to environmental signals (i.e., growth conditions), resulting in the coexistence of a certain subset of dehydrogenases and reductases/oxidases under environmental Cr(VI) conditions.
Thioredoxins, which have a dithiol/disulfide active site (CGPC), comprise part of a ubiquitous disulfide-reducing system that is involved in a number of important cellular functions, including maintenance of the intracellular redox state and redox regulation of protein function, defense against oxidative stress, and prevention of disulfide bond formation, and function as efficient oxidoreductases of disulfides (for a review, see reference 4). Glutaredoxins catalyze the reduction of disulfides via reduced glutathione in a reaction coupled to glutathione reductase and NADPH. Chromate toxicity is associated with the intracellular reduction of Cr(VI) to lower oxidation states, a process that generates free radicals and imposes oxidative damage on DNA (16). In response to prolonged Cr(VI) exposure, two putative thioredoxin-coding transcripts (SO0476 and SO3117) were upregulated two- to threefold, whereas a gene (SO2100) encoding a thioredoxin family protein was downregulated under the same growth conditions (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). S. oneidensis glutaredoxin (SO2745) also was transcriptionally induced by 6.9-fold under Cr(VI) stress. Other investigators have demonstrated that, while thioredoxin protein products are partially redundant in terms of biochemical function, their gene expression can be differentially regulated depending on the type of stressor (32).
The vast majority of MR-1 proteins downregulated under our experimental conditions are annotated as hypotheticals (Table (Table1).1). Besides poorly characterized proteins, other downregulated proteins belonged to the functional categories of cellular process proteins, transport and binding proteins, cell envelope proteins, and energy metabolism proteins. The relative abundance levels of seven proteins were found to be decreased under 24-h Cr(VI) exposure compared to the control conditions: two chemotaxis proteins (SO1144 and SO3207), a cation transport ATPase (SO1689), an AcrB/AcrD/AcrF family protein (SO3103), and three proteins involved in motility (FliD, FlgE, and SO3314) (Table (Table3).3). Downregulation of proteins involved in motility and chemotaxis was consistent with confocal laser scanning microscopy observations, which indicated a prevalence of nonmotile cells under prolonged Cr(VI) exposure (data not shown). Additional chemotaxis genes displayed transcriptional repression under Cr stress and includedcheY1 (0.4-fold), cheA (0.5-fold), cheW (0.3-fold), and cheB1 (0.5-fold) (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental).
Our previous work showed that a major feature characterizing the S. oneidensis response to acute chromate stress was the high induction of genes encoding components of a TonB1 iron transport system (tonB1-exbB1-exbD1), hemin ATP-binding cassette (ABC) transporters (hmuTUV), and various outer membrane TonB-dependent receptors (14). Almost one-half of the proteins upregulated under chromate shock were predicted to have transport and binding functions. Prolonged exposure to Cr(VI) elicited a very different molecular response. In general, iron sequestration and other transport/binding genes were either repressed in response to 24-h Cr(VI) exposure or induced at substantially lower levels. Genes encoding a heme transport protein (hugA; 0.50-fold) and a TonB1 protein (tonB1; 0.30-fold), for example, were downregulated under the Cr(VI) conditions used in this study, in contrast to the 120- and 99-fold inductions, respectively, observed in response to a 90-min acute exposure (14) (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). Both AlcA (SO3030, a siderophore biosynthesis protein) and IrgA (SO4523, an iron-regulated outer membrane virulence protein) were upregulated under acute chromate stress (14) but were downregulated in response to a 24-h exposure (Table (Table3).3). Additionally, both bfr-1 and bfr-2, encoding bacterioferritin subunits 1 and 2, respectively, were downregulated twofold under our Cr(VI) treatment conditions but were shown to be upregulated during active Cr(VI) reduction by Bencheikh-Latmani et al. (7).
Genes involved in phosphate transport also were downregulated. This was observed only at the transcript level with the repression of two distinct gene clusters with annotated functions in phosphate assimilation: (i) the first cluster (SO1723 to -26) is composed of three genes encoding phosphate ABC transporters (SO1723, 0.30-fold; SO1724, 0.30-fold; pstB1, 0.40-fold) and phoU (0.50-fold), a regulatory gene for phosphate transport that is involved in downregulating the phosphate (Pho) regulon under conditions of excess phosphate; (ii) a second gene cluster (SO1557 to -60) encodes a putative outer membrane porin (SO1557, 0.30-fold), PhoB (SO1558, 0.30-fold), PhoR (SO1559, 0.40-fold), and a phosphate-binding protein (SO1560, 0.20-fold) (see supplemental Table S1 at http://compbio.ornl.gov/shewanella_metal_stress/chronic/supplemental). The outer membrane porin (SO1557) may play a role in phosphate uptake. The signal-transducing proteins of the PhoB-PhoR two-component regulatory system are required for transcriptional activation of the E. coli Pho regulon in response to external inorganic phosphate limitation (reviewed in reference 62). Some of these genes, namely, SO1557, phoB, and SO1723, were observed to be repressed as well in response to anaerobic Cr(VI)-reducing conditions (7). Other transporter genes displaying repression included SO0986 (0.50-fold) for chromate transport and several efflux family genes (SO1918, SO2045, and SOA0159) for multidrug and cation transport. Similarly, gene SO0986, encoding a putative chromate transporter, was not found to be upregulated under anaerobic Cr(VI)-reducing conditions (7).
Combined transcriptomic and proteomic analyses provided a global view of the molecular response induced upon 24-h Cr(VI) challenge of S. oneidensis MR-1. Analytical assays indicated the complete disappearance of Cr(VI) from the growth medium at the 24-h time point of exposure, when cells were harvested for analysis, suggesting the intracellular accumulation of chromate and/or its derivatives. As revealed by integrated transcriptome and proteome analyses, the predominant molecular response of MR-1 to 24-h Cr(VI) challenge was the strong upregulation of numerous prophage-related (LambdaSo, MuSo1, and MuSo2) genes, suggesting that prophage activation may be the primary factor leading to cell lysis under these environmental conditions. These findings were in marked contrast to our previous global examination of the MR-1 temporal response to an acute treatment of 1 mM chromate over a 90-min period (14). Other differences in the gene and protein expression profiles included the downregulation or low-level induction (two- to fourfold) of genes encoding iron binding and transport functions (e.g., putative operons alcA-SO3031-SO3032, tonB1-exbB1-exbD1, hmuTUV), a DNA-binding response regulator (SO2426), and a putative azoreductase (SO3585) under 24-h Cr(VI) exposure (see supplemental Tables S1 and S2 at http://compbio.ornl.gov-/shewanella_metal_stress/chronic/supplemental), whereas these genes and their corresponding proteins were upregulated at high expression levels (in some cases, >50-fold) in response to an acute chromate challenge (14). Transcriptional induction of cell wall biosynthesis and cell division genes (i.e., the putative mur and fts operon), as well as rodA, mreD, and era, appeared to correlate with the extreme filamentous morphology observed using confocal microscopy and AFM topographic imaging. The induction of genes encoding molecular chaperones, classical heat shock proteins, oxidative stress protection proteins, thioredoxins, and glutaredoxins correlated with the susceptibility of S. oneidensis MR-1 to prolonged Cr(VI) exposure.
We thank Liyou Wu for whole-genome S. oneidensis MR-1 microarrays. We are grateful to Tingfen Yan and Xueduan Liu for technical assistance in the construction of whole-genome MR-1 microarrays. M. R. Thompson acknowledges support from the University of Tennessee (Knoxville)-ORNL Graduate School of Genome Science and Technology.
This research was supported by the U.S. Department of Energy, Office of Biological and Environmental Research, Environmental Remediation Sciences Program. Oak Ridge National Laboratory is managed by University of Tennessee-Battelle LLC for the Department of Energy under contract DOE-AC05-00OR22725.