Physiological and genomic characterization of two novel marine thaumarchaeal strains indicates niche differentiation.

Ammonia-oxidizing Archaea (AOA) are ubiquitous throughout the oceanic water column; however, our knowledge on their physiological and ecological diversity in different oceanic regions is rather limited. Here, we report the cultivation and characterization of two novel Nitrosopumilus strains, originating from coastal surface waters of the Northern Adriatic Sea. The combined physiological and genomic information revealed that each strain exhibits different metabolic and functional traits, potentially reflecting contrasting life modes. Strain NF5 contains many chemotaxis-related genes and is able to express archaella, suggesting that it can sense and actively seek favorable microenvironments such as nutrient-rich particles. In contrast, strain D3C is non-motile and shows higher versatility in substrate utilization, being able to use urea as an alternative substrate in addition to ammonia. Furthermore, it encodes a divergent, second copy of the AmoB subunit of the key enzyme ammonia monooxygenase, which might have an additional catalytic function and suggests further metabolic versatility. However, the role of this gene requires further investigation. Our results provide evidence for functional diversity and metabolic versatility among phylogenetically closely related thaumarchaeal strains, and point toward adaptations to free-living versus particle-associated life styles and possible niche differentiation among AOA in marine ecosystems.

Introduction

Ammonia oxidation represents the first and rate-limiting step in nitrification and the only known oxidative biological process in nature that converts reduced inorganic nitrogen species into their oxidized form (Prosser, 1989). For more than a century, ammonia oxidation was thought to be mainly performed by certain lineages of Beta- and Gammaproteobacteria (Purkhold ); however, the discovery of ammonia-oxidizing archaea (AOA) (Könneke ; Treusch ), which now comprise the class Nitrososphaeria (Stieglmeier ), within the phylum Thaumarchaeota (Brochier-Armanet ; Spang ) has dramatically changed our perception of microbial nitrification and nitrogen cycling.

During the past decade, several studies based on 16 S rRNA and ammonia monooxygenase subunit A (amoA) gene sequences of AOA have shown that this archaeal group appears to be extremely diverse and abundant in various marine and terrestrial environments (Francis ; Leininger ; Zhang ), constituting up to 20% of the ocean's picoplankton (Karner ; Teira ). Cultivated thaumarchaeal ammonia oxidizers have been shown to be adapted to oligotrophic conditions, specifically to low ammonium concentrations (nm–μm range) (Martens-Habbena ; Prosser and Nicol, 2012), as generally found in surface waters of the global oceans.

The majority of AOA present in marine ecosystems belongs to a specific subgroup referred to as group 1.1a or Marine Group I, recently defined as the candidate order Nitrosopumilales (Stieglmeier ) which includes the first cultured representative, 'Candidatus Nitrosopumilus maritimus' SCM1, isolated from a tropical seawater aquarium (Könneke ). Organisms closely related to 'Ca. N. maritimus' SCM1 have now been successfully enriched from arctic marine and estuarine sediments (Mosier ; Park ) and were isolated from coastal surface waters (Qin ). So far, all cultured representatives of group 1.1a couple ammonia oxidation to autotrophic carbon fixation, but vary in other metabolic traits such as the ability to utilize urea or the dependence on small amounts of organic compounds to sustain growth (Park ; Qin ). Such distinctive differences can support niche partitioning of AOA populations, potentially influencing nitrification rates in various marine habitats.

Despite the global distribution and ecological importance of members of group 1.1a, studies on their metabolic diversity and detailed genomic analyses are relatively sparse. In this study, two novel thaumarchaeal ammonia oxidizers were obtained from Northern Adriatic coastal waters. The two enriched organisms were designated 'Candidatus Nitrosopumilus sp.' strains NF5 and D3C, and their complete genomes were assembled. Physiological and genomic characterization revealed that both strains differ in their metabolic capabilities as well as in their functional traits, potentially reflecting contrasting life modes.

Materials and methods

Sampling and establishment of enrichment cultures

Seawater was collected in the Northern Adriatic Sea from approx. 0.5 m depth off Piran, Slovenia (45.518 N, 13.568 E) in two consecutive months (November and December 2011). At the time of sampling, seawater temperatures ranged between 13.1 and 15.5 °C and ammonium concentrations varied between 0.12 and 0.48 μm.

Enrichment cultures were initiated by inoculating 10 ml of surface seawater to 90 ml of modified Synthetic Crenarchaeota Medium (SCM) supplemented with KH2PO4, NaHCO3, NH4Cl, Fe-NaEDTA and non-chelated trace element solution, as described by Könneke . The medium also contained 5–10% of filter-sterilized (0.2 μm mesh) and autoclaved Adriatic seawater to reduce isolation stress. The pH of the medium was adjusted to 7.6 by adding 10 ml l−1 of HEPES buffer (1 m HEPES, 0.6 m NaOH). Cultures were incubated aerobically at 17 and 25 °C in the dark without shaking.

Ammonia-oxidizing activity in the enrichment cultures was monitored by measuring NH4+ consumption and NO2− production, applying fluorescence and absorbance spectroscopy methods using the OPA reagent (Holmes ) and Griess reagent (Griess, 1879), respectively. Additionally, archaeal growth was monitored by flow cytometry (as described in Supplementary material).

