Transcriptional changes underlying elemental stoichiometry shifts in a marine heterotrophic bacterium.

Marine bacteria drive the biogeochemical processing of oceanic dissolved organic carbon (DOC), a 750-Tg C reservoir that is a critical component of the global C cycle. Catabolism of DOC is thought to be regulated by the biomass composition of heterotrophic bacteria, as cells maintain a C:N:P ratio of ∼50:10:1 during DOC processing. Yet a complicating factor in stoichiometry-based analyses is that bacteria can change the C:N:P ratio of their biomass in response to resource composition. We investigated the physiological mechanisms of resource-driven shifts in biomass stoichiometry in continuous cultures of the marine heterotrophic bacterium Ruegeria pomeroyi (a member of the Roseobacter clade) under four element limitation regimes (C, N, P, and S). Microarray analysis indicated that the bacterium scavenged for alternate sources of the scarce element when cells were C-, N-, or P-limited; reworked the ratios of biomolecules when C- and P- limited; and exerted tighter control over import/export and cytoplasmic pools when N-limited. Under S limitation, a scenario not existing naturally for surface ocean microbes, stress responses dominated transcriptional changes. Resource-driven changes in C:N ratios of up to 2.5-fold and in C:P ratios of up to sixfold were measured in R. pomeroyi biomass. These changes were best explained if the C and P content of the cells was flexible in the face of shifting resources but N content was not, achieved through the net balance of different transcriptional strategies. The cellular-level metabolic trade-offs that govern biomass stoichiometry in R. pomeroyi may have implications for global carbon cycling if extendable to other heterotrophic bacteria. Strong homeostatic responses to N limitation by marine bacteria would intensify competition with autotrophs. Modification of cellular inventories in C- and P-limited heterotrophs would vary the elemental ratio of particulate organic matter sequestered in the deep ocean.

Introduction

Bacterioplankton control the flux of dissolved organic carbon (DOC) into the microbial food web and influence the release of carbon to atmospheric, offshore, and deep sea reservoirs. Many studies suggest that catabolism of DOC is regulated by the biomass stoichiometry of heterotrophic bacteria with respect to N and other nutrients (del Giorgio and Cole, 1998). Thus bacteria incorporate a higher percentage of metabolized DOC into biomass when C is the limiting element and a lower percentage when other elements (typically N or P) limit growth, resulting in decreased growth efficiency as the substrate C:N ratio increases (Goldman et al., 1987). While the molar stoichiometry of C:N:P in marine seston averages ∼106:16:1 (the Redfield ratio; Redfield, 1934), bacteria from a variety of freshwater and marine environments have higher N and P requirements relative to C and typically must attain a biomass stoichiometry closer to 50:10:1 during DOC processing (Fagerbakke et al., 1996; Cotner et al., 2010).

A confounding factor in stoichiometry-based analyses is that heterotrophic bacteria can change the C:N:P ratio of their biomass in response to the substrate composition (Martinussen and Thingstad, 1987; Tezuka, 1990; Fagerbakke et al., 1996; Gundersen et al., 2002). This plasticity occurs either by uncoupling catabolism and anabolism to overproduce metabolites rich in the excess element (which can be stored or excreted), or by reducing the requirement for the limiting element. Resource-driven shifts in C:N:P ratios of bacterial biomass can affect the efficiency of DOC transformation, modify C flux through the marine food web (Elser et al., 1995), buffer mismatches between bacterial requirements and ecosystem resource availability, and alter the composition of particulate organic matter sequestered in the deep ocean (Thingstad et al., 2008).

While it is clear that any modifications in elemental stoichiometry must occur within the basic macromolecular constraints of a functioning heterotrophic bacterial cell, there is yet only rudimentary knowledge of the physiological mechanisms that might allow C:N:P ratio variability in marine bacteria. Proposed avenues for changes in stoichiometry include alteration of the cellular inventory of P-rich rRNA (Elser et al., 1996), trade-offs in macromolecular composition to maximize growth (Franklin et al., 2011), substitution of S for P in bacterial membrane lipids (Van Mooy et al., 2009), depletion of P reserves stored in the form of polyphosphate (Kornberg, 1999), depletion of C reserves stored in the form of polyhydroxybutyrate (Dawes and Senior, 1973), and reduction of the cell quota of the limiting element (i.e., the required cell concentration; Harder and Dijkhuizen, 1983; Rivkin and Anderson, 1997).

Here we investigate the global transcripn>tional response underlying resource-driven shifts in biomass stoichiometry in the marine heterotrophic bacterium Ruegeria pomeroyi DSS-3, a member of the ubiquitous Roseobacter clade, under four different element limitation regimes (C, N, P, and S). Our experimental approach equalized growth rates across treatments using continuous cultures (Ferenci, 2008), thereby disentangling effects of the limiting element from any growth rate-driven differences in cell composition (Elser et al., 1995). R. pomeroyi exhibited significant flexibility in biomass content of some elements but not others, and gene transcription patterns deduced from whole-genome microarrays revealed element-specific metabolic strategies underlying the stoichiometric shifts.

Materials and Methods

Cell growth conditions

Ruegeria pomeroyi DSS-3 was grown in 320 ml custom-made chemostats with a culture volume of 200 ml. Steady-state biomass was limited either by carbon (glucose, 1 mmol l−1), nitrogen (NH4Cl, 0.26 mmol l−1), phosphorus (KH2PO4, 9.2 μmol l−1), or sulfur (Na2SO4, 25 μmol l−1). The appropriate concentrations of limiting nutrients needed to produce similar biomass were initially approximated from batch culture experiments that measured growth yields of R. pomeroyi under a range of concentrations. The range was then narrowed in test runs of the chemostats. Once the correct element concentration was determined, six independent chemostat cultures for each macroelement limitation run at two different times were used for elemental analysis and triplicate RNA sampling.

Chemostat culture medium (Table A1 in Appendix) modified from Henriksen (2008) was buffered with 10 mmol l−1 1,3-bis(tris(hydroxymethyl)methyl amino)propane (Bis–tris propane), a C- and N-containing compound; batch cultures of R. pomeroyi did not grow with 10 mmol l−1 Bis–tris propane when it was added as the sole N or C source. Vitamins and trace metals were added to the culture medium (Table A1 in Appendix). The feed medium was added to each chemostat at a rate of 8.4 ml h−1, equivalent to a dilution rate of 0.042 h−1, chosen to approach the rate of marine bacterial growth in situ (∼1–2 day−1; Ducklow and Hill, 1985) while maintaining sufficient cell yield for biological and chemical analysis. During the incubation, cell cultures were mixed by constant stirring, and temperature was maintained at 30°C using a circulating water bath. Air was bubbled into the culture at a flow rate of 2 ml min−1. Salinity was constant at 25 and pH ranged from 6.6 to 7.0 during all growth regimes. Cell cultures were considered to be at steady-state when the change in OD600 was ≤10% between two successive measurements. At steady-state, the feed medium dilution rate equals cell growth rate, which was equivalent to a 16.5-h doubling time for R. pomeroyi DSS-3. All cultures reached steady-state after three volume exchanges, and cells were harvested after five volume exchanges

Table A1

Chemostat culture medium.

Compound	Concentration
MACRONUTRIENTa
Glucose (C6H12O6)	4.50 mmol l−1 (1.00 mmol l−1)
Phosphate (KH2PO4)	0.50 mmol l−1 (9.20 μmol l−1)
Ammonium (NH4Cl)	2.80 mmol l−1 (0.26 mmol l−1)
Sulfate (Na2SO4)	2.50 mmol l−1 (25.00 μmol l−1)
BUFFER, IRON, AND SALTb
Bis–tris propane (C11H26N2O6)c	9.91 mmol l−1
CaCl2·2H2O	7.42 mmol l−1
MgCl2	106.21 mmol l−1
Iron-EDTA (C10H12FeN2NaO8)d	67.52 μmol l−1
KCl	10.64 mmol l−1
NaCl	198.49 mmol l−1
TRACE ELEMENT
C6H9NO6	12.95 μmol l−1
CoCl2·6H2O	0.80 μmol l−1
Na2SeO3	1.87 μmol l−1
Na2WO4·2H2O	0.55 μmol l−1
BaCl2·2H2O	0.49 μmol l−1
CuSO4	55.57 nmol l−1
MnSO4·H2O	0.77 μmol l−1
ZnSO4·7H2O	0.54 μmol l−1
Na2MoO4·2H2O	0.57 μmol l−1
Na2SiO3·9H2O	0.45 μmol l−1
SrCl2·6H2O	0.19 mmol l−1
NiCl2·6H2O	96.32 nmol l−1
VITAMIN
Thiamine (B1)	0.15 μmol l−1
Nicotinic acid (B3)	0.40 μmol l−1
Pyridoxine-HCl (B6)	0.48 μmol l−1
Folic acid (B9)	44.86 nmol l−1
Cyanocobalamin (B12)	7.30 nmol l−1
Riboflavin (B2)	0.13 μmol l−1
Pantothenic acid (B5)	0.23 μmol l−1
Biotin (B7)	81.05 nmol l−1
p-aminobenzoic acid (B10)	0.36 μmol l−1

