Molecular Dynamics of the Shewanella oneidensis Response to Chromate Stress* □ S

Temporal genomic profiling and whole-cell proteomic analyses were performed to characterize the dynamic molecular response of the metal-reducing bacterium Shewanella oneidensis MR-1 to an acute chromate shock. The complex dynamics of cellular processes demand the integration of methodologies that describe biological systems at the levels of regulation, gene and protein expression, and metabolite production. Genomic microarray analysis of the transcriptome dynamics of midexponential phase cells subjected to 1 m M potassium chromate (K 2 CrO 4 ) at exposure time intervals of 5, 30, 60, and 90 min revealed 910 genes that were differentially expressed at one or more time points. Strongly induced genes included those encoding components of a TonB1 iron transport system ( tonB1 - exbB1 - exbD1 ), hemin ATP-binding cas-sette transporters ( hmuTUV ), TonB-dependent receptors as well as sulfate transporters ( cysP , cysW-2 , and cysA-2 ), and enzymes involved in assimilative sulfur metabolism ( cysC , cysN , cysD , cysH , cysI , and cysJ ). Transcript levels for genes with annotated functions in DNA repair ( lexA , recX , recA , recN , dinP , and umuD ), cellular detoxification ( so1756 , so3585 , and so3586 ), and two-component signal transduction systems ( so2426 ) were also significantly up-regulated ( p < 0.05) in Cr(VI)-exposed cells relative to untreated cells. By contrast, genes with functions linked to energy metabolism, particularly electron transport ( e.g. so0902 - 03 - 04 , mtrA , omcA ,

Temporal genomic profiling and whole-cell proteomic analyses were performed to characterize the dynamic molecular response of the metal-reducing bacterium Shewanella oneidensis MR-1 to an acute chromate shock. The complex dynamics of cellular processes demand the integration of methodologies that describe biological systems at the levels of regulation, gene and protein expression, and metabolite production. Genomic microarray analysis of the transcriptome dynamics of midexponential phase cells subjected to 1 mM potassium chromate (K 2 CrO 4 ) at exposure time intervals of 5, 30, 60, and 90 min revealed 910 genes that were differentially expressed at one or more time points. Strongly induced genes included those encoding components of a TonB1 iron transport system (tonB1-exbB1-exbD1), hemin ATP-binding cassette transporters (hmuTUV), TonB-dependent receptors as well as sulfate transporters (cysP, cysW-2, and cysA-2), and enzymes involved in assimilative sulfur metabolism (cysC, cysN, cysD, cysH, cysI, and cysJ). Transcript levels for genes with annotated functions in DNA repair (lexA, recX, recA, recN, dinP, and umuD), cellular detoxification (so1756, so3585, and so3586), and two-component signal transduction systems (so2426) were also significantly upregulated (p < 0.05) in Cr(VI)-exposed cells relative to untreated cells. By contrast, genes with functions linked to energy metabolism, particularly electron transport (e.g. so0902-03-04, mtrA, omcA, and omcB), showed dramatic temporal alterations in expression with the majority exhibiting repression. Differential proteomics based on multidimensional HPLC-MS/MS was used to complement the transcriptome data, resulting in comparable induction and repression patterns for a subset of corresponding proteins. In total, expression of 2,370 proteins were confidently verified with 624 (26%) of these annotated as hypothetical or conserved hypothetical proteins. The initial response of S. oneidensis to chromate shock appears to require a combination of different regulatory networks that involve genes with annotated functions in oxidative stress protection, detoxification, protein stress protection, iron and sulfur acquisition, and SOS-controlled DNA repair mechanisms. ) species, is a serious anthropogenic pollutant found in the discharged effluents from various manufacturing and military industries (1,2) and is one of several risk-driving contaminants found at unacceptable levels in groundwater as well as soil and sediment at the United States Department of Energy waste sites (3). The most chemically stable and common forms of chromium in the environment are the trivalent Cr(III) and hexavalent Cr(VI) species, which differ significantly in their physicochemical properties and biological responsiveness (i.e. toxicity). Cr(VI) is highly soluble and a toxic non-essential metal for most organisms with chronic chromate exposure leading to mutagenesis and carcinogenesis, whereas Cr(III) is sparingly soluble and relatively innocuous (4 -6).
In situ microbial bioreduction of Cr(VI) to Cr(III) is of particular interest because it may serve as a potential strategy for the detoxification and immobilization of chromate compared with more cost-prohibitive physical and chemical treatment methods (7). Microorganisms exhibiting Cr(VI)-reducing activities and resistance have been detected in chromate-contaminated sites as well as natural, uncontaminated ecosystems (8 -10). Previous studies have indicated that species of the ␥-proteobacterial genus Shewanella are capable of both direct (enzymatic) dissimilatory Cr(VI) reduction (11)(12)(13) and indirect (chemical) Cr(VI) reduction driven by the reduction of Fe(III) to Fe(II) (14 -16). Myers et al. (12) demonstrated that formate-dependent Cr(VI) reductase activity was localized to the cytoplasmic membrane of anaerobically grown Shewanella oneidensis MR-1 (previously named Shewanella putrefaciens MR-1). Similarly chromate reductase activity was shown to be preferentially associated with the membrane fraction of Enterobacter cloacae cells (17), whereas soluble Cr(VI)-reducing enzymes have been purified to varying degrees from Pseudomonas putida MK1 (18), P. putida PRS2000 (19), and Pseudomonas ambigua G-1 (20). Reductases with other primary cellular functions have been shown to be efficient in chromate reduction (21), suggesting that a variety of enzymes may participate in the transfer of electrons to Cr(VI). Studies investigating the effect of nitrite on hexavalent chromium reduction indicated that S. oneidensis MR-1 might possess multiple nonspecific Cr(VI) reduction mechanisms, as well as metal resistance mechanisms, that are dependent on physiological growth conditions (13).
Oxidatively induced DNA damage is considered to be the basis of chromate genotoxicity (22)(23)(24) with active transport of chromate across biological membranes being mediated by the sulfate transport system in prokaryotic and eukaryotic cells (7,25,26). The chromate resistance mechanisms displayed by microorganisms are diverse and 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 less toxic Cr(III) (for a review, see Ref. 7). Plasmid-determined resistance to chromate has been shown to occur in bacteria, including species of the genera Pseudomonas (27)(28)(29) and Alcaligenes (25). A hydrophobic protein, designated ChrA, was found to be responsible for the plasmid-specified resistance phenotype in these organisms (30,31) and appears to function as a secondary transport system for the extrusion of chromate ions (32). Although Cr(VI) reduction and toxicity resistance mechanisms are considered to be unlinked cellular processes, the biotransformation of Cr(VI) to Cr(III) likely contributes to the detoxification of chromate (33,34).
The goal of the study described here was to obtain global insight into temporal alterations in mRNA expression and protein synthesis that occur in response to toxic acute levels of chromate and thus to gain an understanding of the potential molecular mechanisms enabling Cr(VI) detoxification under aerobic respiratory conditions. For the most part, the observed alterations in the protein patterns reflected the changes identified at the transcriptomic level. In addition, differential proteomics revealed post-transcriptional levels of regulation that cannot be captured by microarray analysis. The combined transcriptome and proteome analyses suggested a close relationship between chromate stress and cellular iron requirement (or limitation) as indicated by the dramatic induction of genes involved in iron sequestration and uptake. Sulfate transport and assimilatory pathways for cysteine production, DNA damage repair systems, and oxidative stress protection were also major features of the initial response of S. oneidensis MR-1 to acute Cr(VI) exposure. The present study identified a number of genes and their encoded products that have been shown previously to be important for the oxidative stress response of MR-1 (35) and other bacteria. Moreover a number of other genes/proteins not described previously were implicated in the cellular detoxification of chromate. This study represents an essential first step toward a global molecular characterization of the cellular response to acute chromate exposure in a metal-reducing bacterium and provides information necessary for future examination of metal stress-linked gene regulatory networks.

