ORIGINAL RESEARCH article
Niche Differentiation of Aerobic and Anaerobic Ammonia Oxidizers in a High Latitude Deep Oxygen Minimum Zone
- 1Department of Limnology and Bio-Oceanography, Center of Functional Ecology, University of Vienna, Vienna, Austria
- 2NIOZ, Department of Marine Microbiology and Biogeochemistry, Royal Netherlands Institute for Sea Research, Utrecht University, Den Burg, Netherlands
- 3Research and Development Center for Marine Biosciences, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Yokosuka, Japan
- 4Ecosystem Oceanography Group (GRECO), Instituto Español de Oceanografía, Centro Oceanográfico de Baleares, Palma, Spain
To elucidate the potential for nitrification and denitrification processes in a high latitude deep oxygen minimum zone (OMZ) we determined the abundance and community composition of the main microbial players in the aerobic and anaerobic (anammox) ammonium oxidation and denitrification processes in the Gulf of Alaska throughout the water column. Within the dominant bacterial groups, Flavobacterales, Rhodobacterales, Actinomarinales, and SAR86 were more abundant in epipelagic waters and decreased with depth, whereas SAR11, SAR324, Marinimicrobia, and Thiomicrospirales increased their contribution to the bacterial community with depth. Nitrosopumilaceae also increased with depth and dominated the OMZ and bathypelagic archaeal communities. Euryarchaeota Marine Group II exhibited an opposite depth pattern to Nitrosopumilaceae, whereas Marine Group III and Woesearchaeota were more abundant in the bathypelagic realm. Candidatus Brocadia contributed 70–100% of the anammox bacterial community throughout the water column. Archaeal ammonia oxidizers (AOA) dominated the microbial community involved in the nitrogen cycle. Two AOA ecotypes, the high ammonia (HAC) and low ammonia (LAC)-AOA, characterized by distinct genes for aerobic ammonia oxidation (amoA) and for denitrification (nirK), exhibited a distinct distribution pattern related to depth and ammonia concentrations. HAC-AOA dominated in epipelagic (80.5 ± 28.3% of total AOA) oxygenated and ammonia-rich waters, and LAC-AOA dominated in the OMZ (90.9 ± 5.1%) and bathypelagic waters (85.5 ± 13.5%), characterized by lower oxygen and ammonia concentrations. Bacterial denitrifiers (3.7 ± 6.9 bacterial nirK gene mL−1) and anaerobic ammonia oxidizers (78 ± 322 anammox 16S rRNA genes L−1) were low in abundance under the oxygen conditions in the Gulf of Alaska throughout the water column. The widespread distribution of bacterial denitrifiers and anaerobic ammonia oxidizers in low abundances reveals a reservoir of genetic and metabolic potential ready to colonize the environment under the predicted increase of OMZs in the ocean. Taken together, our results reinforce the niche partitioning of archaeal ammonia oxidizers based on their distinct metabolic characteristics resulting in the dominance of LAC-AOA in a high latitude deep OMZ. Considering the different ecological roles and functions of the two archaeal ecotypes, the expansion of the zones dominated by the LAC-ecotype might have implications for the nitrogen cycle in the future ocean.
Microorganisms mediate most of the biogeochemical transformations in the global nitrogen (N) cycle (Kuypers et al., 2018). Over the last decade, fundamental pathways and key microbial players in the N cycle have been discovered (Kuypers et al., 2018). However, major uncertainties still exist on the extent and connection of nitrification and denitrification processes, especially in the open ocean (Ward and Jensen, 2014). The coupling of these processes, which affects the flow of N in the ecosystems, requires close interaction between nitrifying and denitrifying microorganisms, both spatially and/or temporally (Ward, 1996; Zehr and Ward, 2002). Oxygen minimum zones (OMZs) play an important role in the global ocean nitrogen cycle (Lam and Kuypers, 2011) providing an array of niches inhabited by metabolically diverse microorganisms (Bertagnolli and Stewart, 2018).
Nitrification is the two-step aerobic oxidation of ammonia (NH3) via nitrite () to nitrate (), mediated by ammonia-oxidizing Archaea and Bacteria and nitrite-oxidizing Bacteria, respectively (Francis et al., 2005; Ward, 2011). Initially, the first stage of nitrification, the aerobic ammonia oxidation, was thought to be performed by ammonia-oxidizing bacteria (AOB) of the phylum Proteobacteria (Purkhold et al., 2000). However, the discovery of archaeal homologs of the bacterial genes encoding for ammonia monooxygenase (amo) in marine and terrestrial metagenomes (Venter et al., 2004; Treusch et al., 2005) and in an archaeal isolate (Könneke et al., 2005) led to the conclusion that members of the archaeal phylum Thaumarchaeota (formerly known as Marine Group I Crenarchaeota; Brochier-Armanet et al., 2008; Spang et al., 2010) are the main ammonia oxidizers in the ocean (Stahl and De La Torre, 2012). In contrast, AOB are present only in low abundances in the oceanic water column (Agogue et al., 2008). Recently, archaeal ammonia oxidizers (AOA) were differentiated into two vertically segregated clusters (Hallam et al., 2006; Beman et al., 2008), water cluster A (WCA) and cluster B (WCB). Later, two ecotypes were distinguished according to the environmental ammonia supply rates (Sintes et al., 2013). The high ammonia concentration (HAC) archaeal ammonia oxidizing ecotype, corresponding to WCA, is dominant in epipelagic and upper mesopelagic waters, especially at high latitudes, while the low ammonia concentration (LAC) (<10 nM) ecotype, corresponding to WCB, is dominant in the oligotrophic gyres and in meso- and bathypelagic waters (Sintes et al., 2013, 2016; Santoro et al., 2017). This general depth distribution of these subgroups might vary as described for the HAC ecotype or WCA, with a surface and a deep-water group (Sintes et al., 2016; Bertagnolli and Ulloa, 2017). Ammonia and nitrite oxidation have been observed in OMZs, even at low or non-detectable oxygen concentrations (Fussel et al., 2012; Kalvelage et al., 2013). Nitrite oxidation can exceed ammonia oxidation rates in these oxygen-deficient ecosystems (Fussel et al., 2012; Kalvelage et al., 2013), suggesting an uncoupling between these two processes. Bacterial nitrite oxidizers (NOB) are distributed throughout the OMZs and can account for a substantial proportion of the prokaryotic community in OMZs (Fussel et al., 2012).
The anaerobic ammonia oxidation (anammox) (Thamdrup and Dalsgaard, 2002), in which is oxidized using as electron acceptor and producing N2 gas (Jetten et al., 1998), has been suggested to be an important sink for fixed inorganic nitrogen in the ocean (Codispoti and Christensen, 1985). Globally, this process might be responsible for up to 50% of the N2 gas produced in the oceans (Dalsgaard et al., 2005). Several studies, based on phylogenetic analyses on the 16S rRNA gene (Brandes and Devol, 2002; Kuypers et al., 2005; Hamersley et al., 2007), indicate that this process is performed by members of the bacterial order Planctomycetales (Strous et al., 1999), represented by five Candidatus genera: Candidatus Brocadia (Strous et al., 1998), Candidatus Kuenenia (Schmid et al., 2000), Candidatus Scalindua (Schmid et al., 2003), Candidatus Anammoxoglobus (Kartal et al., 2007) and Candidatus Jettenia (Quan et al., 2008). Hydrazine-oxidoreductase (hzo) is the key enzyme of the anammox process (Shimamura et al., 2007), catalyzing the oxidation of hydrazine (N2H4) to N2. The hzo genes can be grouped into three clusters based on their phylogeny: hzo cluster 1, cluster 2a and 2b (Schmid et al., 2008). Bacteria of cluster 2 exhibit lower ammonia oxidation rates as compared to members of cluster 1 (Kartal et al., 2011). Up to now, few studies have focused on the abundance and community composition of anammox bacteria in the marine water column (Pitcher et al., 2011; Kong et al., 2013) and sediments (Dang et al., 2010, 2013; Shao et al., 2014) based on hzo genes (Schmid et al., 2008). Moreover, the dominance or co-existence of anammox and denitrification processes in oxygen-depleted environments is under debate (Dalsgaard et al., 2005; Bulow et al., 2010; Russ et al., 2014; Bristow et al., 2016). Overall, the balance between anammox and denitrification processes seems to be constrained by the flux and the C/N ratio of available organic matter (Babbin et al., 2014). The bioavailability of organic carbon and organic matter stoichiometry plays a key role in determining the relative contributions of anammox and denitrification to fixed nitrogen removal in the ocean (Dang and Chen, 2017).
Denitrification is an anaerobic respiratory process found in both autotrophic and heterotrophic prokaryotes (Ward et al., 2009), energetically less favorable than aerobic respiration (Deutsch et al., 2001) and predominantly taking place under anaerobic or low oxygen conditions. and are sequentially respired to NO and N2O, resulting in the production of N2. N2O gas has been reported to accumulate at the oxyclines in OMZs as a result of nitrification and denitrification and organic matter oxidation processes (Bange et al., 2001; Kock et al., 2016). Similarly, N2O concentrations in the Gulf of Alaska increase concurrently to the decrease in oxygen concentrations down to 600 m depth (Grundle et al., 2012). Several key enzymes are characteristic for denitrification processes and thus, can be used to assess the community composition and distribution of denitrifiers. Membrane-bound nitrate reductase (nar), encoded by the narH and narG genes, catalyzes the first step of dissimilatory nitrate reduction (Simon and Klotz, 2013). Nitrate reductase has been found in Firmicutes, Enterobacteria, Betaproteobacteria, and Gammaproteobacteria (Petri and Imhoff, 2000) in anoxic environments. The periplasmatic nitrate reductase, encoded by napA, is expressed in the presence of nitrate (Imhoff, 2016) and is highly conserved in chemoautotrophic nitrate-reducing Epsilonproteobacteria (Vetriani et al., 2014). In general, denitrifying Bacteria harbor one of two types of nitrite reductases, either the cytochrome-containing enzyme encoded by nirS or the copper-containing enzyme encoded by nirK (Zumft, 1997; Petri and Imhoff, 2000). Previous studies have indicated a high diversity of Bacteria harboring nir genes within the Proteobacteria phylum (Braker et al., 1998; Huang et al., 2011). Interestingly, the nirK gene has also been found in bacterial and archaeal ammonia oxidizers (Casciotti and Ward, 2001; Treusch et al., 2005; Lund et al., 2012), and the nirS in anammox bacteria (Imhoff, 2016). Thus, it has been suggested that nirK and amo have co-evolved and that carrying both genes allow the cells to adapt to both aerobic and low oxygen environmental conditions (Treusch et al., 2005; Imhoff, 2016).
The aim of this study was to expand our knowledge on the coupling between aerobic and anaerobic ammonia oxidizing and denitrifying communities throughout the water column of the Gulf of Alaska (GoA). The GoA is characterized by the cyclonic Subarctic Alaskan Gyre, which comprises three major current systems (Stabeno et al., 2004). The southern boundary of the gyre is the West Wind Drift flowing eastward, which splits into the equatorward California Current and the poleward Alaska Current. The Alaska Current flows counterclockwise northward from ~48°N along the west coast of North America (Hickey and Royer, 2001) and becomes the narrow, fast flowing Alaskan Stream at the head of the Gulf, following the shelf-break south-westward. The Alaska Current is interspersed with frequent mesoscale anticyclonic eddies and meanders (Stabeno et al., 2004). Three types of eddies, which can persist for years, carry nutrient-rich warm waters offshore into regions with low ambient nutrient levels, and thus, can have a great impact on the chemistry and biology of the GoA ecosystem (Ladd et al., 2009). The basin is also well-known for its oxygen depleted deep waters (Hood and Zimmerman, 1986) and its seasonal deep OMZ (Paulmier and Ruiz-Pino, 2009).
We hypothesized that high latitude deep-ocean oxygen-limited areas might be characterized by bacterial and archaeal nitrifiers and denitrifiers different from the shallow tropical and subtropical OMZs due to distinctly different environmental conditions. To address the coupling or segregation of the different microbial nitrifiers and denitrifies, we used phylogenetic and functional marker genes to assess the distribution of two previously described ecotypes of AOA, amoA-HAC, and amoA-LAC (Sintes et al., 2013), and nirK harboring Archaea and Bacteria by q-PCR. Moreover, we characterized the community composition of Bacteria and Archaea, ammonia oxidizers, denitrifiers, and anammox in epi-, meso-, and bathypelagic waters. Additionally, we determined the environmental factors shaping the abundance and diversity of these communities. Our results support the previously suggested niche-separation of the two AOA ecotypes revealing a strong correlation between those ecotypes and two archaeal denitrification genes, indicating their co-occurrence within the same genome. Moreover, our results point to a larger contribution of the LAC-ecotype in deep ocean OMZs than in other, low latitude OMZs (Bertagnolli and Ulloa, 2017).
Materials and Methods
Study Area and Sampling
Sampling was conducted in the Gulf of Alaska during the DORC (Deep Ocean Refractory Carbon) cruise on board R/V Melville in August 2013. Water samples were collected at 12 stations (Figure 1) at six depths for physico-chemical and biological parameters (bacterial abundance, DNA): the surface layer (5 m), the deep chlorophyll maximum (~50 m), the OMZ (located in the mesopelagic, ~1,000 m), the bathypelagic (at 2,000 and 3,000 m), and the bottom layer (max. depth 5,100 m). Water samples were collected with 12-L Niskin bottles mounted on a CTD (conductivity-temperature-depth) rosette sampler equipped with specific sensors for pressure (Paroscientific Digiquartz pressure transducer), temperature (SBE3plus), conductivity (SBE4C), and dissolved oxygen (SBE43). Salinity was determined on Guildline models 8400B and 8410 Portasal at the Marine Chemistry Laboratory at the University of Washington's School of Oceanography to calibrate the conductivity sensors (Hansell, 2016). Oxygen sensors were calibrated against water samples measured by Winkler titration (Langdon, 2010). Environmental parameters (depth, potential temperature, salinity, oxygen, AOU, nitrite, nitrate, and ammonia) are summarized in Table S1.
Figure 1. Map of the study area in the Gulf of Alaska. Sampling sites occupied during the DORC research cruise in 2013 are indicated by numbered dots.
Inorganic Nutrient Concentrations
Dissolved inorganic nutrient concentrations (, , , , and ) were analyzed as described elsewhere (Hansell, 2016). The full dataset of measured inorganic nutrients can be downloaded from the Biological and Chemical Oceanography Data Management Office website (https://www.bco-dmo.org/dataset/527121).
Abundance of Prokaryotes
Prokaryotic abundance was determined by flow cytometry following an established protocol (Brussaard, 2004). Water samples were fixed with glutaraldehyde (0.5% final concentration), flash-frozen in liquid nitrogen and stored at −80°C until further analysis. The samples were subsequently enumerated on a FACSAria II Cell sorter (Becton Dickinson) after staining with SybrGreen I based on their signature in green fluorescence vs. side scatter (Brussaard, 2004).
DNA Extraction and PCR Amplification
Four to 9 L of water, depending on the expected prokaryotic abundance, were filtered onto 0.22 μm pore-size membrane filter (Millipore, GTTP), flash-frozen in liquid nitrogen and stored at −80°C until further processing in the lab. DNA extraction was performed using the UltraClean Soil DNA Isolation Kit (MoBio Laboratories) following the manufacturer's Advanced User's protocol. DNA extracts were stored at −80°C.
The primers used for amplification of phylogenetic and functional marker genes to generate q-PCR standards and for sequencing are listed in Tables S2, S3, respectively. The thermocycling conditions were chosen as previously described (Tables S2, S3) with slight modifications. The amplification reaction consisted of PicoMaxx Polymerase (2.5 U), 1x final concentration of PicoMaxx Reaction Buffer (Agilent Technologies), 1 μL DNA extract, primers (0.25 μM of amoA primers, 0.5 μM of archaeal nirK primers, and 0.4 μM of bacterial nirK primers, 0.5 μM of each primer set for sequencing purposes), 2 μg bovine serum albumin, 0.25 mM deoxyribonucleoside triphosphate equimolar solution mix, 5 mM MgCl2, made-up to 25 μL with PCR-grade water (Roche). The PCR products were checked by gel electrophoresis on a 2% agarose gel and visualized after staining with SYBR® Gold Nucleic Acid Gel Stain (ThermoFisher) for the right band size. The PCR product was purified using PCRExtract MiniKit (5-PRIME).
Q-PCR analysis was conducted on samples collected at the 12 stations and 6 depths per station (Figure 1). Q-PCR was used to quantify gene abundance of the bacterial recA, the 16S rRNA of Archaea and the 16S rRNA of anammox Bacteria, as well as the abundance of the key functional genes of nitrification and denitrification pathways (archaeal amoA-HAC and amoA-LAC, archaeal nirK-a and nirK-b, and bacterial nirK) using specific primer sets (Table S2). The bacterial nirS gene was amplified with the primer set nirS-1F and nirS-6R (Braker et al., 1998), however, due to the ambiguous PCR product, the data have been excluded from further analysis. The recombinase A (recA) was chosen as a proxy for bacterial abundance, as it is a single-copy gene in Bacteria (Lin et al., 2006). Standards for the different genes were prepared by amplifying a deep-sea sample or previously cloned PCR fragments (in the case of bacterial and archaeal nirK) with the specific primer sets as previously described (Sintes et al., 2013). For archaeal amoA-HAC, DNA from Nitrosopumilus cultures was used to prepare the standard (Sintes et al., 2013). Briefly, the amplified PCR products were checked on a 1% agarose gel stained with SYBR® Gold (Invitrogen) for the right size-band. Subsequently, the target PCR product was purified with the PCR Extract Mini Kit (5Prime). After the quantification of the purified products using a Nanodrop® spectrophotometer, the gene abundance was calculated from the concentration of the purified DNA and the corresponding fragment length. Triplicate 10-fold serial dilutions ranging from 100 to 107 specific gene copies were used to generate an internal quantification standard for each gene. The detection limit ranged between 1 gene copy μL−1 of extract (anammox 16S rRNA), 10 gene copies μL−1 (recA, archaeal nirK-a, archaeal 16S rRNA, bacterial nirK), 100 gene copies μL−1 (HAC-amoA, LAC-amoA), and 1,000 gene copies μL−1 (archaeal nirK-b).
Q-PCR analyses were performed using the LightCycler 480 thermocycler (Roche) equipped with LightCycler 480 gene scanning software (version 1.5, Roche) as previously described (Sintes et al., 2016). The reaction mixtures for each gene contained 1 × LightCycler 480 SYBRGreen® I Master (Roche), primers (0.1 μM of amoA primers, 0.2 μM of archaeal nirK primers, and 0.15 μM of bacterial nirK primers), 1 μL of DNA extract (with nucleic acid concentration ranging between 1.8 and 44.6 ng μL−1) or 1μL of the corresponding standard dilution, and was made-up to 10 μL with PCR-grade water (Roche). Details on used primers and q-PCR thermocycling conditions are listed in Table S2. All standards, environmental samples, and negative controls were run in triplicate in 96-well q-PCR plates (Roche) with optical tape. The specificity of the q-PCR reaction was tested by gel electrophoresis (2% agarose) and by melting curve analysis (65–95°C) to identify unspecific PCR products. PCR efficiency was on average 78.1% for bacterial recA, 83.9% for archaeal 16S rRNA, 92.6% for anammox 16S rRNA, 86.9% for archaeal amoA-HAC, 78.8% for archaeal amoA-LAC, 85.1% for archaeal nirK-a, 98.1% for archaeal nirK-b, and 74.4% for bacterial nirK. Gene abundance was normalized per mL of seawater considering the seawater volume filtered and the extracted DNA volume for each sample assuming 100% filtration and extraction efficiency.
Primer coverage of the HAC- and LAC-amoA was checked against the marine sequences from the most recent curated non-redundant database (Alves et al., 2018) using ARB (Ludwig et al., 2004). Allowing two mismatches, 35% of the amoA marine sequences were targeted by both primer sets, while ~35% of the sequences were only targeted by the HAC- primer set and 9% by the LAC- primer set, resulting in a total primer coverage of 79%.
The correlation between q-PCR quantification of the recA gene abundance and prokaryotic cell enumeration assessed by flow cytometry indicated that recA gene abundance was an appropriate proxy for prokaryotic abundance. A similar correlation between prokaryotic abundance and recA gene abundance as found in this study was also found for the Atlantic (Steiner, 2013).
Sequencing and Bioinformatic Analyses
Next generation sequencing was performed at IMGM Laboratories GmbH (Martinsried, Germany) on an Illumina MiSeq using v3 chemistry. Sequence libraries from three stations (Station 5, 10, and 11) and three depth layers (epi- meso- and bathypelagic) were constructed for the 16S rRNA gene of Bacteria, Archaea and anammox Bacteria, archaeal amoA, archaeal nirK-a and nirK-b, and bacterial nirK genes following the pipeline described previously (Sintes et al., 2016). These stations and depths were selected to cover distinct locations and contrasting environmental conditions (Figure 1, Figure S1). Primers and amplification conditions are listed in Table S3. All samples were barcoded using multiplex identifiers and sequenced together in one run.
The bioinformatics analysis of the 16S rRNA and functional genes followed the standard operating procedure pipeline (https://mothur.org/wiki/MiSeq_SOP) of Mothur (Schloss et al., 2009). The SILVA database (release 132) was used to identify the phylogenetic affiliation of the 16S rRNA gene. Briefly, the reads were quality checked using make.conting script implemented in Mothur. Sequences below 250 bp (below 200 for 16S rRNA of anammox bacteria) or larger than the corresponding fragment size (Tables S2, S3) were excluded. Sequences with ambiguities or more than 8 homopolymers were also excluded. The obtained sequences were further screened for chimeras using vsearch (Rognes et al., 2016) as implemented in the Mothur script chimera.vsearch with default settings. The sensitivity was decreased, however, for 16S rRNA of Bacteria (minimum score to report chimera increased to 5). Operational taxonomic units (OTUs) were defined as sequences with 97% similarity using the average linkage method cutoff implemented in Mothur. Phylogenetic diversity, Chao1, OTU richness, Shannon index of diversity, and the Simpson evenness index were calculated with Mothur using the subsampled OTU table to ensure equal number of OTUs for all samples for each of the different genes sequenced.
The sequence data of the 100 most abundant OTUs from each functional gene were aligned together with environmental and reference sequences from NCBI using MEGA-6 (Tamura et al., 2013). Phylogenetic trees of the functional genes were constructed in MEGA-6 (Tamura et al., 2013) using the Neighbor-Joining method (Saitou and Nei, 1987). Subsequently, the resulting trees were established using iTOL (Letunic and Bork, 2016). Obtained sequences were deposited at the sequence read archive at NCBI, accession number PRJNA507511.
Spearman's rank-order correlation analysis was performed using SPSS v20.0 (SPSS, Inc., Chicago, IL) to determine correlations between the community composition of the different phylogenetic and functional groups with environmental parameters (depth, potential temperature, salinity, apparent oxygen utilization, ammonia, nitrate, and nitrite concentration). Correlations were considered statistically significant at −0.5 > rs > 0.5 and p ≤ 0.05 for all analyses. Single correlation analysis was conducted in SPSS to test the relationship between prokaryotic abundance measured by flow cytometry and the recA gene abundance assessed by q-PCR. Significant differences within phylogenetic and functional genes were assessed by the one-way ANOVA on ranks (Shapiro–Wilk normality test and the Mann–Whitney U-test) and t-test, both performed in Sigmaplot 12.0 (Systat Software, Chicago, IL, USA). The Canonical Correspondence Analysis (CCA) was performed using Past v.3.20 (Hammer et al., 2001) to test the relationship between the prokaryotic communities (bacterial and archaeal) and the environmental variables.
Physico-Chemical Characteristics of the Water Column
Dissolved oxygen concentrations decreased with depth ranging between 225.10 and 298.70 μmol kg−1 in epipelagic waters to 11.80–23.60 μmol kg−1 at 1,000 m depth (Table S1). Consequently, the 1,000 m depth layer complies with the CRIO (CRIterion on O2) definition of an OMZ core (O2 < ~20 μM) (Paulmier and Ruiz-Pino, 2009). Below 1,000 m depth, oxygen concentrations increased toward the near bottom waters ranging between 111.00 and 152.40 μmol kg−1 corresponding to an apparent oxygen utilization (AOU) of 188.36 and 228.80 μmol kg−1. Nitrite concentrations were rather constant throughout the water column (≤0.89 μmol kg−1; Table S1). Nitrate concentrations generally increased with depth (0.23–46.47 μmol kg−1) reaching a maximum within the OMZ (38.88–46.47 μmol kg−1) and decreasing toward the deeper layers (≥33.35 μmol kg−1). Ammonia concentrations were typically highest in the epipelagic waters (≤3.03 μmol kg−1) and decreased with depth reaching minimum concentrations within the OMZ (≤0.01 μmol kg−1), and increasing again toward the bottom waters (≤3.19 μmol kg−1; Table S1).
Total prokaryotic abundance determined by flow cytometry decreased with depth by one order of magnitude from an average of 4.7 × 105 cells mL−1 in the epipelagic to 3.0 × 104 cells mL−1 in bathypelagic waters (ANOVA on Ranks, P < 0.001). Total prokaryotic abundance was positively correlated with temperature, O2 concentration and nitrite concentration (r > 0.5, p = 0.000 for all three variables) and negatively correlated with depth, AOU, and salinity (r > −0.55, p = 0.000 for both variables).
Bacterial abundance, assessed as the recA gene abundance determined by q-PCR, decreased with depth from 2.9 × 106 genes mL−1 in the epipelagic to 1 × 104 genes mL−1 in the bathypelagic waters (Figure 2A, Dataset S1). Archaeal 16S rRNA gene abundance increased from 1 gene mL−1 in the surface waters to 9.5 × 103 genes mL−1 in the bathypelagic realm (Figure 2B). Highest archaeal 16S rRNA gene abundance was found in the OMZ (1.4 × 104 genes mL−1; Figure 2B, Dataset S1). Archaeal 16S rRNA gene abundance was positively correlated with depth, salinity (both r > 0.5, p = 0.000) and AOU (r = 0.850, p = 0.000) and negatively correlated with temperature (r > −0.5, p = 0.000). Anammox bacterial 16S rRNA gene abundance was very low (highest abundance 57 × 10−3 genes mL−1) and without a clear depth distribution (Figure 2C). Prokaryotic abundance determined by flow cytometry was related to recA gene abundance determined by q-PCR (y = 0.05x1.15, r = 0.80, p < 0.001).
Figure 2. Depth profiles of (A) bacterial recA, (B) archaeal 16S rRNA, and (C) anammox 16S rRNA gene abundance throughout the water column of the Gulf of Alaska.
Abundance and Distribution of Archaeal Ammonia Oxidizers
Archaeal amoA-HAC (high ammonia concentration amoA) gene abundance decreased from 0.5 × 102-2.4 × 105 genes mL−1 in epipelagic waters to 0.4–30.9 × 102 genes mL−1 in the OMZ, and 1–12 × 103 genes mL−1 in bathypelagic waters (Figure 3D, Dataset S1). In contrast, archaeal amoA-LAC (low ammonia concentration amoA) gene abundance increased from 0.3 to 12 × 102 genes mL−1 in epipelagic waters up to 7.5–142.7 × 102 genes mL−1 in the OMZ and 1.1–112.5 × 102 genes mL−1 in bathypelagic waters (Figure 3E, Dataset S1). The ratio of amoA-HAC to amoA-LAC gene was ≤ 1 throughout the water column, with the exception of the 50 m depth layer (deep chlorophyll maximum, DCM), where it reached 10–2,837 (data not shown). The ratio of total amoA (i.e., HAC plus LAC) to archaeal 16S rRNA gene abundance showed a similar pattern, with most values close to 1 throughout the water column, but reaching a maximum value of 427 in the DCM (data not shown). The amoA-LAC gene abundance was significantly correlated to archaeal 16S rRNA gene abundance (y = 4.29x0.86, r = 0.88, p = 0.000). Furthermore, archaeal amoA-LAC gene abundance was positively correlated with depth, salinity, (r ≤ 0.53, p < 0.001 for all these variables), and AOU (r = 0.74, p < 0.001). Moreover, archaeal amoA-LAC abundance was negatively correlated with potential temperature (r = −0.51, p < 0.001) and (r = −0.54, p < 0.001). The archaeal amoA-HAC: amoA-LAC ratio was positively correlated with O2 and concentrations (r ≥ 0.5, p < 0.001) and negatively correlated with AOU, and (−0.5 ≥ r, p < 0.001). No significant correlations were found between the abundance of archaeal amoA-HAC genes and the physico-chemical parameters (data not shown).
Figure 3. Depth profiles of functional gene abundances: archaeal nirK-a (A) and nirK-b (B) and bacterial nirK (C) archaeal HAC-amoA (D) and LAC-amoA (E) nitrification genes measured by q-PCR. amoA-HAC, “high-ammonia concentration” archaeal amoA; amoA-LAC, “low-ammonia concentration” archaeal amoA; nirK-a and nirK-b, archaeal nitrite reductase K type a and b; bac nirK, bacterial nitrite reductase K.
Abundance and Distribution of nirK Containing Prokaryotes
Archaeal nirK gene abundance dominated over bacterial nirK genes (Figure 3, Dataset S1). Archaeal nirK-a gene abundance varied from 1 to 1 × 104 genes mL−1 in epipelagic waters and from 1 to 294 genes mL−1 in the OMZ and 1–1 × 103 genes mL−1 in bathypelagic waters (Figure 3A). Archaeal nirK-b gene abundance ranged between 1.9 and 78.7 × 102 genes mL−1 in epipelagic waters, between 1.6 and 65.7 × 102 genes mL−1 in the OMZ and between 0.5 and 49.5 × 102 genes mL−1 in bathypelagic waters (Figure 3B). Overall, the abundance of nirK-a and nirK-b genes was significantly different (P < 0.001, t-test) throughout the water column. Bacterial nirK gene abundance was highest in surface and near-bottom layers, and lowest in the OMZ. Bacterial nirK gene abundance was ≤ 33 genes mL−1 in epipelagic waters, ≤ 3 genes mL−1 in the OMZ and ≤ 30 genes mL−1 in bathypelagic waters (Figure 3C).
The two variants of archaeal denitrification genes co-varied strongly with the two archaeal nitrification genes. Archaeal nirK-a gene abundance strongly correlated with archaeal amoA-HAC gene abundance (Figure S2A) and the abundance of archaeal nirK-b correlated with the amoA-LAC gene abundance (Figure S2B). Total archaeal nirK vs. total amoA (calculated as the sum of nirK-a and nirK-b, and amoA-HAC and amoA-LAC gene abundance, respectively) were significantly correlated (y = 2.38 × 0.99, r = 0.67, p < 0.001). Archaeal nirK-b gene distribution was related to archaeal 16S rRNA gene abundance (r = 0.60, p < 0.001). No significant (r ≥ 0.5, p ≤ 0.01) correlation between denitrification gene abundance and environmental parameters was found.
Prokaryotic Community Composition
Within the Bacteria domain, Alphaproteobacteria (Rhodobacterales, Rhodospirillales, and SAR11 clade) dominated throughout the water column (Figure 4), ranging between 40.9% in the epipelagic (station 10) and 13.6% in the bathypelagic (station 11). Bacteroidetes (Flavobacteriales) decreased with depth, from 20.7% in epipelagic waters (station 11) to 2.9% in bathypelagic waters (station 11). The abundance of Gammaproteobacteria increased with depth from 11.2% in the epipelagic to 28.7% in the OMZ and 26.4% in the bathypelagic (at station 11). Some Gammaproteobacteria orders (Alteromonadales, HOC36, Oceanospirillales, Thiomicrospirales) increased with depth, whereas SAR86 clade, Thiotrichales, Cellvibrionales, and Betaproteobacteriales decreased with depth. Within Actinobacteria, the Actinomarinales decreased with depth (Figure 4) from 11.1% in the epipelagic (station 11) to undetectable in the bathypelagic (station 5 and 10), whereas Microtrichales were more abundant at OMZ (4.0% at station 10) and bathypelagic (3.7% at station 5) than in epipelagic waters. Members of Verrucomicrobiales and Synechococcales (Cyanobacteria) were more abundant in epipelagic waters, up to 1.8 and 1.9%, respectively, at station 5. SAR324 clade (Deltaproteobacteria) (8.0–13.3%), members of the Phycisphaerales from Planctomycetes (up to 1.2%) and SAR202 (Chloroflexi) (up to ~4.1%) and Marinimicrobia (5.4–10.8%) contributed more to the bacterial communities in the OMZ and bathypelagic waters than to epipelagic waters, where SAR324, Phycisphaerales, and SAR202 contributed <0.9% and Marinimicrobia contributed 0.9–3.0%. The abundance of Betaproteobacteria was very low (≤1%) throughout the water column, and betaproteobacterial ammonia oxidizers ranged between undetectable and 0.06% in the 16S rRNA libraries. Nitrospinae bacteria, including NOB, accounted for up to 1.3% of the bacterial community composition in the epipelagic waters at station 5 (Figure 4).
Figure 4. Relative abundance of major bacterial groups at the order level in specific stations and depth layers based on bacterial 16S rRNA gene sequencing. Orders that contribute ≥1% to the community are shown. Members of orders contributing ≤1% to the community are grouped under “Other Bacteria,” “Other Proteobacteria,” “Other Alphaproteobacteria,” and “Other Gammaproteobacteria.” Epi, epipelagic zone; OMZ, oxygen minimum zone; Bathy, bathypelagic zone; St., station.
Candidatus Brocadia (family Brocadiae) were the most abundant members of the anammox Bacteria throughout the water column (Table 1) ranging between 100 and 70% of the total anammox bacterial sequences, while Candidatus Scalindua ranged between 0 and 30% throughout the water column. The Planctomycetales order, to which Brocadia and Scalindua belong, represented between 0 and 0.3% of the 16S rRNA bacterial communities.
Table 1. Contribution of different archaeal families and groups and of different bacterial anammox genera to the archaeal and anammox communities, respectively.
Members of the thaumarchaeal family Nitrosopumilaceae, belonging to Marine group I (MGI), were the dominant Archaea in the OMZ and bathypelagic waters (Table 1) ranging from 54.3 to 66.9% of total archaeal abundance. The abundance of members of Marine group II (MGII) Euryarchaeota decreased with depth from 40.8 to 55.4% in epipelagic waters to 26.6–30.1% in bathypelagic waters (Table 1). Members of the Marine group III (MGIII) Euryarchaeota increased from 0.6 to 0.9% in epipelagic waters to 5.9–6.7% in bathypelagic waters. Marine group IV (MGIV; Halobacteria) and Marine benthic group A (MBGA) were found in the OMZ and bathypelagic waters, ranging between 0.03–0.2% and 0.2–0.5%, respectively (Table 1). Members of Nanoarchaeota were more abundant at depth, reaching 1.5% at station 10.
Spearman's rank correlations of members of the prokaryotic community and the measured physico-chemical parameters are summarized in Table S4. Briefly, members of “other” Alphaproteobacteria and Verrucomicrobiales were positively correlated with nitrite (r ≥ 0.71, p ≤ 0.05). Flavobacteriales, Parvibaculales were positively correlated with ammonia concentrations (r ≥ 0.68, p ≤ 0.05; Figure S3). Members of SAR202 clade, Marinimicrobia, HOC36 correlated negatively only with ammonia concentrations (r ≤ −0.69, p ≤ 0.05). In contrast, members of SAR324, Alteromonadales, and UBA10353 marine group were negatively correlated with nitrite (r ≤ −0.71, p ≤ 0.01) and ammonia concentrations (r ≤ −0.69, p ≤ 0.05). Members of Deltaproteobacteria (SAR324), Gammaproteobacteria (Thiomicrospirales) and “other” Alphaproteobacteria were positively correlated with nitrate (r ≥ 0.71, p ≤ 0.05; Figure S3). In contrast, Actinobacteria, SAR86, and Verrumicrobiales were negatively correlated with nitrate concentrations (r ≤ −0.71, p ≤ 0.05). Members of Halobacteria, unclassified Nanoarchaeaeota, MBGA, and unclassified Thaumarchaeota were negatively correlated with nitrite (r ≤ −0.71, p ≤ 0.05). Halobacteria, Woesearchaeia, and MBGA correlated with ammonia concentrations negatively (r ≤ −0.73, p ≤ 0.05), whereas MGII correlated positively (r ≥ 0.67, p ≤ 0.05; Figure S3). No significant correlations were found between anammox members with nitrite, nitrate nor ammonia concentrations.
Archaeal richness and diversity indexes increased with depth (p ≤ 0.004, Kruskal–Wallis one-way ANOVA on ranks), while remaining fairly stable for Bacteria and anammox across the different depth layers (Table S5).
Composition of Denitrifying and Ammonia Oxidizing Communities
Phylogenetic analysis of the archaeal nirK genes revealed two main clusters for each of the two gene types (Figure 5). Three OTUs dominated the archaeal nirK-a community, A-Otu0001 and A-Otu0002 mainly included sequences from epipelagic and OMZ waters, while A-Otu0003 was present in bathypelagic waters and the OMZ. All representative archaeal nirK-a OTUs were closely related to the Nitrosopumilus genus nirK. Archaeal nirK-b revealed two main clusters dominated by five OTUs. One cluster consisted mainly of sequences from bathypelagic and OMZ waters (with very few sequences from epipelagic layers), while the other cluster consisted exclusively of epipelagic sequences. The OTUs from archaeal nirK-b were only closely related to environmental sequences from the oxygenated epi-, meso- and bathypelagic water of the open ocean and coastal environments but not to isolates (Figure 5).
Figure 5. Phylogenetic tree of the 100 most abundant archaeal nirK OTUs type a and b sampled in the Gulf of Alaska. One representative of sequence group >97% identical is shown; the bar size shows the number of sequences represented by an OTU, and the colors indicate the depth layer. Green tones: epipelagic; light blue tones: oxygen minimum zone; dark blue tones: bathypelagic. Reference sequences as well as environmental sequences from archaeal nirK are included.
Sequences of bacterial nirK genes were only successfully amplified in four samples (Figure S4). Phylogenetic analysis of the bacterial nirK showed a clear stratification with depth. Two OTUs, closely related to the alphaproteobacterial genus Phaeobacter, dominated the epipelagic waters (Otu0002 and Otu0003). The dominant bathypelagic nirK harboring Bacteria were closely related to the genera Nitrosomonas (Otu0001 and Otu0004) and Pseudomonas (Otu0006), or to environmental sequences (Otu0005, Otu0007, Otu0008, and Otu0009; Figure S4).
Archaeal and bacterial nirK richness and diversity generally increased with depth (Table S5). Archaeal nirK-b diversity parameters were significantly lower in epipelagic waters than in the OMZ (number of observed OTUs, p = 0.004; Shannon index, p = 0.032; Chao1 index, p = 0.007; Kruskal–Wallis ANOVA on ranks) and bathypelagic waters (Inverse Simpson index, p = 0.015).
The amoA harboring community was dominated by three OTUs (Figure 6). The two most abundant OTUs (Otu0001 and Otu0002) were affiliated to the amoA-HAC ecotype, closely related to the Nitrosopumilus genus, and dominated in epipelagic waters. Two OTUs (Otu0003 and Otu0006) related to amoA-LAC ecotype dominated in deep waters, and were only related environmental sequences. The diversity of the ammonia oxidizing community (Inverse Simpson index) increased with depth, while richness (Chao1 index, number of observed OTUs) decreased with depth (Table S5). The number of OTUs found in epipelagic waters was significantly larger than in bathypelagic waters (p = 0.011, Kruskal–Wallis ANOVA on rank).
Figure 6. Phylogenetic tree of the 100 most abundant archaeal amoA OTUs recovered from the Gulf of Alaska. One representative of sequence group >97% identical is shown; the bar shows the number of sequences represented by an OTU, and the colors indicate the depth layer (Green tones: epipelagic; light blue tones: oxygen minimum zone; dark blue tones: bathypelagic). Reference sequences as well as environmental sequences from HAC- and LAC-amoA are included.
The eastern subtropical North Pacific is characterized by a pronounced OMZ at 700–1,300 m depth, with oxygen concentration below 20 μM in its core. Currently, this OMZ seasonally extends northwards into the Gulf of Alaska (Paulmier and Ruiz-Pino, 2009), being more intense in spring and fall. It is forecasted that OMZs are expanding and intensifying as a response to anthropogenic impacts (Stramma et al., 2008; Keeling et al., 2010). OMZs are hotspots of microbial activity and particularly of processes involved in the nitrogen cycle (Lam and Kuypers, 2011). Microbial players of the nitrogen cycle have been thoroughly studied in permanent open ocean OMZs, such as the Arabian Sea and the eastern Tropical South Pacific (Pitcher et al., 2011; De Brabandere et al., 2014) and shallow seasonal OMZs (Molina et al., 2010; Galan et al., 2012). Seasonal OMZs, characterized by alternation of environmental conditions stimulate or inhibit several processes involved in the nitrogen cycle (Galan et al., 2017). These seasonal OMZs are particularly suited to study the coupling and dynamics between different microbial players in the nitrogen cycling and to understand the microbial response to expanding OMZs (Wright et al., 2012; Hallam et al., 2017).
Depth Distribution of Bacterial and Archaeal Phylotypes
The relative abundance of Thaumarchaeota in the Gulf of Alaska increased in meso- and bathypelagic waters, while the relative bacterial abundance decreased with depth in agreement with previous studies in the Pacific, Atlantic, and the Mediterranean Sea (Agogue et al., 2008; De Corte et al., 2009; Church et al., 2010; Santoro et al., 2010).
Overall, the bacterial community composition in the Gulf of Alaska exhibited a similar depth stratification in relation to various environmental parameters as previously reported for the Atlantic and Pacific (Delong et al., 2006; Treusch et al., 2009; Agogue et al., 2011). Briefly, the epipelagic waters are inhabited mainly by Rhodobacterales and Flavobacteriales previously reported to respond to phytoplankton blooms (Teeling et al., 2012, 2016; Taylor et al., 2014). The OMZ and bathypelagic waters are characterized by an increase in the contribution of chemoautotrophs such as SAR324, and Planctomycetes, to which anammox Bacteria belong, as well as groups involved in the sulfur cycle (e.g., Marinimicrobia, Thiotrichales, Thiomicrospirales), similarly to coastal oxygen-depleted environments (Aldunate et al., 2018). These groups correlated negatively with oxygen concentrations (Figure S3) and might benefit from gradients in oxygen and redox potential in the OMZ providing different electron acceptors and donors for microbial metabolism (Stevens and Ulloa, 2008).
The archaeal community is dominated by MGI Thaumarchaeota (Nitrosopumilaceae) and MGII Euryarchaeota followed by MGIII Euryarchaeota, Nanoarchaeota, Marine Benthic Group A, and members of the Halobacteria (MGIV). MGI Thaumarchaeota and MGII Euryarchaeota have been reported to be the dominant groups of the marine archaeal community from epipelagic (Karner et al., 2001; Martin-Cuadrado et al., 2015) to bathypelagic waters (López-García et al., 2001; Teira et al., 2006; Li et al., 2015). In contrast, MGIII and MGIV Euryarchaeota are present in low abundance in bathypelagic communities (Fuhrman and Davis, 1997; Galand et al., 2009).
Potential for Nitrification and Denitrification of Archaeal Ecotypes
Thaumarchaeota have been estimated to represent ~20% of the prokaryotic cells in the ocean (Karner et al., 2001), however, they are present in only low abundance in open ocean surface waters (Tolar et al., 2016b). Besides light inhibition (Merbt et al., 2012), it has been suggested that Thaumarchaeota are inhibited by hydrogen peroxide (H2O2), a by-product of photochemical and biological processes in surface waters (Tolar et al., 2016a). Thaumarchaeal cells lack detoxifying enzymes for this highly reactive compound (Tolar et al., 2016b; Bayer et al., 2019a). Additionally, other environmental factors such as ammonia concentration (Sintes et al., 2013) and temperature (Groussin and Gouy, 2011) support the niche differentiation of marine Thaumarchaeota into two ecotypes (Francis et al., 2005; Beman et al., 2008; Sintes et al., 2013, 2016; Luo et al., 2014). The two ecotypes of ammonia oxidizing Archaea, HAC- and LAC-AOA, showed a depth distribution in the Gulf of Alaska similar to that previously reported for the Atlantic (Sintes et al., 2016) and Pacific (Santoro et al., 2017). HAC-AOA dominates in the epipelagic, characterized by micromolar concentrations of ammonium (Table S1), and LAC-AOA in meso- and bathypelagic waters, where ammonium concentrations are in the nanomolar range to below the detection limit (Table S1). These data are in agreement with the findings of Sintes et al. (2013, 2016), Smith et al. (2016), and Santoro et al. (2017) who suggested that these ecotypes are adapted to different ammonium concentrations and supply rates. The location of the OMZ, at around 1,000 m depth, coincides with the increase in the LAC-AOA, as indicated in previous studies in other shallow OMZs (Bertagnolli and Ulloa, 2017). The deeper location of the OMZ in the Gulf of Alaska results in a larger dominance of the LAC-AOA, which comprised 90.9 ± 5.1% of the AOA community as compared with 42% in shallower OMZs (Bertagnolli and Ulloa, 2017), suggesting that oxygen concentration is not the main driver of the distribution of these two ecotypes. However, the environmental factors behind this depth distribution pattern remain unknown.
Notably, the two variants of thaumarchaeal nirK genes (nirK-a and nirK-b) also showed a distinct depth distribution pattern (Figure 3), in agreement with previous findings from the Pacific off Monterey Bay and from the California Current (Lund et al., 2012). Archaeal nirK-a containing cells dominate in epi- and upper mesopelagic waters while archaeal nirK-b containing cells dominate in meso- and bathypelagic waters, similar to previous findings (Lund et al., 2012). This distribution pattern tentatively suggests that archaeal nirK-a is mainly found in the HAC-amoA containing cells, and that archaeal cells containing nirK-b are also harboring LAC-amoA. This notion is further supported by the close affiliation of nirK-a and amoA-HAC to nirK and amoA gene from Nitrosopumilus maritimus, Cand. Nitrosopumilus piranensis, and Cand. Nitrosopumilus adriaticus, while archaeal nirK-b and amoA-LAC sequences are solely associated with environmental sequences (Figures 5, 6). Moreover, the covariation of nirK-a and HAC-amoA abundance, on the one hand, and nirK-b and LAC-amoA abundance, on the other hand, explain the variation in archaeal and bacterial communities over the depth profile in the Gulf of Alaska OMZ (Figure S3). Collectively, these findings reveal that the AOA ecotypes not only vary in their ammonia monooxygenase, but also in their nitrite reductase (Figure S3). Furthermore, this indicates a resource-based niche-partitioning not only based on the affinity to ammonia (Sintes et al., 2013; Smith et al., 2016), but also by their potential to use alternative sources of energy (Qin et al., 2014; Smith et al., 2016). The potentially wider substrate range of LAC-AOA (Water Cluster B, WCB) is further supported by the mismatch between their relatively high abundance in deep waters and the low nitrification rates at these depths, which otherwise would require an extremely long turnover time based on ammonia chemoautotrophy (Newell et al., 2011; Smith et al., 2016). All cultured isolates of AOA (Könneke et al., 2005; Santoro et al., 2015; Bayer et al., 2019b) belong to the HAC-AOA (or water cluster A, WCA), proposed to be mostly obligate chemoautotrophs. However, less is known about LAC-AOA, which have eluded cultivation thus far. Some AOA can use alternative energy sources such as urea (Bayer et al., 2019b) for inorganic carbon fixation, or could exhibit mixotrophic or heterotrophic lifestyles. Potential mixotrophy or heterotrophy have been suggested for Thaumarchaeota cells based on the observation of organic acids incorporation by Archaea (Varela et al., 2008; Clifford et al., 2019) and expression of transporter proteins for organic acids (Bergauer et al., 2018). Recent findings also indicate differences in ecological interactions of the two ecotypes. Preferential association of bacterial phylotypes with one of the AOA ecotypes has been described, such as the specific association between Nitrospina and members of HAC (WCA-like cluster; Reji et al., 2019). This linkage supports a major role of the HAC-ecotype in nitrification (Smith et al., 2014) suggesting that reciprocal feeding between NOB and AOA (Pachiadaki et al., 2017) would be predominantly between NOB and HAC-AOA, with implications for the nitrogen cycle in zones where LAC-AOA are dominant, such as in the deep ocean OMZ. However, further research is needed to elucidate the contribution of heterotrophic or mixotrophic Thaumarchaeota in the ocean.
The role of the nitrite reductase in aerobic ammonia oxidizers is still unclear. Copper containing nitrite reductase, nirK, has been found in both, bacterial and archaeal ammonia-oxidizers (Casciotti and Ward, 2001; Treusch et al., 2005; Bartossek et al., 2010; Lund et al., 2012). It has been suggested that the nirK enables bacterial ammonia oxidizers to tolerate environments with high ammonia and low O2 concentrations (Kozlowski et al., 2014). Recently, a three-step ammonia oxidation pathway was proposed for Bacteria (Caranto and Lancaster, 2017) with two obligate intermediates, hydroxylamine (NH2OH) and nitric oxide (NO). In this proposed pathway, nirK would catalyze the oxidation of NO produced by hydroxylamine oxidoreductase (HAO) in ammonia-oxidizing Bacteria (AOB) to form (Caranto and Lancaster, 2017). Despite the absence of a HAO homolog in Thaumarchaeota (Kerou et al., 2016), the production of these two intermediates in AOA (Vajrala et al., 2013; Kozlowski et al., 2016) suggests that ammonia oxidation in Archaea might also occur via a three-step process (Carini et al., 2018). The reported accumulation of N2O in this zone (Grundle et al., 2012) implies active nitrification/denitrification processes regardless of the low abundance of Bacteria involved in the nitrogen cycle, thus suggesting an important role of AOA.
Bacterial Players on the Nitrogen Cycle in the Gulf of Alaska
Bacterial ammonia oxidizers, including Nitrosomonas and Nitrosococcales, were only present at very low abundances in the 16S rRNA libraries (0–0.06%), suggesting a minor role of bacterial nitrifiers throughout the water column of the Gulf of Alaska. Nitrite oxidizing Bacteria only accounted for ~1% of the bacterial community, in contrast to other OMZs (Fussel et al., 2012).
Denitrifying Bacteria include AOB, which under oxygen-limited conditions can substitute O2 with as an alternative electron acceptor in the “nitrifier denitrification” pathway (Stein, 2011). However, other members of the Bacteria are capable of denitrification. The alphaproteobacterium Phaeobacter gallaeciensis is the closest cultured relative to the two most abundant epipelagic OTUs of denitrifying Bacteria in this study. Phaeobacter gallaeciencis is adapted to colonize surfaces (Thole et al., 2012; Freese et al., 2017). The distribution of bacterial nirK in the Gulf of Alaska does not seem to be related to dissolved oxygen concentrations. However, marine snow and the zooplankton gut provide diffusion limited microenvironments where aerobic respiration results in oxygen depletion (Alldredge and Cohen, 1987) and thus, expand the niche of anoxic metabolism such as denitrification to otherwise well oxygenated ocean waters (Dang and Lovell, 2016; Bianchi et al., 2018). Ganesh et al. (2015) reported an enrichment of transcripts of genes related to denitrification in particle-associated communities relative to free-living communities in the OMZ. In addition, AOB, such as Nitrosomonas eutropha (Betaproteobacteria), a close relative to the dominant denitrifier OTUs in the bathypelagic realm in this study (Figure S4), have been related to a particle-associated life style (Phillips et al., 1999). Marine snow could provide not only an oxygen-depleted microenvironment but also a regular supply of ammonia associated to the decomposition of organic material (Thole et al., 2012). The particle-associated lifestyle of the denitrifying Bacteria could help explaining their low abundance in our samples (<33 nirK genes mL−1) and their increase in the vicinity to the seafloor where resuspension might introduce particles into the water column. However, the bacterial denitrifiers abundance and diversity in the Gulf of Alaska are most probably underestimated, as we could not successfully amplify the nirS-containing organisms, which have been reported as abundant and diverse nitrifiers (Dang et al., 2009; Jayakumar et al., 2013).
Anaerobic Ammonia Oxidizers
Anammox Bacteria were present in low abundance in the Gulf of Alaska, consisting of Candidatus Brocadia followed by Candidatus Scalindua. This is in contrast to previous observations in marine and freshwater environments where Cand. Scalindua was found to be relatively abundant (Penton et al., 2006; Schmid et al., 2007; Van De Vossenberg et al., 2008; Woebken et al., 2008; Villanueva et al., 2014). Cand. Brocadia has previously been reported as the dominant anammox bacterium in waste water treatment plants (Taylor, 2012; Taylor et al., 2014) corresponding to the general preference of anammox Bacteria to environments with high and concentrations (Oshiki et al., 2016). Anammox bacteria can be active at low levels of O2 ranging between 1 and ~10 μmol O2 L−1 (Jensen et al., 2011) and up to ~20 μmol O2 L−1 (Kalvelage et al., 2011; Dalsgaard et al., 2014). Bristow et al. (2016) demonstrated anammox activity at O2 levels as low as 5–30 nM using a highly sensitive oxygen sensor. However, the minimum O2 concentration we measured in our study was 11.8 μmol L−1 thus, in the upper range of the oxygen concentration for anammox activity. Although 16S rRNA genes affiliated to anammox Bacteria were detected in the GoA waters in agreement with previous findings in zones with low but measurable oxygen concentrations (Dalsgaard et al., 2012), anammox Planctomycetes are highly oxygen-sensitive.
Taken together, our results support the broad distribution of archaeal ammonia oxidizers in the ocean's OMZs and the diversity of microbial niches present in OMZs (Bertagnolli and Stewart, 2018). Bacterial denitrifiers harboring the nirK gene are present throughout the water column potentially exhibiting a particle-associated life style, which can provide anoxic or microaerobic microenvironments facilitating oxygen-sensitive metabolisms (Dang and Lovell, 2016). This niche-specific association explains the generally low abundance of bacterial denitrifiers due to the patchiness of marine snow (Silver et al., 1978), and their relatively higher abundances in the epipelagic realm, where marine snow is formed by phytoplankton activity and their exudates, and in bathypelagic waters, where aggregation can lead to the formation of large particles (Bochdansky et al., 2016). Anammox Bacteria are confined to an even narrower niche than denitrifying Bacteria. Their low abundance throughout the water column suggests that they are members of the rare biosphere (Sogin et al., 2006) and potentially, members of a seed bank (Lennon and Jones, 2011), ready to increase in abundance once favorable environmental conditions are established, i.e., the intensification and areal expansion of the OMZ during the fall in this area. Consequently, archaeal ammonia oxidizers play a major role in the nitrogen cycle throughout the water column in the Gulf of Alaska. However, AOA display a distinct ecotype distribution with amoA-HAC/nirK-a and amoA-LAC/nirK-b dominating in epi- and meso- to bathypelagic waters, respectively. Archaeal ammonia oxidizing ecotypes differ not only in their amoA gene but also in their nirK gene, suggesting that this niche separation implicates a profound change in the substrate preferences of these ecotypes.
ES and GH designed the work. SM, DD, and EC performed the research. BB provided Nitrosopumilus cultures. SM, DD, and ES performed the data analysis. SM, DD, GH, and ES wrote the study. All coauthors approved the final version of the manuscript.
Shiptime and laboratory work was supported by the Austrian Science Fund (FWF) to ES (P27696-B22) and to GH (I486-B09, Z194, and P23234-B11) and by the European Research Council under the European Community's Seventh Framework Program (FP7/2007-2013)/ERC (grant agreement No. 268595) to GH. This work was also supported by the FWF DK project on microbial nitrogen cycling (W1257-B20) and the project ARTEMIS (P28781-B21) to GH.
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.
We thank the captain and crew of RV Melville for their support and great atmosphere on board. Thanks to the chief scientist Dennis A. Hansell for cruise planning and support. Thanks to J. Sinninghe Damsté and M. van der Meer from the Department of Marine Microbiology and Biogeochemistry, Royal Netherlands Institute for Sea Research (NIOZ) for granting SM guest status at the MMB. R. Alves provided the filtered marine Thaumarchaeota database used to test primer coverage.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02141/full#supplementary-material
Table S1. Environmental parameters measured in the different stations sampled. Pot. Temp., potential temperature; AOU, apparent oxygen utilization. Dashes indicate no data available.
Table S2. Primer sets and amplification parameters used for the phylogenetic and functional marker genes amplification by qPCR. HAC, “high-ammonia concentration” archaeal amoA gene; LAC, “low-ammonia concentration” archaeal amoA gene; Thau, Thaumarchaeota; nirK-a and nirK-b, archaeal nitrate reductase K type a and type b, respectively.
Table S3. Primer sets and cycling conditions used to amplify bacterial, archaeal, and anammox phylogenetic and functional genes for next-generation sequencing.
Table S4. Spearman's rank correlation coefficients between the different bacterial and archaeal groups and the environmental parameters measured. Only significant results of the bacterial and archaeal community are shown (−0.5 > rs > 0.5). No significant correlations were found for anammox bacteria. Numbers in bold are p ≤ 0.01. Dashes indicate no significant correlation for the specific parameter. AOU, Apparent oxygen utilization; Pot. Temp., potential temperature.
Table S5. Diversity parameters for phylogenetic (based on 16S rRNA) and functional prokaryotic groups. nitrif., nitrification; amoA, archaeal amoA genes; nirK-a and nirK-b, archaeal nitrite reductase K variants a and b; bac nirK, bacterial nitrite reductase K; epi, epipelagic zone; OMZ, oxygen minimum zone; bathy, bathypelagic zone; nseqs, number of sequences sampled for OTU definition; inv. Simpson, inverse of the Simpson index; sobs, number of OTUs observed. Numbers in bold are significantly different (Kruskal–Wallis ANOVA on rank, p < 0.004).
Dataset S1. Average gene abundances and standard deviation determined by q-PCR.
Figure S1. Salinity-temperature diagram of the different stations sampled in the Gulf of Alaska (A). Depth profiles of dissolved oxygen (B), dissolved organic carbon (C), and ammonia (D) concentration in different stations.
Figure S2. Correlation between archaeal denitrification (nirK-a, nirK-b) and nitrification (amoA HAC, amoA LAC) genes. (A) “High-ammonia concentration” archaeal amoA (amoA HAC) vs. archaeal nitrate reductase K variant a (nirK-a); (B) “low-ammonia concentration” archaeal amoA (amoA LAC) vs. archaeal nitrate reductase K variant b (nirK-b).
Figure S3. Canonical correspondence analysis of the bacterial and archaeal phylotypes using the environmental parameters and the phylogenetic and functional gene abundance.
Figure S4. Phylogenetic tree of the 100 most abundant bacterial nirK OTUs sampled in the Gulf of Alaska. One representative of each sequence group >97% identical is shown; the bar shows the number of sequences represented by the OTU, and the color indicates the depth layer. Green tones: epipelagic; dark blue tones: bathypelagic. Reference and environmental sequences are included.
Agogue, H., Brink, M., Dinasquet, J., and Herndl, G. J. (2008). Major gradients in putatively nitrifying and non-nitrifying Archaea in the deep North Atlantic. Nature 456, 788–791. doi: 10.1038/nature07535
Agogue, H., Lamy, D., Neal, P. R., Sogin, M. L., and Herndl, G. J. (2011). Water mass-specificity of bacterial communities in the North Atlantic revealed by massively parallel sequencing. Mol. Ecol. 20, 258–274. doi: 10.1111/j.1365-294X.2010.04932.x
Aldunate, M., De La Iglesia, R., Bertagnolli, A. D., and Ulloa, O. (2018). Oxygen modulates bacterial community composition in the coastal upwelling waters off central Chile. Deep-Sea Res. Part II Top. Stud. Oceanogr. 156, 68–79. doi: 10.1016/j.dsr2.2018.02.001
Alves, R. J. E., Minh, B. Q., Urich, T., Von Haeseler, A., and Schleper, C. (2018). Unifying the global phylogeny and environmental distribution of ammonia-oxidising archaea based on amoA genes. Nat. Commun. 9:1517. doi: 10.1038/s41467-018-03861-1
Bartossek, R., Nicol, G. W., Lanzen, A., Klenk, H. P., and Schleper, C. (2010). Homologues of nitrite reductases in ammonia-oxidizing archaea: diversity and genomic context. Environ. Microbiol. 12, 1075–1088. doi: 10.1111/j.1462-2920.2010.02153.x
Bayer, B., Pelikan, C., Bittner, M. J., Reinthaler, T., Könneke, M., Herndl, G. J., et al. (2019a). Proteomic response of three marine ammonia-oxidizing archaea to hydrogen peroxide and their metabolic interactions with a heterotrophic alphaproteobacterium. mSystems 4:e00181-19. doi: 10.1128/mSystems.00181-19
Bayer, B., Vojvoda, J., Reinthaler, T., Reyes, C., Pinto, M., and Herndl, G. J. (2019b). Nitrosopumilus adriaticus sp. nov. and Nitrosopumilus piranensis sp. nov., two ammonia-oxidizing archaea from the Adriatic Sea and members of the class Nitrososphaeria. Int. J. Syst. Evol. Microbiol. 7, 1892–1902. doi: 10.1099/ijsem.0.003360
Beman, J. M., Popp, B. N., and Francis, C. A. (2008). Molecular and biogeochemical evidence for ammonia oxidation by marine Crenarchaeota in the Gulf of California. ISME J. 2, 429–441. doi: 10.1038/ismej.2007.118
Bergauer, K., Fernandez-Guerra, A., Garcia, J. L., Sprenger, R. R., Stepanauskas, R., Pachiadaki, M. G., et al. (2018). Organic matter processing by microbial communities throughout the Atlantic water column as revealed by metaproteomics. Proc. Natl. Acad. Sci. U.S.A. 115, E400–E408. doi: 10.1073/pnas.1708779115
Bertagnolli, A. D., and Ulloa, O. (2017). Hydrography shapes community composition and diversity of amoA-containing Thaumarchaeota in the coastal waters off central Chile. Environ. Microbiol. Rep. 9, 717–728. doi: 10.1111/1758-2229.12579
Bianchi, D., Weber, T. S., Kiko, R., and Deutsch, C. (2018). Global niche of marine anaerobic metabolisms expanded by particle microenvironments. Nat. Geosci. 11, 263–268. doi: 10.1038/s41561-018-0081-0
Bochdansky, A. B., Clouse, M. A., and Herndl, G. J. (2016). Dragon kings of the deep sea: marine particles deviate markedly from the common number-size spectrum. Sci. Rep. 6:22633. doi: 10.1038/srep22633
Braker, G., Fesefeldt, A., and Witzel, K. P. (1998). Development of PCR primer systems for amplification of nitrite reductase genes (nirK and nirS) to detect denitrifying bacteria in environmental samples. Appl. Environ. Microbiol. 64, 3769–3775.
Brandes, J. A., and Devol, A. H. (2002). A global marine-fixed nitrogen isotopic budget: implications for Holocene nitrogen cycling. Glob. Biogeochem. Cycles 16, 67-61–67-14. doi: 10.1029/2001GB001856
Bristow, L. A., Callbeck, C. M., Larsen, M., Altabet, M. A., Dekaezemacker, J., Forth, M., et al. (2016). N2 production rates limited by nitrite availability in the Bay of Bengal oxygen minimum zone. Nat. Geosci. 10:24. doi: 10.1038/ngeo2847
Brochier-Armanet, C., Boussau, B., Gribaldo, S., and Forterre, P. (2008). Mesophilic crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nat. Rev. Microbiol. 6:245. doi: 10.1038/nrmicro1852
Bulow, S. E., Rich, J. J., Naik, H. S., Pratihary, A. K., and Ward, B. B. (2010). Denitrification exceeds anammox as a nitrogen loss pathway in the Arabian Sea oxygen minimum zone. Deep-Sea Res. Part I 57, 384–393. doi: 10.1016/j.dsr.2009.10.014
Caranto, J. D., and Lancaster, K. M. (2017). Nitric oxide is an obligate bacterial nitrification intermediate produced by hydroxylamine oxidoreductase. Proc. Natl. Acad. Sci. U.S.A. 114, 8217–8222. doi: 10.1073/pnas.1704504114
Carini, P., Dupont, C. L., and Santoro, A. E. (2018). Patterns of thaumarchaeal gene expression in culture and diverse marine environments. Environ. Microbiol. 20, 2112–2124. doi: 10.1111/1462-2920.14107
Casciotti, K. L., and Ward, B. B. (2001). Dissimilatory nitrite reductase genes from autotrophic ammonia-oxidizing bacteria. Appl. Environ. Microbiol. 67, 2213–2221. doi: 10.1128/AEM.67.5.2213-2221.2001
Church, M. J., Wai, B., Karl, D. M., and Delong, E. F. (2010). Abundances of crenarchaeal amoA genes and transcripts in the Pacific Ocean. Environ. Microbiol. 12, 679–688. doi: 10.1111/j.1462-2920.2009.02108.x
Clifford, E. L., Varela, M. M., De Corte, D., Bode, A., Ortiz, V., Herndl, G. J., et al. (2019). Taurine is a major carbon and energy source for marine prokaryotes in the North Atlantic Ocean off the Iberian Peninsula. Microb. Ecol. 78, 299–312. doi: 10.1007/s00248-019-01320-y
Codispoti, L. A., and Christensen, J. P. (1985). Nitrification, denitrification and nitrous oxide cycling in the eastern tropical South Pacific ocean. Mar. Chem. 16, 277–300. doi: 10.1016/0304-4203(85)90051-9
Dalsgaard, T., Stewart, F. J., Thamdrup, B., De Brabandere, L., Revsbech, N. P., Ulloa, O., et al. (2014). Oxygen at nanomolar levels reversibly suppresses process rates and gene expression in anammox and denitrification in the oxygen minimum zone off northern Chile. MBio 5:e01966. doi: 10.1128/mBio.01966-14
Dalsgaard, T., Thamdrup, B., Farías, L., and Revsbech, N. P. (2012). Anammox and denitrification in the oxygen minimum zone of the eastern South Pacific. Limnol. Oceanogr. 57, 1331–1346. doi: 10.4319/lo.2012.57.5.1331
Dang, H. Y., and Chen, C. T. A. (2017). Ecological energetic perspectives on responses of nitrogen-transforming chemolithoautotrophic microbiota to changes in the marine environment. Front. Microbiol. 8:1246. doi: 10.3389/fmicb.2017.01246
Dang, H. Y., Chen, R. P., Wang, L., Guo, L. Z., Chen, P. P., Tang, Z. W., et al. (2010). Environmental factors shape sediment anammox bacterial communities in hypernutrified Jiaozhou Bay, China. Appl. Environ. Microbiol. 76, 7036–7047. doi: 10.1128/AEM.01264-10
Dang, H. Y., Wang, C. Y., Li, J., Li, T. G., Tian, F., Jin, W., 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. doi: 10.1007/s00248-008-9469-5
Dang, H. Y., Zhou, H. X., Zhang, Z. N., Yu, Z. S., Hua, E., Liu, X. S., et al. (2013). Molecular detection of Candidatus scalindua pacifica and environmental responses of sediment anammox bacterial community in the Bohai Sea, China. PLoS ONE 8:e61330. doi: 10.1371/journal.pone.0061330
De Brabandere, L., Canfield, D. E., Dalsgaard, T., Friederich, G. E., Revsbech, N. P., Ulloa, O., et al. (2014). Vertical partitioning of nitrogen-loss processes across the oxic-anoxic interface of an oceanic oxygen minimum zone. Environ. Microbiol. 16, 3041–3054. doi: 10.1111/1462-2920.12255
De Corte, D., Yokokawa, T., Varela, M. M., Agogue, H., and Herndl, G. J. (2009). Spatial distribution of Bacteria and Archaea and amoA gene copy numbers throughout the water column of the Eastern Mediterranean Sea. ISME J. 3, 147–158. doi: 10.1038/ismej.2008.94
Delong, E. F., Preston, C. M., Mincer, T., Rich, V., Hallam, S. J., Frigaard, N.-U., et al. (2006). Community genomics among stratified microbial assemblages in the ocean's interior. Science 311, 496–503. doi: 10.1126/science.1120250
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. U.S.A. 102, 14683–14688. doi: 10.1073/pnas.0506625102
Freese, H. M., Methner, A., and Overmann, J. (2017). Adaptation of surface-associated Bacteria to the open ocean: a genomically distinct subpopulation of Phaeobacter gallaeciensis colonizes Pacific mesozooplankton. Front. Microbiol. 8:1659. doi: 10.3389/fmicb.2017.01659
Galan, A., Molina, V., Belmar, L., and Ulloa, O. (2012). Temporal variability and phylogenetic characterization of planktonic anammox bacteria in the coastal upwelling ecosystem off central Chile. Prog. Oceanogr. 92–95, 110–120. doi: 10.1016/j.pocean.2011.07.007
Galan, A., Thamdrup, B., Saldias, G. S., and Farias, L. (2017). Vertical segregation among pathways mediating nitrogen loss (N-2 and N2O production) across the oxygen gradient in a coastal upwelling ecosystem. Biogeosciences 14, 4795–4813. doi: 10.5194/bg-14-4795-2017
Galand, P. E., Casamayor, E. O., Kirchman, D. L., Potvin, M., and Lovejoy, C. (2009). Unique archaeal assemblages in the Arctic Ocean unveiled by massively parallel tag sequencing. ISME J. 3, 860–869. doi: 10.1038/ismej.2009.23
Ganesh, S., Bristow, L. A., Larsen, M., Sarode, N., Thamdrup, B., and Stewart, F. J. (2015). Size-fraction partitioning of community gene transcription and nitrogen metabolism in a marine oxygen minimum zone. ISME J. 9, 2682–2696. doi: 10.1038/ismej.2015.44
Groussin, M., and Gouy, M. (2011). Adaptation to environmental temperature is a major determinant of molecular evolutionary rates in Archaea. Mol. Biol. Evol. 28, 2661–2674. doi: 10.1093/molbev/msr098
Grundle, D. S., Maranger, R., and Juniper, S. K. (2012). Upper water column nitrous oxide distributions in the Northeast Subarctic Pacific Ocean. Atmos. Ocean 50, 475–486. doi: 10.1080/07055900.2012.727779
Hallam, S. J., Mincer, T. J., Schleper, C., Preston, C. M., Roberts, K., Richardson, P. M., et al. (2006). Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine Crenarchaeota. PLoS Biol. 4:e95. doi: 10.1371/journal.pbio.0040095
Hallam, S. J., Torres-Beltran, M., and Hawley, A. K. (2017). Comment: monitoring microbial responses to ocean deoxygenation in a model oxygen minimum zone. Sci. Data 4:170158. doi: 10.1038/sdata.2017.158
Hamersley, M. R., Lavik, G., Woebken, D., Rattray, J. E., Lam, P., Hopmans, E. C., et al. (2007). Anaerobic ammonium oxidation in the Peruvian oxygen minimum zone. Limnol. Oceanogr. 52, 923–933. doi: 10.4319/lo.2007.52.3.0923
Hansell, D. (2016). CTD and Bottle Data: TOM and Nutrients From R/V Melville Cruise MV1310 in the North Pacific Gulf of Alaska; 48N to 59N and 129W to 153W in 2013 (North Pacific RDOC Project). Biological and Chemical Oceanography Data Management Office (BCO-DMO). Dataset Version 2016-11-28. Available online at: http://lod.bco-dmo.org/id/dataset/527121
Hickey, B., and Royer, T. C. (2001). “California and alaska currents,” in Encyclopedia of Ocean Sciences, eds J. H. Steele, S. H. Thorpe, and K. A. Turekian (San Diego, CA: Academic Press), 368–379. doi: 10.1006/rwos.2001.0352
Hood, D. W., and Zimmerman, S. T. (1986). The Gulf of Alaska: Physical Environment and Biological Resources. Anchorage, AK: U.S. Department of Commerce and U.S. Department of the Interior. doi: 10.5962/bhl.title.60759
Huang, S., Chen, C., Yang, X., Wu, Q., and Zhang, R. (2011). Distribution of typical denitrifying functional genes and diversity of the nirS-encoding bacterial community related to environmental characteristics of river sediments. Biogeosciences 8, 3041–3051. doi: 10.5194/bg-8-3041-2011
Imhoff, J. F. (2016). New dimensions in microbial ecology-functional genes in studies to unravel the biodiversity and role of functional microbial groups in the environment. Microorganisms 4:e19. doi: 10.3390/microorganisms4020019
Jayakumar, A., Peng, X. F., and Ward, B. B. (2013). Community composition of bacteria involved in fixed nitrogen loss in the water column of two major oxygen minimum zones in the ocean. Aquat. Microb. Ecol. 70, 245–259. doi: 10.3354/ame01654
Jensen, M. M., Lam, P., Revsbech, N. P., Nagel, B., Gaye, B., Jetten, M. 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. doi: 10.1038/ismej.2011.44
Jetten, M. S., Strous, M., Van De Pas-Schoonen, K. T., Schalk, J., Van Dongen, U. G., Van De Graaf, A. A., et al. (1998). The anaerobic oxidation of ammonium. FEMS Microbiol. Rev. 22, 421–437. doi: 10.1111/j.1574-6976.1998.tb00379.x
Kalvelage, T., Jensen, M. M., Contreras, S., Revsbech, N. P., Lam, P., Gunter, 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
Kalvelage, T., Lavik, G., Lam, P., Contreras, S., Arteaga, L., Loscher, C. R., et al. (2013). Nitrogen cycling driven by organic matter export in the South Pacific oxygen minimum zone. Nat. Geosci. 6, 228–234. doi: 10.1038/ngeo1739
Kartal, B., Kuypers, M. M., Lavik, G., Schalk, J., Op Den Camp, H. J., Jetten, M. S., et al. (2007). Anammox bacteria disguised as denitrifiers: nitrate reduction to dinitrogen gas via nitrite and ammonium. Environ. Microbiol. 9, 635–642. doi: 10.1111/j.1462-2920.2006.01183.x
Kerou, M., Offre, P., Valledor, L., Abby, S. S., Melcher, M., Nagler, M., et al. (2016). Proteomics and comparative genomics of Nitrososphaera viennensis reveal the core genome and adaptations of archaeal ammonia oxidizers. Proc. Natl. Acad. Sci. U.S.A. 113, E7937–E7946. doi: 10.1073/pnas.1601212113
Kock, A., Arevalo-Martinez, D. L., Loscher, C. R., and Bange, H. W. (2016). Extreme N2O accumulation in the coastal oxygen minimum zone off Peru. Biogeosciences 13, 827–840. doi: 10.5194/bg-13-827-2016
Kong, L., Jing, H., Kataoka, T., Buchwald, C., and Liu, H. (2013). Diversity and spatial distribution of hydrazine oxidoreductase (hzo) gene in the oxygen minimum zone off Costa Rica. PLoS ONE 8:e78275. doi: 10.1371/journal.pone.0078275
Könneke, M., Bernhard, A. E., De La Torre, J. R., Walker, C. B., Waterbury, J. B., and Stahl, D. A. (2005). Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature 437, 543–546. doi: 10.1038/nature03911
Kozlowski, J. A., Price, J., and Stein, L. Y. (2014). Revision of N2O-producing pathways in the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC 19718. Appl. Environ. Microbiol. 80, 4930–4935. doi: 10.1128/AEM.01061-14
Kozlowski, J. A., Stieglmeier, M., Schleper, C., Klotz, M. G., and Stein, L. Y. (2016). Pathways and key intermediates required for obligate aerobic ammonia-dependent chemolithotrophy in bacteria and Thaumarchaeota. ISME J. 10, 1836–1845. doi: 10.1038/ismej.2016.2
Kuypers, M. M. M., Lavik, G., Woebken, D., Schmid, M., Fuchs, B. M., Amann, R., 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. doi: 10.1073/pnas.0502088102
Ladd, C., Crawford, W. R., Harpold, C. E., Johnson, W. K., Kachel, N. B., Stabeno, P. J., et al. (2009). A synoptic survey of young mesoscale eddies in the Eastern Gulf of Alaska. Deep-Sea Res. Part II 56, 2460–2473. doi: 10.1016/j.dsr2.2009.02.007
Langdon, C. (2010). Determination of Dissolved Oxygen in Seawater by Winkler Titration Using the Amperometric Technique. International Ocean Carbon Coordination Project Report N° 14 International CLIVAR Project Office Publication Series N° 134, Version 1, UNESCO, 1–18.
Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245. doi: 10.1093/nar/gkw290
Li, M., Baker, B. J., Anantharaman, K., Jain, S., Breier, J. A., and Dick, G. J. (2015). Genomic and transcriptomic evidence for scavenging of diverse organic compounds by widespread deep-sea archaea. Nat. Commun. 6:8933. doi: 10.1038/ncomms9933
Lin, Z., Kong, H., Nei, M., and Ma, H. (2006). Origins and evolution of the recA/RAD51 gene family: evidence for ancient gene duplication and endosymbiotic gene transfer. Proc. Natl. Acad. Sci. U.S.A. 103, 10328–10333. doi: 10.1073/pnas.0604232103
López-García, P., López-López, A., Moreira, D., and Rodríguez-Valera, F. (2001). Diversity of free-living prokaryotes from a deep-sea site at the Antarctic Polar Front. FEMS Microbiol. Ecol. 36, 193–202. doi: 10.1016/S0168-6496(01)00133-7
Lund, M. B., Smith, J. M., and Francis, C. A. (2012). Diversity, abundance and expression of nitrite reductase (nirK)-like genes in marine thaumarchaea. ISME J. 6, 1966–1977. doi: 10.1038/ismej.2012.40
Luo, H., Tolar, B. B., Swan, B. K., Zhang, C. L., Stepanauskas, R., Ann Moran, M., et al. (2014). Single-cell genomics shedding light on marine Thaumarchaeota diversification. ISME J. 8, 732–736. doi: 10.1038/ismej.2013.202
Martin-Cuadrado, A. B., Garcia-Heredia, I., Molto, A. G., Lopez-Ubeda, R., Kimes, N., Lopez-Garcia, P., et al. (2015). A new class of marine Euryarchaeota group II from the Mediterranean deep chlorophyll maximum. ISME J. 9, 1619–1634. doi: 10.1038/ismej.2014.249
Merbt, S. N., Stahl, D. A., Casamayor, E. O., Marti, E., Nicol, G. W., and Prosser, J. I. (2012). Differential photoinhibition of bacterial and archaeal ammonia oxidation. FEMS Microbiol. Lett. 327, 41–46. doi: 10.1111/j.1574-6968.2011.02457.x
Molina, V., Belmar, L., and Ulloa, O. (2010). High diversity of ammonia-oxidizing archaea in permanent and seasonal oxygen-deficient waters of the eastern South Pacific. Environ. Microbiol. 12, 2450–2465. doi: 10.1111/j.1462-2920.2010.02218.x
Pachiadaki, M. G., Sintes, E., Bergauer, K., Brown, J. M., Record, N. R., Swan, B. K., et al. (2017). Major role of nitrite-oxidizing bacteria in dark ocean carbon fixation. Science 358, 1046–1050. doi: 10.1126/science.aan8260
Penton, C. R., Devol, A. H., and Tiedje, J. M. (2006). Molecular evidence for the broad distribution of anaerobic ammonium-oxidizing bacteria in freshwater and marine sediments. Appl. Environ. Microb. 72, 6829–6832. doi: 10.1128/AEM.01254-06
Petri, R., and Imhoff, J. F. (2000). The relationship of nitrate reducing bacteria on the basis of narH gene sequences and comparison of narH and 16S rDNA based phylogeny. Syst. Appl. Microbiol. 23, 47–57. doi: 10.1016/S0723-2020(00)80045-4
Phillips, C. J., Smith, Z., Embley, T. M., and Prosser, J. I. (1999). Phylogenetic differences between particle-associated and planktonic ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria in the Northwestern Mediterranean Sea. Appl. Environ. Microbiol. 65, 779–786.
Pitcher, A., Villanueva, L., Hopmans, E. C., Schouten, S., Reichart, G. J., and Sinninghe Damste, J. S. (2011). Niche segregation of ammonia-oxidizing archaea and anammox bacteria in the Arabian Sea oxygen minimum zone. ISME J. 5, 1896–1904. doi: 10.1038/ismej.2011.60
Purkhold, U., Pommerening-Röser, A., Juretschko, S., Schmid, M. C., Koops, H.-P., and Wagner, M. (2000). Phylogeny of all recognized species of ammonia oxidizers based on comparative 16S rRNA and amoA sequence analysis: implications for molecular diversity surveys. Appl. Environ. Microbiol. 66, 5368–5382. doi: 10.1128/AEM.66.12.5368-5382.2000
Qin, W., Amin, S. A., Martens-Habbena, W., Walker, C. B., Urakawa, H., Devol, A. H., et al. (2014). Marine ammonia-oxidizing archaeal isolates display obligate mixotrophy and wide ecotypic variation. Proc. Natl. Acad. Sci. U.S.A. 111, 12504–12509. doi: 10.1073/pnas.1324115111
Quan, Z. X., Rhee, S. K., Zuo, J. E., Yang, Y., Bae, J. W., Park, J. R., et al. (2008). Diversity of ammonium-oxidizing bacteria in a granular sludge anaerobic ammonium-oxidizing (anammox) reactor. Environ. Microbiol. 10, 3130–3139. doi: 10.1111/j.1462-2920.2008.01642.x
Reji, L., Tolar, B. B., Smith, J. M., Chavez, F. P., and Francis, C. A. (2019). Differential co-occurrence relationships shaping ecotype diversification within Thaumarchaeota populations in the coastal ocean water column. ISME J. 13, 1144–1158. doi: 10.1038/s41396-018-0311-x
Russ, L., Speth, D. R., Jetten, M. S. M., Op Den Camp, H. J. M., and Kartal, B. (2014). Interactions between anaerobic ammonium and sulfur-oxidizing bacteria in a laboratory scale model system. Environ. Microbiol. 16, 3487–3498. doi: 10.1111/1462-2920.12487
Santoro, A. E., Casciotti, K. L., and Francis, C. A. (2010). Activity, abundance and diversity of nitrifying archaea and bacteria in the central California Current. Environ. Microbiol. 12, 1989–2006. doi: 10.1111/j.1462-2920.2010.02205.x
Santoro, A. E., Dupont, C. L., Richter, R. A., Craig, M. T., Carini, P., Mcilvin, M. R., et al. (2015). Genomic and proteomic characterization of “Candidatus Nitrosopelagicus brevis”: an ammonia-oxidizing archaeon from the open ocean. Proc. Natl. Acad. Sci. U.S.A. 112, 1173–1178. doi: 10.1073/pnas.1416223112
Santoro, A. E., Saito, M. A., Goepfert, T. J., Lamborg, C. H., Dupont, C. L., and Ditullio, G. R. (2017). Thaumarchaeal ecotype distributions across the equatorial Pacific Ocean and their potential roles in nitrification and sinking flux attenuation. Limnol. Oceanogr. 62, 1984–2003. doi: 10.1002/lno.10547
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., and Hollister, E. B. (2009) Introducing MOTHUR: open-source, platform-independent, community-supported software for describing comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09.
Schmid, M., Twachtmann, U., Klein, M., Strous, M., Juretschko, S., Jetten, M., et al. (2000). Molecular evidence for genus level diversity of bacteria capable of catalyzing anaerobic ammonium oxidation. Syst. Appl. Microbiol. 23, 93–106. doi: 10.1016/S0723-2020(00)80050-8
Schmid, M., Walsh, K., Webb, R., Rijpstra, W. I., Van De Pas-Schoonen, K., Verbruggen, M. J., et al. (2003). Candidatus “Scalindua brodae”, sp. nov., Candidatus “Scalindua wagneri”, sp. nov., two new species of anaerobic ammonium oxidizing Bacteria. Syst. Appl. Microbiol. 26, 529–538. doi: 10.1078/072320203770865837
Schmid, M. C., Hooper, A. B., Klotz, M. G., Woebken, D., Lam, P., Kuypers, M. M., et al. (2008). Environmental detection of octahaem cytochrome c hydroxylamine/hydrazine oxidoreductase genes of aerobic and anaerobic ammonium-oxidizing bacteria. Environ. Microbiol. 10, 3140–3149. doi: 10.1111/j.1462-2920.2008.01732.x
Schmid, M. C., Risgaard-Petersen, N., Van De Vossenberg, J., Kuypers, M. M. M., Lavik, G., Petersen, J., et al. (2007). Anaerobic ammonium-oxidizing bacteria in marine environments: widespread occurrence but low diversity. Environ. Microbiol. 9, 1476–1484. doi: 10.1111/j.1462-2920.2007.01266.x
Shao, S. D., Luan, X. W., Dang, H. Y., Zhou, H. X., Zhao, Y. K., Liu, H. T., et al. (2014). Deep-sea methane seep sediments in the Okhotsk Sea sustain diverse and abundant anammox bacteria. FEMS Microbiol. Ecol. 87, 503–516. doi: 10.1111/1574-6941.12241
Shimamura, M., Nishiyama, T., Shigetomo, H., Toyomoto, T., Kawahara, Y., Furukawa, K., et al. (2007). Isolation of a multiheme protein with features of a hydrazine-oxidizing enzyme from an anaerobic ammonium-oxidizing enrichment culture. Appl. Environ. Microbiol. 73, 1065–1072. doi: 10.1128/AEM.01978-06
Silver, M. W., Shanks, A. L., and Trent, J. D. (1978). Marine snow: microplankton habitat and source of small-scale patchiness in pelagic populations. Science 201, 371–373. doi: 10.1126/science.201.4353.371
Simon, J., and Klotz, M. G. (2013). Diversity and evolution of bioenergetic systems involved in microbial nitrogen compound transformations. Bioenergetics 1827, 114–135. doi: 10.1016/j.bbabio.2012.07.005
Sintes, E., Bergauer, K., De Corte, D., Yokokawa, T., and Herndl, G. J. (2013). Archaeal amoA gene diversity points to distinct biogeography of ammonia-oxidizing Crenarchaeota in the ocean. Environ. Microbiol. 15, 1647–1658. doi: 10.1111/j.1462-2920.2012.02801.x
Sintes, E., De Corte, D., Haberleitner, E., and Herndl, G. J. (2016). Geographic distribution of archaeal ammonia oxidizing ecotypes in the Atlantic Ocean. Front. Microbiol. 7:77. doi: 10.3389/fmicb.2016.00077
Smith, J. M., Casciotti, K. L., Chavez, F. P., and Francis, C. A. (2014). Differential contributions of archaeal ammonia oxidizer ecotypes to nitrification in coastal surface waters. ISME J. 8, 1704–1714. doi: 10.1038/ismej.2014.11
Smith, J. M., Damashek, J., Chavez, F. P., and Francis, C. A. (2016). Factors influencing nitrification rates and the abundance and transcriptional activity of ammonia-oxidizing microorganisms in the dark northeast Pacific Ocean. Limnol. Oceanogr. 61, 596–609. doi: 10.1002/lno.10235
Sogin, M. L., Morrison, H. G., Huber, J. A., Mark Welch, D., Huse, S. M., Neal, P. R., et al. (2006). Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc. Natl. Acad. Sci. U.S.A. 103, 12115–12120. doi: 10.1073/pnas.0605127103
Spang, A., Hatzenpichler, R., Brochier-Armanet, C., Rattei, T., Tischler, P., Spieck, E., et al. (2010). Distinct gene set in two different lineages of ammonia-oxidizing archaea supports the phylum Thaumarchaeota. Trends Microbiol. 18, 331–340. doi: 10.1016/j.tim.2010.06.003
Stabeno, P. J., Bond, N. A., Hermann, A. J., Kachel, N. B., Mordy, C. W., and Overland, J. E. (2004). Meteorology and oceanography of the Northern Gulf of Alaska. Cont. Shelf Res. 24, 859–897. doi: 10.1016/j.csr.2004.02.007
Strous, M., Fuerst, J. A., Kramer, E. H. M., Logemann, S., Muyzer, G., Van De Pas-Schoonen, K. T., et al. (1999). Missing lithotroph identified as new planctomycete. Nature 400:446. doi: 10.1038/22749
Strous, M., Heijnen, J. J., Kuenen, J. G., and Jetten, M. S. M. (1998). The sequencing batch reactor as a powerful tool for the study of slowly growing anaerobic ammonium-oxidizing microorganisms. Appl. Microbiol. Biotechnol. 50, 589–596. doi: 10.1007/s002530051340
Taylor, J. D., Cottingham, S. D., Billinge, J., and Cunliffe, M. (2014). Seasonal microbial community dynamics correlate with phytoplankton-derived polysaccharides in surface coastal waters. ISME J. 8, 245–248. doi: 10.1038/ismej.2013.178
Teeling, H., Fuchs, B. M., Becher, D., Klockow, C., Gardebrecht, A., Bennke, C. M., et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611. doi: 10.1126/science.1218344
Teeling, H., Fuchs, B. M., Bennke, C. M., Kruger, K., Chafee, M., Kappelmann, L., et al. (2016). Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. Elife 5:e11888. doi: 10.7554/eLife.11888
Teira, E., Lebaron, P., Van Aken, H., and Herndl, G. J. (2006). Distribution and activity of Bacteria and Archaea in the deep water masses of the North Atlantic. Limnol. Oceanogr. 51, 2131–2144. doi: 10.4319/lo.2006.51.5.2131
Thamdrup, B., and Dalsgaard, T. (2002). Production of N2 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
Thole, S., Kalhoefer, D., Voget, S., Berger, M., Engelhardt, T., Liesegang, H., et al. (2012). Phaeobacter gallaeciensis genomes from globally opposite locations reveal high similarity of adaptation to surface life. ISME J. 6, 2229–2244. doi: 10.1038/ismej.2012.62
Tolar, B. B., Powers, L. C., Miller, W. L., Wallsgrove, N. J., Popp, B. N., and Hollibaugh, J. T. (2016a). Ammonia oxidation in the ocean can be inhibited by nanomolar concentrations of hydrogen peroxide. Front. Mar. Sci. 3:237. doi: 10.3389/fmars.2016.00237
Tolar, B. B., Ross, M. J., Wallsgrove, N. J., Liu, Q., Aluwihare, L. I., Popp, B. N., et al. (2016b). Contribution of ammonia oxidation to chemoautotrophy in Antarctic coastal waters. ISME J. 10, 2605–2619. doi: 10.1038/ismej.2016.61
Treusch, A. H., Leininger, S., Kletzin, A., Schuster, S. C., Klenk, H. P., and Schleper, C. (2005). Novel genes for nitrite reductase and Amo-related proteins indicate a role of uncultivated mesophilic crenarchaeota in nitrogen cycling. Environ. Microbiol. 7, 1985–1995. doi: 10.1111/j.1462-2920.2005.00906.x
Treusch, A. H., Vergin, K. L., Finlay, L. A., Donatz, M. G., Burton, R. M., Carlson, C. A., et al. (2009). Seasonality and vertical structure of microbial communities in an ocean gyre. ISME J. 3, 1148–1163. doi: 10.1038/ismej.2009.60
Vajrala, N., Martens-Habbena, W., Sayavedra-Soto, L. A., Schauer, A., Bottomley, P. J., Stahl, D. A., et al. (2013). Hydroxylamine as an intermediate in ammonia oxidation by globally abundant marine archaea. Proc. Natl. Acad. Sci. U.S.A. 110, 1006–1011. doi: 10.1073/pnas.1214272110
Van De Vossenberg, J., Rattray, J. E., Geerts, W., Kartal, B., Van Niftrik, L., Van Donselaar, E. G., et al. (2008). Enrichment and characterization of marine anammox bacteria associated with global nitrogen gas production. Environ. Microbiol. 10, 3120–3129. doi: 10.1111/j.1462-2920.2008.01643.x
Varela, M. M., Van Aken, H. M., Sintes, E., and Herndl, G. J. (2008). Latitudinal trends of Crenarchaeota and Bacteria in the meso- and bathypelagic water masses of the Eastern North Atlantic. Environ. Microbiol. 10, 110–124.
Venter, J. C., Remington, K., Heidelberg, J. F., Halpern, A. L., Rusch, D., Eisen, J. A., et al. (2004). Environmental genome shotgun sequencing of the Sargasso Sea. Science 304, 66–74. doi: 10.1126/science.1093857
Vetriani, C., Voordeckers, J. W., Crespo-Medina, M., O'brien, C. E., Giovannelli, D., and Lutz, R. A. (2014). Deep-sea hydrothermal vent Epsilonproteobacteria encode a conserved and widespread nitrate reduction pathway (Nap). ISME J. 8, 1510–1521. doi: 10.1038/ismej.2013.246
Villanueva, L., Speth, D. R., Van Alen, T., Hoischen, A., and Jetten, M. S. M. (2014). Shotgun metagenomic data reveals significant abundance but low diversity of “Candidatus Scalindua” marine anammox bacteria in the Arabian Sea oxygen minimum zone. Front. Microbiol. 5:31. doi: 10.3389/fmicb.2014.00031
Ward, B. B., Devol, A. H., Rich, J. J., Chang, B. X., Bulow, S. E., Naik, H., et al. (2009). Denitrification as the dominant nitrogen loss process in the Arabian Sea. Nature 461:78. doi: 10.1038/nature08276
Woebken, D., Lam, P., Kuypers, M. M. M., Naqvi, S. W. A., Kartal, B., Strous, M., et al. (2008). A microdiversity study of anammox bacteria reveals a novel Candidatus Scalindua phylotype in marine oxygen minimum zones. Environ. Microbiol. 10, 3106–3119. doi: 10.1111/j.1462-2920.2008.01640.x
Keywords: ammonia oxidizers, denitrifiers, anammox, archaea, OMZ, Gulf of Alaska
Citation: Muck S, De Corte D, Clifford EL, Bayer B, Herndl GJ and Sintes E (2019) Niche Differentiation of Aerobic and Anaerobic Ammonia Oxidizers in a High Latitude Deep Oxygen Minimum Zone. Front. Microbiol. 10:2141. doi: 10.3389/fmicb.2019.02141
Received: 10 December 2018; Accepted: 30 August 2019;
Published: 13 September 2019.
Edited by:Osvaldo Ulloa, Universidad de Concepción, Chile
Reviewed by:Veronica Molina, Universidad de Playa Ancha, Chile
Hongyue Dang, Xiamen University, China
Copyright © 2019 Muck, De Corte, Clifford, Bayer, Herndl and Sintes. 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) and the copyright owner(s) 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: Eva Sintes, firstname.lastname@example.org