.

.

.

.

Inoculum for the chemostat cultures was from a frozen glycerol stock of R. pomeroyi DSS-3, revived by streaking onto 1/2 strength YTSS agar (yeast extract, 4 g l−1; tryptone, 2.5 g l−1; sea salts, 20 g l−1; agar, 20 g l−1; González et al., 2003; pH adjusted to 6.8), and incubating at 30°C in the dark. After colonies appeared (∼3 days after streaking), a patch of cells was transferred into 1/2 YTSS liquid medium and grown at 30°C with shaking at 150 rpm. After an overnight incubation, cells were washed two times and resuspended in fresh medium lacking glucose, ammonium, phosphate, and sulfate (Table A1 in Appendix). Cells were inoculated into 200 ml chemostat medium to an OD600 of 0.05 (∼7 × 106 cells ml−1), and cultured initially with the outflow pump turned off. After ∼16 h, the flow carrying the feed medium was started and the medium supply was maintained at a constant rate during the chemostat run.

Cell counts and protein measurements

Chemostat cultures were routinely monitored for OD600, cell counts, and cell protein concentration. For cell counts, aliquots taken from the chemostats were fixed with glutaraldehyde (2% vol/vol, final concentration), incubated at room temperature for 10 min, and frozen at −20°C. One freeze–thaw cycle did not cause a change in the cell count (data not shown). Thawed cells were stained with SYBR Green II (Molecular Probes/Invitrogen, Carlsbad, CA, USA), incubated in the dark for 10 min, and analyzed with a Beckman Colter Cyan Flow Cytometer. Bacterial cells were quantified using a combination of forward light scatter and fluorescence detection. To reduce the contribution of background fluorescence, cells were quantified within a gate encompassing ≥94% of total fluorescence. Light scatter values provided an index of cell size. The cell protein concentration was measured from cells flash-frozen in a dry ice–ethanol bath according to the method of Bradford (1976) using Bovine Serum Albumin as a standard.

Spent medium analysis

For each treatment regime, concentrations of glucose, and were measured in one sample (one or two technical duplicates) of cell-free spent medium pooled from three independent chemostat cultures. Cells were pelleted by centrifugation at 16,250 × g. The supernatant from three replicates was combined and passed through a 0.22-μm pore-size cellulose acetate filter to remove remaining cells and stored at −20°C until analysis. An estimate of ammonium in the supernatant was obtained by assaying the spent medium using the Ammonia Assay Kit (BioVision, Mountain View, CA, USA; detection limit of 20 μmol l−1) according to the manufacturer’s instructions and correcting for the interference from pyruvate based on a control prepared in parallel with the sample. Attempts to measure ammonium with the more sensitive phenol–hypochlorite method (Solorzano, 1969) failed due to inhibiting agents in the culture medium. Glucose concentration was measured with the Glucose Assay Kit (Sigma, St. Louis, MO, USA; detection limit of 50 μmol l−1), and control samples were prepared in parallel with the supplied glucose standard. The concentration of was determined with a Shimadzu UV-1601 spectrophotometer using the molybdenum blue reaction (Strickland and Parsons, 1972) and the concentration of was determined with a DIONEX ICS-2000 ion chromatography system. Internal standards were used in and measurements.

Growth efficiency calculations

Carbon assimilatory efficiency was calculated as Cb/Cipt, where Cb is the C present in cell biomass per chemostat volume (200 ml), calculated as 3,072, 1,007, 2,928, and 511 μg C for the C-, N-, P-, and S-limited conditions, respectively, based on measures of biomass %C and dry weight; and Cipt is the C input to the chemostat per cell doubling minus any C remaining in the spent medium, calculated as 9,988 μg C for C-limited conditions (8.4 ml h−1 of a 1-mmol l−1 solution) and 44,946 μg C minus spent medium C for N-, P-, and S-limited conditions (8.4 ml h−1 of a 4.5-mmol l−1 solution).

Cell elemental analysis

Ruegeria pomeroyi DSS-3 cells from 240 ml of culture, pooled from three independent chemostat replicates (80 ml per replicate), were pelleted by centrifugation at 16,250 × g, washed twice in fresh medium without glucose, ammonium, phosphate, and sulfate (Table A1 in Appendix), immediately frozen in a dry ice–ethanol bath, and stored at −70°C. Frozen cell samples were lyophilized for dry weight measurements and then used in elemental analysis with a Perkin-Elmer 2400 CHN analyzer (C and N) and an inductively coupled plasma mass spectrometer (P).

Free amino acid analysis

Cells from 180 ml of culture pooled from three independent chemostat replicates (60 ml per replicate) were washed and pelleted as described above. A portion was saved for cell counts by flow cytometry. To the remainder, 9.7 nmol of β-amino isobutyric acid was added as an internal standard (Hendrickson et al., 2008) and cells were rapidly frozen in a dry ice–ethanol bath and stored at −80°C until analysis. Frozen cell pellets were sent to the Molecular Structure Facility at the University of California Davis for assay with a Hitachi L-8900 Amino Acid Analyzer according to the facility’s protocol. Briefly, cell pellets were lysed with acetonitrile (50% vol/vol) and formic acid (5% vol/vol), and amino-ethyl cysteine was added as an internal standard for quality control. The final amount of amino acids was corrected for the standard and converted to units of pmoles (mg dry weight)−1.

Microarray hybridization

Cell culture from each replicate chemostat (45 ml) was rapidly fixed with 5 ml of phenol–ethanol solution (5% vol/vol). Cells were collected by centrifugation at 4,500 × g, the supernatant was removed, and the loose cell pellets were frozen at −70°C. Total RNA was purified using the RNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to procedures described previously (Poretsky et al., 2009). DNA was removed by the TURBO DNA-free kit (Applied Biosystems/Ambion, Austin, TX, USA), ribosomal RNA was removed with the MicrobeExpress kit (Applied Biosystems/Ambion), and mRNA-enriched RNA was amplified with the MessageAmp II-Bacteria Kit (Applied Biosystems/Ambion).

A non-compn>etitive, single-color hybridization strategy was used to hybridize the mRNA samples to microarray slides. Amplified mRNA was labeled with AlexaFluor 647 dye (Invitrogen), with modifications to the manufacturer’s protocol as described previously (Bürgmann et al., 2007). Amplified and labeled mRNA was purified with the MEGAclear kit (Applied Biosystems/Ambion) and concentrated by ethanol precipitation. Samples were hybridized to a whole-genome microarray manufactured by Combimatrix (Mukilteo, WA, USA) and imaged according to the procedures described in the manufacturer’s protocol. The CustomArray 12 K microarray slides for R. pomeroyi DSS-3 used in this study were described in detail by Bürgmann et al. (2007). A total of 8,143 probes, each present up to three times, were represented on the array (Bürgmann et al., 2007). Hybridized slides were stripped according to Combimatrix’s protocol and re-used up to three times. Arrays were re-imaged after each stripping, and re-stripped as needed until minimal fluorescence signal remained.

Microarray data analysis