Subsequent sub-culturing was performed by inoculating fresh medium (5–10% vol/vol inoculum) with cells in mid-exponential phase, which corresponds to a residual NH4+ concentration of approx. 300 μm (Figure 3). As faster growth was observed at 25 °C than at 17 °C, sub-cultures were grown at 25 °C. In the latest stages of the enrichment process, addition of Adriatic Sea water was omitted from the medium. Different antibiotics were alternately added to the medium over several cell passages to further purify the enrichments. The following antibiotics were used: kanamycin (100 μg ml−1), streptomycin (100 μg ml−1), carbenicillin (100 μg ml−1), spectinomycin (100 μg ml−1), ofloxacin (25 μg ml−1) and chloramphenicol (15 μg ml−1).

This procedure resulted in the enrichment of two AOA strains thereafter referred to as NF5 and D3C (strain NF5 derived from samples taken in November 2011 whereas strain D3C derived from samples taken in December 2011). Temperature and pH optima were determined for both strains (see Supplementary material).

Sequencing, genome assembly and analyses

Genomic DNA was extracted from 1 l batch cultures in early stationary phase, and sequenced using a 454/FLX-Titanium sequencer (Roche, MWG Eurofins, Ebersberg, Germany). Both cultures were >90% enriched at the time of sampling for genomic sequencing as indicated by catalyzed reporter deposition fluorescence In Situ hybridization analyses. The total number of obtained reads was 114 359 and 80 867 and the average read length was 622 and 623 bp for the D3C and NF5 enrichment culture, respectively. The obtained reads were length- and quality-filtered (min. length: 150 bp, min. average quality: 25) using Prinseq (Schmieder and Edwards, 2011), and draft genomes were de novo assembled with Newbler 2.6 (Margulies ). The draft genome assembly of the NF5 enrichment culture resulted in 23 contigs (N50=239 081bp) containing 63 618 reads. Ten contigs (500–1500 bp) of the NF5 genome could not be linked to the rest of the genome, suggesting that these are from bacterial contaminants, which was supported by the fact that they only gave bacterial hits in nucleotide homology (Blastn) searches. For the D3C enrichment, 89 704 reads were assembled into 3 contigs (N50=1 617 100bp). The remaining gaps in both genomes were closed via PCR, using specific primers targeting the start and end positions of each contig. The complete genomes were deposited under the accession numbers CP010868 and CP011070.

Annotations of the draft genomes were performed using the MicroScope genome platform (Vallenet ) and manually curated as previously described (Spang ). Synteny statistics were generated by the respective tools of MaGe (Vallenet ) (blastP alignment threshold=min. 35% sequence identity on min. 80% of the length of the smallest protein, gap parameter=5 genes). The circular genome plot (Supplementary Figure S4) was generated using the Circos software (version 0.54) (Krzywinski ). Average nucleotide identity (ANI) values were calculated based on whole-genome sequences using nucleotide homology (Blastn) searches (Altschul ).

Detailed descriptions of all methodological procedures used in this study can be found in the Supplementary material.

Results and Discussion

Enrichment of two novel marine AOA and their morphology

The two novel AOA strains described in this study were enriched from Northern Adriatic coastal surface waters (approx. 0.5 m depth) sampled off Piran, Slovenia. Ammonia consumption in the initial enrichment cultures was first observed after 8 weeks of incubation at 25 °C. Consecutive cell passages in medium containing 1 mm NH4Cl over a period of 2 years, and the addition of various antibiotics resulted in two stable, highly enriched cultures (>99%) (Supplementary Figure S1A and B). A total of 30 archaeal ammonia monooxygenase subunit A (amoA) gene fragments were cloned from each of the two cultures, and the corresponding sequences were nearly identical (>99.5% nucleotide identity) within each culture but differed slightly (94% nucleotide identity) between the two enrichments, suggesting that two different, but closely related AOA strains had been enriched, thereafter referred to as D3C and NF5.

Strains D3C and NF5 are affiliated with the genus Nitrosopumilus, within the provisional order Nitrosopumilales (also known as Thaumarchaeota group I.1a). Based on the concatenated 16 S–23 S rRNA gene phylogeny of all cultivated Thaumarchaeota with sequenced genomes, strain D3C forms a sub-cluster with 'Ca. N. maritimus' SCM1, isolated from a tropical marine aquarium (Könneke ), and 'Ca. N. koreensis' AR1, enriched from Arctic marine sediment (Park ), whereas strain NF5 branches separately from all other cultivated Nitrosopumilus species (Figure 1). This topology is consistent with that of an amoA gene phylogeny including 'Ca. Nitrosopumilus sp.' strains PS0 and HCA1 (Qin ), for which full rRNA gene sequences are not available (Supplementary Figure S2). The amoA sequence of strain NF5 is almost identical (>99% id) to an amoA sequence (KR072695) that was retrieved from a marine snow particle sampled at the same location off the coast of Piran (Vojvoda et al., unpublished).

Figure 1

