PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(5): e10487.
Published online 2010 May 5. doi:  10.1371/journal.pone.0010487
PMCID: PMC2864753

Transcriptional and Translational Regulatory Responses to Iron Limitation in the Globally Distributed Marine Bacterium Candidatus Pelagibacter ubique

Francisco Rodriguez-Valera, Editor

Abstract

Iron is recognized as an important micronutrient that limits microbial plankton productivity over vast regions of the oceans. We investigated the gene expression responses of Candidatus Pelagibacter ubique cultures to iron limitation in natural seawater media supplemented with a siderophore to chelate iron. Microarray data indicated transcription of the periplasmic iron binding protein sfuC increased by 16-fold, and iron transporter subunits, iron-sulfur center assembly genes, and the putative ferroxidase rubrerythrin transcripts increased to a lesser extent. Quantitative peptide mass spectrometry revealed that sfuC protein abundance increased 27-fold, despite an average decrease of 59% across the global proteome. Thus, we propose sfuC as a marker gene for indicating iron limitation in marine metatranscriptomic and metaproteomic ecological surveys. The marked proteome reduction was not directly correlated to changes in the transcriptome, implicating post-transcriptional regulatory mechanisms as modulators of protein expression. Two RNA-binding proteins, CspE and CspL, correlated well with iron availability, suggesting that they may contribute to the observed differences between the transcriptome and proteome. We propose a model in which the RNA-binding activity of CspE and CspL selectively enables protein synthesis of the iron acquisition protein SfuC during transient growth-limiting episodes of iron scarcity.

Introduction

The importance of iron as a nutrient in the oceans was first recognized by Martin [1] and later experiments verified that iron limits primary production over broad regions of the marine environment [2][4]. A variety of biological processes such as photosynthesis, N2 fixation, methanogenesis, respiration, oxygen transport, gene regulation, and DNA synthesis all depend on iron-containing proteins [5]. In pelagic surface waters, planktonic communities must cope with iron concentrations that average just 70 picomolar [6]. The inhibitory effect that this has on growth was most clearly illustrated by a series of iron fertilization experiments in which iron was added to large swaths of the ocean, resulting in a marked increase in nutrient utilization [2], [3], [7].

Bacteria commonly have specialized systems for responding to iron limitation. Genes for iron uptake and utilization are primarily regulated by the Fur protein [8], [9]. When complexed with Fe(II) cations, Fur binds the “Fur box” recognition sequence, which is made of several GATAAT hexamers [10][14]. In some bacteria, this single transcription factor can directly repress or activate more than 100 genes in response to iron scarcity [9]. Irr is a similar transcription factor that couples intracellular heme levels to expression of many different iron-related pathways [15][21]. Small RNAs [22][25] and mRNA-binding proteins [26], [27] can also regulate nonessential iron-utilizing proteins at the post-transcriptional level by selectively targeting their transcripts for degradation. To improve their chances of encountering Fe(III), many bacteria secrete siderophores [28][30]. These chelating agents help dissolve the poorly soluble particles and sequester them in a form that is unusable by competing microorganisms. Due to the spontaneous reactivity of iron ions, cells often encapsulate these atoms inside containers made of ferritin proteins to better modulate redox reactions [31], [32].

Candidatus Pelagibacter ubique was selected as an iron limitation model for two reasons. First, this alphaproteobacterium is regularly the most numerically abundant microorganism in surveys of marine microbial diversity. Second, its proteome of just 1,354 genes is possibly the simplest of any free-living heterotrophic organism [33]. Ca. Pelagibacter ubique's genome encodes Fur and Irr, but not ferritin or siderophore-related proteins, raising questions about how or if this bacterium can cope with iron stress. Investigating how this organism's relatively small genetic repertoire produces thriving populations in the variable ocean environment has been impeded by the lack of a genetic system able to create knockouts or other genetic modifications. Thus, observing how the entire transcriptome or proteome changes in response to growth conditions has become a primary approach for elucidating metabolic and regulatory schemes [34][36]. A comparison of cultures in exponential and stationary phase did not reveal a major remodeling of the proteome nor evidence of a global regulatory mechanism [34], suggesting that this organism may continue to benefit from temporary nutrient availability regardless of overall cellular activity. That study, along with a follow-up using a metaproteomics approach on environmental samples [35], found Ca. Pelagibacter ubique's proteome to be consistently composed of an unusually high proportion of transport-related proteins.

Arguably the most important characteristic of organisms is their ability to express the right proteins in the right amounts at the right times. The interplay between stimuli, sensors, and regulators precisely optimizes the combination of mRNA transcripts and protein products present in the cell. Several known and putative transcriptional regulators have been identified in the Ca. Pelagibacter ubique genome, as well as cis-acting riboswitches capable of modulating mRNA translation based on the concentration of particular metabolites. This method of decoupling production of mRNA from protein synthesis has been found on glycine metabolism genes in Ca. Pelagibacter ubique [37] and sequence motif searches [38] found additional candidate riboswitches for s-adenylmethionine [39], [40] and thiamine pyrophosphate [40], [41]. Meyer and colleagues also identified homologs to ribosomal proteins capable of regulating their own translation, as well as regions in the genome with riboswitch-like characteristics but lacking homologous annotated motifs. One of these putative structural RNA regions is located immediately upstream of the sfuA–C operon, which encodes an iron-acquisition system. Post-transcriptional regulation schemes allow the cell to conserve amino acids while still rapidly providing ephemeral enzymes. The success of these characteristics are evidenced by direct observation; Ca. Pelagibacter ubique is the most abundant heterotroph in the oceans, accounting for one-third of surface water bacteria [42], [43] and consuming up to half of some dissolved organic matter compounds [44].

