Skip to main content


Front. Microbiol., 25 March 2014
Sec. Aquatic Microbiology
Volume 5 - 2014 |

The microbial nitrogen cycling potential is impacted by polyaromatic hydrocarbon pollution of marine sediments

  • 1Institute of Genomic and Systems Biology, Argonne National Laboratory, Lemont, IL, USA
  • 2Department of Ecology and Evolutionary Biology, University of Chicago, Chicago, IL, USA
  • 3Energy and Efficiency Division, Chemical and Biological Process Development Group, Pacific Northwest National Laboratory, Richland, WA, USA
  • 4Systems Microbiology and Biotechnology Group, Washington State University, Richland, WA, USA
  • 5Ecology Department, Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
  • 6Earth, Ocean and Atmospheric Science, Florida State University, Tallahassee, FL, USA

During hydrocarbon exposure, the composition and functional dynamics of marine microbial communities are altered, favoring bacteria that can utilize this rich carbon source. Initial exposure of high levels of hydrocarbons in aerobic surface sediments can enrich growth of heterotrophic microorganisms having hydrocarbon degradation capacity. As a result, there can be a localized reduction in oxygen potential within the surface layer of marine sediments causing anaerobic zones. We hypothesized that increasing exposure to elevated hydrocarbon concentrations would positively correlate with an increase in denitrification processes and the net accumulation of dinitrogen. This hypothesis was tested by comparing the relative abundance of genes associated with nitrogen metabolism and nitrogen cycling identified in 6 metagenomes from sediments contaminated by polyaromatic hydrocarbons from the Deepwater Horizon (DWH) oil spill in the Gulf of Mexico, and 3 metagenomes from sediments associated with natural oil seeps in the Santa Barbara Channel. An additional 8 metagenomes from uncontaminated sediments from the Gulf of Mexico were analyzed for comparison. We predicted relative changes in metabolite turnover as a function of the differential microbial gene abundances, which showed predicted accumulation of metabolites associated with denitrification processes, including anammox, in the contaminated samples compared to uncontaminated sediments, with the magnitude of this change being positively correlated to the hydrocarbon concentration and exposure duration. These data highlight the potential impact of hydrocarbon inputs on N cycling processes in marine sediments and provide information relevant for system scale models of nitrogen metabolism in affected ecosystems.


Petroleum hydrocarbon spills considerably alter the composition and functional dynamics of marine microbial communities (Hazen et al., 2010; Gutierrez et al., 2013). Ultimately, microorganisms that can respond to complex hydrocarbon mixtures are preferentially enriched by hydrocarbons provided during an oil spill; growing from small natural seed populations (Röling et al., 2002, 2004). A number of studies have demonstrated changes in the spatiotemporal ecology of microbial communities in the presence of oil contamination (Bordenave et al., 2007; Paissé et al., 2008; Hazen et al., 2010; Bælum et al., 2012; Lu et al., 2012; Mason et al., 2012), which has concomitant impacts on microbially-mediated biogeochemistry (Atlas, 1991). Proliferation of heterotrophs in the presence of hydrocarbons as the sole electron donor can result in the rapid depletion of oxygen in the ecosystem through aerobic respiration (Riser-Roberts, 1992).

One potential limiting factor for hydrocarbon degradation in marine sediments is nitrogen availability (Herbert, 1999). Marine sediments are a key location for global nitrogen cycling (Engström et al., 2009; Lam and Kuypers, 2011; Hamme and Emerson, 2013), because they provide both a physical location for organic matter remineralization through ammonification, and also an anoxic environment for denitrification processes (Laverock et al., 2011). At the water-sediment interface, ammonium generated through organic matter remineralization, is converted to nitrite and then nitrate by nitrification (Henriksen and Kemp, 1988). Available nitrite or nitrate may diffuse to the water column or be consumed in various biotic reactions depending on oxygen availability. Under oxygen limited conditions, deeper within the sediment, nitrate may be reduced via dissimilatory nitrate reduction to ammonium (DNRA), which can be anaerobically oxidized, with nitrite, to dinitrogen, through anammox processes, also the experimental evidence for this mechanism remains limited (Brandes et al., 2007; Kalvelage et al., 2011). It has been estimated that dinitrogen production in deep benthic sediments (1000–3000 m) is responsible for ~7–16% of the total nitrogen loss in the marine ecosystem (Engström et al., 2009).

How oil contamination affects nitrogen cycling processes in situ is still not well understood (Deni and Penninckx, 1999; Bell et al., 2011; Trimmer et al., 2013). Although petroleum oil is typically rich in nitrogen, most of it is bound in aromatic heterocyclic compounds whose carbon to nitrogen bonds are difficult to break (Snyder, 1970). Thus, bioremediation of oil contamination often requires the addition of inorganic nutrients including nitrogen, phosphorus, and/or iron, to increase enzyme activity (Brook et al., 2001; Head et al., 2006; Bell et al., 2011). In the case of high carbon and low nitrogen environments, there is evidence of increased diazotrophy within microbial communities (Karl et al., 2002), although this has rarely been observed in the presence of hydrocarbons despite nitrogen limiting conditions (Laguerre et al., 1987; Eckford et al., 2002; Musat et al., 2006). A more recent study exploring the relationship between microbial nitrogen cycling dynamics and oil contamination found evidence of DNRA being coupled to nitrogen cycling in a suboxic hydrocarbon contaminated subsurface well (Yagi et al., 2010).