Phylogenetic tree of concatenated 16S–23S rRNA genes from cultivated Thaumarchaeota with sequenced genomes (16 organisms) showing the affiliation of strains D3C and NF5. Full-length gene sequences were aligned with MAFFT (L-INS-I method) based on Archaea-specific structurally accurate seed alignments (1474 and 3026 aligned positions for 16S and 23S rRNA genes, respectively), and unreliable positions were filtered from the resulting alignments with TCS before concatenation (4443 aligned positions in the concatenated alignment). The tree was calculated by maximum likelihood with IQ-Tree based on the GTR+I+Γ4 model, with ultrafast bootstrap (UFBoot) and SH-aLRT support values inferred from 1000 replicates each (see Supplementary materials and methods for details). Support values ⩾85% are represented on the respective branches by semi-circles color-coded as indicated on the figure.

In electron microscopy, cells of strain NF5 appeared as small, regular rods with a diameter of 0.20–0.25 μm (mean±s.d.: 0.23±0.015 μm; n=50) and a length of 0.59–1.74 μm (mean±s.d.: 0.93±0.26 μm; n=50) and cells of strain D3C ranged in their diameter from 0.20 to 0.25 μm (mean±s.d.: 0.22±0.008; n=50) and in their length from 0.49 to 2.00 μm (mean±s.d.: 0.85±0.28 μm; n=50) (Figure 2). In agreement with previous observations of all other described Nitrosopumilus strains, we did not observe flagella on cells of strain D3C. In contrast, strain NF5 showed, at least occasionally, flagella-like appendages (21 out of 300 cells) with a diameter of 11–14 nm and hence, within the typical size range reported for archaeal flagella (archaella) (Jarrell ). 'Ca. Nitrosopumilus sp.' NF5 is the first representative of the Nitrosopumilus genus shown to possess archaella. Thus far, within marine representatives of the Nitrosopumilales, archaella have only been observed for 'Ca. Nitrosoarchaeum limnia' strains SFB1 and BG20 (Blainey , Mosier ).

Figure 2

Transmission electron microscopy (a, b) and scanning electron microscopy (c, d) images of strains NF5 (left) and D3C (right).

In Northern Adriatic coastal waters, marine snow is common and represents a nutrient-rich environment in the euphotic layer (Kaltenböck and Herndl, 1992), characterized by a patchy distribution of nutrients. Hence, the potential ability of flagellated microbes to sense and actively move toward favorable microhabitats might be beneficial for microorganisms in Northern Adriatic coastal waters, suggesting that strain NF5 could possibly be associated with such particles. Under culture conditions, nutrients are likely to be homogenously distributed, suggesting that the expression of archaella might be less important, potentially also explaining the low number of flagellated cells observed by TEM.

Physiological characterization

Both enrichment cultures converted ammonia into nitrite stoichiometrically, which was accompanied by growth of AOA as measured by flow cytometry (Figure 3). Both cultures showed very similar growth patterns. Typically after 7–12 days of incubation (starting from a 10% inoculum), 1 mm ammonium was consumed and cultures reached cell abundances of 5 × 107–108 cells ml−1. Ammonium uptake was also observed for very low ammonium concentrations ranging between 50 and 100 nm (limit of detection, data not shown), consistent with previous reports of low kinetic constants (Km 132 nm) of ammonium uptake of 'Ca. Nitrosopumilus maritimus' SCM1 (Martens-Habbena ).

Figure 3

Typical growth curve of strains D3C (a) and NF5 (b) in medium containing 1 mm NH4Cl. Triangles: NH4+ concentrations, full circles: NO2− concentrations, open circles: cell abundances. Error bars represent standard deviations of measurements from triplicate cultures.

Generally, cells grew at temperatures between 14 and 34 °C (Supplementary Figure S3A), which is similar to the temperature range reported for 'Ca. N. maritimus' SCM1, and strains HCA1 and PS0 (Qin ). The temperature optimum of strain NF5 was between 30 and 32 °C, whereas cells of strain D3C grew faster at slightly higher temperatures of 32–34 °C. When cultures were exposed to 37 °C, strain D3C maintained both ammonia oxidation and growth whereas strain NF5 maintained ammonia oxidation in the first days of exposure but rapidly decreased in cell abundance over time. In contrast to terrestrial thaumarchaeal enrichment and pure cultures (Tourna ; Lebedeva ), marine representatives of group I.1a have so far not been shown to be able to sustain growth at a temperature above 35 °C (Qin ; Santoro ). Unlike for strains HCA1 and PS0 (Qin ), growth at 10 °C could not be observed for strains D3C and NF5, indicating a potential adaptation to Adriatic surface waters, where temperatures are typically well above 10 °C for most of the year.

Both strains grew well between pH 6.8 and 8.0, although their highest specific growth rates were observed at approx. pH 7.1 (Supplementary Figure S3B). Despite the average pH range in surface waters of 8.0–8.2, all cultured marine AOA described so far have typical growth optima at around pH 7 (Könneke ; Qin ). However, it should also be considered that the process of ammonia oxidation (uptake of NH3/NH4+ and generation of NO2−) possibly reduces the pH outside the cell, potentially creating microenvironments with lower pH that could facilitate growth in natural environments. Additionally, it has recently been hypothesized that some AOA populations have the capacity to adapt to lower pH and that ocean acidification might promote growth of AOA (Qin ).