EXPERIMENTAL PROCEDURES
Bacterial Strains, Culture Conditions, and Assays-S. oneidensis strain MR-1 (36) was used for both transcriptomic and proteomic analyses and was propagated in Luria-Bertani (LB) 1 medium (pH 7.2) at 30°C under aerobic conditions. A Bioscreen C microbiological culture system (Growth Curves USA, Piscataway, NJ) was used to monitor the growth of MR-1 in response to different concentrations of K 2 CrO 4 over 48 h. Chromate reduction was measured using the 1,5-diphenylcarbazide method essentially as described previously (18).
RNA Isolation and Preparation of Fluorescein-labeled cDNA-For the time series microarray experiment, a 1:1000 dilution of a fresh overnight culture (16 h) of S. oneidensis MR-1 was used to inoculate prewarmed LB medium in six 250-ml sidearm Pyrex flasks. Batch cultures were grown to midexponential phase (A 600 , 0.5) followed by the addition of prewarmed 2 M K 2 CrO 4 to a final concentration of 1 mM for three of the six cultures. The remaining three cultures served as the reference (control) samples and were grown in parallel with the treated cells. For microarray hybridization, the control and treated cells were harvested in parallel for total cellular RNA extraction at 5, 30, 60, and 90 min post-K 2 CrO 4 addition. RNA isolation, fluorescent labeling reactions, probe purification, and microarray hybridization were performed as described previously (37).
Microarray Hybridization, Scanning, Image Quantification, and Data Analysis-S. oneidensis MR-1 microarray construction, hybridization, scanning, image quantification, and data analyses were performed as described previously (38,39). Briefly temporal gene expression analysis was performed using six microarrays for each time point (three biological replicates ϫ two dye reversal reactions) with each slide containing two spots representing each gene at different array locations for a total of 12 signal intensity measurements per gene per time point. The two separately labeled cDNA pools (i.e. the K 2 CrO 4 -treated and the corresponding control time point) to be compared were mixed together in a hybridization solution containing 50% (v/v) formamide. Microarray hybridization signals were quantified using Im-aGene Version 5.5 (Biodiscovery, Inc., Los Angeles, CA) followed by data transformation and normalization using GeneSite Light (Biodiscovery, Inc.). ArrayStat TM (Imaging Research, Inc., Ontario, Canada) was used to determine the common error of these values, remove outliers, and determine statistical significance via a z test for two independent conditions and the false discovery rate method (nominal ␣, p Ͻ 0.05). Genes exhibiting significant changes in expression at a degree of 2-fold or greater (40) were further analyzed using the program Hierarchical Clustering Explorer Version 3.0 (www.cs.umd. edu/hcil/multi-cluster/).
Preparation of Proteomes for HPLC-MS/MS Analysis-All chemical reagents were obtained from Sigma unless stated otherwise. Modified sequencing grade trypsin (Promega, Madison, WI) was used in all digestions. HPLC grade water and acetonitrile were acquired from Burdick & Jackson (Muskegon, MI), and 99% formic acid was purchased from EM Science (Darmstadt, Germany).
For large scale proteomic characterization, 500-ml cultures of S. oneidensis MR-1 in 4-liter flasks (a total of 1 liter of culture for treatment and control) were grown to midexponential phase (A 600 , 0.5) under the same conditions as described above for the microarray studies and then either exposed to a final K 2 CrO 4 concentration of 1 mM or allowed to continue growing in the absence of added chromate. At time points of 45 and 90 min after treatment with K 2 CrO 4 , cells were harvested from each of the following conditions for HPLC-MS/MS analysis: 1) Control 1 (untreated midlog phase cells after 45 min of further growth), 2) Treatment 1 (45 min post-Cr addition), 3) Control 2 (untreated midlog phase cells after 90 min of further growth), and 4) Treatment 2 (90 min post-Cr addition). For this, cells were pelleted by centrifugation (5,000 ϫ g for 5 min), resuspended in ice-cold LB medium, washed two times in 50 mM Tris, 10 mM EDTA (pH 7.6), and centrifuged at 5,000 ϫ g for 10 min. The S. oneidensis cells were then placed on ice and lysed by sonication using a microprobe at high power with 30-s pulses for five times with a 30-s cooling period between each sonication. Cellular debris were removed by centrifugation at 5,000 ϫ g for 10 min. The supernatant was centrifuged at 100,000 ϫ g for 60 min in an ultracentrifuge to separate a soluble fraction from a pellet. The pellet (membrane fraction) was washed with 50 mM Tris, 10 mM EDTA (pH 7.6) and centrifuged at 100,000 ϫ g for 60 min; this fraction was then resuspended in 50 mM Tris, 10 mM EDTA (pH 7.6) by brief sonication. Both proteome fractions were quantified using bicinchoninic acid (BCA) analysis, aliquoted, and stored at Ϫ80°C until ready for digestion. Approximately 2 mg of each proteome fraction (soluble and membrane) prepared for the two growth conditions was denatured and reduced in 6 M guanidine and 10 mM DTT (60°C for 1 h). The denatured/reduced proteome mixture was diluted 6-fold with 50 mM Tris, 10 mM CaCl 2 (pH 7.8), and sequencing grade trypsin was added at 1:100 (protease/protein (w/ w)). The digestions were run with gentle shaking at 37°C for 18 h followed by a second addition of trypsin at 1:100 and an additional 5 h incubation. The samples were treated with 20 mM DTT for 1 h at 37°C as a final reduction step to remove remaining disulfide bonds and then immediately desalted using Sep-Pak Plus C 18 solid phase extraction (Waters, Milford, MA). All samples were concentrated and solvent-exchanged into 0.1% formic acid in water by centrifugal evaporation to ϳ10 g/l starting material, filtered, aliquoted, and stored at Ϫ80°C until ready for LC-MS/MS analysis.
LC/LC-MS/MS Analysis-The proteome fractions (soluble and membrane) prepared from control and chromate-treated samples were analyzed in duplicate via two-dimensional (2-D) LC-MS/MS experiments using an Ultimate HPLC system (LC Packings, a division of Dionex, San Francisco, CA) coupled to a linear trapping quadrupole (LTQ) mass spectrometer (ThermoFinnigan, San Jose, CA). The HPLC pump provided a flow rate of ϳ100 l/min that was split precolumn to achieve a final flow rate of ϳ200 nl/min at the nanospray tip. A split phase column (150-m inner diameter fused silica) was packed via a pressure cell as follows: ϳ3.5 cm of strong cation exchange (Luna SCX, 5 m, 100 Å; Phenomenex, Torrance, CA) followed by ϳ3.5 cm of C 18 reverse phase (Aqua C 18 , 5 m, 200 Å; Phenomenex). Subse-quently ϳ500 g of sample was loaded onto the split phase column via a pressure cell. The loaded split phase column was then inserted behind a PicoFrit tip (100-m inner diameter, 15-m inner diameter at the tip; New Objective, Woburn, MA) packed via a pressure cell with ϳ15 cm of C 18 reverse phase (Jupiter C 18 , 5 m, 300 Å; Phenomenex). This entire column system was positioned in front of the LTQ on a nanospray source (ThermoFinnigan).
All samples were analyzed via a 24-h 12-step 2-D analysis consisting of increasing concentration (0 -500 mM) salt pulses of ammonium acetate followed by 2-h reverse phase gradients from 100% aqueous solvent (95% H 2 O, 5% ACN, 0.1% formic acid) to 50% organic solvent (30% H 2 O, 70% ACN, 0.1% formic acid). During the entire chromatographic process, the LTQ was operated in a data-dependent MS/MS mode detailed below. The chromatographic methods and HPLC columns were identical for all analyses. The LC-MS/MS system was fully automated and under direct control of the Xcalibur software system (ThermoFinnigan). The LTQ was operated with a nanospray voltage of 2.6 kV, heated capillary temperature of 200°C, and a full scan m/z range of 400 -1700. Data-dependent MS/MS mode was operated as follows. Five MS/MS spectra were acquired following every full scan, two microscans were averaged for every full MS and MS/MS spectra, a 3 m/z isolation width was used, and 35% collision energy was used for fragmentation. The dynamic exclusion was set to 1 with an exclusion duration of 3 min.
Proteome Bioinformatics-A protein database was created by combining the most recent version of the S. oneidensis MR-1 database (Version 8; www.tigr.org/) containing a total of 4,798 predicted proteins with 36 common contaminants (trypsin, keratin, etc.). The database can be downloaded from the website compbio.ornl.gov/ shewanella_chromium_stress/databases/. For all database searches, the MS/MS spectra were searched using SEQUEST (Ref. 41; Ther-moFinnigan) with the following parameters: enzyme type, trypsin; parent mass tolerance, 3.0; fragment ion tolerance, 0.5; up to four missed cleavages allowed; fully tryptic peptides only. The MS/MS spectra from individual RAW files were first converted to mzXML format by using ReAdW software written at the Institute for Systems Biology in Seattle, WA (www.systemsbiology.org) and can be downloaded from the SourceForge repository (sashimi.sourceforge.net). Individual spectra were then converted to DTA files by mzXML2Other, also from the Institute for Systems Biology. DTA files are the required format for input into SEQUEST (see Ref. 52). The output data files were then filtered and sorted with the DTASelect algorithm (42) using the following parameters: fully tryptic peptides only with ⌬CN of at least 0.08 and cross-correlation scores (Xcorrs) of at least 1.8 (ϩ1), 2.5 (ϩ2), and 3.5 (ϩ3). These threshold scores have been tested rigorously in our laboratory and provide a high confidence of identification (see Refs. 43 and 52 for more discussion) with a maximum false-positive rate of 1-2%. Post-translational modifications and other fixed modifications (i.e. due to addition of iodoacetamide) were not included in the search parameters. DTASelect files are available on the analysis page (compbio.ornl.gov/shewanella_chromium_ stress/ms_analysis) under the corresponding dataset and are filtered at one peptide and two peptides per protein. The files are presented in a text format or a viewable html version where every identified spectrum can be viewed by clicking on the spectral number (first column, labeled by filename). The DTASelect results from all control and chromate-treated samples were then compared with the Contrast program (42) for each time point. These results are located under the global contrast heading on the analysis page. A list was made of all proteins showing a reproducible significant change of at least 40% sequence coverage, five or more unique peptides, and/or a reproducible spectral count difference of 2ϫ between the control and chromate-treated samples at each time point (adapted from Refs. 39 and 43). The analysis page also contains inter-run contrast files (compares duplicate runs on same sample) as well as fractionation comparisons (compares replicate runs of the same proteome broken down by fraction).

