- Research article
- Open Access
Global transcriptional responses of Pseudomonas syringae DC3000 to changes in iron bioavailability in vitro
© Bronstein et al; licensee BioMed Central Ltd. 2008
- Received: 17 September 2008
- Accepted: 02 December 2008
- Published: 02 December 2008
Pseudomonas syringae pv tomato DC3000 (DC3000) is a Gram-negative model plant pathogen that is found in a wide variety of environments. To survive in these diverse conditions it must sense and respond to various environmental cues. One micronutrient required for most forms of life is iron. Bioavailable iron has been shown to be an important global regulator for many bacteria where it not only regulates a wide variety of genes involved in general cell physiology but also virulence determinants. In this study we used microarrays to study differential gene regulation in DC3000 in response to changes in levels of cell-associated iron.
DC3000 cultures were grown under highly controlled conditions and analyzed after the addition of iron citrate or sodium citrate to the media. In the cultures supplemented with iron, we found that cell-associated iron increased rapidly while culture densities were not significantly different over 4 hours when compared to cultures with sodium citrate added. Microarray analysis of samples taken from before and after the addition of either sodium citrate or iron citrate identified 386 differentially regulated genes with high statistical confidence. Differentially regulated genes were clustered based on expression patterns observed between comparison of samples taken at different time points and with different supplements. This analysis grouped genes associated with the same regulatory motifs and/or had similar putative or known function.
This study shows iron is rapidly taken up from the medium by iron-depleted DC3000 cultures and that bioavailable iron is a global cue for the expression of iron transport, storage, and known virulence factors in DC3000. Furthermore approximately 34% of the differentially regulated genes are associated with one of four regulatory motifs for Fur, PvdS, HrpL, or RpoD.
- Sigma Factor
- Regulatory Motif
- Hypothetical Gene
- Bioavailable Iron
- Tomato DC3000
Pseudomonas syringae pv tomato DC3000 (DC3000), a pathogen of tomato and Arabidopsis, is a Gram-negative bacterium that is found in a wide variety of environmental niches. This bacterium must to be able to interpret environmental cues to determine if it is in the soil, in water drops, on a plant surface, or in the plant apoplast. The pathogen adapts to each local environment by altering its internal physiology, including changing mRNA, protein content, and internal metabolite stores. By altering its behavior DC3000 is better able to scavenge available nutrients, withstand various environmental stresses and, if necessary, coordinate expression and deployment of appropriate virulence factors. For growth within plants DC3000 can express two major virulence factors: the phytotoxin coronatine, which causes chlorosis in plants , and a type III secretion system (T3SS), whose main function is to help the bacteria to evade host defenses by interfering with intracellular signaling pathways of plant cells (for recent review ). Several regulatory proteins have been identified that control T3SS and coronatine production but the precise nature of the chemical and physical factors used as environmental cues involved in the expression of these pathways are unknown.
One environmental signal detected by many bacterial pathogens, including Pseudomonas aeruginosa, is iron. Iron homeostasis in P. aeruginosa has been investigated due to its essential role in opportunistic human disease, especially the expression of extracellular virulence determinants including an exotoxin and proteases [3, 4]. Many iron-responsive regulatory pathways have been studied in this species. One key regulator is the ferric uptake regulator, Fur [5, 6]. Fur forms a homodimer complex in the presence of Zn2+ and Fe2+ ions [7–9] enabling it to bind to specific DNA motifs, "Fur boxes". These motifs are overlapping inverted repeat sequences that often occur in close proximity to, or overlap with, promoters. Fur binds to DNA in iron-replete conditions leading to repression of downstream transcriptional units. The characteristic repressor behavior of Fur is clearly evident in microarray studies of patterns of expression of iron-responsive genes in P. aeruginosa PAO1 (PAO1) [10, 11]. Recently the Fur regulon has been studied using computational methods to identify putative Fur regulated genes that specifically contain cooperative upstream binding sites .
The regulatory networks that respond to iron are complex. For instance in PAO1 there are several extracytoplasmic function (ECF) sigma factors and small RNAs (sRNAs) that are associated with iron homeostasis [11, 13, 14]. Twelve sigma factor genes have confirmed or putative upstream Fur binding sites, and nine of these genes exhibit iron-responsive expression under one or more conditions [10, 11]. Sigma factors modulate transcription by recruiting core RNA polymerase to promoters, which in turn regulate many genes. A recent study in PAO1 characterized several ECF sigma factors that control TonB-dependent receptors and a variety of putative metal transport systems .
Previous studies in PAO1 have implicated more than 300 genes that are differentially expressed in response in iron [10, 11]. Many of the regulatory elements and downstream targets found in PAO1 have homologs in other pseudomonads. However, relatively little direct experimentation has been reported for P. syringae pathovars regarding iron homeostasis. The iron-dependent production of toxins in P. syringae pv syringae has been known for some time [15, 16]. Loper and Lindow, using P. syringae 31-derived strains showed limited iron bioavailability in colonized regions of leaf surfaces . In addition, Joyner and Lindow found considerable variation in iron bioavailability on plant leaves . Analysis of the DC3000 genome sequence revealed the existence of two extracellular siderophores, yersiniabactin and pyoverdin, along with multiple putative TonB-dependent siderophore receptors, putative iron uptake pathways, and proteins involved in iron storage . Recent studies show that both yersiniabactin and pyoverdin are produced in DC3000 culture and that yersiniabactin is also produced in planta . Our group has also recently characterized the PvdS regulon, which is the regulon of an ECF type sigma factor that controls genes responsible for the production of pyoverdin along with other iron responsive genes . These studies imply that iron acquisition is important for growth in planta and that iron may be an important environmental signal sensed by DC3000. However, since a strain lacking in the production of yersiniabactin has no growth defect in the plant and an increased growth rate in low-iron media, the role of environmental iron signaling and acquisition is likely to be a complex phenomenon .
In this study we evaluated the effect of bioavailable iron on bacterial gene transcription using carefully controlled growth conditions to maximize differential cell associated iron levels and reduce confounding environmental and growth effects. We measured cell-associated iron concentrations in DC3000 at multiple time points after the addition of iron citrate or sodium citrate and evaluated the transcriptome to identify differentially regulated genes. RNA levels from each sample were analyzed using microarrays. We found that cell-associated iron levels are strongly associated with the differential expression of a variety of pathways in DC3000 including siderophore production, iron transport systems, sigma factors, T3SS and coronatine production. The differentially expressed genes were grouped by their patterns of regulation and we found that genes that clustered together had similar known or putative functions and shared bioinformatically derived regulatory motifs. Based on these results, we believe that many previously uncharacterized genes can be assigned putative functions and we present an initial systems view of iron-dependent regulation in DC3000.
Kinetics of iron uptake by iron-depleted P. syringae cultures
Reproducibility of culture growth and iron uptake
Global expression changes determined via microarray analysis
Clustering of differentially expressed genes
Evaluating gene clusters with regulatory motifs
After clustering we correlated the results with putative regulatory motifs generated by bioinformatics from other available data sets. Four motifs were used to search the DC3000 genome: the HrpL associated motif from Ferreira, et al., the PvdS motif from Swingle, et al., a Fur binding motif derived from PAO1 microarrays, and a RpoD motif generated by our group from unpublished proteomics data from DC3000 [10, 11, 21, 25]. Genes were considered associated with the HrpL or PvdS motifs if they were reported as such in the previous publications [21, 25]. For the RpoD and Fur motifs, genes were considered associated with these motifs if i) the motif was within 220 base pairs of the putative start of the first gene in the transcript (operon) containing the gene, as defined by a curated set of operons derived from VIMSS Operon Predictions , and ii) in the case of RpoD, if the motif and the putative transcript were in the same orientation (see Additional File 6 for RpoD and Fur motifs). Our analyses predicted 98 HrpL associated, 70 PvdS associated, 90 Fur associated, and 448 RpoD associated genes in the DC3000 genome.
Approximately 34% of the differentially expressed genes (140 of 386) were associated with one of these four motifs, and many genes associated with a particular regulatory motif grouped closely in the microarray-derived clusters. These include 16 differentially expressed operons (21 genes) that were associated with the Fur motif, 8 operons (13 genes) with the PvdS motif, 46 operons (74 genes) with HrpL, and 21 operons (32 genes) with the putative RpoD motif. As expected, Fur motifs were mostly associated with genes that showed differential expression at early time point comparisons after the addition of iron (8 out of 16 operons in the b-/b+ or a/b+ comparisons), while genes associated with the PvdS motif were associated with differential expression at later time points in response to iron (12 genes were differentially expressed in the c-/c+ comparison while none were identified in the b-/b+ comparison and only 1 was identified in the a/b+ comparison). Most of the HrpL associated genes are more highly expressed in iron supplemented media (a/c+ or c-/c+ comparisons). There appears to be a subgroup of HrpL associated genes that are also responsive to iron at earlier time points. The early set is diverse and includes some genes that encode for the Hrp secretion apparatus and Hrp secreted proteins. Additionally, the early set does not completely overlap with genes that were induced early in previous HrpL regulatory studies . Finally the RpoD motif is generally associated with genes that have lower expression at later time points (higher culture density). Many of these (27 of 33) have less than 2-fold differential expression, including genes that encode for ribosomal proteins, all of which show less than 1.2 fold change.
HrpL regulon cluster
Cluster I from Figure 5A contains a total of 91 differentially expressed genes that are induced in high iron conditions, 61 of which are associated with an HrpL regulatory motif. Other genes encode for bacterioferritins, superoxide dismutase, catalases, and peroxidase. There are also 6 hypothetical genes in this cluster and 4 known regulators: hrpRS, hrpL, and corR [27–29]. Additionally there are 3 putative response regulators and 1 putative sensor histidine kinase that may be involved in expression of virulence factors, iron storage or coping with free oxygen radicals.
While it is reported that iron influences T3SS in Shigella and Salmonella [30, 31], this is the first report of bioavailable iron regulating T3SS in the plant pathogen DC3000. The mechanism by which iron modulates T3SS seems to differ between Shigella and Salmonella. In Shigella a Fur regulated small RNA, RyhB, is implicated in regulating T3SS transcription ; while in Salmonella Fur, through an unknown mechanism, affects the activity of HilD, an AraC-like regulator, while in turn regulates expression of T3SS . A ryhB homolog in DC3000, prrF2, is associated with a Fur motif and was differentially regulated in response to iron in our experiments and experiments in P. aerginosa . If PrrF2 acts as a negative regulator of T3SS, which has been reported for RyhB , then we would expect prrF2 to co-cluster with T3SS regulatory genes in analyses using pair-wise average linkage with absolute un-centered correlation distance (Figure 5B). In this analysis we found that prrF2 clustered closely with T3SS genes including hrpL and hrpS, both known regulators of T3SS in DC3000. However, further experiments are required to characterize the role of prrF2 in the regulation of T3SS genes.
Coronatine production clusters
Cluster II contains 29 genes that are up-regulated in response to iron after 4 hours. Some encode known or putative enzymes that contain iron as a co-factor. These include genes involved in the production of succinate semialdehyde, a precursor for coronatine . While there are no putative regulatory proteins found in this cluster, 9 genes are described as hypothetical.
The 61 genes in cluster III are differentially expressed in the comparison between cultures supplemented with iron citrate at the 0 and 4 hour time points. Several T3SS substrates appear in this cluster as well as 4 genes involved in the formation of coronafacic acid, a precursor to coronatine, and one (gabD-3) that is involved in formation of succinate semialdehyde . Additionally there are genes encoding for porins, heat shock proteins, and many putative ABC transporters. About a third of the genes in this cluster are differentially expressed less than 1.5 fold. There are 15 hypothetical genes in the cluster, 3 of which have homology to functional domains of transcriptional regulators.
Growth phase dependent clusters
Four distinct clusters appear to be regulated in a growth phase dependent manner. The first of these, cluster IV, is up-regulated between 0 hour and 4 hour time points and/or 30 minute and 4 hour time points in the presence of iron citrate. Of the 41 genes in this cluster, 8 are differentially expressed less than 1.5 fold in both comparisons. These include several involved in general metabolism. Additionally there are two putative regulators in this set; one encodes for a PrtN-like regulator found only in Pseudomonas syringae pathovars (PSPTO_0570) and the other encodes for a protein homologous to OmpR. PrtN regulates pyocins in Pseudomonas aeruginosa [34, 35] and PSPTO_0570 is found within a region contains a putative F-type pyocin and at least one other gene, PSPTO_0569, clusters closely in PSPTO_0570 in this analysis. There are also two known regulators, OmpR, and RpoS. Since RpoS is reported to up-regulate genes at the onset of stationary phase [36, 37], we hypothesize that these regulators control the expression of many of the genes found in this cluster. However, we have not been able to derive regulatory motifs from the upstream regions of the genes in this cluster or other known RpoS regulated genes in other pseudomonads.
Cluster V contains 20 genes and 1 regulator, glnK, which controls expression of genes involved in nitrogen metabolism. Cluster V also includes two genes thought to be involved in nitrogen metabolism and several genes predicted to be involved in mannitol or amino acid transport [38, 39]. We hypothesize that GlnK is responsible for the differential regulation for many of these genes.
The next cluster (VI) contains 25 genes that are repressed at later time points. Seven show less than a 1.5 fold change and eight are associated with the RpoD motif. Many of these are involved in general metabolism. One interesting gene, PSPTO_4442, is homologous to bolA in E.coli. The BolA protein is hypothesized to be involved in the switch between cell elongation and septation and is induced during the transition into early stationary phase .
The last growth phase dependent gene cluster (VIII) contains 35 genes that are repressed between 0 hours and 4 hours after the addition of iron citrate. Almost half (15) show less than a 1.5 change between these two conditions and nine of those exhibit 1.2 fold or less change. Many are involved in ribosomal synthesis. This suggests that the cultures were beginning to transition from exponential growth and early stationary phase at this time point. 14 of the genes in this cluster are associated with the predicted RpoD motif. Because our RpoD motif contains a weak -35 box, a strong -10 box and a conserved motif directly downstream of the -10 box, this cluster may represent a subset of RpoD regulated genes that are subjected to additional regulation by other factors. Two other regulatory genes in this cluster, algU (rpoE) and sigX, are involved in outer membrane modifications and exhibit less than a 1.5 fold change in this comparison [40, 41].
Citrate uptake cluster
This cluster (VII) of 8 genes includes 3 hypothetical genes and is up-regulated starting at 30 minutes in response to the addition of either iron citrate or sodium citrate. Four genes are homologous to genes involved in citrate uptake (tctABC, and a citrate permease, citM) [42, 43]. While we cannot identify any regulatory proteins among these genes, we speculate that the 3 hypothetical genes that cluster in this group are also involved in citrate uptake.
Iron uptake and storage cluster
Cluster IX contains genes that are down-regulated in response to the addition of iron citrate. Thirteen of the twenty-three genes in this cluster are associated with a Fur motif. These include three genes that encode for putative lipoproteins, two hypothetical genes, and other genes whose homologs have been implicated as iron-responsive in other bacteria. Cluster IX also includes two regions in the genome predicted to contain small RNAs, both of which were included in the microarray analysis. One sRNA designated here as cobalamin-4, was identified by the Rfam database as a cobalamin riboswitch . The cobalamin riboswitch terminates transcription of downstream genes if B12 is bound to the riboswitch . In DC3000, cobalamin-4 is upstream of an operon containing PSPTO_3256-PSPTO_3258 that is not differentially regulated in our analyses, even at p = 0.05. The genes are predicted to encode for members of an iron ABC transport system  but without further analysis it is impossible to rule out that these genes are involved in the transport of cobalamin. Therefore, we hypothesize that the promoter upstream of cobalamin-4 is differentially expressed in response to iron but since cobalamin is still present in the cell in sufficient amounts (data not shown), B12 is binding to the RNA and attenuating the transcription of the downstream genes.
The other sRNA found in this cluster is prrF2. As stated previously prrF2 is a homolog of ryhB and has been characterized in P. aeruginosa for its role in iron homeostasis [13, 45]. In Figure 5B prrF2 clusters with genes that are similar to genes in P. aeruginosa that are affected by prrF1/2 including genes involved in catalase production, bacterioferritins, and superoxide dismutase . This leads us to believe that, in addition to possibly regulating T3SS, prrF1/2 play similar roles for iron homeostasis in DC3000 as they do in P. aeruginosa.
Also included within cluster IX are genes that encode for the TonB transport system (exbB-1, exbD-1,tonB-1) and putative TonB-dependent siderophore receptors (PSPTO_2152, 3294, 3574). Because of the co-clustering with the TonB transport system and the strong association with the Fur motif, we propose that the putative lipoproteins and hypothetical proteins encoded by genes in this cluster are involved in iron-uptake processes.
Siderophore production cluster
Cluster X contains 25 genes that are repressed in response to iron after 4 hours. Twenty of these genes are known or hypothesized to be associated with either yersiniabactin (8 genes) or pyoverdine (12 genes) production or siderophore transport/uptake [19, 20, 46]. Additionally the regulator for pyoverdine production, pvdS, is found within this cluster. There are also 2 genes that encode for putative TonB-dependent siderophore receptors (PSPTO_2605, 3462). PSPTO_2605 is located on the chromosome next to the yersiniabactin synthesis genes and is probably responsible for the uptake of that siderophore. Finally, fecB, involved in iron dicitrate transport, is also found within this group[47, 48]. We hypothesize that genes within this cluster are involved in iron acquisition and that intracellular iron levels at these time points are close to "iron-saturated". If so, we expect genes involved in iron uptake to be down-regulated by regulators downstream from Fur regulation, such sigma factors like PvdS.
In order to be a successful pathogen, a bacterium must sense and respond to a diverse array of environmental signals. Many signals cause the differential regulation of hundreds of genes by primary and downstream regulatory events. In this study we have investigated the connection between multiple bacterial regulons and iron availability using a systems biology approach by integrating global expression analysis with computational biology. By analyzing samples taken from cultures with different cell associated iron concentrations at multiple time points, we have attempted to unravel complex regulatory pathways. We found that clustering differentially expressed genes based on their patterns of expression grouped of genes with like function together. We also used regulatory motifs derived from other data sets to show that many closely grouped genes also share common regulatory features. This global study has allowed us to hypothesize on functions for many previously uncharacterized genes base on clustering and has given us an initial systems level view of gene regulation in response to bioavailable iron of Pseudomonas syringae pv tomato DC3000.
Media preparation and Bacterial growth in bioreactors
A Sixfors bioreactor system (Infors, Sweden) was used for culturing bacteria for microarray experiments. Reactors were soaked overnight in 20% nitric acid to removed residual bound iron. The 500 ml reactors were thoroughly washed and 400 ml of defined minimal medium [Mannitol-Glutamate (MG) media (10 g/L of mannitol, 2 g/L of L-glutamic acid, 0.5 g/L of KH2PO4, 0.2 g/L of NaCl, 0.2 g/L of MgSO4, final pH of 7)] was added . When available, Sigma Ultrapure components were used to minimize the amount of iron contamination in the media. Reactors were autoclaved and allowed to oxygenate for at least 4 hours prior to inoculation. Running conditions were as follows: 25°C, 1 L/min of air supplied via sparging, and a Rushton impeller spinning at 500 RPM for additional perturbation of the media. Vessels were inoculated to an OD600 of 0.01 with bacteria that had been grown to confluency on LM agar  and then resuspended in MG medium prior to inoculation.
Sample collection and Isolation of RNA
When cultures reached an OD600 of 0.3 (~16 hours after inoculation), samples were taken (t = 0 h) and iron citrate (Sigma, St Louis, MO) or sodium citrate (Sigma) was added to a final concentration of 50 μM. At each time point (t = 0 h, t = 0.5 h, and t = 4 h) 35 ml of culture was taken from each reactor via an aseptic method. Twenty-five ml of the cultures was centrifuged to separate bacteria from the supernatant and each fraction was frozen at -20°C for further analysis of iron levels. Five ml of culture was pelleted by centrifugation at room temperature for 5 min at 10,000 × g and the supernatant was removed. RNA was isolated using the RNeasy kit (Qiagen, Carlsbad, CA) following manufacture's instructions; with the exception that lyzozyme was used at a concentration of 5 mg/ml. RNA was treated with DNase I (Ambion, Austin, TX) to remove residual DNA and then cleaned and concentrated using the MinElute kit (Qiagen). Removal of DNA was verified by qRT-PCR . Integrity of the RNA was assessed using the Agilent Bioanalyzer (Microarray Core Facility, Cornell University).
Measurement of iron concentrations
Bacterial pellets from 25 ml of culture were digested with 1.0 ml of concentrated nitric acid at 120°C until dry, then 1.0 ml of a 1:1 mixture of concentrated nitric acid and perchloric acid was added and heated at 220°C until dry. The ash was dissolved in 20.0 ml of 5% Nitric acid and analyzed on an axially viewed ICP trace analyzer emission spectrometer (model ICAP 61E trace analyzer, Thermo Electron, Waltham Ma). The transfer optics were replaced with a short depth of field transfer optics to reduce matrix effects, (2000) US Patent No. 6,122050. A specialized spray chamber and desolvation system (2005) US Patent No. 6859272 was also used.
The microarrays used in this study were NimbleExpress Made-to-Order Prokaryotic Arrays for P. syringae pv tomato DC3000 (part# 530131) (NimbleGen Systems Inc., Madison, WI, USA) purchased through Affymetrix (Santa Clara, CA, USA). These microarrays consisted of 190,000 probes, with probe length of 24 nucleotides, representing 5,608 genes, as well as several intergenic regions, with a minimum of 17 probes/gene, including perfect match and mismatch probes. The NimbleExpress arrays are packaged in an Affymetrix GeneChip cartridge (49 format), and can be used with GeneChip reagents and processed on the GeneChip Instrument System. cDNA synthesis, labeling, and fragmentation were performed at the Microarray Core Facility at Cornell University using standard Affymetrix protocols. Hybridization, staining, washing, and chip scanning were also performed at the Microarray Core Facility using the Flexmidi_Euk2v3 fluidics protocol from Affymetrix with a modification in wash temperature B on a Fluidic station FS450. Microarray chips were scanned on a GeneChip scanner 3000-7G with 4-color upgrade and array autoloader.
Analysis of Microarray data
Microarray data were analyzed using Bioconductor . Quality assessment of microarray data was performed as previously described . Redundant and unwanted probes were removed from data sets prior to final statistical analysis. Samples were taken from 5 biological replicates for all time points. All data was normalized using gcrma and pair wise comparisons were made using genefilter requiring at least 3 chips where probe sets had intensity above 10. Subsequent FDR analysis was done using multtest at different p values (0.5, 0.1, 0.02, 0.004). Differentially expressed genes were clustered using the Bio.Cluster module in BioPython (version 1.44). The resulting dendograms were visualized using Java TreeView software .
Our group previously published HrpL and PvdS motifs [21, 25]. The microarray data on iron-responsive genes in PAO1 from Ochsner et al.  and Palma et al.  was used to identify putative Fur binding motifs in all pseudomonads. The process of identifying putative Fur binding sites is a variant of computational methods used previously by our group [25, 50, 53]. Briefly, training sets were constructed from the manually annotated set from Ochsner et al. and by performing Gibbs sampling with various constraints on motif length and symmetry on regions upstream of iron-responsive genes in the microarray experiments of Ochsner et al. and Palma et al. The HMM was calibrated and used to scan each available pseudomonad genome sequence. The RpoD promoter model and putative regulatory targets were derived from unpublished proteomics data. From the set of proteins that were present in a preliminary proteomics survey, a list of candidate genes and upstream regulatory elements were assembled, based on operon predictions used and reported in . Intergenic regions up to 90 bases upstream of identified operons were included in an input set, excluding regions identified in  to contain a HrpL promoter element (Hrp box) or that were upstream of an annotated transposase gene. Gibbs sampling using PhyloGibbs  was performed to identify promoter-size (approximately 30 nt), conserved motifs, and a motif was consistently found that resembled the known E. coli RpoD promoter motif . The aligned sequences were extracted and used to training an HMM for scanning genomic sequences to identify candidate promoters.
For Figure 6 and the associated p-values quoted in the text, the following computations were performed. For every pair of genes within the set of 386 differentially expressed genes shown in Figure 5A, the uncentered correlation (c) of the expression pattern of the genes in that pair was computed, using the same metric with which clustering was performed. A histogram of these correlations forms the blue background distribution labeled "all" in Figure 6(A–D). For each regulatory factor, all possible pairs were formed from the subset of genes predicted to contain an upstream regulatory motif for that factor. The histogram of correlations within each subset is shown in red for each factor. In addition, p-values were computed to estimate the likelihood that the subset of genes predicted to be regulated by a given factor could have been drawn at random from the background distribution of all pair-wise distances. For each factor, 100,000 random subsets of genes of the same size as the predicted subset were drawn from the full set of 386 genes, and all pair-wise correlations were computed for each random subset, from which the mean pair-wise correlation was extracted. The number of random samples with a mean pair-wise correlation that exceeded the mean pair-wise correlation of the predicted gene set for each factor was used to estimate a p-value for that factor. Given the trimodal nature of the background distribution, the mean pair-wise correlation is not a particularly discriminating test statistic, and we therefore expect the quoted p-values to be somewhat conservative, i.e., more discriminating test statistics would be expected to give smaller p-values than those quoted here for the means.
Microarray data accession
The microarray data from this study is available on the GEO database at http://www.ncbi.nlm.nih.gov/geo with the accession number GES13500.
We would like to thank Monica Moll for qRT-PCR analysis of selected targets, the Cornell University Life Sciences Core Laboratories Center/DNA microarray facility for processing our microarray chips, and Jenny Drnevich and Ariel Chernomoretz for portions of R scripts to modify probe sets for microarray analysis. The authors also thank Bryan Swingle for reviewing the manuscript prior to submission.
- Gnanamanickam SS, Starratt AN, Ward EWB: Coronatine Production Invitro and Invivo and Its Relation to Symptom Development in Bacterial-Blight of Soybean. Canadian of Bot. 1982, 60 (5): 645-650.View ArticleGoogle Scholar
- McCann HC, Guttman DS: Evolution of the type III secretion system and its effectors in plant-microbe interactions. New Phytol. 2008, 177 (1): 33-47.View ArticlePubMedGoogle Scholar
- Prince RW, Cox CD, Vasil ML: Coordinate regulation of siderophore and exotoxin A production: molecular cloning and sequencing of the Pseudomonas aeruginosa fur gene. J Bacteriol. 1993, 175 (9): 2589-2598.PubMed CentralPubMedGoogle Scholar
- Wilderman PJ, Vasil AI, Johnson Z, Wilson MJ, Cunliffe HE, Lamont IL, Vasil ML: Characterization of an endoprotease (PrpL) encoded by a PvdS-regulated gene in Pseudomonas aeruginosa. Infect Immun. 2001, 69 (9): 5385-5394. 10.1128/IAI.69.9.5385-5394.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Escolar L, Perez-Martin J, de Lorenzo V: Opening the iron box: transcriptional metalloregulation by the Fur protein. J Bacteriol. 1999, 181 (20): 6223-6229.PubMed CentralPubMedGoogle Scholar
- Litwin CM, Calderwood SB: Role of iron in regulation of virulence genes. Clin Microbiol Rev. 1993, 6 (2): 137-149.PubMed CentralPubMedGoogle Scholar
- Bagg A, Neilands JB: Ferric uptake regulation protein acts as a repressor, employing iron (II) as a cofactor to bind the operator of an iron transport operon in Escherichia coli. Biochemistry. 1987, 26 (17): 5471-5477. 10.1021/bi00391a039.View ArticlePubMedGoogle Scholar
- Jacquamet L, Aberdam D, Adrait A, Hazemann JL, Latour JM, Michaud-Soret I: X-ray absorption spectroscopy of a new zinc site in the fur protein from Escherichia coli. Biochemistry. 1998, 37 (8): 2564-2571. 10.1021/bi9721344.View ArticlePubMedGoogle Scholar
- Althaus EW, Outten CE, Olson KE, Cao H, O'Halloran TV: The ferric uptake regulation (Fur) repressor is a zinc metalloprotein. Biochemistry. 1999, 38 (20): 6559-6569. 10.1021/bi982788s.View ArticlePubMedGoogle Scholar
- Palma M, Worgall S, Quadri LE: Transcriptome analysis of the Pseudomonas aeruginosa response to iron. Arch Microbiol. 2003, 180 (5): 374-379. 10.1007/s00203-003-0602-z.View ArticlePubMedGoogle Scholar
- Ochsner UA, Wilderman PJ, Vasil AI, Vasil ML: GeneChip expression analysis of the iron starvation response in Pseudomonas aeruginosa: identification of novel pyoverdine biosynthesis genes. Mol Microbiol. 2002, 45 (5): 1277-1287. 10.1046/j.1365-2958.2002.03084.x.View ArticlePubMedGoogle Scholar
- van Oeffelen L, Cornelis P, van Delm W, De Ridder F, De Moor B, Moreau Y: Detecting cis-regulatory binding sites for cooperatively binding proteins. Nucleic Acids Res. 2008, 36 (8): e46-10.1093/nar/gkn140.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilderman PJ, Sowa NA, FitzGerald DJ, FitzGerald PC, Gottesman S, Ochsner UA, Vasil ML: Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis. Proc Natl Acad Sci USA. 2004, 101 (26): 9792-9797. 10.1073/pnas.0403423101.PubMed CentralView ArticlePubMedGoogle Scholar
- Llamas MA, Mooij MJ, Sparrius M, Vandenbroucke-Grauls CM, Ratledge C, Bitter W: Characterization of five novel Pseudomonas aeruginosa cell-surface signalling systems. Mol Microbiol. 2008, 67 (2): 458-472.View ArticlePubMedGoogle Scholar
- Morgan MK, Chatterjee AK: Genetic organization and regulation of proteins associated with production of syringotoxin by Pseudomonas syringae pv. syringae. J Bacteriol. 1988, 170 (12): 5689-5697.PubMed CentralPubMedGoogle Scholar
- Gross DC: Regulation of syringomycin synthesis in Pseudomonas syringae pv. syringae and defined conditions for its production. J Appl Bacteriol. 1985, 58 (2): 167-174.View ArticlePubMedGoogle Scholar
- Loper JE, Lindow SE: A Biological Sensor for Iron Available to Bacteria in Their Habitats on Plant Surfaces. Appl Environ Microbiol. 1994, 60 (6): 1934-1941.PubMed CentralPubMedGoogle Scholar
- Joyner DC, Lindow SE: Heterogeneity of iron bioavailability on plants assessed with a whole-cell GFP-based bacterial biosensor. Microbiology. 2000, 146 (Pt 10): 2435-2445.View ArticlePubMedGoogle Scholar
- Buell CR, Joardar V, Lindeberg M, Selengut J, Paulsen IT, Gwinn ML, Dodson RJ, Deboy RT, Durkin AS, Kolonay JF, Madupu R, Daugherty S, Brinkac L, Beanan MJ, Haft DH, Nelson WC, Davidsen T, Zafar N, Zhou LW, Liu J, Yuan QP, Khouri H, Fedorova N, Tran B, Russell D, Berry K, Utterback T, Van Aken SE, Feldblyum TV, D'Ascenzo M, Deng WL, Ramos AR, Alfano JR, Cartinhour S, Chatterjee AK, Delaney TP, Lazarowitz SG, Martin GB, Schneider DJ, Tang XY, Bender CL, White O, Fraser CM, Collmer A: The complete genome sequence of the Arabidopsis and tomato pathogen Pseudomonas syringae pv. tomato DC3000. Proc Natl Acad Sci USA. 2003, 100 (18): 10181-10186. 10.1073/pnas.1731982100.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones AM, Lindow SE, Wildermuth MC: Salicylic acid, yersiniabactin, and pyoverdin production by the model phytopathogen Pseudomonas syringae pv. tomato DC3000: synthesis, regulation, and impact on tomato and Arabidopsis host plants. J Bacteriol. 2007, 189 (19): 6773-6786. 10.1128/JB.00827-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Swingle B, Thete D, Moll M, Myers CR, Schneider DJ, Cartinhour S: Characterization of the PvdS-regulated promoter motif in Pseudomonas syringae pv. tomato DC3000 reveals regulon members and insights regarding PvdS function in other pseudomonads. Mol Microbiol. 2008, 68 (4): 871-889. 10.1111/j.1365-2958.2008.06209.x.View ArticlePubMedGoogle Scholar
- Keane PJ, Kerr A, New PB: Crown Gall of Stone Fruit .2. Identification and Nomenclature of Agrobacterium Isolates. Austn J Biol Sci. 1970, 23 (3): 585-&.Google Scholar
- Heber S, Sick B: Quality assessment of Affymetrix GeneChip data. Omics. 2006, 10 (3): 358-368. 10.1089/omi.2006.10.358.View ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR: Rfam: an RNA family database. Nucleic Acids Res. 2003, 31 (1): 439-441. 10.1093/nar/gkg006.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferreira AO, Myers CR, Gordon JS, Martin GB, Vencato M, Collmer A, Wehling MD, Alfano JR, Moreno-Hagelsieb G, Lamboy WF, DeClerck G, Schneider DJ, Cartinhour SW: Whole-genome expression profiling defines the HrpL regulon of Pseudomonas syringae pv. tomato DC allows de novo reconstruction of the Hrp cis clement, and identifies novel coregulated genes. Mol Plant Microbe Interact. 3000, 19 (11): 1167-1179. 10.1094/MPMI-19-1167.View 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 (3): 880-892. 10.1093/nar/gki232.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao Y, Heu S, Yi J, Lu Y, Hutcheson SW: Identification of a putative alternate sigma factor and characterization of a multicomponent regulatory cascade controlling the expression of Pseudomonas syringae pv. syringae Pss61 hrp and hrmA genes. J Bacteriol. 1994, 176 (4): 1025-1036.PubMed CentralPubMedGoogle Scholar
- Rahme LG, Mindrinos MN, Panopoulos NJ: Plant and environmental sensory signals control the expression of hrp genes in Pseudomonas syringae pv. phaseolicola. J Bacteriol. 1992, 174 (11): 3499-3507.PubMed CentralPubMedGoogle Scholar
- Ullrich M, Penaloza-Vazquez A, Bailey AM, Bender CL: A modified two-component regulatory system is involved in temperature-dependent biosynthesis of the Pseudomonas syringae phytotoxin coronatine. J Bacteriol. 1995, 177 (21): 6160-6169.PubMed CentralPubMedGoogle Scholar
- Murphy ER, Payne SM: RyhB, an iron-responsive small RNA molecule, regulates Shigella dysenteriae virulence. Infect Immun. 2007, 75 (7): 3470-3477. 10.1128/IAI.00112-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Ellermeier JR, Slauch JM: Fur regulates expression of the Salmonella pathogenicity island 1 type III secretion system through HilD. J Bacteriol. 2008, 190 (2): 476-486. 10.1128/JB.00926-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Masse E, Gottesman S: A small RNA regulates the expression of genes involved in iron metabolism in Escherichia coli. Proc Natl Acad Sci USA. 2002, 99 (7): 4620-4625. 10.1073/pnas.032066599.PubMed CentralView ArticlePubMedGoogle Scholar
- Rangaswamy V, Jiralerspong S, Parry R, Bender CL: Biosynthesis of the Pseudomonas polyketide coronafacic acid requires monofunctional and multifunctional polyketide synthase proteins. Proc Natl Acad Sci USA. 1998, 95 (26): 15469-15474. 10.1073/pnas.95.26.15469.PubMed CentralView ArticlePubMedGoogle Scholar
- Shinomiya T, Shiga S, Kikuchi A, Kageyama M: Genetic determinant of pyocin R2 in Pseudomonas aeruginosa PAO. II. Physical characterization of pyocin R2 genes using R-prime plasmids constructed from R68.45. Mol Gen Genet. 1983, 189 (3): 382-389. 10.1007/BF00325899.View ArticlePubMedGoogle Scholar
- Matsui H, Sano Y, Ishihara H, Shinomiya T: Regulation of pyocin genes in Pseudomonas aeruginosa by positive (prtN) and negative (prtR) regulatory genes. J Bacteriol. 1993, 175 (5): 1257-1263.PubMed CentralPubMedGoogle Scholar
- Fujita M, Tanaka K, Takahashi H, Amemura A: Transcription of the principal sigma-factor genes, rpoD and rpoS, in Pseudomonas aeruginosa is controlled according to the growth phase. Mol Microbiol. 1994, 13 (6): 1071-1077. 10.1111/j.1365-2958.1994.tb00498.x.View ArticlePubMedGoogle Scholar
- Miksch G, Dobrowolski P: Growth phase-dependent induction of stationary-phase promoters of Escherichia coli in different gram-negative bacteria. J Bacteriol. 1995, 177 (18): 5374-5378.PubMed CentralPubMedGoogle Scholar
- van Heeswijk WC, Stegeman B, Hoving S, Molenaar D, Kahn D, Westerhoff HV: An additional PII in Escherichia coli: a new regulatory protein in the glutamine synthetase cascade. FEMS Microbiol Lett. 1995, 132: 1-2. 10.1016/0378-1097(95)00271-6.View ArticleGoogle Scholar
- Atkinson MR, Ninfa AJ: Role of the GlnK signal transduction protein in the regulation of nitrogen assimilation in Escherichia coli. Mol Microbiol. 1998, 29 (2): 431-447. 10.1046/j.1365-2958.1998.00932.x.View ArticlePubMedGoogle Scholar
- Martin DW, Holloway BW, Deretic V: Characterization of a locus determining the mucoid status of Pseudomonas aeruginosa: AlgU shows sequence similarities with a Bacillus sigma factor. J Bacteriol. 1993, 175 (4): 1153-1164.PubMed CentralPubMedGoogle Scholar
- Brink S, Bogsch EG, Edwards WR, Hynds PJ, Robinson C: Targeting of thylakoid proteins by the Delta pH-driven twin-arginine translocation pathway requires a specific signal in the hydrophobic domain in conjunction with the twin-arginine motif. Febs Letters. 1998, 434 (3): 425-430. 10.1016/S0014-5793(98)01028-X.View ArticlePubMedGoogle Scholar
- Winnen B, Hvorup RN, Saier MH: The tripartite tricarboxylate transporter (TTT) family. Res Microbiol. 2003, 154 (7): 457-465. 10.1016/S0923-2508(03)00126-8.View ArticlePubMedGoogle Scholar
- Boorsma A, Rest van der ME, Lolkema JS, Konings WN: Secondary transporters for citrate and the Mg(2+)-citrate complex in Bacillus subtilis are homologous proteins. J Bacteriol. 1996, 178 (21): 6216-6222.PubMed CentralPubMedGoogle Scholar
- Nahvi A, Barrick JE, Breaker RR: Coenzyme B12 riboswitches are widespread genetic control elements in prokaryotes. Nucleic Acids Res. 2004, 32 (1): 143-150. 10.1093/nar/gkh167.PubMed CentralView ArticlePubMedGoogle Scholar
- Oglesby AG, Farrow JM, Lee JH, Tomaras AP, Greenberg EP, Pesci EC, Vasil ML: The influence of iron on Pseudomonas aeruginosa physiology: a regulatory link between iron and quorum sensing. J Biol Chem. 2008, 283 (23): 15558-15567. 10.1074/jbc.M707840200.PubMed CentralView ArticlePubMedGoogle Scholar
- Bultreys A, Gheysen I, de Hoffmann E: Yersiniabactin production by Pseudomonas syringae and Escherichia coli, and description of a second yersiniabactin locus evolutionary group. Appl Environ Microbiol. 2006, 72 (6): 3814-3825. 10.1128/AEM.00119-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Pressler U, Staudenmaier H, Zimmermann L, Braun V: Genetics of the iron dicitrate transport system of Escherichia coli. J Bacteriol. 1988, 170 (6): 2716-2724.PubMed CentralPubMedGoogle Scholar
- Staudenmaier H, van Hove B, Yaraghi Z, Braun V: Nucleotide sequences of the fecBCDE genes and locations of the proteins suggest a periplasmic-binding-protein-dependent transport mechanism for iron(III) dicitrate in Escherichia coli. J Bacteriol. 1989, 171 (5): 2626-2633.PubMed CentralPubMedGoogle Scholar
- Hanahan D: Studies on Transformation of Escherichia-Coli with Plasmids. J Mol Biol. 1983, 166 (4): 557-580. 10.1016/S0022-2836(83)80284-8.View ArticlePubMedGoogle Scholar
- Vencato M, Tian F, Alfano JR, Buell CR, Cartinhour S, DeClerck GA, Guttman DS, Stavrinides J, Joardar V, Lindeberg M, Bronstein PA, Mansfield JW, Myers CR, Collmer A, Schneider DJ: Bioinformatics-enabled identification of the HrpL regulon and type III secretion system effector proteins of Pseudomonas syringae pv. phaseolicola 1448A. Mol Plant Microbe Interact. 2006, 19 (11): 1193-1206. 10.1094/MPMI-19-1193.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMedGoogle Scholar
- Saldanha AJ: Java Treeview – extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- Fouts DE, Abramovitch RB, Alfano JR, Baldo AM, Buell CR, Cartinhour S, Chatterjee AK, D'Ascenzo M, Gwinn ML, Lazarowitz SG, Lin NC, Martin GB, Rehm AH, Schneider DJ, van Dijk K, Tang X, Collmer A: Genomewide identification of Pseudomonas syringae pv. tomato DC3000 promoters controlled by the HrpL alternative sigma factor. Proc Natl Acad Sci USA. 2002, 99 (4): 2275-2280. 10.1073/pnas.032514099.PubMed CentralView ArticlePubMedGoogle Scholar
- Siddharthan R, Siggia ED, van Nimwegen E: PhyloGibbs: a Gibbs sampling motif finder that incorporates phylogeny. PLoS Comput Biol. 2005, 1 (7): e67-10.1371/journal.pcbi.0010067.PubMed CentralView ArticlePubMedGoogle Scholar
- Wagner R: Transcription regulation in prokaryotes. 2000, Oxford; New York: Oxford University PressGoogle 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.