Abstract
A pronounced deficit of nitrogen (N) in the oxygen minimum zone (OMZ) of the Arabian Sea suggests the occurrence of heavy N-loss that is commonly attributed to pelagic processes. However, the OMZ water is in direct contact with sediments on three sides of the basin. Contribution from benthic N-loss to the total N-loss in the Arabian Sea remains largely unassessed. In October 2007, we sampled the water column and surface sediments along a transect cross-cutting the Arabian Sea OMZ at the Pakistan continental margin, covering a range of station depths from 360 to 1430 m. Benthic denitrification and anammox rates were determined by using 15N-stable isotope pairing experiments. Intact core incubations showed declining rates of total benthic N-loss with water depth from 0.55 to 0.18 mmol N m−2 day−1. While denitrification rates measured in slurry incubations decreased from 2.73 to 1.46 mmol N m−2 day−1 with water depth, anammox rates increased from 0.21 to 0.89 mmol N m−2 day−1. Hence, the contribution from anammox to total benthic N-loss increased from 7% at 360 m to 40% at 1430 m. This trend is further supported by the quantification of cd1-containing nitrite reductase (nirS), the biomarker functional gene encoding for cytochrome cd1-Nir of microorganisms involved in both N-loss processes. Anammox-like nirS genes within the sediments increased in proportion to total nirS gene copies with water depth. Moreover, phylogenetic analyses of NirS revealed different communities of both denitrifying and anammox bacteria between shallow and deep stations. Together, rate measurement and nirS analyses showed that anammox, determined for the first time in the Arabian Sea sediments, is an important benthic N-loss process at the continental margin off Pakistan, especially in the sediments at deeper water depths. Extrapolation from the measured benthic N-loss to all shelf sediments within the basin suggests that benthic N-loss may be responsible for about half of the overall N-loss in the Arabian Sea.
Introduction
The Arabian Sea is the semi-enclosed, north-western part of the Indian Ocean. Connected with the Red Sea and the Persian Gulf, it also receives discharge from some of the largest rivers in the world and is fringed by amongst the densest human populations. Although covering only 1% of the ocean surface, the Arabian Sea accounts for ∼5% of the global phytoplankton production, which has characteristic seasonal variability driven by two monsoons each year (Marra and Barber, ; Wiggert et al., ). Owing to the high seasonal surface production, high respiration in subsurface waters along with slow ventilation produces a pronounced oxygen minimum zone (OMZ) at depths between ∼100 and 1000 m. This OMZ is associated with a high nitrogen deficit (Codispoti et al., ; Deutsch et al., ) and a strong secondary nitrite maximum found at similar depths, which have been attributed to high pelagic N-loss activities therein (Naqvi, ; Naqvi et al., ; Ward et al., ). Due to its size, the Arabian Sea OMZ is assumed to be one of the biggest pelagic N-sinks, with annual estimated rates varying between ∼30 and 60 Tg N year−1 (Bange et al., ; Codispoti et al., ; Devol et al., ).
On the other hand, N-loss processes also occur in marine sediments. In fact, benthic N-loss is believed to contribute ∼50–70% of global oceanic N-loss (Codispoti et al., ; Galloway et al., ; Gruber, ). Because of the geographical configuration of the Arabian Sea, OMZ waters therein impinge on the sediments along the continental margins off the coasts of India, Pakistan as well as Oman. Consequently, any in situ N-transformations within the OMZ waters would undoubtedly affect the N-budget of the sediments, and vice versa. Nevertheless, despite the obvious importance of benthic N-loss in the Arabian Sea, benthic N-loss activities have hardly been assessed, and thus estimates of the benthic contribution to the N-deficit and overall N-loss in the Arabian Sea remain poorly constrained. Based on depth-integrated primary production rates (Seitzinger and Giblin, ), Bange et al. () estimated that shelf and margin sediments may account for 17% of the N-loss in the Arabian Sea; or up to 26% estimated by Schwartz et al. () from the changes in N2:Ar and nitrate consumption rates in Arabian Sea sediments. No direct measurements have been made, however, to distinguish the benthic N-loss pathways, nor have the potential interactions with overlying OMZ waters been much considered.
In general, two processes are known to remove nitrogen from marine systems: the N2O and N2 production via canonical denitrification and the N2 production by anaerobic ammonium oxidation (anammox, ). In marine environments, anammox activities were first detected in sediments (Dalsgaard and Thamdrup, ; Thamdrup and Dalsgaard, ), and later in the suboxic water columns of the Black Sea (Kuypers et al., ) and Golfo Dulce, Costa Rica (Dalsgaard et al., ). Since then anammox bacteria have been found in marine habitats ranging from the Arctic sea ice (Rysgaard et al., ) to deep sea hydrothermal vents (Byrne et al., ). In sediments, anammox has been shown to contribute up to 80% to the N2 production (Dalsgaard et al., ), but anammox rates measured by 15N stable isotope pairing experiments in sediments underlying a major OMZ to our knowledge have never been made before.
The reduction of nitrite to nitric oxide is an essential step in both anammox and denitrification, and is mediated by specific nitrite reductases (Nir; Schalk et al., ; Strous et al., ; Kartal et al., ). In general, two different types of nitrite reductases are known to occur, the copper-(NirK), and the cd1-containing nitrite reductase (NirS), but organisms harbor either of the reductases. The nirK genes are not only present in denitrifiers, but also known to occur in nitrifiers and therefore not suitable for the quantification of denitrifiers. Hence, genes encoding for the (nirS) are more commonly used as biomarkers for denitrifiers (Jayakumar et al., ; Castro-Gonzalez et al., ; Tiquia et al., ; Dang et al., ) and found to be more abundant in general and in an estuary system (Abell et al., ). Meanwhile, anammox bacteria also use a NirS, which is phylogenetically distinct from denitrifier NirS. Thus, nirS can be a useful biomarker to distinguish between denitrifiers and anammox bacteria, as evidenced by studies in the Peruvian and Arabian Sea pelagic OMZs (Lam et al., ).
In this study, we determined N-loss rates of denitrification and anammox in surface sediments at the continental margin off Pakistan via 15N-stable isotope experiments in intact core and slurry incubations. The relative abundances of denitrifying and anammox bacteria in the sediment were quantified based on their respective nirS genes and their phylogenies were further evaluated to characterize the benthic microbial communities at various station depths. In order to explore the potential interaction between benthic and OMZ N-loss rates, stations with water depths between 360 and 1430 m were sampled. Accordingly, sediments at one station lay below the OMZ, while the others were within OMZ waters.
Materials and Methods
Sampling procedures and chemical analyses
Sampling was conducted during the R/V Meteor cruise M74/2, on 7th to 28th October 2007, in the Arabian Sea over the Pakistan shelf (Makran region, Figure 1). Four stations ranging from 360 to 1430 m were selected for detailed sampling and sediment incubations. (Please note, that original station names have been shortened for simplicity, from, e.g. GeoB12204 to station 04). Dissolved oxygen and temperature of the water column were measured with a conductivity-temperature-depth (CTD) probe, equipped with an oxygen sensor (Sea Bird Electronics). The oxygen concentration was calibrated against Winkler titration. Water samples were taken with a CTD-rosette. On ship board, concentrations of ammonium and nitrite were measured fluorometrically (Holmes et al., ) and photometrically (Grasshoff and Johannsen, ), respectively. Additional subsamples were stored at −20°C for later analyses for ammonium, nitrate, nitrite, and phosphate in a shore based laboratory using an autoanalyzer (TRAACS 800, Bran & Luebbe).
Figure 1
Sediment cores were taken with a multicorer (MUC) equipped with eight acrylic liners (10 cm diameter). Subsamples for molecular analyses were taken directly from MUC cores at 2 cm intervals from the surface down to 8 cm. DNA samples were stored at −80°C, shipped on dry ice and kept at −80°C until DNA extraction. Pore water extraction was conducted on board, sediment cores were sliced in a resolution of 0.5 cm (sediment depth 0–1 cm), 1 cm (sediment depth 2–5 cm), and 4 cm (sediment depth 5–0 cm), and pore water was squeezed out of the sediment slices with a pore water press (Schlüter, ). Pore water samples for nitrate and nitrite were kept frozen until measurement with a chemiluminescence NOx analyzer (Thermo Environmental Instruments, Inc; Braman and Hendrix, ) in a shore based laboratory. For analyses of dissolved iron and sulfide, sediment cores were sampled on board with rhizomes at 1 cm (sediment depth 0–5 cm) and 2 cm intervals (sediment depth 5–30 cm). Subsequently the obtained pore water was analyzed for Fe2+ and HS− on board according to Grasshoff et al. () and Cline (). Concentrations of organic carbon and nitrogen were determined by combustion/gas chromatography (Carlo Erba NA-1500 CNS analyzer) of dried sediment samples after acidification with 3 mol l−1 phosphoric acid in a shore based laboratory.
Incubation experiments
Benthic denitrification and anammox rates were determined from N2 production of 15N-labeled slurry and intact core incubations. Rates from slurry incubations were used to calculate the contribution of anammox and denitrification to the total N-loss. Furthermore, volumetric rates from slurry incubations were integrated over the nitrate penetration depth to derive areal N-loss rates. Areal rates were also estimated from intact core incubations according to the revised isotope pairing technique (rIPT) detailed in Risgaard-Petersen et al. ().
Intact core incubations
Sediment cores (10 cm diameter) were subsampled with 3.6 cm diameter liners and the overlying water was adjusted to a height of 12.5 cm above sediment surface. (Campro Scientific GmbH) was added to a final concentration of 50 μmol l−1 in the overlying water, which was constantly mixed with magnetic stirrers. After pre-incubation for 8–12 h, the cores were sealed with rubber stoppers and incubated without gas phase in the dark at in situ temperature (6–16°C). Five time points were taken at ∼0, 2, 6, 10, and 15 h after the cores have been sealed. At each time point, three cores were randomly selected and sacrificed by first removing the rubber stopper and injecting 1 ml of 50% (w/v) zinc chloride to the overlying water to precipitate any free sulfide. Then the first 6 cm of the sediment were mixed with the overlying water. A subsample of the slurry was transferred into 12 ml gas tight sterile glass vials (Exetainer®, Labco), poisoned with 100 μl of saturated HgCl2 solution to stop biological activity and kept at room temperature in the dark until further processing.
Slurry incubations
Vertical distributions of denitrification and anammox rates within the sediment were estimated from slurry incubation experiments in gas tight bags made of plastic-laminated aluminum-foil (Gao et al., ). Briefly, MUC sediment cores were sliced in 2 cm intervals between 0 and 8 cm depth. Each slice (volume of ∼160 cm3) was transferred into a gas tight bag that was subsequently heat-sealed from all sides. To prepare the slurry, 200 ml of degassed bottom water, taken from the overlying water in the MUC cores, was injected through a gas tight port into the bag. The residual air was removed from the bag and the slurry was thoroughly mixed. After pre-incubating the bags for 2 h, to remove potential air-contamination introduced by the sub-sampling, 15N-labeled substrates were injected into the bags and the slurries were again thoroughly mixed. Incubations were performed in the dark at in situ temperatures. In Experiment 1, and were added to the slurries to final concentrations of 200 and 100 μmol l−1, respectively. Furthermore, allylthiourea (ATU) was added to a final concentration of 86 μmol l−1 (Ginestet et al., ) to inhibit possible bacterial ammonia oxidation. In Experiment 2, was added to the slurries to a final concentration of 200 μmol l−1. For both experiments, a subsample of 6 ml was drawn from the bags immediately after tracer addition, transferred into sterile gas tight glass vials (Exetainer®, Labco) and fixed with 100 μl of saturated HgCl2 solution. Between five and seven subsamples were drawn from the bags during the subsequent 26–28 h. The exetainers containing the subsamples were kept and shipped upside down in the dark at room temperature.
In the laboratory, a 2 ml helium headspace was introduced into the 12 ml exetainer of the whole core incubations while a headspace of 1 ml was used for the 6 ml exetainer of the slurries. The exetainers were shaken vigorously to allow N2 to equilibrate between the headspace and the liquid phase. The N2 isotope ratio (28N2, 29N2, and 30N2) of the headspace was determined by gas chromatography-isotopic ratio mass spectrometry (VG Optima, Micromass) by direct injections from the headspace according to Kuypers et al. (). Concentrations of 30N2 and 29N2 were normalized to 28N2 and calculated as excess relative to air according to Holtappels et al. (). N2 production rates were calculated from the 29N2 and 30N2 increase over time. Only production with a significant linear slope (p < 0.05) over time without delay was used for further calculations.
Calculation of N-loss in the sediment slurries
In Experiment 1, the anammox pathway combines either or with to form 28N2 and 29N2. Anammox activity was indicated when the production of 29N2 (p29N2) was measured without any production of 30N2 (p30N2). The production of 30N2 was not detected in our samples, only a small amount of 30N2 production was measured at station 16, depth 2–4 cm. The total N2 production via anammox in Experiment 1 [A(Ex1)] was calculated from: where is the labeling percentage of the 15N-substrate . For Experiment 1, was calculated from the measured -concentrations in bottom waters and pore waters and the known addition of
In Experiment 2, the addition of to background concentrations of and would produce 28N2 and 29N2 via anammox and 28N2, 29N2, and 30N2 via denitrification. Thus, the production of 30N2 (p30N2) indicates active denitrification. The total N2 production by denitrification in Experiment 2 [D(Ex2)] was calculated according to Thamdrup and Dalsgaard () from p30N2: where is the labeling percentage of nitrate In Experiment 2, both, anammox and denitrification produce 29N2. To calculate anammox from Experiment 2, Eq. 1 is modified to: where p29N2(Den) is the 29N2 production via denitrification. With (Thamdrup and Dalsgaard, ), we derive:
Results from slurry incubations indicated the presence of -storing organisms releasing intracellular in the course of the experiment (for further details, see discussion). An estimate of on the basis of measured bottom water and pore water concentrations was therefore not possible. Instead, we calculated from Eq. 3 by inserting the measured p29N2 and p30N2 and by assuming A(Ex1) = A(Ex2). The derived value, in the following referred to as was then used to estimate denitrification according to Eq. 2. For sediments without the release of stored nitrate we expected equals , whereas indicated an additional source of which was not dissolved initially in the pore water. We denoted the additional nitrate as excess that was calculated from , and the known concentration of in the slurry:
Calculation of N-loss in intact sediment cores
From the slurry incubation, the contribution of anammox to the total N-loss was estimated as ra = A(Ex1)/[A(Ex1) + D(Ex2)]. The total N-loss due to denitrification and anammox was calculated according to Risgaard-Petersen et al. () from ra and the production of 30N2 and 29N2 in the core incubations: where R29 is the ratio between the 29N2 and 30N2 production.
Detection and phylogenetic analyses of denitrifier and anammox nirS genes
The biomarker functional gene nirS, encoding the cd1-containing nitrite reductase, for both denitrifiers and marine anammox bacteria were targeted using qualitative and quantitative analyses. Nucleic acids were extracted from the sediment layers corresponding to those used for rate measurements (0–2, 2–4, 4–6, and 6–8 cm, respectively), by applying the UltraClean™ Soil DNA Isolation Kit (MO BIO Labratories, Inc.) according to the manufacturer’s instructions. Triplicate DNA extractions were made for each sample to reduce bias through the extraction procedure and pooled together through purification with the Wizard® Genomic DNA Purification Kit (Promega GmbH). DNA was stored in 10 mM Tris-HCl at −80°C until further analyses. The concentrations of the DNA in the samples were measured spectrophotometrically with a NanoDrop instrument (Thermo Fisher Scientific Inc.).
Denitrifier nirS gene fragments were PCR amplified with the primers cd3aF/R3cd (5′-GTSAACGTSAAGGARACSGG-3′ (Michotey et al., )/5′-GASTTCGGRTGSGTCTTGA-3′; Throback et al., ). The primers Scnir372F/Scnir845R (5′-TGTAGCCAGCATTGTAGCGT-3′/5′-TCAAGCCAGACCCATTTGCT-3′; Lam et al., ) were used to target the specific nirS for marine anammox bacteria, so far believed to all fall into the Candidatus Scalindua clade. PCR reactions were performed with the Master Taq Kit (5 Prime) on a thermal cycler (Eppendorf AG) and were examined with gel electrophoresis on 1% LE agarose gels (Biozym Scientific GmbH).
Subsequently, clone libraries were constructed from PCR amplicons of correct sizes. The PCR products were purified with the QIAquick PCR Purification Kit (QIAGEN) and the cloning reactions were performed with the TOPO TA Cloning® Kit for sequencing (pCR4 vector) with One Shot® TOP10 chemically competent E. coli cells (Invitrogen GmbH). Clones were screened for correct inserts by performing PCR with the primers M13F/M13R (5′-GTAAAACGACGGCCAG-3′/5′-CAGGAAACAGCTATGAC-3′), the number of non-nirS sequences was ≤2 per library. PCR products of the correct size were sequenced using the dye terminator sequencing method (Sanger et al., ) with the BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Inc.) and the T7 primer (5′-TAATACGACTCACTATAGGG-3′). Sequencing was performed on an ABI3730 capillary sequencer system (ABI) according to the manufacturer’s protocol. For the primers cd3aF/R3cd, sequence length used for phylogenetic analyses was ∼400 bp, while for primers Scnir372F/Scnir845R sequences had a minimum length of 440 bp.
Sequences were initially processed using BioEdit (Hall, ), aligned with ClustalW (Thompson et al., ) and screened for NirS encoding genes in the GenBank using the BLAST searches (Altschul et al., ). Mothur (Schloss et al., ) was used to calculate a similarity cut-off for operational taxonomic units (OTUs) based on nucleic acids of ≥95% and rarefaction curves. The screened sequences were imported into the ARB software package for phylogenetic analyses (Ludwig et al., ). Phylogenetic analyses were performed according to the amino acid sequences translated from the obtained sequences, together with some related sequences retrieved from GenBank. Phylogenetic trees were calculated based on the algorithms of maximum likelihood and maximum parsimony. Bootstrapped analyses of 100 resamplings were conducted. The sequences were deposited in GenBank under the accession numbers KC111208 to KC111421.
Quantitative PCR
Both denitrifier- and Scalindua-specific nirS genes were further quantified with real-time PCR, using the primers cd3aF/R3cd (Michotey et al., ) and Scnir372F/Scnir845R (Lam et al., ), which result in amplicons of 425 and 473 bp, respectively. The reactions were performed on an iQ5 cycler (Bio-Rad Laboratories GmbH) with the PowerSYBR® Green Master Mix (Applied Biosystems Inc.), as previously described (Lam et al., ; Jensen et al., ). All samples and non-template controls were analyzed as triplicate and the standards were analyzed in every qPCR run. The specificities of PCR amplicons were checked with subsequent melt curve analyses, as well as with 2% agarose gel electrophoresis.
Results
Hydro- and geochemistry
The compiled oxygen concentration profiles of the four investigated stations revealed an OMZ with a vertical expanse of ∼900 m (Figure 2). Within the oxycline (50–100 m), oxygen concentrations decreased from ∼200 to ∼5 μmol O2 l−1. From 200 to 300 m, an intrusion of Persian Gulf Water, identified by higher salinity (Shetye et al., ), led to increased oxygen concentrations of up to 16 μmol O2 l−1. At 300 m, oxygen concentrations dropped below the detection limit (∼1 μmol O2 l−1) and increased again below ∼900 m water depth. Bottom water oxygen concentrations of 23 μmol O2 l−1 were measured at the deepest station (1430 m), whereas no oxygen was detectable in the bottom water of the three shallower stations. Concentrations of ammonium were low throughout the water column (, data not shown). Nitrite concentrations were close to the detection limit of but peaked at distinct depths to maximum concentrations of at 30 m depth and 0.33 μmol l−1 between 400 and 600 m (Figure 2). Nitrate was depleted in the surface waters but increased below the oxycline (Figure 2) so that bottom water concentrations increased from at the shallowest station to at the deepest station. The nitrogen deficit, calculated according to Gruber and Sarmiento () as , was zero in surface waters (Figure 2), then decreased to −11 μmol N l−1 between 300 and 600 m depth and rose slightly to −7 μmol N l−1 below 800 m depth.
Figure 2
Within the sediments, the pore water was analyzed for the upper ∼30 cm (Figure 3). Nitrate concentrations in the first 0.5 cm of the sediment ranged from and dropped sharply to below. Similar to nitrate, nitrite generally declined within the upper centimeters from at the surface to below 2 cm. Significant subsurface maxima of nitrate ( at 17.5 cm, station 04) and nitrite (up to , stations 10 and 04) were sometimes found at the shallower stations. Concentrations of dissolved Fe2+ increased from 0.2 μmol l−1 at the surface to maximum concentrations ranging from 55 to 133 μmol l−1 at 5–9 cm depth and decreased within the layers below (Figure 3). Sulfide was detected only at the shallowest station 10 below 23 cm sediment depth where it increased with depth to a maximum of ∼400 μmol HS− l−1 at the lowermost sampled layer (31 cm, Figure 3).
Figure 3

