Multiple nutrient stresses at intersecting Pacific Ocean biomes detected by protein biomarkers

See allHide authors and affiliations

Science  05 Sep 2014:
Vol. 345, Issue 6201, pp. 1173-1177
DOI: 10.1126/science.1256450

This article has a correction. Please see:

Protein markers of cyanobacterial stress

Nutrients including iron, nitrogen, and phosphorus limit primary productivity in the oceans. Determining how abundant cyanobacteria such as Prochlorococcus adapt to nutrient stress across marine settings requires accurate molecular assays. Saito et al. developed a proteomic and metaproteomic approach capable of targeting specific metabolic biomarkers from mixed communities in seawater (see the Perspective by Moore). Prochlorococcus proteins are indicative of a major limiting nutrient across a wide transect in the Pacific Ocean; however, they also show that the limitation of multiple nutrients at overlapping biomes is an additional source of stress.

Science, this issue p. 1173; see also p. 1120


Marine primary productivity is strongly influenced by the scarcity of required nutrients, yet our understanding of these nutrient limitations is informed by experimental observations with sparse geographical coverage and methodological limitations. We developed a quantitative proteomic method to directly assess nutrient stress in high-light ecotypes of the abundant cyanobacterium Prochlorococcus across a meridional transect in the central Pacific Ocean. Multiple peptide biomarkers detected widespread and overlapping regions of nutritional stress for nitrogen and phosphorus in the North Pacific Subtropical Gyre and iron in the equatorial Pacific. Quantitative protein analyses demonstrated simultaneous stress for these nutrients at biome interfaces. This application of proteomic biomarkers to diagnose ocean metabolism demonstrated Prochlorococcus actively and simultaneously deploying multiple biochemical strategies for low-nutrient conditions in the oceans.

Marine photosynthetic activity in the oceans is largely controlled by nutrient and micronutrient availability. Although this process is critical to marine ecosystem structure, ocean-climate interactions, and nutrient cycling processes, the nutrient-addition incubation experiments typically used to assess which nutrients are limiting are time-intensive and subject to artifacts (1), making their deployment to broad geographic regions difficult. Nutrient limitation is commonly parameterized in marine ecosystem and biogeochemistry models so that growth is controlled by the single scarcest nutrient (relative to cellular requirements) (2); however, multiple scarce nutrients may also influence phytoplankton community structure (3, 4). Molecular methodologies have shown potential for detecting in situ biomarkers for nutrient stress by measurement of up-regulated transcripts or proteins associated with nutrient scarcity (5, 6), but they can have relatively broad biological specificity and can be difficult to deploy to a wide geographic region (5). With advances in mass spectrometry, microbial proteins have begun to be directly measured in complex natural environments; recent metaproteomic surveys of oceanic environments have identified and examined the relative abundance of biogeochemically relevant proteins (7, 8).

In this study, we conducted calibrated quantitative mass spectrometry–based protein biomarker measurements using multiple reaction monitoring (MRM) to characterize nutrient limitation patterns for multiple nutrients on the abundant marine cyanobacterium Prochlorococcus (9). Samples were collected along a ~4500-km meridional transect through the central Pacific and equatorial Pacific Ocean initiating south of the Hawaiian Islands and ending near the Samoan Islands in October of 2011 (Fig. 1A). The equatorial Pacific is an ideal region to examine the range of nutrient stresses affecting phytoplankton because it crosses several biogeochemical biomes (10). Equatorial upwelling caused by Ekman divergence supplies cold, nitrogen-, and phosphorus-rich waters (Fig. 1B and fig. S1) (11), resulting in iron limitation of phytoplankton communities (1). In contrast, the subtropical gyres to the north and south are highly depleted in multiple nutrients (N, P, and micronutrients; Fig. 1B and fig. S1). The scarcity of multiple nutrients has made it difficult to experimentally determine a single primary limiting nutrient of primary productivity in oligotrophic regions (3), and as an alternative, marine ecosystem models have been used to predict iron and nitrogen limitation of the equatorial and subtropical gyre regions of the Pacific, respectively, with a sharp transition between them (2, 12). Prochlorococcus, an oxygenic photoautotrophic cyanobacterium, is a major contributor to marine photosynthesis in tropical and subtropical oceanic regions (1, 13) and was observed throughout the transect by its unique divinyl chlorophyll a pigment, particularly in the equatorial upwelling region when total chlorophyll a was also highest (Fig. 1C and fig. S1).

