Search tips
Search criteria 


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

Metagenomic Profiling of a Microbial Assemblage Associated with the California Mussel: A Node in Networks of Carbon and Nitrogen Cycling

Robert DeSalle, Editor


Mussels are conspicuous and often abundant members of rocky shores and may constitute an important site for the nitrogen cycle due to their feeding and excretion activities. We used shotgun metagenomics of the microbial community associated with the surface of mussels (Mytilus californianus) on Tatoosh Island in Washington state to test whether there is a nitrogen-based microbial assemblage associated with mussels. Analyses of both tidepool mussels and those on emergent benches revealed a diverse community of Bacteria and Archaea with approximately 31 million bp from 6 mussels in each habitat. Using MG-RAST, between 22.5–25.6% were identifiable using the SEED non-redundant database for proteins. Of those fragments that were identifiable through MG-RAST, the composition was dominated by Cyanobacteria and Alpha- and Gamma-proteobacteria. Microbial composition was highly similar between the tidepool and emergent bench mussels, suggesting similar functions across these different microhabitats. One percent of the proteins identified in each sample were related to nitrogen cycling. When normalized to protein discovery rate, the high diversity and abundance of enzymes related to the nitrogen cycle in mussel-associated microbes is as great or greater than that described for other marine metagenomes. In some instances, the nitrogen-utilizing profile of this assemblage was more concordant with soil metagenomes in the Midwestern U.S. than for open ocean system. Carbon fixation and Calvin cycle enzymes further represented 0.65 and 1.26% of all proteins and their abundance was comparable to a number of open ocean marine metagenomes. In sum, the diversity and abundance of nitrogen and carbon cycle related enzymes in the microbes occupying the shells of Mytilus californianus suggest these mussels provide a node for microbial populations and thus biogeochemical processes.


In many locales in coastal oceans, nitrogen has been demonstrated to be the limiting nutrient, with large-scale circulation patterns (such as upwelling) being the primary determinant of coastal productivity. Although circulation patterns that drive upwelling can import substantial amounts of nitrate into coastal areas, regeneration of nitrogen in situ can also contribute to local productivity [1][3]. Regenerated nitrogen is mostly due to the metabolism and excretion of animals, while marine plants, seaweeds and microbes utilize the nitrogenous waste. Although the response of some coastal eukaryotic primary producers to nitrogen production by animals has been described [4][6], microbial population abundance and diversity in response to nitrogen is less studied. Nonetheless, there is ample evidence that microbes are ubiquitous consumers of nitrogeneous byproducts from animals, chemolithotrophy is well-established, and there is a great potential for regenerated nitrogen availability to drive enhanced carbon dioxide fixation.

Despite the importance of nitrate delivery with upwelling along the margins of northeast Pacific Ocean, ammonium excretion by animals is detectable [7][10] and has been shown to contribute to local productivity [5], [6], [10] and diversity [11]. Although marine mammals, seabirds, fishes and dense aggregations of invertebrates all may contribute to regenerated nitrogen in coastal areas, mussels (Mytilus californianus, henceforth mussels) have only recently been recognized as significant contributors [6], [10]. Experimental manipulation of the presence of mussels demonstrated that ammonium excretion by invertebrates not only boosts the productivity of macroalgae, but also drives microbial productivity via nitrification [6]. The use of animal-regenerated nitrogen for chemolithotrophy by marine microbes has been relatively ignored in these well-studied rocky shores; arguably, their abundance and function is probably better understood in the open ocean [12],[13] and deep sea environs [14]. To date, we know relatively little about the identity or function of rocky shore microbes and their importance to nitrogen and carbon cycling. Marine benthic nearshore microbes may play an important role mediating the abundance of different forms of nitrogen via nitrification, ammonification, detnitrification and potentially all aspects of the nitrogen cycle. Additionally, they are likely providing increased opportunities for microbial CO2 fixation, while also competing with other primary producers, including the ecologically important macroalgae, for nitrogen. Here we describe shotgun metagenomic-based analysis of the microbes associated with mussels including analyses of their function in rocky shore ecosystems.