Pore water profiles for nitrate, nitrite, iron, and sulfide for the investigated stations. The gray zone indicates the layers sampled for slurry incubations and DNA extraction, while the zone immediately above in the first panel represents the water column or bottom water. Samples for nitrate concentrations in the bottom water were retrieved from the bottom-most CTD sample. Sulfide was only measurable at and thus shown for station 10.
Organic carbon and nitrogen contents were measured in the sediment layers corresponding to the slurry incubations (Table 1). Within the OMZ, the organic carbon content (% of dry weight) in the surface sediment layer increased from 1.6% at 360 m to 2.4% at 970 m, but decreased again to 1.7% below the OMZ at 1430 m. Although organic carbon and nitrogen contents decreased within sediment depth at all stations, there was no clear trend for C:N ratios with sediment depth. However, the C:N ratios were slightly enhanced with station depth within the OMZ (C:N = 8–9 at 680 and 970 m), compared to the shallowest and deepest stations (C:N = 7–8).
Table 1
| Station | Sediment depth [cm] | Organic carbon [% dry wt] | Organic nitrogen [% dry wt] | C:N [mol:mol] | Excess nitrate [nmol (cm sed)−3] | DNA [ng DNA (mg dry sed)−1] |
|---|---|---|---|---|---|---|
| 10 | 0–2 | 1.6 | 0.26 | 7.2 | 95.8 | 4.92 ± 0.03 |
| 2–4 | 1.5 | 0.23 | 7.8 | 22.2 | 3.69 ± 0.23 | |
| 4–6 | 1.5 | 0.21 | 8.2 | 12.3 | 3.47 ± 0.23 | |
| 6–8 | 1.2 | 0.21 | 7.1 | n.d. | 3.05 ± 0.24 | |
| 04 | 0–2 | 2.2 | 0.31 | 8.4 | 110.6 | 9.18 ± 0.65 |
| 2–4 | 2.0 | 0.28 | 8.1 | 23.6 | 5.04 ± 0.59 | |
| 4–6 | 1.9 | 0.25 | 8.8 | 3.1 | 6.42 ± 0.84 | |
| 6–8 | 1.8 | 0.26 | 8.1 | n.d. | 3.96 ± 0.33 | |
| 16 | 0–2 | 2.4 | 0.31 | 9.0 | 222.0 | 5.08 ± 0.14 |
| 2–4 | 2.1 | 0.28 | 8.6 | n.d. | 5.61 ± 0.49 | |
| 4–6 | 2.0 | 0.27 | 8.7 | n.d. | 4.43 ± 0.43 | |
| 6–8 | 1.3 | 0.19 | 8.0 | n.d. | 3.12 ± 0.41 | |
| 08 | 0–2 | 1.7 | 0.24 | 8.3 | 111.6 | 6.48 ± 0.41 |
| 2–4 | 1.6 | 0.23 | 7.8 | n.d. | 5.13 ± 0.21 | |
| 4–6 | 1.4 | 0.21 | 7.9 | n.d. | 5.62 ± 0.40 | |
| 6–8 | 1.2 | 0.19 | 7.4 | n.d. | 2.70 ± 0.08 |
Organic carbon and nitrogen, C:N ratios, N-loss rates, and gene copy numbers of the investigated sediment layers.
n.s., not significant; n.d., not detectable; sed, sediment.
Benthic N-loss rates
Benthic N-loss activity was detected in both sediment slurries and intact sediment cores. In the intact core incubations, total benthic N-loss rates increased within the OMZ waters from 0.39 mmol N m−2 day−1 at 360 m to a maximum of 0.52 mmol N m−2 day−1 at 680 m (Figure 5A). At the lower boundary of the OMZ, rates decreased to 0.22 mmol N m−2 day−1 (970 m) and were the lowest at 1430 m (0.18 mmol N m−2 day−1). The relative contribution of denitrification and anammox to the total N-loss was estimated from slurry incubations. Denitrification rates in intact sediment cores ranged between 0.11 and 0.46 mmol N m−2 day−1, while anammox rates increased from 0.03 mmol N m−2 day−1 at the shallowest station to 0.07 mmol N m−2 day−1 at the deepest station (Figure 5A).
There were strong indications of the release of intracellular during the slurry incubations. The release of stored was most apparent in the measurements in the HgCl2-fixed subsamples from the initial time point (T0). -concentrations at T0 were significantly above the total sum of in the bottom water, pore water, and 15N-amendment combined, thus indicating an excess of in the slurry. Unfortunately, the true labeling percentage during the slurry incubation could not be determined from these subsamples, since any residual intracellular nitrate would have been released after poisoning with HgCl2. For this reason, was calculated from Eq. 3 (see Materials and Methods) and subsequently the excess concentrations of were calculated according to Eq. 4. Excess nitrate was calculated for all depths with denitrification rates (Table 1) and generally decreased with sediment depth. Excess nitrate ranged between 222 nmol N (cm3 sediment)−1 in the surface at station 16 and 3.1 nmol N (cm3 sediment)−1 in 4–6 cm at station 04.
In slurry incubations, both denitrification and anammox rates generally decreased with increasing sediment depth (Figures 4A,B). Due to insignificant 29N2 and 30N2 production, denitrification rates could not be obtained for 6–8 cm at all stations and 4–6 cm at stations 16 and 08. Denitrification rates at the sediment surface (0–2 cm layer) decreased with increasing water depth, from 136 nmol N cm−3 day−1 at 360 m to 73 nmol N cm−3 day−1 at 1430 m (Figure 4A). Anammox rates in surface sediments were lower than denitrification rates. However, in contrast to denitrification rates, anammox rates increased with water depth from 10 nmol N cm−3 day−1 at 360 m to 45 nmol N cm−3 day−1 at 1430 m (Figure 4B). Anammox and denitrification rates from slurry incubations were integrated down to the nitrate penetration depth of 2 cm (Figure 5B), which represents a rather conservative estimate, given that nitrate was found deeper in the sediment at some stations. Integrated denitrification rates decreased from 2.7 (±0.07) mmol m−2 day−1 at 360 m to 1.5 (±0.17) mmol m−2 day−1 at 1430 m. Anammox rates on the other hand increased with water depth from 0.21 (±0.03) mmol m−2 day−1 at 360 m to 0.89 (±0.04) mmol m−2 day−1 at the deepest station. As a result, the relative contribution of anammox to total N-loss increased with water depth from 7% at the shallowest station to 38% at the deepest station (Figure 5C).
Figure 4