Fig. 1 Station locations, nitrate distributions, and vertical profiles of iron and nitrogen stress biomarkers.

(A) Cruise track for the METZYME expedition in the central Pacific Ocean aboard the R/V Kilo Moana in October 2011 from Hawaii to Samoa (station coordinates in table S1). (B) Nitrate distributions along the transect (units of micromolar). Vertical profiles of (C) nitrogen and (D) iron stress biomarkers P-II and IdiA, respectively, from the abundant marine cyanobacterium Prochlorococcus marinus as measured by targeted proteomics along the transect. With southward progression, a shift in nitrogen to iron stress was observed. Divinyl chlorophyll a distributions are also shown for comparison (in gray). Station 5 was located on the equator. Zero values below the photic zone, where Prochlorococcus is absent, demonstrate a lack of false positives.

Peptide biomarkers were selected from proteins of interest identified within global metaproteome analyses of the microbial samples taken from this transect. These metaproteome analyses used numerous cyanobacterial genomes and Pacific metagenomes to enable protein identification specific to this region (table S2). Targeted methods were designed using peptide-specific parent to fragment ion transitions (table S3), and synthesized isotope-labeled peptides served as internal standards using MRM on a triple quadrupole mass spectrometer (14). Several proteins had been shown in laboratory culture studies of marine cyanobacteria to have potential use as biomarkers of iron and nitrogen stress, including iron deficiency protein IdiA (5); flavodoxin, which replaces iron-requiring ferredoxin under iron stress (15, 16); and the global nitrogen response regulator NtcA that serves as a transcriptional activator for alternate forms of nitrogen (such as urea) and is under negative control by ammonium (17, 18). We also targeted the nitrogen regulatory protein P-II and the urea transporter of Prochlorococcus as additional nitrogen stress biomarkers based on their presence in the metaproteome. One challenge in deploying targeted proteomics to the oceans is the reliance of the MRM technique on exact sequence matches, which, in a complex environmental community, can result in the measurement of a subset of a protein’s functional sequence diversity (fig. S2). Yet, genomic analysis of cyanobacterial genomes demonstrated that peptides targeted for Prochlorococcus nitrogen and iron stress were specific not only to the Prochlorococcus species, but also showed subspecies specificity to the two currently known high-light ecotypes of Prochlorococcus in the upper water column [table S4, the low-B/A I and low-B/A II clades (19)], with the exception of one of the two nitrogen regulatory proteins (NtcA and P-II), where NtcA also targeted the related cyanobacterium Synechococcus, whereas P-II was specific to high-light Prochlorococcus. Hence, by selecting protein targets identified in metaproteomes that were generated with regional metagenomes and microbial genomes, we could target biomarkers that were representative and abundant for this region. Owing to this subspecies-level resolution, the data presented here are not normalized to Prochlorococcus population abundance, because defining and quantifying ecotypes can be challenging. Yet, the unnormalized protein concentration data are of biogeochemical value, as they can be directly related to ecological stoichiometry, nutrient uptake, and enzymatic activity.