Quality control of microarray hybridizations was performed as described previously (Bürgmann et al., 2007). Each hybridized spot was manually examined in GenePix Pro 6.0 Software (Molecular Devices, Sunnyvale, CA, USA), and the background fluorescence was subtracted (i.e., fluorescence values of the five closest empty spots). Spots were flagged as bad with Acuity 4.0 (Molecular Devices) when they were unevenly hybridized (determined manually by checking the hybridization circularity), when the detected fluorescence minus the background fluorescence was less than 2 SDs of the background fluorescence [F635–B635 < 2 SD (B635)], or when the signal to noise ratio was <3. Probes were removed when the flagged spots were present in >25% of arrays. For the remaining probes with good hybridization signal, those targeting 5S rRNA, 16S rRNA, 23S rRNA, mismatch probes, and manufacturer’s control probes were removed. Each gene was represented by up to two probes and each probe was spotted up to three times on the array (Bürgmann et al., 2007). After probe removal, 5,145 probes (covering 3,324 genes out of the total 4,252 genes in the R. pomeroyi DSS-3 genome) remained for analysis. The complete microarray dataset, reported according to the Minimum Information About a Microarray Experiment guidelines (Brazma et al., 2001), was appended to the Gene Expression Omnibus (GEO; The National Center for Biotechnology Information), platform GPL4067, under series accession number GSE27032.

Florescence values from 14 microarray hybridizations (three biological replicates from each macroelement-limited condition and one technical replicate each for the N and S limitations) were normalized globally, as in Rinta-Kanto et al. (2010). Data for the technical replicates were averaged and combined (Figure A1 in Appendix). Spearman rank correlation coefficients between replicate arrays ranged between 0.82 and 0.96, while those between arrays from different treatments were ≤0.80 (Table A2 in Appendix). The exception was between P and S limitation arrays, for which the correlation coefficients ranged between 0.71 and 0.85 (Table A2 in Appendix). Statistical analysis was performed in MultiExperiment Viewer (MeV, Dana-Farber Cancer Institute, Boston, MA, USA; Saeed et al., 2006). Significance analysis of microarray (SAM; Tusher et al., 2001), available through the MeV interface, was used for statistical analysis. Pairwise comparisons were made between treatments, with the following settings: two-class unpaired, 5,000 permutations, and false discovery rate (FDR; Tusher et al., 2001) ≤2%. Pairwise SAM analyses were performed as follows: C limitation compared to C excess conditions (where N, P, or S was limiting); N limitation compared to N excess conditions (where C, P, or S was limiting); P limitation compared to P excess conditions (where N, C, or S was limiting); and S limitation compared to S excess conditions (where N, C, or P were limiting; Figure A1 in Appendix).

Figure A1

Schematic of microarray data analysis. Probe counts are shown in italicized font and gene counts are in bold font. (A) Significant Analysis of Microarray (SAM) was used to analyze the microarray data (pairwise SAM comparisons with 5,000 permutations and a false discovery rate ≤ 2%). From these results, probes significantly different in all three comparisons to the other treatments (that is, falling into the three-way overlap region of each Venn diagram) were selected for further analyses. (B) From the normalized fluorescence values, the fold-change for each probe was calculated as the median value between all limitation-to-excess pairwise comparisons and genes whose fold-change was ≥ 3 were retained.

Table A2

Spearman’s rank correlation coefficients of the normalized fluorescence values in 5,145 probes between replicate arrays.

C-limited		N-limiteda			P-limited			S-limiteda
0.91	0.92	0.80	0.78	0.79	0.74	0.77	0.70	0.76	0.76	0.74	C-limited
	0.96	0.70	0.75	0.74	0.67	0.66	0.64	0.73	0.69	0.69
		0.72	0.78	0.79	0.66	0.69	0.63	0.71	0.71	0.69
			0.91	0.88	0.76	0.77	0.69	0.67	0.69	0.73	N-limiteda
				0.92	0.75	0.75	0.69	0.63	0.69	0.70
					0.68	0.75	0.60	0.58	0.69	0.63
						0.91	0.96	0.75	0.81	0.82	P-limited
							0.89	0.71	0.85	0.79
								0.76	0.81	0.84
									0.82	0.90	S-limiteda
										0.87

.

Transcripn>tionally responsive genes were identified as those with a ≥3-fold enrichment in all three pairwise comparisons across treatments (Figure A1 in Appendix). Four genes (SPO0968, SPO1031, SPO1795, and SPO3300) were transcriptionally enriched under one treatment but transcriptionally depleted under another, due to different responses of two probes for a single gene; these were not considered further. Designation of enriched or depleted genes was conservative. In many cases genes had significant responses in two out of three pairwise comparisons or had a ≥2-fold (but ≤3-fold) transcriptional change, but they were not considered further.

Gene annotation and metabolic pathway construction

Ruegeria pomeroyi DSS-3 gene products were generally assigned according to the annotations in GenBank (Accession number NC_003911.11) with updated gene calls as appropriate (SPO0781, phnD; SPO1948, pstS; SPO1979, OlsA-like protein; SPO1980, OlsB-like protein; SPO0372 and SPO2287, LuxI-type proteins). Operons were as predicted in Roseobase, and metabolic pathways were constructed manually from BioCyc (SRI International, Menlo Park, CA, USA). For analysis of the global ocean sampling (GOS) metagenome (Rusch et al., 2007; Yooseph et al., 2007; data obtained from Community Cyberinfrastructure for Advanced Microbial Ecology Research and Analysis), BLASTp was used to identify homologous sequences using an E-value of <10−5. The percent of bacterioplankton cells in GOS samples (as of May 2011; Sargasso Sea, northwest Atlantic, Pacific, and Indian Ocean sequences; Rusch et al., 2007; Yooseph et al., 2007) harboring a homolog was calculated as (number of homologs × 100)/number of recA, assuming only one homolog and one recA per genome and 10,196 recA homologs in the analyzed portion of the GOS dataset (Moran et al., 2011). For analysis of other microbial genomes, BLASTp was used to identify homologs using an E-value < 10−30. In both cases, candidate homologs were subjected to phylogenetic analysis using pplacer (Matsen et al., 2010) to confirm homology with genes of known function.

Cell stoichiometry

The stoichiometry rules that best reproduced observed C:N:P ratios under elemental limitation were explored. A simple model that assumed limitation shifted the concentration of the limiting element away from the concentration found under balanced growth (when all elements are available in non-limiting conditions) by a factor α, while the other two (non-limiting) elements were unaffected by the limitation gave the best fit to the observed C:N:P ratios. Thus the C:N:P ratio for R. pomeroyi under balanced growth conditions was approximated as:

where αX is the factor by which element x changes when it is limiting, and C, N, and P are the C, N, and P content, respectively, under limiting element x.

Results and Discussion

Attributes of steady-state cultures

Measurements of cell-free spent media indicated that limiting elements were below the detection limit while non-limiting elements were in excess for chemostat-grown R. pomeroyi DSS-3 (Table 1). Cell dry weights for the C, N, and S treatments averaged 486 fg cell−1, similar to those reported for heterotrophic marine bacteria during exponential growth (322–512 fg cell−1, Vrede et al., 2002; Table A3 in Appendix), but P-limited cells had an unusually high dry weight (1,388 fg cell−1). C-limited cells were smaller and had higher growth efficiency (i.e., amount of new bacterial biomass produced per unit of organic C assimilated) than N-, P-, or S- limited cells (31 versus 1–8%; Table 1). By way of comparison, del Giorgio and Cole (1998) report 9–47% bacterial growth efficiencies (interquartile range) for heterotrophic bacteria in oceans and estuaries. The protein: dry weight ratio of S-limited cells was atypically high, likely reflecting an unnatural condition for marine bacterioplankton. Molar C:N ratios in biomass varied up to 2.5-fold and C:P ratios varied up to sixfold among the different limiting element regimes (Table 1).

Table 1

Biological and chemical status of .

Limitation	Spent medium (mmol l−1)				Steady-state culture
	Glucosea	NH4+a	PO43-b	SO42-b	Protein (fg cell−1)c	Dry weight (fg cell−1)b	Protein: dry weight (%)	Flow scatter (mode)c	C growth efficiency	Molar ratiob
										C:N	C:P	N:P
C	BDd	1.02	0.40	1.76	140.8	454.0	31.0	1,792	30.8	4.4	55.6	12.7
N	0.72	BD	0.46	1.52	353.9	757.9	46.7	2,048	2.7	11.2	166.7	14.8
P	0.67	0.54	BD	1.81	371.5	1388.4	26.8	2,048	7.7	9.5	333.3	35.0
S	0.66	NDe	0.43	0.01	234.6	245.8	95.4	2,099 ± 280	1.3	4.9	62.5	12.7

.

.

.

.

.

