- Research article
- Open Access
Fnr (EtrA) acts as a fine-tuning regulator of anaerobic metabolism in Shewanella oneidensisMR-1
BMC Microbiologyvolume 11, Article number: 64 (2011)
EtrA in Shewanella oneidensis MR-1, a model organism for study of adaptation to varied redox niches, shares 73.6% and 50.8% amino acid sequence identity with the oxygen-sensing regulators Fnr in E. coli and Anr in Pseudomonas aeruginosa, respectively; however, its regulatory role of anaerobic metabolism in Shewanella spp. is complex and not well understood.
The expression of the nap genes, nrfA, cymA and hcp was significantly reduced in etrA deletion mutant EtrA7-1; however, limited anaerobic growth and nitrate reduction occurred, suggesting that multiple regulators control nitrate reduction in this strain. Dimethyl sulfoxide (DMSO) and fumarate reductase gene expression was down-regulated at least 2-fold in the mutant, which, showed lower or no reduction of these electron acceptors when compared to the wild type, suggesting both respiratory pathways are under EtrA control. Transcript analysis further suggested a role of EtrA in prophage activation and down-regulation of genes implicated in aerobic metabolism.
In contrast to previous studies that attributed a minor regulatory role to EtrA in Shewanella spp., this study demonstrates that EtrA acts as a global transcriptional regulator and, in conjunction with other regulators, fine-tunes the expression of genes involved in anaerobic metabolism in S. oneidensis strain MR-1. Transcriptomic and sequence analyses of the genes differentially expressed showed that those mostly affected by the mutation belonged to the "Energy metabolism" category, while stress-related genes were indirectly regulated in the mutant possibly as a result of a secondary perturbation (e.g. oxidative stress, starvation). We also conclude based on sequence, physiological and expression analyses that this regulator is more appropriately termed Fnr and recommend this descriptor be used in future publications.
Due to its respiratory versatility, Shewanella oneidensis strain MR-1 serves as a model organism for studying the regulation of aerobic and anaerobic growth [1–3]. In contrast to Escherichia coli, the regulatory systems that control transcription of genes responsible for different respiratory processes are poorly understood in environmentally relevant Shewanella spp. [4–7]. In E. coli, the transition from aerobic to anaerobic metabolism is primarily regulated by Fnr (fumarate and nitrate reduction regulator) and by the two-component regulatory system ArcAB (aerobic respiration control) [8–11]. A gene expression study in E. coli K12 indicated that one-third of its 4,290 genes were differentially expressed during aerobic versus anaerobic growth . Among the differentially expressed genes, 712 (49%) genes were directly or indirectly affected by Fnr. Fnr possesses a [4Fe-4S]2+ cluster that acts as an oxygen sensory domain . Fnr in its active dimeric form binds to target DNA sequences inducing or repressing transcription [14, 15]. Under aerobic conditions, or when oxygen levels increase, an Fe2+ atom in the [4Fe-4S]2+ cluster is oxidized resulting in the formation of a [2Fe-2S]2+ cluster via a [3Fe-4S]1+ intermediate. This oxidation causes a conformation change in Fnr, thus altering its affinity to DNA and regulatory control of transcription [14, 15].
Studies with strain MR-1 mutants have identified three important regulators of anaerobic metabolism: EtrA (electron transport regulator protein), ArcA and CRP (cyclic AMP receptor protein) [4, 5, 16–19]. Sequence alignment of the protein encoded by etrA reveal that the four cysteine residues that form the [4Fe-4S]2+ cluster in Fnr are conserved in EtrA . In a gene replacement study, etrA of strain MR-1 restored wild type physiology of an E. coli fnr deletion mutant . EtrA shares 73.6% and 50.8% of amino acid sequence identity with Fnr in E. coli and Anr (arginine deaminase and nitrate reductase anaerobic regulator) in Pseudomonas aeruginosa, respectively. This high degree of similarity suggests that EtrA has a regulatory function in MR-1, possibly by sensing oxygen. Despite the lack of physiological evidence to support a regulatory role of EtrA in the anaerobic metabolism of strain MR-1 , a gene expression study using a partial microarray (691 ORFs) of the strain MR-1's genome suggested involvement of EtrA in the regulation of the transcription of genes associated with aerobic and anaerobic metabolism . Growth experiments with an etrA deletion mutant in S. oneidensis strain DSP10 (a spontaneous rifampicin resistant mutant of MR-1) implicated EtrA in the regulation of genes related to aerobic and anaerobic metabolism, similar to what has been observed for Fnr in E. coli [12, 20]. Unfortunately, the implications of these findings cannot be interpreted unambiguously since the rifampicin resistance of strain DSP10 influences electron transport .
To examine the regulatory role of EtrA in strain MR-1 in more detail, we generated an etrA knockout mutant EtrA7-1 in a wild type background. Growth and phenotypic characterization of this mutant combined with a whole genome transcriptome analysis confirms that EtrA regulates nitrate and fumarate reduction, plus provides experimental evidence for its positive regulatory role in DMSO reduction. Our genome-wide expression analysis shows differential expression of 612 genes for which sequence analysis recognized a EtrA motif for 72 of the operons encoding 118 genes, suggesting that their regulation is via direct interaction of EtrA with its promoters. Most of these genes are associated with metabolic functions.
Genotypic and phenotypic characterization of a ΔetrA::loxPmutant
The growth of the etrA knockout mutant EtrA7-1 with nitrate was significantly impaired as cultures reached a maximum OD600 of 0.02, at least 5-fold lower than the wild type strain (Figure 1). In addition, the doubling time for the mutant under these conditions was approximately 10 h compared to a doubling time of 2 h for the wild type. Plasmid pCCG03 carrying etrA, but not the parental pCM62 vector lacking etrA, restored near wild type growth to the EtrA7-1 mutant, which confirms that the observed phenotype was attributable to the deletion of etrA. After 10 h of incubation, nitrate was reduced in wild type and complemented EtrA7-1 cultures though less nitrate was reduced in the latter consistent with its slightly slower growth (Figure 2). After 23 h, both the wild type and EtrA7-1 complement had completely converted nitrate and nitrite to ammonia, while nitrite was the major product in cultures of the EtrA7-1 mutant and the control strain harboring pCM62 (Figure 2).
Anaerobic cultures of the mutant and the wild type strain were analyzed for the reduction of different electron acceptors with lactate as the electron donor. No growth of the EtrA7-1 mutant was observed with fumarate as electron acceptor whereas the wild type strain reached an OD600 of 0.053 ± 0.005. Limited growth (approximately 50% lower OD600 compared with the wild type cultures) was observed in mutant cultures amended with trimethylamine N-oxide (TMAO) or thiosulfate (data not shown). No OD increases with the mutant and the wild type were measured with DMSO provided as electron acceptor at 2 and 10 mM; however, HPLC analyses of cultures with 2 mM DMSO revealed that DMSO was completely consumed in wild type cultures, whereas no DMSO consumption was evident in the mutant cultures (Figure 3). No changes in DMSO concentrations were observed in cultures with 10 mM DMSO. No significant differences in Fe(III), Mn(IV) and sulfite reduction rates were observed between the wild type and the EtrA7-1 deletion mutant (Figure 3). Anaerobic cultures of the mutant and the wild type strains grown with pyruvate instead of lactate as electron donor showed similar results, i.e., the mutant showed limited or no growth with nitrate, fumarate and DMSO provided as electron acceptors compared to the wild type (Figure 4). Similar to the lactate-amended cultures, the rates of nitrate, fumarate and DMSO reduction in wild type cultures exceeded those measured in cultures of the mutant strain (Table 1). Resting cell assays corroborated these findings and nitrate reduction and ammonium production occurred at higher rates in assays with wild type cells. Complete stoichiometric conversion to ammonium also occurred in the assays with mutant cells, although lower rates and a 3-fold longer incubation were required for complete reduction (i.e., 24 h for the EtrA7-1 versus 8 h for the wild type) (Figure 5).
Effects of etrAdeletion on transcription
The global transcriptome profile of mutant strain EtrA7-1 grown anaerobically with nitrate as the sole electron acceptor was compared to that of the wild type under the same growth conditions. A complete list of all the genes differentially expressed two-fold or higher is provided as supplemental information (Additional file 1). Out of 612 differentially transcribed genes in the EtrA7-1 mutant relative to the wild type, 289 were up-regulated and 323 were down-regulated. The differentially transcribed genes were classified in 19 functional "TIGR Role" categories (Additional file 2) based on the MR-1 genome annotation (GenBank accession number AE014299) . Genes with unknown functions represented the largest category of up-regulated (14.8%) and the second most common category of down-regulated genes (17.3%). Genes associated with energy metabolism were the largest category (17.6%) of down-regulated genes (Additional file 2). Among the up-regulated genes, the "Protein synthesis" category ranked second (12.5%) and the "Other categories" ranked third (11.4%). This latter category included phage-, transposase- and plasmid-related genes. The "Energy metabolism" category represented 9.7%, ranking fourth.
Identification of putative EtrA binding sites
The promoters of the differentially expressed genes were examined for putative EtrA binding sites in order to identify those genes that were likely directly regulated by EtrA from the many genes whose expression changes were most likely due to secondary effects. For example, the up-regulation of phage-related genes is likely a response to stress, and not a direct result of the etrA deletion. Putative EtrA binding sites were identified for those genes that showed at least 2.5 fold change in the EtrA7-1 mutant relative to the wild type by extracting the upstream intergenic sequences and applying the Gibbs centroid sampler to identify a putative binding site motif. A motif was identified (Additional file 3) that displays similarities to the E. coli Fnr and Crp binding sites motifs; this motif was present upstream of 44 operons that encode a total of 78 genes. The largest proportion of these genes is in the "Energy metabolism" category (Table 2 and 3, Additional file 2). Binding sites were detected upstream of an additional 28 operons when the detected motif (Additional file 3) was used to scan the upstream intergenic regions of all genes listed in Additional file 1.
Regulatory role of EtrA in energy metabolism
Since the "Energy metabolism" category contained the largest group of genes responsive to EtrA, these genes were analyzed in more detail. Up-regulated genes (Table 2) in this group included genes encoding a cytochrome c oxidase (ccoPQN [SO2361-2362, SO2364]), proteins involved in gluconeogenesis such as PckA (SO0162), and nqrABCDEF-2 genes (SO1103-1108) encoding NADH:ubiquinone oxidoreductases. From this group, only the nqr gene clusters had a putative EtrA binding site.
While the nqr-2 gene cluster was up-regulated in the etrA knockout mutant, the nqr-1 gene cluster (SO0903-0907) was down-regulated. Nqr is a Na+ pump that during respiration generates a sodium motive force to mediate solute transport, flagellar motility and ATP synthesis . Both nqr gene clusters had putative EtrA binding sites. The microarray data indicated that EtrA affects the transcription pattern of these genes differently. Similarly, the etrA deletion had a distinct effect on the expression of the fdh gene clusters encoding a formate dehydrogenase. The fdh-1 genes (SO4508-4511) were up-regulated whereas the fdh-2 gene cluster (SO4512-4515) was down-regulated. An EtrA binding site was only identified for the fdh-2 cluster and not for the fdh-1 cluster, indicating EtrA affects both clusters differently.
Other up-regulated genes in the "Energy metabolism" category included the succinate dehydrogenase gene sdhC (SO1927), the succinyl-CoA synthase operon sucABCD (SO1930-1933), the butyryl-CoA:acetate CoA-transferase and the acetyl CoA-synthase genes (SO1891-1892).
Down-regulated genes in the "Energy metabolism" category (Table 3), which were predicted for direct regulation by EtrA, included genes involved in anaerobic metabolism, especially in nitrate respiration like the napBHGAD operon (SO0845-0849, strongly repressed with all subunit genes except napD with ratios ≤ 0.18), and the nrfA (SO3980) genes. cymA (SO4591; ratio 0.39), the prismane protein hcp gene (SO1363), and neighboring protein hcr gene (SO1364), both of which were strongly repressed (ratios ≤ 0.13) and have been associated with the nitrate reduction pathway [24–27], did not show evidence of EtrA binding sites. Also indirectly down-regulated were the fumarate reductase genes frdAB (SO0398-0399) and fccA (SO0970), the ackA and the pta (SO2915-16) genes involved in acetate production and the ppc (SO0274) gene encoding an acetate phosphoenol pyruvate carboxylase. The hyaCBA (SO2097-2099) genes encoding a quinone-reactive Ni/Fe hydrogenase were highly indirectly repressed (ratio ≤ 0.11). Among the genes identified as directly down-regulated are all the genes in the operon that encodes the anaerobic DMSO reductase (dmsAB) (SO1428-32), the cydAB genes (SO3285-3286) encoding a cytochrome d oxidase complex, as well as genes involved in metabolism of organic compounds such as the pflAB (SO2912-2913).
Other down-regulated genes grouped in different categories included genes encoding ABC transporters (cydCD [SO3779-3780], SO4446-4448), TonB-dependent receptors (nosA [SO0630]), and L-lactate permease (lldP [SO0827]) and a putative lactate permease (SO1522). The only gene directly down-regulated from this later group is lldP (SO0827), for which an EtrA binding site was predicted (Table 3). As expected, the cDNA for etrA, shows no significant hybridization signal in EtrA7-1 mutant (ratio 0.05).
Stress response caused by the etrAdeletion
We detected induction of genes from various categories, which have been associated with stress response i.e., starvation, phage infection and oxidative stress, possibly due to accumulation of nitrogen oxide reactive species. Up-regulated genes (Additional file 1) were dominated by genes grouped in "Other categories". The majority of up-regulated genes were phage-related. For example, 25 genes of the LambdaSo phage (SO2940-2974), a gene encoding a viral capsid protein of the MuSo1 phage (SO0675), and genes of MuSo2 phage (SO2684-2685, SO2687, SO2702) were up-regulated. In contrast, the gene encoding the LambdaSo phage transcriptional regulator of the Cro/CI family (SO2990) was down-regulated (ratio 0.43). Transcriptional changes of most of these genes are likely indirect effects due to the deletion of the etrA gene and only for the LambdaSo phage genes S02957-2962 was an EtrA binding site predicted.
The category "Transport and binding proteins" contains a large number of genes associated with stress response. Up-regulated genes in this group contain genes encoding heavy metal efflux pumps (SO0520, SO4596-4598, SOA0153) and genes for phosphate transport (SO1560, SO1723-1724 including pstB-1 [SO1725], pstB-2 [SO4289] and pstA [SO4290]). There was up-regulation of phoBR (SO1558-59) and phoU (SO1726) genes, which regulate the phosphate transporters genes during phosphate starvation [28–32].
Up-regulated genes in response to stress conditions i.e., starvation, phage infection, oxidative stress, include a stringent starvation protein encoded by the sspAB genes (SO0611-0612), and a phage shock protein operon pspABC (SO1807-1809). Other up-regulated stress-related genes were the RNA polymerase sigma-70 factor rpoD (SO1284)[32, 35], a GTP-binding protein that regulates the TCA cycle and responds to starvation (era [SO1349]), and a DNA repair protein (recO [SO1350]).
The results of this study demonstrate that EtrA positively regulates dissimilatory nitrate, fumarate and DMSO reduction pathways in S. oneidensis MR-1. The generation of etrA knockout mutant EtrA7-1 in the wild type strain MR-1 background eliminated any possible secondary effects on the phenotype, such as the electron transfer perturbation suspected with the rifampicin resistant DSP10 strain . Similar to other etrA mutants of strain MR-1, EtrA7-1 retained its ability to reduce nitrate [6, 7, 16]; however, our results show that the anaerobic growth of the mutant was significantly impaired compared to the wild type when nitrate was the only electron acceptor. Likewise, the etrA deletion mutant lost its ability to reduce fumarate and DMSO with both lactate and pyruvate as electron donor. Regulation of DMSO reduction by EtrA in strain MR-1 was suggested previously  however this study provides physiological evidence that confirm its role. The ability of the EtrA7-1 mutant to reduce TMAO and thiosulfate also decreased; however the reduction of Fe(III) citrate, Mn(IV) and sulfite was not affected by the deletion. No differences in growth performance between the wild type and the mutant were observed under aerobic conditions (data not shown).
The transcriptome analysis provides a genome-wide expression profile of S. oneidensis MR-1 instead of the partial genome array that was previously evaluated (691 ORFs  vs 4,648 genes in this study). We observed in 612 (13%) differentially expressed genes represented though some are likely due to differences in growth rate between the mutant strain and the wild type strain. Nonetheless, the expression patterns of genes are consistent with the physiological data and with the transcription data reported for Fnr in E. coli [11, 12, 20] and with the more limited data by Beliaev et al. . Genes involved in nitrate reduction (napDAHGB, nrfA, and hcp) were significantly down-regulated by the etrA deletion as well as those encoding the fumarate reduction (frdAB, fccA) and all the genes encoding for the DMSO reductases (dmsAB). All of these genes have been considered candidates for EtrA regulation in previous studies; however, results were not conclusive [5–7, 16]. Sequence analysis of the regulatory regions of the differentially expressed genes, indicated possible EtrA recognition sites for most genes in the "Energy metabolism" category, some of which i.e. napDAHGB, nrfA, frdAB and dmsAB, confirms previous results  and further suggests that regulation of these genes is via direct interaction of EtrA with their promoters. Putative recognition sites for EtrA were also identified for the two nqr gene clusters, which had not been identified previously. Also, the regulatory regions for fdh gene clusters were evaluated and an EtrA binding site was recognized for only fdhA-1. The fdh-2 cluster does not possess an EtrA binding site, suggesting a different regulatory system.
Our data indicate that EtrA is a global regulator acting in cooperation with other regulatory proteins to control anaerobic metabolic processes in strain MR-1 [6, 7, 16], therefore, the expression of these genes cannot be expected to be under an "all or none" regulatory mechanism. Rather, these global regulators respond to multiple stimuli (e.g., oxygen levels, substrates) and fine-tune regulation via transcriptional control and interactions between regulatory proteins. Studies in S. oneidensis and in other Shewanella species that indicate the combined action of transcriptional regulators for the anaerobic metabolism in this organism [4, 17–19]. For example, recent studies showed that CRP, EtrA and the product of the cya genes act as expression regulators of several anaerobic respiratory systems, including nitrate reduction in S. oneidensis MR-1 and Shewanella sp. strain ANA-3 [4, 17–19]. In E. coli, Fnr and NarP positively regulate the nap and nrf genes [12, 20, 38, 39]. MR-1 possesses the genes for a homolog of the two-component regulatory system in E. coli NarQ/NarP (SO3981-3982). The presence of alternate regulators that partially fulfill the function of EtrA can explain why nitrate reduction even though impaired, still occurred in the EtrA7-1 knockout mutant.
Down-regulation of genes for lactate transport was also observed. Since lactate was the source of reducing equivalents and carbon, a lack of electron donor and carbon may have contributed to the impaired growth of the EtrA7-1 mutant. Induction of transport proteins for carbon sources and electron acceptors has also been credited to Fnr in E. coli [12, 20], and a putative EtrA binding site was predicted for the gene encoding a lactate permease (SO0827) in MR-1.
Impaired growth of EtrA7-1 could also be due to stress factors caused or enhanced by the deletion (e.g. accumulation of nitrogen oxide reactive species and starvation). The expression of phage-related genes induced in response to irradiation in strain MR-1 has been reported . Up-regulation of the genes involved in activation of the strain MR-1 prophages LambdaSo, MuSo1 and MuSo2 in the EtrA7-1 mutant was observed, suggesting phage activity. Induction of bacterial genes (e.g., nusAG) required to stabilize the Lambda protein antitermination complex in E. coli was also shown [41, 42]. Conversely, there is repression of the LambdaSo transcription regulator Cro/CI family, which represses the transcription of the Lambda genes in E. coli . Also, the induction of genes associated with starvation, i.e., a condition that could activate the lytic cycle of prophages , was confirmed in the expression analysis.
The involvement of several regulatory controls has complicated the interpretation of gene expression patterns and functions in Shewanella spp. Results from the above etrA deletion mutant studies suggest a global regulatory role for EtrA, but one which works in conjunction with other regulators to fine-tune the expression of key genes in anaerobic metabolic pathways in S. oneidensis strain MR-1. Besides confirming and clarifying previous reports on Fnr regulation, we also provide experimental evidence for a positive regulatory role of EtrA in the DMSO reduction pathway of strain MR-1. Furthermore, our whole-genome transcriptional profile shows the effects of EtrA on the expression of genes not previously evaluated (e.g. nqr, fdh-1, phage- and stress-related genes), and differences in the expression pattern of genes previously analyzed (e.g. cydAB and sdhC)[6, 12]. These observations are consistent with results obtained by Gralnick et al.  suggesting a distinctive regulatory system, although very similar to Fnr in E. coli. A stringent sequence analysis of the regulatory region of the genes affected by the mutation suggest direct interaction of EtrA to those in the "Energy metabolism" category, while stress- and phage-related genes are up-regulated indirectly as a consequence of a secondary perturbation. This and previous work taken together suggest that this regulator is more properly termed Fnr.
Bacterial strains and culture conditions
The bacterial strains, plasmids, primers and, probes used in this study are described in Table 4. S. oneidensis strain MR-1 and its mutant strains were grown in HEPES medium as described . The medium was supplemented with 20 mM lactate and KNO3 was added as electron acceptor in concentrations specified below. Oxygen was removed from the medium by boiling and purging with helium . Cultures of E. coli strain β2155 (auxotroph of diaminopimelic acid [DAP]) were grown in Luria-Bertani (LB) medium supplemented with 100 μg/ml of DAP at 37°C. S. oneidensis strain MR-1 was cultivated in aerobic LB medium at 30°C during the mutagenesis process. Antibiotics used for the selection of MR-1 transformants were added in the following concentrations: 25 μg/ml of kanamycin, 7.5 μg/ml of gentamycin, and 10 μg/ml of tetracycline. Vessels that received no inoculum or no KNO3 served as negative controls.
Construction of a chromosomal ΔetrA::loxPallele
PCR primers were designed with Vector NTI® software (InforMax, Inc., Frederick, MD) and synthesized at Integrated DNA Technologies (http://www.idtdna.com).
A deletion mutant of the etrA gene (SO2356) was constructed by allelic replacement as described . Primers etrAN Rev (SacI) and etrAN Fwd (SacII) generated a 520 bp fragment containing about two-thirds of the upstream SO2357 gene, the SO2357-etrA intergenic region, and the first four base pairs of the 753-bp etrA gene (Table 4). Primers etrAC Rev (NdeII) and etrAC Fwd (EcoRI) generated a 526 bp fragment containing the last five base pairs of the etrA gene, the etrA-SO2355 intergenic region and about half of the downstream SO2355 gene (Table 4). The resulting 'SO2357-etrA'-loxP-kan-loxP-'etrA-SO2355' assembly was cloned into the conditionally-replicating plasmid pKNOCK-Gm (Table 4), which encodes resistance to gentamycin, to generate plasmid pCCG02. S. oneidensis MR-1 Kmr Gms colonies were screened by PCR using primers etrAScreenout Fwd and etrAScreenout Rev (Table 4) to identify recombinants in which the etrA gene was replaced by the loxP-kan-loxP cassette.
The kan gene was subsequently removed from the mutant MR-1 genome by Cre recombination [44, 46]. PCR and DNA sequence analyses confirmed the presence of the chromosomal ΔetrA::loxP allele. DNA sequencing was performed at the Genomics Technical Support Facility at Michigan State University.
Plasmid pCM62 (Table 4) was used as the vector for the expression of the etrA gene in a ΔetrA::loxP mutant (strain EtrA7-1). The etrA gene (SO2356) was PCR amplified from S. oneidensis MR-1 genomic DNA using the etrAcomp Fwd (BamHI) and etrAcomp Rev (EcoRI)(Table 4). The amplicon was double digested with BamHI and EcoRI and ligated to the multiple cloning site in pCM62. This construct (pCCG03) was transformed into EtrA7-1 by conjugation from E. coli β2155. Ligation, electroporation into E. coli β2155, and conjugation in strain EtrA7-1 were performed as described . Plasmid pCM62 was also transformed into EtrA7-1 via conjugation from E. coli β2155 and used as a control for any plasmid effects. Transformants were selected by streaking on LB plates with tetracycline. EtrA7-1 Tcr colonies were diagnosed by PCR using the etrAcomp primers (Table 4) and subsequently sequenced to verify the deletion of the etrA gene.
Phenotypic characterization of the ΔetrA::loxPmutant
Cultures of the wild type, EtrA7-1, EtrA7-1 complement and EtrA7-1 harboring pCM62 were grown anaerobically with 3 mM KNO3 in HEPES medium. Growth was monitored periodically by OD measurements at 600 nm. Samples (2 mL) were periodically withdrawn for analysis of nitrate, nitrite and ammonium concentrations as described [44, 47].
Cultures of the wild type and EtrA7-1 were also cultivated anaerobically with ferric citrate (10 mM), fumarate (10 mM), disodium thiosulfate (10 mM), trimethylamine N-oxide (TMAO; 10 mM), manganese dioxide (1 mM, nominal concentration), dimethyl sulfoxide (DMSO; 2 and 10 mM) and disodium sulfite (1 mM), as electron acceptors. The ferric citrate and the manganese dioxide were prepared as described . Evidence of growth via reduction of TMAO, thiosulfate and fumarate was determined by OD600 measurements. Fe(III) reduction was determined by the ferrozine assay following HCl extraction [49, 50]. Mn(IV) reduction was assayed colorimetrically . Cultures supplied with DMSO as the terminal electron acceptor were analyzed by high-performance liquid chromatography (HPLC) for lactate consumption and acetate formation . Sulfite consumption was measured using a DX-100 ion chromatograph (Dionex Corp., Sunnyvale, CA) equipped with an IonPac AS14A Column.
To determine the effects of lactate on the reduction of DMSO, nitrate and fumarate, cultures of the wild type and the EtrA7-1 mutant strain were grown anaerobically with 20 mM sodium pyruvate as the electron donor and dimethyl sulfoxide (DMSO; 1 mM), fumarate (10 mM) or nitrate (2 mM) as electron acceptors. DMSO and fumarate reduction were monitored as mentioned above. Nitrate reduction was measured using a Dionex ICS-3000 ion chromatograph (Dionex Corp., Sunnyvale, CA) equipped with an IonPac AS14 Column.
Resting cell assays
Cells of the wild type strain MR-1 and the EtrA7-1 mutant were grown aerobically in HEPES medium with 20 mM lactate as the electron donor to OD600 = 0.2. Cells were harvested by centrifuging for 10 min at 3,000 g, washed twice with 50 mM saline phosphate buffer (pH 7.0) [52, 53], and resuspended in the same buffer to an OD600 = 1.0 under anoxic conditions. The cells were incubated with 5 mM KNO3, and after 0, 1, 2, 4, 8 and 24 h were removed by centrifugation and three 1-mL replicate samples of the supernatant were assayed to determine nitrate and nitrite reduction rates. Assays with autoclaved wild type cells served as negative controls. Nitrate, nitrite and ammonium concentrations were determined as described .
Total RNA preparations
Total RNA was extracted from triplicate cultures of strains MR-1 and EtrA7-1 grown with 2 mM nitrate as the sole electron acceptor. The RNA was extracted with RNAwiz Solution following the instructions of the manufacturer (Ambion, Inc., Austin, TX). RNA samples were treated with RNase-free DNaseI (Roche Pharmaceuticals, Basel, Switzerland) and purified by phenol:chloroform (1:1) and chloroform extractions , and stored in ethanol at -80°C until use. Quality of the RNA was verified using the RNA 6000 Pico LabChip kit and the 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA).
Global expression analyses
A S. oneidensis strain MR-1 whole genome microarrays  were provided by Liyou Wu and Jizhong Zhou (Oak Ridge National Laboratory, Oak Ridge, TN). cDNA preparation and labeling were performed as described  using a 2:3 ratio of 5-(3-aminoallyl)-dUTP and dTTP. Hybridization and post-hybridization washes were done as described . Three biological replicates per treatment were used for the hybridization of six microarray slides including technical duplicates (dye-swap). Data analysis was performed using the GeneSpring 6.0 software (Silicon Genetics, Redwood City, CA). The data were normalized per chip and per gene (Lowess Normalization) and the spots with less than 55% pixel intensity above background plus two standard deviations were eliminated from the analyses . The data were filtered using the Benjamini and Hochberg false discovery rate with 95% confidence and only those genes with a > 2-fold change in expression were considered significant.
Microarray data accession number
The raw microarray intensity data has been deposited in the GenBank Gene Expression Omnibus (GEO) database under the accession number GSE26935.
Identification of putative EtrA binding sites
Regulatory motifs were predicted in the intergenic regions of differentially expressed genes using the Gibbs centroid sampler . Intergenic regions were extracted, based on the S. oneidensis MR-1 genome annotation, that were at least 50 bp in length and upstream of differentially expressed genes or operons whose change in expression (average ± one standard deviation) was at least 2.5-fold. A total of 170 genes met these requirements: 55 had expression levels minus one standard deviation of ≥ 2.5 (up-regulated genes) and 115 had expression levels plus one standard deviation of ≤ 0.4 (down-regulated genes). For these genes, 118 upstream intergenic regions were extracted for analysis, after accounting for multiple genes within an operon. The parameters for the Gibbs centroid sampler used on these sequences were the following: up to two motif models were allowed, where each model was specified to be palindromic and 16-24 bases long, a maximum of three sites per intergenic was allowed, a position-specific background model  was employed, and centroid sampling was performed with 1000 burn-in iterations, 5000 sampling iterations and 10 random seeds. The results from four independent runs were compared, and the subset of 47 intergenic regions extracted that contained a predicted regulatory motif in at least one of those runs. These 47 intergenic sequences were analyzed with the Gibbs centroid sampler, using the same parameters as above, except that only one motif model was specified. Additional binding sites were detected using dscan (http://ccmbweb.ccv.brown.edu/cgi-bin/dscan.pl) to search the set of promoters for all the genes that exhibited ≥ 2-fold change in expression (Additional file 1). This set included a total of 424 intergenic regions.
Holden M, Bentley S, Sebaihia M, Thompson N, Cerdeño-Tárraga A, Parkhill J: The magnificent seven. Trends Microbiol. 2003, 11: 12-14.22. 10.1016/S0966-842X(02)00009-4.
Tiedje JM: Shewanella-the environmentally versatile genome. Nat Biotechnol. 2002, 20: 1093-1094. 10.1038/nbt1102-1093.
Heidelberg JF, Paulsen IT, Nelson KE, Gaidos EJ, Nelson WC, Read TD, Eisen JA, Seshadri R, Ward N, Methe B, Clayton RA, Meyer T, Tsapin A, Scott J, Beanan M, Brinkac L, Daugherty S, DeBoy RT, Dodson RJ, Durkin , Haft DH, Kolonay JF, Madupu R, Peterson JD, Umayam LA, White O, Wolf AM, Vamathevan J, Weidman J, Impraim M, Lee K, Berry K, Lee C, Mueller J, Khouri H, Gill J, Utterback TR, McDonald LA, Feldblyum TV, Smith HO, Venter JC, Nealson KH, Fraser CM: Genome sequence of the dissimilatory metal ion-reducing bacterium Shewanella oneidensis. Nat Biotechnol. 2002, 20: 1118-1123. 10.1038/nbt749.
Gralnick JA, Brown CT, Newman DK: Anaerobic regulation by an atypical Arc system in Shewanella oneidensis. Mol Microbiol. 2005, 56: 1347-1357. 10.1111/j.1365-2958.2005.04628.x.
Saffarini DA, Schultz R, Beliaev A: Involvement of cyclic AMP (cAMP) and cAMP receptor protein in anaerobic respiration of Shewanella oneidensis. J Bacteriol. 2003, 185: 3668-3671. 10.1128/JB.185.12.3668-3671.2003.
Beliaev AS, Thompson DK, Fields MW, Wu L, Lies DP, Nealson KH, Zhou J: Microarray transcription profiling of a Shewanella oneidensis etrA mutant. J Bacteriol. 2002, 184: 4612-4616. 10.1128/JB.184.16.4612-4616.2002.
Maier TM, Myers CR: Isolation and characterization of a Shewanella putrefaciens MR-1 electron transport regulator etrA mutant: reassessment of the role of EtrA. J Bacteriol. 2001, 183: 4918-4926. 10.1128/JB.183.16.4918-4926.2001.
Darwin AJ, Ziegelhoffer EC, Kiley PJ, Stewart V: Fnr, NarP, and NarL regulation of E. coli K-12 napF (periplasmic nitrate reductase) operon transcription in vitro. J Bacteriol. 1998, 180: 4192-4198.
Darwin AJ, Stewart V: Nitrate and nitrite regulation of the Fnr-dependent aeg-46.5 promoter of E. coli K-12 is mediated by competition between homologous response regulators (NarL and NarP) for a common DNA-binding site. J Mol Biol. 1997, 251: 15-29. 10.1006/jmbi.1995.0412.
Choe M, Reznikoff WS: Anaerobically expressed Escherichia coli genes identified by operon fusion techniques. J Bacteriol. 1991, 173: 6139-6146.
Shalel-Levanon S, San KY, Bennett GN: Effect of ArcA and FNR on the expression of genes related to the oxygen regulation and the glycolysis pathway in Escherichia coli under microaerobic growth conditions. Metab Eng. 2005, 7: 445-456. 10.1016/j.ymben.2005.07.001.
Kang Y, Weber KD, Qiu Y, Kiley PJ, Blattner FR: Genome-wide expression analysis indicates that FNR of Escherichia coli K-12 regulates a large number of genes of unknown function. J Bacteriol. 2005, 187: 1135-1160. 10.1128/JB.187.3.1135-1160.2005.
Crack J, Green J, Thomson AJ: Mechanism of oxygen sensing by the bacterial transcription factor fumarate-nitrate reduction (FNR). J Biol Chem. 2004, 279: 9278-9286. 10.1074/jbc.M309878200.
Outten FW: Iron-sulfur clusters as oxygen-responsive molecular switches. Nat Chem Biol. 2007, 3: 206-207. 10.1038/nchembio0407-206.
Moore LJ, Mettert EL, Kiley PJ: Regulation of FNR dimerization by subunit charge repulsion. J Biol Chem. 2006, 281: 33268-33275. 10.1074/jbc.M608331200.
Saffarini DA, Nealson KH: Sequence and genetic characterization of etrA, an fnr analog that regulates anaerobic respiration in Shewanella putrefaciens MR-1. J Bacteriol. 1993, 175: 7938-7944.
Charania MA, Brockman KL, Zhang Y, Banerjee A, Pinchuk GE, Fredrickson JK, Beliaev AS, Saffarini DA: Involvement of a membrane-bound class III adenylate cyclase in regulation of anaerobic respiration in Shewanella oneidensis MR-1. J Bacteriol. 2009, 191: 4298-4306. 10.1128/JB.01829-08.
Murphy JN, Saltikov CW: The ArsR repressor mediates arsenite-dependent regulation of arsenate respiration and detoxification operons of Shewanella sp. Strain ANA-3. J Bacteriol. 2009, 191: 6722-6731. 10.1128/JB.00801-09.
Murphy JN, Durbin KJ, Saltikov CW: Functional roles of arcA, etrA, cyclic AMP (cAMP)-cAMP receptor protein, and cya in the arsenate respiration pathway in Shewanella sp. strain ANA-3. J Bacteriol. 2009, 191: 1035-1043. 10.1128/JB.01293-08.
Overton TW, Griffiths L, Patel MD, Hobman JL, Penn CW, Cole JA, Constantinidou C: Microarray analysis of gene regulation by oxygen, nitrate, nitrite, FNR, NarL and NarP during anaerobic growth of Escherichia coli: new insights into microbial physiology. Biochem Soc T. 2006, 34: 104-107. 10.1042/BST0340104.
Myers CR, Myers JM: Role of the tetraheme cytochrome CymA in anaerobic electron transport in cells of Shewanella putrefaciens MR-1 with normal levels of menaquinone. J Bacteriol. 2000, 182: 67-75. 10.1128/JB.182.1.67-75.2000.
Fredrickson JK, Romine MF: Genome-assisted analysis of dissimilatory metal-reducing bacteria. Curr Opin Biotechnol. 2005, 16: 269-274. 10.1016/j.copbio.2005.04.001.
Hayashi M, Hirai K, Unemoto T: Sequencing and the alignment of structural genes in the nqr operon encoding the Na(+)-translocating NADH-quinone reductase from Vibrio alginolyticus. FEBS Lett. 1995, 363: 75-77. 10.1016/0014-5793(95)00283-F.
Filenko NA, Browning DF, Cole JA: Transcriptional regulation of a hybrid cluster (prismane) protein. Biochem Soc Trans. 2005, 33: 195-197. 10.1042/BST0330195.
Cabello P, Pino C, Olmo-Mira MF, Castillo F, Roldán MD, Moreno-Vivián C: Hydroxylamine assimilation by Rhodobacter capsulatus E1F1. requirement of the hcp gene (hybrid cluster protein) located in the nitrate assimilation nas gene region for hydroxylamine reduction. J Biol Chem. 2004, 279: 45485-45494. 10.1074/jbc.M404417200.
Wolfe MT, Heo J, Garavelli JS, Ludden PW: Hydroxylamine reductase activity of the hybrid cluster protein from Escherichia coli. J Bacteriol. 2002, 184: 5898-5902. 10.1128/JB.184.21.5898-5902.2002.
van den Berg WA, Hagen WR, van Dongen WM: The hybrid-cluster protein ('prismane protein') from Escherichia coli. Characterization of the hybrid-cluster protein, redox properties of the [2Fe-2S] and [4Fe-2S-2O] clusters and identification of an associated NADH oxidoreductase containing FAD and [2Fe-2S]. Eur J Biochem. 2000, 267: 666-676. 10.1046/j.1432-1327.2000.01032.x.
White AK, Metcalf WW: The htx and ptx operons of Pseudomonas stutzeri WM88 are new members of the pho regulon. J Bacteriol. 2004, 186: 5876-5882. 10.1128/JB.186.17.5876-5882.2004.
White AK, Metcalf WW: Isolation and biochemical characterization of hypophosphite/2-oxoglutarate dioxygenase. A novel phosphorus-oxidizing enzyme from Pseudomonas stutzeri WM88. J Biol Chem. 2002, 277: 38262-38271. 10.1074/jbc.M204605200.
Jiang ZD, Greenberg D, Natarro JP, Stephen R, DuPont HL: Rate of occurrence and pathogenic effect of enteroaggregative Escherichia coli virulence factors in international travelers. J Clin Microbiol. 2002, 40: 4185-4190. 10.1128/JCM.40.11.4185-4190.2002.
Jiang W, Metcalf WW, Lee KS, Wanner BL: Molecular cloning, mapping, and regulation of Pho regulon genes for phosphonate breakdown by the phosphonatase pathway of Salmonella typhimurium LT2. J Bacteriol. 1995, 177: 6411-6421.
Kim SK, Makino K, Amemura M, Nakata A, Shinagawa H: Mutational analysis of the role of the first helix region 4.2 of the sigma 70 subunit of Escherichia coli RNA polymerase in transcriptional activation by activator protein PhoB. Mol Gen Genet. 1995, 248: 1-8. 10.1007/BF02456607.
Hansen AM, Gu Y, Li M, Andrykovitch M, Waugh DS, Jin DJ, Ji X: Structural basis for the function of stringent starvation protein a as a transcription factor. J Biol Chem. 2005, 280: 17380-17391. 10.1074/jbc.M501444200.
Jovanovic G, Weiner L, Model P: Identification, nucleotide sequence, and characterization of PspF, the transcriptional activator of the Escherichia coli stress-induced psp operon. J Bacteriol. 1996, 178: 1936-1945.
Klancnik A, Bottledoorn N, Herman L, Mozina SS: Survival and stress induced expression of groEL and rpoD of Campylobacter jejuni from different growth phases. Int J Food Microbiol. 2006, 112: 200-207. 10.1016/j.ijfoodmicro.2006.03.015.
Pillutla RC, Ahnn J, Inouye M: Deletion of the putative effector region of Era, an essential GTP-binding protein in Escherichia coli, causes a dominant-negative phenotype. FEMS Microbiol Lett. 1996, 143: 47-55. 10.1111/j.1574-6968.1996.tb08460.x.
Luisi-DeLuca C, Kolodner R: Purification and characterization of the Escherichia coli RecO protein. J Mol Biol. 1994, 236: 124-138. 10.1006/jmbi.1994.1123.
Cotter PA, Gunsalus RP: Oxygen, nitrate and molybdenum regulation of dmsABC genes expression in Escherichia coli. J Bacteriol. 1989, 171: 3817-3823.
Stewart V, Bledsoe PJ, Williams SB: Dual overlapping promoters control napF (periplasmic nitrate reductase) operon expression in Escherichia coli K-12. J Bacteriol. 2003, 185: 5862-5870. 10.1128/JB.185.19.5862-5870.2003.
Qiu X, Sundin GW, Wu L, Zhou J, Tiedje JM: Comparative analysis of differentially expressed genes in Shewanella oneidensis MR-1 following exposure to UVC, UVB, and UVA radiation. J Bacteriol. 2005, 187: 3556-3564. 10.1128/JB.187.10.3556-3564.2005.
Bonin I, Muhlberger R, Bourenkov GP, Huber R, Bacher A, Richter G, Wahl WC: Structural basis for the interaction of Escherichia coli NusA with protein N of phage lambda. Proc Natl Acad Sci USA. 2004, 101: 13762-13767. 10.1073/pnas.0405883101.
Torres M, Balada JM, Zellars M, Squires C, Squires C: In vivo effect of NusB and NusG on rRNA transcription antitermination. J Bacteriol. 2004, 186: 1304-1310. 10.1128/JB.186.5.1304-1310.2004.
Voyles BA: The biology of viruses. 1993, Mosby-Year Book, Inc., St. Louis, MO
Cruz-García C, Murray AE, Klappenbach NJA, Stewart V, Tiedje JM: Respiratory Nitrate Ammonification by Shewanella oneidensis MR-1. J Bacteriol. 2007, 189: 656-662.
Balch WE, Wolfe RS: New approach to the cultivation of methanogenic bacteria: 2-mercaptoethanesulfonic acid (HS-CoM)-dependent growth of Methanobacterium ruminantium in a pressureized atmosphere. Appl Environ Microbiol. 1976, 32: 781-791.
Marx CJ, Lidstrom ME: Broad-host-range cre-lox system for antibiotic marker recycling in gram-negative bacteria. Biotechniques. 2002, 33: 1062-1067.
Nelson DW: Determination of ammonium in KCl extracts of soils by the salicylate method. Comm Soil Sci Plant Anal. 1983, 14: 1051-1062. 10.1080/00103628309367431.
Burlage RS, Atlas R, Stahl D, Gessey G, Sayler G: Techniques in Microbial Ecology. 1988, Oxford University Press US, New York, NY
Lovley DR, Phillips EJP: Novel Mode of Microbial Energy Metabolism: Organic Carbon oxidation Coupled to Dissimilatory Reduction of Iron or Manganese. Appl Environ Microbiol. 1988, 54: 1472-1480.
Lovley DR, Phillips EJP: Availability of Ferric Iron Microbial Reduction in Bottom Sediments of the Freshwater Tidal Potomac River. Appl Environ Microbiol. 1986, 52: 751-757.
He J, Ritalahti RM, Aiello MR, Löffler FE: Complete Detoxification of Vinyl Chloride by an Anaerobic Enrichment Culture and Identification of the Reductively Dechlorinating Population as a Dehalococcoides species. Appl Environ Microbiol. 2003, 69: 996-1003. 10.1128/AEM.69.2.996-1003.2003.
Rodrigues JLM, Maltseva OV, Tsoi TV, Helton RR, Quensen JF, Fukuda M, Tiedje JM: Development of a Rhodococcus recombinant strain for degradation of products from anaerobic dechlorination of PCBs. Environ Sci Technol. 2001, 35: 663-668. 10.1021/es001308t.
Löffler FE, Champine JE, Ritalahti KM, Sprague SJ, Tiedje JM: Complete reductive dechlorination of 1,2-dichloropropane by anaerobic bacteria. Appl Environ Microbiol. 1997, 63: 2870-2875.
Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY
Gao H, Wang Y, Liu X, Yan T, Wu L, Alm E, Arkin A, Thompson DK, Zhou J: Global transcriptome analysis of the heat shock response of Shewanella oneidensis. J Bacteriol. 2004, 186: 7796-7803. 10.1128/JB.186.22.7796-7803.2004.
Schroeder RG, Peterson LM, Fleischmann RD: Improved quantitation and reproducibility in Mycobacterium tuberculosis DNA microarrays. J Mol Microbiol Biotechnol. 2002, 4: 123-126.
Hedge P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Earlehughes J, Snesrud E, Lee N, Quackenbush J: A concise guide to cDNA microarray analysis. BioTechniques. 2000, 29: 548-562.
Murray AE, Lies D, Li G, Nealson KH, Zhou J, Tiedje JM: DNA/DNA hybridization to microarrays reveals gene-specific differences between closely related microbial genomes. Proc Natl Acad Sci. 2001, 98: 9853-9858. 10.1073/pnas.171178898.
Thompson WA, Newberg LA, Conlan S, McCue L, Lawrence CE: The Gibbs centroid sampler. Nucleic Acids Res. 2007, 35: W232-237. 10.1093/nar/gkm265.
Liu JS, Lawrence CE: Bayesian inference on biopolymer models. Bioinformatics. 1999, 15: 38-52. 10.1093/bioinformatics/15.1.38.
Demarre G, Guérout AM, Matsumoto-Mashimo C, Rowe-Magnus DA, Marlière P, Mazel D: A new family of mobilizable suicide plasmids based on broad host range R388 plasmid (IncW) and RP4 plasmid (IncPα) conjugative machineries and their cognate Escherichia coli host strains. Res Microbiol. 2005, 156: 245-255.
Myers CR, Nealson KH: Bacterial manganese reduction and growth with manganese oxide as the sole electron acceptor. Science. 1988, 240: 1319-1321. 10.1126/science.240.4857.1319.
Marx CJ, Chistoserdova L: Development of versatile broad-host-range vectors for use in methylotrophs and other gram-negative bacteria. Microbiology. 2001, 147: 2065-2075.
Alexeyev MF: The pKNOCK series of broad-host-range mobilizable suicide vectors for gene knockout and targeted DNA insertion into the chromosome of gram-negative bacteria. BioTechniques. 1999, 26: 824-828.
We thank Xiaoyun Qiu for advice on the DNA microarray work, Valley Stewart and Joel Klappenbach for advice and discussion. We thank Benjamin K. Amos, Jed Costanza, Qingzhong Wu and Sara H. Thomas for technical assistance in the phenotypic characterization of the EtrA7-1 strain. We also acknowledge members of the Shewanella Federation for helpful discussions.
This study was supported by Department of Energy grants DE-FG02-02ER63342 from the Genomics Program, Office of Biological and Environmental Research (awarded to JMT), DE-FG02-04ER63718.25 from the Environmental Remediation Science Division, Biological and Environmental Research (awarded to FEL) and DE-FG02-04ER63942 from the Genomes to Life Program, Office of Biological and Environmental Research (awarded to LAM). Contributions by MFR and LAM were performed at Pacific Northwest National Laboratory, which is operated by Battelle for the United States Department of Energy under Contract DE-AC05-76RL01830.
All authors contributed in the organization and design of experiments as well as data interpretation and manuscript preparation. CCG, FEL, and JMT wrote the paper. CCG designed and carried out the majority of the experimental work including mutant construction, cDNA microarray experiments and analysis, and growth studies. AEM, MFR and LAM contributed in experimental design and cDNA microarray data analysis and interpretation. JLMR performed resting cell assays. JAG replicated DMSO physiological experiments. LAM performed EtrA binding site identification. MFR provided updated genome sequence annotation. FEL provided laboratory equipment, materials, and funding and supervision for the phenotypic characterization work. JMT supervised experimental work. All authors read and approved the final version of the manuscript.