The vertical profiles of these biomarker proteins identified a transition between regions of iron and nitrogen stress along the transect (Figs. 1 and 2). Biomarkers for IdiA and flavodoxin went from nondetectable to highly abundant upon entering the equatorial upwelling region as the integrated dissolved iron inventory (from 20 to 100 m) became depleted, from 24 μmol m−2 at station 1 to 3 μmol m−2 by station 7 (Fig. 2F). This decrease in integrated dissolved iron was observed despite a vertical structure that was more complex than that of nitrate owing to the regional hydrography and long-distance zonal transport of particulate iron (table S5) (11, 20). In contrast, P-II was abundant in the oligotrophic gyre region of the North Pacific, but became scarce in the equatorial Pacific as the integrated nitrate inventory (20 to 100 m) increased from 4 to 565 mmol m−2 (Fig. 2C). Nitrite and ammonia concentrations both increased at the equator and southward (fig. S6), consistent with P-II distributions. Several criteria corroborate this method for diagnosing the nutritional status of environmental populations of Prochlorococcus. First, multiple independent protein biomarkers provided a consistent indication of nutritional stress (P-II and urea transporter for nitrogen, and IdiA and flavodoxin for iron) and were oceanographically consistent with cyanobacterial proteins being nondetectable below the photic zone. Second, within two of these protein biomarkers, the urea transporter and flavodoxin, multiple peptides were targeted and yielded very similar geographical trends, with variations in abundance likely reflecting the diversity of protein sequence in the natural environment (table S2 and fig. S2). Third, inspection of individual fragmentation patterns (ms2 spectra) from metaproteomic analyses showed consistent fragment patterns and a lack of interfering ions for the MRM analyses (figs. S3 and S4).

Fig. 2 Ocean sections of Prochlorococcus nitrogen and iron metabolic protein biomarkers.

(A and B) P-II and urea transporter for nitrogen stress in units of fmol/liter (inset: transect station map; peptide IDs 34 and 41 from tables S3 and S4), and (D and E) IdiA and flavodoxin for iron stress in units of fmol/liter (peptide IDs 27 and 31). (C and F) Depth-integrated nitrate and dissolved iron for the 20- to 100-m region. At stations 3 and 4, samples were not available at 20 m, so the 40-m values were used as a surface depleted estimate. When no 100-m sample was available, a value was extrapolated from adjacent depths (see table S3).

These protein biomarkers also provided insight into the forms of nitrogen used by Prochlorococcus. Although phytoplankton, including Prochlorococcus, are often considered to prefer ammonia as a source of nitrogen (21), nitrite (21), nitrate (22), urea (8, 22, 23), and cyanate (24) have also been observed to be available to Prochlorococcus. In this study, the Prochlorococcus urea transporter protein (UrtA) was among the most abundant proteins that we detected with concentrations of 1000 to 2500 fmol liter−1 in the North Pacific Subtropical Gyre (NPSG), comparable to that of the iron transporter IdiA in the iron-limited regions to the south, implying that dissolved organic nitrogen was an important nutritional source for Prochlorococcus. Biomarkers for a subunit of the urease enzyme (UreC) and an ammonia transporter (Amt1) from high-light Prochlorococcus showed meridional distributions similar to that of P-II and the urea transporter, but both were found at much lower abundances compared to the urea transporter (fig. S5 and tables S3 and S4). Low urease abundance is likely related to the enzyme’s efficiency relative to the challenge of acquiring scarce urea. This urea utilization could be a useful chemical and ecological strategy for high-light Prochlorococcus in ammonia-limiting waters, allowing direct acquisition of this zooplankton waste product without the need for a heterotrophic microbial intermediary with urease capability. The low concentrations of ammonia and urea transporters on the equator (station 5, fig. S5) could reflect the use of the more abundant nitrate and nitrite (25) (nitrite was 1.8 to 7.5 times more abundant than ammonia, and nitrate was 17 times more abundant than nitrite in the upper 100 m; figs. S6 and S15). Alternatively, the abundant ammonia concentrations could have induced low ammonia transporter concentrations, or another ammonia transporter sequence was in use that was not being targeted in our study since the targeted peptide was only identified in one Prochlorococcus isolate genome (table S4).

