<title>Background</title> <p><italic>Shewanella oneidensis </italic>MR-1, a facultatively anaerobic member of the γ-class of <italic>Proteobacteria</italic>, can respire anaerobically by reducing a wide variety of alternative electron acceptors including fumarate, nitrate, elemental sulfur, and such oxidized metals as Fe(III), Mn(IV), and U(VI) [##REF##15518832##1##, ####UREF##0##2##, ##REF##11536738##3##, ##REF##17815852##4##, ##REF##2172208##5####2172208##5##]. In addition, MR-1 is capable of reducing hexavalent chromium, or chromate [Cr(VI)], to less toxic and sparingly soluble trivalent chromium [Cr(III)] under both aerobic and anaerobic conditions [##REF##12889027##6##, ####REF##10735248##7##, ##REF##12001169##8####12001169##8##]; however, there have been no reports to date of the ability of MR-1 to generate energy for growth using Cr(VI) as the sole terminal electron acceptor for anaerobic respiration [##REF##16269787##9##]. This plasticity in the use of alternative electron acceptors for anaerobic respiration has engendered interest in <italic>S. oneidensis </italic>MR-1 as a model environmental organism with potential utility in the bioremediation of dissolved metal ions, and as a consequence, the complete MR-1 genome was sequenced to advance understanding of <italic>Shewanella </italic>biology [##REF##12368813##10##].</p> <p>To expand knowledge of metal stress responses in particular, we previously initiated genome-based studies focused on identifying the molecular components involved in the <italic>S. oneidensis </italic>response to chromate [##REF##16524964##11##, ####REF##16957260##12##, ##REF##17385904##13####17385904##13##], an anthropogenic pollutant widely distributed in the environment due to its prevalent use in manufacturing and military industries [##UREF##1##14##,##UREF##2##15##]. Chromate toxicity is associated with the generation of reactive oxygen species during the intracellular partial reduction of Cr(VI) to the highly reactive radical Cr(V) by various <italic>in vivo </italic>nonspecific reductants or cellular one-electron reductases [##REF##11348688##16##]. Our previous work employing DNA microarrays and multidimensional liquid chromatography-tandem mass spectrometry (LC-MS/MS) demonstrated that a functionally unknown DNA-binding response regulator (designated SO2426) in MR-1 was reproducibly and significantly upregulated at both the transcript and protein levels in response to acute chromate challenge [##REF##16524964##11##,##REF##17385904##13##].</p> <p>In accordance with its ability to adapt to various environmental conditions, the <italic>S. oneidensis </italic>MR-1 genome encodes a relatively large repertoire of transcriptional regulators, including 88 predicted two-component regulatory system proteins consisting of 23 histidine protein kinase (HK) genes, 57 response regulators (RR), and 8 HK-RR hybrids [##REF##12368813##10##]. Prototypical two-component systems, which constitute the predominant mechanism used by bacteria for coupling environmental signals to specific adaptive responses, comprise a sensor histidine kinase and a cognate response regulator [##REF##10745001##17##, ####UREF##3##18##, ##REF##10966457##19####10966457##19##]. Fundamental to signal transduction is a phosphotransfer mechanism that forms the core of both mechanistically simple as well as more complex signaling pathways. In this scheme, a membrane-bound sensor kinase catalyzes ATP-dependent autophosphorylation at a conserved histidine residue in the HK cytoplasmic autokinase domain in response to a specific environmental cue. The signal is transduced through transfer of the phosphoryl group from the phosphoHis of the HK to a conserved aspartate residue within the receiver module of a cytosolic RR. Phosphorylation of the RR alters its affinity for DNA, leading to the transcriptional activation or repression of specific genes. An essential feature of two-component systems is that specific histidine kinases and response regulators function as cognate pairs and are often arranged within the same operon. The <italic>S. oneidensis </italic>SO2426 protein appears to be an orphan RR as defined by the fact that there is no obvious cognate histidine kinase flanking the <italic>so2426 </italic>gene or proximally located based on the genome annotation [##REF##12368813##10##].</p> <p>The primary aim of the present study was to investigate the function of the uncharacterized SO2426 regulator within the context of chromate stress by performing a genome-wide screen for potential target genes of this RR in <italic>S. oneidensis </italic>MR-1. Toward this end, we created an <italic>so2426 </italic>in-frame deletion mutant and compared its global transcriptional profiles with that of wild-type MR-1 in response to chromate challenge using time-resolved microarray experiments. Comparative transcriptome analysis suggested that SO2426 was associated with the activation of genes involved in siderophore-mediated Fe uptake, Fe storage, and possibly other cation transport processes.</p>
<title>Methods</title> <title>Bacterial strains, plasmids, and growth conditions</title> <p>All bacterial strains and plasmids used in this study are listed in Table ##TAB##1##2##. <italic>S. oneidensis </italic>strain MR-1 [##REF##17815852##4##,##REF##10319494##53##] was used as the wild type. The Δ<italic>so2426 </italic>strain, a derivative of MR-1, contains an in-frame deletion of the <italic>so2426 </italic>locus on the chromosome (see below for a detailed description of mutant construction). For growth studies, wild-type and mutant strains were cultivated aerobically with shaking (250 rpm) at 30°C in Luria-Bertani (LB) medium (pH 7.2) alone or in media amended with 0.3 mM (final concentration) of K<sub>2</sub>CrO<sub>4</sub>. Optical density was monitored at a wavelength of 600 nm (OD<sub>600</sub>) 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 [##REF##16391131##36##].</p> <p>For transcriptome analyses, three separate batch cultures each of wild-type MR-1 and the Δ<italic>so2426 </italic>mutant were grown to mid-exponential phase (OD<sub>600</sub>, 0.5) in 100 ml of LB broth, followed by addition of 2 M K<sub>2</sub>CrO<sub>4 </sub>to a final concentration of 0.3 mM. The two different sets of chromate-challenged cultures (wild type and mutant), each consisting of three biological replicates, were grown aerobically in parallel at 30°C. Cells were harvested for RNA extraction at 5, 30, 60, 90, 180 min, and 24 h post-chromate exposure. For RT-PCR and 5' RACE experiments, wild-type <italic>S. oneidensis </italic>MR-1 cells were grown to mid-log phase (OD<sub>600</sub>, 0.5) in LB broth at 30°C and then exposed to chromate (final concentration, 1 mM) for 30 min before harvesting cells for RNA extraction. Untreated MR-1 cells were grown in parallel and used as the control.</p> <title>Chromate reduction and CAS siderophore assays</title> <p>Disappearance of extracellular Cr(VI) was quantified spectrophotometrically using the 1,5-diphenylcarbazide (DPC) method as described elsewhere [##REF##10788340##34##]. Cultures were assayed at different time points (5, 30, 60, 90, 180 min, and 24 h) post chromate-challenge 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). Cell-free LB medium containing 0.3 mM K<sub>2</sub>CrO<sub>4 </sub>served as the abiotic negative control and was monitored in parallel with the experimental samples.</p> <p>The chrome azurol S (CAS) assay for detection of siderophore production was performed essentially as described by Schwyn and Neilands [##REF##2952030##54##]. Siderophore biosynthesis and excretion by the Δ<italic>so2426 </italic>mutant was compared with those by the wild-type MR-1 strain in LB broth with or without the addition of 50 μM FeCl<sub>3 </sub>or 0.3 mM K<sub>2</sub>CrO<sub>4</sub>. Cultures of these strains were grown aerobically to stationary phase (OD<sub>600 </sub>&gt; 1.0) at 30°C for 24 h. Cell-free supernatants of <italic>S. oneidensis </italic>cultures were mixed 1:1 with the CAS assay solution and equilibrated at room temperature for 2 h before the absorbance at 630 nm was measured. The relative siderophore production was calculated as the ratio of the control (uninoculated LB medium) <italic>A</italic><sub>630 </sub>to the wild-type or mutant strain supernatant <italic>A</italic><sub>630</sub>. All CAS measurements were performed in triplicate, and at least two independent determinations were conducted.</p> <title>Construction of an <italic>S. oneidensis so2426 </italic>deletion strain and complementation</title> <p>An <italic>S. oneidensis </italic>MR-1 in-frame deletion mutant lacking the <italic>so2426 </italic>locus was generated using the <italic>cre-lox </italic>recombination system as described elsewhere [##REF##16391095##32##,##REF##12449384##33##]. The application of this mutagenesis strategy to targeted <italic>S. oneidensis </italic>MR-1 genes has been described previously [##REF##16740962##55##]. The primer sequences used in the generation and verification of the Δ<italic>so2426 </italic>mutant are listed in Table ##TAB##2##3##. PCR was used to amplify a 900-bp region upstream and a 840-bp region downstream of the <italic>so2426 </italic>gene using primer pairs Del900-F/Del900-R and Del840-F/Del840-R, respectively. The amplified regions were cloned into the kanamycin-resistant (Km<sup>r</sup>) plasmid pJK100 [##REF##16391095##32##]. The resulting plasmid, pJK102, carrying a correct construct was introduced into the <italic>E. coli </italic>strain WM3064 [##REF##16391095##32##], a diaminopimelic acid (DAP) auxotroph, and then subsequently moved into MR-1 via conjugation with WM3064. The resultant MR-1 strain was Km<sup>r</sup>/Tet<sup>s </sup>and contained the <italic>loxP-Km</italic><sup><italic>r</italic></sup>-<italic>loxP </italic>cassette. Removal of the kanamycin cassette and the helper plasmid was performed as described elsewhere [##REF##12449384##33##]. A cre-recombinase enzyme-producing plasmid, pCM157/Tet<sup>r</sup>, was introduced into MR-1 to resolve the <italic>loxP-Km</italic><sup><italic>r</italic></sup>-<italic>loxP </italic>cassette through conjugation. The Tet<sup>r </sup>plasmid pCM157 was cured from MR-1 by continuous culturing in non-selective LB medium. The in-frame deletion of <italic>so2426 </italic>was confirmed by PCR amplification using several sets of primers (ISU2426/ISD2426, OP840, OP900, kanF/kanR in Table ##TAB##2##3##) and DNA sequencing.</p> <p>The Δ<italic>so2426 </italic>mutant was complemented by reintroduction of the wild-type <italic>so2426 </italic>allele on a low-copy-number plasmid. For this, the <italic>so2426 </italic>open reading frame with 100 bp of upstream sequence was amplified from <italic>S. oneidensis </italic>MR-1 genomic DNA using primers 2426com-F and 2426com-R (Table ##TAB##2##3##). The PCR product was gel purified, digested, and ligated into the <italic>Bam</italic>H1 and <italic>Hind</italic>III sites of the plasmid pBBR1MCS-5 [##REF##8529885##56##], which contains a gentamicin resistance (Gm<sup>r</sup>) cassette. The resultant construct harboring the complete <italic>S. oneidensis so2426 </italic>gene with its endogenous promoter was designated pBB2426. The pBB2426 plasmid was electroporated into <italic>E. coli </italic>strain WM3064 and then introduced into the <italic>S. oneidensis </italic>Δ<italic>so2426 </italic>mutant and wild-type MR-1 strains via conjugation. In addition, plasmid pBBR1MCS5 without the insert (empty vector) was also transferred into the mutant and wild-type strains via conjugation with WM3064. Gm<sup>r </sup><italic>Shewanella </italic>colonies were verified for the presence of the <italic>so2426 </italic>gene by PCR analysis and DNA sequencing.</p> <title>RNA isolation</title> <p>For microarray profiling and real-time RT-PCR experiments, total cellular RNA was isolated from <italic>S. oneidensis </italic>cultures using the TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. RNA preparations were treated with RNase-free DNase I (Ambion Applied Biosystems, Foster City, CA) to remove residual genomic DNA and subsequently further purified using the RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's RNA cleanup protocol. Total RNA was quantitated as described previously [##REF##16391131##36##]. For RT-PCR and 5'-RACE analyses, total RNA was extracted from <italic>S. oneidensis </italic>cultures using the RNeasy Mini Kit (Qiagen). Chromosomal DNA contamination was removed by incubating total purified RNA with 3 U of RNase-free DNase (Ambion Applied Biosystems) for 30 min at 37°C. The DNA digestion reaction was performed twice for each RNA sample. The quantity and purity of the RNA was assessed as described previously [##REF##16391131##36##]. The integrity of all RNA samples was assessed visually using 1% agarose gel electrophoresis and ethidium bromide staining.</p> <title>Microarray hybridizations and data analysis</title> <p>Transcriptome analyses were performed with a microarray containing 4,197 PCR amplicons and 451 specific 50-mer oligonucleotides, covering approximately 94% of the total predicted gene content of <italic>S. oneidensis </italic>MR-1. Fabrication of the MR-1 arrays has been described in detail elsewhere [##REF##15516594##57##]. Synthesis of the two differentially labeled cDNA pools (wild-type MR-1 and the Δ<italic>so2426 </italic>mutant) to be compared, microarray prehybridization and hybridization, and post-hybridization washings were performed as described previously [##REF##16391131##36##]. Temporal gene expression analysis was performed using six independent microarray hybridizations (three biological replicates × two dye-swap reactions) for each of six time points, 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. Image quantification, data normalization, and analysis of gene expression data for statistical significance were conducted as described by Brown <italic>et al</italic>. [##REF##16391131##36##]. The time-series microarray expression profiles of the Δ<italic>so2426 </italic>mutant relative to the parental strain were clustered using Hierarchical Clustering Explorer (HCE) [##REF##16418236##58##]. During the clustering process, only genes with an expression value of at least ≥ 2-fold or ≤ 0.5-fold in one or more of 6 expression profiling time points were included in the analyses. As a result, a dataset of 841 genes was clustered based on average linkage using Euclidean distance.</p> <title>Access to microarray data</title> <p>The microarray data reported in this study have been deposited in MIAME-compliant format at Gene Expression Omnibus on the NCBI website [##UREF##6##59##] under series accession number GSE12129. The statistically analyzed microarray output is provided in Additional file ##SUPPL##0##1##.</p> <title>Quantitative real-time RT-PCR</title> <p>Reverse transcriptase, quantitative real-time PCR (qRT-PCR) was used to provide an independent assessment of gene expression for five selected genes (<italic>so0404</italic>, <italic>so2426</italic>, <italic>so3670</italic>, <italic>so1826</italic>, and <italic>so3585</italic>), which exhibited different expression patterns (<italic>i.e</italic>., down-regulated or no change). These genes are predicted to encode a hypothetical protein (SO0404), a DNA-binding response regulator (SO2426), proteins involved in transport and binding of cations (SO3670 and SO1826), 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 <ext-link ext-link-type="uri" xlink:href="http://www-genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi"/>: SO0404, 5'-AGTATAACCAAGCGCCAGTA and 3'-GCATCGGTATTAACTTGCTC; SO2426, 5'-GCAGAAGGATTTAGGTCGAT and 3'-CGGTGTTGATTAAAGTACGC; SO3670, 5'-TCTAAACAGTCGCAGGAGCA and 3'-GCGCCATATTGCTATCCATT; SO1826, 5'-GGGTGTCCCAAGCTAGT CAA and 3'-GAGCATTACTCGTCCCCTGA; SO3585, 5'-CGAGGCTATCCATC ACTTAG and 3'-TGGAAAACACGATAAAGACC. Reverse transcription was performed on the 180-min time point samples with 3 μg of total cellular RNA and 2.5 mM random hexamers using Superscript™ II RNase H<sup>- </sup>Reverse Transcriptase (Invitrogen) as described previously [##REF##15576789##31##]. Quantitative PCR was carried out in an iCycler iQ<sup>® </sup>real-time PCR system (BioRad, Hercules, CA) in 50-μl reaction mixtures containing 1 μl cDNA, 600 nM forward and reverse primers, and iQ SYBR green supermix (BioRad) according to the manufacturer's instructions and conducted under the following conditions: 30 sec at 95°C, followed by 40 cycles of 15 sec at 95°C, 30 sec at the specific annealing temperature, and 30 sec at 72°C. The qRT-PCR reactions were performed in triplicate for each of the three biological replicates tested. Standards for each gene of interest were included in the analysis as described previously [##REF##15576789##31##]. The final Δ<italic>so2426</italic>/WT ratio for each target gene and the standard error were calculated as described elsewhere [##REF##15576789##31##], and the linear correlation between the qRT-PCR and microarray data was determined based on the log mean values using Sigma-Plot version 9.0 (SPSS Inc., Chicago, IL).</p> <title>RT-PCR and 5'-RACE analyses</title> <p>Reverse transcription-PCR (RT-PCR) was performed to investigate the transcription of the <italic>so2427</italic>-<italic>so2422 </italic>gene cluster on the MR-1 chromosome. First-strand cDNA was synthesized from 1 μg of DNase-treated total RNA isolated from MR-1 cultures grown in LB broth under either non-stress (no chromate added) or chromate stress (1 mM, 30-min exposure) conditions. Random hexamers (250 ng), dNTPs (10 mM), 5× First-Strand buffer (4 μl), 0.1 M DTT (1 μl), RNaseOUT (1 μl), recombinant RNase inhibitor (40 U), and Superscript™ II Reverse Transcriptase (200 U; Invitrogen) were added to the reaction mixture, which was incubated at 25°C for 10 min, 42°C for 50 min, and then 70°C for 15 min in accordance with the recommendations of the supplier. Synthesized cDNA was used as the template in subsequent PCR amplification reactions with gene-specific oligonucleotide primers P1–P9 (Table ##TAB##2##3##) spanning the <italic>so2427</italic>-<italic>so2422 </italic>gene region. The TripleMaster PCR System (Eppendorf, Westbury, NY) was used for amplification, with primers added to each reaction mixture at a final concentration of 250 nmol. The reaction mixtures were subjected to 30 cycles of denaturation at 94°C for 60 s, annealing at 44–52°C (depending on the primer T<sub>m</sub>) for 60 s, and an extension at 72°C for 1–2 min in an Eppendorf Mastercycler ep gradient thermocycler. Each amplified product was analyzed by 1% agarose gel electrophoresis with ethidium bromide staining. In addition, each PCR product obtained was cloned into the pGEM-T vector (Promega, Madison, WI) and sequenced at Purdue University's Low-Throughput DNA Sequencing Laboratory to confirm that they were the expected loci of each respective gene region.</p> <p>The transcriptional start site of the <italic>so2426 </italic>gene was localized using the 5' RACE System for Rapid Amplification of cDNA Ends version 2.0 (Invitrogen) according to the manufacturer's instructions. Briefly, cDNA was generated using Superscript™ II reverse transcriptase (42°C, 1 h; Invitrogen) in a reaction containing 2 μg of total cellular RNA (from non-stressed or chromate-stressed MR-1 cells) and a <italic>so2426</italic>-specific primer (2426-GSP1: 5'-TAACCGATTGAATTGTT-3'). Nested PCR was performed using two <italic>so2426</italic>-specific primers (2426-GSP2: 5'-TGCACTAACCGTCGCTCTATGGCCTGCAAA-3' and 2426-GSP3: 5'-GCACCAAGTTCATAGCTTCGTATCCTGTCT-3') and the manufacturer-supplied abridged anchor and abridged universal anchor primers. The 5'-RACE products were analyzed by agarose gel electrophoresis and cloned into the pGEM-T vector (Promega) prior to DNA sequencing.</p>
<title>Results and discussion</title> <title>Sequence analysis of the response regulator SO2426</title> <p>Based on its nucleotide sequence, the <italic>so2426 </italic>gene of <italic>S. oneidensis </italic>MR-1 is predicted to encode a 237-amino-acid protein with a calculated molecular mass of 27.4 kDa [##REF##12368813##10##]. Comparative sequence analysis using ClustalW [##REF##7984417##20##] revealed that SO2426 shares sequence similarity to CpxR from <italic>Vibrio cholerae </italic>(36% identity) and <italic>E. coli </italic>(34%), and to OmpR from <italic>V</italic>. <italic>cholerae </italic>(29%) and <italic>E. coli </italic>(27%) (Figure ##FIG##0##1##). In addition, SO2426 orthologs are present and conserved among many of the <italic>Shewanella </italic>species sequenced to date, including <italic>Shewanella </italic>sp. MR-7 (88% sequence identity), <italic>Shewanella </italic>sp. MR-4 and sp. ANA-3 (87%), <italic>S. baltica </italic>OS195 (74%), <italic>S. putrefaciens </italic>200 and CN-32 (73%), <italic>S. woodyi </italic>ATCC 51908 (60%), <italic>S. amazonensis </italic>SB2B (61%), and <italic>S. frigidimarina </italic>NCIMB 400 (49%).</p> <p>The deduced SO2426 protein consists of an N-terminal CheY-like receiver domain and a C-terminal winged-helix DNA-binding domain. The N-terminal portion of SO2426 (amino acids 13–124; see Figure ##FIG##0##1##) shares homology with receiver domains and contains the highly conserved signature residues that are thought to constitute the phosphorylated acid-pocket active site [##REF##10966457##19##,##REF##2689446##21##]. In CheY, these residues correspond to Asp12, Asp13, and Asp57, with Thr87 and Lys109 completing the cluster of conserved residues surrounding the phosphorylation site [##REF##10966457##19##]. Thr87 and Lys109, although not absolutely required for phosphorylation, have been implicated in contributing to the phosphorylation-induced conformational change [##REF##9657998##22##, ####REF##1902474##23##, ##REF##9030562##24####9030562##24##]. <italic>S. oneidensis </italic>SO2426 contains three contiguous aspartate residues near its N-terminus at positions 18, 19 and 20, and there is a conserved substitution to serine at the equivalent position of 87. The phosphate-accepting aspartic residue Asp57 and Lys109 (relative to CheY residue numbering) are retained at equivalent positions in the SO2426 primary sequence. It is presently not known whether a cognate histidine kinase interacts with SO2426 or whether phosphorylation plays a regulatory role in the cellular activity of SO2426. The fact that SO2426 contains the highly conserved residues constituting the active phosphorylation pocket (D18, D19, D62, and K109) suggests that the protein might be differentially controlled by phosphorylation. Furthermore, it remains to be determined whether the regulatory activity of SO2426 is modulated in response to environmental or intracellular stimuli.</p> <p>The C-terminal portion of SO2426 (amino acids 158–235; see Figure ##FIG##0##1##) contains a predicted winged helix-turn-helix (wHTH) motif indicative of the DNA-binding domains of response regulators in the OmpR/PhoB subfamily [##REF##8374950##25##, ####REF##11934608##26##, ##REF##9199401##27####9199401##27##]. The basic structure of the carboxy-terminal domains of OmpR/PhoB subfamily response regulators is characterized by an amino-terminal four-stranded β sheet, a central three-helical bundle, and a C-terminal β-strand hairpin (the wing), which provides an additional interface for DNA contact [##REF##11934608##26##].</p> <title>The <italic>so2426 </italic>gene is co-transcribed in an operon</title> <p>Based on the MR-1 genome annotation, SO2426 is an apparent orphan response regulator given that the <italic>so2426 </italic>locus is not linked with a gene that encodes a potential cognate histidine kinase. As schematically represented in Figure ##FIG##1##2A##, <italic>so2426 </italic>is tightly clustered with four functionally unknown downstream genes, all possessing the same transcriptional polarity. Located immediately downstream of <italic>so2426 </italic>with a short intergenic spacer of 10 bp is a small ORF (132 bp) predicted to encode a soluble hypothetical protein (SO2425), which was shown previously to be induced at the transcriptional level in response to acute chromate challenge [##REF##16524964##11##]. The gene cluster is also characterized by ORFs encoding a zinc carboxypeptidase domain protein (SO2424), a hypothetical protein (SO2423), and a conserved hypothetical protein (SO2422) (Figure ##FIG##1##2A##). The <italic>so2423 </italic>ORF overlaps <italic>so2422 </italic>by 3 bp based on the MR-1 genome annotation (<ext-link ext-link-type="uri" xlink:href="http://www.tigr.org"/>). The <italic>so2426 </italic>region also includes a gene encoding a putative TonB-dependent receptor (SO2427) positioned approximately 170 bp upstream of <italic>so2426</italic>. Putative ρ-independent transcription terminators were identified downstream of the <italic>so2427 </italic>and <italic>so2422 </italic>ORFs (<ext-link ext-link-type="uri" xlink:href="http://cmr.tigr.org"/>; [##UREF##4##28##]).</p> <p>Reverse transcription-PCR (RT-PCR) experiments were performed to determine whether the <italic>so2427</italic>-<italic>so2422 </italic>gene cluster is transcribed as a polycistronic mRNA. Primers P1–P9 (Figure ##FIG##1##2A##) used in the PCR reactions following reverse transcription with random hexamers were designed to amplify regions spanning two or more ORFs in order to establish the transcriptional organization. To assess the reproducibility of the results, multiple independent RT-PCR experiments were carried out using different preparations of total RNA isolated from untreated or chromate-treated MR-1 cells. Figure ##FIG##1##2B## shows a representative gel of the RT-PCR results using the entire suite of primers. Each PCR product was confirmed by DNA sequencing. No products were amplified when using either total RNA as the template for PCR amplification (Figure ##FIG##1##2B##) or primers complementary to an intergenic region located outside of the <italic>so2427</italic>-<italic>so2422 </italic>cluster (results not shown), thus indicating the absence of contaminating genomic DNA in the total RNA preparations. <italic>S. oneidensis </italic>MR-1 genomic DNA was used as a positive control for the PCR conditions (results not shown). Amplicons of the expected sizes, <italic>i</italic>.<italic>e</italic>., 0.85 kb for primers P1/P5, 1.9 kb for primers P1/P6, 2.3 kb for primers P1/P7, 2.5 kb for primers P2/P8, and 3.3 kb for primers P9/P3, were obtained (Figure ##FIG##1##2B##), indicating that the <italic>so2427</italic>-<italic>so2422 </italic>region is co-transcribed as a polycistronic mRNA under no-metal and chromate-amended conditions. However, given the presence of a putative ρ-independent terminator immediately downstream of <italic>so2427</italic>, we cannot rule out the possibility that the product generated from the <italic>so2427</italic>-<italic>so2426 </italic>transcript might be due to read-through transcription.</p> <title>Characterization of the <italic>so2426 </italic>promoter region</title> <p>The transcription start site of the <italic>so2426 </italic>gene was localized using 5' RACE analysis of mRNA transcripts isolated from wild-type cells grown in either the absence or presence of chromate. Two clear 5' termini, both mapping to A residues and corresponding to positions 25 and 27 bp upstream of a second in-frame ATG (amino acid residue M11) (Figure ##FIG##2##3##), were obtained by DNA sequencing of the 5' RACE amplicons. These sites are located at positions 4 and 6 bp downstream of the initial ATG (M1) start codon, indicating that the original sequence annotation for this gene is incorrect and that the second ATG (M11) most likely constitutes the translation start codon for SO2426. The <italic>so2426 </italic>transcription start point indicated potential -10 and -35 basal promoter elements, 5'-TAatAT and 5'-gTGACA, respectively, with a 15-bp spacing between them. The putative -10 promoter sequence (5'-TAatAT) has a 4/6 match to the consensus σ<sup>70</sup>-dependent -10 hexamer element (5'-TATAAT) from <italic>E. coli</italic>. Interestingly, located one base upstream of the -10 hexamer is the extended -10 element 5'-TaTG-3' (Figure ##FIG##2##3##), identified by Mitchell <italic>et al</italic>. [##REF##12907708##29##] and shown to be an important determinant for promoter activity. Finally, the identified <italic>so2426 </italic>promoter region contains the putative -35 element 5'-gTGACA, which is a strong match (5/6) to the <italic>E. coli </italic>consensus -35 hexamer (5'-TTGACA).</p> <p>Further analysis of the <italic>so2426 </italic>upstream regulatory region revealed the presence of a putative recognition site for the ferric uptake regulator (Fur), the iron-binding global transcriptional regulator that serves as the dominant sensor of Fe availability in both gram-positive and gram-negative bacteria (for a review, see [##REF##12829269##30##]). A putative palindromic Fur box, 5'-AAATGAtATTgATTcTCgTTT-3', which closely matched our previously derived predicted Fur-binding motif for <italic>S. oneidensis </italic>MR-1 [##REF##15576789##31##], was identified at positions -5 to +16 immediately downstream of the -10 promoter element and overlapped the mapped transcription start sites for <italic>so2426 </italic>(Figure ##FIG##2##3##). Identification of a Fur box is consistent with our previous microarray studies demonstrating increased <italic>so2426 </italic>mRNA levels in a <italic>S. oneidensis </italic>MR-1 <italic>fur </italic>knockout strain [##REF##15576789##31##] and suggests that expression of the SO2426 response regulator might be subject to Fur-mediated repression.</p> <title>Growth and chromate reduction kinetics of the <italic>so2426 </italic>deletion mutant</title> <p>A <italic>cre</italic>-<italic>lox </italic>recombination method [##REF##16391095##32##,##REF##12449384##33##] was used to create a nonpolar in-frame deletion of the <italic>so2426 </italic>locus in MR-1 for the purpose of analyzing in depth the functional role of this response regulator in metal stress responses (see Methods for a detailed description). The deleted segment was 721 bp in length and corresponded to positions -77 to +645 (<italic>S. oneidensis </italic>MR-1 genome coordinates 2534690 and 2533969) from the annotated translational start site, leaving only 69 nucleotides of the <italic>so2426 </italic>gene intact at the C-terminal end. The growth behavior of the Δ<italic>so2426 </italic>mutant in response to Cr(VI) was examined using a Bioscreen C reader to monitor culture turbidity (OD<sub>600</sub>) at 30-min intervals under aerobic conditions over 72 h. As shown in Figure ##FIG##3##4A##, growth of the <italic>S. oneidensis </italic>Δ<italic>so2426 </italic>strain was comparable to the wild type in LB broth without added chromate. By contrast, growth of the Δ<italic>so2426 </italic>mutant was inhibited to a greater extent than the wild type in the presence of 0.3 mM chromate as evident from the longer time interval required before the commencement of exponential growth.</p> <p>To determine whether deletion of <italic>so2426 </italic>affected the ability of MR-1 to remove extracellular Cr(VI) from the medium, the wild-type and mutant strains were challenged with 0.3 mM chromate (the same final concentration used for the transcriptomic studies) when cells reached mid-log phase (OD<sub>600</sub>, 0.5). Residual hexavalent chromium was measured spectrophotometrically at 3-h intervals using the chromogenic 1,5-diphenylcarbazide (DPC) method [##REF##10788340##34##]. The chromate disappearance profiles in Figure ##FIG##3##4B## show that the parental MR-1 strain removed 100% of the external Cr(VI) within 9 h, while the Δ<italic>so2426 </italic>mutant transformed approximately 50% within the same time period. No abiotic conversion of chromate was observed for the cell-free LB control (data not shown).</p> <p>Complementation experiments were performed to ascertain whether the enhanced chromate sensitivity and slower chromate transformation rate for the <italic>so2426 </italic>null mutant were due to the deletion of the response regulator gene <italic>so2426 </italic>and not to secondary-site mutations introduced during the mutagenesis procedure. A complemented Δ<italic>so2426 </italic>strain was constructed by cloning the full-length <italic>so2426 </italic>wild-type ORF plus 100 bp of contiguous upstream sequence into the broad-host-range vector pBBR2426 as described in the Methods. Transformation of pBBR2426 into the Δ<italic>so2426 </italic>mutant resulted in a complemented strain, designated Δ<italic>so2426</italic>(<italic>so2426</italic><sup>+</sup>). Genetic complementation of the <italic>so2426 </italic>gene in the mutant strain restored growth (data not shown) and Cr(VI) reduction capacity to near wild-type levels (Figure ##FIG##3##4B##), thus reaffirming the integral role of SO2426 in the cellular response to chromate.</p> <title>Temporal transcriptome analysis of Δ<italic>so2426 </italic>and wild-type strains</title> <p>Prior to this study, nothing was known about the regulatory targets of the OmpR family response regulator SO2426. To investigate the effect of the <italic>so2426 </italic>deletion on global gene expression and to identify candidate target genes under SO2426 control, the chromate-perturbed transcriptomes of the Δ<italic>so2426 </italic>mutant and the <italic>S. oneidensis </italic>MR-1 wild type were compared using whole-genome microarray analysis. Time-series microarray experiments were performed in which global transcriptional profiles of the wild type and the Δ<italic>so2426 </italic>mutant strains were compared at 5, 30, 60, 90, 180 min, and 24 hours following the addition of K<sub>2</sub>CrO<sub>4 </sub>(0.3 mM final concentration) to mid-exponentially growing cells (see Additional file ##SUPPL##0##1## for a complete list of the processed microarray data). Genes with altered expression ratios of ≥ 2-fold for at least one of the time points and having a <italic>p </italic>value of &lt; 0.05 were considered to be significantly changed and were further analyzed using clustering methods. Overall, microarray expression profiling revealed that 82 (5 min), 90 (30 min), 80 (60 min), 109, (90 min), 443 (180 min), and 465 (24 h) genes were induced in the Δ<italic>so2426 </italic>mutant under chromate exposure, whereas 125 (5 min), 81 (30 min), 56 (60 min), 105 (90 min), 314 (180 min), and 364 (24 h) genes were repressed. The total number of differentially expressed genes was substantially less at the shorter, acute exposure times (<italic>i.e</italic>., 5, 30, 60, and 90 min) compared to the longer exposure periods (<italic>i.e</italic>., 180 min and 24 h), at which point other effects such as accumulation of intracellular chromium (hexavalent and/or reduced forms) and stationary phase growth are likely impacting global expression profiles.</p> <p>Quantitative real-time RT-PCR (qRT-PCR) analysis was employed to assess the general quality of the microarray data by providing an independent measurement of gene expression for a small subset of genes. The five genes (<italic>so0404</italic>, <italic>so2426</italic>, <italic>so3670</italic>, <italic>so1826</italic>, and <italic>so3585</italic>) selected for comparative qRT-PCR displayed downregulated or unaltered expression patterns at the 180-min time point of chromate exposure based on microarray hybridization. Comparison of gene expression levels as determined by microarray and qRT-PCR analyses indicated a high degree of concordance, with a Pearson's correlation coefficient (<italic>r</italic>) value of 0.99 (data not shown). For the mutant strain, both microarray hybridization and qRT-PCR indicated that <italic>so2426 </italic>expression levels were not detectable above background, thus providing additional confirmation of deletion of the <italic>so2426 </italic>gene.</p> <title>Functional classification of genes differentially expressed in the Δ<italic>so2426 </italic>mutant</title> <p>Pairwise complete-linkage clustering was used to identify groups of potentially co-ordinately regulated genes among the 841 ORFs (18% of the arrayed genome) showing a twofold or greater statistically significant (<italic>p </italic>&lt; 0.05) change in expression for at least one time point during Cr(VI) exposure. One particularly noteworthy cluster of similarly expressed genes (Cluster A in Figure ##FIG##4##5##) comprises 46 genes with functions distributed across the following role categories: cellular processes (2 genes), energy metabolism (9 genes), hypothetical proteins (17 genes), protein synthesis (1 gene), regulatory functions/signal transduction (4 genes), and transport and binding proteins (13 genes). These clustered genes displayed a distinct pattern of downregulated expression in the Δ<italic>so2426 </italic>mutant relative to the wild-type MR-1 strain, particularly at the earlier time points (5, 30, 60, and 90 min) of chromate exposure (see Table ##TAB##0##1##). This is consistent with our previous studies demonstrating that the wild-type <italic>so2426 </italic>gene was induced at the transcriptional level early in Cr(VI) exposure, <italic>i.e</italic>. 90 min or less, but not at 24 h post-treatment [##REF##16524964##11##,##REF##16957260##12##]. Closer inspection of the 46 potentially co-regulated genes (Table ##TAB##0##1##) revealed that many of the transport/binding, hypothetical, and regulatory genes were shown previously to be upregulated in wild-type MR-1 cells in response to chromate treatment [##REF##16524964##11##,##REF##17385904##13##], thus suggesting that these ORFs are candidate gene targets of positive control by SO2426. The genes constituting Cluster A are discussed in greater detail in the sections below.</p> <p>Cluster B (Figure ##FIG##4##5##) includes 20 genes, three of which encode AcrB/AcrD/AcrF family proteins (SO1882, SO3484, SO4692) presumably involved in drug efflux and resistance. All of these 20 genes showed decreased expression in the Δ<italic>so2426 </italic>mutant except at the 24-h time point, where mRNA levels increased (Figure ##FIG##4##5## and Additional file ##SUPPL##0##1##). Notably, Cluster B also contains genes <italic>so0392 </italic>(putative lipoprotein) and <italic>so4688 </italic>(glycosyl transferase, group 2 family protein), both with annotated functions associated with the cell envelope, as well as <italic>cpxA </italic>(<italic>so4478</italic>). In <italic>E. coli</italic>, CpxA is the sensory inner membrane kinase that functions with its cognate partner, a cytoplasmic response regulator (CpxR), in a prototypical two-component system to sense and respond to perturbations in the bacterial cell envelope (reviewed in [##REF##15802241##35##]). <italic>S. oneidensis cpxR </italic>also displayed downregulated temporal expression profiles in the Δ<italic>so2426 </italic>mutant under chromate challenge but grouped with genes in Cluster A.</p> <title>Genes encoding transport and binding functions</title> <p>Thirteen of the genes in Cluster A have annotated functions in metal transport and binding, in particular Fe acquisition and homeostasis: <italic>ftn </italic>(ferritin), <italic>so1580 </italic>(TonB-dependent heme receptor), <italic>so1771 </italic>(GntP family permease), <italic>so2045 </italic>(cation efflux family protein), <italic>alcA </italic>(siderophore biosynthesis protein), <italic>so3031 </italic>and <italic>so3032 </italic>(both putative siderophore biosynthesis proteins encoded immediately downstream of <italic>alcA</italic>), <italic>so3033 </italic>(ferric alcaligin siderophore receptor), <italic>so3063 </italic>(sodium:alanine symporter family protein), <italic>so4150 </italic>(putative transporter), <italic>viuA </italic>(ferric vibriobactin receptor), <italic>so4712 </italic>(putative ABC transporter, ATP-binding protein), and <italic>so4743 </italic>(putative TonB-dependent receptor). In contrast to the TonB-dependent receptor genes <italic>so1580 </italic>and <italic>so4743</italic>, either no significant change or a maximal twofold decrease at the 180-min time point was observed for <italic>so2427 </italic>(putative TonB-dependent receptor), located directly upstream of the <italic>so2426 </italic>response regulator gene. Of the 13 genes with transport and binding functions in Cluster A, ten ORFs (<italic>ftn</italic>, <italic>so1580</italic>, <italic>so2045</italic>, <italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032</italic>, <italic>so3033</italic>, <italic>so3063</italic>, <italic>viuA</italic>, <italic>so4743</italic>) were shown in an earlier study to be induced in chromate-challenged wild-type MR-1 cells compared to untreated cells [##REF##16524964##11##,##REF##17385904##13##]. Temporal expression patterns for these genes demonstrated that they were downregulated 2- to 20-fold in the Δ<italic>so2426 </italic>strain over the 180-min time course (Table ##TAB##0##1##), suggesting that SO2426 acts as a direct or indirect positive regulator of a subset of Fe uptake and storage genes in <italic>S. oneidensis</italic>.</p> <p>The <italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032 </italic>operon encodes proteins required for siderophore biosynthesis in MR-1 [##REF##16391131##36##], and the product of gene <italic>so3033 </italic>is predicted to allow for the cellular utilization of the structurally undetermined MR-1 siderophore. The differential profiles for these four genes were characterized by a peak in down-regulated expression (ranging from ~7- to 17-fold) at the 180-min time interval post chromate addition, followed by no significant change in expression at 24 h (Table ##TAB##0##1##). The first gene in the siderophore biosynthetic operon shows 48% sequence identity to <italic>Bordetella pertussis alcA</italic>, which is required for alcaligin production [##REF##8759851##37##], while the two downstream genes (<italic>so3031 </italic>and <italic>so3032</italic>, respectively) diverge from <italic>alcB </italic>and <italic>alcC </italic>and are annotated as putative siderophore biosynthesis genes. While the structural identity of the MR-1 siderophore has not been elucidated, we predict that MR-1 likely produces a siderophore similar to putrebactin [##UREF##5##38##], a novel cyclic dihydroxamate siderophore characterized from <italic>Shewanella putrefaciens </italic>strain 200 and structurally similar to alcaligin.</p> <p>Consistent with the transcriptomic data, further physiological evidence for the involvement of the SO2426 response regulator in controlling siderophore-dependent iron acquisition was obtained by performing semiquantitative liquid CAS assays in which relative siderophore production levels in supernatants of the Δ<italic>so2426 </italic>mutant were compared to those in wild-type MR-1 cultures. Following 24-h growth in LB medium in the absence of added FeCl<sub>3 </sub>or chromate, the Δ<italic>so2426 </italic>mutant (<italic>A</italic><sub>630</sub>, 0.739 ± 0.02) exhibited essentially no detectable siderophore production over that of the cell-free control (<italic>A</italic><sub>630</sub>, 0.726 ± 0.005), whereas wild-type MR-1 cultures (<italic>A</italic><sub>630</sub>, 0.108 ± 0.06) produced approximately 7-fold more siderophore than the Δ<italic>so2426 </italic>mutant. As expected, addition of 50 μM FeCl<sub>3 </sub>to the culture medium reduced siderophore production by the wild-type to near background levels (an approximately 6-fold reduction) but did not affect siderophore accumulation in the Δ<italic>so2426 </italic>mutant. In the presence of 0.3 mM chromate, replicate cultures of wild-type MR-1 exhibited increasing levels in relative siderophore production over time (Figure ##FIG##5##6##), which was consistent with previous microarray expression data showing induction of siderophore biosynthesis genes in response to chromate exposure [##REF##16524964##11##]. By contrast, the temporal profile for the Δ<italic>so2426 </italic>mutant showed an initial moderate increase and then a dramatic reduction in siderophore excretion compared to the wild type (Figure ##FIG##5##6##). Collectively, these data support the hypothesis that SO2426 is a positive regulator of siderophore-mediated Fe uptake.</p> <p>Cluster A also includes an <italic>S. oneidensis </italic>gene that is involved in Fe storage, <italic>i.e</italic>., the <italic>ftn </italic>gene (<italic>so0139</italic>) encoding ferritin. The <italic>ftn </italic>mRNA expression levels were down-regulated 4- to 20-fold in the SO2426 null mutant across the entire temporal range, with the peak in repression occurring at the 180-min time point (Table ##TAB##0##1##). Although iron is an essential micronutrient for most organisms, free excess Fe in the presence of oxygen is potentially toxic because of its tendency to participate in Fenton chemistry reactions, which generate cell-damaging reactive oxygen intermediates. By scavenging and storing intracellular free Fe in a remobilizable form, Ftn protects against Fe(II)-mediated formation of hydroxyl radicals [##REF##9889981##39##,##REF##10620317##40##]. Studies have indicated that ferritins make important contributions to both Fe storage as well as cellular protection against oxidative stress [##REF##9782498##41##,##REF##8809765##42##] and that <italic>ftn </italic>transcription is co-ordinately regulated by iron and redox stress in some bacteria [##REF##15256555##43##]. <italic>S. oneidensis ftn </italic>was shown previously to be induced at a low level (~2-fold) in chromate-exposed wild-type MR-1 cells in contrast to untreated cells [##REF##16524964##11##]. The induction of <italic>ftn </italic>in <italic>S. oneidensis </italic>might act, at least partially, to mitigate chromate-generated oxidative stress by sequestering excess free Fe available for the catalysis of Fenton reactions and the formation of additional reactive hydroxyl radicals.</p> <title>Genes of unknown function</title> <p>The majority of coordinately repressed genes in Cluster A belong to the functional category of hypothetical and conserved hypothetical proteins (Table ##TAB##0##1##). Most notably, genes comprising the conserved hypothetical <italic>so1188</italic>-<italic>so1189</italic>-<italic>so1190 </italic>operon exhibited the highest degree of down-regulated expression in the chromate-exposed Δ<italic>so2426 </italic>mutant relative to the wild type, with decreases in mRNA ratios ranging from 45- to &gt; 300-fold over the 180-min period. No significant changes in expression, however, were observed for these genes at the 24-h time point. In addition, the hypothetical gene <italic>so2425</italic>, which is co-transcribed with <italic>so2426 </italic>in a polycistronic transcript (shown in this study), was repressed as much as 11-fold over the 180-min period, pointing to the possibility that the <italic>so2426 </italic>operon may be under positive autoregulatory control. Surprisingly, the functionally unknown genes <italic>so2424</italic>, <italic>so2423</italic>, and <italic>so2422</italic>, which also form part of the <italic>so2426 </italic>operon, showed either a less than 2-fold change or no significant differential expression over the time course.</p> <title>Regulatory genes</title> <p>In addition to metal transport and hypothetical genes described above, the temporal expression profiles for four transcriptional regulatory genes (<italic>so0916</italic>, <italic>so0544</italic>, <italic>cpxR</italic>, <italic>so4567</italic>) were altered in the Δ<italic>so2426 </italic>mutant. The function of the regulators encoded by <italic>so0916</italic>, <italic>so0544</italic>, and <italic>so4567 </italic>in <italic>S. oneidensis </italic>as well as their target genes are not known. However, gene <italic>so4477 </italic>encodes a homolog of CpxR, the response regulator component of the CpxAR envelope stress response system, which has been well characterized in <italic>E. coli </italic>(reviewed in [##REF##15802241##35##,##REF##16487683##44##]). The annotated <italic>cpxR </italic>gene in <italic>S. oneidensis </italic>MR-1 was downregulated as much as threefold in response to the <italic>so2426 </italic>deletion under chromate stress conditions (Table ##TAB##0##1##). This is in contrast to the two- to three-fold induction measured for <italic>cpxR </italic>in chromate-treated wild-type MR-1 cells [##REF##16524964##11##]. Located immediately downstream of <italic>cpxR </italic>in an apparent operon is <italic>cpxA</italic>, which was downregulated as much as 3-fold in the mutant but grouped in a separate cluster (Cluster B in Figure ##FIG##4##5##).</p> <p>In <italic>E. coli</italic>, the CpxAR two-component pair constitutes a stress response system the main function of which is to sense and respond to cell envelope distress, principally misfolded periplasmic proteins, by activating genes encoding proteases and folding catalysts [##REF##15802241##35##,##REF##7883164##45##,##REF##11544368##46##]. However, the complexity, extent, and response overlap of the Cpx system have become progressively apparent, as evident from its involvement in outer membrane porin expression [##REF##16077119##47##], stationary-phase survival [##REF##10542180##48##], and the response to high pH stress [##REF##9473036##49##] to name a few. Although the specific cellular role of the CpxR homolog in <italic>S. oneidensis </italic>has not been established, it is intriguing to consider the possibility that interplay between CpxR and SO2426 may occur in response to metal stress. It is conceivable that the CpxA-CpxR two-component system may respond to oxidative damage imposed on biomolecules as a consequence of reactive oxygen species generated during intracellular partial Cr(VI) reduction.</p> <title>Other differentially regulated genes in the <italic>so2426 </italic>deletion mutant</title> <p>While hypothetical and metal transport/binding genes are dominant in Cluster A (representing 65%), other notable genes with ≥ 2-fold-repressed mRNA levels in the Δ<italic>so2426 </italic>mutant include ORFs codifying a putative bicyclomycin resistance (Bcr) protein (SO2280) and a ferric iron reductase protein (SO3034). Bcr proteins are members of the major facilitator superfamily class of membrane efflux pumps and play important roles in drug resistance [##REF##11104814##50##]. In a previous work, the <italic>so2280 </italic>gene in a wild-type MR-1 genetic background exhibited 2- to 4-fold upregulated expression at the transcript level in response to short-term chromate exposure, <italic>i.e</italic>., up to 90 min [##REF##16524964##11##], but was shown in this study to be repressed as much as 4-fold in Cr(VI)-challenged Δ<italic>so2426 </italic>cells (Table ##TAB##0##1##), suggesting that the <italic>so2280 </italic>gene is a direct or indirect target of SO2426-mediated activation under chromate conditions. Although the cellular function of SO2280 is presently unknown, its co-ordinately regulated expression with the siderophore biosynthesis operon (<italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032</italic>) and gene <italic>so3033</italic>, encoding a ferric alcaligin siderophore receptor, is intriguing. Studies in other bacterial systems have shown that certain proteins resembling Bcr at the amino acid level function in siderophore export from the cell [##REF##15901687##51##].</p> <p>Downregulated expression of gene <italic>so3034 </italic>(a putative ferric iron reductase) peaked at the 90-min time point with a 4-fold decrease in mRNA levels (Table ##TAB##0##1##). Ferric iron reductases are thought to be involved in reducing siderophore-complexed Fe(III) to Fe(II), thus enabling dissociation of the high-affinity Fe<sup>3+</sup>-chelating siderophore from iron [##REF##12829269##30##]. The entire <italic>so3030</italic>-<italic>so3034 </italic>gene cluster, which is predicted to be involved in siderophore production and utilization, exhibits similar temporal expression profiles in the mutant, suggesting a role for the SO2426 response regulator in coordinating these functions. It is not clear whether the other differentially expressed energy metabolism genes (<italic>fdhG</italic>, <italic>fdnH</italic>, <italic>fdhE</italic>, <italic>so4503</italic>, <italic>so4506</italic>, <italic>so4509</italic>) contribute directly to chromate stress response or show altered transcription as a result of secondary effects associated with the <italic>so2426 </italic>deletion. These genes were not observed to be induced in response to chromate exposure in our previous transcriptomic studies using wild-type <italic>S. oneidensis </italic>MR-1 [##REF##16524964##11##].</p>
<title>Results and discussion</title> <title>Sequence analysis of the response regulator SO2426</title> <p>Based on its nucleotide sequence, the <italic>so2426 </italic>gene of <italic>S. oneidensis </italic>MR-1 is predicted to encode a 237-amino-acid protein with a calculated molecular mass of 27.4 kDa [##REF##12368813##10##]. Comparative sequence analysis using ClustalW [##REF##7984417##20##] revealed that SO2426 shares sequence similarity to CpxR from <italic>Vibrio cholerae </italic>(36% identity) and <italic>E. coli </italic>(34%), and to OmpR from <italic>V</italic>. <italic>cholerae </italic>(29%) and <italic>E. coli </italic>(27%) (Figure ##FIG##0##1##). In addition, SO2426 orthologs are present and conserved among many of the <italic>Shewanella </italic>species sequenced to date, including <italic>Shewanella </italic>sp. MR-7 (88% sequence identity), <italic>Shewanella </italic>sp. MR-4 and sp. ANA-3 (87%), <italic>S. baltica </italic>OS195 (74%), <italic>S. putrefaciens </italic>200 and CN-32 (73%), <italic>S. woodyi </italic>ATCC 51908 (60%), <italic>S. amazonensis </italic>SB2B (61%), and <italic>S. frigidimarina </italic>NCIMB 400 (49%).</p> <p>The deduced SO2426 protein consists of an N-terminal CheY-like receiver domain and a C-terminal winged-helix DNA-binding domain. The N-terminal portion of SO2426 (amino acids 13–124; see Figure ##FIG##0##1##) shares homology with receiver domains and contains the highly conserved signature residues that are thought to constitute the phosphorylated acid-pocket active site [##REF##10966457##19##,##REF##2689446##21##]. In CheY, these residues correspond to Asp12, Asp13, and Asp57, with Thr87 and Lys109 completing the cluster of conserved residues surrounding the phosphorylation site [##REF##10966457##19##]. Thr87 and Lys109, although not absolutely required for phosphorylation, have been implicated in contributing to the phosphorylation-induced conformational change [##REF##9657998##22##, ####REF##1902474##23##, ##REF##9030562##24####9030562##24##]. <italic>S. oneidensis </italic>SO2426 contains three contiguous aspartate residues near its N-terminus at positions 18, 19 and 20, and there is a conserved substitution to serine at the equivalent position of 87. The phosphate-accepting aspartic residue Asp57 and Lys109 (relative to CheY residue numbering) are retained at equivalent positions in the SO2426 primary sequence. It is presently not known whether a cognate histidine kinase interacts with SO2426 or whether phosphorylation plays a regulatory role in the cellular activity of SO2426. The fact that SO2426 contains the highly conserved residues constituting the active phosphorylation pocket (D18, D19, D62, and K109) suggests that the protein might be differentially controlled by phosphorylation. Furthermore, it remains to be determined whether the regulatory activity of SO2426 is modulated in response to environmental or intracellular stimuli.</p> <p>The C-terminal portion of SO2426 (amino acids 158–235; see Figure ##FIG##0##1##) contains a predicted winged helix-turn-helix (wHTH) motif indicative of the DNA-binding domains of response regulators in the OmpR/PhoB subfamily [##REF##8374950##25##, ####REF##11934608##26##, ##REF##9199401##27####9199401##27##]. The basic structure of the carboxy-terminal domains of OmpR/PhoB subfamily response regulators is characterized by an amino-terminal four-stranded β sheet, a central three-helical bundle, and a C-terminal β-strand hairpin (the wing), which provides an additional interface for DNA contact [##REF##11934608##26##].</p> <title>The <italic>so2426 </italic>gene is co-transcribed in an operon</title> <p>Based on the MR-1 genome annotation, SO2426 is an apparent orphan response regulator given that the <italic>so2426 </italic>locus is not linked with a gene that encodes a potential cognate histidine kinase. As schematically represented in Figure ##FIG##1##2A##, <italic>so2426 </italic>is tightly clustered with four functionally unknown downstream genes, all possessing the same transcriptional polarity. Located immediately downstream of <italic>so2426 </italic>with a short intergenic spacer of 10 bp is a small ORF (132 bp) predicted to encode a soluble hypothetical protein (SO2425), which was shown previously to be induced at the transcriptional level in response to acute chromate challenge [##REF##16524964##11##]. The gene cluster is also characterized by ORFs encoding a zinc carboxypeptidase domain protein (SO2424), a hypothetical protein (SO2423), and a conserved hypothetical protein (SO2422) (Figure ##FIG##1##2A##). The <italic>so2423 </italic>ORF overlaps <italic>so2422 </italic>by 3 bp based on the MR-1 genome annotation (<ext-link ext-link-type="uri" xlink:href="http://www.tigr.org"/>). The <italic>so2426 </italic>region also includes a gene encoding a putative TonB-dependent receptor (SO2427) positioned approximately 170 bp upstream of <italic>so2426</italic>. Putative ρ-independent transcription terminators were identified downstream of the <italic>so2427 </italic>and <italic>so2422 </italic>ORFs (<ext-link ext-link-type="uri" xlink:href="http://cmr.tigr.org"/>; [##UREF##4##28##]).</p> <p>Reverse transcription-PCR (RT-PCR) experiments were performed to determine whether the <italic>so2427</italic>-<italic>so2422 </italic>gene cluster is transcribed as a polycistronic mRNA. Primers P1–P9 (Figure ##FIG##1##2A##) used in the PCR reactions following reverse transcription with random hexamers were designed to amplify regions spanning two or more ORFs in order to establish the transcriptional organization. To assess the reproducibility of the results, multiple independent RT-PCR experiments were carried out using different preparations of total RNA isolated from untreated or chromate-treated MR-1 cells. Figure ##FIG##1##2B## shows a representative gel of the RT-PCR results using the entire suite of primers. Each PCR product was confirmed by DNA sequencing. No products were amplified when using either total RNA as the template for PCR amplification (Figure ##FIG##1##2B##) or primers complementary to an intergenic region located outside of the <italic>so2427</italic>-<italic>so2422 </italic>cluster (results not shown), thus indicating the absence of contaminating genomic DNA in the total RNA preparations. <italic>S. oneidensis </italic>MR-1 genomic DNA was used as a positive control for the PCR conditions (results not shown). Amplicons of the expected sizes, <italic>i</italic>.<italic>e</italic>., 0.85 kb for primers P1/P5, 1.9 kb for primers P1/P6, 2.3 kb for primers P1/P7, 2.5 kb for primers P2/P8, and 3.3 kb for primers P9/P3, were obtained (Figure ##FIG##1##2B##), indicating that the <italic>so2427</italic>-<italic>so2422 </italic>region is co-transcribed as a polycistronic mRNA under no-metal and chromate-amended conditions. However, given the presence of a putative ρ-independent terminator immediately downstream of <italic>so2427</italic>, we cannot rule out the possibility that the product generated from the <italic>so2427</italic>-<italic>so2426 </italic>transcript might be due to read-through transcription.</p> <title>Characterization of the <italic>so2426 </italic>promoter region</title> <p>The transcription start site of the <italic>so2426 </italic>gene was localized using 5' RACE analysis of mRNA transcripts isolated from wild-type cells grown in either the absence or presence of chromate. Two clear 5' termini, both mapping to A residues and corresponding to positions 25 and 27 bp upstream of a second in-frame ATG (amino acid residue M11) (Figure ##FIG##2##3##), were obtained by DNA sequencing of the 5' RACE amplicons. These sites are located at positions 4 and 6 bp downstream of the initial ATG (M1) start codon, indicating that the original sequence annotation for this gene is incorrect and that the second ATG (M11) most likely constitutes the translation start codon for SO2426. The <italic>so2426 </italic>transcription start point indicated potential -10 and -35 basal promoter elements, 5'-TAatAT and 5'-gTGACA, respectively, with a 15-bp spacing between them. The putative -10 promoter sequence (5'-TAatAT) has a 4/6 match to the consensus σ<sup>70</sup>-dependent -10 hexamer element (5'-TATAAT) from <italic>E. coli</italic>. Interestingly, located one base upstream of the -10 hexamer is the extended -10 element 5'-TaTG-3' (Figure ##FIG##2##3##), identified by Mitchell <italic>et al</italic>. [##REF##12907708##29##] and shown to be an important determinant for promoter activity. Finally, the identified <italic>so2426 </italic>promoter region contains the putative -35 element 5'-gTGACA, which is a strong match (5/6) to the <italic>E. coli </italic>consensus -35 hexamer (5'-TTGACA).</p> <p>Further analysis of the <italic>so2426 </italic>upstream regulatory region revealed the presence of a putative recognition site for the ferric uptake regulator (Fur), the iron-binding global transcriptional regulator that serves as the dominant sensor of Fe availability in both gram-positive and gram-negative bacteria (for a review, see [##REF##12829269##30##]). A putative palindromic Fur box, 5'-AAATGAtATTgATTcTCgTTT-3', which closely matched our previously derived predicted Fur-binding motif for <italic>S. oneidensis </italic>MR-1 [##REF##15576789##31##], was identified at positions -5 to +16 immediately downstream of the -10 promoter element and overlapped the mapped transcription start sites for <italic>so2426 </italic>(Figure ##FIG##2##3##). Identification of a Fur box is consistent with our previous microarray studies demonstrating increased <italic>so2426 </italic>mRNA levels in a <italic>S. oneidensis </italic>MR-1 <italic>fur </italic>knockout strain [##REF##15576789##31##] and suggests that expression of the SO2426 response regulator might be subject to Fur-mediated repression.</p> <title>Growth and chromate reduction kinetics of the <italic>so2426 </italic>deletion mutant</title> <p>A <italic>cre</italic>-<italic>lox </italic>recombination method [##REF##16391095##32##,##REF##12449384##33##] was used to create a nonpolar in-frame deletion of the <italic>so2426 </italic>locus in MR-1 for the purpose of analyzing in depth the functional role of this response regulator in metal stress responses (see Methods for a detailed description). The deleted segment was 721 bp in length and corresponded to positions -77 to +645 (<italic>S. oneidensis </italic>MR-1 genome coordinates 2534690 and 2533969) from the annotated translational start site, leaving only 69 nucleotides of the <italic>so2426 </italic>gene intact at the C-terminal end. The growth behavior of the Δ<italic>so2426 </italic>mutant in response to Cr(VI) was examined using a Bioscreen C reader to monitor culture turbidity (OD<sub>600</sub>) at 30-min intervals under aerobic conditions over 72 h. As shown in Figure ##FIG##3##4A##, growth of the <italic>S. oneidensis </italic>Δ<italic>so2426 </italic>strain was comparable to the wild type in LB broth without added chromate. By contrast, growth of the Δ<italic>so2426 </italic>mutant was inhibited to a greater extent than the wild type in the presence of 0.3 mM chromate as evident from the longer time interval required before the commencement of exponential growth.</p> <p>To determine whether deletion of <italic>so2426 </italic>affected the ability of MR-1 to remove extracellular Cr(VI) from the medium, the wild-type and mutant strains were challenged with 0.3 mM chromate (the same final concentration used for the transcriptomic studies) when cells reached mid-log phase (OD<sub>600</sub>, 0.5). Residual hexavalent chromium was measured spectrophotometrically at 3-h intervals using the chromogenic 1,5-diphenylcarbazide (DPC) method [##REF##10788340##34##]. The chromate disappearance profiles in Figure ##FIG##3##4B## show that the parental MR-1 strain removed 100% of the external Cr(VI) within 9 h, while the Δ<italic>so2426 </italic>mutant transformed approximately 50% within the same time period. No abiotic conversion of chromate was observed for the cell-free LB control (data not shown).</p> <p>Complementation experiments were performed to ascertain whether the enhanced chromate sensitivity and slower chromate transformation rate for the <italic>so2426 </italic>null mutant were due to the deletion of the response regulator gene <italic>so2426 </italic>and not to secondary-site mutations introduced during the mutagenesis procedure. A complemented Δ<italic>so2426 </italic>strain was constructed by cloning the full-length <italic>so2426 </italic>wild-type ORF plus 100 bp of contiguous upstream sequence into the broad-host-range vector pBBR2426 as described in the Methods. Transformation of pBBR2426 into the Δ<italic>so2426 </italic>mutant resulted in a complemented strain, designated Δ<italic>so2426</italic>(<italic>so2426</italic><sup>+</sup>). Genetic complementation of the <italic>so2426 </italic>gene in the mutant strain restored growth (data not shown) and Cr(VI) reduction capacity to near wild-type levels (Figure ##FIG##3##4B##), thus reaffirming the integral role of SO2426 in the cellular response to chromate.</p> <title>Temporal transcriptome analysis of Δ<italic>so2426 </italic>and wild-type strains</title> <p>Prior to this study, nothing was known about the regulatory targets of the OmpR family response regulator SO2426. To investigate the effect of the <italic>so2426 </italic>deletion on global gene expression and to identify candidate target genes under SO2426 control, the chromate-perturbed transcriptomes of the Δ<italic>so2426 </italic>mutant and the <italic>S. oneidensis </italic>MR-1 wild type were compared using whole-genome microarray analysis. Time-series microarray experiments were performed in which global transcriptional profiles of the wild type and the Δ<italic>so2426 </italic>mutant strains were compared at 5, 30, 60, 90, 180 min, and 24 hours following the addition of K<sub>2</sub>CrO<sub>4 </sub>(0.3 mM final concentration) to mid-exponentially growing cells (see Additional file ##SUPPL##0##1## for a complete list of the processed microarray data). Genes with altered expression ratios of ≥ 2-fold for at least one of the time points and having a <italic>p </italic>value of &lt; 0.05 were considered to be significantly changed and were further analyzed using clustering methods. Overall, microarray expression profiling revealed that 82 (5 min), 90 (30 min), 80 (60 min), 109, (90 min), 443 (180 min), and 465 (24 h) genes were induced in the Δ<italic>so2426 </italic>mutant under chromate exposure, whereas 125 (5 min), 81 (30 min), 56 (60 min), 105 (90 min), 314 (180 min), and 364 (24 h) genes were repressed. The total number of differentially expressed genes was substantially less at the shorter, acute exposure times (<italic>i.e</italic>., 5, 30, 60, and 90 min) compared to the longer exposure periods (<italic>i.e</italic>., 180 min and 24 h), at which point other effects such as accumulation of intracellular chromium (hexavalent and/or reduced forms) and stationary phase growth are likely impacting global expression profiles.</p> <p>Quantitative real-time RT-PCR (qRT-PCR) analysis was employed to assess the general quality of the microarray data by providing an independent measurement of gene expression for a small subset of genes. The five genes (<italic>so0404</italic>, <italic>so2426</italic>, <italic>so3670</italic>, <italic>so1826</italic>, and <italic>so3585</italic>) selected for comparative qRT-PCR displayed downregulated or unaltered expression patterns at the 180-min time point of chromate exposure based on microarray hybridization. Comparison of gene expression levels as determined by microarray and qRT-PCR analyses indicated a high degree of concordance, with a Pearson's correlation coefficient (<italic>r</italic>) value of 0.99 (data not shown). For the mutant strain, both microarray hybridization and qRT-PCR indicated that <italic>so2426 </italic>expression levels were not detectable above background, thus providing additional confirmation of deletion of the <italic>so2426 </italic>gene.</p> <title>Functional classification of genes differentially expressed in the Δ<italic>so2426 </italic>mutant</title> <p>Pairwise complete-linkage clustering was used to identify groups of potentially co-ordinately regulated genes among the 841 ORFs (18% of the arrayed genome) showing a twofold or greater statistically significant (<italic>p </italic>&lt; 0.05) change in expression for at least one time point during Cr(VI) exposure. One particularly noteworthy cluster of similarly expressed genes (Cluster A in Figure ##FIG##4##5##) comprises 46 genes with functions distributed across the following role categories: cellular processes (2 genes), energy metabolism (9 genes), hypothetical proteins (17 genes), protein synthesis (1 gene), regulatory functions/signal transduction (4 genes), and transport and binding proteins (13 genes). These clustered genes displayed a distinct pattern of downregulated expression in the Δ<italic>so2426 </italic>mutant relative to the wild-type MR-1 strain, particularly at the earlier time points (5, 30, 60, and 90 min) of chromate exposure (see Table ##TAB##0##1##). This is consistent with our previous studies demonstrating that the wild-type <italic>so2426 </italic>gene was induced at the transcriptional level early in Cr(VI) exposure, <italic>i.e</italic>. 90 min or less, but not at 24 h post-treatment [##REF##16524964##11##,##REF##16957260##12##]. Closer inspection of the 46 potentially co-regulated genes (Table ##TAB##0##1##) revealed that many of the transport/binding, hypothetical, and regulatory genes were shown previously to be upregulated in wild-type MR-1 cells in response to chromate treatment [##REF##16524964##11##,##REF##17385904##13##], thus suggesting that these ORFs are candidate gene targets of positive control by SO2426. The genes constituting Cluster A are discussed in greater detail in the sections below.</p> <p>Cluster B (Figure ##FIG##4##5##) includes 20 genes, three of which encode AcrB/AcrD/AcrF family proteins (SO1882, SO3484, SO4692) presumably involved in drug efflux and resistance. All of these 20 genes showed decreased expression in the Δ<italic>so2426 </italic>mutant except at the 24-h time point, where mRNA levels increased (Figure ##FIG##4##5## and Additional file ##SUPPL##0##1##). Notably, Cluster B also contains genes <italic>so0392 </italic>(putative lipoprotein) and <italic>so4688 </italic>(glycosyl transferase, group 2 family protein), both with annotated functions associated with the cell envelope, as well as <italic>cpxA </italic>(<italic>so4478</italic>). In <italic>E. coli</italic>, CpxA is the sensory inner membrane kinase that functions with its cognate partner, a cytoplasmic response regulator (CpxR), in a prototypical two-component system to sense and respond to perturbations in the bacterial cell envelope (reviewed in [##REF##15802241##35##]). <italic>S. oneidensis cpxR </italic>also displayed downregulated temporal expression profiles in the Δ<italic>so2426 </italic>mutant under chromate challenge but grouped with genes in Cluster A.</p> <title>Genes encoding transport and binding functions</title> <p>Thirteen of the genes in Cluster A have annotated functions in metal transport and binding, in particular Fe acquisition and homeostasis: <italic>ftn </italic>(ferritin), <italic>so1580 </italic>(TonB-dependent heme receptor), <italic>so1771 </italic>(GntP family permease), <italic>so2045 </italic>(cation efflux family protein), <italic>alcA </italic>(siderophore biosynthesis protein), <italic>so3031 </italic>and <italic>so3032 </italic>(both putative siderophore biosynthesis proteins encoded immediately downstream of <italic>alcA</italic>), <italic>so3033 </italic>(ferric alcaligin siderophore receptor), <italic>so3063 </italic>(sodium:alanine symporter family protein), <italic>so4150 </italic>(putative transporter), <italic>viuA </italic>(ferric vibriobactin receptor), <italic>so4712 </italic>(putative ABC transporter, ATP-binding protein), and <italic>so4743 </italic>(putative TonB-dependent receptor). In contrast to the TonB-dependent receptor genes <italic>so1580 </italic>and <italic>so4743</italic>, either no significant change or a maximal twofold decrease at the 180-min time point was observed for <italic>so2427 </italic>(putative TonB-dependent receptor), located directly upstream of the <italic>so2426 </italic>response regulator gene. Of the 13 genes with transport and binding functions in Cluster A, ten ORFs (<italic>ftn</italic>, <italic>so1580</italic>, <italic>so2045</italic>, <italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032</italic>, <italic>so3033</italic>, <italic>so3063</italic>, <italic>viuA</italic>, <italic>so4743</italic>) were shown in an earlier study to be induced in chromate-challenged wild-type MR-1 cells compared to untreated cells [##REF##16524964##11##,##REF##17385904##13##]. Temporal expression patterns for these genes demonstrated that they were downregulated 2- to 20-fold in the Δ<italic>so2426 </italic>strain over the 180-min time course (Table ##TAB##0##1##), suggesting that SO2426 acts as a direct or indirect positive regulator of a subset of Fe uptake and storage genes in <italic>S. oneidensis</italic>.</p> <p>The <italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032 </italic>operon encodes proteins required for siderophore biosynthesis in MR-1 [##REF##16391131##36##], and the product of gene <italic>so3033 </italic>is predicted to allow for the cellular utilization of the structurally undetermined MR-1 siderophore. The differential profiles for these four genes were characterized by a peak in down-regulated expression (ranging from ~7- to 17-fold) at the 180-min time interval post chromate addition, followed by no significant change in expression at 24 h (Table ##TAB##0##1##). The first gene in the siderophore biosynthetic operon shows 48% sequence identity to <italic>Bordetella pertussis alcA</italic>, which is required for alcaligin production [##REF##8759851##37##], while the two downstream genes (<italic>so3031 </italic>and <italic>so3032</italic>, respectively) diverge from <italic>alcB </italic>and <italic>alcC </italic>and are annotated as putative siderophore biosynthesis genes. While the structural identity of the MR-1 siderophore has not been elucidated, we predict that MR-1 likely produces a siderophore similar to putrebactin [##UREF##5##38##], a novel cyclic dihydroxamate siderophore characterized from <italic>Shewanella putrefaciens </italic>strain 200 and structurally similar to alcaligin.</p> <p>Consistent with the transcriptomic data, further physiological evidence for the involvement of the SO2426 response regulator in controlling siderophore-dependent iron acquisition was obtained by performing semiquantitative liquid CAS assays in which relative siderophore production levels in supernatants of the Δ<italic>so2426 </italic>mutant were compared to those in wild-type MR-1 cultures. Following 24-h growth in LB medium in the absence of added FeCl<sub>3 </sub>or chromate, the Δ<italic>so2426 </italic>mutant (<italic>A</italic><sub>630</sub>, 0.739 ± 0.02) exhibited essentially no detectable siderophore production over that of the cell-free control (<italic>A</italic><sub>630</sub>, 0.726 ± 0.005), whereas wild-type MR-1 cultures (<italic>A</italic><sub>630</sub>, 0.108 ± 0.06) produced approximately 7-fold more siderophore than the Δ<italic>so2426 </italic>mutant. As expected, addition of 50 μM FeCl<sub>3 </sub>to the culture medium reduced siderophore production by the wild-type to near background levels (an approximately 6-fold reduction) but did not affect siderophore accumulation in the Δ<italic>so2426 </italic>mutant. In the presence of 0.3 mM chromate, replicate cultures of wild-type MR-1 exhibited increasing levels in relative siderophore production over time (Figure ##FIG##5##6##), which was consistent with previous microarray expression data showing induction of siderophore biosynthesis genes in response to chromate exposure [##REF##16524964##11##]. By contrast, the temporal profile for the Δ<italic>so2426 </italic>mutant showed an initial moderate increase and then a dramatic reduction in siderophore excretion compared to the wild type (Figure ##FIG##5##6##). Collectively, these data support the hypothesis that SO2426 is a positive regulator of siderophore-mediated Fe uptake.</p> <p>Cluster A also includes an <italic>S. oneidensis </italic>gene that is involved in Fe storage, <italic>i.e</italic>., the <italic>ftn </italic>gene (<italic>so0139</italic>) encoding ferritin. The <italic>ftn </italic>mRNA expression levels were down-regulated 4- to 20-fold in the SO2426 null mutant across the entire temporal range, with the peak in repression occurring at the 180-min time point (Table ##TAB##0##1##). Although iron is an essential micronutrient for most organisms, free excess Fe in the presence of oxygen is potentially toxic because of its tendency to participate in Fenton chemistry reactions, which generate cell-damaging reactive oxygen intermediates. By scavenging and storing intracellular free Fe in a remobilizable form, Ftn protects against Fe(II)-mediated formation of hydroxyl radicals [##REF##9889981##39##,##REF##10620317##40##]. Studies have indicated that ferritins make important contributions to both Fe storage as well as cellular protection against oxidative stress [##REF##9782498##41##,##REF##8809765##42##] and that <italic>ftn </italic>transcription is co-ordinately regulated by iron and redox stress in some bacteria [##REF##15256555##43##]. <italic>S. oneidensis ftn </italic>was shown previously to be induced at a low level (~2-fold) in chromate-exposed wild-type MR-1 cells in contrast to untreated cells [##REF##16524964##11##]. The induction of <italic>ftn </italic>in <italic>S. oneidensis </italic>might act, at least partially, to mitigate chromate-generated oxidative stress by sequestering excess free Fe available for the catalysis of Fenton reactions and the formation of additional reactive hydroxyl radicals.</p> <title>Genes of unknown function</title> <p>The majority of coordinately repressed genes in Cluster A belong to the functional category of hypothetical and conserved hypothetical proteins (Table ##TAB##0##1##). Most notably, genes comprising the conserved hypothetical <italic>so1188</italic>-<italic>so1189</italic>-<italic>so1190 </italic>operon exhibited the highest degree of down-regulated expression in the chromate-exposed Δ<italic>so2426 </italic>mutant relative to the wild type, with decreases in mRNA ratios ranging from 45- to &gt; 300-fold over the 180-min period. No significant changes in expression, however, were observed for these genes at the 24-h time point. In addition, the hypothetical gene <italic>so2425</italic>, which is co-transcribed with <italic>so2426 </italic>in a polycistronic transcript (shown in this study), was repressed as much as 11-fold over the 180-min period, pointing to the possibility that the <italic>so2426 </italic>operon may be under positive autoregulatory control. Surprisingly, the functionally unknown genes <italic>so2424</italic>, <italic>so2423</italic>, and <italic>so2422</italic>, which also form part of the <italic>so2426 </italic>operon, showed either a less than 2-fold change or no significant differential expression over the time course.</p> <title>Regulatory genes</title> <p>In addition to metal transport and hypothetical genes described above, the temporal expression profiles for four transcriptional regulatory genes (<italic>so0916</italic>, <italic>so0544</italic>, <italic>cpxR</italic>, <italic>so4567</italic>) were altered in the Δ<italic>so2426 </italic>mutant. The function of the regulators encoded by <italic>so0916</italic>, <italic>so0544</italic>, and <italic>so4567 </italic>in <italic>S. oneidensis </italic>as well as their target genes are not known. However, gene <italic>so4477 </italic>encodes a homolog of CpxR, the response regulator component of the CpxAR envelope stress response system, which has been well characterized in <italic>E. coli </italic>(reviewed in [##REF##15802241##35##,##REF##16487683##44##]). The annotated <italic>cpxR </italic>gene in <italic>S. oneidensis </italic>MR-1 was downregulated as much as threefold in response to the <italic>so2426 </italic>deletion under chromate stress conditions (Table ##TAB##0##1##). This is in contrast to the two- to three-fold induction measured for <italic>cpxR </italic>in chromate-treated wild-type MR-1 cells [##REF##16524964##11##]. Located immediately downstream of <italic>cpxR </italic>in an apparent operon is <italic>cpxA</italic>, which was downregulated as much as 3-fold in the mutant but grouped in a separate cluster (Cluster B in Figure ##FIG##4##5##).</p> <p>In <italic>E. coli</italic>, the CpxAR two-component pair constitutes a stress response system the main function of which is to sense and respond to cell envelope distress, principally misfolded periplasmic proteins, by activating genes encoding proteases and folding catalysts [##REF##15802241##35##,##REF##7883164##45##,##REF##11544368##46##]. However, the complexity, extent, and response overlap of the Cpx system have become progressively apparent, as evident from its involvement in outer membrane porin expression [##REF##16077119##47##], stationary-phase survival [##REF##10542180##48##], and the response to high pH stress [##REF##9473036##49##] to name a few. Although the specific cellular role of the CpxR homolog in <italic>S. oneidensis </italic>has not been established, it is intriguing to consider the possibility that interplay between CpxR and SO2426 may occur in response to metal stress. It is conceivable that the CpxA-CpxR two-component system may respond to oxidative damage imposed on biomolecules as a consequence of reactive oxygen species generated during intracellular partial Cr(VI) reduction.</p> <title>Other differentially regulated genes in the <italic>so2426 </italic>deletion mutant</title> <p>While hypothetical and metal transport/binding genes are dominant in Cluster A (representing 65%), other notable genes with ≥ 2-fold-repressed mRNA levels in the Δ<italic>so2426 </italic>mutant include ORFs codifying a putative bicyclomycin resistance (Bcr) protein (SO2280) and a ferric iron reductase protein (SO3034). Bcr proteins are members of the major facilitator superfamily class of membrane efflux pumps and play important roles in drug resistance [##REF##11104814##50##]. In a previous work, the <italic>so2280 </italic>gene in a wild-type MR-1 genetic background exhibited 2- to 4-fold upregulated expression at the transcript level in response to short-term chromate exposure, <italic>i.e</italic>., up to 90 min [##REF##16524964##11##], but was shown in this study to be repressed as much as 4-fold in Cr(VI)-challenged Δ<italic>so2426 </italic>cells (Table ##TAB##0##1##), suggesting that the <italic>so2280 </italic>gene is a direct or indirect target of SO2426-mediated activation under chromate conditions. Although the cellular function of SO2280 is presently unknown, its co-ordinately regulated expression with the siderophore biosynthesis operon (<italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032</italic>) and gene <italic>so3033</italic>, encoding a ferric alcaligin siderophore receptor, is intriguing. Studies in other bacterial systems have shown that certain proteins resembling Bcr at the amino acid level function in siderophore export from the cell [##REF##15901687##51##].</p> <p>Downregulated expression of gene <italic>so3034 </italic>(a putative ferric iron reductase) peaked at the 90-min time point with a 4-fold decrease in mRNA levels (Table ##TAB##0##1##). Ferric iron reductases are thought to be involved in reducing siderophore-complexed Fe(III) to Fe(II), thus enabling dissociation of the high-affinity Fe<sup>3+</sup>-chelating siderophore from iron [##REF##12829269##30##]. The entire <italic>so3030</italic>-<italic>so3034 </italic>gene cluster, which is predicted to be involved in siderophore production and utilization, exhibits similar temporal expression profiles in the mutant, suggesting a role for the SO2426 response regulator in coordinating these functions. It is not clear whether the other differentially expressed energy metabolism genes (<italic>fdhG</italic>, <italic>fdnH</italic>, <italic>fdhE</italic>, <italic>so4503</italic>, <italic>so4506</italic>, <italic>so4509</italic>) contribute directly to chromate stress response or show altered transcription as a result of secondary effects associated with the <italic>so2426 </italic>deletion. These genes were not observed to be induced in response to chromate exposure in our previous transcriptomic studies using wild-type <italic>S. oneidensis </italic>MR-1 [##REF##16524964##11##].</p>
<title>Conclusion</title> <p><italic>S. oneidensis </italic>SO2426 is an apparent orphan response regulator with an N-terminal domain that contains an aspartate residue at the canonical receiver phosphorylation site and a predicted C-terminal wHTH motif indicative of response regulators in the OmpR subfamily. Although previous work in this laboratory demonstrated enhanced mRNA levels for the <italic>so2426 </italic>gene in a <italic>S. oneidensis </italic>Fur-deficient strain [##REF##15576789##31##] and in chromate-challenged wild-type MR-1 cells [##REF##16524964##11##], the basic molecular function of this predicted transcriptional regulator remained undefined. In this study, by integrating functional genomics and genetic approaches, we have shown that SO2426 is a direct or indirect positive regulator of siderophore-mediated Fe transport systems, the Fe storage gene <italic>ftn</italic>, and other as-yet uncharacterized cation transport genes. At this time, however, we cannot discount the possibility of negative regulation by SO2426 as well. Clustering analysis of temporal transcriptome data revealed a distinctive subset of 46 coordinately expressed genes that were consistently downregulated in the <italic>so2426 </italic>deletion mutant during the initial 180-min period of chromate exposure. Many of these genes – namely, <italic>ftn</italic>, the <italic>alcA</italic>-<italic>so3031</italic>-<italic>so3032 </italic>siderophore biosynthesis operon, <italic>so1580 </italic>(TonB-dependent heme receptor), <italic>so4743 </italic>(TonB-dependent receptor), <italic>viuA</italic>, <italic>so2280 </italic>(<italic>bcr</italic>), as well as a number of hypothetical genes – constituted the predominant molecular response to acute chromate stress in wild-type MR-1 [##REF##16524964##11##]. The transcriptomic analyses described here identified a number of possible directly SO2426-regulated genes, including those of unknown function, that may serve as targets for future studies.</p> <p>One possible hypothesis that can be generated from the collective data described here is that <italic>S. oneidensis </italic>MR-1 employs at least two control mechanisms for regulating Fe acquisition via siderophores under conditions of external Fe insufficiency: derepression as a result of decreased Fur activity and activation by SO2426 to enhance needed expression of Fe(III) transport systems. CAS siderophore assays clearly demonstrated that the Δ<italic>so2426 </italic>mutant was deficient in siderophore production compared to wild-type MR-1 cultures in the presence and absence of chromate, indicating that SO2426 plays a positive role in coordinating the regulation of siderophore-dependent Fe uptake mechanisms in <italic>S. oneidensis</italic>. Furthermore, the impaired chromate reduction rate of the <italic>so2426 </italic>deletion strain points to a connection between Fe transport and chromate stress response that requires further delineation at the molecular level. Our previous global studies characterizing the transcriptome and proteome of wild-type MR-1 in response to sub-lethal acute doses of chromate demonstrated a dramatic and unexpected induction of genes involved in Fe sequestration and uptake in contrast to unchallenged cells [##REF##16524964##11##,##REF##17385904##13##]. Two possible explanations for this observation are (i) Fe limitation under external chromate conditions, and (ii) direct or indirect displacement of the Fe(II) cofactor from Fur-binding sites under elevated Cr(VI), resulting in decreased DNA-binding activity of the Fur repressor and subsequent derepression of target genes. Exposure to other metals has been shown to perturb intracellular Fe pools and hence influence Fur-mediated repression [##REF##12950915##52##]. Future studies will focus on experimentally confirming direct gene targets of SO2426 control and further elucidating the linkage between regulation of Fe homeostasis and chromate stress response at molecular and physiological levels.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p><italic>Shewanella oneidensis </italic>MR-1 exhibits diverse metal ion-reducing capabilities and thus is of potential utility as a bioremediation agent. Knowledge of the molecular components and regulatory mechanisms dictating cellular responses to heavy metal stress, however, remains incomplete. In a previous work, the <italic>S. oneidensis so2426 </italic>gene, annotated as a DNA-binding response regulator, was demonstrated to be specifically responsive at both the transcript and protein levels to acute chromate [Cr(VI)] challenge. To delineate the cellular function of SO2426 and its contribution to metal stress response, we integrated genetic and physiological approaches with a genome-wide screen for target gene candidates comprising the SO2426 regulon.</p> <title>Results</title> <p>Inactivation of <italic>so2426 </italic>by an in-frame deletion resulted in enhanced chromate sensitivity and a reduced capacity to remove extracellular Cr(VI) relative to the parental strain. Time-resolved microarray analysis was used to compare transcriptomic profiles of wild-type and SO2426-deficient mutant <italic>S. oneidensis </italic>under conditions of chromate exposure. In total, 841 genes (18% of the arrayed genome) were up- or downregulated at least twofold in the Δ<italic>so2426 </italic>mutant for at least one of six time-point conditions. Hierarchical cluster analysis of temporal transcriptional profiles identified a distinct cluster (n = 46) comprised of co-ordinately regulated genes exhibiting significant downregulated expression (<italic>p </italic>&lt; 0.05) over time. Thirteen of these genes encoded proteins associated with transport and binding functions, particularly those involved in Fe transport and homeostasis (<italic>e.g</italic>., siderophore biosynthetic enzymes, TonB-dependent receptors, and the iron-storage protein ferritin). A conserved hypothetical operon (<italic>so1188</italic>-<italic>so1189</italic>-<italic>so1190</italic>), previously identified as a potential target of Fur-mediated repression, as well as a putative bicyclomycin resistance gene (<italic>so2280</italic>) and cation efflux family protein gene (<italic>so2045</italic>) also were repressed in the <italic>so2426 </italic>deletion mutant. Furthermore, the temporal expression profiles of four regulatory genes including a <italic>cpxR </italic>homolog were perturbed in the chromate-challenged mutant.</p> <title>Conclusion</title> <p>Our findings suggest a previously unrecognized functional role for the response regulator SO2426 in the activation of genes required for siderophore-mediated Fe acquisition, Fe storage, and other cation transport mechanisms. SO2426 regulatory function is involved at a fundamental molecular level in the linkage between Fe homeostasis and the cellular response to chromate-induced stress in <italic>S. oneidensis</italic>.</p>
<title>Authors' contributions</title> <p>KC created the deletion mutant analyzed in this work; performed the microarray experiments, qRT-PCR verifications, and the initial data analysis; and carried out the growth studies and Cr(VI) reduction assays. WW performed the operon structure analysis, transcription start site determination, and CAS siderophore assays. X–FW carried out statistical and clustering analyses of the microarray data. KC and X–FW both contributed to the manuscript writing. DKT conceived and coordinated the study and performed the majority of the manuscript writing. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Dr. Joel Klappenbach for plasmids and <italic>E. coli </italic>strain WM3064 used in site-directed mutagenesis and Dr. Steven Brown for assistance with the construction of the <italic>so2426 </italic>deletion mutant. We also thank Drs. Jizhong Zhou and Liyou Wu for <italic>S. oneidensis </italic>MR-1 microarrays, and Dr. Gene Wickham for helpful comments concerning the manuscript. This research was supported in part by the Office of Science (BER), U. S. Department of Energy, Grant No. DE-FG02-06ER64163, to DKT.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Sequence alignment of <italic>S. oneidensis </italic>SO2426 with other two-component response regulators</bold>. ClustalW [##REF##7984417##20##] was used to perform a multiple sequence alignment consisting of <italic>S. oneidensis </italic>SO2426 (GenBank ID <ext-link ext-link-type="gen" xlink:href="24348419">24348419</ext-link>), CpxR from <italic>E</italic>. <italic>coli </italic>(GenBank ID <ext-link ext-link-type="gen" xlink:href="16131752">16131752</ext-link>) and <italic>Vibrio cholerae </italic>(GenBank ID <ext-link ext-link-type="gen" xlink:href="147675245">147675245</ext-link>), and OmpR from <italic>E</italic>. <italic>coli </italic>(GenBank ID <ext-link ext-link-type="gen" xlink:href="16131282">16131282</ext-link>) and <italic>V</italic>. <italic>cholerae </italic>(GenBank ID <ext-link ext-link-type="gen" xlink:href="147673571">147673571</ext-link>). The region underlined with "=" is the aligned regulator receiver domain with predicted domain (SO2426: positions 13–124), and the region denoted with "~" is the aligned C-terminal domain containing the wHTH DNA-binding motif (SO2426: positions 158–235). Boldface letters highlighted in grey indicate conserved signature residues of receiver domains [##REF##10966457##19##]. Residue D62 is predicted as 4-aspartylphosphate, the putative phosphorylation site (highlighted in yellow). The star, colon, and dot notations rank the sequence conservation from high to low.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Structural organization of the <italic>S. oneidensis </italic>MR-1 <italic>so2426 </italic>locus</bold>. (A) Schematic representation of the <italic>so2426 </italic>gene region. ORFs located upstream and downstream of <italic>so2426 </italic>(black arrow) are indicated by open arrows, which also indicate the direction of transcription. The deduced proteins of the ORFs flanking <italic>so2426 </italic>have the following annotations based on the J. Craig Venter Institute's Comprehensive Microbial Resource <ext-link ext-link-type="uri" xlink:href="http://cmr.tigr.org"/>: SO2421 is a L-asparaginase I; SO2422 is a conserved hypothetical protein; SO2423 is a hypothetical protein; SO2424 is a zinc carboxypeptidase domain protein; SO2425 is a hypothetical protein; SO2426 is a DNA-binding response regulator; SO2427 is a putative TonB-dependent receptor; and SO2428 is a hypothetical protein. The locations of the oligonucleotide primers (P1–P9) used for RT-PCR are shown by the small solid arrows. The expected PCR product sizes for the different primer pairs are given below. (B) Agarose gel (1%) electrophoresis of amplified DNA fragments derived from <italic>S</italic>. <italic>oneidensis </italic>MR-1 cDNA templates under conditions of chromate stress or no stress. Lane designations: (1) 100-bp DNA ladder; (2) PCR primer pair P1/P4, (3) P1/P5, (4) P1/P6, (5) P1/P7, (6) P2/P8, (7) P9/P3, (8) blank, (9) 100-bp DNA ladder, (10) P1/P4, (11) P1/P5, (12) P1/P6, (13) P1/P7, (14) P2/P8, (15) P9/P3, (16) chromate-treated total RNA used as the template in PCR amplification with P1/P4 (negative control), and (17) non-stressed total RNA used as the template in PCR amplification with P1/P4 (negative control).</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Nucleotide sequence of the <italic>so2426 </italic>promoter region and N-terminal coding region</bold>. The transcriptional start site (+1) is indicated by arrows, and putative -10 and -35 basal promoter elements are presented in boldface and underlined. Closed circles indicate conserved residues matching a predicted Fur-binding recognition site. The annotated (original) translation start codon (M1 residue) of <italic>so2426 </italic>is shown in boldface. 5'-RACE results described here indicate that the annotated translation start codon is incorrect, and residue M11 constitutes the actual start codon. The 5'-RACE primers used to identify the transcription start site are shown in italics and underlined.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Growth and chromate reduction kinetics of <italic>S. oneidensis </italic>MR-1 wild type and the Δ<italic>so2426 </italic>mutant</bold>. (A) Wild-type (●) and mutant (○) cultures were cultivated in LB broth without the addition of chromate under aerobic conditions and compared to wild-type (▼) and mutant (△) cultures grown in LB broth amended with a final K<sub>2</sub>CrO<sub>4 </sub>concentration of 0.3 mM. (B) Cr(VI) removal rates for the <italic>S</italic>. <italic>oneidensis </italic>MR-1 (●), Δ<italic>so2426 </italic>(○), and complemented Δ<italic>so2426</italic>(<italic>so2426</italic><sup>+</sup>) (▼) strains in LB broth supplemented with 0.3 mM chromate. The different cultures were challenged with 0.3 mM K<sub>2</sub>CrO<sub>4 </sub>when cells reached an OD<sub>600 </sub>of 0.5 (mid-log point). The error bars indicate standard deviations of triplicate measurements.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Complete linkage clustering analysis of genes with altered expression profiles in the Δ<italic>so2426 </italic>mutant</bold>. The 841 genes (out of the 4,648 represented on the microarray) exhibiting at least a twofold statistically significant (<italic>p </italic>&lt; 0.05) change in expression for at least one of the time points during Cr(VI) exposure were analyzed by pairwise complete-linkage clustering. Transcriptional profiles are shown at 5, 30, 60, 90, 180 min, and 24 h post chromate addition. Individual genes are represented by a single row, and each exposure time interval is represented by a single column. Red represents induction, while green represents repression. Two noteworthy clusters (A and B) are indicated, with their respective expression heat maps enlarged.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Relative siderophore production by wild-type MR-1 and Δ<italic>so2426 </italic>mutant cells under chromate conditions</bold>. A final chromate concentration of 0.3 mM was added to mid-log-phase (OD<sub>600</sub>, 0.5) wild-type MR-1 (●) and Δ<italic>so2426 </italic>mutant (○) cultures cultivated in LB broth, and relative siderophore synthesis by each culture was monitored over time using the CAS assay. The relative siderophore production was calculated by subtracting the supernatant <italic>A</italic><sub>630 </sub>for the wild type or mutant from the uninoculated control and then determining the ratio of corrected supernatant <italic>A</italic><sub>630 </sub>to control <italic>A</italic><sub>630</sub>. The graphed values represent mean ratios ± standard errors (bars) for three replicate CAS measurements.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Selected genes with down-regulated expression profiles in the Δ<italic>so2426 </italic>mutant relative to wild-type MR-1</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="center" colspan="6"><bold>Mean (Δ<italic>so2426</italic>/WT) mRNA ratio<sup>†</sup></bold></td></tr><tr><td/><td/><td colspan="6"><hr/></td></tr><tr><td/><td/><td align="center" colspan="6"><bold>Time (min) post Cr(VI) addition</bold></td></tr><tr><td/><td/><td colspan="6"><hr/></td></tr><tr><td align="left"><bold>Gene</bold></td><td align="left"><bold>Description</bold></td><td align="left"><bold>5</bold></td><td align="left"><bold>30</bold></td><td align="left"><bold>60</bold></td><td align="left"><bold>90</bold></td><td align="left"><bold>180</bold></td><td align="left"><bold>1440</bold></td></tr></thead><tbody><tr><td align="left" colspan="8"><bold>Cellular processes</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SO2280</td><td align="left">Bicyclomycin resistance protein</td><td align="left">0.64</td><td align="left">0.36</td><td align="left">0.43</td><td align="left">0.23</td><td align="left">0.24</td><td align="left">0.39</td></tr><tr><td align="left">SO4274</td><td align="left">Undecaprenol kinase, putative</td><td align="left">0.67</td><td align="left">0.7*</td><td align="left">0.57</td><td align="left">0.56*</td><td align="left">0.8*</td><td align="left">0.36</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left" colspan="8"><bold>Energy metabolism</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SO0101</td><td align="left">Formate dehydrogenase (<italic>fdnG</italic>)</td><td align="left">0.3</td><td align="left">0.32</td><td align="left">0.55</td><td align="left">0.56</td><td align="left">0.29</td><td align="left">0.36</td></tr><tr><td align="left">SO0102</td><td align="left">Formate dehydrogenase (<italic>fdnH</italic>)</td><td align="left">0.37</td><td align="left">0.31</td><td align="left">0.72</td><td align="left">0.53</td><td align="left">0.25</td><td align="left">0.67*</td></tr><tr><td align="left">SO0104</td><td align="left">FdhE protein (<italic>fdhE</italic>)</td><td align="left">0.32</td><td align="left">0.74*</td><td align="left">0.45</td><td align="left">0.61*</td><td align="left">0.43</td><td align="left">0.76*</td></tr><tr><td align="left">SO0809</td><td align="left">Azurin precursor (<italic>azu</italic>)</td><td align="left">0.44</td><td align="left">0.69</td><td align="left">0.28</td><td align="left">0.84</td><td align="left">0.3</td><td align="left">0.4</td></tr><tr><td align="left">SO3034</td><td align="left">Ferric iron reductase protein, putative</td><td align="left">1.09*</td><td align="left">0.61</td><td align="left">0.91*</td><td align="left">0.25</td><td align="left">0.45</td><td align="left">0.42</td></tr><tr><td align="left">SO4151</td><td align="left">Polysaccharide deacetylase family protein</td><td align="left">0.29</td><td align="left">0.3</td><td align="left">0.35</td><td align="left">0.6</td><td align="left">0.24</td><td align="left">0.23</td></tr><tr><td align="left">SO4503</td><td align="left">Formate dehydrogenase accessory protein FdhD</td><td align="left">1.13*</td><td align="left">0.56</td><td align="left">0.53</td><td align="left">0.5</td><td align="left">0.28</td><td align="left">0.36</td></tr><tr><td align="left">SO4506</td><td align="left">Iron-sulfur cluster-binding protein</td><td align="left">0.42</td><td align="left">0.28</td><td align="left">0.77</td><td align="left">0.45</td><td align="left">0.47</td><td align="left">0.39</td></tr><tr><td align="left">SO4509</td><td align="left">Formate dehydrogenase, alpha subunit</td><td align="left">0.4</td><td align="left">0.29</td><td align="left">0.63</td><td align="left">0.38</td><td align="left">0.23</td><td align="left">0.37</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left" colspan="8"><bold>Hypothetical</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SOA0058</td><td align="left">Hypothetical protein</td><td align="left">1.19</td><td align="left">0.82*</td><td align="left">1.0*</td><td align="left">0.37*</td><td align="left">0.29</td><td align="left">0.27</td></tr><tr><td align="left">SO0496</td><td align="left">Conserved hypothetical protein</td><td align="left">0.59</td><td align="left">0.4</td><td align="left">0.48</td><td align="left">0.45</td><td align="left">0.46</td><td align="left">0.66*</td></tr><tr><td align="left">SO1188</td><td align="left">Conserved hypothetical protein</td><td align="left">0.01</td><td align="left">0.013</td><td align="left">0.007</td><td align="left">0.006</td><td align="left">0.004</td><td align="left">1.52*</td></tr><tr><td align="left">SO1189</td><td align="left">Conserved hypothetical protein</td><td align="left">0.003</td><td align="left">0.022</td><td align="left">0.006</td><td align="left">0.004</td><td align="left">0.005</td><td align="left">1.68*</td></tr><tr><td align="left">SO1190</td><td align="left">Conserved hypothetical protein</td><td align="left">0.01</td><td align="left">0.014</td><td align="left">0.02</td><td align="left">0.009</td><td align="left">0.02</td><td align="left">1.03*</td></tr><tr><td align="left">SO1770</td><td align="left">Glycerate kinase, putative</td><td align="left">0.56</td><td align="left">0.48</td><td align="left">0.35</td><td align="left">0.18</td><td align="left">0.22</td><td align="left">0.72*</td></tr><tr><td align="left">SO1967</td><td align="left">Hypothetical protein</td><td align="left">0.97*</td><td align="left">0.78*</td><td align="left">0.34</td><td align="left">0.55*</td><td align="left">0.72</td><td align="left">0.42</td></tr><tr><td align="left">SO2128</td><td align="left">Hypothetical protein</td><td align="left">0.38</td><td align="left">0.43</td><td align="left">0.32</td><td align="left">1.21</td><td align="left">0.54</td><td align="left">0.89*</td></tr><tr><td align="left">SO2425</td><td align="left">Hypothetical protein</td><td align="left">0.21</td><td align="left">0.29</td><td align="left">0.23</td><td align="left">0.09</td><td align="left">0.14</td><td align="left">1.63*</td></tr><tr><td align="left">SO2469</td><td align="left">Conserved hypothetical protein</td><td align="left">0.35</td><td align="left">0.53</td><td align="left">0.68</td><td align="left">0.62</td><td align="left">0.46</td><td align="left">0.96*</td></tr><tr><td align="left">SO3025</td><td align="left">Conserved hypothetical protein</td><td align="left">0.17</td><td align="left">0.16</td><td align="left">0.16</td><td align="left">0.2</td><td align="left">0.23</td><td align="left">0.73*</td></tr><tr><td align="left">SO3062</td><td align="left">Hypothetical protein</td><td align="left">0.48</td><td align="left">0.15</td><td align="left">0.16</td><td align="left">0.14</td><td align="left">0.21</td><td align="left">1.21*</td></tr><tr><td align="left">SO4502</td><td align="left">Conserved domain protein</td><td align="left">0.93*</td><td align="left">0.74*</td><td align="left">0.54</td><td align="left">0.54*</td><td align="left">0.26</td><td align="left">0.3</td></tr><tr><td align="left">SO4504</td><td align="left">Conserved hypothetical protein</td><td align="left">0.2</td><td align="left">0.06</td><td align="left">0.15</td><td align="left">0.14</td><td align="left">0.23</td><td align="left">0.13</td></tr><tr><td align="left">SO4505</td><td align="left">Conserved hypothetical protein</td><td align="left">0.23</td><td align="left">0.16</td><td align="left">0.45</td><td align="left">0.23</td><td align="left">0.29</td><td align="left">0.16</td></tr><tr><td align="left">SO4689</td><td align="left">Conserved hypothetical protein</td><td align="left">0.51</td><td align="left">0.43</td><td align="left">0.4</td><td align="left">0.45</td><td align="left">0.54</td><td align="left">0.76*</td></tr><tr><td align="left">SO4719</td><td align="left">Conserved hypothetical protein</td><td align="left">0.23</td><td align="left">0.36</td><td align="left">0.36</td><td align="left">0.31</td><td align="left">0.18</td><td align="left">0.62*</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left" colspan="8"><bold>Protein synthesis</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SO0106</td><td align="left">Selenocysteine-specific translation elongation factor (<italic>selB</italic>)</td><td align="left">0.38</td><td align="left">0.72*</td><td align="left">0.57</td><td align="left">0.76*</td><td align="left">0.46</td><td align="left">0.83*</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left" colspan="8"><bold>Regulatory functions</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SO0916</td><td align="left">Transcriptional regulator, MarR family</td><td align="left">0.64</td><td align="left">0.42</td><td align="left">0.4</td><td align="left">0.34</td><td align="left">0.33</td><td align="left">0.84*</td></tr><tr><td align="left">SO0544</td><td align="left">Sensory box histidine kinase</td><td align="left">0.53</td><td align="left">0.31</td><td align="left">0.35</td><td align="left">0.51</td><td align="left">0.52</td><td align="left">0.58*</td></tr><tr><td align="left">SO4477</td><td align="left">Transcriptional regulatory protein CpxR (<italic>cpxR</italic>)</td><td align="left">0.54</td><td align="left">0.32</td><td align="left">0.37</td><td align="left">0.35</td><td align="left">0.43</td><td align="left">0.89*</td></tr><tr><td align="left">SO4567</td><td align="left">Transcriptional regulator, AsnC family</td><td align="left">0.41</td><td align="left">0.39</td><td align="left">0.44</td><td align="left">0.79</td><td align="left">0.34</td><td align="left">0.45</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left" colspan="8"><bold>Transport and binding proteins</bold></td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="left">SO0139</td><td align="left">Ferritin (<italic>ftn</italic>)</td><td align="left">0.2</td><td align="left">0.22</td><td align="left">0.11</td><td align="left">0.08</td><td align="left">0.05</td><td align="left">0.19</td></tr><tr><td align="left">SO1580</td><td align="left">TonB-dependent heme receptor</td><td align="left">0.32</td><td align="left">0.39</td><td align="left">0.37</td><td align="left">0.19</td><td align="left">0.23</td><td align="left">0.92*</td></tr><tr><td align="left">SO1771</td><td align="left">Permease, GntP family</td><td align="left">0.88*</td><td align="left">0.23</td><td align="left">0.43</td><td align="left">0.37</td><td align="left">0.14</td><td align="left">0.2</td></tr><tr><td align="left">SO2045</td><td align="left">Cation efflux family protein</td><td align="left">0.57</td><td align="left">0.24</td><td align="left">0.3</td><td align="left">0.22</td><td align="left">0.13</td><td align="left">1.19*</td></tr><tr><td align="left">SO3030</td><td align="left">Siderophore biosynthesis protein AlcA (<italic>alcA</italic>)</td><td align="left">0.67</td><td align="left">0.27</td><td align="left">0.81*</td><td align="left">0.21</td><td align="left">0.15</td><td align="left">0.69*</td></tr><tr><td align="left">SO3031</td><td align="left">Siderophore biosynthesis protein, putative</td><td align="left">0.52</td><td align="left">0.34</td><td align="left">0.3</td><td align="left">0.16</td><td align="left">0.08</td><td align="left">1.37*</td></tr><tr><td align="left">SO3032</td><td align="left">Siderophore biosynthesis protein, putative</td><td align="left">0.5</td><td align="left">0.33</td><td align="left">0.25</td><td align="left">0.14</td><td align="left">0.06</td><td align="left">0.73*</td></tr><tr><td align="left">SO3033</td><td align="left">Ferric alcaligin siderophore receptor</td><td align="left">0.73</td><td align="left">0.34</td><td align="left">0.29</td><td align="left">0.17</td><td align="left">0.07</td><td align="left">0.77*</td></tr><tr><td align="left">SO3063</td><td align="left">Sodium:alanine symporter family protein</td><td align="left">0.26</td><td align="left">0.18</td><td align="left">0.12</td><td align="left">0.15</td><td align="left">0.18</td><td align="left">0.34</td></tr><tr><td align="left">SO4150</td><td align="left">Transporter, putative</td><td align="left">0.07</td><td align="left">0.12</td><td align="left">0.11</td><td align="left">0.34</td><td align="left">0.27</td><td align="left">0.3</td></tr><tr><td align="left">SO4516</td><td align="left">Ferric vibriobactin receptor (<italic>viuA</italic>)</td><td align="left">0.59</td><td align="left">0.54</td><td align="left">0.37</td><td align="left">0.21</td><td align="left">0.16</td><td align="left">0.67*</td></tr><tr><td align="left">SO4712</td><td align="left">ABC transporter, ATP-binding protein, putative</td><td align="left">0.49</td><td align="left">0.84*</td><td align="left">0.59</td><td align="left">0.77*</td><td align="left">0.54</td><td align="left">0.59*</td></tr><tr><td align="left">SO4743</td><td align="left">TonB-dependent receptor, putative</td><td align="left">0.08</td><td align="left">0.16</td><td align="left">0.09</td><td align="left">0.06</td><td align="left">0.05</td><td align="left">0.43</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Bacterial strains and plasmids used in this study</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Strain or plasmid</bold></td><td align="left"><bold>Relevant characteristics</bold></td><td align="left"><bold>Source or reference</bold></td></tr></thead><tbody><tr><td align="left" colspan="3"><bold>Strains</bold></td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><italic>S. oneidensis</italic></td><td/><td/></tr><tr><td align="left"> MR-1 (ATCC 700550)</td><td align="left">Wild type</td><td align="left">Laboratory Stock</td></tr><tr><td align="left"> MR-1/Δ<italic>so2426</italic></td><td align="left">In-frame deletion of the <italic>so2426 </italic>locus</td><td align="left">This study</td></tr><tr><td align="left"> MR-1/Δ<italic>so2426</italic>(<italic>so2426</italic><sup>+</sup>)</td><td align="left">Δ<italic>so2426 </italic>mutant complemented with pBBRso2426</td><td align="left">This study</td></tr><tr><td align="left"><italic>E. coli </italic>WM3064</td><td align="left">DAP-auxotrophic wild type</td><td align="left">[##REF##16391095##32##]</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left" colspan="3"><bold>Plasmids</bold></td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"> pJK100</td><td align="left">Allelic-exchange vector</td><td align="left">[##REF##16391095##32##]</td></tr><tr><td align="left"> pJK102</td><td align="left">pJK100 with <italic>so2426 </italic>upstream and downstream regions</td><td align="left">This study</td></tr><tr><td align="left"> pCM157</td><td align="left">Cre recombinase expression vector</td><td align="left">[##REF##16391095##32##]</td></tr><tr><td align="left"> pBBR1MCS-5</td><td align="left">Gm<sup>r</sup>, broad-host range vector, <italic>lacPOZ'</italic></td><td align="left">[##REF##8529885##56##]</td></tr><tr><td align="left"> pBBRso2426</td><td align="left">pBBR1MCS5-based construct harboring entire <italic>so2426 </italic>gene</td><td align="left">This study</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Oligonucleotide primers used in this study</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Primer</bold></td><td align="left"><bold>DNA sequence (5' → 3')</bold></td><td align="left"><bold>Reference or source</bold></td></tr></thead><tbody><tr><td align="left" colspan="3"><bold>Construction of Δ<italic>so2426 </italic>strain</bold></td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">Del840-F</td><td align="left">CTTGGTTACCGGCTAGTGAAC</td><td align="left">This study</td></tr><tr><td align="left">Del840-R</td><td align="left">GGCAGGTATTGATAACAATGA</td><td align="left">This study</td></tr><tr><td align="left">Del900-F</td><td align="left">GGTTCACACCAATCGCATTAG</td><td align="left">This study</td></tr><tr><td align="left">Del900-R</td><td align="left">TGGCCAATACCCGCTTACCGC</td><td align="left">This study</td></tr><tr><td align="left">ISU2426</td><td align="left">GCCTAAGATGCCATCAGT</td><td align="left">This study</td></tr><tr><td align="left">ISD2426</td><td align="left">TCTTCAAGATTTAGCTTATCC</td><td align="left">This study</td></tr><tr><td align="left">OP840</td><td align="left">CACATAAGGCAGACCTTCGTC</td><td align="left">This study</td></tr><tr><td align="left">OP900</td><td align="left">ATGGTCCGTACTGTGGCCGC</td><td align="left">This study</td></tr><tr><td align="left">kanF</td><td align="left">ATTGTTGATGCGCTGGCAGT</td><td align="left">[##REF##16391095##32##]</td></tr><tr><td align="left">kanR</td><td align="left">TCCGGTGAGAATGGCAAAAG</td><td align="left">[##REF##16391095##32##]</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left" colspan="3"><bold>Construction of complementation plasmid</bold></td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">2426com-F</td><td align="left">ACACACAAGCTTGCGCTTTTCTTTTAGGTACAA</td><td align="left">This study</td></tr><tr><td align="left">2426com-R</td><td align="left">ACACACGGATCCGACTCACAGAGGGCGCTTA</td><td align="left">This study</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left" colspan="3"><bold>RT-PCR analysis</bold></td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">P1</td><td align="left">ATGATATTGATTCTCGTTTG</td><td align="left">This study</td></tr><tr><td align="left">P2</td><td align="left">CCGAGTGTGGTTAACTGATG</td><td align="left">This study</td></tr><tr><td align="left">P3</td><td align="left">CGCCGAGTATTACTGATATGC</td><td align="left">This study</td></tr><tr><td align="left">P4</td><td align="left">AAGCGCTAAAACTGTATCC</td><td align="left">This study</td></tr><tr><td align="left">P5</td><td align="left">TTAACATGCATCTACTTTTA</td><td align="left">This study</td></tr><tr><td align="left">P6</td><td align="left">GATCTTGCAGGTTGTTGTT</td><td align="left">This study</td></tr><tr><td align="left">P7</td><td align="left">TCATACACTCTTTCGCTTAT</td><td align="left">This study</td></tr><tr><td align="left">P8</td><td align="left">CTGTTTCTTCAACTCAGCCT</td><td align="left">This study</td></tr><tr><td align="left">P9</td><td align="left">CAGTCGTTAGCTCAATTGCT</td><td align="left">This study</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Complete processed microarray dataset for the <italic>S. oneidensis </italic>Δ<italic>so2426 </italic>mutant compared to the wild type under chromate conditions</bold>. This file provides the complete statistically analyzed microarray output for all six time points (5, 30, 60, 90, 180 min, and 24 h) during chromate (0.3 mM) exposure starting with the Gene ID, mean Δ<italic>so2426 </italic>mutant/WT expression ratio, <italic>n </italic>(number of signal values out of 12 total per gene included in the statistical analysis), statistical significance, gene name, gene product, main role category, and subrole category. ArrayStat™ (Imaging Research, Inc., Ontario, Canada) was used to determine the common error of normalized expression values, remove outliers, and determine statistical significance via a <italic>z </italic>test for two independent conditions and the false discovery rate method (nominal α, confidence cutoff of <italic>p </italic>&lt; 0.05).</p></caption></supplementary-material>
<table-wrap-foot><p><sup>†</sup>The mRNA ratios (Δ<italic>so2426 </italic>mutant/wild type) per time interval represent mean values derived from six DNA microarray experiments performed with total RNA isolated from three independent cultures (plus two dye reversal reactions per biological replicate) in chromate-amended LB medium.</p><p>*Expression ratio values determined to be statistically non-significant (<italic>p </italic>&gt; 0.05).</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-395-1"/> <graphic xlink:href="1471-2164-9-395-2"/> <graphic xlink:href="1471-2164-9-395-3"/> <graphic xlink:href="1471-2164-9-395-4"/> <graphic xlink:href="1471-2164-9-395-5"/> <graphic xlink:href="1471-2164-9-395-6"/>
<media xlink:href="1471-2164-9-395-S1.zip" mimetype="application" mime-subtype="x-zip-compressed"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Lovley", "Phillips", "Gorby", "Landa"], "given-names": ["DR", "EJP", "YA", "ER"], "article-title": ["Microbial reduction of uranium"], "source": ["Nature"], "year": ["1991"], "volume": ["350"], "fpage": ["413"], "lpage": ["416"]}, {"surname": ["James"], "given-names": ["BR"], "article-title": ["The challenge of remediating chromium-contaminated soil"], "source": ["Environ Sci Technol"], "year": ["1996"], "volume": ["30"], "fpage": ["A248"], "lpage": ["A251"]}, {"surname": ["Langard", "Waldron HA"], "given-names": ["S"], "article-title": ["Chromium"], "source": ["Metals in the Environment"], "year": ["1980"], "publisher-name": ["New York: Academy Press Inc"], "fpage": ["111"], "lpage": ["132"]}, {"surname": ["Hoch", "Silhavy"], "given-names": ["JA", "TJ"], "source": ["Two-Component Signal Transduction"], "year": ["1995"], "publisher-name": ["Washington, DC: American Society for Microbiology Press"]}, {"surname": ["Wan", "Xu"], "given-names": ["XF", "D"], "article-title": ["Intrinsic terminator prediction and its application in "], "italic": ["Synechococcus "], "source": ["J Comp Sci Tech"], "year": ["2005"], "volume": ["20"], "fpage": ["465"], "lpage": ["482"]}, {"surname": ["Ledyard", "Butler"], "given-names": ["KM", "A"], "article-title": ["Structure of putrebactin, a new dihydroamate siderophore produced by "], "italic": ["Shewanella putrefaciens"], "source": ["J Biol Inorg Chem"], "year": ["1997"], "volume": ["2"], "fpage": ["93"], "lpage": ["97"]}, {"article-title": ["Gene Expression Omnibus"]}]
BMC Genomics. 2008 Aug 21; 9:395