A number of studies have evaluated how well the relative abundance of genes encoding nitrogen metabolic enzymes correlate with biogeochemical measurements of N metabolism. For instance, the abundance of nirS, nrfA, narG, and napA genes involved in nitrite and nitrate reduction significantly predicted the rates of denitrification and DNRA in the Colne estuary (Dong et al., 2009). However, nirK, another nitrite reduction gene, has been shown to be a poor predictor of functional traits relevant for denitrification (Salles et al., 2012). The metabolic pathways themselves are a network of interactions, so we hypothesize that the non-linear relationship between gene abundance and metabolite turnover is best evaluated as a function of compound changes in the relative abundance of many different related genes, rather than any single gene abundance. Predicted relative metabolic turnover (PRMT) quantifies relative changes in the metabolic potential as a network of predicted metabolic reactions inferred from the relative abundance of genes annotated from a metagenome (Larsen et al., 2011). Each predicted metabolite is then a function of the predicted enzymes and their “metabolic community of reactions” rather than simply the relative abundance of just the single gene that codes for the enzyme responsible for the metabolism. PRMT has previously been used to accurately predict seasonal variation in metabolites in marine surface waters (Larsen et al., 2011).

Here we hypothesize that increasing exposure to elevated hydrocarbon concentrations will positively correlate with predicted metabolic shifts toward denitrification in anaerobic zones in sediments. This is based on the premise that previously challenged and constantly exposed hydrocarbon samples are more likely to be “primed” for hydrocarbon response (Deni and Penninckx, 1999; Labbé et al., 2007; Taketani et al., 2010). To test this hypothesis we analyzed metagenomic sequence data from 17 sediment samples from the Gulf of Mexico and the Santa Barbara Channel, which represent sites with short term exposure to oil contamination [those from the Deepwater Horizon (DWH) oil spill], sites with a long history of exposure to hydrocarbons (those from the natural oil seeps), and sites unaffected by hydrocarbon contamination.

Materials and Methods