Table A3

Bacterial dry weights.

Organism	Growth conditionsa	Dry weight (fg cell−1)	Reference
Ruegeria pomeroyi	Steady-state culture, C-limited (dt = 16.5 h)	454	This study
	Steady-state culture, N-limited (dt = 16.5 h)	758	This study
	Steady-state culture, P-limited (dt = 16.5 h)	1388	This study
	Steady-state culture, S-limited (dt = 16.5 h)	246	This study
Escherichia coli	Batch culture in low potassium medium (growing)	710b	Fagerbakke et al. (1996)
	Batch culture in low potassium medium (stationary)	180b	Fagerbakke et al. (1996)
	Batch culture in glucose minimal medium (dt = 40 min)	280c	Neidhardt et al. (1996)
Vibrio natriegens	Batch culture in brain heart infusion medium (growing)	850b	Fagerbakke et al. (1996)
	Batch culture in brain heart infusion medium (stationary)	145b	Fagerbakke et al. (1996)
Natural marine bacterioplankton isolate	Batch culture (exponential)	427d	Vrede et al. (2002)
	Batch culture, C-limited (stationary)	110d	Vrede et al. (2002)
	Batch culture, N-limited (stationary)	276d	Vrede et al. (2002)
	Batch culture, P-limited (stationary)	270d	Vrede et al. (2002)
	Native isolates from various aquatic environments	21–60e	Fagerbakke et al. (1996)

.

.

.

.

.

.

.

.

.

.

Transcripn>ts for 190 genes (4.5% of the 4,252 genes in the R. pomeroyi DSS-3 genome) were responsive exclusively to C, N, P, or S limitation, with 134 genes enriched (3.2%) and 56 genes depleted (1.3%; Figure A1 in Appendix; see Tables A4 and A5 in Appendix for a complete list of significantly enriched and depleted genes).

Table A4

Transcriptionally enriched .

Locus tag	Gene	Product	Fold-change ratioa
			Probe 1	Probe 2
C-LIMITED (68 GENES)
SPO0090	–	Hypothetical protein	5.3	–
SPO0107	recQ	ATP-dependent DNA helicase	1.0	6.3
SPO0120	–	Hypothetical protein	3.9	–
SPO0149	dnaA	Chromosomal replication initiation protein	16.5	1.2
SPO0171	fliR	Flagellar biosynthetic protein	1.3	6.7
SPO0183	fliI	H+-transporting two-sector ATPase, flagellum-specific	6.3	1.0
SPO0265	–	Auxin efflux carrier family protein	2.1	4.4
SPO0276	–	LuxR family transcriptional regulator	0.8	3.0
SPO0590	–	LacI family transcription regulator (carbon catabolite repression domain-containing)	14.5	–
SPO0875	gap-2	Glyceraldehyde-3-phosphate dehydrogenase, type I	8.4	1.3
SPO0911	proC	Pyrroline-5-carboxylate reductase	3.2	3.2
SPO0939	–	M48 family peptidase	1.1	8.5
SPO0940	–	Hypothetical protein	3.4	–
SPO0970	apqZ	Aquaporin Z	3.1	1.3
SPO0988	–	YeeE/YedE family protein	1.1	3.0
SPO0999	soxD	Diheme cytochrome c	1.4	3.8
SPO1001	soxF	Sulfur oxidation F protein	3.2	–
SPO1063	–	Hypothetical protein	1.2	3.3
SPO1092	–	LysR family transcriptional regulator	3.0	1.8
SPO1132	–	Glycine betaine/proline ABC transporter, ATP-binding protein	3.1	3.0
SPO1138	–	AsnC family transcriptional regulator	6.2	4.7
SPO1156	hisF	Imidazole glycerol phosphate synthase subunit	3.0	3.5
SPO1337	–	Hypothetical protein	3.9	–
SPO1446	–	Cyclase, putative	3.8	–
SPO1463	–	TRAP dicarboxylate transporter, DctM subunit	8.3	–
SPO1499	–	Alpha/beta fold family hydrolase	5.2	1.4
SPO1592	–	Aminomethyl transferase family protein	2.3	3.4
SPO1645	–	Oligopeptide/dipeptide ABC transporter, permease protein	4.6	–
SPO1702	–	Hypothetical protein	3.6	–
SPO1774	–	3-hydroxyanthranilate 3,4-dioxygenase	3.3	–
SPO1796	–	Formate dehydrogenase, alpha subunit, putative	3.4	4.7
SPO1807	–	DNA polymerase III epsilon subunit family exonuclease	3.7	1.0
SPO2029	–	Glutamine amidotransferase class-II	3.4	–
SPO2166	–	Putative lipoprotein	0.6	3.3
SPO2205	–	ErfK/YbiS/YcfS/YnhG family protein	7.2	0.5
SPO2301	–	OsmC-like family protein	5.2	–
SPO2306	–	Hypothetical protein	3.4	–
SPO2325	serS	Seryl-tRNA synthetase	3.2	–
SPO2370	–	Sodium:alanine symporter family protein	3.4	–
SPO2371	–	Universal stress protein family protein	4.4	1.8
SPO2422	–	d-isomer specific 2-hydroxyacid dehydrogenase family protein	6.4	–
SPO2644	–	VWA domain CoxE-like family protein	11.3	–
SPO2647	–	Hypothetical protein	3.0	–
SPO2668	–	LysR family transcriptional regulator	3.2	1.1
SPO2747	–	Diguanylate cyclase, putative	3.7	–
SPO2845	–	Hypothetical protein	3.8	–
SPO2881	–	Xanthine dehydrogenase family protein, large subunit	5.0	1.0
SPO2989	–	Cytochrome b562	3.1	1.4
SPO3121	–	MATE efflux family protein	4.3	0.9
SPO3291	–	Branched-chain amino acid ABC transporter, periplasmic	5.5	4.0
SPO3407	–	Hypothetical protein	1.7	5.7
SPO3461	–	Flagellar protein FlgJ, putative	8.5	–
SPO3552	–	Oxidoreductase, FAD-binding	3.2	–
SPO3562	tauR	Transcriptional regulator	5.2	1.5
SPO3577	–	Hypothetical protein	3.8	–
SPO3589	–	Hypothetical protein	3.0	0.8
SPO3641	–	ABC transporter, permease protein	1.2	3.7
SPO3734	–	MerR family transcriptional regulator	4.2	–
SPO3735	–	Major facilitator family protein	1.4	6.2
SPO3737	–	Pyridine nucleotide-disulfide oxidoreductase family protein	4.4	–
SPO3764	–	Glutathione S-transferase family protein	1.2	6.8
SPO3783	–	Sugar ABC transporter, ATP-binding protein	8.0	–
SPOA0097	–	Branched-chain amino acid ABC transporter, permease protein	3.6	1.4
SPOA0160	–	TRAP dicarboxylate transporter, DctM subunit	8.7	–
SPOA0247	–	Hypothetical protein	3.0	–
SPOA0289	–	AraC family transcriptional regulator	4.1	–
SPOA0299	–	Branched-chain amino acid ABC transporter, permease protein	0.9	3.7
SPOA0317	–	AsnC family transcriptional regulator	10.1	3.0
N-LIMITED (35 GENES)
SPO0021	–	Hpt domain-containing protein	3.3	–
SPO0161	–	LuxR family DNA-binding response regulator	3.1	0.8
SPO0296	–	Hypothetical protein	4.6	3.0
SPO0322	–	Hypothetical protein	3.9	6.5
SPO0327	–	EAL domain-containing protein	0.4	4.9
SPO0438	–	ErfK/YbiS/YcfS/YnhG family protein	1.5	3.0
SPO0804	ihfB	Integration host factor subunit beta	3.0	1.4
SPO0860	–	Xylose repressor, putative	3.2	3.2
SPO0919	–	MarR family transcriptional regulator	3.3	2.8
SPO1078	–	Chromosome replication initiation inhibitor protein	3.0	–
SPO1275	–	Cold shock family protein	2.5	3.5
SPO1286	–	Hypothetical protein	0.5	3.7
SPO1384	–	MarR family transcriptional regulator	4.4	1.4
SPO1707	–	Urea transporter, ATP-binding protein (urtE)	3.8	3.0
SPOA0447b	–	Urea transporter, ATP-binding protein (urtD)	5.5	–
SPO1708	–	Urea transporter, permease protein (urtC)	11.6	–
SPO1808	–	Hypothetical protein	4.6	–
SPO1853	–	TetR family transcriptional regulator	3.0	–
SPO1872	–	LysR family transcriptional regulator	6.2	2.1
SPO1911	–	Hypothetical protein	1.1	4.3
SPO2067	–	Hypothetical protein	6.5	–
SPO2087	ntrC	Nitrogen regulation protein	1.1	4.3
SPO2224	–	Hypothetical protein	2.2	4.1
SPO2294	glnB-1	Nitrogen regulatory protein PII	9.0	6.6
SPO2295	glnA	Glutamine synthetase, type I	4.7	5.1
SPO2364	–	Amino acid ABC transporter, periplasmic binding protein	1.3	4.1
SPO2602	–	RpiR family transcriptional regulator	5.1	2.2
SPO2778	–	Hypothetical protein	3.0	1.5
SPO2818	–	Hypothetical protein	4.9	3.0
SPO3607	–	Hypothetical protein	4.6	0.6
SPO3723	amt-2	Ammonium transporter	28.0	–
SPO3724	glnB-2	Nitrogen regulatory protein PII	124.3	39.1
SPOA0280	–	TRAP dicarboxylate transporter, DctP subunit, putative	7.9	4.5
SPOA0300	–	Branched-chain amino acid ABC transporter, periplasmic binding protein	6.6	–
SPOA0399	–	R body protein RebB-like protein	5.6	–
P-LIMITED (20 GENES)
SPO0147	–	Enoyl-CoA hydratase	4.3	3.0
SPO0304	–	Putative lipoprotein	7.0	5.4
SPO0468	phnG	Alkylphosphonate utilization protein	5.2	–
SPO0472	phnK	Phosphonate C-P lyase system protein	10.9	–
SPO0505	rplO	Ribosomal protein L15	3.2	5.7
SPO0781	phnD	Phosphonate ABC transporter, periplasmic binding protein	103.5	–
SPO1227	–	Hypothetical protein	3.2	–
SPO1504	pqqA	Coenzyme PQQ biosynthesis protein A	5.1	4.7
SPO1860	–	Twin-arginine translocation pathway signal sequence domain-containing protein	7.4	–
SPO1928	–	Tat pathway signal sequence domain-containing protein	3.5	0.7
SPO1948	pstS	Phosphate ABC transporter, periplasmic binding protein	47.9	–
SPO1949	pstC	Phosphate ABC transporter, permease protein	4.7	–
SPO1951	pstB	Phosphate transporter ATP-binding protein	7.4	–
SPO1953	phoB	Phosphate regulon transcriptional regulatory protein	10.7	26.1
SPO2626	–	TRAP transporter, DctM subunit	1.8	4.2
SPO2627	–	TRAP transporter, DctQ subunit	3.3	–
SPO3198	rnc	Ribonuclease III	4.1	1.6
SPO3625	cspA	Cold shock protein	2.6	4.2
SPO3868	–	Hypothetical protein	1.6	3.6
SPOA0294	pmtA	Phosphatidylethanolamine N-methyltransferase	3.6	1.2
S-LIMITED (11 GENES)
SPO0412	–	Hypothetical protein	13.2	–
SPO0636	–	EF hand domain-containing protein	3.9	–
SPO1256	ppk2	Polyphosphate kinase 2	5.0	4.0
SPO1330	hflC	HflC protein	4.0	3.5
SPO1409	rpoH-2	RNA polymerase factor sigma-32	12.5	–
SPO2596	hemA-1	5-amino-levulinate synthase	1.8	4.6
SPO2632	cobA-1	Uroporphyrin-III C-methyltransferase	3.4	5.9
SPO2634	–	Sulfite reductase, putative	5.7	5.3
SPO3383	–	Thiol-specific antioxidant protein	4.9	4.5
SPO3527	–	Universal stress protein family protein	3.7	5.3
SPO3532	hemN	Coproporphyrinogen III oxidase	7.3	–