Together these biomarker results reveal a gradation of adaptive responses between the NPSG and equatorial regions, with nitrogen and iron stress proteins co-occurring in the transition regions between biomes (Fig. 2). Comparing biomarker abundances with coherent and distinct ratios of P-II:IdiA and UrtA:IdiA peptide abundances at each station (Fig. 3, A and B, and table S6) showed the primary nutrient stress, consistent with previous incubation observations (26), yet also revealed co-stresses at biome transitions. Because the protein measurements were calibrated on an absolute scale, we calculated the cellular nitrogen being dedicated to both of these enhanced nitrogen and iron acquisition strategies for five Prochlorococcus proteins across this region (Fig. 3C and fig. S7; IdiA, flavodoxin, P-II, urea transporter, glutamine synthetase). Stations 3 and 9 had a considerable fraction of the nitrogen within this suite of five targeted proteins dedicated to the acquisition of both nutrients simultaneously (Fig. 3D and fig. S7, A and B). Although this estimate using selected proteins serves as a minimum nutritional estimate for the Prochlorococcus nitrogen and iron scarcity response, this analysis demonstrates that the community was deploying multiple nutrient acquisition strategies simultaneously, and hence was carrying the burden of an increased nutritional and energetic cost for maintaining both systems simultaneously as a potentially colimiting scenario. Moreover, although a previous study in a similar geographic region found that a novel low-iron ecotype dominated Prochlorococcus surface populations throughout (27), it is also possible that these biomarkers reflect an aggregated community signal of two or more coexisting high-light Prochlorococcus clades present in similar abundances. These results arguably present a challenge to current ecosystem parameterization influenced by Liebig’s law of the minimum, where the required nutrient that is lowest in abundance relative to cellular requirements (e.g., cellular stoichiometry) solely controls growth rate by applying a minimum function to multiple Monod growth expression terms (3, 4). Implicit in this minimum approach is an assumption that there is no metabolic or energetic penalty for the occurrence of more than one scarce nutrient simultaneously. Yet, these observations point to nutritional costs associated with maintaining the cellular machinery required to obtain a second scarce nutrient at the biomes’ intersections.

Fig. 3

Comparison of distributions of nitrogen and iron stress biomarkers across the transect. (A) Relationships between IdiA and P-II biomarkers, with shallower depths tending to have increased protein abundances. (B) Comparison of distributions of urea transporter abundances and IdiA iron stress biomarker. The coherent ratios for these pairs of biomarkers was indicative of the unique combination of environmental stimuli tuning the protein expression of nutrient-scarcity responses (ratios of biomarkers presented in table S6). (C) The nitrogen content within each targeted protein was calculated and summed for selected nitrogen and iron metabolism proteins along the transect (P-II, urea transporter, and glutamine synthetase for nitrogen metabolism and IdiA and flavodoxin for iron metabolism; peptide IDs 34, 41, 145, 27, 31) and (D) represented as a fraction of the total nitrogen for these five proteins. Color bars between (C) and (D) correspond to station symbols in (A) and (B).