It is thought that many microbial taxa cannot be cultured outside of their natural environment; thus, microbial diversity remains poorly described [15], [16]. The metagenome techniques developed recently have therefore greatly extended our knowledge of microbial genetic diversity [17][19]. Because they are acclimated to high energy waves and cold temperatures, many rocky shore species, including microbes, are difficult to accommodate in laboratory environs. The recent findings of the previously undescribed nitrifying Archaea in a diversity of habitats [20][21], suggest that there is much microbial diversity yet to be described. Additionally, the ability to analyze vast numbers of genomes allows probable metabolic functions to be determined [22]. Because we had strong experimental evidence that microbial nitrification was present in tidepools with abundant mussels [6], we hypothesized that these microbes would live in close proximity to a reliable source of both habitat and ammonium – the shells of the mussels themselves. We further hypothesized microbial assemblages would be common to mussels in a variety of habitats on rocky shores, due to their dominance and abundance [23], [24]. Indeed, mussels average densities of mussels are 4661 per m2 on Tatoosh Island [25], the site of the work reported here. We thus report metagenome analyses of the microbial community obtained from shells of mussels, including separate analyses of the community from tidepool mussels versus those from mussels that reside on rock that is emergent at low tide. Specifically, we ask about the taxonomic affiliations of these microbial communities as well as the likely function of these microbes given their affiliations and their sequence homology with enzymes of known function in nitrogen metabolism.

The increasing public availability of environmental metagenomes has further allowed us to compare our mussel microbial assemblage both in terms of taxonomy and metabolism to other ecosystems. We further use results from other marine ecosystems to test whether mussel-associated microbes have similar nitrogen-based metabolism.

Materials and Methods

Mussels were collected from the Main Beach site of Tatoosh Island (48.32°N, 124.74°W), located in the eastern Pacific 0.7 km off the northwestern tip of Washington State, USA. Six mussel shells were collected from among 6 tidepools, while 6 more were collected at a distance of approximately 5 m apart on an adjacent exposed bench on 10 April 2008 and immediately cleaned of all soft tissue. The shells (mean length 4.47 cm and 4.42 cm for tidepool and bench mussels respectively) were put on ice and brought to Argonne National Labs.

DNA was extracted and purified using Ultraclean Mega Prep Soil DNA Isolation Kit and following directions therein (MO BIO Laboratories,Inc.) and the two extractions are referred to as tidepool versus bench mussels. The tidepool sample yielded 4320 ng in 108 µL (Invitrogen Qubit fluorometer dsDNA HS Kit), while bench mussels had 168 ng in 350 µL and required use of the GenomiPhi V2 DNA Amplification Kit (GE Healthcare). We followed the Roche GS-FLX (454) shotgun library preparation protocol; the tidepool sample used 2.4 µg and the bench sample used 5.0 µg for library preparation. Both samples had a mean fragment size of 750 bp after library preparation. All sequencing was performed with the 454 GS-FLX instrument and LR70 sequencing chemistry (Roche Applied Science).

We analyzed the taxonomic composition of our two metagenome sample sets with the MG-RAST server [26] using similarity to a large non-redundant protein database. Using the same non-redundant database, we also tested the affinities of our sequences for known metabolic function against both SEED subsystems [27] and KEGG metabolic pathways [28] using a maximum e-value of e<10−5. Although there are a number of metabolic functions that can be tested, our specific interest in microbial contributions to the nitrogen cycle focused our efforts on both nitrogen metabolism and carbon dioxide fixation. Thus, we probed particularly for enzymes related to the components of nitrogen and CO2 use.