Combining transcriptomic and proteomic data offers a perspective on cellular activity that cannot be obtained from either method individually. Numerous studies have shown that changes in the transcriptome poorly correlate with changes to the proteome, except for very highly expressed genes [45][48]. Although much of the disparity between these two types of datasets has been attributed to measurement inaccuracy [49] and differences in protein degradation rates [50], some studies have revealed systematic post-transcriptional regulatory schemes. For instance, in the eukaryotic protozoan Plasmodium falciparum, mRNAs were often upregulated an entire life phase before the one in which the encoded protein was needed [51]. Additionally, translation in Escherichia coli was found to be partially regulated by mRNA secondary structure [52]. Therefore, it is evident that the transcriptome does not necessarily represent the current state of the proteome, but is rather a mixture of transcripts being actively translated and others that are standing by, awaiting activation by post-transcriptional regulatory mechanisms. This study integrates both transcriptomic and proteomic analyses in order to attain a more complete understanding of the cellular response to iron limitation in Ca. Pelagibacter ubique. The results strongly suggest that transcription and translation are not always tightly coupled in this bacterium.

Results

Reaction to the Siderophore

Two iron-sequestering siderophores were tested on cultures of Ca. Pelagibacter ubique to determine their feasibility for creating iron-limiting conditions. Ferrichrome (Sigma #F8014) and deferoxamine mesylate salt (Sigma #D9533) were both found to arrest batch culture growth within 1/3 of a doubling – an inhibition which could be reversed by addition of iron (Figure 1A). Sufficient bioavailable iron was present in the natural seawater media collected from the Oregon coast to enable cultures to grow when supplemented with 10 nM siderophore, but not when supplemented with 100 nM siderophore.

Figure 1
Growth of Candidatus Pelagibacter ubique cells was arrested by iron-sequestering siderophores.

Microarrays

Six 20 L carboys inoculated with Ca. Pelagibacter ubique were grown to near-maximum density, then randomly selected for treatment with either ferrichrome or ferrichrome plus excess iron (Figure 1B). To measure the amount of messenger RNA transcripts present in cells, mRNA from each carboy was hybridized to separate microarray chips containing probes for all Ca. Pelagibacter ubique genes. Microarray data was deposited in the NCBI GEO database under accession number GSE20962. Of the three time points where mRNA abundance was measured, the greatest difference in expression of known iron-related genes was observed 24 hours after the siderophore amendment. Table 1 lists the 23 transcripts that were expressed at least 50 percent higher in the iron-limited culture relative to the control. Two-thirds of the genes in this list come from two operons: the first containing iron-sulfur center assembly proteins including sufA–E and the second made up of iron uptake proteins such as sfuA–C. The four other genes with a known function are: rubrerythrin, an iron-binding protein that is postulated to act as a ferroxidase for converting Fe(II) to Fe(III), hslU and hslV which together form a protease complex, and dimethylglycine dehydrogenase – an enzyme that is necessary for converting betaine to glycine.

Table 1
All 23 Ca. Pelagibacter ubique mRNA transcripts that were at least 50 percent more abundant in the iron-limited cultures compared to the control cultures, 24 hours after addition of an iron-chelator.

A modified radial coordinate visualization plot (Figure 2) of the microarray data shows four distinct clusters of genes: exponential growth, stationary phase, early iron stress, and late iron stress. As detailed in Supplementary Table S1, the early iron stress cluster is dominated by genes from three genomic loci: the sfu iron uptake operon, the suf iron-sulfur center assembly operon, and the functionally unclear loci SAR11_1157, SAR11_1158, SAR11_1163, and SAR11_1164. The late iron stress cluster also contains different genes that are located in or adjacent to the sfu and suf operons, but is better characterized by lexA, recA, and mucA – three genes involved in the SOS response. The iron response regulators fur and irr cluster with stationary phase genes, indicating that the abundances of these two transcripts are more affected by stationary phase than by iron limitation.

Figure 2
Genes transcribed during iron limitation were different from stationary phase genes.

Proteomics

Cellular protein fractions from each treatment were isolated and digested before being separated with liquid chromatography and injected into a tandem mass spectrometer. An Accurate Mass and Time Tag library, developed previously [34], was used to make quantitative comparisons of the abundance of individual peptides between samples. This dataset is available at http://omics.pnl.gov/. Of the 216 proteins detected with high certainty in this study, 18 were observed to be at least 50% more abundant in the iron-limited cultures: four on day 18, and 17 on day 28 (Table 2). The proteins SfuC, CspL, and GroES were higher in the iron-limited cultures at both timepoints. The iron-binding SfuC is unique in that it was the only one of these 18 proteins to increase in both protein and mRNA abundance by at least 50%. CspL was originally annotated as a DNA-binding protein, however, similar proteins have been found to modulate the accessibility of mRNA binding sites by selectively melting secondary RNA structures [53]. The third protein, GroES, forms a complex with GroEL to mediate protein folding. Because the required GroEL subunit was much less abundant in iron-limited cultures, and since the three largest GroES peptide spectra (out of 9) were less pronounced in the iron-limited cultures, GroES may be a false positive. Mass spectrometry measurements did not reveal a signficant change in Fur or Irr abundance between treatments or timepoints.

Table 2
All 17 Ca. Pelagibacter ubique proteins that were at least 50 percent more abundant in the iron-limited cultures compared to the iron replete cultures, two and 12 days after addition of an iron chelator.

Iron-limitation had a marked impact on the overall proteome. Two days after addition of an iron-chelator, 181 of the 216 proteins were significantly (P≤0.05) less abundant in the iron-limited cultures relative to the control cultures. Using the same criteria, only 32 of the 216 proteins were found to significantly decrease in the control cultures between days 18 and 28 as Ca. Pelagibacter ubique cells entered stationary phase due to an unknown, non-iron, limitation.

Comparing Changes in mRNA and Protein Abundances

Aside from the highly expressed iron-binding protein SfuC, the abundances of individual proteins appeared to be independent of the amount of mRNA encoding them (Figure 3).

Figure 3
Protein abundances were largely decoupled from transcript abundances.

Discussion

We are studying keystone microbial plankton species such as Ca. Pelagibacter ubique in culture to provide a basis for interpreting data emerging from molecular ecology studies. In an era of rapid environmental change, metagenomics, and allied technologies such as metaproteomics and metatranscriptomics, are being used to monitor the structure and health of natural ecosystems and to identify ecological processes that impact biogeochemistry. Interpretations of these data depend on understanding how complex cellular systems respond to environmental factors. We focused on a microorganism, Ca. Pelagibacter ubique, that produces the largest signal in most environmental studies of marine macromolecules, and a process, iron limitation, that impacts marine ecology on very large geographical scales.

Upregulation of sfuC during iron limitation

The only gene to clearly increase in both mRNA and protein abundance during iron limitation was sfuC. This protein localizes to the periplasmic space and binds dissolved Fe(III) with high affinity. The SfuC-Fe complex associates with the ATPase (SfuA) and permease (SfuB) components of the tripartite ABC transporter complex to actively transport iron into the cell. The fact that sfuA and sfuB were not observed to increase in protein abundance is not wholly unexpected – SfuA–SfuB complexes only interact with iron-bound SfuC proteins, which are a very small fraction of the total SfuC pool in an iron-limited environment. Additionally, integral membrane proteins such as SfuB are particularly challenging to recover in proteomic studies because they are not readily soluble. This likely contributed to the complete absence of SfuB peptides in all mass spectrometry studies of Ca. Pelagibacter ubique to date.

The identification of sfuC expression as a readily quantifiable iron limitation marker is particularly useful for ecological surveys. As its name suggests, Ca. Pelagibacter ubique's genome, transcriptome, and proteome regularly dominate bacterial surveys throughout the pelagic environment. Future oceanographic studies seeking evidence of iron availability limiting bacterioplankton growth may use metatranscriptomic or metaproteomic analyses to assess the expression of sfuC in the local Ca. Pelagibacter ubique population.

Transcriptome distinct from proteome

Protein abundance was generally uncorrelated with changes in mRNA abundance, suggesting that post-transcriptional mechanisms might be acting at the RNA level to suppress translation. As reviewed in the introduction, previous studies have shown that disparities between a cell's transcriptome and proteome are the norm rather than the exception. However, the observation that iron-related genes such as sufA–E increased in mRNA but not protein indicates that expression of these proteins are controlled at both the level of transcription and at the level of translation.

Cold-shock proteins correlated with iron stress

CspL was significantly more abundant in iron-limited cultures (Figure 4), leading us to closely examine the biological activity of this protein as well as the inversely expressed homolog CspE. The first discovered member of the cold-shock protein (CSP) family, E. coli's CspA, is highly upregulated under cold stress; it is believed to associate with and melt double-stranded RNA complexes as a mechanism to prevent spurious stem loop structures from interfering with transcription and translation [54][59]. Despite their homology to CspA, many CSP variants are not cold-inducible, but rather are involved in regulating cellular processes [54], [60][62] and can even target their activity to specific RNA sequences [63]. A growing body of literature has described mRNAs which modulate their own expression via temperature- (RNA thermometer) or ligand-sensitive (riboswitch) secondary structures [64], [65]. Due to the episodic nature of iron deposition into ocean surface waters [66] and the resulting selective pressure favoring rapid response systems for this limiting nutrient [67], we speculate that Ca. Pelagibacter ubique CspE and/or CspL affects a reversible inhibition of translation by facilitating an mRNA secondary structure unfavorable for ribosome processing, thereby maintaining the transcriptome in a state of cell growth readiness during times of stress such as iron limitation.

Figure 4
Translation of Ca. Pelagibacter ubique's cold shock and iron-binding genes are influenced by iron availability.

This is the first report describing the general suppression of translation across the entire transcriptome of the cell. In this case, the apparent adaptive significance of protein synthesis suppression is related to urgent cellular requirements to acquire an essential nutrient. The model we propose to explain this phenomenon incorporates activity previously observed in cold-shock proteins, however, the essence of our model assigns cold-shock proteins a new systemic role in Ca. Pelagibacter ubique cells with the apparent result of focusing protein synthesis on transporters that target a missing essential nutrient. The validation of this model is beyond the scope of this study. Future work may more precisely identify interactions between cold-shock proteins and specific RNA motifs.

Summary

Census information has left little doubt that Ca. Pelagibacter ubique plays its role in biological oceanography on a vast scale. To understand this role, we turned inward, investigating the mechanisms used by these cells to respond to a common form of nutritional stress. One motivation for this study can be described with a term borrowed from satellite remote sensing: the term “ground truth” was coined to describe the validation, by direct measurements, of remotely sensed observations. Metatranscriptomic and metaproteomic measurements are being widely adopted by microbial ecologists anticipating that these approaches will reveal the metabolic status of cells in microbial communities, providing information that can be extrapolated to interpret broader levels of ecosystem function. Essential to this vision is an understanding of how cells respond to environmental variables. Our findings indicate that the periplasmic iron binding protein sfuC is uniquely suitable for assessing the iron limitation status of Ca. Pelagibacter ubique cells. We anticipate that ecologists will use this data for interpreting the nutritional status of Ca. Pelagibacter ubique cells in nature.

This study, one of the few to simultaneously examine both transcriptional and translational responses in a bacteria cell, uncovered evidence suggesting that cspL might play a role in the cellular response to iron limitation. We offer the model that this protein controls translation in response to environmental conditions for a specific subset of genes present in the transcriptome. We hypothesize that this activity might serve an emergency function, limiting the synthesis of proteins to those that are critical for survival. This finding is consistent with previous reports of post-transcriptional regulation of the iron stress response in which a protein was found to facilitate the degradation of specific mRNAs which encoded nonessential iron-consuming pathways [26], [27].

Not only is Ca. Pelagibacter ubique one of the most successful cells known, it is also one of the simplest, giving it value as a model for understanding bacterial cell responses. Indeed, numerous new structural RNAs, some widely distributed among bacteria, have been discovered and described in Ca. Pelagibacter ubique [37], [38]. It is perhaps hubris to imagine that the concept of systems biology might one day be extended from the machinery of cells to the machinery of microbial ecosystems at work on the scale of oceans. But, if that vision has a chance, it will be by combining studies that cross scales and disciplines to understand the keystone species of the oceans.

Materials and Methods

Growth Media and Harvesting

Seawater was collected on 6/14/08 at the Newport Hydroline station NH5 (44°39.1′N, 124°10.6′W) from a depth of 10 m. The water was then filtered through a 0.2 µM filter, autoclaved, and sparged with CO2 for 24 hours followed by air for 24–48 hours as previously described [68], [69]. Immediately prior to inoculation with Ca. Pelagibacter ubique HTCC1062, the media was amended with 50 µM pyruvate, 50 µM glucose, 10 µM nitrogen, 1 µM methionine, 1 µM glycine, 1 µM phosphate, and vitamins. Cells were grown at 20°C (flasks) or 16°C (carboys) with intermittent light and sparging with air. On day 16, three 20 L control cultures were amended with 100 nM ferrichrome and 1 µM FeCl3, and three 20 L treatment cultures were amended with 100 nM ferrichrome only. On day 18, 8 L from each carboy was harvested. On day 28, the remaining ~10L from each was harvested. Prior to each harvest, and on day 17, three 40 mL samples of culture were removed from each culture for microarrays. Water from the three replicate cultures were then combined and growth was arrested using 0.01g chloramphenical and 0.1 mL protease inhibitor cocktail II (CalBiochem #539132) per liter of culture. Tangential flow filtration, followed by centrifugation produced cell pellets for the mass-spectrometry analysis. All samples were kept at −80°C until analysis.

Messenger RNA Preparation

Ca. Pelagibacter ubique cells used in microarray experiments were grown in batch cultures as described above. Cells (40 ml for each biological replicate) were collected via centrifugation, and RNA was extracted using RNeasy Mini kits (Qiagen), followed by amplification with MessageAmp-II Bacteria RNA amplification kit (Ambion). The resulting aRNA was then screened for length and quality using a Bioanalyzer 2100 (Agilent) and quantified utilizing a Nanodrop 1000 spectrophotometer (Thermo Fisher Scientific). 5.5 µg of biotinylated aRNA from each sample was then fractionated and hybridized (45°C) overnight to custom Ca. Pelagibacter ubique Affymetrix GeneChip arrays that contained probes for strains HTCC1002, HTCC1062 and HTCC7211 (Pubiquea520471f) using Affymetrix GeneChip Fluidics Station 450, and Affymetrix GeneChip Hybridization Oven 640. Arrays were then washed as per the manufacturer's instructions and the resulting images were analyzed using an Affymetrix GeneChip Scanner 3000. Fluorescence measurements were normalized over all 18 microarray chips.

Microarray Clustering

A modified radial coordinate visualization plot was used for illustrating mRNA expression in a manner that accentuated condition-specific preferential transcription. In Figure 2, dimensional anchors (DA) representing each of the six microarray samples were positioned manually around the circumference of a circle such that iron-limited samples are on the left, samples with excess iron are on the right, and the vertical placement corresponds to the culture's transition from exponential growth (bottom) to stationary phase (top). Each gene is represented by a single point, positioned according to the relative abundances between every sample pair, and sized according to the largest observed change in expression level. PTg is the point for gene g, with attributes x, y, and s describing its x-axis position, y-axis position, and size, respectively. Si,g is the log base-10 average fluorescence for gene g in sample i. DAi is sample i's dimensional anchor positioned at (DAi,x, DAi,y) on the graph.

equation image
equation image
equation image

This type of graph is ideal for revealing if a given gene's transcript abundance is changing as a result of iron limitation or as a result of the stationary phase transcriptome remodeling induced by iron limitation.

Global TFE Protein Preparation

Four samples were prepared using the TFE (2,2,2-Trifluoroethanol) digestion method. The cell pellets were reconstituted in 100 mM NH4HCO3, pH 8.4 buffer and transferred to a siliconized 0.6 mL microcentrifuge tube. 0.1 mm Zirconia/Silica Beads were added to the top of the tube and bead beat at maximum speed for 3 minutes and immediately placed on ice. A hole was poked in the base of the 0.6 mL siliconized eppendorf tube and placed in a 1.5 mL siliconized eppendorf tube. The sample was then centrifuged for 5 minutes at 14,000 rpm at 4°C. The cell lysis was mixed to a homogenized state and the volume was determined using a pipette. The sample concentration was determined with a Coomassie protein assay and read on a microplate reader. TFE was added to a concentration of 50%. The sample was then homogenized by sonication for one minute in an ice bath followed by incubation at 60°C for two hours with gentle shaking (300 rpm). Proteins were reduced by adding DTT to a final concentration of 2 mM, sonicated for one minute in an ice bath and incubated at 37°C for one hour with gentle shaking. Samples were then diluted 5-fold with 100 mM NH4HCO3 to reduce the salt concentration, and CaCl2 was added to a final concentration of 1 mM. The sample was digested for 3 hours with Trypsin (Promega, Madison WI) at 37°C at a concentration of 1 unit trypsin/50 units protein. After trypsin incubation, a BCA protein assay was performed on the sample to determine the final concentration and vialed for mass spectrometer analysis.

Capillary LC-MS Analysis

The custom HPLC system was configured using 65-mL Isco Model 65D syringe pumps (Isco, Inc., Lincoln, NE), 2-position Valco valves (Valco Instruments Co., Houston, TX), and a PAL autosampler (Leap Technologies, Carrboro, NC), allowing for fully automated sample analysis across four separate HPLC columns. Reversed-phase capillary HPLC columns were manufactured in-house by slurry packing 3-µm Jupiter C18 stationary phase (Phenomenex, Torrence, CA) into a 70-cm length of 360 µm o.d.×75 µm i.d. fused silica capillary tubing (Polymicro Technologies Inc., Phoenix, AZ) that incorporated a 0.5-µm retaining screen in a 1/16” custom laser-bored 75 µm i.d. union (screen and union – Valco Instruments Co., Houston, TX; laser bore - Lenox Laser, Glen Arm, MD). Mobile phases consisted of 0.2% acetic acid and 0.05% TFA in water (A) and 0.1% TFA in 90% acetonitrile/10% water (B). The mobile phase flowed through an in-line Degassex DG4400 degasser (Phenomenex, Torrance, CA). The HPLC system was equilibrated at 10 k psi with 100% mobile phase A. Fifty minutes after sample injection the mobile phase was switched to 100% B, which created a near-exponential gradient as mobile phase B displaced A in a 2.5 mL active mixer. A 30-cm length of 360 µm o.d.×15 µm i.d. fused silica tubing was used to split ~20 µL/min of flow before it reached the injection valve (5 µL sample loop). The split flow controlled the gradient speed under conditions of constant pressure operation (10 k psi). Flow through the capillary HPLC column when equilibrated to 100% mobile phase A was ~400 nL/min.

MS analysis was performed using a ThermoFinnigan LTQ-Orbitrap mass spectrometer (Thermo Scientific, San Jose, CA) with electrospray ionization (ESI). The HPLC columns were coupled to the mass spectrometer by using an in-house manufactured interface. Chemically etched electrospray emitters, 150 um o.d.×20 um i.d, were used [70]. The heated capillary temperature and spray voltage were 200°C and 2.2 kV, respectively. Data was acquired for 100 min, beginning 65 min after sample injection (15 min into gradient). Orbitrap spectra (AGC 1×106) were collected from 400–2000 m/z at a resolution of 100k followed by data dependant ion trap MS/MS spectra (AGC 1×104) of the six most abundant ions using a collision energy of 35%. A dynamic exclusion time of 60 sec was used to discriminate against previously analyzed ions. Three technical replicates were run on the mass spectrometer for each cell pellet.

Quantitative Proteomics

Quantitative estimates of peptide abundances, calculated from the area under the isotopic profile, were obtained by using a previously developed accurate mass and time (AMT) tag library [34] to search the mass spectra generated by the 12 runs for the four samples. After deisotoping and calculating monoisotopic mass, mass spectrometric features were matched to database peptides with a mass tolerance window of +/−6ppm and an elution time window of +/−0.1% after alignment in both dimensions. Peptide abundances were reported for those which had observations in at least 2 of the 3 technical replicates. Linear regression normalization was used to normalize each set of technical replicates as described elsewhere [71]. Briefly, the abundance of peptide x in sample i was transformed into minus versus average space using the following formulas:

equation image
equation image

Next, the transformed value was corrected based on a linear regression:

equation image

where mi* is the value for mi calculated from the m vs a regression equation. Lastly, the computed values were deconvoluted to yield the normalized abundances:

equation image

Peptides were excluded from further analysis if the standard deviation exceeded the average measurement value among the three technical replicates for a sample. A final filter was applied to exclude the lowest third of peptides for a given protein, when sorted by the peptides' maximum PeptideProphet F-Score.

Protein abundance was calculated only if a protein had three or more peptides which passed the above filters. Calculating the difference in protein abundance between two samples was a three step process. First, the three replicate peptide abundance measurements were averaged together. Next, the peptide average from sample 1 was divided by the peptide average from sample 2, then log10 transformed. Finally, all log10 peptide ratios from the same protein were averaged together.

To represent the likelihood that a protein was equally abundant in both samples, the multiple peptide measurements were combined into a single statistic as previously described [72]. Briefly, p-values for individual peptides were calculated using a one-tailed Student's t-test on the technical replicates' x′ values. A two-tailed Student's t-test was not used because p-values reflecting a large increase would be indistinguishable from p-values reflecting a large decrease. Instead, peptides which changed in the opposite direction from the protein average were assigned a p-value of 1 for their one-tailed Student's t-test. All peptide p-values for a single protein were then combined into a single chi-square statistic using Fisher's method:

equation image

Supporting Information

Table S1

Coordinates of All Genes Plotted in Figure 2.

(0.25 MB XLS)

Acknowledgments

The authors would like to thank the crew of the Elaka for their help in seawater collection, Zanna Chase for measuring the iron concentration in that water, and Mark Wells for advice on siderophore selection and protocols. Input from Paul Carini greatly helped in the design of Figure 3.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by Marine Microbiology Initiative award number 607_01 from the Gordon and Betty Moore Foundation (http://www.moore.org/). Portions of this research were performed at the W.R. Wiley Environmental Molecular Science Laboratory, a national scientific user facility sponsored by the U.S. Department of Energy's Office of Biological and Environmental Research, located at Pacific Northwest National Laboratory (http://www.pnl.gov/). PNNL is operated by Battelle Memorial Institute under DOE contract DE-AC05-76RL01830. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Martin JH, Gordon RM, Fitzwater SE. Iron in Antarctic waters. Nature. 1990;345:156–158. doi: 10.1038/345156a0.
2. Coale KH, Johnson KS, Fitzwater SE, Gordon RM, Tanner S, et al. A massive phytoplankton bloom induced by an ecosystem-scale iron fertilization experiment in the equatorial Pacific Ocean. Nature. 1996;383:495–501. [PubMed]
3. Boyd PW, Watson AJ, Law CS, Abraham ER, Trull T, et al. A mesoscale phytoplankton bloom in the polar Southern Ocean stimulated by iron fertilization. Nature. 2000;407:695–702. [PubMed]
4. Fung IY, Meyn SK, Tegen I, Doney SC, John JG, et al. Iron supply and demand in the upper ocean. Global Biogeochem Cycles. 2000;14:281–295.
5. Andrews SC, Robinson AK, Rodriguez-Quinones F. Bacterial iron homeostasis. FEMS Microbiol Rev. 2003;27:215–237. doi: 10.1016/S0168-6445(03)00055-X. [PubMed]
6. Johnson KS, Gordon RM, Coale KH. What controls dissolved iron concentrations in the world ocean? Mar Chem. 1997;57:137–161.
7. Martin JH, Coale KH, Johnson KS, Fitzwater SE, Gordon RM, et al. Testing the iron hypothesis in ecosystems of the equatorial Pacific Ocean. Nature. 1994;371:123–129. doi: 10.1038/371123a0.
8. Escolar L, Pérez-Martín J, de Lorenzo V. Opening the iron box: transcriptional metalloregulation by the Fur protein. J Bacteriol. 1999;181:6223–6229. [PMC free article] [PubMed]
9. Ducey TF, Carson MB, Orvis J, Stintzi AP, Dyer DW. Identification of the iron-responsive genes of Neisseria gonorrhoeae by microarray analysis in defined medium. J Bacteriol. 2005;187:4865–4874. doi: 10.1128/JB.187.14.4865-4874.2005. [PMC free article] [PubMed]
10. Baichoo N, Helmann JD. Recognition of DNA by Fur: a reinterpretation of the Fur box consensus sequence. J Bacteriol. 2002;184:5826–5832. doi: 10.1128/JB.184.21.5826-5832.2002. [PMC free article] [PubMed]
11. Fuangthong M, Helmann JD. Recognition of DNA by three ferric uptake regulator (Fur) homologs in Bacillus subtilis. J Bacteriol. 2003;185:6348–6357. [PMC free article] [PubMed]
12. Grifantini R, Sebastian S, Frigimelica E, Draghi M, Bartolini E, et al. Identification of iron-activated and -repressed Fur-dependent genes by transcriptome analysis of Neisseria meningitidis group B. Proc Natl Acad Sci USA. 2003;100:9542–9547. doi: 10.1073/pnas.1033001100. [PubMed]
13. Lavrrar JL, McIntosh MA. Architecture of a fur binding site: a comparative analysis. J Bacteriol. 2003;185:2194–2202. [PMC free article] [PubMed]
14. Tiss A, Barre O, Michaud-Soret I, Forest E. Characterization of the DNA-binding site in the ferric uptake regulator protein from Escherichia coli by UV crosslinking and mass spectrometry. FEBS Lett. 2005;579:5454–5460. doi: 10.1016/j.febslet.2005.08.067. [PubMed]
15. Yang J, Sangwan I, Lindemann A, Hauser F, Hennecke H, et al. Bradyrhizobium japonicum senses iron through the status of haem to regulate iron homeostasis and metabolism. Mol Microbiol. 2006;60:427–437. doi: 10.1111/j.1365-2958.2006.05101.x. [PMC free article] [PubMed]
16. Rudolph G, Semini G, Hauser F, Lindemann A, Friberg M, et al. The Iron control element, acting in positive and negative control of iron-regulated Bradyrhizobium japonicum genes, is a target for the Irr protein. J Bacteriol. 2006;188:733–744. doi: 10.1128/JB.188.2.733-744.2006. [PMC free article] [PubMed]
17. Martínez M, Ugalde RA, Almirón M. Dimeric Brucella abortus Irr protein controls its own expression and binds haem. Microbiology (Reading, Engl.) 2005;151:3427–3433. doi: 10.1099/mic.0.28213-0. [PubMed]
18. Yang J, Ishimori K, O'Brian MR. Two heme binding sites are involved in the regulated degradation of the bacterial iron response regulator (Irr) protein. J Biol Chem. 2005;280:7671–7676. doi: 10.1074/jbc.M411664200. [PubMed]
19. Friedman YE, O'Brian MR. A novel DNA-binding site for the ferric uptake regulator (Fur) protein from Bradyrhizobium japonicum. J Biol Chem. 2003;278:38395–38401. doi: 10.1074/jbc.M306710200. [PubMed]
20. Hamza I, Chauhan S, Hassett R, O'Brian MR. The bacterial irr protein is required for coordination of heme biosynthesis with iron availability. J Biol Chem. 1998;273:21669–21674. [PubMed]
21. Rodionov DA, Gelfand MS, Todd JD, Curson ARJ, Johnston AWB. Computational reconstruction of iron- and manganese-responsive transcriptional networks in alpha-proteobacteria. PLoS Comput Biol. 2006;2:e163. doi: 10.1371/journal.pcbi.0020163. [PubMed]
22. Massé E, Arguin M. Ironing out the problem: new mechanisms of iron homeostasis. Trends Biochem Sci. 2005;30:462–468. doi: 10.1016/j.tibs.2005.06.005. [PubMed]
23. Massé E, Salvail H, Desnoyers G, Arguin M. Small RNAs controlling iron metabolism. Curr Opin Microbiol. 2007;10:140–145. doi: 10.1016/j.mib.2007.03.013. [PubMed]
24. Gaballa A, Antelmann H, Aguilar C, Khakh SK, Song K, et al. The Bacillus subtilis iron-sparing response is mediated by a Fur-regulated small RNA and three small, basic proteins. Proc Natl Acad Sci USA. 2008;105:11927–11932. doi: 10.1073/pnas.0711752105. [PubMed]
25. Vecerek B, Moll I, Bläsi U. Control of Fur synthesis by the non-coding RNA RyhB and iron-responsive decoding. EMBO J. 2007;26:965–975. doi: 10.1038/sj.emboj.7601553. [PubMed]
26. Vergara SV, Thiele DJ. Post-transcriptional regulation of gene expression in response to iron deficiency: co-ordinated metabolic reprogramming by yeast mRNA-binding proteins. Biochem Soc Trans. 2008;36:1088–1090. doi: 10.1042/BST0361088. [PubMed]
27. Dubrac S, Touati D. Fur positive regulation of iron superoxide dismutase in Escherichia coli: functional analysis of the sodB promoter. J Bacteriol. 2000;182:3802–3808. [PMC free article] [PubMed]
28. Vraspir JM, Butler A. Chemistry of Marine Ligands and Siderophores. Annu Rev Mar Sci. 2009;1:43–63. doi: 10.1146/annurev.marine.010908.163712. [PMC free article] [PubMed]
29. Venturi V, Ottevanger C, Bracke M, Weisbeek P. Iron regulation of siderophore biosynthesis and transport in Pseudomonas putida WCS358: involvement of a transcriptional activator and of the Fur protein. Mol Microbiol. 1995;15:1081–1093. [PubMed]
30. Wilhelm SW, Maxwell DP, Trick CG. Growth, iron requirements, and siderophore production in iron-limited Synechococcus PCC 7002. Limnol Oceanogr. 1996;41:89–97.
31. Liu X, Theil EC. Ferritins: dynamic management of biological iron and oxygen chemistry. Acc Chem Res. 2005;38:167–175. doi: 10.1021/ar0302336. [PubMed]
32. Hintze KJ, Theil EC. Cellular regulation and molecular interactions of the ferritins. Cell Mol Life Sci. 2006;63:591–600. doi: 10.1007/s00018-005-5285-y. [PubMed]
33. Giovannoni SJ, Tripp HJ, Givan S, Podar M, Vergin KL, et al. Genome streamlining in a cosmopolitan oceanic bacterium. Science. 2005;309:1242–1245. doi: 10.1126/science.1114057. [PubMed]
34. Sowell SM, Norbeck AD, Lipton MS, Nicora CD, Callister SJ, et al. Proteomic analysis of stationary phase in the marine bacterium “Candidatus Pelagibacter ubique”. Appl Environ Microbiol. 2008;74:4091–4100. doi: 10.1128/AEM.00599-08. [PMC free article] [PubMed]
35. Sowell SM, Wilhelm LJ, Norbeck AD, Lipton MS, Nicora CD, et al. Transport functions dominate the SAR11 metaproteome at low-nutrient extremes in the Sargasso Sea. ISME J. 2009;3:93–105. doi: 10.1038/ismej.2008.83. [PubMed]
36. Schwalbach MS, Tripp HJ, Steindler L, Smith DP, Giovannoni SJ. The presence of the glycolysis operon in SAR11 genomes is positively correlated with ocean productivity. Environ Microbiol. 2009 Available: http://www.ncbi.nlm.nih.gov/pubmed/19889000. [PubMed]
37. Tripp HJ, Schwalbach MS, Meyer MM, Kitner JB, Breaker RR, et al. Unique glycine-activated riboswitch linked to glycine-serine auxotrophy in SAR11. Environ Microbiol. 2009;11:230–238. doi: 10.1111/j.1462-2920.2008.01758.x. [PMC free article] [PubMed]
38. Meyer MM, Ames TD, Smith DP, Weinberg Z, Schwalbach MS, et al. Identification of candidate structured RNAs in the marine organism ‘Candidatus Pelagibacter ubique’. BMC Genomics. 2009;10:268. doi: 10.1186/1471-2164-10-268. [PMC free article] [PubMed]
39. Corbino KA, Barrick JE, Lim J, Welz R, Tucker BJ, et al. Evidence for a second class of S-adenosylmethionine riboswitches and other regulatory RNA motifs in alpha-proteobacteria. Genome Biol. 2005;6:R70. doi: 10.1186/gb-2005-6-8-r70. [PMC free article] [PubMed]
40. Worden AZ, Lee J, Mock T, Rouzé P, Simmons MP, et al. Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes Micromonas. Science. 2009;324:268–272. doi: 10.1126/science.1167222. [PubMed]
41. Winkler W, Nahvi A, Breaker RR. Thiamine derivatives bind messenger RNAs directly to regulate bacterial gene expression. Nature. 2002;419:952–956. doi: 10.1038/nature01145. [PubMed]
42. Morris RM, Rappé MS, Connon SA, Vergin KL, Siebold WA, et al. SAR11 clade dominates ocean surface bacterioplankton communities. Nature. 2002;420:806–810. doi: 10.1038/nature01240. [PubMed]
43. Malmstrom RR, Kiene RP, Cottrell MT, Kirchman DL. Contribution of SAR11 bacteria to dissolved dimethylsulfoniopropionate and amino acid uptake in the North Atlantic ocean. Appl Environ Microbiol. 2004;70:4129–4135. doi: 10.1128/AEM.70.7.4129-4135.2004. [PMC free article] [PubMed]
44. Malmstrom RR, Cottrell MT, Elifantz H, Kirchman DL. Biomass production and assimilation of dissolved organic matter by SAR11 bacteria in the Northwest Atlantic Ocean. Appl Environ Microbiol. 2005;71:2979–2986. doi: 10.1128/AEM.71.6.2979-2986.2005. [PMC free article] [PubMed]
45. Chen G, Gharib TG, Huang C, Taylor JMG, Misek DE, et al. Discordant protein and mRNA expression in lung adenocarcinomas. Mol Cell Proteomics. 2002;1:304–313. [PubMed]
46. Griffin TJ, Gygi SP, Ideker T, Rist B, Eng J, et al. Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol Cell Proteomics. 2002;1:323–333. [PubMed]
47. Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, et al. Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science. 2001;292:929–934. doi: 10.1126/science.292.5518.929. [PubMed]
48. de Sousa Abreu R, Penalva LO, Marcotte EM, Vogel C. Global signatures of protein and mRNA expression levels. Mol Biosyst. 2009;5:1512–1526. doi: 10.1039/b908315d. [PubMed]
49. Nie L, Wu G, Culley DE, Scholten JCM, Zhang W. Integrative analysis of transcriptomic and proteomic data: challenges, solutions and applications. Crit Rev Biotechnol. 2007;27:63–75. doi: 10.1080/07388550701334212. [PubMed]
50. Wu G, Nie L, Zhang W. Integrative analyses of posttranscriptional regulation in the yeast Saccharomyces cerevisiae using transcriptomic and proteomic data. Curr Microbiol. 2008;57:18–22. doi: 10.1007/s00284-008-9145-5. [PubMed]
51. Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, et al. Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004;14:2308–2318. doi: 10.1101/gr.2523904. [PubMed]
52. Kudla G, Murray AW, Tollervey D, Plotkin JB. Coding-sequence determinants of gene expression in Escherichia coli. Science. 2009;324:255–258. doi: 10.1126/science.1170160. [PubMed]
53. Ermolenko DN, Makhatadze GI. Bacterial cold-shock proteins. Cell Mol Life Sci. 2002;59:1902–1913. [PubMed]
54. Bae W, Xia B, Inouye M, Severinov K. Escherichia coli CspA-family RNA chaperones are transcription antiterminators. Proc Natl Acad Sci USA. 2000;97:7784–7789. [PubMed]
55. Stülke J. Control of transcription termination in bacteria by RNA-binding proteins that modulate RNA structures. Arch Microbiol. 2002;177:433–440. doi: 10.1007/s00203-002-0407-5. [PubMed]
56. Graumann PL, Marahiel MA. A superfamily of proteins that contain the cold-shock domain. Trends Biochem Sci. 1998;23:286–290. [PubMed]
57. Henkin TM. Transcription termination control in bacteria. Curr Opin Microbiol. 2000;3:149–153. [PubMed]
58. Hofweber R, Horn G, Langmann T, Balbach J, Kremer W, et al. The influence of cold shock proteins on transcription and translation studied in cell-free model systems. FEBS J. 2005;272:4691–4702. doi: 10.1111/j.1742-4658.2005.04885.x. [PubMed]
59. Phadtare S, Severinov K, Inouye M. Assay of transcription antitermination by proteins of the CspA family. Meth Enzymol. 2003;371:460–471. doi: 10.1016/S0076-6879(03)71034-9. [PubMed]
60. Hu KH, Liu E, Dean K, Gingras M, DeGraff W, et al. Overproduction of three genes leads to camphor resistance and chromosome condensation in Escherichia coli. Genetics. 1996;143:1521–1532. [PubMed]
61. Xia B, Ke H, Inouye M. Acquirement of cold sensitivity by quadruple deletion of the cspA family and its suppression by PNPase S1 domain in Escherichia coli. Mol Microbiol. 2001;40:179–188. [PubMed]
62. Phadtare S, Inouye M. Role of CspC and CspE in regulation of expression of RpoS and UspA, the stress response proteins in Escherichia coli. J Bacteriol. 2001;183:1205–1214. doi: 10.1128/JB.183.4.1205-1214.2001. [PMC free article] [PubMed]
63. Phadtare S, Inouye M. Sequence-selective interactions with RNA by CspB, CspC and CspE, members of the CspA family of Escherichia coli. Mol Microbiol. 1999;33:1004–1014. [PubMed]
64. Geissmann T, Marzi S, Romby P. The role of mRNA structure in translational control in bacteria. RNA Biol. 2009;6:153–160. [PubMed]
65. Marzi S, Fechter P, Chevalier C, Romby P, Geissmann T. RNA switches regulate initiation of translation in bacteria. Biol Chem. 2008;389:585–598. [PubMed]
66. Jickells TD, An ZS, Andersen KK, Baker AR, Bergametti G, et al. Global iron connections between desert dust, ocean biogeochemistry, and climate. Science. 2005;308:67–71. doi: 10.1126/science.1105959. [PubMed]
67. Goldman JC. Conceptual Role for Microaggregates in Pelagic Waters. Bulletin of Marine Science. 1984;35:462–476.
68. Connon SA, Giovannoni SJ. High-throughput methods for culturing microorganisms in very-low-nutrient media yield diverse new marine isolates. Appl Environ Microbiol. 2002;68:3878–3885. [PMC free article] [PubMed]
69. Rappé MS, Connon SA, Vergin KL, Giovannoni SJ. Cultivation of the ubiquitous SAR11 marine bacterioplankton clade. Nature. 2002;418:630–633. doi: 10.1038/nature00917. [PubMed]
70. Kelly RT, Page JS, Luo Q, Moore RJ, Orton DJ, et al. Chemically etched open tubular and monolithic emitters for nanoelectrospray ionization mass spectrometry. Anal Chem. 2006;78:7796–7801. doi: 10.1021/ac061133r. [PMC free article] [PubMed]
71. Callister SJ, Barry RC, Adkins JN, Johnson ET, Qian W, et al. Normalization approaches for removing systematic biases associated with mass spectrometry and label-free proteomics. J Proteome Res. 2006;5:277–286. doi: 10.1021/pr050300l. [PMC free article] [PubMed]
72. Hess A, Iyer H. Fisher's combined p-value for detecting differentially expressed genes using Affymetrix expression arrays. BMC Genomics. 2007;8:96. doi: 10.1186/1471-2164-8-96. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science