Neither growth nor nitrite production was observed when cells of strain NF5 were incubated in ammonium-free medium supplemented with 0.5 mm urea, whereas growth of strain D3C was comparable to that with an equivalent amount of ammonium. Growth of strain D3C on 0.5 mm urea was accompanied by the generation of approx. 1 mm nitrite, indicating that this strain can efficiently utilize urea as an alternative substrate for ammonia oxidation (Figure 4).

Figure 4

Growth of strains D3C and NF5 in medium supplemented with 0.5 mm urea. Open circles: cell abundances, full circles: NO2− production. Error bars represent standard deviations of measurements from triplicate cultures.

Archaeal urease genes appear to be predominantly present in marine polar waters (Alonso-Saez ; Connelly ) and were also identified in 'Ca. Nitrosopumilus sediminis' AR2, enriched from Arctic sediments (Park ), and in the psychrophilic sponge symbiont 'Ca. Cenarchaeum symbiosum' (Hallam ). Together with the recently described 'Ca. Nitrosopumilus sp.' PS0 that was enriched from Puget Sound coastal surface waters and has been shown to grow on urea (Qin ), our findings indicate that urea metabolism might also be important in coastal waters where urea concentrations are typically high (Remsen, 1971; Painter ).

Mixotrophy and potential interaction with heterotrophic bacteria

Amendment of the growth medium with alpha-ketoglutaric acid, pyruvic acid and oxaloacetic acid did not have any effect on the growth of either strain, suggesting that strains D3C and NF5 do not depend on organic molecules to sustain growth. However, the dependence on small organic molecules cannot be completely excluded due to the presence of remaining bacterial contaminants that could potentially produce required metabolites.

The co-cultured bacteria belong to the Alphaproteobacteria and share 99% 16 S rRNA gene sequence with Oceanicaulis alexandrii C116-3 (AJ309863). Oceanicaulis spp. belong to the family Hyphomonadaceae, which consists of chemoorganotrophic prosthecate marine bacteria (Strompl ) known to persist under extremely oligotrophic conditions (Abraham and Rohde, 2014), as typically found in most parts of the ocean.

Recently, co-occurrence patterns of Thaumarchaeota with chemoheterotrophic Proteobacteria have been described in the deep chlorophyll maximum of the Southern California Bight (Beman ). It was suggested that in concert, these groups might be involved in the remineralization of organic material and hence, nutrient cycling.

Physiological studies on Nitrososphaera viennensis EN76T revealed that pure cultures require pyruvate to reach growth rates comparable to those of enrichment cultures containing co-cultured Hyphomicrobium vulgare and Mesorhizobium sp. of the alphaproteobacterial order Rhizobiales (Tourna ). A similar observation was made for the marine strains PS0 and HCA1, which require organic carbon for growth (Qin ), indicating possible symbiotic interactions with heterotrophic bacteria in natural environments. Beneficial trophic interactions between autotrophic AOA and heterotrophic bacteria could potentially also explain the apparent ecological success of Thaumarchaeota in oligotrophic environments such as the open ocean.

Genomic characterization

Genomic plasticity

The complete genomes of the two newly characterized AOA strains D3C and NF5 were assembled. When compared with each other, both genomes share an ANI of 84%. Compared with the genome of 'Ca. Nitrosopumilus maritimus' SCM1, the genomes of strains D3C and NF5 show 88% and 81% nucleotide identity, respectively. These ANI values are far below the previously defined threshold of ⩾95% for species definition (Goris ; Mende ; Kim ). Therefore, we propose two new species names for strains NF5 and D3C: 'Candidatus Nitrosopumilus adriaticus' NF5 and 'Candidatus Nitrosopumilus piranensis' D3C.

Both genomes contain several unique genomic regions that are absent in 'Ca. N. maritimus' SCM1 (Walker ). The genome of 'Ca. N. adriaticus' NF5 contains 375 and 395 coding sequences that are not present in 'Ca. N. piranensis' D3C and 'Ca. N. maritimus' SCM1, whereas 'Ca. N. piranensis' D3C contains 319 and 329 coding sequences that are absent in the genomes of 'Ca. N. adriaticus' NF5 and 'Ca. N. maritimus' SCM1, respectively (Supplementary Figure S4). When compared with each other and with 'Ca. N. maritimus' SCM1, ⩽75% of the genes are conserved in syntenic gene clusters (NF5 vs D3C: 73%, D3C vs NF5: 75%, NF5 vs SCM1: 72%, D3C vs SCM1: 74%), despite their highly similar 16 S rRNA gene sequences.

Energy metabolism

Both genomes encode a complete set of genes that have been proposed to be involved in ammonia oxidation and electron transfer (Walker ; Spang ; Stahl and de la Torre, 2012), including ammonia monooxygenase genes amoAXCB and genes for several blue-copper domain-containing proteins and multicopper oxidases. Additionally, two genes for putative ammonium transporters and genes encoding putative copper-containing nitrite reductases (NirK) have been identified in both genomes.