.

.

Light gray shading indicates genes included in Table 2.

Table A5

Transcriptionally depleted .

Locus tag	Gene	Product	Fold-change ratioa
			Probe 1	Probe 2
C-LIMITED (12 GENES)
SPO0416	mreD	Rod shape-determining protein	0.3	–
SPO0956	–	Hypothetical protein	0.3	–
SPO1244	efp	Elongation factor P	0.3	0.3
SPO1293	phaP	Phasin	0.1	0.1
SPO1756	–	Polysaccharide biosynthesis/export protein	0.3	0.7
SPO1868	–	Glutaredoxin-related protein	0.4	0.3
SPO2135	–	Hypothetical protein	0.3	–
SPO3022	valS	Valyl-tRNA synthetase	0.3	–
SPO3253	rpsP	30S ribosomal protein S16	2.0	0.2
SPO3658	–	d-alanyl-d-alanine carboxypeptidase family protein	0.4	0.3
SPO3840	rpsO	30S ribosomal protein S15	0.2	0.2
SPOA0031	nqrD	Na+-translocating NADH-quinone reductase subunit D	0.5	0.3
N-LIMITED (34 GENES)
SPO0227	–	PaxA, putative	0.3	0.2
SPO0406	rpoH-2	RNA polymerase factor sigma-32	0.6	0.3
SPO0862	xylH	Xylose ABC transporter, permease protein	0.3	0.3
SPO1060	–	Hypothetical protein	0.1	–
SPO1079	–	Lysine exporter, putative	0.1	2.0
SPO1308	–	Hypothetical protein	0.3	–
SPO1336	ispZ	Intracellular septation protein A	0.6	0.3
SPO1685	–	tRNA synthetase, class I family protein	0.1	–
SPO1743	–	Glutamate dehydrogenase	0.2	0.2
SPO1747	–	Sarcosine oxidase gamma subunit family protein	1.2	0.2
SPO1954	–	LysR family transcriptional regulator	1.7	0.2
SPO2071	–	Bcr/CflA subfamily drug resistance transporter	0.3	0.3
SPO2271	fabF	3-oxoacyl-(acyl carrier protein) synthase II	0.5	0.3
SPO2373	–	Hypothetical protein	0.1	–
SPO2414	–	Mandelate racemase/muconate lactonizing enzyme family protein	0.3	–
SPO2556	–	Allantoate amidohydrolase	0.2	–
SPO2592	–	Beta-lactamase, putative	1.0	0.3
SPO2616	tgt	Queuine tRNA-ribosyltransferase	0.7	0.3
SPO3012	–	Inositol-1-monophosphatase, putative	0.1	0.2
SPO3156	–	l-threonine aldolase, low-specificity, putative	1.3	0.3
SPO3419	–	UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase	1.9	0.3
SPO3437	–	Mechanosensitive ion channel family protein	0.2	0.3
SPO3484	–	HSP20 family protein	0.5	0.2
SPO3712	asd	Aspartate-semialdehyde dehydrogenase	0.2	0.3
SPO3720	tyrB	Aromatic amino acid aminotransferase	0.3	0.9
SPO3731	–	Glycerophosphoryl diester phosphodiesterase, putative	0.3	–
SPO3752	–	Radical SAM domain-containing protein	0.2	0.7
SPOA0011	metK	S-adenosylmethionine synthetase	0.2	0.4
SPOA0057	gcvT	Glycine cleavage system T protein	0.2	0.3
SPOA0094	–	Hypothetical protein	0.2	–
SPOA0099	–	Branched-chain amino acid ABC transporter, ATP-binding protein	0.1	–
SPOA0115	gtdA-2	Gentisate 1,2-dioxygenase	0.8	0.1
SPOA0215	norQ	Nitric oxide reductase Q protein	2.2	0.1
SPOA0220	nirS	Cytochrome cd1 nitrite reductase	0.2	–
P-LIMITED (NO GENES)
–	–	–	–	–
S-LIMITED (10 GENES)
SPO0371	luxR-1	Autoinducer-binding transcriptional regulator	0.0	0.2
SPO0490	–	Calcium-binding domain–containing protein	0.3	0.3
SPO0491	–	Hypothetical protein	0.2	0.7
SPO1054	–	Hypothetical protein	0.3	–
SPO1221	–	Hypothetical protein	0.1	0.3
SPO1679	ctrA	DNA-binding response regulator	0.2	0.2
SPO2580	–	Hypothetical protein	0.1	0.1
SPO3223	–	Response regulator	0.3	0.2
SPO3673	–	Type I secretion target repeat-containing protein	0.3	0.2
SPOA0337	–	Hypothetical protein	0.1	–