Denitrification (A) and anammox (B) rates determined from slurry incubation experiments, along with the corresponding gene copy numbers of denitrifier nirS (A) and anammox nirS (B) quantified by real-time qPCR in the sediments. The asterisks mark the incubations that showed insignificant 15N–N2 production. Abbreviations: “n.d.” refers to non-detectable anammox nirS gene copy numbers.
Figure 5

(A) Anammox and denitrification rates measured in intact core incubation experiments, using the rIPT method described by Risgaard-Petersen et al. (
Detection of nirS genes from denitrifiers and anammox bacteria
The presence of microorganisms mediating the denitrification and anammox processes was verified by the detection of their respective biomarker functional genes nirS. Altogether, 225 denitrifier nirS sequences were obtained, and they formed 114 OTUs that could be grouped into seven clusters (Figure 6; Table A1 in Appendix). The nirS sequences from the Pakistan continental margin are diverse, and show clustering pattern that seems to be depth-related: certain clusters are dominated by sequences from the two shallow stations (10 and 04), while others are dominated by sequences from the two deeper stations (16 and 08). The majority of the sequences derived from stations 10 and 04 are found in clusters D2 (33 sequences) and D3 (49 sequences), to which the contributions from the deep stations (08 and 16) are considerably lower (only 13 sequences for D2 and 6 for D3). Meanwhile, clusters D1, D4, D5, and D6 seemed to be dominated by sequences from the deeper stations (08 and 16), with 41, 29, 6, and 9 sequences, respectively. Although the cd3aF–R3cd primer pair amplified predominantly denitrifier nirS genes, two sequences (OTU 04nir375) obtained from station 04 (680 m) were found to be more closely affiliated with the freshwater anammox bacterium Candidatus “Kuenenia stuttgartiensis” in cluster D7, with a similarity of 73% based on amino acid sequence. It should be noted that cluster D7 also includes cultured species like Halomonas campisalis and Methylomirabilis oxyfera, which share up to 59 and 69% amino acid sequence similarity, respectively, to the Arabian Sea D7 sequences obtained in this study.
Figure 6

Phylogenetic reconstruction of the denitrifier NirS based on amino acids sequences translated from gene sequences. Sequences were retrieved from clone libraries constructed for sediments of all stations based on maximum likelihood and maximum parsimony algorithms. Bootstrapped values of >50% are shown for maximum likelihood/maximum parsimony. Indicated in black are the related reference sequences obtained from GenBank. Labeling of sequences: station = 04, 08, 10, or 16, “nir” = amplicons from primers cd3aF/R3cd, or “sc” = amplicons from primers Scnir372F/Scnir845R, followed by unique sequence number. Numbers in parentheses are the numbers of sequences represented by the same OTU with ≥95% nucleic acids sequence similarity. OTUs from the shallow stations are in orange red, while OTUs from deep stations are in purple. D1–D7 indicate the different clusters identified in this study, while the ratio below gives the ratio of sequences from shallow stations (10 and 04) to the deeper stations (08 and 16).
A total of 109 OTUs from 241 anammox nirS sequences were retrieved from the Pakistani margin sediments (Figure 7), and they formed three clusters that might also carry some water-depth-related pattern, though not as obvious as for the denitrifier nirS sequences. Cluster S1 and S3, closely related to OTUs from the Arabian Sea water column, were dominated by sequences from deep stations (16 and 08) with 33 sequences compared to 17 and 22 sequences from the shallow stations (10 and 04). In contrast, cluster S2 affiliated with OTUs from the Peruvian water column seemed to have similar contributions from both shallow and deep stations.
Figure 7

Phylogenetic reconstruction of the anammox (Scalindua) NirS based on amino acids sequences. In blue are the OTUs from the deep stations (08 and 16), while the green sequences were obtained from the shallow stations (10 and 04). Please refer to Figure 6 for additional information regarding the sequence labels in tree.
Quantification of denitrifier- and anammox-nirS genes
Consistent with benthic N-loss rate measurements, the anammox nirS genes were generally less abundant than denitrifier nirS genes (Figures 4A,B). Both nirS gene copy numbers showed a decreasing trend with sediment depth. Amongst all stations, the highest denitrifier nirS gene abundance of 7245 ± 813 gene copies (mg dry sediment)−1 was detected in the surface sediment layer at station 16 (970 m), whereas the lowest denitrifier nirS abundance of 439 (±90) gene copies (mg dry sediment)−1 was detected in the uppermost 2 cm at the deepest station 08.
The abundance of anammox nirS genes was usually an order of magnitude lower than that of the denitrifier nirS (Figure 4B), and was often found to be close to the detection limit. Similar to the denitrifier nirS genes, the highest numbers of anammox nirS genes were also detected at station 16, ranging from 1728 ± 198 gene copies (mg dry sediment)−1 in the surface to undetectable at 6–8 cm. Although the highest rates of anammox were measured in the slurry incubation experiment at station 08, only low gene copy numbers of anammox nirS, in the range of 93 ± 44 to 203 ± 44 gene copies (mg dry sediment)−1, were detected.
The relative contribution of the anammox nirS to the total nirS gene copy numbers in the uppermost 2 cm increased with water depth from 1% at 360 m to 19% at station 16 (Figure 5C). These results are consistent with depth-integrated rates, which show an increase of anammox contribution to total N-loss with increasing water depth.
Discussion
Benthic N-loss due to denitrification
Consistent with previous benthic N-loss studies from other continental slopes, e.g., the North Atlantic (Trimmer and Nicholls,
Denitrification rates have been determined for the continental shelf sediments off central Chile, where seasonal hypoxia develops each year (Farías et al.,
Several lines of observations collectively indicate the likely presence of nitrate-storing organisms in the sediments of the Pakistani margin. Firstly, high production of 29N2 relative to 30N2 was measured in the incubations, which did not agree with the calculated labeling percentage and the measured anammox rates. Secondly, nitrate concentrations in the T0 subsamples of the slurry incubations exceeded the sum of bottom water, pore water and 15N-nitrate. Thirdly, subsurface maxima of pore water nitrate/nitrite, similar to those previously observed at the Pakistani margin (Woulds et al.,
Nitrate-storing sulfur bacteria, such as Thioploca spp. and Beggiotoa spp., have been associated with high pore water nitrate concentrations (Fossing et al.,
Denitrification from the stored by foraminifera would lead to false denitrification estimates if the intracellular labeling percentage was not known. However, the increased concentrations in the slurry subsamples at T0 suggest that the stored was released into the pore water when the slurry was mixed at the start of the experiment. Thus, a subsequent uptake of from the pore water would lead to an intracellular that is close to the pore water . Furthermore, the linear increase of 29N2 and 30N2 with time indicates that either intracellular did not change over time or that the N2 production by foraminifera was minor, as was also observed in other regions (Glud et al.,
The dominance of denitrification in benthic N-loss in the Pakistan margin sediments is strongly supported by the high abundance of denitrifier nirS genes. Moreover, the gene copy numbers generally followed similar decreasing trends as the rates measured in slurry incubation within the sediments (Figure 4). Exceptions were noted particularly in the topmost layer(s) at the deepest station (station 8), and these could potentially be due to nucleic acid extraction efficiency or biases, and/or the presence of PCR inhibitors. In addition, the primers used only target nirS, while any occurrence of the nirK genes would not have been accounted for. Although there are also primers designed for nirK, those currently available may also target those of nitrifiers. Consequently, quantification of nirK in addition to that of nirS would likely overestimate denitrifier abundance instead. Future refinement of primer designs, or the assessment of multiple biomarker genes in parallel, may help shed light on the true quantitative distribution of denitrifiers in the environment. Compared to previous studies in various sediments, most of which also focused on denitrifier nirS and found gene copy numbers ranging from ca. 0.6 × 103 copies (mg sediment)−1 at the mouth of the Colne estuary (Smith et al.,
In agreement with studies addressing nirS genes in the water column of the Arabian Sea (Jayakumar et al.,
It is particularly interesting to find an OTU amplified with the primers targeting denitrifier nirS genes, to be related to the Ca. “K. stuttgartiensis” (73% similarity, Figure 6). Ca. “K. stuttgartiensis” is an anammox bacterium known to occur in freshwater (Jetten et al.,
Benthic N-loss via anammox
This study provides the first direct measurement of anammox activity in the sediments of the Arabian Sea, or any OMZs. The very recent study by Bohlen et al. (
Because the nirS gene is present as a single copy in anammox bacteria, according to the sequenced genomes of both the freshwater Ca. “K. stuttgartiensis” (Strous et al.,
According to the phylogenetic reconstruction of Scalindua NirS (Figure 6), three different clusters could be identified and the diversity of the community (Chao1 richness estimate = 275), though lower than the diversity of the denitrifier NirS (Figures A1A,C in Appendix), seems to be higher compared to those found in the water column OMZ of the Arabian Sea (Chao1 richness estimate = 8; Jensen et al.,
Anammox contribution increased with water depth
In agreement with other studies (Engström et al.,
Although anammox rates and cell abundance increase with water depth, it is unlikely that water depth or rather pressure itself is a direct regulating factor for the anammox contribution, since bacterial communities and denitrifiers in particular are able to cope with high pressure very well (Tamegai et al.,
Organic carbon concentrations usually decrease with water depth and therefore it is hypothesized in some studies (Nicholls and Trimmer,
Anammox activity depends on sufficient supplies of (Dalsgaard and Thamdrup,
Contribution of benthic N-loss to the N-deficit in the Arabian Sea
In general, the water column of the central Arabian Sea is believed to be an important sink for fixed nitrogen in global oceans as indicated by a prominent N-deficit (Naqvi,
Naqvi et al. (
Furthermore, N-loss rates measured in the central Arabian Sea of 0.3–0.6 mmol N m−2 day−1 (Jensen et al.,
Conclusion
Benthic N-loss due to anammox increased with water depth on the Pakistan margin and the contribution of anammox to total N-loss seemed to co-vary with temperature and nitrate concentrations in the bottom water. Compared to shallow sediments, anammox bacteria seem to be more successful in deeper sediments, as anammox accounted for almost 40% to the total benthic N-loss at 1430 m water depth. The shift from a denitrifier-dominated heterotrophic system in shallow sediments, to a system in which the autotrophic anammox process plays a more important role in sediments at deeper water depths, could also be coupled to the availability of labile organic carbon. Owing to their chemolithoautotrophic lifestyle, anammox bacteria could have a competitive advantage over denitrifiers in deeper sediments due to their lesser dependence on the often seasonally fluctuating primary production in surface waters for sources of electron donor and carbon. Extrapolation from our data suggests that benthic N-loss could account for up to half of the total N-loss in the Arabian Sea as a whole, and may thus have contributed to the N-deficits in the water column, though further investigations during different seasons are necessary to fully evaluate the role of sediments in the annual marine N-loss. Since human populations and anthropogenic atmospheric N deposition (Duce et al.,
Statements
Acknowledgments
We are very thankful for the technical support and analyses by G. Klockgether, D. Franzke, M. Meier, I. Vieweg, A. Schipper, S. Kühn, and for the excellent cooperation with captain and crew of R/V Meteor M74/2. We thank K. Zonneveld and S. Kasten for providing the equipment and fruitful discussions. Funding from DFG-Research Center/Excellence Cluster “The Ocean in the Earth System” (MARUM) and the Max Planck Society are gratefully acknowledged.
Conflict of interest
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.
References
1
AbellG. C. J.RevillA. T.SmithC.BissettA. P.VolkmanJ. K.RobertS. S. (2010). Archaeal ammonia oxidizers and nirS-type denitrifiers dominate sediment nitrifying and denitrifying populations in a subtropical macrotidal estuary. ISME J.4, 286–300.10.1038/ismej.2009.105
2
AltschulS. F.MaddenT. L.SchafferA. A.ZhangJ. H.ZhangZ.MillerW.et al (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res.25, 3389–3402.10.1093/nar/25.17.3389
3
BangeH. W.RixenT.JohansenA. M.SiefertR. L.RameshR.IttekkotV.et al (2000). A revised nitrogen budget for the Arabian Sea. Global Biogeochem. Cycles14, 1283–1297.10.1029/1999GB001228
4
BohlenL.DaleA. W.SommerS.MoschT.HensenC.NoffkeA.et al (2011). Benthic nitrogen cycling traversing the Peruvian oxygen minimum zone. Geochim. Cosmochim. Acta75, 6094–6111.10.1016/j.gca.2011.08.010
5
BramanR. S.HendrixS. A. (2002). Nanogram nitrite and nitrate determination in environmental and biological materials by vanadium (III) reduction with chemiluminescence detection. Anal. Chem.61, 2715–2718.10.1021/ac00199a007
6
BulowS. E.FrancisC. A.JacksonG. A.WardB. B. (2008). Sediment denitrifier community composition and nirS gene expression investigated with functional gene microarrays. Environ. Microbiol.10, 3057–3069.10.1111/j.1462-2920.2008.01765.x
7
ByrneN.StrousM.CrepeauV.KartalB.BirrienJ.-L.SchmidM.et al (2009). Presence and activity of anaerobic ammonium-oxidizing bacteria at deep-sea hydrothermal vents. ISME J.3, 117–123.10.1038/ismej.2008.72
8
Castro-GonzalezM.BrakerG.FariasL.UlloaO. (2005). Communities of nirS-type denitrifiers in the water column of the oxygen minimum zone in the eastern South Pacific. Environ. Microbiol.7, 1298–1306.10.1111/j.1462-2920.2005.00809.x
9
ClineJ. D. (1969). Spectrophotometric determination of hydrogen sulfide in natural waters. Limnol. Oceanogr.14, 454–458.10.4319/lo.1969.14.3.0454
10
CodispotiL. A.BrandesJ. A.ChristensenJ. P.DevolA. H.NaqviS. W. A.PaerlH. W.et al (2001). The oceanic fixed nitrogen and nitrous oxide budgets: moving targets as we enter the anthropocene?Sci. Mar.65, 85–105.10.3989/scimar.2001.65s285
11
DalsgaardT.CanfieldD. E.PetersenJ.ThamdrupB.Acuna-GonzalezJ. (2003). Anammox is a significant pathway of N2 production in the anoxic water column of Golfo Dulce, Costa Rica. Nature422, 606–608.10.1038/nature01526
12
DalsgaardT.ThamdrupB. (2002). Factors controlling anaerobic ammonium oxidation with nitrite in marine sediments. Appl. Environ. Microbiol.68, 3802–3808.10.1128/AEM.68.8.3802-3808.2002
13
DalsgaardT.ThamdrupB.CanfieldD. E. (2005). Anaerobic ammonium oxidation (anammox) in the marine environment. Res. Microbiol.156, 457–464.10.1016/j.resmic.2005.01.011
14
DangH.ChenR.WangL.GuoL.ChenP.TangZ.et al (2010). Environmental factors shape sediment anammox bacterial communities in hypernutrified Jiaozhou Bay, China. Appl. Environ. Microbiol.76, 7036–7047.10.1128/AEM.01264-10
15
DangH.WangC.LiJ.LiT.TianF.JinW.et al (2009). Diversity and distribution of sediment NirS-encoding bacterial assemblages in response to environmental gradients in the eutrophied Jiaozhou Bay, China. Microb. Ecol.58, 161–169.10.1007/s00248-008-9469-5
16
DeutschC.GruberN.KeyR.SarmientoJ. L.GanachaudA. (2001). Denitrification and N2 fixation in the Pacific Ocean. Global Biogeochem. Cycles15, 483–506.10.1029/2000GB001291
17
DevolA. H.NaqviS. W. A.CodispotiL. A. (2006). “Nitrogen cycling in the suboxic waters of the Arabian sea,” in NATO Science Series IV Earth and Environmental Sciences: 64, Past and Present Water Column Anoxia, ed. NeretinL. (Amsterdam: IOS Press and Kluwer Academic Publishers in conjunction with the NATO Scientific Affairs Division), 283–310.
18
DeVriesT.DeutschC.PrimeauF.ChangB.DevolA. (2012). Global rates of water-column denitrification derived from nitrogen gas measurements. Nat. Geosci.5, 547–550.10.1038/ngeo1515
19
DuceR. A.LarocheJ.AltieriK.ArrigoK. R.BakerA. R.CaponeD. G.et al (2008). Impacts of atmospheric anthropogenic nitrogen on the open ocean. Science320, 893–897.10.1126/science.1150369
20
EngströmP.DalsgaardT.HulthS.AllerR. C. (2005). Anaerobic ammonium oxidation by nitrite (anammox): implications for N2 production in coastal marine sediments. Geochim. Cosmochim. Acta69, 2057–2065.10.1016/j.gca.2004.09.032
21
EngströmP.PentonC. R.DevolA. H. (2009). Anaerobic ammonium oxidation in deep-sea sediments off the Washington margin. Limnol. Oceanogr.54, 1643–1652.10.4319/lo.2009.54.5.1643
22
EttwigK. F.ButlerM. K.Le PaslierD.PelletierE.MangenotS.KuypersM. M. M.et al (2010). Nitrite-driven anaerobic methane oxidation by oxygenic bacteria. Nature464, 543–548.10.1038/nature08883
23
FaríasL.GracoM.UlloaO. (2004). Temporal variability of nitrogen cycling in continental-shelf sediments of the upwelling ecosystem off central Chile. Deep Sea Res. Part II Top. Stud. Oceanogr.51, 2491–2505.10.1016/j.dsr2.2004.07.029
24
FossingH.GallardoV. A.JorgensenB. B.HuttelM.NielsenL. P.SchulzH.et al (1995). Concentration and transport of nitrate by the mat-forming sulfur bacterium Thioploca. Nature374, 713–715.10.1038/374713a0
25
GallowayJ. N.DentenerF. J.CaponeD. G.BoyerE. W.HowarthR. W.SeitzingerS. P.et al (2004). Nitrogen cycles: past, present, and future. Biogeochemistry70, 153–226.10.1007/s10533-004-0370-0
26
GaoH.SchreiberF.CollinsG.JensenM. M.KostkaJ. E.LavikG.et al (2009). Aerobic denitrification in permeable Wadden sea sediments. ISME J.4, 417–426.10.1038/ismej.2009.127
27
GilloolyJ. F.BrownJ. H.WestG. B.SavageV. M.CharnovE. L. (2001). Effects of size and temperature on metabolic rate. Science293, 2248–2251.10.1126/science.1061967
28
GilloolyJ. F.CharnovE. L.WestG. B.SavageV. M.BrownJ. H. (2002). Effects of size and temperature on developmental time. Nature417, 70–73.10.1038/417070a
29
GinestetP.AudicJ. M.UrbainV.BlockJ. C. (1998). Estimation of nitrifying bacterial activities by measuring oxygen uptake in the presence of the metabolic inhibitors allylthiourea and azide. Appl. Environ. Microbiol.64, 2266–2268.
30
GludR. N.ThamdrupB.StahlH.WenzhoeferF.GludA.NomakiH.et al (2009). Nitrogen cycling in a deep ocean margin sediment (Sagami Bay, Japan). Limnol. Oceanogr.54, 723–734.10.4319/lo.2009.54.3.0723
31
GrasshoffK.JohannsenH. (1972). New sensitive and direct method for automatic determination of ammonia in sea-water. ICES J. Mar. Sci.34, 516–521.10.1093/icesjms/34.3.516
32
GrasshoffK.KremlingK.EhrhardtM. (1999). Methods of Seawater Analysis. Weinheim: Wiley-VCH Verlag GmbH.
33
GruberN. (2004). “The dynamics of the marine nitrogen cycle and its influence on atmospheric CO2 variations,” in NATO Science Series IV Earth and Environmental Sciences: 40, The Ocean Carbon Cycle and Climate, eds FollowsM.OguzT. (Dordrecht: Kluwer Academic Publishers), 97–148.
34
GruberN.SarmientoJ. L. (1997). Global patterns of marine nitrogen fixation and denitrification. Global Biogeochem. Cycles11, 235–266.10.1029/97GB00077
35
HallT. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser.41, 95–98.
36
HarhangiH. R.Le RoyM.Van AlenT.HuB.-L.GroenJ.KartalB.et al (2012). Hydrazine synthase, a unique phylomarker with which to study the presence and biodiversity of anammox bacteria. Appl. Environ. Microbiol.78, 752–758.10.1128/AEM.07113-11
37
HiraD.TohH.MigitaC. T.OkuboH.NishiyamaT.HattoriM.et al (2012). Anammox organism KSU-1 expresses a NirK-type copper-containing nitrite reductase instead of a NirS-type with cytochrome cd(1). FEBS Lett.586, 1658–1663.10.1016/j.febslet.2012.04.041
38
HolmesR. M.AminotA.KerouelR.HookerB. A.PetersonB. J. (1999). A simple and precise method for measuring ammonium in marine and freshwater ecosystems. Can. J. Fish. Aquat. Sci.56, 1801–1808.10.1139/f99-128
39
HoltappelsM.LavikG.JensenM. M.KuypersM. M. M. (2011). “(15)N-labeling experiments to dissect the contributions of heterotrophic denitrification and anammox to nitrogen removal in the OMZ waters of the ocean,” in Methods in Enzymology: Research on Nitrification and Related Processes, Vol. 486, Part A, ed. KlotzM. (San Diego: Elsevier Academic Press Inc.), 223–251.
40
JayakumarD. A.FrancisC. A.NaqviS. W. A.WardB. B. (2004). Diversity of nitrite reductase genes (nirS) in the denitrifying water column of the coastal Arabian Sea. Aquat. Microb. Ecol.34, 69–78.10.3354/ame034069
41
JensenM. M.LamP.RevsbechN. P.NagelB.GayeB.JettenM. S. M.et al (2011). Intensive nitrogen loss over the Omani Shelf due to anammox coupled with dissimilatory nitrite reduction to ammonium. ISME J.5, 1660–1670.10.1038/ismej.2011.44
42
JettenM. S. M.SliekersO.KuypersM.DalsgaardT.van NiftrikL.CirpusI.et al (2003). Anaerobic ammonium oxidation by marine and freshwater planctomycete-like bacteria. Appl. Microbiol. Biotechnol.63, 107–114.10.1007/s00253-003-1422-4
43
KartalB.KolevaM.ArsovR.van der StarW.JettenM. S. M.StrousM. (2006). Adaptation of a freshwater anammox population to high salinity wastewater. J. Biotechnol.126, 546–553.10.1016/j.jbiotec.2006.05.012
44
KartalB.MaalckeW. J.de AlmeidaN. M.CirpusI.GloerichJ.GeertsW.et al (2011). Molecular mechanism of anaerobic ammonium oxidation. Nature479, 127–130.10.1038/nature10453
45
KuypersM. M. M.LavikG.WoebkenD.SchmidM.FuchsB. M.AmannR.et al (2005). Massive nitrogen loss from the Benguela upwelling system through anaerobic ammonium oxidation. Proc. Natl. Acad. Sci. U.S.A.102, 6478–6483.10.1073/pnas.0502088102
46
KuypersM. M. M.SliekersA. O.LavikG.SchmidM.JorgensenB. B.KuenenJ. G.et al (2003). Anaerobic ammonium oxidation by anammox bacteria in the Black Sea. Nature422, 608–611.10.1038/nature01472
47
LamP.JensenM. M.KockA.LettmannK. A.PlancherelY.LavikG.et al (2011). Origin and fate of the secondary nitrite maximum in the Arabian sea. Biogeosciences8, 1565–1577.10.5194/bg-8-1565-2011
48
LamP.LavikG.JensenM. M.van de VossenbergJ.SchmidM.WoebkenD.et al (2009). Revising the nitrogen cycle in the Peruvian oxygen minimum zone. Proc. Natl. Acad. Sci. U.S.A.106, 4752–4757.10.1073/pnas.0812444106
49
LiM.FordT.LiX.GuJ.-D. (2011). Cytochrome cd1-containing nitrite reductase encoding gene nirS as a new functional biomarker for detection of anaerobic ammonium oxidizing (Anammox) bacteria. Environ. Sci. Technol.45, 3547–3553.10.1021/es103830z
50
LiuX. D.TiquiaS. M.HolguinG.WuL. Y.NoldS. C.DevolA. H.et al (2003). Molecular diversity of denitrifying genes in continental margin sediments within the oxygen-deficient zone off the Pacific coast of Mexico. Appl. Environ. Microbiol.69, 3549–3560.10.1128/AEM.69.4.2110-2115.2003
51
LudwigW.StrunkO.WestramR.RichterL.MeierH.Yadhukumaret al (2004). ARB: a software environment for sequence data. Nucleic Acids Res.32, 1363–1371.10.1093/nar/gkh293
52
MarraJ.BarberR. T. (2005). Primary productivity in the Arabian sea: a synthesis of JGOFS data. Prog. Oceanogr.65, 159–175.10.1016/j.pocean.2005.03.004
53
MichoteyV.MejeanV.BoninP. (2000). Comparison of methods for quantification of cytochrome cd1-denitrifying bacteria in environmental marine samples. Appl. Environ. Microbiol.66, 1564–1571.10.1128/AEM.66.4.1564-1571.2000
54
MormileM. R.RomineM. F.GarciaT.VentosaA.BaileyT. J.PeytonB. M. (1999). Halomonas campisalis sp nov., a denitrifying, moderately haloalkaliphilic bacterium. Syst. Appl. Microbiol.22, 551–558.10.1016/S0723-2020(99)80008-3
55
NaqviS. W. A. (1994). Denitrification processes in the Arabian sea. Proceedings of the Indian Academy of Sciences – Earth and Planetary Sciences103, 279–300.
56
NaqviS. W. A.NaikH.PratiharyA.D’SouzaW.NarvekarP. V.JayakumarD. A.et al (2006). Coastal versus open-ocean denitrification in the Arabian sea. Biogeosciences3, 621–633.10.5194/bg-3-621-2006
57
NichollsJ. C.TrimmerM. (2009). Widespread occurrence of the anammox reaction in estuarine sediments. Aquat. Microb. Ecol.55, 105–113.10.3354/ame01285
58
Piña-OchoaE.HøgslundS.GeslinE.CedhagenT.RevsbechN. P.NielsenL. P.et al (2010). Widespread occurrence of nitrate storage and denitrification among Foraminifera and Gromiida. Proc. Natl. Acad. Sci. U.S.A.107, 1148–1153.10.1073/pnas.0908440107
59
Risgaard-PetersenN.LangezaalA. M.IngvardsenS.SchmidM. C.JettenM. S. M.Op den CampH. J. M.et al (2006). Evidence for complete denitrification in a benthic foraminifer. Nature443, 93–96.10.1038/nature05070
60
Risgaard-PetersenN.NielsenL. P.RysgaardS.DalsgaardT.MeyerR. L. (2003). Application of the isotope pairing technique in sediments where anammox and denitrification coexist. Limnol. Oceanogr. Methods1, 63–73.10.4319/lom.2011.1.63
61
RysgaardS.GludR. N.Risgaard-PetersenN.DalsgaardT. (2004). Denitrification and anammox activity in Arctic marine sediments. Limnol. Oceanogr.49, 1493–1502.10.4319/lo.2004.49.1.0086
62
SangerF.NicklenS.CoulsonA. R. (1977). DNA sequencing with chain-terminating inhibitors. Proc. Natl. Acad. Sci. U.S.A.74, 5463–5467.10.1073/pnas.74.12.5463
63
SavageV. M.GilloolyJ. F.BrownJ. H.WestG. B.CharnovE. L. (2004). Effects of body size and temperature on population growth. Am. Nat.163, 429–441.10.1086/381872
64
SchalkJ.de VriesS.KuenenJ. G.JettenM. S. M. (2000). Involvement of a novel hydroxylamine oxidoreductase in anaerobic ammonium oxidation. Biochemistry39, 5405–5412.10.1021/bi992721k
65
SchlossP. D.WestcottS. L.RyabinT.HallJ. R.HartmannM.HollisterE. B.et al (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol.75, 7537–7541.10.1128/AEM.01541-09
66
SchlüterM. (1990). “Zur Frühdiagenese von organischem Kohlenstoff und Opal in Sedimenten des südlichen und östlichen Weddellmeeres: geochemische Analyse und Modellierung=Early diagenesis of organic carbon and opal in sediments of the southern and eastern Weddell Sea: geochemical analysis and modelling,” in Berichte zur Polarforschung (Reports on Polar Research), Bremerhaven: Alfred Wegener Institute for Polar and Marine Research.
67
SchottW.von StackelbergU.EckhardtF. J.MattiatB.PetersJ.ZobelB. (1970). Geologische Untersuchungen an Sedimenten des indisch-pakistanischen Kontinentalrandes (Arabisches Meer). Int. J. Earth Sci.60, 264–275.
68
SchumacherS.JorissenF. J.DissardD.LarkinK. E.GoodayA. J. (2007). Live (Rose Bengal stained) and dead benthic foraminifera from the oxygen minimum zone of the Pakistan continental margin (Arabian Sea). Mar. Micropaleontol.62, 45–73.10.1016/j.marmicro.2006.07.004
69
SchwartzM. C.WouldsC.CowieG. L. (2009). Sedimentary denitrification rates across the Arabian Sea oxygen minimum zone. Deep Sea Res. Part II Top. Stud. Oceanogr.56, 324–332.10.1016/j.dsr2.2008.05.028
70
SeitzingerS.GiblinA. (1996). Estimating denitrification in North Atlantic continental shelf sediments. Biogeochemistry35, 235–260.10.1007/BF02179829
71
ShetyeS.GouveiaA.ShenoiS. (1994). Circulation and water masses of the Arabian Sea. J. Earth Syst. Sci.103, 107–123.
72
SmithC. J.NedwellD. B.DongL. F.OsbornA. M. (2007). Diversity and abundance of nitrate reductase genes (narG and napA), nitrite reductase genes (nirS and nrfA), and their transcripts in estuarine sediments. Appl. Environ. Microbiol.73, 3612–3622.10.1128/AEM.01278-07
73
StrousM.PelletierE.MangenotS.RatteiT.LehnerA.TaylorM. W.et al (2006). Deciphering the evolution and metabolism of an anammox bacterium from a community genome. Nature440, 790–794.10.1038/nature04647
74
TamegaiH.LiL.MasuiN.KatoC. (1997). A denitrifying bacterium from the deep sea at 11,000-m depth. Extremophiles1, 207–211.10.1007/s007920050035
75
ThamdrupB.DalsgaardT. (2002). Production of N2 through anaerobic ammonium oxidation coupled to nitrate reduction in marine sediment. Appl. Environ. Microbiol.68, 1312–1318.10.1128/AEM.68.3.1312-1318.2002
76
ThompsonJ. D.HigginsD. G.GibsonT. J. (1994). Clustal-W – improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res.22, 4673–4680.10.1093/nar/22.22.4673
77
ThrobackI. N.EnwallK.JarvisA.HallinS. (2004). Reassessing PCR primers targeting nirS, nirK and nosZ genes for community surveys of denitrifying bacteria with DGGE. FEMS Microbiol. Ecol.49, 401–417.10.1016/j.femsec.2004.04.011
78
TiquiaS. M.MassonS. A.DevolA. (2006). Vertical distribution of nitrite reductase genes (nirS) in continental margin sediments of the Gulf of Mexico. FEMS Microbiol. Ecol.58, 464–475.10.1111/j.1574-6941.2006.00173.x
79
TrimmerM.EngströmP. (2011). “Distribution, activity, and ecology of anammox bacteria in aquatic environments,” in Nitrification, eds WardB.ArpD.KlotzM. (Washingtopn: American Society for Microbiology), 201–235.
80
TrimmerM.NichollsJ. C. (2009). Production of nitrogen gas via anammox and denitrification in intact sediment cores along a continental shelf to slope transect in the North Atlantic. Limnol. Oceanogr.54, 577–589.10.4319/lo.2009.54.2.0577
81
van de VossenbergJ.RattrayJ. E.GeertsW.KartalB.van NiftrikL.van DonselaarE. G.et al (2008). Enrichment and characterization of marine anammox bacteria associated with global nitrogen gas production. Environ. Microbiol.10, 3120–3129.10.1111/j.1462-2920.2008.01643.x
82
WardB. B.DevolA. H.RichJ. J.ChangB. X.BulowS. E.NaikH.et al (2009). Denitrification as the dominant nitrogen loss process in the Arabian Sea. Nature461, 78–81.10.1038/nature08276
83
WiggertJ. D.HoodR. R.BanseK.KindleJ. C. (2005). Monsoon-driven biogeochemical processes in the Arabian Sea. Prog. Oceanogr.65, 176–213.10.1016/j.pocean.2005.03.008
84
WouldsC.SchwartzM. C.BrandT.CowieG. L.LawG.MowbrayS. R. (2009). Pore water nutrient concentrations and benthic nutrient fluxes across the Pakistan margin OMZ. Deep Sea Res. Part II Top. Stud. Oceanogr.56, 333–346.10.1016/j.dsr2.2008.10.008
Appendix
Table A1
| Denitrifier or Kuenenia-like | Scalindua-like | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| D1 | D2 | D3 | D4 | D5 | D6 | D7 | Kuenenia stuttgartiensis | Methylomirabilis oxyfera | S1 | S2 | S3 | Scalindua profunda | |
| D1 | |||||||||||||
| D2 | 82–54 | ||||||||||||
| D3 | 82–60 | 76–56 | |||||||||||
| D4 | 70–47 | 71–53 | 70–54 | ||||||||||
| D5 | 74–51 | 70–55 | 70–56 | 72–55 | |||||||||
| D6 | 59–38 | 60–47 | 58–44 | 62–47 | 63–48 | ||||||||
| D7 | 58–38 | 58–43 | 59–41 | 58–46 | 61–45 | 66–48 | |||||||
| Kuenenia stuttgartiensis | 61–48 | 60–53 | 61–54 | 62–55 | 62–53 | 60–51 | 73–38 | ||||||
| Methylomirabilis oxyfera | 58–42 | 58–49 | 55–45 | 59–52 | 60–54 | 62–47 | 69–42 | 56 | |||||
| S1 | 62–39 | 61–42 | 57–39 | 63–44 | 66–47 | 66–46 | 60–44 | 61–51 | 64–56 | ||||
| S2 | 61–41 | 62–44 | 57–40 | 65–43 | 68–48 | 65–49 | 63–44 | 61–50 | 66–54 | 94–57 | |||
| S3 | 66–38 | 62–37 | 59–36 | 65–38 | 67–44 | 65–41 | 63–42 | 59–44 | 64–49 | 85–52 | 84–50 | ||
| Scalindua profunda | 60–49 | 61–52 | 57–47 | 64–54 | 65–57 | 64–56 | 62–45 | 51 | 60 | 94–77 | 95–77 | 82–64 | |
Range of similarities (%) between the different clusters identified in the NirS phylogenetic trees, based on amino acid sequences.
Figure A1

Rarefaction curves for denitrifier and anammox (Scalindua) nirS sequences based on 95% and 90% nucleotide sequences cut offs. (A) Shows rarefaction curves over all sampled stations, (B) shows the denitrifier nirS rarefaction curves separately for each sampled station, and (C) the corresponding curves for anammox (Scalindua) nirS.
Summary
Keywords
Arabian Sea, benthic N-loss, anammox, denitrification, qPCR, nirS
Citation
Sokoll S, Holtappels M, Lam P, Collins G, Schlüter M, Lavik G and Kuypers MMM (2012) Benthic Nitrogen Loss in the Arabian Sea Off Pakistan. Front. Microbio. 3:395. doi: 10.3389/fmicb.2012.00395
Received
09 July 2012
Accepted
29 October 2012
Published
28 November 2012
Volume
3 - 2012
Edited by
Hongyue Dang, China University of Petroleum, China
Reviewed by
Susannah Green Tringe, DOE Joint Genome Institute, USA; Zhe-Xue Quan, Fudan University, China
Copyright
© 2012 Sokoll, Holtappels, Lam, Collins, Schlüter, Lavik and Kuypers.
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Sarah Sokoll, Biogeochemistry Department, Max Planck Institute for Marine Microbiology, Celsiusstrasse 1, D-28359 Bremen, Germany. e-mail: ssokoll@mpi-bremen.de
This article was submitted to Frontiers in Aquatic Microbiology, a specialty of Frontiers in Microbiology.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.