The genomes of 'Ca. Nitrosopumilus adriaticus' NF5 and 'Ca. Nitrosopumilus maritimus' SCM1 contain two highly similar nirK gene copies (95% and 91% identity, respectively), in contrast to all other sequenced genomes of members of the Nitrosopumilales order, which contain one nirK gene copy each. The exact function of NirK in nitrifiers is not yet clear, although nitric oxide (NO), one of the reduction products of NO2−, might possibly have a central role in the ammonia oxidation process of AOA (Klotz and Stein, 2008; Shen Martens-Habbena ). So far, there are no apparent physiological differences related to ammonia oxidation of 'Ca. N. adriaticus' NF5 and 'Ca. N. maritimus' SCM1 as compared with other strains, suggesting a potential additional role of a second NirK copy in those organisms. Another possible function of NirK in nitrifiers could be nitrite detoxification, as proposed for ammonia-oxidizing bacteria (Cho ; Cantera and Stein, 2007). AOA are typically present in environments where nitrite concentrations are far below toxic levels. However, 'Ca. N. maritimus' SCM1 was isolated from an aquarium where nitrite concentrations can occasionally be high. 'Ca. N. adriaticus' NF5 was enriched from Northern Adriatic coastal waters, where marine snow formation occurs regularly during the summer months (Müller-Niklas ) and nitrite concentrations in these aggregates can be up to 10 times higher than in the surrounding seawater (Kaltenböck and Herndl, 1992). The capability to detoxify nitrite could thus be beneficial in such environments.

'Ca. N. piranensis' D3C harbors a second, highly divergent amoB gene ortholog (amoB2), a unique feature among the genomes of all characterized archaeal and bacterial ammonia oxidizers. Although several ammonia-oxidizing Betaproteobacteria harbor multiple amo gene copies, these are invariably identical or highly similar (Norton ) and only duplications of the amoC gene have been described among characterized AOA genomes (Spang ). The amoB2 gene of 'Ca. N. piranensis' D3C is 567 bp long and encodes a predicted mature protein and a signal peptide of 157 and 31 amino acids, respectively. The full polypeptide chain of 'Ca. N. piranensis' D3C AmoB1 has a similar length (189 aa) and is also predicted to include a signal peptide (35 aa). The AmoB2 shares only 48% amino-acid identity with the AmoB1 of 'Ca. N. piranensis' D3C and between 38% and 50% with all other thaumarchaeal AmoB proteins over the full polypeptide chain (between 42% and 55% with thaumarchaeal AmoB proteins for the predicted mature protein). These overall low amino-acid identities reflect the unusual protein sequence of AmoB2, which comprises regions of high similarity with other AmoB proteins, interspaced by unique amino-acid motifs (Supplementary Figure S5). The amoB2 gene is distantly located from the typical amoAXCB gene cluster of marine AOA and is present within a region of low similarity to other Nitrosopumilus genomes. This cluster comprises several genes encoding secreted and/or membrane-bound proteins of unknown function, suggesting that amoB2 might be a component of an unknown membrane protein complex.

Phylogenetic analysis of available full-length archaeal amoB genes shows that the amoB2 gene of 'Ca. N. piranensis' D3C originates from within the order Nitrosopumilales, but inference of more detailed relationships was likely limited by the short length of thaumarchaeal amoB genes (between 567 and 576 bp, including predicted signal peptides) (Supplementary Figure S6). It is unclear whether amoB2 originates from a duplication event within the genome of 'Ca. N. piranensis' D3C, followed by rapid diversification, or results from horizontal gene transfer from another member of the Nitrosopumilales. Previous studies have indeed shown evidence for frequent horizontal gene transfers in marine Thaumarchaeota (Lopez-Garcia ; Brochier-Armanet ; Deschamps ). The long branch of the amoB2 gene in the phylogenetic tree (Supplementary Figure S6) reflects the highly divergent codons encoding amino-acid motifs unique to this protein, despite the several conserved residues shared with other AmoB protein sequences (Supplementary Figure S5).

Extensive biochemical characterization of particulate methane monooxygenases, the only structurally characterized copper-dependent membrane-bound monooxygenase complexes (which include the archaeal and bacterial AMOs), has shown that the particulate methane monooxygenases active site is a dicopper center coordinated by three histidines in the PmoB subunit (Balasubramanian ; Culpepper and Rosenzweig, 2012). Furthermore, recent crystallization of the AmoB of the thaumarchaeon 'Ca. Nitrosocaldus yellowstonii' HL72 indicated that these three conserved histidines coordinate copper in this protein, possibly constituting the catalytic center of the thaumarchaeal AMO (Lawton ). The AmoB2 protein of 'Ca. N. piranensis' D3C also contains the three histidine residues conserved in all other AmoB and nearly all PmoB proteins (Supplementary Figure S5) (Lawton ), suggesting that it might also bind copper and has a catalytic role. Nevertheless, given the lack of functional evidence, the location of the active site of AMO complexes remains uncertain, as earlier studies indirectly suggested that AmoA might be the catalytic subunit of the bacterial AMO (Hyman and Wood, 1985).