.

Light gray shading indicates genes included in Table 2.

Transcriptional response to C limitation

The C-limited transcriptome was enriched in transporters for dicarboxylic acids, peptides, branched-chain amino acids, sugars, and glycine betaine/proline (Table 2), indicating increased investment in substrate acquisition. Several cell division transcripts were also present in higher proportion, including those for a chromosome initiation replication protein (dnaA), DNA helicase (recQ), DNA polymerase, and a purine biosynthesis enzyme (hisF; Table 2), coinciding with the smaller and more numerous cells observed under C limitation (Table 1). Transcripts for the rod shape-determining mreD, a mutated version of which causes Escherichia coli cells to take on a spherical shape (Wachi et al., 1989), were depleted (Table 2). Spherical shape and small cell size are characteristic of marine bacterioplankton in oligotrophic conditions, and are hypothesized to afford a more favorable surface-to-volume ratio to compete for scarce substrates (Azam et al., 1983). Relative increases in transcripts for three proteins involved in flagellum synthesis indicated increased motility of C-limited cells. Together, these changes represent a strong C scavenging response by C-limited R. pomeroyi cells.

Table 2

Selected .

Locus tag	Product (gene name, if available)	Fold-change ratioa		Response strategy
		Probe 1	Probe 2	Scavenging	Quota Change	Gate-keeping
C-LIMITED
SPO0107	ATP-dependent DNA helicase (recQ)	1.0	6.3	•
SPO0149	Chromosomal replication initiation protein (dnaA)	16.5	1.2	•
SPO0171	Flagellar biosynthetic protein	1.3	6.7	•
SPO0183	H+-transporting two-sector ATPase, flagellum-specific	6.3	1.0	•
SPO0590	LacI family transcriptional regulator (carbon catabolite repression domain-containing)	14.5	–	•
SPO0970	Aquaporin Z (apqZ)	3.1	1.3	•
SPO1132	Glycine betaine/proline ABC transporter, ATP-binding protein	3.1	3.0	•
SPO1156	Imidazole glycerol phosphate synthase subunit (hisF)	3.0	3.5	•
SPO1463	TRAP dicarboxylate transporter, DctM subunit	8.3	–	•
SPO1645	Oligopeptide/dipeptide ABC transporter, permease protein	4.6	–	•
SPO1807	DNA polymerase III epsilon subunit family exonuclease	3.7	1.0	•
SPO2370	Sodium:alanine symporter family protein	3.4	–	•
SPO3291	Branched-chain amino acid ABC transporter, periplasmic binding protein, putative	5.5	4.0	•
SPO3461	Flagellar protein FlgJ, putative	8.5	–	•
SPO3783	Sugar ABC transporter, ATP-binding protein	8.0	–	•
SPOA0097	Branched-chain amino acid ABC transporter, permease protein	3.6	1.4	•
SPOA0160	TRAP dicarboxylate transporter, DctM subunit	8.7	–	•
SPOA0299	Branched-chain amino acid ABC transporter, permease protein	0.9	3.7	•
SPO0416	Rod-shaped determining protein (mreD)	0.3	–	•
SPO1293	Phasin (phaP)	0.1	0.1		•
SPO1756	Polysaccharide biosynthesis/export protein	0.3	0.7		•
N-LIMITED
SPO0860	Xylose repressor, putative	3.2	3.2			•
SPO1707	Urea transporter, ATP-binding protein (urtE)	3.8	3.0	•
SPOA0447	Urea transporter, ATP-binding protein (urtD)	5.5	–	•
SPO1708	Urea transporter, permease protein (urtC)	11.6	–	•
SPO2087	Nitrogen regulation protein (ntrC)	1.1	4.3	•
SPO2294	Nitrogen regulatory protein PII (glnB-1)	9.0	6.6	•
SPO2295	Glutamine synthetase, type I (glnA)	4.7	5.1	•
SPO2364	Amino acid ABC transporter, periplasmic binding protein	1.3	4.1	•
SPO3723	Ammonium transporter (amt-2)	28.0	–	•
SPO3724	Nitrogen regulatory protein PII (glnB-2)	124.3	39.1	•
SPOA0300	Branched-chain amino acid ABC transporter, periplasmic binding protein	6.6	–	•
SPO0862	Xylose ABC transporter, permease protein (xylH)	0.3	0.3			•
SPO1079	Lysine exporter, putative	0.1	2.0			•
SPO1743	Glutamate dehydrogenase	0.2	0.2			•
SPO3156	l-threonine aldolase, low-specificity, putative	1.3	0.3
SPO3712	Aspartate-semialdehyde dehydrogenase (asd)	0.2	0.3			•
SPO3720	Aromatic amino acid aminotransferase (tyrB)	0.3	0.9			•
SPOA0011	S-adenosylmethionine synthetase (metK)	0.2	0.4			•
SPOA0057	Glycine cleavage system T protein (gcvT)	0.2	0.3			•
SPOA0215	Nitric oxide reductase Q protein (norQ)	2.2	0.1			•
SPOA0220	Cytochrome cd1 nitrite reductase (nirS)	0.2	–			•
P-LIMITED
SPO0147	Enoyl-CoA hydratase	4.3	3.0		•
SPO0304	Lipoprotein, putative	7.0	5.4		•
SPO0468	Alkylphosphonate utilization protein (phnG)	5.2	–	•
SPO0472	Phosphonate C-P lyase system protein (phnK)	10.9	–	•
SPO0781	Phosphonate ABC transporter, periplasmic phosphonate-binding (phnD)	103.5	–	•
SPO1860	Alkaline phosphatase (phoX)	7.4	–	•
SPO1948	Phosphate ABC transporter, periplasmic binding protein (pstS)	47.9	–	•
SPO1949	Phosphate ABC transporter, permease protein (pstC)	4.7	–	•
SPO1951	Phosphate transporter ATP-binding protein (pstB)	7.4	–	•
SPO1953	Phosphate regulon transcriptional regulatory protein (phoB)	10.7	26.1	•
SPOA0294	Phosphatidylethanolamine N-methyltransferase (pmtA)	3.6	1.2		•
S-LIMITED
SPO1256	Polyphosphate kinase 2 (ppk2)	5.0	4.0
SPO1409	RNA polymerase factor sigma-32 (rpoH-2)	12.5	–
SPO2596	5-amino-levulinate synthase (hemA-1)	1.8	4.6
SPO2632	Uroporphyrinogen-III C-methyltransferase (cobA-1)	3.4	5.9
SPO2634	Sulfite reductase, putative	5.7	5.3
SPO3383	Thiol-specific antioxidant protein	4.9	4.5
SPO3527	Universal stress protein family protein	3.7	5.3
SPO3532	Coproporphyrinogen III oxidase (hemN)	7.3	–
SPO0371	Autoinducer-binding transcriptional regulator (luxR-1)	0.0	0.2
SPO1679	DNA-binding response regulator (ctrA)	0.2	0.2