Other nutrients, in addition to iron and nitrogen, were also scarce in this region and could influence primary productivity and community composition (3, 4). Phosphorus and cobalt, in particular, are required for growth of Prochlorococcus (28, 29) and were both highly depleted in the upper photic zone of the NPSG (Fig. 4, B and C). Two peptide biomarkers for the sulfolipid uridine 5-diphosphate (UDP) sulfoquinovose enzyme from high-light Prochlorococcus were found to be present at stations 1 and 3, but were diminished by station 5 as phosphate approached 1 μM in surface waters (Fig. 4A). Prochlorococcus uses this enzyme to replace phospholipids for sulfolipids as an adaptive response to phosphate scarcity (30). Deeper (low light) cyanobacterial communities responded to nutrient distributions along the transect. A protein homologous to the two-component phosphate response regulator PhoP in Bacillus subtilis (31) was identified during the metaproteome analyses (table S3) that corresponded to low-light Prochlorococcus ecotype and Synechococcus genomes (table S4). This two-component response regulator is one of seven identified in Prochlorococcus (32), yet was not previously observed to be responsive to phosphate scarcity (29). In this transect, this putative PhoP showed a similarity to the photic zone vertical distributions of phosphate and cobalt (Fig. 4), consistent with a role in regulation of phosphate metabolism near the nutricline. Notably, the expression of PhoP was opposite that of the two-component regulatory system PhoR identified in high-light Prochlorococcus ecotypes that is abundant under phosphate-scarcity rather than phosphate-replete conditions. PhoR is only present as a pseudogene in the low-light Prochlorococcus ecotype (29), implying that alternate response regulators such as PhoP are needed for regulation of low-light Prochlorococcus communities.

Fig. 4 Phosphorus-related proteins and phosphate and cobalt distributions.

(A and B) Vertical distributions of the sulfolipid UDP-sulfoquinovose protein and the putative two-component regulator PhoP at stations 1, 3, and 5 (peptide IDs 134 and 171; circle and triangle, respectively), in relation to (B) soluble reactive phosphate and (C) total dissolved cobalt and oxygen distributions (solid line).

The distributions of the putative PhoP also resembled those of dissolved cobalt, whose distribution often correlates with that of phosphate (Fig. 4). Although cobalt is an essential micronutrient for Prochlorococcus and is the rarest of the metallic micronutrients in the sea (28), there are as yet no cyanobacterial biomarkers known to represent cobalt stress. The similarity to cobalt is also consistent with a metal requirement for the alkaline phosphatase enzyme, globally observed photic zone correlations of cobalt and phosphate (33), and several potential metal-related proteins within the gene neighborhood of phoP (including a transporter, protease, and carbonic anhydrase). Utilization of the micronutrient nickel also varied across the transect, yet in contrast to the low-abundance nickel-containing urease mentioned above, the nickel superoxide dismutase enzyme was highly abundant and increased in concentration toward the southern half of the transect, presumably in response to oxidative stress associated with iron-limiting physiology (fig. S6). Consistent with nickel superoxide dismutase being much more abundant than urease, the nickel measured in filtered particulate material was more abundant on the equatorial and southern regions of the transect (fig. S6).

These results show the utility of targeted metaproteomics in diagnosing nutrient stress of a major phytoplankton group, studying changes in ecological stoichiometry, and identifying novel regulatory components and their potential role in the systems biology of the oceans. Future climate is predicted to alter the oceanic biogeochemical provinces (34), with marine cyanobacterial populations expected to be affected, although it has been difficult to predict the nature of these impacts (35, 36). Diagnostic measurements deployed over large geographic regions, such as the targeted metaproteomics methods shown here, could be valuable in characterizing the response to multiple influences.

Supplementary Materials

Materials and Methods

Supplementary Text

Figs. S1 to S18

Tables S1 to S6

References (3763)

References and Notes

  1. Materials and Methods are available on Science online.
  2. Acknowledgments: This work was made possible by NSF grants OCE-1031271, 1155566, and 1233261 Gordon Betty Moore Foundation grants 2724 and 3782, the Center for Microbial Research and Education, and the Woods Hole Oceanographic Institution Ocean Life Institute. We are grateful to the captain and crew of the R/V Kilo Moana, particularly V. Polidoro for his resourcefulness. We thank N. Hawco, D. Wang, and P. Balcom for sampling assistance and J. Dusenberry for programming assistance. We thank A. Santoro, J. Waterbury, and P. Chivers for helpful discussions and J. Blethrow and V. Zabrouskov (Thermo Scientific) for analyses on the prototype Fusion. Data for the Metyme Expedition are available at
View Abstract

Stay Connected to Science

Editor's Blog

Navigate This Article