In addition to genes involved in ammonia oxidation, 'Ca. N. piranensis' D3C harbors a ure gene cluster consisting of three urease-encoding subunits (ureABC), four urease accessory proteins (ureDGFE) and two urea transporters. The presence of urease confirms our observations that 'Ca. N. piranensis' D3C is able to grow on urea as a sole energy source. Genes for urea metabolism have only been found in a few members of group I.1a (Hallam ; Park ), but appear to be rather conserved among different species (Supplementary Figure S7). Although urea might represent an important alternative substrate for ammonia oxidation, AOA could scavenge ammonia/nitrogen also through cleavage of peptides. Thus far, all described Nitrosopumilales members with sequenced genomes, including 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C, contain an ABC transport system for di-/oligopeptides and genes encoding membrane-bound aminopeptidases (alanyl-, methionyl-, and leucyl-) as well as various endopeptidases.

Central carbon metabolism

As for all Thaumarchaeota that have been characterized thus far, both genomes contain all the enzymes required for a modified 3-hydroxypropionate/4-hydroxybutyrate (3HP/4HB) pathway for the fixation of inorganic carbon as well as all gluconeogenic enzymes for the generation of glucose-6-phosphate during autotrophic growth. The 3HP/4HB cycle has recently been proposed to be the most energy-efficient aerobic pathway for CO2 fixation (Könneke ). Both 'Ca. Nitrosopumilus adriaticus' NF5 and 'Ca. Nitrosopumilus piranensis' D3C contain one putative bicarbonate transporter (Offre ), providing the ability to actively take up inorganic carbon in the form of bicarbonate, which is used to form malonyl-CoA/methylmalonyl-CoA from acetyl-/propionyl-CoA by the acetyl-CoA/propionyl-CoA carboxylase in the first and seventh steps of the 3HP/4HB cycle (Berg ).

In addition to the genomic repertoire supporting autotrophy, genes encoding candidate enzymes for all steps of an oxidative tricarboxylic acid (TCA) cycle are present in 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C. However, as recently described in Spang , it remains unclear whether Thaumarchaeota are able to exploit classical TCA bypass reactions or use an incomplete TCA cycle exclusively for biosynthetic purposes. There is no conclusive evidence for the operation of a reverse TCA cycle, as Thaumarchaeota, including 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C, harbor a putative citrate/citryl-coA lyase but lack citryl-coA synthetase which together catalyze the ATP-dependent cleavage of citrate into acetyl-CoA and oxaloacetate, representing the key step of the rTCA cycle (Aoshima , 2004b).

'Ca. N. adriaticus' NF5 contains a gene encoding a NAD(P)-dependent malic enzyme, which is present in the genomes of 'Ca. Nitrosopumilus salaria' BD31 (Mosier ), 'Ca. Nitrosoarchaeum limnia' SFB1 (Blainey ) and members of the order Nitrososphaerales (Spang ), but appears to be absent in all other Nitrosopumilus species. This alternative route for PEP (phosphoenolpyruvate) synthesis using a NAD(P)-dependent malic enzyme and pyruvate phosphate dikinase is energetically less favorable than the direct conversion from oxaloacetate via ATP-dependent PEP carboxykinase (Siebers and Schönheit, 2005). The presence of this alternative pathway, however, likely provides higher flexibility in the utilization of the central metabolite pyruvate.

Motility and chemotaxis

In contrast to all other characterized representatives of the Nitrosopumilus genus, 'Ca. Nitrosopumilus adriaticus' NF5 encodes a number of genes associated with archaeal flagella (archaella) assembly and chemotaxis. The thaumarchaeal flagella apparatus shares genetic features with those of Euryarchaeota and Crenarchaeota. As typically found in Crenarchaeota, Thaumarchaeota lack flaCDE genes (Blainey ; Spang ; Lebedeva ), which are thought to be responsible for switching rotation of the archaella. In contrast, both Thaumarchaeota and Euryarchaeota contain the preflagellin peptidase FlaK, which has been shown to specifically process pre-archaellins, not being able to process pre-pilins like the more versatile prepilin-like peptidase PibD found in Crenarchaeota (Szabo ; Jarrell ).

The archaella-associated gene cluster in 'Ca. N. adriaticus' NF5 comprises five flagellin genes (flaB1–5), the putative flagellar accessory genes flaFGHIJ and one gene encoding the preflagellin peptidase FlaK. The identified chemotaxis-related genes include the chemotaxis proteins cheCABW, three methyltransferases (cheR), and several response regulators (cheB and cheY), most of which are located adjacent to the fla gene cluster.

'Ca. N. adriaticus' NF5 lacks a gene encoding the deamidase CheD, which is located between cheC and cheR in 'Ca. Nitrosoarchaeum limnia' SFB1 (Blainey ) and 'Ca. Nitrososphaera gargensis' Ga9.2 (Spang ) (Supplementary Figure S7). Deamidation was proposed to be necessary for receptor activation, and without CheD, receptors are undermethylated resulting in poor kinase activation (Rosario ). Although Archaea usually harbor cheD genes, the function of CheD in some Bacteria, such as Escherichia coli, is carried out by CheB (Szurmant and Ordal, 2004). Additionally, the 'Ca. N. limnia' SFB1 genome only contains a fragment of cheD (Blainey ), suggesting that some Thaumarchaeota may not require CheD for chemotaxis, and possess alternative methylation pathways, or have different requirements for methylation.