The enriched genes included here (unshaded rows) represent 36% of all enriched genes (Table .

.

Black dots indicate response strategy assignment.

Selepan class="Chemical">cted .

.

In the C-limited cultures, the R. pomeroyi transcriptome was depleted in transcripts encoding the synthesis and export of high C-content surface-associated polysaccharides, a strategy that would lower the cellular C quota (Table 2). Transcripts for phasin (phaP), a protein that binds to the surface of C storage inclusions composed of polyhydroxyalkanoate (PHA; Anderson and Dawes, 1990), were depleted (Table 2), indicating degradation of internal C reservoirs. Other genes involved in PHA metabolism (phbA, phbB, phaZ, phaR) did not have altered transcription levels (Figure A2 in Appendix), but this is consistent with data from the related alphaproteobacterium Rhodobacter capsulatus showing PHA synthesis genes to be constitutively expressed, with control occurring post-translationally (Kranz et al., 1997). Since PHA storage molecules have been observed in R. pomeroyi cells (González et al., 2003) and can occupy 40–50% of cell volume in some C-replete bacteria (Yurkov, 2006), this finding coincides with the smaller cell size observed under C limitation (Table 1). Together, this set of transcriptional changes comprise strategies for a reduced cell quota for C.

Figure A2

Biosynthesis of polyhydroxyalkanoate (PHA) in . (A) Pathway inferred from BioCyc; (B) arrangement of genes for PHA synthesis (phbAB, phaC), degradation (phaZ), and regulation (phaR) on R. pomeroyi DSS-3 chromosome; (C) transcriptional ratios for relevant genes; values are shown for genes listed in the pathway (solid red box) and genes that are not in the pathway, but are located proximal to PHA synthesis genes (dotted red box). Genes in the same operon are identically color-coded. *, transcriptional change is significant (≥3-fold, all limitation-to-excess pairwise comparisons).

Transcriptional response to N limitation

N limitation induced a strong scavenging response for acquiring available N compounds mediated through a conserved bacterial regulatory system involving PII regulators (glnB-1 and glnB-2) and the N response regulator NtrC (Ikeda et al., 1996; Gyaneshwar et al., 2005; Walter et al., 2007; Table 2). Two genes previously found to be controlled by these regulators in E. coli (Zimmer et al., 2000; Gyaneshwar et al., 2005) were also part of the R. pomeroyi scavenging response: an ammonium transporter (amt-2, located adjacent to glnB-2) and a glutamine synthetase (glnA, mediating ammonium incorporation; located adjacent to glnB-1). Two other N scavenging activities included enrichment of transcripts for urea transport and amino acid transport (Table 2).

There was no evidence of a role for N storage compounds in the R. pomeroyi response to N limitation. The genome does not contain homologs to the genes required for synthesis of cyanophycin granule polypeptide (CGP), the only known bacterial mechanism for N storage (Füser and Steinbüchel, 2007). Instead, transcriptional responses focused on tightening of the import/export balance of N through better “gate-keeping.” In one response, depletion of a xylose ABC transporter (xylH) and enrichment of a putative xylose uptake repressor (Table 2) would work to decrease the import of high C:N content carbohydrates. In another, depletion of transcripts for a putative lysine exporter (Table 2) involved in cytoplasmic amino acid regulation (Vrljic et al., 1999) would suppress the export of low C:N content amino acids. Two genes mediating dissimilatory nitrite reduction (norQ and nirS, which are part of a partial denitrification pathway; Moran et al., 2004) were also depleted (Table 2), suggesting another mechanism for N retention through a decrease in N-based respiration. The underrepresentation of transcripts for amino acid metabolism, including glutamate dehydrogenase, threonine aldolase, asd (aspartate-semialdehyde dehydrogenase), metK (S-adenosylmethionine synthetase), tyrB (aromatic amino acid aminotransferase), and gcvT (glycine cleavage system T protein; Table 2) was consistent with the substantial change in free amino acid concentrations observed under N-limited conditions (0.5% of cell N content compared to 1.3 and 4.2% under C and P limitation; Table A6 in Appendix), suggesting a reallocation of cytoplasmic N. Collectively, the transcriptional data support strong N homeostasis as the major response to N limitation, enabled through intense scavenging and tighter controls over N import/export and cytoplasmic pools.

Table A6

N content in free amino acids measured in steady-state .

Amino acida	N (ng) in free amino acid (mg dry weight)−1
	C-limited	N-limited	P-limited	S-limited
Alanine	17.1	34.5	53.6	665.6
Arginine	0.0	0.0	52.0	550.3
Glutamate	1405.6	123.3	2310.3	18142.8
Glycine	46.7	29.2	57.7	1196.5
Isoleucine	0.0	11.7	0.0	425.4
Leucine	5.0	13.1	22.7	219.4
Lysine	10.6	28.5	37.2	391.2
Phenylalanine	4.7	12.0	0.0	111.5
Serine	0.0	0.0	24.8	204.5
Threonine	3.2	14.9	17.7	686.4
Valine	4.5	14.1	19.6	395.6
Sum	1497.4	281.3	2595.7	22989.2
%Nb	1.3	0.5	4.2	25.5

.

.

Transcriptional response to P limitation

A phoB-mediated response involving scavenging for inorganic and organic sources of P dominated the P-limited transcriptome (Table 2). Transcripts were enriched for transport or utilization of phosphate (pstB, -C, -S; Yuan et al., 2006), phosphonate (phnD, -G, and -K; Metcalf and Wanner, 1991; Parker et al., 1999), and phosphate esters (phoX; Sebastian and Ammerman, 2011; Table 2). While inorganic P storage is common among bacteria, R. pomeroyi does not have a close homolog to the widespread P storage enzyme polyphosphate kinase I (ppk1; Newton et al., 2010) or to 1,3-diphosphoglycerate-polyphosphate phosphotransferase or polyphosphate:AMP-phosphotransferase (Kulaev et al., 1971; Ishige and Noguchi, 2000). It has been suggested that SPO0224 functions as a ppk1 in R. pomeroyi (Nahálka and Pätoprstý, 2009), but neither that gene nor the two ppk2 genes purported to degrade polyphosphates (SPOP1256 and SPO1727; Nahálka and Pätoprstý, 2009) were transcriptionally modified under P limitation. Instead, there was evidence for lowering the P cell quota by reworking of membrane phospholipids, including enrichment of transcripts for an enoyl-CoA hydratase mediating fatty acid degradation, a phosphatidyl ethanolamine N-methyltransferase (pmtA), and a lipoprotein likely representing a membrane-bound degradative enzyme or stress sensor (Table 2; Figure A3 in Appendix). R. pomeroyi possesses genes for the synthesis of ornithine-containing lipids (olsAB-like genes), which incorporate N rather than P into lipid membranes. The transcription of olsB was highest under P limitation (Figure A3 in Appendix), although the probes did not meet the significance cut off of ≥3-fold. These genes are important under P limitation in closely related Alphaproteobacteria (Weissenmayer et al., 2002; Gao et al., 2004) and in marine bacterioplankton communities (Van Mooy et al., 2009), and we leave open the possibility that this pathway is involved in P content reworking in R. pomeroyi.

Figure A3

Biosynthesis of phospholipid and ornithine lipid in . Strain DSS-3 does not encode genes to synthesize other bacterial P-free lipids (betaine lipid, btaAB; glycolipid/sulfolipid, sqdBCDX). (A) Pathway inferred from BioCyc; (B) arrangement of lipid synthesis genes on the chromosome; (C) transcription ratios for relevant genes; values shown for genes listed in (A,B) (solid red box). Genes in the same operon are identically color-coded. Symbols for ornithine lipid synthesis genes (yellow and orange) are bordered in black. *, transcriptional change is significant (≥3-fold, all limitation-to-excess pairwise comparisons). Abbreviation in (A,B) as follows: G3P, glyceraldehyde 3-phosphate.

Transcriptional response to S limitation

Unlike responses to C, N, and P limitation, none of the 21 transcriptionally responsive genes to S limitation were indicative of scavenging (Tables A4 and A5 in Appendix). Instead, response took the form of increased relative expression of enzymes requiring S-containing substrates or cofactors (Figure A4 in Appendix). Changes in transcript abundance for the synthesis of 5-amino-levulinate (hemA-1), protoporphyrinogen IX (hemN), sirohydrochlorin (cobA-1), and the quorum sensor regulator luxR-1 (Table 2; Figure A4 in Appendix) may all be linked to a depleted pool of S-adenosyl-methionine (Figure A4 in Appendix). Evidence for a cellular stress response was strongest under S limitation, and this included enrichment of transcripts for a thiol-containing antioxidant protein, a universal stress protein, and a σ32 factor (rpoH-2), along with the gene ppk2 (Table A4 in Appendix), which degrades polyphosphate to GTP and is implicated in bacterial stress responses (Brown and Kornberg, 2008; Gangaiah et al., 2010).

Figure A4

Relevant S-genes in . (A) Pathway inferred from BioCyc; (B) arrangement of relevant genes on the chromosome; (C) transcriptional ratios for relevant genes; ratios are shown for genes listed in the pathway (A) (solid red box) and genes that are not in the pathway, but are located in proximity to relevant S-genes [dotted red box in (B)]. Genes located in the same operon are identically color-coded. Two sets of orthologous LuxR–LuxI-type proteins are boxed in (B). Symbols for autoinducer-binding or synthesis genes (yellow, orange, brown) are bordered in black in (C). *Transcriptional change is significant (≥3-fold, all limitation-to-excess pairwise comparisons). Abbreviations in (A) as follows: THF, tetrahydrofolate; 5-mthglu, 5-methyltetrahydropteroyltri-l-glutamate; thglu, tetrahydropteroyltri-l-glutamate.

The S-limited transcripn>tome provided a useful perspective on transcriptional responses to limitation by an element that is not naturally limiting in the ocean (seawater SO4 concentrations are ∼28 mmol l−1) and therefore for which there is little selective pressure to evolve responses to scarcity. The biasing of the S-limited transcriptome toward biosynthetic pathways with depleted end products, while a successful strategy for addressing metabolic imbalances in bacteria (Goyal et al., 2010), cannot solve a cell-wide elemental deficit. In E. coli, scavenging responses occur through the CysB system that upregulates transporters for sulfate, cysteine, and alternate sources of S (Gyaneshwar et al., 2005). However, the R. pomeroyi DSS-3 genome does not contain a CysB ortholog.

The fact that S limitation was the only condition that did not invoke a scavenging response for obtaining more of the limiting element indicates that transcriptional regulation in R. pomeroyi is not evolutionarily attuned for S limitation. By analogy, the absence of an N storage response may indicate that conditions under which excess N is available for transport and storage by bacterioplankton are rare in the ocean. Indeed, a search of the marine bacterioplankton genes captured in the Global Ocean Sampling (GOS) dataset (Rusch et al., 2007; Yooseph et al., 2007) indicated that fewer than 3% of surface ocean bacterioplankton have homologs for cysB. Similarly, fewer than 1% of bacterioplankton carry homologs for either of the two genes required for N storage (CPG synthetase and cyanophycinase; Table A7 in Appendix). Transcriptional responses of R. pomeroyi may therefore provide insights into the role of seawater chemistry in molding the evolution of marine bacterioplankton genome content.

Table A7

Abundance of homologs for .

Gene function	Gene name	Number of homologsa	Cells with geneb (%)
S regulation	cysB	281	2.8
Cyanophycin granule synthesis	CPG synthetase	53	0.5
Cyanophycin granule synthesis	Cyanophycinase	24	0.2

.

.

Explaining cell stoichiometry changes

To link transcripn>tional responses to cell stoichiometry under different element limitations, we estimated the composition of R. pomeroyi biomass under balanced growth (i.e., when all elements are available in excess and the C:N:P ratio is optimal); this is the conceptual starting point from which the cells’ transcriptional responses produced the observed element-limited ratios (Table 1). The simplest model, one in which C:N:P ratios under element limitation are reached through a proportional change in only the limiting element, fit our data well. This does not necessarily require that C, N, and P are unlinked in the cell, but rather that the net result of resource limitation is a decrease in only the limiting element, where αX is the factor by which element x (c, n, or p) changes when it is limiting. The best solution to the elemental ratio data emerged with a C:N:P of ∼154:15:1 under balanced growth (Figure 1). From this initial ratio, C-limited cells (56:13:1; Table 1) would have undergone a relative C decrease to 38% of balanced growth levels (αC = 1/0.38), P-limited cells (333:35:1) would have undergone a relative P decrease to 40% of balanced growth (αP = 1/0.4), but N-limited cells (167:15:1) would have undergone almost no change in N content (αN = 1/1.0; Figure 1).

Figure 1

Modeled C:P and N:P atomic ratios of . The model assumes that limitation decreases the concentration of the limiting element relative to its value under balanced growth by a factor α, but the non-limiting elements are not affected. Measured C:N:P ratios in the element-limited chemostats are shown as symbols with bold borders, and modeled balanced growth ratios assuming different values of α are shown as solid lines. Dotted lines show elemental ratios if one of the non-limiting elements also changes during limitation, but by a factor 10-fold less than for the limiting element. The region where the three lines intersect (yellow shading) indicates the simultaneous solution for the C:N:P ratio (∼154:15:1). Other plotted C:P and N:P ratios are as follows: gray circle with black border, S-limited R. pomeroyi DSS-3 cells; purple X, Redfield ratio (Redfield, 1934); orange X, marine bacteria (Cotner et al., 1997); brown X, aquatic bacteria (Fagerbakke et al., 1996); blue X, freshwater bacteria (Cotner et al., 2010).

The transcripn>tional changes in R. pomeroyi strongly support this model. Only the C- and P-limited cells exhibited reworking strategies to decrease cell quotas (degradation of stored reserves, decreases in production or secretion of extracellular products, restructuring of membranes; Table 2). Only the N-limited cells showed a substantial gate-keeping response and tightening of control over cytoplasmic pools (Table 2). It is not yet evident whether strong N homeostasis compared to C and P is a feature that can be generalized to other heterotrophic marine bacterioplankton, although R. pomeroyi represents a ubiquitous bacterial taxon in surface ocean waters (Moran et al., 2004) and its growth rate in the continuous cultures was within the range measured for marine bacterioplankton in situ (Ducklow and Hill, 1985; Whitman et al., 1998). Four unidentified marine bacterial isolates studied by Vrede et al. (2002) did not show a similar N homeostasis (Figure A5 in Appendix), although they were not grown in continuous culture. Elemental ratios measured for natural bacterioplankton typically do not fall within the balanced growth region for R. pomeroyi (Figure 1) and are a better match to the ratios found under C limitation, suggesting that marine bacteria may be most often C-limited in the ocean.

Figure A5

C:N:P ratios from various studies. Black symbols, steady-state R. pomeroyi DSS-3; the unlimited C:N:P ratio was modeled as ∼154:15:1 (this study); blue, batch-cultured marine bacterioplankton (Vrede et al., 2002); purple, Redfield ratio, marine seston (Redfield, 1934); brown, natural bacteria (Fagerbakke et al., 1996); orange, natural bacteria (Cotner et al., 1997).

A central concept in stoichiometric ecology is that growth rate changes drive differences in biomass elemental composition through shifts in the content of P-rich RNA (Sterner et al., 2008). For marine microbes, this growth rate hypothesis (GRH) will be most relevant when bacterioplankton are growing at different rates, such as during phytoplankton blooms or mixing events relative to non-bloom conditions or dormancy. By comparison, R. pomeroyi was studied under constant growth rate but differing stoichiometry of external resources, a scenario more relevant when the element limiting heterotrophic growth varies over time or space (e.g., Thingstad, 1987; Kirchman, 1990; Pomeroy et al., 1995; Cherrier et al., 1996; Cotner et al., 1997; Kirchman and Rich, 1997; Zohary and Robarts, 1998; Torréton et al., 2000; Van Wambeke et al., 2002; Obernosterer et al., 2003; Thingstad et al., 2005). In exploring stoichiometric changes predicted under the GRH, Loladze and Elser (2011) calculated a microbial protein:rRNA ratio of 3 under optimal growth conditions (cellular N:P ratio = 16). Our balanced growth-stoichiometry calculation for R. pomeroyi agrees with this (N:P ratio = 15; Figure 1). Under non-optimal conditions however, which are commonplace in the ocean, transcriptional responses of heterotrophic bacteria can drive biomass C:N:P ratios from this ideal value through physiological mechanisms that store certain elements when available in excess and strongly retain others when scarce. The consistency of transcriptional reactions among the diverse taxa that make up marine bacterial assemblages will determine the aggregate community response to element limitation in the ocean (Hall et al., 2011; Scott et al., 2012).

Conclusion

Bacterioplankton transcriptional changes reveal the cellular basis for biomass stoichiometry outcomes whose effects are manifested at the ecosystem level. Flexible biomass ratios provide a buffer in the development of nutrient limitation (Thingstad et al., 2008) and, at least in the short term, delay intensification of bacterial–phytoplankton competition. Stoichiometric flexibility also broadens the range of C:N:P ratios in bacterial biomass produced in ocean surface waters, with consequences for basin-scale differences in particulate material and sequestration of organic matter deviating from Redfield proportions (Geider and La Roche, 2002; Hessen et al., 2004; Sterner et al., 2008; Weber and Deutsch, 2010). For model marine heterotrophic bacterium R. pomeroyi, cellular biomass ratios reflect a net transcriptional balance between element scavenging responses, gate-keeping activities, and cell quota changes, with the balance between these expression responses dictated by which element is limiting its growth.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.