RESULTS
Physiological Effect of Chromate during Aerobic Growth of S. oneidensis-Initially we investigated the effects of chromate-induced stress on the growth of aerobically growing cultures of S. oneidensis MR-1 to identify adequate concentrations for transcriptome and proteome analyses. Three independent experiments based on end point culture turbidity determinations indicated that MR-1 responded to K 2 CrO 4 in a dose-dependent manner. MR-1 cells that were grown in the presence of 0.06 and 0.125 mM Cr(VI) attained A 600 values similar to those attained for cells grown in the absence of heavy metal in the medium after 48 h (Fig. 1A). MR-1 cultured with 0.25, 0.5, and 1 mM K 2 CrO 4 showed ϳ17, 50, and 66% reductions in growth, respectively. Culture turbidity did not increase for inocula in medium containing 2 mM K 2 CrO 4 after 48 h.
To provide a complete description of the growth behavior exhibited by MR-1 in response to different amounts of Cr(VI), culture turbidity was measured at 30-min intervals under aerobic conditions over the course of 48 h using a Bioscreen C reader. As shown in Fig. 1B, growth rates of MR-1 in the presence of 0.06 and 0.125 mM K 2 CrO 4 were decreased slightly relative to growth in the absence of added K 2 CrO 4 , but similar maximal turbidity measurements were obtained after 48 h. Extended lag periods, ranging from 20 -40 h, followed by biphasic growth with a period of growth arrest and then by complete or partial recovery were observed for MR-1 grown in LB medium containing higher Cr(VI) concentrations (i.e. 0.25, 0.5, and 1 mM; Fig. 1B). However, MR-1 was unable to grow after 48 h in 2 mM K 2 CrO 4 under the conditions tested in the present study (Fig. 1B).
Chromate Reduction-A chromate concentration of 1 mM was selected for global gene and protein expression profiling. To determine whether chromate reduction occurs in the presence of exponentially growing S. oneidensis cells upon addition of 1 mM K 2 CrO 4 , a spectrophotometric method using 1,5-diphenylcarbazide (18) was used to monitor chromate conversion over the entire time course of the transcriptome and proteome analyses. No chromate disappearance was observed for Cr(VI)-shocked midlog phase cells over a 150min period (see Supplemental Fig. S1). Similarly no abiotic conversion of chromate was detected in the LB broth-only control (Supplemental Fig. S1), thus eliminating the possible contribution of Fe 3ϩ in the chemical reduction of chromate. After 24 h of incubation in LB medium containing 1 mM K 2 CrO 4 , wild-type MR-1 cells transformed ϳ42% of Cr(VI) (results not shown).
General Transcriptome Dynamic Patterns in Response to Chromate Shock and Array Data Validation-To define the repertoire and temporal expression patterns of MR-1 genes responding to acute Cr(VI) exposure (1 mM K 2 CrO 4 ), transcriptome dynamics were examined based on time series DNA microarray experiments at 5, 30, 60, and 90 min postshock using whole-genome microarrays. Temporal gene expression profiles of cells exposed to Cr(VI) were compared with those of the untreated control cells grown in parallel. Approximately 20% (n ϭ 910) of the total predicted S. oneidensis genes (n ϭ 4,648) represented on the microarray showed at least a 2-fold statistically significant (p Ͻ 0.05) change in expression for at least one of the time points during acute Cr(VI) stress exposure. Generally the number of up-and down-regulated genes increased with time of exposure to chromate. At all four time points, as many as 54 genes were induced more than 2-fold, whereas only 12 genes were repressed at all time points   2). The early response (5 min) was characterized predominantly by up-regulation of genes whose products mostly belonged to the functional categories of energy metabolism, hypothetical proteins, and transport/binding proteins ( Fig.  2A). General expression patterns during the middle time intervals (30 and 60 min) of K 2 CrO 4 exposure were characterized by a marked increase in the number of down-regulated genes in the following functional categories: biosynthesis of cofactors, prosthetic groups, and carriers; cellular processes; energy metabolism; hypothetical proteins; protein fate; purines, pyrimidines, nucleosides, and nucleotides; regulatory functions; transport and binding proteins; and proteins of unknown function (Fig. 2, B and C).  interval, we observed a large increase in the number of repressed genes involved in energy metabolism. Real time quantitative RT-PCR was used to provide an independent assessment of gene expression for selected genes across the entire time course of chromate treatment. The six genes (tonB1, tonB2, exbB1, so2426, so3032, and so3585) selected for comparative real time QRT-PCR analysis displayed up-regulated or unchanged expression patterns as identified by microarray analysis. A comparison of gene expression measurements determined by QRT-PCR and microarray analyses indicated that the datasets were highly correlated with an r 2 value of 0.81 (p Ͻ 0.0001) (see Supplemental Fig. S3).
Induction of Iron Binding and Transport Genes-Of the 910 genes shown to be significantly differentially regulated, we identified a small subgroup comprising 16 genes that exhibited substantially greater expression levels (in some cases Ͼ100-fold) compared with the other induced Cr(VI)-responsive genes (Table II). Subgroup IIC (see Supplemental Table  S1 and Table II) was dominated by genes encoding products with functions in iron binding and transport: a putative siderophore biosynthesis protein (so3032), a ferric alcaligin siderophore receptor (so3033), a probable heme transport operon (hugA-hugX-hugZ) (44), three TonB1 complex proteins (tonB1, exbB1, and exbD1), and a hemin ABC transporter (hmuTUV) ( Table II). The expression levels of these affected genes remained induced throughout the time course of acute Cr(VI) exposure. Other genes in the siderophore biosynthesis operon (so3030-31-32) and those encoding TonB-dependent receptors (so1482 and so1580), a ferric vibriobactin receptor (so4516), and an iron-regulated outer membrane virulence protein (so4523) were also up-regulated during the 90-min treatment period but at expression levels much lower in mag-nitude compared with the genes comprising the putative tonB1-exbB1-exbD1 operon (Tables II and III).
A number of genes located in different subgroups, whose products are also predicted to be involved in sequestering iron, were induced but at lower levels (e.g. 2-5-fold) compared with those in Subgroup IIC (see Supplemental Table S1 and Table II). These genes included so3034 (encoding a putative ferric iron reductase protein), ftn (encoding ferritin), and so1158 (encoding a ferritin-like Dps protein). Although the majority of iron binding and transport genes were induced, several genes such as so1112 and so1111 (encoding bacterioferritin subunits 1 and 2, respectively) were repressed about 2-5-fold over the time course, and others like bcp, which is predicted to encode a bacterioferritin comigratory protein, did not change in response to chromate.
Induction of Detoxification, DNA Damage Repair, and Other Stress-related Genes after Chromate Exposure-The biological basis of chromate-induced toxicity is thought to be the generation of reactive oxygen species (ROS) during enzymatic chromate reduction (for reviews, see Refs. 5 and 7). Cr(V) species are transient intermediates in the flavoenzyme-catalyzed one-electron reduction of chromate whose redox cycling results in the formation of ROS and H 2 O 2 generation, which has been shown to be involved in partial chromate reduction (20,46,47). Recently the central mechanism of chromium toxicity in Saccharomyces cerevisiae was shown to involve oxidative damage to cellular proteins (48). As a strategy for scavenging ROS, bacteria commonly modulate gene expression by inducing genes encoding antioxidant enzymes and proteins. Cr(VI)-stressed S. oneidensis MR-1 indicated a 2-5-fold induction of katG-1 (SO0725, catalase/peroxidase hydroperoxidase), katB (SO1070, catalase), and so4640 (an- a Relative gene expression (-fold induction) is presented as the mean ratio of the fluorescence intensity of K 2 CrO 4 -exposed cells to control cells. Each gene showed significant differential expression (p Ͻ 0.05) at least at one time point.
b Time in minutes at which cells were harvested for RNA isolation following addition of 1 mM K 2 CrO 4 to the experimental culture. c Genes so3667 (hugZ) and so3668 (hugX) were recently reannotated and predicted to encode heme iron utilization proteins (45).
tioxidant AhpC/Tsa family protein) at the 60-and 90-min time points (see Supplemental Table S1). Differential expression was not measured for sodB, predicted to encode an ironcofactored superoxide dismutase. The SOS regulon is a genetic network that enables Escherichia coli and related bacteria to maximize their chances of survival when exposed to environmental stresses that dam-age DNA (for a review, see Ref. 49). Key components in homologous recombination and the SOS pathway of DNA repair were up-regulated at the transcription level in response to Cr(VI) exposure. These up-regulated MR-1 repair genes encoded proteins similar to the LexA (SO4603) repressor, regulatory protein RecX (SO3429), RecA (SO3430), and DNA repair protein RecN (SO3462), which in E. coli comprise the a Relative gene expression (-fold induction) is presented as the mean ratio of the fluorescence intensity of K 2 CrO 4 -exposed cells to control cells. Each gene showed significant differential expression (p Ͻ 0.05) at least at one time point.
b Time in minutes at which cells were harvested for RNA isolation following addition of 1 mM K 2 CrO 4 to the experimental culture.
core of the SOS response. Expression of lexA, recX, recA, and recN increased substantially after 5 min of acute Cr(VI) exposure and then remained at similar induced levels (ϳ4 -6-fold) over the rest of the time course (Table III). Other known SOS-controlled genes, like dinP (so1114) and plasmid-borne umuD (soa0013), also displayed expression profiles like that of recA (Table III), thus suggesting the induction of an SOSlike pathway of DNA repair and mutagenesis in S. oneidensis MR-1 in response to chromate. The lexA and so4604 genes are coupled by overlapping stop and start codons and coinduced in response to UV irradiation, suggesting expression as an operon (35). The downstream so4605 gene was in the lexA and so4604 expression cluster in response to both UV and chromate stresses, although the putative initiation codon does not overlap with the stop codon of the so4604 gene. A putative azoreductase gene (so3585) with an annotated function in detoxification, a glyoxalase family gene (so3586) of undefined cellular function, and a gene encoding a hypothetical protein (so3587) (44) were located within the cluster of highly induced iron-sequestering genes (Table II). Maximal transcriptional induction of so3585 (61-fold), so3586 (26-fold), and so3587 (18-fold) occurred at the 30-min time point followed by continued up-regulation at the 60-and 90-min time points but at lower levels (Table II). These genes have not been shown previously to be involved in the cellular response to chromate stress. so3585, so3586, and so3587 are transcribed in the same direction and exhibit co-regulated temporal expression patterns in response to acute Cr(VI) exposure, suggesting that the encoded proteins might function together in a complex. Another gene predicted to encode a glyoxalase family protein (SO1756) exhibited a different temporal expression pattern, characterized by a maximum induction of 12-fold at the 90-min time point (Table III).
Up-regulated Expression of Sulfate Transport and Sulfur Metabolism Genes-The structural similarity of chromate (CrO 4 2Ϫ ) anions to such biologically important inorganic anions as SO 4 2Ϫ and PO 4 3Ϫ most likely constitutes the basis for its active transport across cell membranes via the sulfate transport system (7). Chromate has been shown to be a competitive inhibitor of sulfate transport in Pseudomonas fluorescens (26). Based on our array experiments, another subset of MR-1 ORFs that were up-regulated in response to chromate consisted of those genes encoding sulfate ABC transporters (cysP, sbp, cysW-2, and cysA-2) and enzymes involved in activation and subsequent reduction of sulfate to sulfide (cysC, cysDN, cysH, and cysIJ) (Table III; for a review of sulfur metabolism, see Ref. 50). These genes and predicted operons, along with five hypothetical genes (so3724, so3725, so4650, so4651, and so4656), displayed a distinctive bimodal pattern of temporal expression with an initial induction peak (ϳ4 -10-fold) at the 30-min time interval followed by a marked decrease in transcription and then a second induction (ϳ13-89-fold) at 90 min. The putative serine acetyltransferase gene (cysE) was induced ϳ2-fold at 90 min; however, cysM and cysK, encoding the two enzymes that catalyze the reaction of O-acetylserine with sulfide to yield cysteine, failed to meet our criteria for differential expression. The so2263-69 genes are immediately downstream of cysE and, with the exception of so2268, were induced between 2-and 6-fold at one or more time points. In E. coli, the B2531-2525 proteins are involved in a sulfur transfer cascade required for the biosynthesis of thiamine, NAD, Fe-S clusters, and thionucleosides (51), and the SO2263-69 proteins showed 65-92% sequence identity between the respective protein counterparts.
Genes Implicated in Regulating the Cellular Response to Chromate-Chromate stress had a greater impact on the expression of genes linked to transcriptional regulatory functions than those with signal transduction functions (Fig. 2). The 90-min time point yielded the most differential expression for transcriptional regulators with 13 and 17 genes being induced and repressed, respectively. These genes included a number of functionally undefined transcriptional regulators from the MerR, TetR, GntR, MarR, and DeoR families (Supplemental Table S1). Genes encoding signal transduction functions were not largely affected by acute chromate treatment: at any one time point only one to three signal transduction genes were induced, and three or four genes were repressed at a single point during the entire exposure period. Most noteworthy was the up-regulation of gene so2426, encoding a DNA-binding response regulator of a two-component signaling system. This gene showed a temporal fold induction profile of 3.7 (5 min), 3.3 (30 min), 3.5 (60 min), and 10.7 (90 min) after challenge with 1 mM chromate (Table III). Interestingly this gene was also markedly induced in response to other heavy metal conditions, such as iron overload (39) and nonradioactive strontium stress (37), and may be involved, at least in part, in the coordinate regulation of the cellular response to heavy metal toxicity.
Chromate-mediated Gene Repression-Next to genes encoding hypothetical proteins, energy metabolism showed the most profound temporal changes in terms of the number of genes down-regulated in response to acute chromate exposure (Fig. 2). Although the total number of induced genes (11-19 genes) remained relatively constant over the time course, the number of repressed genes with functions in energy metabolism changed dramatically with 5, 34, 44, and 67 ORFs showing significant decreases in mRNA abundance at the 5-, 30-, 60-, and 90-min time points, respectively. Many of these more highly repressed energy metabolism genes had annotated subrole functions in amino acids and amines (e.g. so1962 and so3774), anaerobic energy metabolism (e.g. fdnG, fdnH, fdnI, fdhE, fdhD, and fdrC), electron transport (e.g. nqrABC-1, mtrA, omcA, omcB, and hydAB), and fermentation (adhE) ( Table IV). As shown in Table IV, transcription of genes functioning in cellular processes (e.g. ccpA), central intermediary metabolism (speF), protein synthesis (selA), transcriptional regulation (e.g. so1422, so3297, and so3627), signal transduction (so4155), and transport and binding (e.g. trkH-2 and modAB) was significantly repressed under chromate stress conditions. Proteome Characterization of the Chromate Stress Response-The transcriptome data were compared with and complemented by proteome studies based on mass spectrometry techniques. Cellular fractions from each of the four S. oneidensis growth conditions (i.e. Control 1, Control 2, Treatment 1, and Treatment 2) were analyzed in duplicate using 2-D LC-ES-MS/MS with a linear trapping quadrupole mass spectrometer (see "Experimental Procedures" for details). All datasets were searched using SEQUEST (41), filtered and sorted with DTASelect and Contrast (42), and can be accessed directly at compbio.ornl.gov/shewanella_chromium_ stress/. The DTASelect files are presented in a viewable html version where every identified spectrum can be viewed by clicking on the spectral number (first column, labeled by filename). The threshold scores used for identification of expressed proteins have been tested rigorously in our laboratory to give a high confidence (see Refs. 43 and 52) and a maximum false-positive rate of 1-2%.
A total of 2,370 of the 4,931 total predicted genes in the S.
oneidensis MR-1 genome were identified with at least two peptides (Table V), representing 48% of the theoretical proteome. Due to the large number of false positives possible at the one-peptide filter level (52), we present a rigorous analysis of the two-peptide dataset only. High stringency filtering was used in this study, giving a maximum false-positive rate of 1-2%. The reproducibility between duplicate protein analyses on the LTQ was as follows: 78.6% (chromate-shocked) and 78.4% (control) for the 45-min poststress time point and 77.7% (chromate-shocked) and 73.5% (control) for the 90min time point. This level of reproducibility is necessary for semiquantification. Variation is likely due to low abundance proteins identified with two peptides in one of the analyses and only one peptide in another, thereby being filtered out. Although previous studies using 2-D PAGE and LC-MS/MS  a Relative gene expression (-fold repression) is presented as the mean ratio of the fluorescence intensity of K 2 CrO 4 -exposed cells to control cells. Each gene showed significant differential expression (p Ͻ 0.05) at least one time point.
b Time in minutes at which cells were harvested for RNA isolation following addition of 1 mM K 2 CrO 4 to the experimental culture. c Gcn5-related N-acetyltransferase.
have been used for global S. oneidensis proteome studies (39,53,54), this study represents, to our knowledge, the largest measurement of the S. oneidensis proteome published to date. The entire list of identified proteins with total sequence coverage, functional categories, pI, and molecular weight information is given in Supplemental Table S2. No major biases were found between the pI and molecular weights of the predicted proteome from the genome and the observed proteome (Supplemental Table S2). The identified proteins with their peptide count (number of identified peptides), spectral count (number of MS/MS spectra identified per protein), and percent sequence coverage (total percentage of the protein sequence covered by tryptic peptides) for the different growth conditions and individual analyses can be found in Supplemental Table S3.
Proteins identified at the two-peptide level were grouped according to the functional categories assigned by The Institute for Genomic Research annotation (Ref. 44; www.tigr.org, Comprehensive Microbial Resource) (Table VI). Proteins with assigned functions in amino acid biosynthesis, cellular processes, protein fate, protein synthesis, nucleotide metabolism, and transcription were found with greater than 70% identified. More than 90% of the proteins comprising the functional classes of protein synthesis, nucleotide metabolism, and transcription were identified, representing an almost complete characterization of these categories at the proteome level. Proteins generally thought to be of lower abundance, such as those with assigned functions in signal transduction and transcriptional regulation, were identified at levels of 64 and 52%, respectively, of the total number of proteins per category. A total of 624 of the 2,039 predicted hypothetical proteins were identified. Of the 624 proteins in this functional category, 209 were annotated as hypothetical proteins, and 415 were annotated as conserved hypothetical proteins. This represents one of the largest identifications of hypothetical protein expression for a microbial proteome to date.
Comparison of Transcriptome and Proteome Data-Although absolute quantification at the global proteome level was not feasible, semiquantification of differentially expressed proteins between the chromate-treated and control samples could be accomplished by using a combination of percent sequence coverage, number of unique peptides, and spectral count from the mass spectra (39, 43, 55) (see "Experimental Procedures" for details). Other methods for relative quantification, such as ICAT, were considered before the method used here was selected. ICAT labels the cysteine residues of a protein, and in S. oneidensis greater than 50% of the predicted proteins have two or fewer cysteine residues present. Moreover almost 20% of the predicted proteome does not contain any cysteine residues. This method would not detect almost 70% of the predicted proteome confidently; thus ICAT was not deemed appropriate for this study.
Comparisons were made between the transcriptome and proteome data to determine the relationship between gene and protein expression. Both up-and down-regulated proteins were measured by comparing the chromate-shocked cell samples to their respective control samples. Using the semiquantitative criteria discussed above, we identified 78 proteins (Supplemental Table S4) as being differentially expressed in response to chromate (Tables VII and VIII and  Supplemental Table S4). Supplemental Table S4 divides the proteins identified at 45 min postshock from the proteins identified at 90 min postshock. We propose that reproducibility of at least 70% between replicate analyses at the protein level is necessary for a successful determination of differentially expressed proteins where a protein must be detected in both of the replicate analyses to be included as a candidate for semiquantification. For the 45-min chromate treatment analysis, a total of 24 proteins were found to be up-regulated with six additional proteins being down-regulated relative to the control sample (Supplemental Table S4). The genes for 23 of these 24 upregulated proteins showed corresponding induction levels at the 30-and 60-min time points based on microarray hybrid-  ization. A putative formate acetyltransferase (encoded by pflB), which was identified as being up-regulated at the protein level, was found to be induced at the transcript level at the 30-min treatment time but repressed at the 60-min time point. The subset of proteins identified as up-regulated were dominated by species (12 total) assigned to the functional category of transport and binding proteins and included TonB-dependent receptors, siderophore biosynthesis proteins, heme transport proteins, and TonB1. Proteins classified as hypothetical were also dominant with a total of six conserved hypothetical proteins identified as up-regulated. For the six proteins down-regulated in the 45-min shocked sample relative to the control condition, five were found to be repressed at the mRNA level by microarray analysis, and the other protein, a hypothetical protein (SOA0141), revealed no change at the transcript level. Proteins that were measured as being down-regulated in the 45-min chromate-shocked samples consisted of three hypothetical proteins (SO1124, SO2929, and SOA0141), a periplasmic nitrate reductase (NapA), a transcriptional regulator (HlyU), and the ␣ subunit of formate dehydrogenase (SO4513).
Proteomic analysis of the 90-min chromate treatment sam-ples revealed 48 up-regulated proteins and 26 down-regulated proteins relative to the control sample (Supplemental Table S4 and Tables VII and VIII   mate addition were also repressed at the transcript level. Of the remaining down-regulated proteins, five showed no change at the mRNA level, one protein of which was found to exhibit no mRNA change at the 90-min time point. Two downregulated proteins (NapG and NapA) were found to be induced at the transcript level. Four of the six proteins downregulated at the 45-min time point were also repressed at 90 min. The transcriptome data revealed an increase in the number of down-regulated energy metabolism genes over time in response to chromate. This trend was also reflected in the proteomic data in which 13 proteins with functions in energy metabolism were identified as repressed at the 90-min time point relative to the 45-min interval.
Differences between the transcriptomic and proteomic responses of S. oneidensis to chromate shock are due to the stringency of filtering used in the proteome study and inherent measurement differences between microarray versus proteome technology. Supplemental Table S5 demonstrates the correlation of the identified and/or differentially expressed proteins to the top 100 mRNAs induced at each time point. The list is composed of a total of 194 mRNAs, which is a concatenated list representing the top 100 mRNAs at each time point (initial list was 400 genes in total before removal of redundant genes). A total of 67% of the genes in Supplemental Table S5 did not meet the differential criteria at the protein level where 45% were not identified by mass spectrometry and another 22% were not considered differentially expressed. The 67% of genes found not to meet differential expression criteria corresponds to 130 genes of 194, which is comparative in percentage to the total global analysis. DISCUSSION A global analysis of the S. oneidensis temporal response to a toxic concentration (1 mM) of chromate was performed by integrating whole-genome microarray hybridization and 2-D HPLC-MS/MS techniques. This study provides the first comprehensive molecular description of the response of a dissimilatory metal-reducing bacterium to an acute challenge of toxic chromate, a serious anthropogenic pollutant found in diverse environmental settings. Five major features characterized the molecular dynamics underlying the S. oneidensis response to acute chromate stress: 1) the induction of genes/ proteins involved in iron transport and sulfur assimilation; 2) co-regulated induction of the gene cluster so3585, so3586, and so3587; 3) up-regulation of DNA repair genes and their encoded protein products; 4) differential expression of a large number of hypothetical genes/proteins; and 5) a dramatic increase in the number of down-regulated energy metabolism genes and proteins over time.
As revealed by transcriptome and proteome analyses, a major feature of the molecular response of S. oneidensis MR-1 to acute chromate challenge was the differential regulation of the TonB1-ExbB1-ExbD1 complex, an integral inner membrane system for iron transport, as well as other genes involved in exogenous iron acquisition. Genes that were highly induced based on time series microarray experiments (Table II) and were identified as up-regulated proteins at the 45-and/or 90-min chromate treatment conditions (Table VII) included two putative siderophore biosynthesis proteins (AlcA and SO3032), a ferric alcaligin siderophore receptor (SO3033), HugA heme transport protein (SO3669), TonB1 (SO3670), a putative TonB-dependent receptor (SO3914), a TonB-dependent heme receptor (SO1580), ViuA (SO4516), and HmuT (SO3673) and HmuV (SO3675) of the hemin ABC transporter. Wang and Newton (56) demonstrated that mutants harboring a deletion of the tonB-trp region of the E. coli chromosome were sensitive to chromic ion (Cr 3ϩ ) due to defective iron transport systems, and residual iron uptake by these strains was shown to be inhibited by chromic ion. Regulation of iron homeostasis is primarily carried out by the Fur protein (for a review, see Ref. 57). It has been suggested previously that iron uptake regulation may not be the only function of Fur but that it may also serve to sequester iron to prevent the generation of highly reactive hydroxyl radicals via Fenton reactions (58). The putative MR-1 ferritin genes, but not bacterioferritin genes, were induced in response to chromate, and these respective iron storage proteins have been suggested to have roles in short term iron flux and long term iron storage in E. coli (59).
In addition to iron transport genes, our global analyses demonstrated that CysP, CysC, CysD, CysN, CysI, CysJ, Sbp, and CysA-2 are up-regulated at both the mRNA and protein levels in response to chromate treatment (Tables III  and VII). The enhanced expression of genes encoding proteins involved in sulfate transport and assimilatory sulfate reduction suggests the possibility of chromate-induced sulfur limitation in S. oneidensis, perhaps through competitive inhibition of sulfate uptake by chromate, as has been shown previously (26,30). Partial reduction of Cr(VI) to Cr(V) produces ROS (6,7,20), leading to chromate-mediated oxidative stress. Researchers working with a number of different bacteria have observed induction of genes involved in sulfur and iron homeostasis following different oxidative stresses (35, 60 -62). A variety of explanations have been proposed including disruption of intracellular redox cycling leading to insufficient sulfite reduction, a reduction in cysteine biosynthesis correlated with cell envelope damage and subsequent leakage of sulfide (63), and increased demand for low molecular weight protective thiol-containing compounds such as glutathione (64). Alternatively induction of genes involved in sulfur metabolism and iron sequestration might represent an adaptive response to sulfur and iron limitation in MR-1 following chromate exposure.
One of the potentially interesting findings to emerge from this integrated global investigation was the co-regulated expression of a cluster of three genes (so3585, so3586, and so3587) at both the mRNA and protein levels. All three genes are transcribed in the same direction on the MR-1 chromo-some and show a similar transcriptional profile in response to chromate with the peak in up-regulated expression occurring at the 30-min time point (Table II). The proteins encoded by two of these genes (so3585 and so3586) were detected in the chromate-treated samples only (Table VII), suggesting that expression of SO3585 and SO3586 was differentially regulated in response to chromate stress conditions. By contrast, hypothetical protein SO3587 was found in both the control and chromate-shocked samples at the two-peptide level (Supplemental Table S3) even though the gene encoding this protein was shown to be up-regulated over the entire time course in response to chromate stress. In addition, SO3587 was found only in the membrane fractions, and a hydrophobicity plot analysis using the computer program SOSUI (65) identified a putative transmembrane domain (IGIALIFADVS-LYLAYFFVGLGV) in SO3587. SO3585 and SO3586 were detected in both the soluble and membrane fractions. Based on their proximity in genome location and co-regulated expression profile within the context of chromate stress, we predict that SO3585, SO3586, and SO3587 function together as a protein complex associated with the cell membrane and play an important role in the response of the cell to chromate toxicity.
S. oneidensis SO3585 and SO3586 are annotated as a putative azoreductase and glyoxalase family protein, respectively (44). Glyoxalase systems are known to serve as key detoxification routes for preventing the intracellular accumulation of methylglyoxal, a natural metabolite with toxic electrophilic properties (for a review, see Ref. 66). Azoreductases are responsible for the reductive cleavage of azo dyes, synthetic organic colorants used extensively in the textile, food, and cosmetics industries. Synthetic azo dyes are not readily reduced under aerobic conditions and are considered pollutants. Protein database searches using BLAST Version 2.2.12 (67) with the derived SO3585 primary sequence revealed ϳ28% sequence identity with P. putida ChrR and E. coli YieF, two soluble flavoproteins that have been demonstrated to exhibit chromate reductase activity (18,46,68). Regions of conservation in the derived amino acid sequence of SO3585 included the characteristic signature, LFVTPEYNX 6 LKNAIDX 2 S (conserved residues in SO3585 are underlined), of the NADH_dh2 family of NAD(P)H oxidoreductases (Ref. 46; results not shown). Recently further investigation demonstrated that the P. putida ChrR functions as a quinone reductase and minimizes oxidative stress induced by intracellular H 2 O 2 , which is generated during the course of chromate reduction (69). An MR-1 strain carrying an in-frame deletion mutation in the so3585 locus has been created, and future studies will characterize it to gain insight into the functional role of SO3585 and to assess the importance of azoreductase in the S. oneidensis response to chromate.
Another primary feature of the molecular response to chromate was the stimulation of genes whose protein products are involved in DNA repair mechanisms, thus indicating a contri-bution of genotoxic damage to chromate toxicity in S. oneidensis. Genes with functions linked to DNA damage, recombination, and repair (dinP, recX, recA, recN, lexA, and umuD) exhibited similar temporal expression profiles when cells were challenged with 1 mM chromate: the transcript levels for all six genes increased only slightly (1.5-2-fold) early in the treatment followed by induction plateaus of 4 -6-fold from the 30-min time point onward (Table III). Only two of these genes (dinP and recN) showed corresponding increases in the abundance levels of their expressed proteins at both the 45-and 90-min time points following chromate treatment (Table VII); RecA and LexA were also identified at the protein level but did not show differential expression. Members of the E. coli SOS regulon, which includes genes encoding the DNA damageinducible (Din) proteins, are controlled by the RecA coprotease activity and LexA repressor and are transcriptionally induced following DNA damage (for a review, see Ref. 49). The induction of homologs for such known SOS-controlled genes as recA, lexA, recN, and umuD suggests chromateinduced DNA damage in S. oneidensis and subsequent RecA activation and LexA cleavage under the selected conditions. Qiu et al. (35) similarly observed induction of the MR-1 SOS response and active detoxification mechanisms following oxidative damage caused by UV radiation. The precise mechanism underlying chromate-induced DNA damage in S. oneidensis has not been determined. However, oxidative DNA damage, triggered by Cr(V) complexes reacting with H 2 O 2 to generate hydroxyl radicals, is considered the basis of chromium genotoxicity (7,(22)(23)(24) and likely produces an SOS-like response in MR-1. The study by Qiu et al. (35) and this investigation revealed two genes (so4604 and so4605) downstream of the lexA gene that were induced following different oxidative stress conditions, and their expression profiles clustered together with lexA. However, further studies are required to determine whether these genes are directly involved in the MR-1 SOS response.
In addition to members comprising a potential SOS regulon in S. oneidensis, we demonstrated the induction of topB (encodes a DNA topoisomerase III) at both the mRNA and protein levels. Topoisomerase III was initially described as exhibiting superhelical DNA relaxing activity (70). A recent study by Nurse et al. (71) demonstrated that topoisomerase III can function as the principal cellular decatenase, capable of uncoupling replicating daughter chromosomes in vivo. Microarray analysis revealed a temporal expression pattern for S. oneidensis topB resembling that of recA and other DNA damage-inducible genes (Supplemental Table S1). A chromate-dependent increase in TopB abundance levels was observed only at the 90-min time point (Table VII). Trivalent chromium has been reported to cause DNA damage and inhibit type II topoisomerase-mediated DNA relaxation activity in bacterial cells (72). Although we did not detect chromate reduction during the treatment period (Supplemental Fig. S1), our transcriptome and proteome data suggest that chromate has an impact on DNA topology leading to adjustments in the expression levels for cellular topoisomerases.
The large portion of genes with unassigned cellular functions among Cr(VI)-induced or Cr(VI)-repressed genes reveals how much of the molecular response of S. oneidensis to chromate toxicity remains to be explored. At each time point after chromate addition, ϳ39 -50% of the total number of induced genes and 31-47% of the total number of repressed genes were functionally classified as hypothetical or conserved hypothetical proteins. Fifty percent of the genes upregulated at the 90-min interval in response to chromate were functionally unknown, whereas poorly characterized genes constituted 31% of the total number of down-regulated genes at that time point. Nine proteins (SO0798, SO0934, SO1045, SO1190, SO3667, SO3907, SO3913, SO4651, and SOA0042) annotated as hypothetical or conserved hypothetical met our criteria of significance for differential expression and were identified as being up-regulated in response to chromate exposure at the 45-and/or 90-min time points (Table VII). For four of these proteins (SO0934, SO1045, SO3907, and SOA0042), we observed no significant change in expression at the mRNA level, suggesting post-transcriptional regulation of these proteins in response to chromate. Our integrated transcriptome and proteome study implicates these differentially regulated proteins of unknown function in the initial response of MR-1 to toxic chromate, thus revealing gene candidates for future functional analysis.
A recent study by Kolker et al. (45) analyzed expression for a subset of 538 hypothetical proteins that were confidently identified in S. oneidensis MR-1 as a result of large scale microarray and proteomic analyses of cell samples generated under different growth conditions. A total of 788 hypothetical proteins have been identified based on the study by Kolker et al. (45) and the present study: 368 of these functionally undefined proteins were found in both studies, 170 were found only by Kolker et al. (45), and 256 were found only in this study (Supplemental Table S6). The 368 hypothetical proteins identified independently by both studies should be considered as expressed proteins, and their annotations should be changed to unknown or conserved unknown (43). Most of these proteins were found under all growth conditions and identified in most of the replicates. Differences between the datasets revealed by this proteomic study and the one reported by Kolker et al. (45) are likely due to differences in the growth conditions used.
Finally the S. oneidensis response to chromate was distinguished by a time-dependent and dramatic increase in the number of down-regulated genes with functions in energy metabolism (Fig. 2). With the exception of the hypothetical protein category, chromate had the greatest impact on energy metabolism in terms of the number of repressed genes: five (5 min), 34 (30 min), 44 (60 min), and 67 (90 min). Genes with annotated functions in anaerobic energy metabolism and electron transport were particularly affected by chromate treatment with ORFs encoding the formate dehydrogenase subunits (fdnGHI), NADH:ubiquinone oxidoreductase subunits (nqrABC-1), MtrA, OmcA, OmcB, and periplasmic iron hydrogenase subunits (hydAB) showing some of the highest repression -fold values at the 60-or 90-min time point (Table  IV). Within this subset of genes, proteins NqrA-1, MtrA, OmcA, OmcB, and HydA were also detected at decreased abundance levels (Table VIII). In addition, MS-based proteomics revealed potential instances of post-transcriptional regulation of gene expression in response to chromate stress. For example, proteome analysis indicated complete repression of the anaerobic dimethyl sulfoxide reductase subunits A (SO1429) and B (SO1430) under chromate treatment conditions, whereas no significant change in expression at the transcript level was observed for dmaA-1 and dmsB-1 (Table VIII).
In summary, we used integrated transcriptome and proteome analyses to reveal a global view of the response of the metal-reducing bacterium S. oneidensis MR-1 to the challenge of acute chromate stress. The molecular response of S. oneidensis to chromate shock elicited a distinctively different transcriptional profile compared with Cr(VI) reduction by MR-1 in a study published recently (73). Approximately 20% of the genes up-regulated during chromate reduction encode putative cytochromes, cytochrome synthesis proteins, or noncytochrome reductases (which were repressed or not differentially expressed in the present study), consistent with cells being potentially limited for iron under Cr(VI)-shocked conditions. The chromate shock response of S. oneidensis requires a combination of different regulatory networks that involve genes with annotated functions in oxidative stress protection, detoxification, protein stress protection, iron and sulfur acquisition, and DNA repair mechanisms.