The gene organization of archaella-associated genes is rather conserved among Thaumarchaeota (Supplementary Figure S7), with the only exception being the number of flagellins/archaellins: one in 'Ca. N. gargensis' Ga9.2 (Spang ), three in N. viennensis EN76 (Kerou et al., in preparation), four in 'Ca. N. limnia' SFB1 (Blainey ) and five in 'Ca. N. adriaticus' NF5 (this study). The number, organization and amino-acid sequences of chemotaxis-related genes are more diverse than archaella-associated genes (Supplementary Figure S7), possibly indicating adaptations in signal sensing in different environments (soil, hot spring, estuarine and marine).

Flagella appear to be favorable attachment sites for bacteriophages (Rakhuba ), and it has recently been shown that the Campylobacter jejuni bacteriophage NCTC 12673 is able to specifically recognize flagellin modified with acetamidino-modified pseudaminic acid (Javed ). In Archaea, post-translational modification of archaellins is widespread and the amount of N-linked glycosylation sites appears to be variable (up to 16 in flaB2 of Methanothermococcus thermolithotrophicus) (Jarrell ).

A higher diversity of flagellins, as present in 'Ca. N. adriaticus' NF5 which contains five different archaellins, may be associated with a potentially higher flexibility of flagellin modifications that could possibly be beneficial for reducing virus infections.

Protection against UV-induced DNA damage

The genomes of 'Ca. Nitrosopumilus adriaticus' NF5 and 'Ca. Nitrosopumilus piranensis' D3C encode putative cyclobutane pyrimidine dimer photolyases (CPD), which may confer greater resistance to ultraviolet radiation (UVR). Photolyases are known to directly repair UVR-induced DNA damage, and, in some cases, are also able to stimulate excision repair nucleases (Ozer ). In Northern Adriatic surface waters, biologically effective UV radiation reaches down to 10 m depth (Herndl ) and UVR may therefore represent a potential stressor for AOA inhabiting surface waters. This observation is supported by single cell genomic data showing that epipelagic Thaumarchaeota, in contrast to mesopelagic populations, harbor putative photolyases (Luo ).

Photolyases belong to the broad cryptochrome/photolyase family, widespread in Bacteria, Archaea and Eukarya, although in Thaumarchaeota they have only been found in 'Ca. Nitrosoarchaeum limnia' strains SFB1 and BG20 (Blainey ; Mosier ), 'Ca. Nitrosopelagicus brevis' CN25 (Santoro ) and 'Ca. Nitrosotenuis uzonensis' N4 (Lebedeva ). Comprehensive phylogenetic analyses show that the putative photolyases of 'Ca. N. piranensis' D3C and 'Ca. N. adriaticus' NF5 cluster within a distinct Thaumarchaeota-specific clade, which, together with Crenarchaeota (Sulfolobales), Euryarchaeota (Thermoplasmatales and Methanosarcinales), Mimivirus and diverse bacterial phyla, represent a previously unrecognized lineage within the cryptochrome/photolyase family (Supplementary Figure S8). The only biochemically and structurally characterized photolyase within this lineage is the CPD class I of Sulfolobus tokodaii (Fujihashi ), which shares 43% and 41% amino-acid identity with that of 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C, respectively. Additionally, the related photolyase of Sulfolobus acidocaldarius functions in vivo in a light-dependent manner and not as part of a default DNA repair mechanism as, at least, some photolyases (Sakofsky ). This further indicates that the related photolyases of 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C are likely involved in the adaptation of the organisms to epipelagic waters. Considering their low similarity with characterized photolyases, it is nevertheless possible that these enzymes might have additional light-dependent regulatory functions, as shown for other CPD I involved in both UVR-damage repair and development regulation (Bayram ).

Provisional classification

We propose the following candidate status:

Nitrososphaeria class, Nitrosopumilales order, Nitrosopumilaceae fam., 'Nitrosopumilus adriaticus' sp. nov. Etymology: L. adj. m. nitrosus: nitrous (nitrite producer); L. adj. m. pumilus: dwarf; L. adriaticus, obtained from the Adriatic Sea.

Nitrososphaeria class, Nitrosopumilales order, Nitrosopumilaceae fam., 'Nitrosopumilus piranensis' sp. nov. Etymology: L. adj. m. nitrosus: nitrous (nitrite producer); L. adj. m. pumilus: dwarf; L. piranensis, obtained off the Piran coast.

Locality: Northern Adriatic coastal surface waters (Piran, Slovenia).

Diagnosis: Chemolithoautotrophic nitrifiers of the domain Archaea, closely related to 'Ca. Nitrosopumilus maritimus' SCM1. 'Ca. Nitrosopumilus piranensis' D3C can alternatively utilize urea as a substrate. Cells appear as straight rods with a diameter of 0.20–0.25 μm and a length of 0.49–2.00 μm ('Ca. Nitrosopumilus piranensis' D3C) and 0.59–1.74 μm ('Ca. Nitrosopumilus adriaticus' NF5). 'Ca. Nitrosopumilus adriaticus' NF5 cells occasionally appear to be flagellated. The species names allude to the habitat of the organisms.

Conclusions

Our results provide new insight into the functional diversity and metabolic versatility of marine Thaumarchaeota. Although 'Ca. Nitrosopumilus adriaticus' NF5 and 'Ca. Nitrosopumilus piranensis' D3C are phylogenetically closely related and were obtained from the same location, both strains appear to exhibit different metabolic traits and functional adaptations (Table 1), potentially reflecting contrasting life strategies.

Table 1

Genomic features of marine Thaumarchaeota with complete genomesa

Organisma	N. adriaticus NF5	N. piranensis D3C	N. maritimus SCM1	N. koreensis AR1	N. sediminis AR2	C. symbiosum A	N. brevis CN25
Habitat	Surface water	Surface water	Aquarium gravel	Arctic sediment	Arctic sediment	Sponge symbiont	Open ocean
Genome size (Mb)	1.80	1.71	1.65	1.64	1.69	2.05	1.23
G+C content (%)	33.4	33.8	34.2	34.2	33.6	57.4	33.2
CDS	2184	2161	1796	1890	1974	2014	1445
16S–23S rRNA genes	1	1	1	1	1	1	1
Number of tRNAsb	43 (45)	43 (44)	44	38	37	45	42
Nitrite reductase genes (nirK)	2	1	2	1	1	0	1
Urease	−	+	−	−	+	+	−
Flagellac	+	−	−	−	−	−	−
Photolyasec	+	+	−	−	−	−	+

Abbreviation: CDS, coding DNA sequences.

'Ca. Nitrosopumilus maritimus' SCM1 (CP000866), 'Ca. Nitrosopumilus koreensis' AR1 (CP003842), 'Ca. Nitrosopumilus sediminis' AR2 (CP003843), 'Ca. Cenarchaeum symbiosum' A (DP000238), 'Ca. Nitrosopelagicus brevis' CN25 (CP007026).

Number in brackets includes putative tRNAs. It should be noted that the number of tRNAs might be underrepresented as archaeal tRNAs sometimes possess non-canonical introns.

Present in the incomplete genomes of 'Ca. Nitrosoarchaeum limnia' strains SFB1 and BG20.

The presence of motility and chemotaxis genes in the genome of 'Ca. N. adriaticus' NF5 indicates that it might actively sense and move toward more favorable microenvironments in marine surface waters such as marine snow. Northern Adriatic coastal waters are known for large marine snow formation (Müller-Niklas ). These aggregates represent a heterogeneous, nutrient-rich environment (Kaltenböck and Herndl, 1992) and potentially, a favorable niche for flagellated microbes. In addition, the amoA sequence of 'Ca. N. adriaticus' NF5 is identical to an amoA sequence retrieved from marine snow sampled at the same location (Vojvoda et al., unpublished), supporting the hypothesis that this organism might be associated with marine snow-type particles.

In contrast, 'Ca. N. piranensis' D3C is non-motile and shows more versatility than 'Ca. N. adriaticus' NF5 in substrate utilization, as shown by its ability to grow on urea as an alternative substrate. Although easily accessible ammonia seems to be preferentially taken up, the potential to utilize urea may be beneficial when ammonium is in short supply. We suggest that urea uptake by marine AOA might not solely be predominant in polar waters as recently proposed (Alonso-Saez ; Connelly ), but might also represent an important alternative energy source for AOA in temperate coastal waters. Moreover, the presence of an additional distinct amoB subunit gene (amoB2) and respective genomic context suggests a potential for the use of additional alternative substrates and/or additional regulatory mechanisms of ammonia oxidation.

It has recently been shown that extremely closely related (>99% ANI) bacterioplankton populations have different behavioral strategies to interact with microscale nutrient patches in the ocean (Yawata ). Whereas one population colonizes particles, the other population rapidly migrates to new nutrient sources, indicating that spatiotemporal structuring contributes to fine-scale ecological differentiation. Our results indicate that 'Ca. N. adriaticus' NF5 and 'Ca. N. piranensis' D3C occupy different niches in marine coastal waters, corresponding to a more particle-associated life mode versus a free-living life mode, respectively. Niche complementation leads to a reduction of competition and improves resource utilization (Cordero and Polz, 2014), which might be beneficial for survival under changing environmental conditions, potentially explaining the remarkable diversity and ecological success of AOA in the global oceans.

Cultivation of marine Thaumarchaeota often follows the culturing conditions established for 'Ca. Nitrosopumilus maritimus' SCM1 (for example, in this study) (Könneke ). The relatively high ammonium concentrations (0.25–1 mm) and/or high temperatures (25–30 °C) may thus be responsible for the selection of AOA strains from the genus Nitrosopumilus in cultivation studies, despite the much greater diversity identified in environmental surveys (Francis ; Sintes ). This is further illustrated by the fact that the recently cultivated 'Ca. Nitrosopelagicus brevis' CN25, which represents a distinct lineage within the order Nitrosopumilales widespread in oceanic waters, grows optimally at ~22 °C (Santoro ).

Given the functional and genomic diversity observed among closely related AOA and the possible bias introduced by current culturing approaches, it is likely that the current range of cultivated AOA greatly underestimates their actual metabolic and functional diversity in the ocean. Thaumarchaeota are ubiquitous and constitute a large fraction of the microbial biomass in the oceans. Therefore, a broader knowledge about the functional adaptations of their populations on different taxonomic levels will be essential to better understand their role in the nutrient dynamics of these ecosystems.