In addition to describing the taxonomic and metabolic features of this microbial community on mussels, we also tested the similarity and differences with other recently described marine microbial assemblages that are public, including those of coastal Georgia [29], 4 tropical Pacific Ocean seawater samples in the Line Islands [19], and the extensive Global Ocean Sampling Expedition [13]. For the latter, we chose for comparison 4 coastal locales that that spanned a wide geography and sampled surface waters, including the Gulf of Maine (GS002, MG-RAST id #4441579.3), Nag's Head, NC (GS013, 4441585.3), Cocos Island, Costa Rica (GS025, 4441593.3), and an upwelling zone off of Fernandina, Galapagos (GS031, 4441597.3). We excluded marine metagenome analyses that had selectively filtered and extracted samples to isolate viruses. We focused our analyses on nitrogen metabolism and CO2 fixation to test the similarities and differences of our mussel-associated microbes. Given the abundance of nitrogen in our mussel-associated waters, we further asked if another nitrogen-rich ecosystem, soils of the agriculture-influenced midwest, showed metabolic similarities. Here, we compared our mussel microbial assemblage to soil samples from Midwestern locales (Waseca farm soil (4441091.3), soybean field (4442657.3), prairie remnant (4442656.3), 2nd yr prairie (4442658.3), 20th year prairie (4442659.3), 33rd year prairie (4441281.3)). For all comparisons, we used a non-redundant protein database with an e-value cut-off of 10−5. We recognize that the ‘discovery rate’ for proteins may depend upon the efficacy of DNA extraction and the length of sequences that result, features that may vary among studies. Although we normalized the number of proteins identified with different metabolic functions by the number of proteins that were found per 100 fragments, we had no means of controlling for the different contiguous sequence lengths that occurred among different studies.

The mussel associated sequences are publicly available in the MG-RAST system under the following project identifiers (IDs 4441185.3 (tidepool), 4441191.3 (emergent,bench). The data in this manuscript and the analyses and comparisons to other public data sets are available via MG-RAST. MIGS/MIMS [30] compliant metadata describing the locations, sampling, data extraction and data is available in GCDML [31] format from within the MG-RAST system as well.


Phylogenetic Analyses

For the tidepool mussel sample, there were 157,599 total DNA fragments with a total sequence size of 30,593,565 and an average sequence length of 194 bp. The bench mussel sample had slightly fewer contiguous sequences (141,293) from a similar sequence size of 31,304,272 and an average sequence length of 222 bp.

The BLASTX analysis against a non-redundant protein database matched 22.5% of the sequences in the tidepool sample of which 74% were bacterial and almost 3% were eukaryotic; the remaining 23% were unidentified. 1% of protein sequences matched to nitrogen metabolism. For the bench mussels, approximately 79% were bacterial with 2% eukaryotic and .4% Archaeal for the 25.6% that could be matched. Our protein ‘discovery rate’ of 22.5 and 25.6% was comparable or greater than other metagenome studies using 454-based sequencing technology [19], [29], but less than studies where direct library construction and sequencing was done [13].

When we compared the taxonomic composition of the 2 mussel samples to each other, they were similar at higher taxonomic organization, but differed slightly in the composition of lower taxonomic groupings (Table 1. Fig. 1a). Cyanobacteria, α-Proteobacteria and γ-Proteobacteria dominated both samples. The Cyanobacteria were more abundant on emergent mussels and were identified primarily as members of the orders Chroococcales and Nostocales (Fig. 2a). Crocosphaera and Synechococcus were identified in both samples, though more on emergent mussels. Both genera are photoautotrophs that are not thought to fix atmospheric nitrogen. The α-Proteobacteria were dominated by Rhizobiales and Rhodobacterales and included the nitrifying Nitrobacter (Rhizobiales) (Fig. 2b). There was an increased incidence of Rhodobacterales on tidepool mussels, including matches with Rhodobacter and Roseobacter, an aerobic anoxygenic phototroph. The β-Proteobacteria were highly similar between tidepool and emergent mussels and the nitrifying Nitrosomonas and Nitrospira were represented in both (Fig. 2c). The γ-Proteobacteria was the taxonomic unit with the greatest membership and was primarily composed of Vibrionales and Alteromonadeles (Fig. 2d). The ammonium oxidizing bacterium Nitrosococcus was represented in both samples. Although relatively few Archaeal proteins were identified, they included representatives of both Crenarchaeota and Euryarchaeota, and Nitrosopumilus, an ammonia-oxidizing chrenarchaeon, was detected in both samples.

Figure 1
Taxonomic composition of surface-associated microbes of tidepool and emergent (bench) mussels.
Figure 2
The relative proportional representation within the most commonly discovered bacterial orders on the surface of the mussel shells.
Table 1
The taxonomic diversity of microbes from the surface of mussels in tidepools and on emergent benches.

When we compared the emergent and tidepool mussels at the finest level for taxonomic affinities, only 7 identities differed between the 2 samples and all were single occurrences within 7 distinct phylogenetic groups (Crenarcheaota, Euryarchaeota, Actinobacteria, Chlorobi, Firmicutes, γ-Proteobacteria). Thus, the two mussel microbe assemblages were highly concordant in their overall composition, despite the fact that they came from different microhabitats.

When we compared the taxonomic composition of mussel shell microbes with other marine metagenomes, the dominance of γ-Proteobacteria in mussels and stromatolites, e.g. the nearshore sites, is apparent (Table 2). α-Proteobacteria were better represented in open ocean waters, though bench mussels had a large representation too. The representation by Cyanobacteria varied among sites with the 2 Line Islands of Fanning and Palmyra having high representation, though primarily by Chroococcales at Fanning (and on the bench mussels) and by Prochlorales at Palmyra. In contrast, the more oceanic Prochlorales were low in incidence on the mussel shells, represented by only 77 Prochlorococcus hits in each of the bench and tidepool samples.

Table 2
The comparative representation of bacterial and archaeal phylogenetic groups on mussel shells versus other marine systems with metagenomes analyzed by shotgun pyrosequencing and using MG-RAST and the SEED subsystem database (e<10−5)).

Metabolic Analyses

Our mussel associated metagenome analysis found many matches to proteins in the non-redundant database relevant to metabolic functions (Table 3). The relevant ranking of metabolic functions was strikingly similar to Dinsdale et al.'s [18] ranking based on the mean of 45 microbial metagenomes from habitats as diverse as the digestive systems of animals to a coral holobiont.

Table 3
Percentage of sequences that matched major metabolic categories (Subsystem Categories [17]) using the SEED non-redundant database for both tidepool mussels and emergent, bench mussels and compared with a mean value for other microbial ...

We hypothesized that enzymes related to ammonium assimilation would be present in mussel shell microbes as a means of utilizing the ammonium excreted by mussels. When we used the protein database to match to metabolic function, we found 1.0% of the sequences in each sample matched to nitrogen metabolism, with a total of 446 sequences found in the tidepool mussels and 445 in the bench mussels. The distribution of sequences associated with different aspects of nitrogen cycling were relatively similar among the 2 samples (Table 4), and included not only ammonium assimilation, but also nitrate and nitrite ammonification, allantoin degradation and nitric oxide synthase as the dominant metabolic components. Enzymes such as ammonium monooxygenase subunit A (amoA) and glutamine synthetase were also detected. Denitrifying enzymes were also present, while nitrogen fixation (as indicated by nitrogenase) was nearly absent.

Table 4
Number of sequences associated with nitrogen metabolism using the SEED database.

The taxonomic affiliations of the enzymes involved in nitrogen metabolism bore strong similarity to the overall representation of Bacteria and Archaea in the samples (Fig. 1b). Thus, all major groups of microbes contributed to nitrogen metabolism in approximate proportion to their abundance, although Delta-Proteobacteria were better represented in the bench mussels and Cyanobacteria and Gamma-proteobacteria were also strongly associated with nitrogen metabolism.

Nitrogen metabolism enzymes in the mussel shell microbes show a strong pattern for much uptake and transformation of inorganic nitrogen especially ammonium uptake and ammonification, a pattern shared with some of the Line Islands metagenomes and also with the waters surrounding the Galapagos upwelling region (Fig. 3a). Other regions were comparatively depauperate in proteins for nitrogen function, including seawater from areas adjacent to Georgia, Maine, North Carolina and Costa Rica. Nitrogen fixation was suggested to be relatively minor in these areas, excepting the Georgia VAN sample. When using MG-RAST to test the hypothesis that our mussel associated microbes would show strong similarity with soil metagenomes from current or former agricultural fields of Illinois, the similarity of enzyme types was marked (Fig. 3b). In these soils, as in association with mussels, enzymes related to ammonium uptake or nitrite and nitrate use were particularly well represented, though soils had an increased incidence of enzymes related to nitrogen fixation.

Figure 3
The number of proteins matched to nitrogen metabolism functions among multiple metagenome studies.

Proteins related to CO2 fixation were also well represented in the mussel shell samples, including enzymes of the Calvin cycle (primarily RUBISCO) and exceeded those in Georgia waters and the 4 Global Ocean samples (Figure 4). We matched 232 CO2 fixation-related proteins in the tidepool sample (0.65% of all proteins) and nearly twice that in the bench sample (456 for 1.26% of all proteins). The Line Islands had a variable amount of CO2 fixation proteins with Palmyra and Kiritimati (Christmas Island) having the greatest number. The percent of proteins identified for CO2 fixation out of the entire protein pool of all metagenomes ranged between 0.0 and 1.38 across all metagenomes.

Figure 4
The number of proteins involved in CO2 fixation, including those of the Calvin-Benson cycle among a selection of marine metagenomes.


The microbes on mussel shells showed both a taxonomic and functional composition that reflects a nitrogen-rich environment. The major nitrifying bacterial genera that are described (Nitrosococcus, Nitrosomonas, Nitrospira, Nitrobacter, [32]) were found in both samples, as well as some Crenarchaeota (Nitrosopumilis, [33]). Ammonium assimilation enzymes were also well represented in both samples. Denitrifying genera were also identified, including Shewanella and Roseobacter, as well as enzymes related to denitrification including nitrite- and nitrate-reductases and all were more prevalent in the tidepool rather than bench mussels. Although the presence of denitrifying enzymes and genera suggest that this process may be occurring in low oxygen microsites in tidepools, its signature feature, the uptake of nitrite is not suggested by water nutrient sampling within tidepools [6]. Whether anaerobic ammonium oxidation (anammox) is important here is unclear. We found no matches to the putative species or enzymes thought to be important to anammox [34], but recognize that our detection may be limited by the relatively little that is known about anammox metabolism. However, as with denitrification, these environments are typically well-oxygenated and we might not expect anammox to be highly important. Although Cyanobacteria and many genera from Proteobacteria that are known nitrogen fixers were represented in our mussel shell microbe samples [35], nitrogenase enzymes were absent, suggesting that these are photoautotrophs that do not fix atmospheric nitrogen. The rich variety of other nitrogen sources in this nearshore environment may select against nitrogen fixation or result in competitive inferiority of nitrogen fixers compared with ammonium or nitrate utilizing microbes, a pattern described in plankton assemblages [36][38]. The potential for relatively high ambient availability of ammonium can be illustrated using Suchanek's [25] estimate of a mean number of 4661 mussels per m2 on Tatoosh Island coupled with per mussel excretion rates [39]. Using both, we estimate >3 g of ammonium (~55 mmol) excreted per day per m2 of mussels, a substantial input of inorganic nitrogen. If rates were known for other invertebrates and vertebrates in this system, such as seabirds and marine mammals, this input would likely be much higher. Nitrate from upwelling is also typically high, and is at concentrations of approximately 20 µM at Tatoosh and nearby sites during spring and summer months [6], [9]. In sum, the nearshore in proximity to M. californianus mussel beds are rich in inorganic nitrogen and may provide an environment that selects against the persistence of nitrogen fixing organisms.

Our comparison of nitrogen metabolism among different systems further substantiated the diverse nitrogen functions on mussel shells. The Line Islands in the South Pacific, the coastal ocean samples of Georgia and the Global Ocean Sampling Expedition, excepting the upwelling region off the Galapagos (Fernandina, GS0031), all had fewer enzyme matches for ammonium assimilation. In terms of the distribution of enzymes related to nitrogen metabolism, mussel shell microbes had much in common with the upwelling region of the Galapagos and Midwestern soils, including a range of different metabolisms with ammonium assimilation and ammonification predominating, though nitrogen fixation was much more common in these soils than indicated for the surface of mussels. The mussels, Fernandina, and the Midwestern soils likely have a rich nitrogen environment that promotes similar microbial opportunities. We note, however, that although all these metagenome studies were analyzed with a common platform (MG-RAST and SEED), comparisons among metagenome studies to date should consider differences among studies in extracting and sequencing. All but the Global Ocean Sampling Expedition used similar shotgun sequencing methodologies. However, the length of reads in the Line Island and Georgia samples were ~100 bp (Roche GS-20) compared with the ~200 bp and greater length in this mussel study and in the soil metagenomes (Roche GS-FLX). The cloning and Sanger-based sequencing approach used by the Global Ocean Sampling program (GS002,GS013,GS025, GS031), however, generated longer sequences (~1000 bp) and thus may have a higher protein discovery rate, although cloning bias would also be a factor.

The microbial assemblages of pool and bench mussels were very similar taxonomically and functionally, indicating that previous results suggesting nitrification in tidepools is probably a phenomenon general to association with mussels regardless of habitat. These few differences have some interesting implications. For example, Cyanobacteria were more abundant on emergent mussels. In the absence of evidence for nitrogenase, it is possible that these are endolithic phototrophs of mussel shells. Research with other mussels (Perna perna) have shown a detrimental effect of Cyanobacteria that are phototrophic endoliths [40], [41]. Although these endolithic forms appear poorly described taxonomically and ecologically, the possible increased incidence on emergent rock suggests that different environmental conditions may affect the distribution of these microbes. The reduced relative composition of Cyanobacteria in tidepools was possibly compensated for with α-Proteobacteria in tidepools, particularly Rhodobacterales, a group referred to as primary surface colonizers [42]. The high density of molluscan grazers in tidepools and their near continuous opportunities to graze might suggest some grazer tolerance or resistance on the part of these α-Proteobacteria, a result supported by experimental grazer removals in marine benthic systems [43]. In terms of nitrogen metabolism the tidepool and emergent bench mussels were very similar with tidepools having only a slightly greater incidence of ammonification.

Some compositional differences between mussel shell microbes and other marine metagenomes was marked. For example, our mussel shell microbes and the Georgia coastal waters were dominated by γ-proteobacteria; dominance by γ-proteobacteria has also been demonstrated in association with the Caribbean coral Porites astreoides [44]. In contrast, open ocean samples such as those from the Sargasso Sea [17] and the locales of the Global Ocean Sampling Expedition (Fig. 3, [13]) were dominated by α-Proteobacteria, a group associated with photoautotrophs of the open ocean. The Cyanobacteria were also variable in abundance, with the bench mussels and the Costa Rican sample (GS0025, Cocos Island) having a relatively large proportion of Cyanobacteria. However, the composition of the Cyanobacteria differed among these sites. The Costa Rican waters were dominated by the order Prochlorales (genus Prochlorococcus), a group that had only 77 contigs per sample in the mussels, a meagre ~.14% of all identifiable sequences. Although Prochlorococcus is known to be an abundant cyanobacterium in the open ocean where it can comprise as much as half of the photosynthetic biomass [45], it did not dominate in this nearshore environment.

Microbial activity related to the carbon cycle also appears to be a strong feature of the assemblage on mussel shells. There were as much or more proteins identified with CO2 fixation for the mussel-associated metagenomes as there were for any of the other marine metagenomes studied, suggesting that microbial nitrogen and carbon cycling are prevalent on these shells. Although there are at least 5 possible microbial pathways for CO2 fixation by microbes [46], the Calvin cycle is likely to be the most prevalent, based on the abundance of Proteobacteria and Cyanobacteria in these samples. The relative low incidence of microbes associated with other carbon fixation pathways (green sulfur bacteria, Chloroflexi) and the aerobic nature of the environment, make anaerobic and anammox pathways less likely.

All the described taxonomic and metabolic diversity came from a surface sample of only 6 mussels in an area where mussels can number in the thousands per square meter, and thus indicates the quantitatively significant role that mussels may play in microbial transformations for the nearshore nitrogen and carbon cycles. We acknowledge, however, that we have no water column censuses nor analyses of other substrates and cannot exclude the possibility that other substrates also serve as nitrogen transforming areas. Further genetic analyses in this system are thus warranted. Whether mussels are alone or not in providing suitable habitat for these microbial populations, the genetic data presented in this study suggests that if nitrogen is continually recycled and transformed by this microbial assemblage, then this provides a significant mechanism for the retention of nitrogen in nearshore areas, thus ameliorating the advection of nitrogen during upwelling events. Whether mussels are a unique node for microbial and biogeochemical activity, or one of several, the threats to their persistence are numerous and include declining ocean pH [47], low oxygen events [48], changing thermal environments [49], anthropogenic nitrogen pollution [50], and toxic algal blooms [51]. Although the work summarized here adds to our understanding of the interaction between macrofauna and their microbial associates, it also underscores how little of this diversity has been described in habitats that are otherwise well characterized in their ecological dynamics.


We thank A. Ammar and M. Domanus for technical expertise in sequencing and R. Edwards and an anonymous reviewer for helpful commentary. The Makah Tribal Nation graciously allowed access to the study areas.


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

Funding: Funding was provided by National Science Foundation OCE-0928232 to CAP and Argonne National Labs. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Dugdale RC, Goering JJ. Uptake of new and regenerated forms of nitrogen in primary productivity. Limnol Oceanogr. 1967;12:196–206.
2. Eppley RW, Reuger EH, Venrick EL, Mullin M. A study of plankton dynamics and nutrient cycling in the Central Gyre of the North Pacific Ocean. Limnol Oceanogr. 1973;18:534–551.
3. Bode A, Barquero S, Gonzalez N, Alvarez-Ossorio M, Varela M. Contribution of heterotrophic plankton to nitrogen regeneration in the upwelling ecosystem of A. Coruña (NW Spain). J Plankton Res. 2004;26:11–28.
4. Wootton JT. Direct and indirect effects of nutrients on intertidal community structure: variable consequences of seabird guano. J Expl Mar Biol Ecol. 1991;151:139–153.
5. Bracken MES. Invertebrate-mediated nutrient loading increases growth of an intertidal macroalga. J Phycol. 2004;40:1032–1041.
6. Pfister CA. Tidepool mussels locally increase nutrients and algal growth. Ecology. 2007;88:1647–1653. [PubMed]
7. Hansen JE. Marine Plants. In: LeBoeuf BJ, editor. Pages 183-204 in The Natural History of Año Nuevo. Pacific Grove, CA: Boxwood Press; 1981.
8. Jensen SL, Muller-Parker G. Inorganic nutrient fluxes in anemone-dominated tide pools. Pac Sci. 1994;48:32–43.
9. Pfister CA, Wootton JT, Neufeld C. The relative roles of coastal and oceanic processes in determining physical and chemical characteristics of an intensively sampled nearshore system. Limnol Oceanogr. 2007;52:1767–1775.
10. Aquilino KM, Bracken MES, Faubel MN, Stachowicz JJ. Local-scale nutrient regeneration facilitates seaweed growth on wave-exposed rocky shores in an upwelling system. Limnol Oceanogr. 2009;54:309–317.
11. Bracken MES, Nielsen KJ. Diversity of intertidal macroalgae increases with nitrogen loading by invertebrates. Ecology. 2004;85:2828–2836.
12. Barber RT. Germany: Ecology Institute; 2001. The Response of Oceanic Ecosystems to the Climate of the 21st Century. Excellence in Ecology 13, Oldendorfe/Luhe,
13. Rusch DB, Halpern AL, Sutton G, Heidelberg KB, Williamson S. The Sorcerer II Global Ocean Sampling expedition: Northwest Atlantic through eastern tropical Pacific. PLoS Biol. 2007;5:e77. doi: 10.1371/journal.pbio.0050077. [PMC free article] [PubMed]
14. Van Dover CL. The Ecology of Deep-Sea Hydrothermal Vents. 2000. Princeton University Press, Princeton, NJ.
15. Pace NR. A molecular view of microbial diversity and the biosphere. Science. 1997;276:734–740. DOI: 10.1126/science.276.5313.734. [PubMed]
16. Sogin ML, Morrison HG, Huber JA, Welch DM, Huse SM, et al. Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc Natl Acad Sci. 2006;103:12115–12120. [PubMed]
17. Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, et al. Environmental genome shotgun sequencing of the Sargasso Sea. Science. 2004;304:66–74. [PubMed]
18. Dinsdale EA, Edwards RA, Hall D, Angly F, Breitbart M, et al. Functional metagenomic profiling of nine biomes. Nature. 2008;452:629–632. doi10.1038/nature0681. [PubMed]
19. Dinsdale EA, Pantos O, Smriga S, Edwards RA, Angly F, et al. Microbial ecology of four coral atolls in the Northern Line Islands. PlosOne. 2008;3:e1584. doi/10.1371/journal.pone.0001584. [PMC free article] [PubMed]
20. DeLong EF. Archaea in coastal marine environments. Proc Natl Acad Sci. 1992;89:5685–5689. [PubMed]
21. Francis CA, Roberts KJ, Beman JM, Santoro AE, Oakley BB. Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc Natl Acad Sci. 2005;102:14683–14688. [PubMed]
22. Hallam SJ, Putnam N, Preston CM, Detter JC, Rokhsar D, et al. Reverse methanogenesis: testing the hypothesis with environmental genomics. Science. 2004;305:1457–1462. [PubMed]
23. Paine RT. Ecological determinism in the competition for space. Ecology. 1984;65:1339–1348.
24. Wootton JT. Field parameterization and experimental test of the neutral theory of biodiversity. Nature. 2005;433:309–312. [PubMed]
25. Suchanek TH. The Mytilus californianus community: studies on the composition, structure, organization, and dynamics of a mussel bed. Dissertation, Univ of Washington 1979
26. Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, et al. The metagenomics RAST server – a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386. doi: 10.1186/1471-2105-9-386. [PMC free article] [PubMed]
27. Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Research. 2005;33:5691–5702. [PMC free article] [PubMed]
28. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Research. 2007. doi: 10.1093/nar/gkm882. [PMC free article] [PubMed]
29. Mou X, Sun S, Edwards RA, Hodson RE, Moran MA. Bacterial carbon processing by generalist species in the coastal ocean. Nature. 2008;451:708–711. doi: 10.1038/nature06513. [PubMed]
30. Field D, Garrity G, Gray T, Morrison N, Selengut J, et al. The minimum information about a genome sequence”(MIGS) specification. Nature Biotech. 2009;26:541–547. [PMC free article] [PubMed]
31. Kottmann R, Gray T, Murphy S, Kagan L, Kravitz S, et al. A standard MIGS/MIMS compliant XML schema: toward the development of the genomic contextual data markup language (GCDML). OMICS: A Journal of Integrative Biology. 2008;12:115–121. [PubMed]
32. Ward BB. Kirchman DL, editor. Nitrification and the Marine Nitrogen Cycle. Pp 427-453 in Microbial Ecology of the Oceans. 2000. Wiley.
33. Könneke M, Bernhard AE, de la Torre JR, Walker CB, Waterbury JB, et al. Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature. 2005;437:543–546. [PubMed]
34. Francis CA, Beman JM, Kuypers MMM. New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation. The International Society for Microbial Ecology Journal. 2007;1:19–27. [PubMed]
35. Paerl HW, Zehr JP. Kirchman DL, editor. Marine nitrogen fixation. Pp 387-426 in Microbial Ecology of the Oceans. 2000. Wiley.
36. Smith VH. Low nitrogen to phosphorus ratios favor dominance by blue-green algae in lake phytoplankton. Science. 1983;221:669–671. [PubMed]
37. Montoya JP, Holl CM, Zehr JP, Hansen A, Villareal TA, et al. High rates of N2 fixation by unicellular diazotrophs in the oligotrophic Pacific Ocean. Nature. 2004;430:1027–1032. [PubMed]
38. Agawin NSR, Rabouille S, Veldhuis MJW, Servatius L, Hol S, et al. Competition and facilitation between unicellular nitrogen-fixing cyanobacteria and non-nitrogen-fixing phytoplankton species. Limnol Oceanogr. 2007;52:2233–2248.
39. Bayne BL, Bayne CJ, Carefoot TC, Thompson RJ. The physiological ecology of Mytilus californianus Conrad. Oecologia. 1976;22:211–228.
40. Kaehler S. Incidence and distribution of phototrophic shell-degrading endoliths of the brown mussel Perna perna. Mar Biol. 1999;135:505–514.
41. Kaehler S, McQuaid CD. Lethal and sub-lethal effects of phototrophic endoliths attacking the shell of the intertidal mussel Perna perna. Mar Biol. 1999;135:497–503.
42. Dang H, Li T, Chen M, Huang G. Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl Environ Microbiol. 2008;74:52–60. [PMC free article] [PubMed]
43. Hillebrand H, Kahlert M, Haglund A, Berninger U, Nagel S, et al. Control of microbenthic communities by grazing and nutrient supply. Ecology. 2002;83:2205–2219.
44. Wegley L, Edwards R, Rodriguez-Brito R, Liu H, Rohwer F. Metagenomic analysis of the microbial community associated with the coral Porites astreoides. Environmental Microbiology. 2007;9:2707–19. [PubMed]
45. Partensky F, Hess WR, Vaulot D. Prochlorococcus, a marine photosynthetic prokaryote of global significance. Microbiol Mol Biol Rev. 1999;63:106–127. [PMC free article] [PubMed]
46. Thauer RK. A fifth pathway of carbon fixation. Science. 2007;318:2009. [PubMed]
47. Wootton JT, Pfister CA, Forester JD. Dynamical patterns and ecological impacts of changing ocean pH in a high-resolution multi-year dataset. Proc Natl Acad Sci. 2008;105:18848–18853. [PubMed]
48. Grantham BA, Chan F, Nielsen KJ, Fox DS, Barth JA, et al. Upwelling-driven nearshore hypoxia signals ecosystem and oceanographic changes in the northeast Pacific. Nature. 2004;429:749–754. doi: 10.1038/nature02605. [PubMed]
49. Helmuth B, Harley CDG, Halpin P, O'Donnell M, Hofmann GE, et al. Climate change and latitudinal patterns of intertidal thermal stress. Science. 2002;298:1015–1017. [PubMed]
50. Howarth R, Anderson D, Cloern J, Elfring C, Hopkinson C, et al. Nutrient Pollution of Coastal Rivers, Bays, and Seas. Issues in Ecology. 2000;7:1–17.
51. Trainer VL, Hickey BM, Lessard EJ, Cochlan WP, Trick CG, et al. Variability of Pseudo-nitzschia and domoic acid in the Juan de Fuca eddy region and its adjacent shelves. Limnol Oceanogr. 2009;54:289–308.

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