Data were downloaded from MG-RAST (Meyer et al., 2008) including the DWH spill project (MGRAST IDs: 4510162.3-4510175.3) and the natural oil seeps study (MGRAST IDs: 4537092.3-4537094.3), which were all annotated using SEED, with a maximum e-value of 1 × 10−3, a minimum identity of 50%, and minimum identity cutoff of 15. Data was also annotated in MG-RAST using Hierarchical Classification subsystems with a maximum e-value cutoff of 10−5, minimum percent identity cutoff of 60%, and a minimum alignment length cutoff of 15; this data was used for looking at functional annotations for gene abundances. For the samples collected from the Gulf of Mexico during the DWH spill, 6 of the samples were from oil-contaminated sites (hereafter referred to as DWH oil spill samples), and 8 samples were from uncontaminated sites (hereafter referred to as uncontaminated samples). This grouping was based on whether the samples clustered based on the normalized gene abundances, and additionally based on whether they exceeded (>1.0 polycyclic aromatic hydrocarbon [PAH] index) or did not exceed EPA (≤1 PAH index) BPA benchmarks for hydrocarbon pollution (Mason et al., 2014) (for more information about how the EPA aquatic benchmarks are calculated please see The Santa Barbara channel oil seep samples included depth, latitude and longitude, and collection date as contextual metadata. This data is summarized in Supplementary Table 1. The Gulf of Mexico samples (Mason et al., 2014) had these contextual metadata in addition to total petroleum hydrocarbons (TPH), polycyclic aromatic hydrocarbons (PAH), dissolved-phosphate (PO4-P), dissolved nitrate (NO3-N), total ammonia nitrogen (NH3-N and NH4-N), dissolved inorganic nitrogen (DIN; NH3-N and NH4-N), total nitrogen (NH3/NH4-N, NO3/NO2-N, and organic nitrogen), total sulfur (S), and total carbon (C). A complete metadata table for the Gulf of Mexico samples is given in Supplementary Table 1 of Mason et al. (2014). For more information about sample collection and the context of these samples please see Hawley et al. (2014), Mason et al. (2014). These values were normalized and log2 transformed before analysis was performed.


The oil seep samples had genetic sequences that annotated to 131 nitrogen metabolism genes that were not present in any of the Gulf samples. Thus for gene annotations, SEED Subsystems-based functional (level 2) annotations were summed and then standardized as a function of total reads within each sample. Predicted Metabolic Turnover Analysis (PRMT) (Larsen et al., 2011) was used to evaluate the community metabolic potential between samples as a function of microbial community gene abundances. PRMT transforms annotated enzyme abundances by a weighted matrix of all possible reactions including those enzymes, their reactions, and associated metabolites as annotated by KEGG (Ogata et al., 1999). Enzyme commission (EC) abundances were gathered from the SEED Subsystems L3 tables, quantile normalized and then log2 transformed before analysis. The EC abundances were compared to a “reference,” which in this analysis was an average of all samples. Positive PRMT score values represent the consumption of a particular metabolite, and negative scores represent the accumulation or production of a particular metabolite. For the nitrogen metabolism pathway (KEGG map00910), the PRMT scores were summed to give either a “net” positive or negative PRMT value. The “net difference” or “pathway flow” was found by adding the net positive and net negative values for PRMT scores for each metabolite in the pathway per sample. For comparisons of sample scores, Kruskal-Wallis rank sum tests were used. Hierarchical annotation for gene abundances were also quantile normalized and log2 transformed. A principal component analysis was performed on the quantile normalized and log2 transformed hierarchical abundances, removing those genes completely absent in the Gulf of Mexico dataset. To pull out hydrazine related gene sequences, bowtie 1.0.0 (Langmead et al., 2009) was used to align reads to custom index of hydrazine hydratase related sequences downloaded from NCBI (Benson et al., 2013). Then we used reads per kilobase per million (RPKM), quantile normalization, and a log2 transformation to normalize the hydrazine hydratase related gene abundances. Pearson correlation coefficients (corr) were used where cited. 10,000 permutations were used to assess empirical p-values.


The contaminated sediments from the Gulf of Mexico were exposed to hydrocarbon contamination from the DWH spill for between 3 and 5 months at the time they were collected. By contrast, the natural oil seeps from the Santa Barbara Channel were estimated to be exposed to hydrocarbons for more than 11,000 years (Hornafius et al., 1999). To the best of our knowledge the uncontaminated sediments from the Gulf of Mexico had not recently been exposed to the amounts of hydrocarbon contamination caused by the DWH Oil spill, although historical presence of temporary natural oil seeps nearby cannot be ruled out (Kvenvolden and Cooper, 2003). The metagenomic data were all generated using the HiSeq2000 platform, with a minimum of 36,851,796 reads and a maximum of 86,321,188 reads per metagenome, with read lengths of ~150 bp.

Oil Seep Samples Maintained a Greater Diversity of Genes Associated with Nitrogen Metabolism

Oil seep samples had sequences that annotated to 131 nitrogen metabolism genes; these genes were not present in the samples from the Gulf of Mexico (oil spill and uncontaminated sites). Of the 11 SEED level 2 annotations within nitrogen metabolism, 4 were present only in the oil seep samples- these included amidase clustered with urea and nitrile hydratase functions, cyanate hydrolysis, citric oxide synthase, and nitrilase.

Nitrogen Metabolism Genes Showed Differential Relative Abundances Between the 3 Different Sample Groups

Anammox pathway specific genes related to hydrazine production (an intermediate of the anammox reaction pathway) were not significantly different between oil seep, oil spill and uncontaminated samples; however they trended toward higher abundance in oil seep samples, followed by oil spill samples and uncontaminated samples. Interestingly, nitrosative stress, was found only in the petroleum-contaminated groups (both oil spill and oil seep), and is involved in response to nitric oxide accumulation (Ridnour et al., 2004). Coincidentally, nitric oxide is predicted to be accumulated by the PRMT analysis (Table 2). Interestingly, only 3 of the 11 SEED level 2 nitrogen pathway annotations showed a significantly different relative abundance between sample types, including dissimilatory nitrite reductase, nitrate and nitrite ammonification, and nitrogen fixation (Table 1, Figure 1).


Table 1. Average proportion of reads annotation to each SEED subsystems-based functional annotation (level 2) function from nitrogen metabolism with their standard deviations in parentheses.


Figure 1. Graphical representation of the changing relative abundance of different key nitrogen pathways and genes within each of the sample types.

A principal component analysis was performed using only the 33 genes that were present across all sample types, with the SEED level 2 hierarchical gene annotation used to demonstrate which genes showed the greatest differentiation between sample types (Figure 2). The first two principal components account for 72.78% of the variance, and while the variance and the distribution of the sample types was due to multiple factors, the influence of key genes in differentiation of sample types was evident. For example, an abundance of nitrogen regulatory protein P-II (glnB) in seep samples, respiratory nitrate reductase (alpha, beta, delta, and gamma chain; narG, narH, narW and narI, respectively) and nitrogenase (alpha and beta chain; nifA and nifB, respectively) in DWH spill samples, and Cu-nitrite reductase (nirK), AnfO protein (anfO), and allantoate amidohydrolase (allC) in uncontaminated sediments played a considerable role in the differentiation of the 3 sample types (Figure 2).


Figure 2. Principal Components Analysis (PCA) of the relative abundance of nitrogen pathway gene annotation shared (33 genes) between the three samples types [oil spill (circles), oil seep (diamonds), and uncontaminated (triangles)]. Genes on the graph are given abbreviated names, but their full names and associated functional groups are given in Supplementary Table 2. The gene abundance loading values are given by the bottom and left axes. Samples from each group shown as shapes in gray with axes scores given by the top and right axes. Gene functions are colored by their level 2 annotation as specified at the right outer margin of the graph. The variance accounted for by each component is listed on each axis.

The Predicted Relative Turnover for Nitrogen Metabolites was Significantly Different Between the 3 Groups

The relative abundance of individual key genes showed differential responses (Figures 1, 2) across the gradient of hydrocarbon exposure time (seep>DWH spill>uncontaminated); therefore PRMT was used to infer how these differential relative abundances could be combined to predict relative metabolite turnover for the different nitrogen pathways. The KEGG nitrogen metabolism reference pathway includes 21 metabolites. The PRMT scores were used to infer whether these metabolites were relatively consumed or accumulated in each group; the PRMT score is positive if the metabolite is being consumed, and negative if it is being accumulated, with the magnitude an indication of the relative level of this metabolism. The overall “pathway flow” (the difference between predicted “net consumption” and “net production” for all metabolites in a pathway) was positive for the three sample types. A positive net pathway flow suggests that overall more nitrogen metabolites were being consumed than accumulated. The difference in “pathway flow” between the sample types was not significantly different; however, the highest “pathway flow” score was found in the oil seeps (PRMT “net difference” score [PRMTdiff] = 11.35), followed by the DWH spill samples (PRMTdiff = 10.77) and uncontaminated sediments (PRMTdiff = 9.56). When the average “net accumulation” (negative PRMT) scores were summed per group, the difference between sample types was significantly different (p = 0.02) and 2-fold higher in the oil seep samples (PRMT “net production” [PRMTp] = −14.39) compared to the DWH spill (PRMTp = −8.08) and uncontaminated sediments (PRMTp = −7.15). For the specific metabolites in the pathway, there were also numerous significant differences between sample types (Table 2). Specifically, nitrate was predicted to be significantly more consumed in the oil seep (mean = 5.1, SD = 2.73) and DWH spill samples (mean = 2.98, SD = 3.41), compared to uncontaminated sediments, where it was predicted to be more significantly accumulated (mean = −0.75, SD = 0.73; p-value < 0.01). Meanwhile, nitrite was predicted to be significantly more accumulated in the oil seep (mean = −3.28, SD = 0.64) and DWH spill samples (mean = −1.99, SD = 1.67), while being relatively consumed in the uncontaminated sediments (mean = 0.37, SD = 0.58; p-value < 0.01). Genes annotating to nitrosative stress were only found in the hydrocarbon contaminated sediments, which is supported by the prediction that nitric oxide was significantly more accumulated in both the oil seep (mean = −2.3, SD = 1.27) and oil spill (mean = −1.1, SD = 0.31) compared to the uncontaminated sediments where it was relatively consumed (mean = 0.45, SD = 0.61). There was also a relative increase in the consumption of ammonia in the oil seep group, although this difference was not statistically significant (Table 1). In addition, nitrile was predicted to be significantly consumed (p-value < 0.001) in the oil seep (mean = 2.11, SD = 0.59) and the DWH spill samples (0.72, SD = 0.14), while being relatively accumulated in the uncontaminated sediments (mean = −0.34, SD = 0.33).


Table 2. Nitrogen metabolism associated metabolites' mean (and standard deviation in parentheses) PRMT scores.

PRMT Scores for Nitrogen Pathway Metabolites Show Significant Correlations with in Situ Biogeochemical Measurements Between Oil Spill and Uncontaminated Sediment Samples in the Gulf of Mexico

The samples collected from the Gulf of Mexico were analyzed in more detail for significant correlations to the available biochemical data (Mason et al., 2014). PRMT scores for dinitrogen showed a significant positive correlation with measured concentrations of in situ total nitrogen (p < 0.05, corr = 0.55). In addition, a number of the other PRMT scores had significant correlations with total sulfur, total carbon, total nitrogen, dissolved nitrate, total ammonium, dissolved inorganic nitrogen, and dissolved phosphate (Table 3). The PRMT scores for nitrite had a significant negative correlation with total carbon (p < 0.01, corr = −0.58), which suggests that when there is more carbon there is a significant accumulation of nitrite. In addition, L-aspartate had a significant correlation (p < 0.05, corr = 0.53) with total hydrocarbon concentration.


Table 3. Pearson correlations of metadata from 14 metagenome samples (Mason et al., 2014) with predicted relative nitrogen metabolite PRMT scores and their associated p-values.


Here we present evidence of the impact of oil contamination, including comparisons of short-term vs. long-term duration of exposure, on nitrogen metabolism in marine sediments. Oil contaminated and uncontaminated sediment samples collected after the DWH spill in the Gulf of Mexico were compared to samples collected from natural oil seeps from the Santa Barbara Channel. Genes and pathways involved in the nitrogen cycle were annotated from metagenomic sequencing data and used to explore differences in the relative abundance of specific genes and to predict relative nitrogen metabolite turnover potential between the 3 sample types. These sample types come from disparate regions (e.g., Gulf of Mexico vs. Santa Barbara Channel), thus numerous other geochemical and physical factors could have played a role in the observed trends in nitrogen metabolism between these environments. However, this study suggests that the selective pressure of oil contamination contributes a significant role toward shaping the functional diversity of these community processes. In addition, we expand on an analysis of metagenome data (Mason et al., 2014) and show that this analysis can be useful for exploring the impacts of hydrocarbon contamination on nitrogen cycling in other contaminated environments.

Studies of the relative abundances of specific genes may not be the best way to study complex, multi-branching metabolic pathways. To overcome this limitation, we used PRMT to better capture the emergent property of the multiphasic gene abundance profiles that make up a metabolic pathway. The PRMT approach captures the relative metabolic changes across an observed assemblage of genes, and therefore the relative abundances of genes and their corresponding metabolic pathways are taken in proportion to each other (Larsen et al., 2011). The predicted “net pathway flow” suggests that overall more nitrogen metabolites were consumed in each sample type than were accumulated, with this value being greatest in the oil seep samples. While this might seem to infer a system mass balance bias, as the data used for predictions is static, these inferences cannot be used to infer mass potential. For those metabolites that are predicted to be accumulated, there was a two-fold increase in the oil seep samples compared to DWH spill and uncontaminated sediment samples. In addition, the specific metabolites that were predicted to accumulate in the contaminated samples were different from those in the uncontaminated samples, which may represent shifts in nitrogen cycling processes in sediments exposed to hydrocarbon saturation.

The metabolites that were significantly different between the three groups, i.e., nitrate, nitrite, and nitric oxide had a common trend in which the oil seep samples had the highest consumption and accumulation, followed by the DWH spill samples, and finally the uncontaminated samples, where the values were often close to 0; suggesting that the pathways involved in consumption and accumulation of nitrogen were balanced. Dinitrogen was an exception in that the uncontaminated sediments had a three-fold higher predicted consumption than in the contaminated samples. This was supported by the relative abundance of genes involved in diazotrophy (nitrogen fixation), which were most abundant in uncontaminated sediments.

The sum of predicted metabolite turnover scores for all nitrogen metabolism pathways evidenced an increase in denitrification processes either through canonical denitrification or anammox, as nitrite and nitric oxide were predicted to significantly accumulate and nitrate was predicted to be consumed by microbial metabolism in contaminated samples. It is more likely that this evidence could be interpreted as relating to canonical denitrification, despite the relatively higher abundance of anammox pathway specific genes related to hydrazine production and ammonium assimilation in seep and spill compared to uncontaminated samples.

The oil spill sediments from the Gulf of Mexico were collected ~3 months after the Deepwater Horizon's Macondo well was capped, thus giving them an active exposure time between 3 and 5 months, if we assume the absence of natural oil seeps near these sites. By contrast, the oil seep samples from the Santa Barbara Channel were actively exposed to petroleum for more than 11,000 years (Hornafius et al., 1999) and samples were taken directly from the seep head. The difference in the time of exposure to hydrocarbon pollution in the oil-contaminated sediments could thus account for differences in the predicted turnover of nitrogen metabolites. Additionally, significant differences in the composition of the oil from both sites (Hornafius et al., 1999; Reddy et al., 2011), may also have influenced the observed compositional differences the microbial communities (Hawley et al., 2014; Mason et al., 2014). Depth of sample collection from the two different environments may also have affected nitrogen cycling as shown in several studies (Engström et al., 2009; Trimmer et al., 2013), however, this trend is probably due to differences in the physicochemical properties of different sites, as has been shown for sites at different distances from shore (Herbert, 1999; Dalsgaard et al., 2005; Zhu et al., 2010). Despite these geographic, depth, and oil composition differences, there were surprising similarities in the response of metabolic turnover to hydrocarbon contamination, suggesting that oil contamination results in a predictable metabolic response despite differences in the affected ecosystems. It is suggested that a topic for future research might be using PRMT on genetic data generated from oil exposed laboratory enrichments or environmental samples through time, to aid in unraveling the relationship between nitrogen cycling and microbial oil remediation.

Nitrite consumption showed a significant negative correlation to the concentration of total carbon in the Gulf of Mexico sediments. The decreased availability of reactive carbon and a high concentration of organic carbon in extremely deep benthic environments would favor the dominance of anammox over denitrification (Thamdrup and Dalsgaard, 2002; Engström et al., 2005), which would therefore lead to direct oxidation of ammonia to dinitrogen reducing nitrite consumption. This would however lead to a nitrogen limited environment, which could be supplanted by the oxidation of organic matter by sulfur reducing bacteria (Canfield et al., 2010); this is potentially supported by significant correlations between the concentration of sulfur and the accumulation of nitrate and nitrite.

To summarize, there is evidence from the PRMT analysis for a shift in the metabolic flow of nitrogen to the denitrification pathways, potentially including the anammox pathway, in hydrocarbon-contaminated sediments (both DWH spill and natural oil seep). Changes in metabolites in the anammox pathway were positively correlated with hydrocarbon concentration, although these were not statistically significant, potentially due to the small sample sizes and confounding environmental factors. The relative abundance of genes related to anammox associated hydrazine metabolism were also greatest in the seep samples that were predicted to have been exposed to hydrocarbons for ~11,000 years.

Marine sediments are very important sites for microbially mediated nitrogen transformation, providing a link between organic matter degradation and nutrient regeneration, essentially supporting primary productivity in the oceans. Exploring factors that significantly influence this process are vital for providing relevant data to propagate system scale models of how basin processes, such as nitrogen cycling and primary productivity in marine sediments, can influence regional and global climate.

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.


This work was supported in part by the US Department of Energy under Contract DE-AC02-06CH11357 and under Contract DE-AC02-05CH11231, and BP/The Gulf of Mexico Research Initiative. The sequence data from the Santa Barbara oils were produced by the US Department of Energy Joint Genome Institute in collaboration with the user community. We also thank Ali Shojaie for his advice about missing data problems in statistical analyses and Ben Rathbone for analytical discussions.

Supplementary Material

The Supplementary Material for this article can be found online at:


Atlas, R. M. (1991). Microbial hydrocarbon degradation—bioremediation of oil spills. J. Chem. Technol. Biotechnol. 52, 149–156. doi: 10.1002/jctb.280520202

CrossRef Full Text

Bælum, J., Borglin, S., Chakraborty, R., Fortney, J. L., Lamendella, R., Mason, O. U., et al. (2012). Deep-sea bacteria enriched by oil and dispersant from the Deepwater Horizon spill. Environ. Microbiol. 14, 2405–2416. doi: 10.1111/j.1462-2920.2012.02780.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bell, T. H., Yergeau, E., Martineau, C., Juck, D., Whyte, L. G., and Greer, C. W. (2011). Identification of nitrogen-incorporating bacteria in petroleum-contaminated arctic soils by using [15N]DNA-based stable isotope probing and pyrosequencing. Appl. Environ. Microbiol. 77, 4163–4171. doi: 10.1128/AEM.00172-11

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Benson, D. A., Cavanaugh, M., Clark, K., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J., et al. (2013). GenBank. Nucleic Acids Res. 41, D36–D42. doi: 10.1093/nar/gks1195

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bordenave, S., Goñi-Urriza, M. S., Caumette, P., and Duran, R. (2007). Effects of heavy fuel oil on the bacterial community structure of a pristine microbial mat. Appl. Environ. Microbiol. 73, 6089–6097. doi: 10.1128/AEM.01352-07

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Brandes, J. A., Devol, A. H., and Deutsch, C. (2007). New developments in the marine nitrogen cycle. Chem. Rev. 107, 577–589. doi: 10.1021/cr050377t

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Brook, T. R., Stiver, W. H., and Zytner, R. G. (2001). Biodegradation of diesel fuel in soil under various nitrogen addition regimes. Soil Sediment Contam. Int. J. 10, 539–553. doi: 10.1080/20015891109428

CrossRef Full Text

Canfield, D. E., Glazer, A. N., and Falkowski, P. G. (2010). The evolution and future of earth's nitrogen cycle. Science 330, 192–196. doi: 10.1126/science.1186120

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Dalsgaard, T., Thamdrup, B., and Canfield, D. E. (2005). Anaerobic ammonium oxidation (anammox) in the marine environment. Res. Microbiol. 156, 457–464. doi: 10.1016/j.resmic.2005.01.011

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Deni, J., and Penninckx, M. J. (1999). Nitrification and autotrophic nitrifying bacteria in a hydrocarbon-polluted soil. Appl. Environ. Microbiol. 65, 4008–4013.

Pubmed Abstract | Pubmed Full Text

Dong, L. F., Smith, C. J., Papaspyrou, S., Stott, A., Osborn, A. M., and Nedwell, D. B. (2009). Changes in benthic denitrification, nitrate ammonification, and anammox process rates and nitrate and nitrite reductase gene abundances along an estuarine nutrient gradient (the Colne Estuary, United Kingdom). Appl. Environ. Microbiol. 75, 3171–3179. doi: 10.1128/AEM.02511-08

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Eckford, R., Cook, F. D., Saul, D., Aislabie, J., and Foght, J. (2002). Free-living heterotrophic nitrogen-fixing bacteria isolated from fuel-contaminated Antarctic soils. Appl. Environ. Microbiol. 68, 5181–5185. doi: 10.1128/AEM.68.10.5181-5185.2002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Engström, P., Penton, C. R., and Devol, A. H. (2009). Anaerobic ammonium oxidation in deep-sea sediments off the Washington margin. Limnol. Oceanogr. 54, 1643–1652. doi: 10.4319/lo.2009.54.5.1643

CrossRef Full Text

Engström, P., Dalsgaard, T., Hulth, S., and Aller, R. C. (2005). Anaerobic ammonium oxidation by nitrite (anammox): implications for N2 production in coastal marine sediments. Geochim. Cosmochim. Acta 69, 2057–2065. doi: 10.1016/j.gca.2004.09.032

CrossRef Full Text

Gutierrez, T., Singleton, D. R., Berry, D., Yang, T., Aitken, M. D., and Teske, A. (2013). Hydrocarbon-degrading bacteria enriched by the Deepwater Horizon oil spill identified by cultivation and DNA-SIP. ISME J. 7, 2091–2104. doi: 10.1038/ismej.2013.98

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hamme, R. C., and Emerson, S. R. (2013). Deep-sea nutrient loss inferred from the marine dissolved N2/Ar ratio. Geophys. Res. Lett. 40, 1149–1153. doi: 10.1002/grl.50275

CrossRef Full Text

Hawley, E. R., Piao, H., Scott, N. M., Malfatti, S., Pagani, I., Huntemann, M., et al. (2014). Metagenomic analysis of microbial consortium from natural crude oil that seeps into the marine ecosystem offshore Southern California. Stand. Genomic Sci. 9, 635–650. doi: 10.4056/sigs.5029016

CrossRef Full Text

Hazen, T. C., Dubinsky, E. A., DeSantis, T. Z., Andersen, G. L., Piceno, Y. M., Singh, N., et al. (2010). Deep-sea oil plume enriches indigenous oil-degrading bacteria. Science 330, 204–208. doi: 10.1126/science.1195979

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Head, I. M., Jones, D. M., and Röling, W. F. M. (2006). Marine microorganisms make a meal of oil. Nat. Rev. Microbiol. 4, 173–182. doi: 10.1038/nrmicro1348

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Henriksen, K., and Kemp, W. M. (1988). “Nitrification in estuarine and coastal marine sediments: methods, patterns and regulating factors,” in Nitrogen Cycling in Coastal Marine Environments, eds T. H. Blackburn and J. Sørensen (Hoboken, NJ: John Wiley and Sons), 207–250.

Herbert, R. A. (1999). Nitrogen cycling in coastal marine ecosystems. FEMS Microbiol. Rev. 23, 563–590. doi: 10.1111/j.1574-6976.1999.tb00414.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hornafius, J. S., Quigley, D., and Luyendyk, B. P. (1999). The world's most spectacular marine hydrocarbon seeps (Coal Oil Point, Santa Barbara Channel, California): quantification of emissions. J. Geophys. Res. Oceans 104, 20703–20711. doi: 10.1029/1999JC900148

CrossRef Full Text

Kalvelage, T., Jensen, M. M., Contreras, S., Revsbech, N. P., Lam, P., Günter, M., et al. (2011). Oxygen sensitivity of anammox and coupled N-cycle processes in oxygen minimum zones. PLoS ONE 6:e29299. doi: 10.1371/journal.pone.0029299

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Karl, D., Michaels, A., Bergman, B., Capone, D., Carpenter, E., Letelier, R., et al. (2002). Dinitrogen fixation in the world's oceans. Biogeochemistry 57–58, 47–98. doi: 10.1023/A:1015798105851

CrossRef Full Text

Kvenvolden, K. A., and Cooper, C. K. (2003). Natural seepage of crude oil into the marine environment. Geo-Mar. Lett. 23, 140–146. doi: 10.1007/s00367-003-0135-0

CrossRef Full Text

Labbé, D., Margesin, R., Schinner, F., Whyte, L. G., and Greer, C. W. (2007). Comparative phylogenetic analysis of microbial communities in pristine and hydrocarbon-contaminated Alpine soils. FEMS Microbiol. Ecol. 59, 466–475. doi: 10.1111/j.1574-6941.2006.00250.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Laguerre, G., Bossand, B., and Bardin, R. (1987). Free-living dinitrogen-fixing bacteria isolated from petroleum refinery oily sludge. Appl. Environ. Microbiol. 53, 1674–1678.

Pubmed Abstract | Pubmed Full Text

Lam, P., and Kuypers, M. M. M. (2011). Microbial nitrogen cycling processes in oxygen minimum zones. Annu. Rev. Mar. Sci. 3, 317–345. doi: 10.1146/annurev-marine-120709-142814

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25. doi: 10.1186/gb-2009-10-3-r25

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Larsen, P. E., Collart, F. R., Field, D., Meyer, F., Keegan, K. P., Henry, C. S., et al. (2011). Predicted Relative Metabolomic Turnover (PRMT): determining metabolic turnover from a coastal marine metagenomic dataset. Microb. Inform. Exp. 1, 4. doi: 10.1186/2042-5783-1-4

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Laverock, B., Gilbert, J. A., Tait, K., Osborn, A. M., and Widdicombe, S. (2011). Bioturbation: impact on the marine nitrogen cycle. Biochem. Soc. Trans. 39, 315–320. doi: 10.1042/BST0390315

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lu, Z., Deng, Y., Van Nostrand, J. D., He, Z., Voordeckers, J., Zhou, A., et al. (2012). Microbial gene functions enriched in the Deepwater Horizon deep-sea oil plume. ISME J. 6, 451–460. doi: 10.1038/ismej.2011.91

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mason, O. U., Hazen, T. C., Borglin, S., Chain, P. S. G., Dubinsky, E. A., Fortney, J. L., et al. (2012). Metagenome, metatranscriptome and single-cell sequencing reveal microbial response to Deepwater Horizon oil spill. ISME J. 6, 1715–1727. doi: 10.1038/ismej.2012.59

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mason, O. U., Scott, N. M., Gonzalez, A., Robbins-Pianka, A., Bælum, J., Kimbrel, J., et al. (2014). Metagenomics reveals sediment microbial community response to Deepwater Horizon oil spill. ISME J. doi: 10.1038/ismej.2013.254. [Epub ahead of print].

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meyer, F., Paarmann, D., D'Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server – a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Musat, F., Harder, J., and Widdel, F. (2006). Study of nitrogen fixation in microbial communities of oil-contaminated marine sediment microcosms. Environ. Microbiol. 8, 1834–1843. doi: 10.1111/j.1462-2920.2006.01069.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ogata, H., Goto, S., Sato, K., Fujibuchi, W., Bono, H., and Kanehisa, M. (1999). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 27, 29–34. doi: 10.1093/nar/27.1.29

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Paissé, S., Coulon, F., Goñi-Urriza, M., Peperzak, L., McGenity, T. J., and Duran, R. (2008). Structure of bacterial communities along a hydrocarbon contamination gradient in a coastal sediment. FEMS Microbiol. Ecol. 66, 295–305. doi: 10.1111/j.1574-6941.2008.00589.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Reddy, C. M., Arey, J. S., Seewald, J. S., Sylva, S. P., Lemkau, K. L., Nelson, R. K., et al. (2011). Composition and fate of gas and oil released to the water column during the Deepwater Horizon oil spill. Proc. Natl. Acad. Sci. U.S.A. 20229–20234. doi: 10.1073/pnas.1101242108

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ridnour, L. A., Thomas, D. D., Mancardi, D., Espey, M. G., Miranda, K. M., Paolocci, N., et al. (2004). The chemistry of nitrosative stress induced by nitric oxide and reactive nitrogen oxide species. Putting perspective on stressful biological situations. Biol. Chem. 385, 1–10. doi: 10.1515/BC.2004.001

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Riser-Roberts, E. (1992). Bioremediation of Petroleum Contaminated Sites. Boca Raton, FL: CRC Press.

Röling, W. F. M., Milner, M. G., Jones, D. M., Fratepietro, F., Swannell, R. P. J., Daniel, F., et al. (2004). Bacterial community dynamics and hydrocarbon degradation during a field-scale evaluation of bioremediation on a mudflat beach contaminated with buried oil. Appl. Environ. Microbiol. 70, 2603–2613. doi: 10.1128/AEM.70.5.2603-2613.2004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Röling, W. F. M., Milner, M. G., Jones, D. M., Lee, K., Daniel, F., Swannell, R. J. P., et al. (2002). Robust hydrocarbon degradation and dynamics of bacterial communities during nutrient-enhanced oil spill bioremediation. Appl. Environ. Microbiol. 68, 5537–5548. doi: 10.1128/AEM.68.11.5537-5548.2002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Salles, J. F., Le Roux, X., and Poly, F. (2012). Relating phylogenetic and functional diversity among denitrifiers and quantifying their capacity to predict community functioning. Front. Microbiol. 3:209. doi: 10.3389/fmicb.2012.00209.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Snyder, L. R. (1970). Petroleum nitrogen compounds and oxygen compounds. Acc. Chem. Res. 3, 290–299. doi: 10.1021/ar50033a002

CrossRef Full Text

Taketani, R. G., Franco, N. O., Rosado, A. S., and van Elsas, J. D. (2010). Microbial community response to a simulated hydrocarbon spill in mangrove sediments. J. Microbiol. 48, 7–15. doi: 10.1007/s12275-009-0147-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Thamdrup, B., and Dalsgaard, T. (2002). Production of N(2) through anaerobic ammonium oxidation coupled to nitrate reduction in marine sediments. Appl. Environ. Microbiol. 68, 1312–1318. doi: 10.1128/AEM.68.3.1312-1318.2002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Trimmer, M., Engström, P., and Thamdrup, B. (2013). Stark contrast in denitrification and anammox across the deep Norwegian trench in the Skagerrak. Appl. Environ. Microbiol. 79, 7381–7389. doi: 10.1128/AEM.01970-13

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yagi, J. M., Suflita, J. M., Gieg, L. M., DeRito, C. M., Jeon, C.-O., and Madsen, E. L. (2010). Subsurface cycling of nitrogen and anaerobic aromatic hydrocarbon biodegradation revealed by nucleic acid and metabolic biomarkers. Appl. Environ. Microbiol. 76, 3124–3134. doi: 10.1128/AEM.00172-10

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zhu, G., Jetten, M. S. M., Kuschk, P., Ettwig, K. F., and Yin, C. (2010). Potential roles of anaerobic ammonium and methane oxidation in the nitrogen cycle of wetland ecosystems. Appl. Microbiol. Biotechnol. 86, 1043–1055. doi: 10.1007/s00253-010-2451-4

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Keywords: nitrogen cycling, marine sediments, denitrification, microbial ecology, metagenomics, deepwater horizon oil spill, oil contamination, oil seeps

Citation: Scott NM, Hess M, Bouskill NJ, Mason OU, Jansson JK and Gilbert JA (2014) The microbial nitrogen cycling potential is impacted by polyaromatic hydrocarbon pollution of marine sediments. Front. Microbiol. 5:108. doi: 10.3389/fmicb.2014.00108

Received: 22 November 2013; Accepted: 03 March 2014;
Published online: 25 March 2014.

Edited by:

Andreas Teske, University of North Carolina at Chapel Hill, USA

Reviewed by:

Mike L. Dyall-Smith, Charles Sturt University, Australia
David Gregory Weissbrodt, ETH Zürich and Eawag, Switzerland

Copyright © 2014 Scott, Hess, Bouskill, Mason, Jansson and Gilbert. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Nicole M. Scott, Department of Ecology and Evolutionary Biology, University of Chicago, 1101 E. 57th Street Chicago, IL 360637, USA e-mail: