Systematic mapping of two component response regulators to gene targets in a model sulfate reducing bacterium
© Rajeev et al.; licensee BioMed Central Ltd. 2011
Received: 11 April 2011
Accepted: 12 October 2011
Published: 12 October 2011
Two component regulatory systems are the primary form of signal transduction in bacteria. Although genomic binding sites have been determined for several eukaryotic and bacterial transcription factors, comprehensive identification of gene targets of two component response regulators remains challenging due to the lack of knowledge of the signals required for their activation. We focused our study on Desulfovibrio vulgaris Hildenborough, a sulfate reducing bacterium that encodes unusually diverse and largely uncharacterized two component signal transduction systems.
We report the first systematic mapping of the genes regulated by all transcriptionally acting response regulators in a single bacterium. Our results enabled functional predictions for several response regulators and include key processes of carbon, nitrogen and energy metabolism, cell motility and biofilm formation, and responses to stresses such as nitrite, low potassium and phosphate starvation. Our study also led to the prediction of new genes and regulatory networks, which found corroboration in a compendium of transcriptome data available for D. vulgaris. For several regulators we predicted and experimentally verified the binding site motifs, most of which were discovered as part of this study.
The gene targets identified for the response regulators allowed strong functional predictions to be made for the corresponding two component systems. By tracking the D. vulgaris regulators and their motifs outside the Desulfovibrio spp. we provide testable hypotheses regarding the functions of orthologous regulators in other organisms. The in vitro array based method optimized here is generally applicable for the study of such systems in all organisms.
Signal transduction to sense and respond to environmental and intracellular changes governs key cellular regulatory functions. In bacteria, two component systems, composed typically of a sensor histidine kinase (HK) and a response regulator (RR), are the primary and best-studied mechanisms for perceiving such changes and controlling downstream responses [1–3]. The regulatory network of an organism is often a reflection of the environments in which it can survive and the signal transduction systems in microbes have been correlated to their sensory IQs . Desulfovibrio vulgaris Hildenborough, an anaerobic sulfate reducing bacterium, occupies a variety of ecological niches and encodes a strikingly large number of these systems with unusual diversity attributed to lineage-specific expansion of existing gene families . Studied since the 1940s, D. vulgaris Hildenborough has come to serve as a model system to evaluate dissimilatory sulfate reduction and hydrogen cycling . However, the function of none of its two component systems, encoded by 72 RRs and 64 HKs, have been characterized to date.
The distribution of RRs in D. vulgaris Hildenborough is considerably different from in other microbes. Of the 72 RRs in D. vulgaris Hildenborough, 29 have a DNA binding domain (DBD), indicating function via gene regulation. Twenty-two of these fall into the NtrC family of σ54-dependent RRs. σ54-dependent response regulators in bacteria typically make up approximately 9% of the total RRs in most organisms  but in D. vulgaris Hildenborough this group constitutes 30% of the total RRs, and 75% of the ones with DBDs. On the other hand, the OmpR family, which typically constitutes the most abundant class of RRs in bacteria, has only two representatives in D. vulgaris Hildenborough. The remaining five RRs fall into the LytR and NarL families (Table S1 in Additional file 1). With the exception of DVU1083, which is an ortholog of the Escherichia coli PhoB , none of the RRs have any characterized orthologs. The targets of these 29 RRs represent the transcriptional portion of the two component regulatory network of this organism and to date remains almost entirely undetermined.
With the exception of a few model organisms such as E. coli [7, 8], Caulobacter crescentus , and Bacillus subtilis [10, 11], genes regulated by these systems remain largely unmapped in most organisms. Even in model bacteria, systematic approaches to delineate gene targets and regulatory networks controlled by two component systems are rare and the available knowledge of their networks represents information compiled from a large body of literature, in silico efforts [8, 12], or by indirect inference of targets based on transcriptomics analysis . While mapping of binding sites via ChIP-on-chip assays is now done routinely for transcription factors, it is effective for two component RRs only if the activating signal or conditions are known. As a result, even in E. coli and B. subtilis the function and targets of some two-component systems remain unmapped.
Here we present a systematic experimental determination of the genes regulated by the transcriptionally acting RRs in D. vulgaris Hildenborough. We optimized an in vitro approach in order to bypass the requirement of using activating conditions that are largely unknown for these two component systems. To our knowledge, this is the first extensive use of an in vitro genome-wide method to map bacterial two component system RR binding sites.
Results and discussion
Gene targets were determined for 24 D. vulgaris Hildenborough RRs
Parameters for RR binding to the sheared genomic DNA were determined using the gel-shift assays and enrichment of the positive control target in the RR-bound DNA fraction was confirmed using quantitative PCR (qPCR; Table S3 in Additional file 1). Successful enrichment of the positive control and no enrichment of a non-specific negative control (Figure 2c) also serves as a validation of the specificity of binding seen in the gel-shift assays, and increases the confidence in the subsequent DAP-chip data set. The chip-based measurements were then conducted as described in Materials and methods. Nimblescan software was used to analyze the tiling array data and rank enriched gene loci for each RR. The top 20 peaks obtained for each DAP-chip are provided in Table S4 in Additional file 1.
Determination of binding site motifs
Our approach proved to be very successful. Binding sites were predicted and confirmed for 15 hitherto uncharacterized RRs (see Figure 8 for the motifs and Figures 4, 5 and 6 for binding site distribution within targets). Experimentally validated motifs further confirm the specificity of the peaks discovered using the DAP-chip method. The majority of the binding sites are palindromic, ranging from 4- to 6-bp inverted repeats separated by 3 to 8 bp in between. In two cases, DVU2394 and DVU1083, the binding site was found to be a direct repeat. Interestingly, RRs DVU0539 and DVU0946, which are paralogs, also recognized the same binding site (Figures 8 and 9).
The confirmed binding motifs were used to assess the general applicability and robustness of this method. First we examined if the binding sequences were present in all the hits for an RR in its DAP-chip data set. Our analyses indicate that for the RRs where a motif was determined, a single motif explained all the targets found in the DAP-chip analysis. The primary exception was DVU1083 (PhoB), where the motif was found upstream of only 14 of the 30 targets. The other exceptions were DVU2394 and DVU0539, where the discovered motif was only present in a subset of the DAP-chip hits (Table S8 in Additional file 1).
A genome-wide scan also revealed that the DAP-chip assay successfully enriched all loci containing the binding site motif for a given RR. An exception to this general observation was the RR DVU1063 (the flagella regulator, see below), for which a sizable number of potential sequences in upstream regions were present that were not enriched as targets in our assay (Table S9 in Additional file 1). Interestingly, many of these are flagella and motility-related genes, suggesting that they are real targets. While it is not clear why these sequences were not enriched in the in vitro assay, one possibility is that the quantity of protein or the provided activated state used was insufficient for interaction with all available loci. DNA modification such as via methylation is another source of regulation that would persist in an in vitro assay and affect motif recognition. Additional experiments would be required to examine these hypotheses.
Functional assignments for response regulators
The DAP-chip gene hits enabled predictions for the function of several RRs (Figure 7). Additional genomics data available for D. vulgaris Hildenborough from transcriptomics studies under different stresses and conditions (MicrobesOnline) and high-resolution tiling arrays conducted to identify transcribed regions and estimate their level of expression  were also used to obtain support for conclusions from our measurements. Using the transcriptomics data, RRs and their target operons were examined for any obvious co-expression patterns while the tiling array data provided support for the expression levels of a gene under routine growth conditions.
Lactate utilization is highly regulated
The lactate/pyruvate utilization genes of the operon DVU3025-3033 are regulated by four RRs - DVU3023, DVU0539, DVU0621, and DVU1083 (PhoB). Lactate is the primary electron donor and carbon source for Desulfovibrio spp., and the genome of D. vulgaris Hildenborough encodes several lactate permeases, lactate dehydrogenases, and pyruvate oxidation genes. Since DVU3025-3033 is so highly regulated, it is likely that it is the primary pathway for lactate utilization. These genes were also observed to be highly expressed in the tiling array data (Figure S1 in Additional file 3) . This operon additionally contains the ack gene encoding acetate kinase, which generates acetate from acetyl phosphate, the energy generating step in dissimilatory sulfate reduction. Since acetyl phosphate can act as a small molecule phosphate donor for RRs , having multiple regulators for this operon may present a mechanism to modulate acetyl phosphate levels inside the cell.
The targets for the RR DVU3023 are lactate-pyruvate oxidation genes DVU3025-3033, and two lactate permeases (Figure 6), suggesting that the corresponding two component system senses lactate and also plays a role in lactate utilization. Further, the two lactate permeases targeted by DVU3023 are expressed differently. Tiling array data indicate that DVU2451 is expressed highly under normal growing conditions, whereas DVU3284 is not (Figure S1 in Additional file 3) . Gene expression correlations show that the two genes are negatively correlated with each other (Figure S2 in Additional file 3). DVU3023 may be activating DVU2451 expression, while repressing DVU3284.
Additionally the lactate permease gene DVU2451 is also targeted by three RRs -DVU3023, DVU0539 and DVU2588 (Figure 6). These three RRs, along with DVU0946, have common targets that may form a network to fine-tune lactate utilization (see below). Thus, lactate consumption appears to be affected in response to multiple stresses or environmental signals.
D. vulgaris encodes paralogous RRs with overlapping functions
The RRs DVU0539 and DVU0946 are paralogous and have identical binding sites (Figures 8 and 9), although their corresponding predicted proximal sensor kinases are not paralogous. The two RRs appear to auto-regulate their own expression, and to regulate each other, and the operon DVU2133-2132. DVU0539 additionally regulates DVU3025-3033, and DVU2451 (Figure 6). Further, within the regulated candidates, genes in the operon DVU2133-2132 appear to be paralogs of those in DVU0943-0944. A gene expression correlation map of these candidates shows that DVU0539-0540 transcript levels positively correlate with DVU0542-0545, and DVU2133-2132, but is negatively correlated with DVU0943-0946 and with DVU3025-3033 and DVU2451 (Figure S2 in Additional file 3). Tiling array data support this observation, where during regular growth DVU0943-0946 genes are well expressed, whereas no expression of DVU0539-0545 genes was measured (Figure S1 in Additional file 3) . Since DVU3025-3033 and DVU2451 are lactate utilization genes, it is likely that the functions of the other target genes encoding hypothetical proteins are also tied in to lactate/pyruvate utilization. Taken together these RRs and their target genes present a highly interconnected and feedback controlled regulatory module to control lactate utilization. Despite identical binding motifs used for DNA binding, our findings suggest that DVU0539 and DVU0946 regulate genes differently. Additional experiments, such as the biochemical evaluation of the phosphotransfer from the respective histidine kinases, may shed more light on the mechanism of such regulation.
Regulation of lipid A biosynthesis
DVU2934 has a single target, the lpxC gene DVU2917 (Figure 4). LpxC is predicted to catalyze the committed step in lipid A biosynthesis. In E. coli and Pseudomonas, regulation appears to be primarily via control of LpxC protein levels. Excess LpxC in these systems is toxic to the cell, although lpxC is an essential gene [21, 22]. Tiling array data for D. vulgaris  suggests that lpxC is highly expressed (Figure S1 in Additional file 3). Regulation by DVU2934 may be an additional mechanism to fine-tune its expression. DVU2934 is part of an operon that encodes the histidine kinase DVU2931 and a second RR with a HD-GYP domain (DVU2933), suggesting that cyclic-di-GMP levels may be used to regulate lipid A biosynthesis in D. vulgaris.
The phosphate starvation response ties in to DNA replication, nitrogen metabolism, and cyclic-di-GMP levels
DVU1083 is annotated as the phosphate starvation response regulator PhoB, and the DAP-chip data support this prediction. Aside from the expected phosphate ABC transporter genes and phosphate transport regulators (PhoU), its targets include DNA polymerase and gyrase genes, amino acid transport genes, glutamate synthase, phosphodiesterase domain (HD/EAL) genes and also some energy metabolism genes, such as the alcohol dehydrogenase and the pyruvate oxidation operon (this also includes the acetate kinase and the phosphotransacetylase (pta-ack) genes) (Figure 6). The Pho regulon in other organisms is known to include members that function outside the direct phosphate starvation response, such as in nitrogen assimilation , DNA replication , and cyclic-di-GMP concentration . The binding site for known PhoB boxes in other bacteria, particularly E. coli, consists of two 7-bp direct repeats (consensus CTGTCAT) separated by a 4-bp spacer [25, 26]. The D. vulgaris PhoB box consensus is a 6-bp repeat with a 4-bp spacer (c/t)GT(n)AC (Figure 8). DVU1083 regulates the DVU2477-2479 operon that encodes the pstS-C-A phosphate ABC transporter genes. The D. vulgaris genome also has another set of phosphate transporter genes, DVU2667-2663 (pstS-C-A-B-ATPase), but interestingly these genes do not have the binding site motif in their upstream regions and were not among the peaks for the PhoB DAP-chip assay.
DVU2114 regulates pili assembly
DVU2114 targets the pilin genes DORF39640-DVU2116, and the operon DVU3342-3345 (Figure 4), which encodes pilus assembly genes . Pili are regulated by σ54-dependent RRs in other species such as Geobacter , Myxococcus xanthus , and Pseudomonas . The type of pilin encoded by DORF39640-DVU2116, and the genes downstream of DVU2116 (DVU2117-2126) appear most similar to the pili assembly machinery (cpaA-F, Tad genes) found in C. crescentus .
DVU1063 appears to be a flagella regulator
DVU1063 has a large number of targets for a σ54-dependent regulator, and among these are some flagella-related genes (Figure 4), suggesting its role as a flagella regulator. DVU1063 is homologous to the flagella regulator FlbD in C. crescentus , and the binding site motif (GGCAxxxxTGCC) resembles that of the C. crescentus FlbD (CCCGGCAxxxxxTGCCGGG), where the italic bases are those that form contacts directly with the RR . Scanning the C. crescentus genome with the D. vulgaris motif identified several of the FlbD-regulated promoters (not shown). Like FlbD, DVU1063 has an atypical receiver domain that lacks some of the active site residues of the phosphorylation pocket, and it may not require activation by phosphorylation. The regulation of cell motility appears to be complex, since other targets for RR DVU1063 include hypothetical proteins, regulatory proteins, membrane proteins and transporters.
Exopolysaccharide and biofilm synthesis is controlled by a pDV1 plasmid encoded regulator
DVUA0057 has been predicted to regulate genes encoding proteins with a PEP-CTERM (Pro-Glu-Pro carboxy-terminal) domain. This domain is predicted to target proteins for export into the exopolysaccharide layer . The RR gene itself is encoded on the native pDV1 plasmid in a ten-gene operon that appears related to exopolysaccharide synthesis. The five predicted PEP-CTERM targets  for this RR were enriched in the DAP-chip assay along with several other hits (Figure 5; Table S5 in Additional file 1). The pDV1- D. vulgaris strain (lacking the native plasmid) produces three-fold less biofilm than the wild type, and the wild-type biofilm also contains less carbohydrate and more protein filaments . It is possible that DVUA0057 is involved in biofilm formation and that the PEP-CTERM proteins are an integral part of the biofilm. Expression correlations show that most of the PEP-CTERM proteins are positively correlated with each other and DVUA0057 (Figure S3 in Additional file 3). The RRs DVU0804 and DVU0110 have a few overlapping targets with DVUA0057 (Figure 5), and may also be involved in similar functions.
A three-component system for modulating general stress responses
DVU3305 may be part of a three-component system that includes a HK, DVU3304, and another RR, DVU3303, that has an unusual domain structure containing a lon protease. DVU3305 regulates its own operon as well as its upstream genes DVU3302-3298, which encode membrane and hypothetical proteins, some of which have UspA type domains (Figure 6; Table S5 in Additional file 1). The two operons have good expression correlation to each other (Figure S3 in Additional file 3). DVU3303-3305 and/or DVU3298-3302 are upregulated during various stresses, such as high pH , heat shock , and stationary phase . DVU3303 may require phosphorylation to be protease-active, which in turn may be part of the general stress response. The other targets for this RR are the operon DVU2822-2825, which encodes a putative dicarboxylate transporter and pyruvate formate lyase genes, and the high affinity potassium transporter DVU3334-3339 (see below).
Potassium uptake genes are controlled by multiple RRs
DVU3334 regulates its own operon, encoding the high affinity potassium uptake genes (kdpFABC). The HK in this operon appears to be split into two ORFs, with DVU3336 having the K+ sensor domain and an UspA domain, and DVU3335 having the HK dimerization and phosphoacceptor domain. Interestingly, two other RRs also target the Kdp operon - DVU0596 and DVU3305 - suggesting that potassium uptake is a response to other stresses as well. The putative binding sites for DVU3305 and DVU3334 upstream of the Kdp genes also overlap (Table S9 in Additional file 1).
DVU0596 has a LytR type DBD, and its other targets are two copies of a putative carbon starvation (cstA) gene that lie downstream of this RR-HK operon. In E. coli the cstA gene is upregulated upon exhaustion of carbon source in the medium, and appears to be a peptide permease. Activated DVU0596 may also function by regulating the putative cstA. Presumably due to the regulation of the DVU3334-3339 Kdp operon by multiple RRs, the genes in this operon do not demonstrate simple correlation patterns with the co-cistronic RR DVU3334 (Figure S3 in Additional file 3).
Regulation of acetyl-CoA levels
A second LytR-type RR, DVU0749, has four genes/operons as its targets (Figure 4). DVU2969 is annotated as an acetyl-CoA synthetase similar to DVU0748, which is in the same operon as the RR itself. DVU2970 has acyl-CoA synthetase and acetyltransferase domains. These seem to suggest that this RR is involved in maintaining acetyl-CoA levels in the cell. Based on the gene annotations for DVU0443-0448, they may constitute a cAMP-dependent membrane transporter or signal transduction protein. Regulon prediction, based on the gene neighbor method , also associates this operon with DVU2969. Additionally, DVU0749 has an atypical receiver domain, with the phosphorylatable aspartic acid replaced by a glycine. Since it also lacks a proximally encoded HK, it is likely that this RR is not activated by phosphorylation.
Two component system involved in energy metabolism regulation
DVU2394 regulates the DVU2405-2397 operon, which encodes the alcohol dehydrogenase and heterodisulfide reductase genes (Figure 6). The adh gene is necessary for growth with ethanol as electron donor, but was also highly expressed during growth on lactate, pyruvate, formate or hydrogen as the electron donor . The operon encoding the RR and DVU2405-2397 are positively correlated, but they appear to be anti-correlated to the other targets DVU3298-3305, which are part of the general stress response (Figure S3 in Additional file 1).
Response to nitrite stress
Nitrates and nitrites are known to inhibit sulfate reduction, and thus pose a stress for sulfate-reducing bacteria. DVU0621 targets its upstream genes DVU0624-0625, which encode a nitrite reductase and are highly upregulated in the presence of nitrite [39, 40]. The nitrite reductase protein has been crystallized and biochemically shown to reduce nitrite as well as sulfite [41, 42]. The tiling array data indicate DVU0624-0625 to be highly expressed in the absence of nitrite (Figure S1 in Additional file 3) , and therefore suggest roles in addition to responding to nitrite stress. DVU0621 also targets the lactate/pyruvate oxidation genes DVU3025-3033, and another hypothetical gene, DVU3384.
The nitrogen regulator is the most conserved among related species
Despite the identification of a target gene via EMSA, no DAP-chip data could be acquired for RR DVU3220, since conditions for qPCR enrichment of the target could not be found. However, DVU3220, along with DVU1083, are the only two RRs that are conserved in all sequenced Desulfovibrio and other related genomes as well (Figure 8). DVU1083 is the phosphate regulator, and it seems likely that DVU3220 is the nitrogen regulator. The identification (by EMSA; Figure 3) of DVU1231-1234, the operon encoding the ammonium transporter, the nitrogen regulatory protein PII, and the P-II uridylyltransferase, as its target supports this hypothesis.
Regulation through sRNAs
The RR DVU0679 has as its target a sRNA that lies downstream of the RR gene (Figure 4) that has been annotated as Dv_sRNA2 (Bender, personal correspondance). Tiling array data  confirm the expression of this sRNA during normal growth (Figure S1 in Additional file 3), but its function is unknown. The validated binding site motif is present in other unique DAP-chip hits, which are not in upstream regions but may be physiologically relevant. It may be that there are as yet undiscovered non-coding RNAs present near these binding sites, or that binding within coding regions or within an operon presents additional ways of modulating expression of the target genes .
Functional validation using reporter assays
For several of the RRs, the corresponding target promoters were tested using in vivo transcriptional reporter analysis. Examination of two component system RRs in D. vulgaris would require the knowledge of the activating signals. To bypass this requirement, we sought to examine the binding specificity in the heterologous host E. coli. The assay utilized two compatible plasmids: pETDEST42, which expressed the RR from an IPTG-inducible promoter; and pBbS2K-RFP, which expressed the red fluorescent protein (RFP) gene under the control of the target D. vulgaris promoter (Materials and methods; Figure S4 in Additional file 3). Expression of RRs in the reporter strains was confirmed using anti-His immunoblots (Figure S4 in Additional file 3).
Binding site motifs are conserved across related species
The presence of orthologous RRs with conserved binding motifs and target genes in other sequenced microbes is a strong indicator of similar function. Orthologs of D. vulgaris RRs could be tracked in several other genomes (Figure 8). We screened 23 sequenced genomes that contained orthologs of any of the D. vulgaris Hildenborough RRs for the presence of conserved binding motifs and target genes. These included six Desulfovibrio genomes, specifically, the closely related D. vulgaris Miyazaki, D. desulfuricans G20, D. salexigens, the magnetotactic D. magneticus, the rumen isolate D. desulfuricans 27774 and the human isolate D. piger. Two RRs that are conserved in all these species were DVU1083, the phosphate regulator, and DVU3220, a possible nitrogen regulator. DVU3220 is also conserved in some of the other sulfate-reducing bacteria we examined, namely Desulfomicrobium baculatum, Desulfohalobium retbaense, Desulfobacterium autotrophicum, Desulfatibacillum alkenivorans, and Desulfotalea psychrophila (Figure 8).
DVU3023, the lactate responsive RR in D. vulgaris, is conserved in all the environmental Desulfovibrio isolates, as well as D. baculatum, D. retbaense and D. psychrophila, where the binding site motif is also conserved upstream of the central lactate-pyruvate oxidation genes. Similarly, the alcohol dehydrogenase regulator DVU2394 is also conserved in many of the sulfate reducers, although the validated motif for D. vulgaris Hildenborough was only found in D. vulgaris Miyazaki (Figure 8).
The flagella regulator DVU1063 is conserved among the Desulfovibrio spp., being absent in only the non-motile D. piger. It is also present in the related pathogen Lawsonia intracellularis and in D. baculatum, and the binding site motif is also conserved upstream of several flagella and motility related genes in these species (Table S9 in Additional file 1). The pili assembly regulator DVU2114 and its binding site motif upstream of pilin genes are conserved in only a few species - D. vulgaris Miyazaki, D. desulfuricans G20 and Syntrophobacter fumaroxidans (also a sulfate reducer with propionate as electron donor). A DVU2114 ortholog is also present in D. magneticus, but this genome lacks the target pilin genes and a similar binding site motif was present upstream of a different gene (Table S9 in Additional file 1).
Orthologs of the paralogous RRs DVU0946 and DVU0539 and their binding site motifs are conserved in all the environmental Desulfovibrio spp. as well as other sulfate reducers such as D. baculatum, D. retbaense and D. alkenivorans (Figure 8). Interestingly, D. vulgaris Miyazaki and D. retbaense also have paralogous copies for this RR. The orthologs in D. desulfuricans G20 and D. vulgaris Miyazaki present a case where the motif occurs upstream of additional genes different from those identified in D. vulgaris Hildenborough. These genomes have the motif upstream of orthologs of the target DVU0943, but they also have it upstream of a putative split soret cytochrome c precursor gene that is not a target in D. vulgaris Hildenborough (Table S9 in Additional file 1). These additional genes may be true targets for the respective RRs, and our genome scans can be a valuable tool for identifying potential targets.
DVU2934, which targets lipid A biosynthesis, has several orthologs, including those found in non-sulfate reducers such as Geobacter lovleyi and Acidobacterium capsulatum. However, the validated motif was present upstream of the target lpxC gene only in D. vulgaris Miyazaki and D. desulfuricans G20 (Figure 8; Table S9 in Additional file 1). It is possible that the other RRs may have evolved to have different functional roles. Alternately, the RR may have diverged to recognize a different motif upstream of the same target. Suggestive of this latter occurrence, in S. fumaroxidans, the lpxC gene is directly upstream of the ortholog RR operon, suggesting that its regulation is likely to be linked to the orthologous RR despite the absence of the motif.
Orthologs for the general stress responsive RR DVU3305 are present in some of the Desulfovibrio species, D. baculatum and D. retbaense, and also in less related species such as Thermodesulfovibrio yellowstonii, Thioalkalivibrio spp. and Dechloromonas aromatica. In most of these genomes, the RR is also associated with the lon-protease RR ortholog of DVU3303, and the binding site is also conserved, indicating that their functions are likely to be related (Figure 8; Table S9 in Additional file 1).
The potassium responsive RR DVU3334 has orthologs only in D. vulgaris Miyazaki and S. fumaroxidans, where the motif is also conserved (Figure 8). Interestingly, other Desulfovibrio species and sulfate reducers do not have the high affinity potassium uptake Kdp genes.
The sRNA regulating RR DVU0679 has orthologs in some of the Desulfovibrio species and in another sulfate reducer, Desulfotomaculum acetoxidans, and in Syymbiobacterium thermophilum. Genome scans show that a target sRNA with a conserved binding site may lie downstream of the RR orthologs in D. vulgaris Miyazaki and D. desulfuricans G20 (Table S9 in Additional file 1), suggesting that sRNA regulation is a function of these orthologous RRs as well. Genome scans of D. acetoxidans and S. thermophilum also showed that the highest scoring hits to the motif were not upstream of coding regions (Table S9 in Additional file 1).
Other RRs such as DVU1419 and DVU3381, which target hypothetical genes and therefore have unknown functions, are conserved along with their binding sites in closely related species. Orthologs of RRs without binding site motif predictions, specifically DVU0621, DVU0653, DVU0804, DVU0744, DVU2675, and DVU2577, are limited to a few of the Desulfovibrio species (Figure 8).
Conservation of RRs in other Deltaproteobacteria, such as Geobacter species and the myxobacteria, are shown in Figure 8. Three RRs do not have any other Desulfovibrio orthologs, although orthologs may be present in more distant species. These include DVU1156 and the pDV1 encoded DVUA0137, with no functional predictions from our study, and DVU2588, which may be involved in lactate utilization. The binding site motif for DVU2588 is also conserved in S. fumaroxidans and A. capsulatus but is present upstream of genes different from those in D. vulgaris (Figure 8; Table S9 in Additional file 1).
Prior to our study, very little was known about the two component regulatory network in sulfate-reducing bacteria. Here we provide a fairly comprehensive map of genes that are transcriptionally regulated by the majority of the two component systems in this model sulfate-reducing organism (Figure 7). Our results include 200 target genes for 24 response regulators and provide strong predictions for the corresponding two component systems that include the regulation of cell motility (flagella and pili), exopolysaccharide production, energy metabolism (lactate utilization, alcohol dehydrogenase regulation, acetyl-CoA levels), lipid A synthesis, nitrogen and phosphate metabolism, and in the responses to stresses such as low potassium, nitrite and carbon starvation. Functions such as lactate utilization and potassium uptake genes are regulated by multiple RRs and appear central to the stress response in D. vulgaris. In addition, the experimentally confirmed binding motifs for several of these RRs could also be used to assess gene targets and function of the orthologs of these regulators in related bacteria. With the exception of a few RRs (for example, DVU1083 (PhoB) and DVU2934), most of the DBD-containing RRs encoded in D. vulgaris Hildenborough appear to be unique to the sequenced Desulfovibrio species and closely related sulfate reducers. This suggests that responses modulated by most two component systems in these sulfate reducing bacteria are unique to the ecological niches they occupy. In addition, the significant numbers of hypothetical proteins and genes with unknown functions amongst the regulated candidates indicate that there remains a lot to be learnt about the environmental stresses faced in these ecosystems. Deeper knowledge of stress response and regulation in these environments is required for robust bioremediation approaches and a better understanding of biogeochemical processes mediated by these bacteria.
Materials and methods
Cloning of response regulator genes in E. coli
D. vulgaris Hildenborough was grown on LS4D medium , and the cell pellet was used to purify genomic DNA with the Qiagen Genomic DNA buffer kit, and the Qiagen midi column (tip 100/G). The response regulator gene was PCR amplified and cloned into the entry vector pENTR™/SD/D-TOPO (Life Technologies, Grand Island, NY, USA) with forward primers carrying the 5' sequence CACC. A list of primers used for cloning the genes is in Table S11 in Additional file 1. The entry clone was transformed into Invitrogen's OneShot Top10 chemically competent cells, and selected on LB-Kanamycin plates. Sequencing was used to confirm the presence of the insert. The expression construct was generated by performing an LR recombination reaction (Gateway LR Clonase II, Life Technologies, Grand Island, NY, USA) with the destination vector pETDEST42 (Invitrogen) such that the RR gene is expressed with a carboxy-terminal V5-epitope and a 6X His-tag. The final construct was transformed into E. coli BL21 Star (DE3) OneShot chemically competent cells (Invitrogen), and selected on LB-Amp plates. Sequencing was used to confirm the presence of the insert.
Protein expression and purification
The E. coli BL21 (DE3) pETDEST42-RR strains were grown on LB-carbenicillin at 37°C. At mid-log phase, the cells were induced with 0.5 mM IPTG (isopropyl-beta-D-thiogalactopyranoside), and then grown at RT for 24 hours. The cells were pelleted, and resuspended in HisTrapFF wash/binding buffer (40 mM imidazole, 500 mM NaCl, and 20 mM sodium phosphate, pH 7.4). Lysozyme (1 mg/ml) and 1X Novagen's benzonase nuclease were added to the cell suspension. The cells were lysed using a French Press at 4°C, and the cell lysate was clarified by spinning at 10, 000 rpm at 4°C. To check for over-expression, a sample was run on 4 to 12% Bis-Tris gel, the gel was transferred onto a nitrocellulose membrane, and a western blot was performed with mouse anti-6X His-tag antibodies. Prior to purification using the AKTA Explorer, the clarified lysate was filtered through a 0.45 μm syringe filter. The lysate was loaded onto a 1 ml HisTrapFF column (GE Healthcare, Piscataway, NJ, USA) that had been equilibrated with HisTrapFF wash/binding buffer (10 ml). The column was washed with 20 ml of the wash buffer, and then eluted with a linear gradient of 0 to 100% elution buffer (500 mM imidazole, 500 mM NaCl, and 20 mM sodium phosphate, pH 7.4). The pooled protein fractions were stripped of imidazole using a HiPrep 26/10 Desalting column (GE Healthcare) and washing with a desalting buffer (500 mM NaCl, 20 mM sodium phosphate, pH 7.4). The protein preps were concentrated using a high molecular weight cutoff spin filter, and glycerol (50%) and DTT (0.1 mM) were added to the preparations for storage. Examples of purified proteins are shown in Figure S6 in Additional file 3.
Preparing substrates for EMSAs
Oligonucleotides (oligos; unlabeled and 5'-biotinylated) were ordered from IDT, San Diego, CA, USA. Full length (200 to 400 bp) upstream regions of target genes (Figure 3) were prepared by PCR amplification from D. vulgaris Hildenborough genomic DNA using one unlabeled oligo and one 5'-biotin-labeled oligo. A list of oligos is given in Table S2 in Additional file 1. The biotin-labeled substrates were then gel-purified from agarose gels using Qiaquick gel extraction kits (Qiagen Inc, Valencia, CA, USA). Substrates for validating binding site motifs (Figure 9) were prepared by annealing oligos carrying the motif to be tested and approximately 10 bp on either side (Table S6 in Additional file 1). The top strand biotinylated oligo was mixed with a slight excess of unlabeled bottom strand oligo in 10 mM Tris HCl, pH 8.0, 1 mM EDTA, and 50 mM NaCl, and heated to 95°C for 5 minutes followed by slow cooling to 25°C. The annealed substrate was then diluted to 1 pmol/μl.
Electrophoretic mobility shift assay
EMSAs were performed using the Pierce Lightshift chemiluminescent kit (Thermo Fisher Scientific, Rockford, IL, USA). The response regulator to be tested was mixed with 50 to 100 fmol of biotinylated DNA in 10 mM Tris HCl pH 7.5, 1 mM DTT, 50 mM KCl, 5 mM MgCl2 and 25% glycerol. Poly dI-dC (1 μg/ml) was added as a non-specific competitor. The reactions were incubated at 25 to 30°C for 20 minutes, and the reactions were loaded on a pre-cast mini 6% polyacrylamide-0.5X TBE gel (Life Technologies), and run at 100 V at room temperature. The gels were transferred to a nylon membrane by semi-dry blotting (BioRad, Hercules, CA), the nylon membranes were UV cross-linked for 3 minutes, and the blot was developed using the Pierce Chemiluminescent Nucleic Acid Detection kit (Thermo Fisher Scientific). The blot was scanned using the Typhoon 8600 Imager (Molecular Dynamics/Amersham Pharmacia (Piscataway, NJ, USA)).
Binding reactions with genomic DNA and DNA affinity purification
D. vulgaris Hildenborough genomic DNA was prepared using a Qiagen genomic tip and genomic buffer set. The DNA was sheared by sonication to an average length of 500 bp. Binding reactions (100 μl) were set up with 2 to 3 μg of sheared genomic DNA, and the appropriate amount of purified response regulator, in 10 mM Tris HCl, pH 7.5, 1 mM DTT, 5 mM MgCl2, 50 mM acetyl phosphate and 25% glycerol. The amount of protein used in the binding reactions was determined based on activity of protein in EMSAs (Table S3 in Additional file 1). Acetyl phosphate was used as a generic method to activate the RRs based on protocols described in the literature . The reactions were incubated at 25°C in a thermal cycler for 30 minutes; 10 μl of the reaction was saved as input DNA. The rest was allowed to bind to 30 μl of Ni-NTA resin that had been washed in the binding/wash buffer (10 mM Tris-HCl, pH 7.5, 5 mM MgCl2, 50 mM KCl, 25% glycerol). After 30 minutes of binding, the tubes were spun down to remove the unbound DNA. The resin was washed three times in 100 μl of the wash buffer, and then the bound DNA was eluted with 35 μl of elution buffer (20 mM sodium phosphate buffer pH 8, 500 mM NaCl, and 500 mM imidazole); 35 μl of this elution buffer was also added to the input DNA. The input and the enriched DNA fractions were cleaned up using Qiaquick PCR purification columns (Qiagen Inc).
Whole genome amplification
The input and enriched DNA samples (10 μl; after clean up) were subjected to whole-genome amplification using Sigma WGA2 kit (Sigma-Aldrich, St. Louis, MO, USA)). Since the starting material was sheared genomic DNA, the first fragmentation heat step was omitted and the number of amplification cycles was increased to 20 as per the manufacturer's suggestions. The amplified samples were cleaned up using Qiagen spin columns and quantified using the nanodrop.
qPCR was performed on the whole genome amplified input and enriched DNA samples using the PerfectA Sybr Green Mix with ROX (Quanta Biosciences, Gaithursburg, MD, USA). The DNA samples were diluted to a concentration of 2 ng/μl, and 5 μl of each sample was used as the template. The primers used for qPCR of target upstream regions are listed in Table S2 in Additional file 1. Each reaction was done in triplicate. Delta CT was calculated as the difference in the CT values of the input and enriched samples (ΔCT = CT (input) - CT (enriched). Fold enrichment was calculated as 2ΔCt. If the target was found to be at higher amounts in the enriched sample, then the samples were Cy3/Cy5 labeled and hybridized to the chip. If no enrichment of the target was observed, then the DAP-whole-genome amplification steps were repeated under different conditions (usually by varying the protein amount) until enrichment was obtained. For each RR tested, the upstream region of a randomly selected gene was also tested to ensure that non-specific gene targets were not enriched. For most RRs, the negative control used was the upstream region of DVU0013. The exceptions were DVU3234 used for RR DVU1083, DVU0599 used for RR DVU1063, and DVU1083 used for RR DVU0946.
DNA labeling and hybridization
Cy3- and Cy5-labeled 9-mers were obtained from Trilink Biotechnologies (San Diego, CA, USA). One microgram each of the input and enriched DNA was labeled with Cy3 and Cy5, respectively, using Klenow polymerase (3'-5' exo-, 50, 000 U/ml, New England Biolabs, Ipswich, MA, USA) at 37°C for 2 hours in the dark. The reactions were stopped by adding EDTA to a final concentration of 50 mM. The DNA was precipitated using NaCl and isopropanol, and the pellet was rinsed in 80% ethanol, and air-dried. The pellets were resuspended in 25 μl water, and quantified using the nanodrop. Six micrograms each of the Cy3- and Cy5-labeled DNA were pooled together, and vacuum dried. The pellet was resuspended in 5 μl water, and hybridized overnight at 42°C to a custom D. vulgaris Hildenborough 385K tiling array (Roche-Nimblegen, Madison, WI) as per the manufacturer's instructions. The tiling array used 50-mer probes with a 46-bp overlap in the intergenic regions, and a 20-bp overlap in the coding regions. The slides were washed, and scanned on a GenePix 4200A scanner. Each slide was scanned at wavelengths of 532 and 635 nm, and the images were saved separately.
Data analysis to generate binding site peaks
Nimblescan software v.2.4 (Roche-Nimblegen) was used to grid the 532 and 635 nm images, and to generate the pair files that contain signal intensities for each probe on the array. Scaled log2R ratios were then computed for the probes. Nimblescan then used a sliding window of 500 bp to search for peaks that have more than four probes above a cutoff value. The cutoff value is a percentage of the hypothetical maximum log2R ratio (mean + 6 standard deviation) that starts at the most stringent value of 90% and decreases in steps of 1 to 15%. The ratio data are randomized 20 times to evaluate the possibility of false positives. As the method to estimate the FDR scores for the peaks, we chose to optimize for strong peaks of varied width. The peak file generated by Nimblescan gives the location of the peak (genomic or plasmid), the start and end loci for the peak, a score that is the log2R ratio for the fourth highest probe in the peak, the cutoff percentage (cutoff_p) at which the peak was identified, and the FDR score for the peak. We used MicrobesOnline  to map the peak loci to the upstream regions of the gene targets. Nimblescan guidelines state that FDR scores of 0 to 0.05 are indicative of a high-confidence binding site, while peaks with scores of up to 0.2 may also be considered as binding sites. Since our DAP-chip data sets had a large number of peaks with FDR scores less than 0.2, a combination of high log2R scores, high cutoff_p values, and low FDR scores was used to identify the most likely targets for each RR. The list of targets presented in Table S5 in Additional file 1 and Figures 4, 5 and 6 is a conservative estimate of the most likely targets. Peaks falling within coding regions, or those that appear in multiple DAP-chip sets, or those with lower scores or higher FDR values were not considered unless they were also found to carry a binding site motif, or appeared to be functionally related to the other top peaks. For DVU1083, where the positive target did not appear in the top hits, DAP-chip was conducted in triplicate and the common hits from the three sets were pooled together to generate a target list. The array data have been deposited in NCBI's Gene Expression Omnibus [GEO:GSE25163] .
Binding site motif predictions for the RRs
MEME was used (through MicrobesOnline  and The Meme Suite ) to predict motifs in the upstream regions. For each RR, motif predictions were run on the upstream regions of the positive target orthologs, and for those RRs that had more than one target MEME was also run on the upstream regions of the target hits obtained from the DAP-chips. Genomes of interest were scanned for the resulting motif using Perl scripts based on the methods as described in Rhodius et al. .
Sigma54 promoter predictions
The rpoN matrix was downloaded from Promscan . Promscan generated an rpoN matrix based on the σ54 promoter sequences in more than 80 bacterial species . We used the Perl script to determine the σ54 promoters of D. vulgaris. Additionally, to test specific upstream regions, we ran the sequence on the Promscan website directly. A list of σ54-promoters that are regulated by RRs is given in Table S10 in Additional file 1.
Construction of pBbS2K-promoter-RFP plasmids
Promoter regions (300 bp upstream) of selected D. vulgaris genes were PCR amplified using primers (Table S12 in Additional file 1) that have a 20-bp overlap with the target vector. The kanamycin-resistant vector pBbS2K-RFP contains the RFP gene under the control of the Ptet promoter and tetR gene. The RFP gene is repressed by TetR, and can be induced using anhydrotetracycline (aTC). The vector pBbS2K-RFP was PCR amplified using promoter-specific primers that have a 20-bp overlap with the promoter of interest (Table S12 in Additional file 1) and the promoters were cloned in by the Gibson method  downstream of the Ptet and upstream of rfp such that they replace the native ribosome binding site. The Gibson reaction was used to transform E. coli Top10 electrocompetent cells and transformants were selected on LB-Kan-aTC plates and red-white screening was used to pick positive colonies. Sequencing was used to confirm the presence of the promoter insert in the pBbS2K-promoter-RFP constructs.
RFP reporter assay in E. coli
E. coli BL21 (DE3) electrocompetent cells were transformed with equal amounts of pBbS2K-promoter-RFP (pSC101 ori, KanR) and pETDEST42-RR (pBR322 ori, AmpR) plasmids (Figure S4 in Additional file 3), and the transformants were selected on LB-Kan-carbenicillin plates. The colonies were adapted to M9 minimal media as described previously . Overnight cultures grown in M9-Kan-Carb were subcultured 1:100 in fresh media and grown at 37°C for 2 hours. The cultures were induced with 50 μM IPTG and 1 ml of each culture (in triplicate) was transferred to a 24-well plate and the plate was incubated with shaking at 25°C in a Tecan Infinite 200 instrument (Tecan, Männedorf, Switzerland) for 24 hours. Cell growth was monitored by measuring absorbance at 590 nm. RFP fluorescence was measured at 535 nm (excitation) and 610 nm (emission). Relative fluorescence units (RFU) were calculated as the ratio of RFP fluorescence to cell growth.
DNA Affinity Purified-chip
electrophoretic mobility shift assay
false discovery rate
open reading frame
polymerase chain reaction
red fluorescent protein
We thank Mario Ouellet (LBNL) and Dr Adam Deutschbauer (LBNL) for helpful discussions about the DAP-chip strategy, Dr Ee-Been Goh (LBNL) and Virgil Rhodius (UCSF) for providing Perl scripts to search for motifs in the various bacterial genomes, Dr Nathan Hillson (LBNL) for help with modifying Perl scripts and Katherine Huang (Broad Institute) for help with generating Figure 7 and Dr Kelly Bender (Southern Illinois University) for sharing unpublished data. We thank Prof. Eric Alm (MIT) and Dr Igor Zwir (Washington University) for reviewing the manuscript. This research is funded by ENIGMA, a Scientific Focus Area Program supported by the US Department of Energy, Office of Science, Office of Biological and Environmental Research, Genomics:GTL Foundational Science through contract DE-AC02-05CH11231 between Lawrence Berkeley National Laboratory and the US Department of Energy.
- Gao R, Stock AM: Biological insights from structures of two-component proteins. Annu Rev Microbiol. 2009, 63: 133-154.PubMedPubMed CentralView ArticleGoogle Scholar
- Galperin MY: Diversity of structure and function of response regulator output domains. Curr Opin Microbiol. 2010, 13: 150-159.PubMedPubMed CentralView ArticleGoogle Scholar
- Laub MT, Goulian M: Specificity in two-component signal transduction pathways. Annu Rev Genet. 2007, 41: 121-145.PubMedView ArticleGoogle Scholar
- Alm E, Huang K, Arkin A: The evolution of two-component systems in bacteria reveals different strategies for niche adaptation. PLoS Comput Biol. 2006, 2: e143-PubMedPubMed CentralView ArticleGoogle Scholar
- Heidelberg JF, Seshadri R, Haveman SA, Hemme CL, Paulsen IT, Kolonay JF, Eisen JA, Ward N, Methe B, Brinkac LM, Daugherty SC, Deboy RT, Dodson RJ, Durkin AS, Madupu R, Nelson WC, Sullivan SA, Fouts D, Haft DH, Selengut J, Peterson JD, Davidsen TM, Zafar N, Zhou L, Radune D, Dimitrov G, Hance M, Tran K, Khouri H, Gill J, et al: The genome sequence of the anaerobic, sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough. Nat Biotechnol. 2004, 22: 554-559.PubMedView ArticleGoogle Scholar
- Han JS, Park JY, Lee YS, Thony B, Hwang DS: PhoB-dependent transcriptional activation of the iciA gene during starvation for phosphate in Escherichia coli. Mol Gen Genet. 1999, 262: 448-452.PubMedView ArticleGoogle Scholar
- Oshima T, Aiba H, Masuda Y, Kanaya S, Sugiura M, Wanner BL, Mori H, Mizuno T: Transcriptome analysis of all two-component regulatory system mutants of Escherichia coli K-12. Mol Microbiol. 2002, 46: 281-291.PubMedView ArticleGoogle Scholar
- Salgado H, Santos-Zavaleta A, Gama-Castro S, Peralta-Gil M, Penaloza-Spinola MI, Martinez-Antonio A, Karp PD, Collado-Vides J: The comprehensive updated regulatory network of Escherichia coli K-12. BMC Bioinformatics. 2006, 7: 5-PubMedPubMed CentralView ArticleGoogle Scholar
- Skerker JM, Prasol MS, Perchuk BS, Biondi EG, Laub MT: Two-component signal transduction pathways regulating growth and cell cycle progression in a bacterium: a system-level analysis. PLoS Biol. 2005, 3: e334-PubMedPubMed CentralView ArticleGoogle Scholar
- Kobayashi K, Ogura M, Yamaguchi H, Yoshida K, Ogasawara N, Tanaka T, Fujita Y: Comprehensive DNA microarray analysis of Bacillus subtilis two-component regulatory systems. J Bacteriol. 2001, 183: 7365-7370.PubMedPubMed CentralView ArticleGoogle Scholar
- Makita Y, Nakao M, Ogasawara N, Nakai K: DBTBS: database of transcriptional regulation in Bacillus subtilis and its contribution to comparative genomics. Nucleic Acids Res. 2004, 32: D75-77.PubMedPubMed CentralView ArticleGoogle Scholar
- Covert MW, Knight EM, Reed JL, Herrgard MJ, Palsson BO: Integrating high-throughput and computational data elucidates bacterial networks. Nature. 2004, 429: 92-96.PubMedView ArticleGoogle Scholar
- Liu X, Noll DM, Lieb JD, Clarke ND: DIP-chip: rapid and accurate determination of DNA-binding specificity. Genome Res. 2005, 15: 421-427.PubMedPubMed CentralView ArticleGoogle Scholar
- Mascher T, Zahner D, Merai M, Balmelle N, de Saizieu AB, Hakenbeck R: The Streptococcus pneumoniae cia regulon: CiaR target sites and transcription profile analysis. J Bacteriol. 2003, 185: 60-70.PubMedPubMed CentralView ArticleGoogle Scholar
- Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E, Volkert TL, Wilson CJ, Bell SP, Young RA: Genome-wide location and function of DNA binding proteins. Science. 2000, 290: 2306-2309.PubMedView ArticleGoogle Scholar
- Dehal PS, Joachimiak MP, Price MN, Bates JT, Baumohl JK, Chivian D, Friedland GD, Huang KH, Keller K, Novichkov PS, Dubchak IL, Alm EJ, Arkin AP: MicrobesOnline: an integrated portal for comparative and functional genomics. Nucleic Acids Res. 2010, 38: D396-400.PubMedPubMed CentralView ArticleGoogle Scholar
- Price MN, Huang KH, Alm EJ, Arkin AP: A novel method for accurate operon predictions in all sequenced prokaryotes. Nucleic Acids Res. 2005, 33: 880-892.PubMedPubMed CentralView ArticleGoogle Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- Price MN, Deutschbauer A, Kuehl JV, Liu H, Witkowska HE, Arkin AP: Evidence-based annotation of transcripts and proteins in the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough. J Bacteriol. 2011, 193: 5716-5727.PubMedPubMed CentralView ArticleGoogle Scholar
- McCleary WR, Stock JB: Acetyl phosphate and the activation of two-component response regulators. J Biol Chem. 1994, 269: 31567-31572.PubMedGoogle Scholar
- Sorensen PG, Lutkenhaus J, Young K, Eveland SS, Anderson MS, Raetz CR: Regulation of UDP-3-O-[R-3-hydroxymyristoyl]-N-acetylglucosamine deacetylase in Escherichia coli. The second enzymatic step of lipid a biosynthesis. J Biol Chem. 1996, 271: 25898-25905.PubMedView ArticleGoogle Scholar
- Fuhrer F, Langklotz S, Narberhaus F: The C-terminal end of LpxC is required for degradation by the FtsH protease. Mol Microbiol. 2006, 59: 1025-1036.PubMedView ArticleGoogle Scholar
- Su Z, Olman V, Xu Y: Computational prediction of Pho regulons in cyanobacteria. BMC Genomics. 2007, 8: 156-PubMedPubMed CentralView ArticleGoogle Scholar
- Pratt JT, McDonough E, Camilli A: PhoB regulates motility, biofilms, and cyclic di-GMP in Vibrio cholerae. J Bacteriol. 2009, 191: 6632-6642.PubMedPubMed CentralView ArticleGoogle Scholar
- Baek JH, Lee SY: Novel gene members in the Pho regulon of Escherichia coli. FEMS Microbiol Lett. 2006, 264: 104-109.PubMedView ArticleGoogle Scholar
- Yuan ZC, Zaheer R, Morton R, Finan TM: Genome prediction of PhoB regulated promoters in Sinorhizobium meliloti and twelve proteobacteria. Nucleic Acids Res. 2006, 34: 2686-2697.PubMedPubMed CentralView ArticleGoogle Scholar
- Bernard CS, Bordi C, Termine E, Filloux A, de Bentzmann S: Organization and PprB-dependent control of the Pseudomonas aeruginosa tad Locus, involved in Flp pilus biology. J Bacteriol. 2009, 191: 1961-1973.PubMedPubMed CentralView ArticleGoogle Scholar
- Juarez K, Kim BC, Nevin K, Olvera L, Reguera G, Lovley DR, Methe BA: PilR, a transcriptional regulator for pilin and other genes required for Fe(III) reduction in Geobacter sulfurreducens. J Mol Microbiol Biotechnol. 2009, 16: 146-158.PubMedView ArticleGoogle Scholar
- Wu SS, Kaiser D: Regulation of expression of the pilA gene in Myxococcus xanthus. J Bacteriol. 1997, 179: 7748-7758.PubMedPubMed CentralGoogle Scholar
- Jin S, Ishimoto KS, Lory S: PilR, a transcriptional regulator of piliation in Pseudomonas aeruginosa, binds to a cis-acting sequence upstream of the pilin gene promoter. Mol Microbiol. 1994, 14: 1049-1057.PubMedView ArticleGoogle Scholar
- Skerker JM, Shapiro L: Identification and cell cycle control of a novel pilus system in Caulobacter crescentus. EMBO J. 2000, 19: 3223-3234.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu J, Benson AK, Newton A: Global regulation of a sigma 54-dependent flagellar gene family in Caulobacter crescentus by the transcriptional activator FlbD. J Bacteriol. 1995, 177: 3241-3250.PubMedPubMed CentralGoogle Scholar
- Haft DH, Paulsen IT, Ward N, Selengut JD: Exopolysaccharide-associated protein sorting in environmental organisms: the PEP-CTERM/EpsH system. Application of a novel phylogenetic profiling heuristic. BMC Biol. 2006, 4: 29-PubMedPubMed CentralView ArticleGoogle Scholar
- Clark ME, Edelmann RE, Duley ML, Wall JD, Fields MW: Biofilm formation in Desulfovibrio vulgaris Hildenborough is dependent upon protein filaments. Environ Microbiol. 2007, 9: 2844-2854.PubMedView ArticleGoogle Scholar
- Stolyar S, He Q, Joachimiak MP, He Z, Yang ZK, Borglin SE, Joyner DC, Huang K, Alm E, Hazen TC, Zhou J, Wall JD, Arkin AP, Stahl DA: Response of Desulfovibrio vulgaris to alkaline stress. J Bacteriol. 2007, 189: 8944-8952.PubMedPubMed CentralView ArticleGoogle Scholar
- Chhabra SR, He Q, Huang KH, Gaucher SP, Alm EJ, He Z, Hadi MZ, Hazen TC, Wall JD, Zhou J, Arkin AP, Singh AK: Global analysis of heat shock response in Desulfovibrio vulgaris Hildenborough. J Bacteriol. 2006, 188: 1817-1828.PubMedPubMed CentralView ArticleGoogle Scholar
- Clark ME, He Q, He Z, Huang KH, Alm EJ, Wan XF, Hazen TC, Arkin AP, Wall JD, Zhou JZ, Fields MW: Temporal transcriptomic analysis as Desulfovibrio vulgaris Hildenborough transitions into stationary phase during electron donor depletion. Appl Environ Microbiol. 2006, 72: 5578-5588.PubMedPubMed CentralView ArticleGoogle Scholar
- Haveman SA, Brunelle V, Voordouw JK, Voordouw G, Heidelberg JF, Rabus R: Gene expression analysis of energy metabolism mutants of Desulfovibrio vulgaris Hildenborough indicates an important role for alcohol dehydrogenase. J Bacteriol. 2003, 185: 4345-4353.PubMedPubMed CentralView ArticleGoogle Scholar
- He Q, Huang KH, He Z, Alm EJ, Fields MW, Hazen TC, Arkin AP, Wall JD, Zhou J: Energetic consequences of nitrite stress in Desulfovibrio vulgaris Hildenborough, inferred from global transcriptional analysis. Appl Environ Microbiol. 2006, 72: 4370-4381.PubMedPubMed CentralView ArticleGoogle Scholar
- Haveman SA, Greene EA, Stilwell CP, Voordouw JK, Voordouw G: Physiological and gene expression analysis of inhibition of Desulfovibrio vulgaris Hildenborough by nitrite. J Bacteriol. 2004, 186: 7944-7950.PubMedPubMed CentralView ArticleGoogle Scholar
- Pereira IA, LeGall J, Xavier AV, Teixeira M: Characterization of a heme c nitrite reductase from a non-ammonifying microorganism, Desulfovibrio vulgaris Hildenborough. Biochim Biophys Acta. 2000, 1481: 119-130.PubMedView ArticleGoogle Scholar
- Rodrigues ML, Oliveira T, Matias PM, Martins IC, Valente FM, Pereira IA, Archer M: Crystallization and preliminary structure determination of the membrane-bound complex cytochrome c nitrite reductase from Desulfovibrio vulgaris Hildenborough. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2006, 62: 565-568.PubMedPubMed CentralView ArticleGoogle Scholar
- Koide T, Reiss DJ, Bare JC, Pang WL, Facciotti MT, Schmid AK, Pan M, Marzolf B, Van PT, Lo FY, Pratap A, Deutsch EW, Peterson A, Martin D, Baliga NS: Prevalence of transcription promoters within archaeal operons and coding sequences. Mol Syst Biol. 2009, 5: 285-PubMedPubMed CentralView ArticleGoogle Scholar
- Mukhopadhyay A, He Z, Alm EJ, Arkin AP, Baidoo EE, Borglin SC, Chen W, Hazen TC, He Q, Holman HY, Huang K, Huang R, Joyner DC, Katz N, Keller M, Oeller P, Redding A, Sun J, Wall J, Wei J, Yang Z, Yen HC, Zhou J, Keasling JD: Salt stress in Desulfovibrio vulgaris Hildenborough: an integrated genomics approach. J Bacteriol. 2006, 188: 4068-4078.PubMedPubMed CentralView ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210.PubMedPubMed CentralView ArticleGoogle Scholar
- The Meme Suite. [http://meme.nbcr.net/meme4_6_1/cgi-bin/meme.cgi]
- Rhodius VA, Suh WC, Nonaka G, West J, Gross CA: Conserved and variable functions of the sigmaE stress response in related genomes. PLoS Biol. 2006, 4: e2-PubMedPubMed CentralView ArticleGoogle Scholar
- Promscan @ Queens. [http://molbiol-tools.ca/promscan/]
- Studholme DJ, Buck M, Nixon T: Identification of potential sigma(N)-dependent promoters in bacterial genomes. Microbiology. 2000, 146: 3021-3023.PubMedView ArticleGoogle Scholar
- Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison CA, Smith HO: Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009, 6: 343-345.PubMedView ArticleGoogle Scholar
- Dunlop MJ, Dossani ZY, Szmidt HL, Chu HC, Lee TS, Keasling JD, Hadi MZ, Mukhopadhyay A: Engineering microbial biofuel tolerance and export using efflux pumps. Mol Syst Biol. 2011, 7: 487-PubMedPubMed CentralView ArticleGoogle Scholar
- Carver T, Berriman M, Tivey A, Patel C, Bohme U, Barrell BG, Parkhill J, Rajandream MA: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676.PubMedPubMed CentralView ArticleGoogle Scholar
- Hertz GZ, Stormo GD: Identifying DNA and protein patterns with statistically significant alignments of multiple sequences. Bioinformatics. 1999, 15: 563-577.PubMedView ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645.PubMedPubMed CentralView ArticleGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009, 26: 1641-1650.PubMedPubMed CentralView ArticleGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14: 1188-1190.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.