# Identifying Potential Mechanisms Enabling Acidophily in the Ammonia-Oxidizing Archaeon "Candidatus Nitrosotalea devanaterra"

Laura E. Lehtovirta-Morley, Luis A. Sayavedra-Soto, Nicolas Gallois, Stefan Schouten, Lisa Y. Stein, James I. Prosser, Graeme W. Nicol

<figure><img src="image_1.png" /></figure>

► To cite this version:

Laura E. Lehtovirta-Morley, Luis A. Sayavedra-Soto, Nicolas Gallois, Stefan Schouten, Lisa Y. Stein, et al.. Identifying Potential Mechanisms Enabling Acidophily in the Ammonia-Oxidizing Archaeon "Candidatus Nitrosotalea devanaterra". Applied and Environmental Microbiology, 2016, 82 (9), pp.2608 - 2619. 10.1128/AEM.04031-15. hal-03906404

HAL Id: hal-03906404

https://hal.science/hal-03906404

Submitted on 5 Jan 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Identifying Potential Mechanisms Enabling Acidophily in the
Ammonia-Oxidizing Archaeon “Candidatus Nitrosotalea devanaterra”

Laura E. Lehtovirta-Morley,ª Luis A. Sayavedra-Soto,ᵇ Nicolas Gallois,ᶜ Stefan Schouten,ᵈ Lisa Y. Stein,ᵉ James I. Prosser,ª Graeme W. Nicolª,f

Institute of Biological and Environmental Sciences, University of Aberdeen, Aberdeen, United Kingdomª; Department of Botany and Plant Pathology, Oregon State University, Corvallis, Oregon, USAᵇ; Université Blaise Pascal, Clermont-Ferrand, Franceᶜ; NIOZ Royal Netherlands Institute for Sea Research, Department of Marine Organic Biogeochemistry, and Utrecht University, Den Burg, Texel, The Netherlandsᵈ; Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canadaᵉ; Laboratoire Ampère, École Centrale de Lyon, L'Université de Lyon, Ecully, Francef

Ammonia oxidation is the first and rate-limiting step in nitrification and is dominated by two distinct groups of microorganisms in soil: ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB). AOA are often more abundant than AOB and dominate activity in acid soils. The mechanism of ammonia oxidation under acidic conditions has been a long-standing paradox. While high rates of ammonia oxidation are frequently measured in acid soils, cultivated ammonia oxidizers grew only at near-neutral pH when grown in standard laboratory culture. Although a number of mechanisms have been demonstrated to enable neutrophilic AOB growth at low pH in the laboratory, these have not been demonstrated in soil, and the recent cultivation of the obligately acidophilic ammonia oxidizer “Candidatus Nitrosotalea devanaterra” provides a more parsimonious explanation for the observed high rates of activity. Analysis of the sequenced genome, transcriptional activity, and lipid content of “Ca. Nitrosotalea devanaterra” reveals that previously proposed mechanisms used by AOB for growth at low pH are not essential for archaeal ammonia oxidation in acidic environments. Instead, the genome indicates that “Ca. Nitrosotalea devanaterra” contains genes encoding both a predicted high-affinity substrate acquisition system and potential pH homeostasis mechanisms absent in neutrophilic AOA. Analysis of mRNA revealed that candidate genes encoding the proposed homeostasis mechanisms were all expressed during acidophilic growth, and lipid profiling by high-performance liquid chromatography–mass spectrometry (HPLC-MS) demonstrated that the membrane lipids of “Ca. Nitrosotalea devanaterra” were not dominated by crenarchaeol, as found in neutrophilic AOA. This study for the first time describes a genome of an obligately acidophilic ammonia oxidizer and identifies potential mechanisms enabling this unique phenotype for future biochemical characterization.

Ammonia oxidation is integral to the global nitrogen cycle and is performed by ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA). AOA are among the most abundant phyla on Earth (1) and frequently outnumber AOB in the environment (2). Ammonia concentration (3) and soil pH (4) have been suggested as drivers of niche differentiation in terrestrial ammonia oxidizers. The effect of soil pH on ammonia oxidation is of particular interest since gross nitrification rates in acidic soils (pH of <5.5), which represent ~50% of Earth's potentially arable land, are similar to those in neutral soils (5). Until recently, growth of all cultivated prokaryotic ammonia oxidizers in batch culture was considered possible only at a pH of >6.5. This view was challenged by the initial cultivation (6) and subsequent isolation (7) of the obligately acidophilic AOA “Candidatus Nitrosotalea devanaterra,” which grows in the pH range 4.0 to 5.5. Members of the *Nitrosotalea* lineage are abundant and widely distributed in acidic soils globally (4).

Inhibition of ammonia oxidizer growth and activity at low pH is poorly understood, and several mechanisms have been suggested, including substrate limitation and product toxicity. Kinetic studies of ammonia oxidation in cell suspensions and extracts of *Nitrosomonas europaea* suggest that ammonia (NH₃), rather than ammonium (NH₄⁺), is the substrate for ammonia monooxygenase (AMO), which catalyzes the first step in ammonia oxidation (8). In contrast, the preferred substrate (NH₃/NH₄⁺) of AOA has not been identified. The pKₐ for NH₄⁺ is 9.25, and ammonia availability decreases with pH (e.g., NH₃/NH₄⁺ ratios are 1:200 and 1:20,000 at pH 7 and 5, respectively). Growth on

NH₃, rather than NH₄⁺, is also supported by lack of growth of AOB in liquid culture below pH 6.5. Although AOB growth occurs under acid conditions in culture (e.g., growth of ureolytic strains on urea [9] or within biofilms [10] or aggregates [11]), there is little evidence that these mechanisms enable AOA growth in acidic soils, nor do they explain acidophilic growth of “Ca. Nitrosotalea devanaterra” in laboratory culture.

All known bacterial and archaeal ammonia oxidizers are autotrophs (12, 13). AOB fix CO₂ via the Calvin cycle (12), whereas AOA fix HCO₃⁻ via the hydroxypropionate/hydroxybutyrate pathway (13). Most inorganic carbon exists as CO₂ at low pH, and the low HCO₃⁻ concentration may require additional adaptations in “Ca. Nitrosotalea devanaterra.”

In addition to low NH₃ and HCO₃⁻ availability, acidic pH

Received 17 December 2015 Accepted 13 February 2016

Accepted manuscript posted online 19 February 2016

Citation Lehtovirta-Morley LE, Sayavedra-Soto LA, Gallois N, Schouten S, Stein LY, Prosser JI, Nicol GW. 2016. Identifying potential mechanisms enabling acidophily in the ammonia-oxidizing archaeon "Candidatus Nitrosotalea devanaterra." *Appl Environ Microbiol* 82:2608–2619. doi:10.1128/AEM.04031-15.

Editor: C. Vieille, Michigan State University

Address correspondence to Laura E. Lehtovirta-Morley, llehtovirta@abdn.ac.uk.

Supplemental material for this article may be found at http://dx.doi.org/10.1128/AEM.04031-15.

Copyright © 2016 Lehtovirta-Morley et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
presents a further challenge to cellular pH homeostasis. Almost all known prokaryotic acidophiles maintain a cytoplasmic pH higher than that in the extracellular environment (14). Extreme acidophiles usually have a reverse membrane potential (inside positive), and the proton motive force is constituted by the transmembrane pH gradient. There are several mechanisms of pH adaptation in acidophiles. In some, reverse membrane potential is achieved by uptake of cations (typically potassium), and most acidophile genomes encode a high number of cation transporters (15). The entry of protons is prevented by altering the membrane composition or modifying the cell surface structures (16, 17). Some acidophiles have proton pumps coupled to the electron transfer chain that actively remove protons from the cytoplasm (18). In addition, protons can be scavenged by buffering, e.g., by carbonic anhydrase in *Helicobacter pylori* (19), or metabolism, e.g., by arginine- and glutamate decarboxylase-based acid resistance mechanisms (20). Since small organic acids can function as uncouplers by passing through the membrane before releasing protons and acidifying the neutral pH cytoplasm, most extreme acidophiles are heterotrophs (15, 21). Acidophiles typically contain a large proportion of secondary transporters coupled to ion gradients (22). In addition, DNA and protein repair genes are prevalent in extreme acidophile genomes and are responsible for the rapid repair required after damage by low pH (15).

"Ca. Nitrosotalea devanaterra" is able to grow and oxidize ammonia within the pH range 4 to 5.5 and is unable to grow at neutral pH (6). This physiology makes "Ca. Nitrosotalea devanaterra" unique among all bacterial and archaeal ammonia oxidizers and indicates adaptations distinct from those of other cultured ammonia oxidizers. There are several major unresolved questions regarding acidophilic ammonia oxidation, as follows. (i) How does "Ca. Nitrosotalea devanaterra" overcome low NH₃ concentrations? Specifically, can this organism utilize NH₄⁺, is the active site of AMO facing the cytoplasm or the periplasm and are there other N metabolism genes that can explain the distinct physiology of "Ca. Nitrosotalea devanaterra"? (ii) How does "Ca. Nitrosotalea devanaterra" fix carbon under acidic conditions where the HCO₃⁻ concentration is low? (iii) Is cytoplasmic pH homeostasis of "Ca. Nitrosotalea devanaterra" achieved by mechanisms similar to those in other acidophiles? The aim of this study was therefore to examine the "Ca. Nitrosotalea devanaterra" genome for evidence of specific adaptations in N and C metabolism and to determine whether the genome contained genes involved in pH homeostasis mechanisms found in other model acidophiles.

## MATERIALS AND METHODS

**DNA isolation.** "Candidatus Nitrosotalea devanaterra" Nd1 was grown in pure culture as previously described (7). Biomass from 10 liters of the culture was concentrated using a Pelicon XL tangential-flow filter cassette (pore size, 0.22 µm) (Merck Millipore, Billerica, MA, USA), followed by pelleting of cells by centrifugation (15 min at 18,000 × g). High-molecular-weight genomic DNA was extracted from a cell pellet by chemical lysis (23). Briefly, cells were lysed in the presence of proteinase K (final concentration, 100 µg ml⁻¹) and SDS (0.5% [wt/vol]), followed by addition of cetyltrimethylammonium bromide (CTAB) in NaCl solution (1% CTAB [wt/vol] and 0.7 M NaCl [final concentrations]) in a final volume of 780 µl. Proteins were removed using an equal volume of phenol-chloroform-isoamyl alcohol (25:24:1) and DNA precipitated in the presence of 7.5 µg linear acrylamide with a 0.7 volume of 100% isopropanol. After centrifugation, the pellet was washed with 70% ethanol and resuspended in Tris-EDTA (TE) buffer. PacBio sequencing was necessary for genome

closure, and as larger quantities of DNA (10 µg) were required, approximately 40 liters of culture was harvested by tangential-flow filtration and the extraction was performed using a DNeasy blood and tissue kit (Qiagen, Venlo, Netherlands). A minor modification to the manufacturer's protocol was made: DNA was eluted fives times with 50 µl prewarmed distilled water (dH₂O) and concentrated using a vacuum concentrator (Eppendorf, Hamburg, Germany). The molecular weight of DNA was assessed by agarose gel electrophoresis, and the concentration was quantified using a NanoDrop spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and a PicoGreen Quant-iT kit (Invitrogen, Carlsbad, CA, USA).

**Genome sequencing, assembly, and deposition.** DNA was sequenced using an Illumina MiSeq sequencer (flow cell v3 with 300-bp paired-end and 2-kb mate pair libraries) and with a PacBio RSII (GTAC, Constance, Germany) that provided 50× and 10× genome coverage, respectively. Assembly of Illumina reads was performed with CLC (CLC bio) and that of PacBio reads with Celera 8.1 (24) and SPAdes 3.1.1 (25). Illumina reads were first quality trimmed, several assembly parameters were tested, and finally the assembly was performed with a word size of 45, a bubble size of 98, and a minimum contig length of 1,000 bp. PacBio reads were error corrected using the PBcr pipeline in Celera and subsequently assembled. Neither approach alone was sufficient to close the genome, and a hybrid assembly of PacBio and Illumina reads was performed with SPAdes. Different assemblies were compared and validated using MAUVE software (26). The entire genome was annotated in MaGe (27) based on a combination of bioinformatics tools, including homology-, structure-, and synteny-based approaches. Other AOA genomes are included in MaGe, which relies on the PkGDB relational database with data from both public data banks and manually curated genomes. Manual curation of the "Ca. Nitrosotalea devanaterra" genome was performed in MaGe, with emphasis on genes absent in other AOA, genes located at synteny breaks, genes characterized in other AOA, and functions involved in major metabolic pathways. Functional categorization into clusters of orthologous genes (COGs) and protein families (Pfams) was performed with COGnitor and HMMPfam software implemented in MaGe, respectively. tRNA predictions were made using tRNAscan-SE (28) and ARAGORN (29). Transporters were annotated using TransAAP implemented in TransportDB (30). Average nucleotide identity (ANI) and tetranucleotide frequency analyses were carried out in JSpecies (31). z scores of the nucleotide word frequency analysis were used to build the matrix for the principal-component analysis in PAST software (32) and SigmaPlot v12 (Systat Software Inc., London, United Kingdom). To construct Venn diagrams, the presence of open reading frames (ORFs) in archaeal genomes was identified using reciprocal BLASTP searches (selected cutoffs were an E value of 10⁻³ and 30% amino acid identity) implemented in the stand-alone BLAST-2.2.25+ suite (33). The cutoff for absence and presence of coding sequences (CDS) in comparison with other archaeal and bacterial genomes was selected as 30% DNA identity with 80% CDS coverage, as recommended in the MaGe guidelines. Synteny analysis was performed in MaGe (27).

**Identification of genes potentially linked to acidophily.** To identify candidate genes enabling acidophily, the "Ca. Nitrosotalea devanaterra" genome was compared to those of previously characterized archaeal and bacterial acidophiles to detect gene homologues conserved between "Ca. Nitrosotalea devanaterra" and acidophiles and absent in other AOA (cut-off, 30% DNA identity with 80% CDS coverage). However, as "Ca. Nitrosotalea devanaterra" represents an uncharacterized genus, many ORFs had no database hits in RefSeq, Swiss-Prot, and TrEMBL. Many ORFs also had homology only to uncharacterized proteins in databases. In addition, some (novel) acid tolerance mechanisms may be a result of convergent evolution and not share genetic homology. Therefore, a literature search was performed on pH adaptation strategies of acidophilic microorganisms, as well as on biochemical limitations of ammonia oxidizers with regard to pH, and potential functional analogues were manually identified from the curated "Ca. Nitrosotalea devanaterra" genome.
**Lipid analysis.** The intact polar and core lipids of "Ca. Nitrosotalea devanaterra" were analyzed as described previously (82). For analysis of intact polar lipids, an aliquot of Bligh-Dyer extracts (BDEs) was analyzed with high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (HPLC-ESI-MS/MS). For analysis of core glycerol dibiphytanyl glycerol tetraether lipids (GDGTs), acid hydrolysis was performed on aliquots of BDEs to cleave off polar head groups and release core GDGTs, which were analyzed by HPLC-atmospheric pressure chemical ionization-MS (HPLC-APCI-MS) using a modified procedure (34).

**Sequence alignments, protein structure predictions, and phylogenetic analysis.** Sequence alignments of AmoB, AmoC, and Amt/Rh genes were generated using ClustalW implemented in BioEdit (35) and used to examine active-site conservation and for phylogenetic analysis. For AMO, all AOA *amoB* and *amoC* sequences deposited in GenBank (July 2014) were used in alignments, and the active site of AMO was modeled against PmoB and PmoC of *Methylococcus capsulatus* Bath and other previously characterized model organisms (36–39). This resulted in an alignment of 115 and 147 sequences for *amoB* and *amoC*, respectively. Transmembrane helices and protein topology were predicted initially using TmHMM (40). However, TmHMM failed to detect some of the helices that experimentally have been indicated to be present in AOA AmoB sequences (39), and TMPred (41) and TopPred (42) were subsequently used. Signal peptides were initially predicted by SignalP (43), which also failed to detect signal peptides of AmoB of AOA (39), and prediction was implemented in Phobius (44). Phobius was used additionally for validation of the transmembrane helix predictions. Phobius uses both signal peptides and transmembrane helices in its topology prediction, which requires that amino acid residues immediately downstream from the signal peptide cleavage site must be extracellular.

For Amt/Rh, the database was compiled from a selected set of sequences (total, 89) from GenBank, with the emphasis on sequences that had been previously characterized electrophysiologically and structurally. Phylogenetic analysis of derived Amt/Rh protein sequences used 228 unambiguously aligned positions. Maximum-likelihood analysis was performed with PhyML (45) with invariable sites and eight variable gamma rates modeled. Bootstrap support was also calculated using parsimony and distance analyses (MEGA6 [46]) with 1,000 replicates.

**Transcriptional and physiological response to pH change.** Exponentially growing "Ca. Nitrosotalea devanaterra" cultures (pH 5) were placed into prerinse d and autoclaved dialysis bags (Spectra/Por 4, molecular mass cutoff, 12 to 14 kDa [Spectrum Labs, Rancho Dominguez, CA, USA]) and transferred twice to fresh medium (pH 5) to remove nitrite to avoid inhibition before transfer to media at pH 4.0, 5.0, and 6.0. The pH was recorded at the start and end of the experiments, and mean pH values were 4.2, 5.5, and 6.2 and 4.2, 5.6, and 6.5, respectively. At each time point, three replicate dialysis bags were destructively sampled. Nitrite was measured inside and outside the bag, and 1 ml of cells was used for transcriptional analysis. Cell and *amoA* abundances were estimated at all time points for all replicates. Phenol-ethanol (1:19) stop solution was added immediately after sampling to preserve transcripts before storing at –80°C until RNA extraction. Nucleic acid extraction was performed by bead beating in the presence of buffer containing 0.5% Triton-X, 0.4% N-lauroylsarcosine, 0.4% SDS, 50 mM Tris-HCl (pH 8), and 100 mM EDTA in a total volume of 500 µl, followed by treatment with an equal volume of phenol-chloroform-isoamyl alcohol (25:24:1). Nucleic acids were precipitated in the presence of 7.5 µg linear acrylamide, 0.1 volume 3 M sodium acetate, and an equal volume of 100% ethanol, followed by DNase treatment and reverse transcription (47). Primers were designed for selected transporters and putative pH homeostasis genes using Primer3 (see Table S1 in the supplemental material) (48), and quantitative PCR (qPCR) was performed with QuantiFast Sybr green mix (Qiagen) for all the assays using a two-step cycle: initial denaturation at 95°C for 15 min, followed by 35 cycles of denaturation at 95°C for 10 s, combined annealing/extension at 60°C for 30 s, and fluorescence reading at

72°C, followed by a final extension at 72°C for 10 min, followed by the melt curve from 55°C to 95°C. The quality of qPCR products was verified by melt curve analysis and agarose gel electrophoresis. *amoA* qPCR was performed as previously described (6). Statistical analysis was performed in SigmaPlot v12 (Systat Software Inc., London, United Kingdom). To determine whether there was a significant increase or decrease in transcript abundance over time, linear regression was performed using individual (rather than averaged) values of transcript/gene data as the dependent variable and time as the independent variable. This analysis was additionally performed for nitrite accumulation, assessed with Griess reagent as previously described (7).

**Nucleotide sequence accession number.** The closed genome sequence with full annotations has been deposited in ENA with accession number LN890280.

## RESULTS AND DISCUSSION

**Genome summary.** The closed and complete genome of "Ca. Nitrosotalea devanaterra" is a single circular 1,805,304-bp chromosome with 2,205 predicted genes (Table 1). The genome was automatically annotated prior to manual examination and curation, and a putative function was successfully assigned to 58.8% of ORFs and seven small RNA (smRNA) genes, 40 tRNA genes, and 1 copy of each rRNA gene (5S, 16S, and 23S), with putative introns in various tRNA sequences and the 23S rRNA gene. Tetranucleotide frequency analysis indicates that the "Ca. Nitrosotalea devanaterra" genome is distinct from those of other AOA genera and that ~26% of ORFs of "Ca. Nitrosotalea devanaterra" are absent from other soil AOA genomes (see Fig. S1 in the supplemental material). Approximately 15% (364) of ORFs of "Ca. Nitrosotalea devanaterra" are unique when all aquatic AOA are included in this comparison. The distinction of genus *Nitrosotalea* from other AOA is supported by the low average nucleotide identity (ANI) scores between "Ca. Nitrosotalea devanaterra" and group 1.1a and 1.1b representatives (66.5 to 68.0% and 62.5 to 62.9%, respectively) (see Fig. S1 in the supplemental material).

"Ca. Nitrosotalea devanaterra" has major metabolic pathways that are conserved in all AOA and discussed elsewhere (49, 50), including biosynthesis of amino acids, lipids, and sugars and central C metabolism (see Table S2 in the supplemental material). In addition, "Ca. Nitrosotalea devanaterra" contains an Ni-Fe hydrogenase and genes for flagellar motility, gas vacuoles, and phosphorus utilization, which are present in some but not all AOA (see Tables S2 and S3 in the supplemental material). It is, however, beyond the scope of this article to discuss all of these genes individually, and this article highlights those features that could be associated with the adaptation of "Ca. Nitrosotalea devanaterra" to low-pH environments.

**Similarities between "Ca. Nitrosotalea devanaterra" and other acidophiles.** The fully annotated "Ca. Nitrosotalea devanaterra" genome contains 364 ORFs without homologues in neutrophilic AOA genomes, of which 45 were shared with characterized bacterial and archaeal acidophiles (see Table S4 in the supplemental material). No homologues were found exclusively in all acidophiles, and only few were present in most acidophiles: Na⁺/solute symporter (NDEV_1297), 2 major facilitator superfamily (MFS) transporters (NDEV_1231 and NDEV_1448), 2 natural resistance-associated macrophage protein (NRAMP) divalent cation transporters (NDEV_1079 and NDEV_1085), an FK506 binding protein (FKBP)-type peptidyl-prolyl *cis-trans* isomerase (NDEV_0529), conserved hypothetical proteins (NDEV_0373 and NDEV_1669), 2 transcriptional regulators (NDEV_0570 and
<table><caption>TABLE 1 Genomic features of selected thaumarchaea</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="5">Value</th></tr><tr><th>"Ca. Nitrosotalea devanaterra" Nd1</th><th>"Ca. Nitrosopumilus koreensis" MY1</th><th>"Ca. Nitrosopumilus maritimus" SCM1</th><th>"Ca. Nitrososphaera viennensis" EN76</th><th>"Ca. Nitrososphaera gargensis" Ga9.2</th></tr></thead><tbody><tr><td>Length (Mb)</td><td>1.81</td><td>1.61</td><td>1.65</td><td>2.53</td><td>2.83</td></tr><tr><td>Coding density (%)</td><td>92.1</td><td>89.99</td><td>91.7</td><td>87.0</td><td>83.4</td></tr><tr><td>No. of predicted CDS</td><td>2,205</td><td>1,957</td><td>1,967</td><td>3,123</td><td>3,997</td></tr><tr><td>% of annotated CDS</td><td>58.8</td><td>62.1</td><td>56.0</td><td>46.7</td><td>43.1</td></tr><tr><td>No. of:</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>16-23S rRNA operons</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Separate 5S rRNAs</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>tRNAs</td><td>40</td><td>42</td><td>44</td><td>44</td><td>40</td></tr><tr><td>Other RNAs</td><td>7</td><td>3</td><td>1</td><td>5</td><td>0</td></tr><tr><td>Primary transporter genes</td><td>51</td><td>38</td><td>45</td><td>66</td><td>59</td></tr><tr><td>Secondary transporter genes</td><td>42</td><td>41</td><td>31</td><td>44</td><td>54</td></tr><tr><td>Other transporter genes</td><td>14</td><td>14</td><td>18</td><td>22</td><td>26</td></tr><tr><td>Ratio of secondary to primary transporter genes<sup>a</sup></td><td>0.8:1</td><td>1.1:1</td><td>0.7:1</td><td>0.7:1</td><td>0.9:1</td></tr></tbody></table>

NDEV_1462), and several subunits of archaeal A-type ATP synthase (complex V) which correspond to the membrane-bound A₀ domain and the central stalk (NDEV_1999 [atpI], NDEV_2002 [atpD], NDEV_2005 [atpF], and NDEV_2006 [atpC]) (51) (see Table S5 in the supplemental material). Of these, cation transporters may function in cation uptake, ATP synthase may function in proton translocation, and FKBP-type peptidyl-prolyl cis-trans isomerase may facilitate protein folding. This modest conservation probably reflects the fact that bacterial and archaeal acidophiles have several strategies for pH adaptation (15). In addition, “Ca. Nitrosotalea devanaterra” may possess mechanisms enabling acidophily that are distinct from those of other acidophiles. The genome, with the emphasis on genes absent in other AOA, was therefore further investigated to identify potential mechanisms contributing to pH homeostasis and facilitating growth of AOA with low concentrations of NH₃ and HCO₃⁻.

**Ion transport.** Extreme acidophiles contain a large number of secondary transporters (22), but “Ca. Nitrosotalea devanaterra,” surprisingly, possesses a similar number of secondary transporters to the number possessed by other AOA (Table 1; see Table S5 in the supplemental material) and has proportionally more primary transporters than secondary transporters. The genome of “Ca. Nitrosotalea devanaterra” contains a variety of cation transporters, with predicted uptake of K⁺, Na⁺, and divalent cations (Fig. 1). Na⁺/solute symporter (NDEV_1297) is absent in other AOA and is one of the few genes consistently found across acidophilic archaeal genomes (see Table S4 in the supplemental material). “Ca. Nitrosotalea devanaterra” lacks the typical CPA2 family Na⁺/H⁺ exchangers of other AOA and has two CPA1 family Na⁺/H⁺ exchangers (NDEV_1447 and NDEV_1587), which share low identity (19.8%) (see Table S6 in the supplemental material). “Ca. Nitrosotalea devanaterra” has two NRAMP family transporters (NDEV_1085 and NDEV_1443) that are absent in other AOA and which may transport divalent cations, e.g., Mn²⁺, Fe²⁺, Zn²⁺, Cd²⁺, and Co²⁺.

The genome carries genes encoding two P-type ATPases, *kdpABCD* (NDEV_1461, NDEV_1460, NDEV_1459, and

NDEV_1455) and *mgtA* (NDEV_1079), predicted for K⁺ and Mg²⁺ transport, respectively (Fig. 1 and 2). While *mgtA* is not found in AOA except “Ca. Nitrosotalea devanaterra” (see Table S4 in the supplemental material), the *kdp* gene cluster is syntenic in many prokaryotes, e.g., proteobacteria, *Nitrososphaera viennensis*, and acidophilic archaea (Fig. 2A). “Ca. Nitrosotalea devanaterra” also lacks the *trk* potassium transporter gene conserved in many AOA genomes (see Table S6 in the supplemental material). Potassium is typically a critically important solute for maintaining the reverse membrane potential in acidophiles (15), and it is therefore surprising that the “Ca. Nitrosotalea devanaterra” genome does not contain genes for potassium transporters which are absent in all neutrophilic AOA.

Several of the “Ca. Nitrosotalea devanaterra” genes potentially involved in pH/ion homeostasis are colocalized on the chromosome (*mgtA*, *mntH*, a carbonic anhydrase gene [NDEV_1082], a universal stress protein gene [NDEV_1080], and an α-acetolactate decarboxylase gene [NDEV_1077]). Interestingly, the region is flanked by a transposase gene, located at a syntenic break, and the gene cluster shares a degree of synteny and homology with bacterial genomes (Fig. 2B; see Table S7 in the supplemental material).

**Scavenging protons by metabolism.** Several bacterial acid resistance mechanisms are based on cytoplasmic proton consumption by various decarboxylases (20). All AOA contain an arginine decarboxylase (NDEV_1165) that is involved in arginine and polyamine metabolism but also in arginine-based acid resistance in bacteria. Due to its presence in neutrophilic AOA, it is an unlikely candidate for pH adaptation of “Ca. Nitrosotalea devanaterra.” In contrast to neutrophilic AOA, “Ca. Nitrosotalea devanaterra” harbors a homologue of α-acetolactate decarboxylase (NDEV_1077), which is upregulated in bacteria during pH stress (52) and is predicted to catalyze the conversion of α-acetolactate to acetoin (Fig. 1 and 2). This is unexpected, as acetoin production is energetically costly for an autotrophic organism and α-acetolactate decarboxylase activity is generally greater under anaerobic conditions (53). How-

<footnote><sup>a</sup> Ratios of secondary to primary transporter genes in other archaeal acidophiles (41): *Thermoplasma acidophilum*, 10:1; *Picrophilus torridus*, 5.6:1; *Sulfolobus solfataricus*, 2.7:1.</footnote>
<figure><img src="image_4.png" /><figcaption>FIG 1 Predicted mechanisms of cytoplasmic pH regulation in "Ca. Nitrosotalea devanaterra" based on the presence of putative functional genes in the genome. (1) Cation influx and proton efflux. kdp, potassium transporting P-type ATPase gene; mgtA, putative magnesium-transporting P-type ATPase gene; mntA, NRAMP-type divalent cation transporter (two copies) gene, Na⁺/solute transporter, Na⁺/hydrogen exchanger (two copies). (2) Proton consumption by metabolism: acetolactate decarboxylase, carbonic anhydrase. (3) Reduced permeability of the cell wall/cell membrane: cell surface glycosylation, GDGT-4-dominated membrane.</figcaption></figure>

ever, α-acetolactate is produced in all AOA from pyruvate as part of valine, leucine, and isoleucine biosynthesis.

**Cell wall and membrane.** The membranes of some acidophiles, e.g., *Picrophilus oshimae*, require a low pH for stability (16), and the lipid composition and proportion of glycolipids change with varying pH (17). “*Ca. Nitrosotalea devanaterra*” contains two large cell surface modification gene clusters (31 and 34 ORFs) that are absent in neutrophilic AOA, with one cluster located next to genes encoding the main S-layer protein. These genes share only modest homology (30 to 61.4% amino acid identity) with bacterial and a limited number of archaeal genes. However, COG- and Pfam-based functional classification and their synteny with bacterial genomes support their involvement in cell surface glycosylation (see Table S8 and Fig. S2 in the supplemental material). The gene clusters of “*Ca. Nitrosotalea devanaterra*” do not belong to the recognized archaeal S-layer glycosylation pathways of *Euryarchaeota* and *Crenarchaeota* (54, 55). In bacteria, these homologues function in synthesis of the polysialic acid capsule (*neu/kps*), pseudoaminic acid glycosylation of flagella (*pse*), and outer spore coat glycosylation of Gram-positive bacteria (*sps*) (but not in bacterial peptidoglycan synthesis or modification) (56–58). There is also no evidence for the presence of peptidoglycan or pseudopeptidoglycan in “*Ca. Nitrosotalea devanaterra*,” which encodes two S-layer proteins (NDEV_0159 and NDEV_0294). It is plausible that the glycosylation genes act by adding nonulosonic acids to the S-layer of “*Ca. Nitrosotalea devanaterra*.”

The membrane lipid composition of “*Ca. Nitrosotalea devanaterra*” was investigated to determine whether it was distinct from that of neutrophilic AOA and therefore a potential adaptation to low pH. Thaumarchaea have glycerol dibiphytanyl glycerol tetraether lipids (GDGTs), with crenarchaeol typically the dominant GDGT in AOA (accounting for >30% of all lipids [24, 59–61]), but its proportion is much lower (11%) in “*Ca. Nitrosotalea devanaterra*.” Instead, GDGT-4 was the dominant membrane lipid in “*Ca. Nitrosotalea devanaterra*” (Fig. 3). The dominance of GDGT-4 over GDGT-0 to -3 in “*Ca. Nitrosotalea devanaterra*” compared to other AOA is in agreement with culture studies of acidophilic archaea which show an increase in the number of cyclopentane moieties in GDGTs with decreasing pH (62, 63). Furthermore, the increase in rings is thought to increase the packing density of the lipid membrane (64), which may reduce membrane permeability. The head groups of the “*Ca. Nitrosotalea devanaterra*” GDGTs are dominated by sugar moieties (glycolipids), in agreement with previous reports on acidophilic archaea which showed an increase in glycolipids versus phospholipids with decreasing pH (16). An increase in sugar units on the outside cell wall has been suggested to provide protection against proton entry (65). Hydroxylated GDGTs were also detected as cores, as in group 1.1a AOA (60), consistent with the 16S rRNA and *amoA* gene phylogenetic placement of *Nitrosotalea* (6).

**Central carbon metabolism and acidophily.** “*Ca. Nitrosotalea devanaterra*” possesses a carbonic anhydrase (NDEV_1082) (Fig. 1 and 2), which is present in many acidophiles (see Table S4 in the supplemental material) and predicted to catalyze rapid re-
<figure><img src="image_5.png" /><figcaption>FIG 2 Gene organization within potassium transporter (A) and pH homeostasis (B) gene clusters. Arrows indicate the relative lengths and directions of ORFs. Values inside the arrows indicate synteny E values. Values above the arrows indicate the gene name. (A) *kdpABC*, potassium-transporting ATPase gene; *kdpD*, signal transduction histidine kinase gene; PAC/PAS, PAS/PAC sensor histidine kinase gene; 1, gene for protein of unknown function; 2, signal transduction response regulator gene; 3, peptidase gene. (B) *mntH*, magnesium/divalent cation transporter gene; *cam*, carbonic anhydrase gene; *uspA*, universal stress protein gene; *mgtA*, magnesium transporting P-type ATPase gene; *alsD*, putative alpha-acetolactate decarboxylase gene; 1, phosphoribosyltransferase-like protein gene; 2, gene for conserved protein of unknown function; 3, Val tRNA gene; 4, adenylate/guanylate cyclase sensor protein gene; 5, gene for protein of unknown function; 6, response regulator receiver protein gene; 7, secreted copper domain-containing protein gene; 8, putative adenylate cyclase gene; 9, putative transposase gene.</figcaption></figure>

versible interconversion of HCO₃⁻ and CO₂. Although found also in *Nitrososphaera viennensis* and *"Candidatus Nitrososphaera evergladensis"*, carbonic anhydrase may have a dual function for carbon accumulation and cytoplasmic buffering in *Ca. Nitrosotalea devanaterra*. Depending on the direction of the reaction, carbonic anhydrase may prevent cytoplasmic acidification by breakdown of HCO₃⁻ or conversion of membrane-permeative CO₂ to HCO₃⁻ for carbon fixation. The reverse mechanism was proposed for carbon assimilation in *N. europaea* (66), with HCO₃⁻ transporter and carbonic anhydrase genes colocalized and HCO₃⁻ taken up, converted to CO₂, and fixed by the Calvin cycle.

**Ammonia oxidation at low pH.** The closed genome of "*Ca. Nitrosotalea devanaterra*" does not carry any genes participating in N metabolism that could explain its adaptation to apparently low substrate concentrations at acid pH (e.g., with urease). As AMO genes and ammonium transporters are found in all AOA (49), this suggests that the adaptations enabling "*Ca. Nitrosotalea devanaterra*" to grow with low ammonia concentrations either are based on differences from the active sites of AMO and transporters or reflect features shared by all AOA.

Since AOB cannot oxidize ammonia at low pH (8, 83), the substrate acquisition systems of "*Ca. Nitrosotalea devanaterra*," AOA, and AOB were compared. The three histidine residues (His³³, His¹³⁷, and His¹³⁹ in *M. capsulatus* Bath [37]) of AmoB/PmoB coordinating a periplasmic dicopper center are conserved in methanotrophs, AOB, and AOA (Fig. 4) but not in the acidophilic methane-oxidizing *Verrucomicrobia* (36, 37). In contrast, the variable metal binding site of AmoC is present in *Verrucomicrobia*, AOB, and AOA (38) (Fig. 4). *In silico* protein topology prediction for AmoB and AmoC favors an extracellular (outward-

facing) location of the active site in "*Ca. Nitrosotalea devanaterra*" (see Fig. S3 and S4 in the supplemental material). The active site of AMO is conserved between "*Ca. Nitrosotalea devanaterra*," AOA, and AOB and thus cannot explain the differences between these organisms.

Ammonium transporters of AOA and AOB are distantly related and belong to two functionally and phylogenetically distinct families: Amt and Rh, respectively (67). The substrate preferences (NH₃ versus NH₄⁺) often differ between Amt and Rh types (68, 69). Rh transporters are thought to perform bidirectional diffusion equilibrating NH₃ (70), and Amt transporters are energy dependent (71) (Fig. 5C). Furthermore, Amt transporters function better or equally well at acidic pH, whereas Rh transporters favor neutral/high pH (68, 72). Although archaeal ammonium transporters were extensively reviewed recently (67), the transported species (NH₃ versus NH₄⁺) of Amt/Rh transporters remains uncertain, and no hypothesis has been put forward for the preferred substrate.

Homology and conservation of coordinating residues with other transporters indicate that NH₄⁺ is bound to the Amt of "*Ca. Nitrosotalea devanaterra*" (Fig. 5; see Table S9 in the supplemental material). The cation binding site (required for NH₄⁺, but not NH₃, transport) is conserved in every Amt of all AOA and is absent in betaproteobacterial AOB (Fig. 5B) (73), as confirmed by the crystal structure of the *N. europaea* Rh protein (74). In addition, ammonium transporters may play a more critical role in AOA than in AOB, as all AOA genomes contain ≥2 Amt transporters, while five of 10 publicly available AOB genomes (*N. eutropha C71*, *Nitrosomas cryotolerans* ATCC 49181, *Nitrosococcus watsonii C113*, *Nitrosococcus oceani* AFC27, and *N. oceani* ATC19707) lack ammonia transporters. Unlike neutrophilic AOA, "*Ca. Nitrosotalea*
<figure><img src="image_6.png" /></figure>

FIG 3 Membrane lipid composition of “Ca. Nitrosotalea devanaterra.” Base peak chromatograms of core glycerol dibiphytanyl glycerol tetraether lipids (GDGTs) (A) and GDGT-based intact polar lipids (B) are shown. Cren' refers to a regioisomer of crenarchaeol. Percentages below “4+Crenarchaeol” indicate the relative contributions of GDGT-4 and crenarchaeol to the peak, respectively. Peaks belong to various core GDGTs from 0 to 4, including GDGTs with hydroxy cores.

lea devanaterra" has two putative low-affinity Amts (Amt1a and Amt1b) and one high-affinity Amt (Amt2) with a conserved signal peptide (75). The conservation of the cation binding site suggests that the mechanism enabling “Ca. Nitrosotalea devanaterra” to acquire ammonium is shared with all AOA but not AOB.

P<sub>II</sub> homologues are one of the most widespread signal trans-
duction proteins in prokaryotes and have a vital and exclusive role
in regulating N metabolism (76). P<sub>II</sub> homologues are found in the
vicinity of low-affinity Amt genes in neutrophilic group 1.1a AOA
(Fig. 4C) (67). In addition, a P<sub>II</sub> homologue is also located near the
high-affinity Amt in “Ca. Nitrosotalea devanaterra” (Fig. 4C).
Some AOB, e.g., *N. europaea*, lack P<sub>II</sub> genes despite the presence of
the Rh transporter (66).

**Transcriptional activity of genes potentially involved in pH homeostasis, ammonia acquisition, and ammonia oxidation.**
Several genes potentially associated with the acidophilic life-
style of “*Ca. Nitrosotalea devanaterra*” were chosen for further
transcriptional analysis: the genes for potassium-transporting
ATPase B chain (*kdpB* [NDEV_1460]), carbonic anhydrase (*cam*
NDEV_1082), α-acetolactate decarboxylase (*alsD* [NDEV_1077]),
or ammonium acquisition and oxidation (*amt1a* [NDEV_1347],
*amt1b* [NDEV_1108], *amt2* [NDEV_1784], and *amoA* [NDEV_1777]). All ammonium transporter genes, *amoA*, and
putative pH homeostasis genes, including those absent in neu-

tral-pH AOA, were transcriptionally active. The transcript/gene
ratios were highest at the optimal growth pH of 5 (see Fig. S5 and
Table S10 in the supplemental material), and transcription rates
were generally associated with nitrification activity rather than
being in response to pH change. Nevertheless, transcription of
ammonium transporter genes, *amoA*, and homeostasis genes at
pH 5 highlights their necessity for growth and metabolism in *Ca.
Nitrosotalea devanaterra.* Transcription of the putative high-af-
finity Amt2 gene was >10-fold higher than that of the low-affinity
Amt1a and Amt1b genes in *Ca. Nitrosotalea devanaterra*, as
reported for nonstarved neutrophilic AOA in culture and meta-
transcriptomic studies (75, 77).

**Obligate acidophily.** Representatives of the *Nitrosotalea* lineage are found overwhelmingly in acidic, rather than neutral, soils (4), and characterized *Nitrosotalea* isolates cannot grow at a pH of >6.1 (7), but the genome analysis provided no explanation for this inability. Ammonia may be toxic at neutral pH due to its higher concentration and its ability to permeate membranes. Neutral pH could also disable secondary transporters coupled to membrane pH gradients and change the bioavailability of many metals compared to that at acid pH (78).

**How does acidophilic ammonia oxidation occur in "Ca. Nitrosotalea devanaterra"?** Several lines of evidence favor an extracellular (outward-facing) orientation of the AMO active site in
<figure><img src="image_7.png" /><figcaption>FIG 4 Ammonia oxidation machinery of “Ca. Nitrosotalea devanaterra.” (A and B) Conservation of the periplasmic active site in PmoB/AmoB (A) and PmoC/AmoC (B) sequences. (C) AMO gene cluster, Amt transporter, and P~II~ gene arrangement. Longer alignments and protein topology predictions for AmoB and AmoC are in Fig. S3 and S4 in the supplemental material, respectively. Illustrated ORFs are drawn to scale, except for genes that are not involved in ammonia metabolism (truncated, white). The asterisks indicate that the “Candidatus Nitrosoarchaeum limnia” SFB1 and “Candidatus Nitrosotenuis uzonensis” genomes are not closed and distance cannot be estimated.</figcaption></figure>

"Ca. Nitrosotalea devanaterra." First, the active site of the pMMO of *M. capsulatus* Bath is periplasmic (36), and there is a high degree of conservation across methanotrophs, AOB, and AOA (38). Second, in silico topology prediction for AmoB and AmoC of "Ca. Nitrosotalea devanaterra" supports the outward orientation of the AMO active site. Third, AMO activity generates reactive intermediates, and it is advantageous to exclude them from the cytoplasm by performing the ammonia oxidation in the pseudoperiplasmic space (79).

Acidophilic growth could be facilitated by high substrate affinity: the whole cells of the marine AOA “Candidatus Nitrosopumilus maritimus” have an apparent half-saturation constant, $K_m$(NH$_3$+NH$_4^+$), of 133 nM at pH 7.5 (2, 80), which is equivalent to ~2 nM NH$_3$ at pH 7.5. At pH 4, 133 nM (NH$_3$+NH$_4^+$) is equivalent to ~0.75 pM, and ammonia oxidation at pH 4 would therefore require a $K_m$(NH$_3$) for AMO in the picomolar range or a mechanism for generating high localized concentrations of ammonia. Import of $NH_4^+$ followed by diffusion of membrane-permeative NH$_3$ into the pseudoperiplasm might generate local ammonia concentrations sufficiently high for oxidation (Fig. 6). $NH_4^+$, rather than NH$_3$, is abundant at low pH, but oxidation of $NH_4^+$ by the AMO would necessitate substantial changes to the enzyme, which would likely be evident in AMO sequence comparisons.

The regulation of pH consumes energy and may explain the lower growth yields of *Nitrosotalea* isolates (7) compared to neutral AOA (80). Unlike extreme acidophiles, "Ca. Nitrosotalea devanaterra" may not have a reverse membrane potential under all conditions. Assuming the same proton motive force as for *N. europaea* (~140 mV) (81) and a protonic pH of 7, membrane potential would be inside positive only at pHs of <4.6 at 25°C (Fig. 6).

The ammonia oxidation machineries of “Ca. Nitrosotalea devanaterra” and other AOA are similar, but most cultivated AOA cannot grow at acidic pH. This may be explained by the absence of the cytoplasmic pH regulation mechanisms found in “Ca. Nitrosotalea devanaterra,” e.g., ion transport and cell envelope modification, from other AOA. Conversely, even if AOB possessed genes for acid pH homeostasis, starvation would result from a lack of ammonia (without ureolytic activity). If AOA are able to use ammonium for transport as our results suggest, then AOA rather than AOB were in a better predisposition to evolve an acidophilic phenotype.

**Conclusions.** Despite its ecological significance, the mechanism of acid-tolerant ammonia oxidation is a long-standing paradox. This investigation has examined, for the first time, how an ammonia-oxidizing organism may perform this process at low pH by examining an archaeal genome with specific emphasis on (i) N metabolism, (ii) C metabolism, and (iii) pH homeostasis mechanisms. We predict the following. (i) Previously proposed mechanisms of acid tolerance in AOB (e.g., urease, biofilm, or aggregate formation) are not required for the growth of the obligate acidophile "Ca. Nitrosotalea devanaterra." Further, predictions suggest that the AMO of "Ca. Nitrosotalea devanaterra" oxidizes NH$_3$ and has an extracellular (outward-facing) active site and that Amt transporters of "Ca. Nitrosotalea devanaterra" bind $NH_4^+$, in contrast to AOB, which transport NH$_3$. (ii) $HCO_3^-$, required for carbon fixation, may be supplied by carbonic anhydrase, which is actively transcribed in "Ca. Nitrosotalea devanaterra." (iii) "Ca. Nitrosotalea devanaterra" shares aspects of pH homeostasis mechanisms with other moderate acidophiles, including having a high number of cation transporters, cytoplasmic proton-scavenging mechanisms, and altered cell membrane and surface compositions compared to neutrophilic AOA. In contrast to the case for
<figure><img src="image_8.png" /><figcaption>FIG 5 Ammonium transporters of "Ca. Nitrosotalea devanaterra." (A) Maximum-likelihood phylogeny of Amt/Rh transporters. The transported substrate (ammonia versus ammonium) has been validated in organisms highlighted with red and blue circles. Red, ammonium (NH₄⁺); blue, ammonia (NH₃). Circles at nodes highlight the level of bootstrap support, and the scale bar represents 0.5 change per amino acid position. (B) Conservation of the cation binding site in Amt/Rh sequences. (C) Transport mechanisms of Amt/Rh transporters previously reported in the literature (see Table S9 in the supplemental material).</figcaption></figure>

extreme acidophiles, secondary transporters are not overrepresented in "Ca. Nitrosotalea devanaterra." We postulate that unlike AOB, all AOA, including "Ca. Nitrosotalea devanaterra," have the genetic potential to oxidize ammonia under acidic conditions and

that acidophilic adaptation of "Ca. Nitrosotalea devanaterra" is achieved by cytoplasmic pH regulation mechanisms not shared by neutrophilic AOA.

The current study (as with all genome studies) is limited by
<figure><img src="image_9.png" /><figcaption>FIG 6 Conceptual model of ammonia acquisition and oxidation mechanisms in “Ca. Nitrosotalea devanaterra.” The model is based on references 70 and 81 and work in this study.</figcaption></figure>

existing sequence databases but proposes hypothetical mechanisms for growth of acidophilic archaea that can provide the basis for future experimental testing. Although the presence and activity of previously uncharacterized genes were associated with the acidophilic growth of “Ca. Nitrosotalea devanaterra,” one caveat of this approach is that “Ca. Nitrosotalea devanaterra” cannot grow at neutral pH. It is therefore impossible to unequivocally confirm that these genes are a prerequisite for acidophily. This question could be further explored by heterologous expression of candidate acidophily genes in a neutrophilic organism (preferably a neutrophilic AOA) to examine whether they confer an acidophilic phenotype on the host. An alternative approach would be site-directed mutagenesis, although unfortunately neither of these techniques has yet been developed for AOA. Further, acetoin production could be determined, transcription of further genes (e.g., Na⁺/solute symporter or Na⁺/H⁺ exchanger genes) studied, and membrane lipid composition determined over a range of pH. Other important future research directions include genome sequencing, lipid profiling, and S-layer glycome profiling of other strains of *Nitrosotalea* (and potentially acid-tolerant urease-positive AOB).

## ACKNOWLEDGMENTS

We thank W. I. C. Rijpstra (NIOZ) for analytical support of the lipid analysis. We thank Marvin Djukic and Heiko Nacke for useful discussion on bioinformatics.

## FUNDING INFORMATION

This work, including the efforts of Luis Sayavedra-Soto, was funded by National Science Foundation (NSF) (0541797). This work, including the efforts of Laura E. Lehtovirta-Morley, James I. Prosser, and Graeme W. Nicol, was funded by Natural Environment Research Council (NERC) (NE/1027835/1). Graeme W. Nicol is also supported by the AXA Research Fund.

## REFERENCES

1. Stahl DA, de la Torre JR. 2012. Physiology and diversity of ammonia-oxidizing archaea. Annu Rev Microbiol 66:83-101. http://dx.doi.org/10.1146/annurev-micro-092611-150128.

2. Prosser JI, Nicol GW. 2012. Archaeal and bacterial ammonia-oxidisers in soil: the quest for niche specialisation and differentiation. Trends Microbiol 20:523-531. http://dx.doi.org/10.1016/j.tim.2012.08.001.

3. Verhamme DT, Prosser JI, Nicol GW. 2011. Ammonia concentration determines differential growth of ammonia-oxidising archaea and bacteria in soil microcosms. ISME J 5:1067-1071. http://dx.doi.org/10.1038/ismej.2010.191.

4. Gubry-Rangin C, Hai B, Quince C, Engel M, Thomson BC, James P, Schloter M, Griffiths RI, Prosser JI, Nicol GW. 2011. Niche specialization of terrestrial archaeal ammonia oxidizers. Proc Natl Acad Sci U S A 108:21206-21211. http://dx.doi.org/10.1073/pnas.1109000108.

5. Booth MS, Stark JM, Rastetter E. 2005. Controls on nitrogen cycling in terrestrial ecosystems: a synthetic analysis of literature data. Ecol Monogr 75:139-157. http://dx.doi.org/10.1890/04-0988.

6. Lehtovirta-Morley LE, Stoecker K, Vilcinskas A, Prosser JI, Nicol GW. 2011. Cultivation of an obligately acidophilic ammonia oxidizer from a nitrifying acid soil. Proc Natl Acad Sci U S A 108:15892-15897. http://dx.doi.org/10.1073/pnas.1107196108.

7. Lehtovirta-Morley LE, Ge C, Ross J, Yao H, Nicol GW, Prosser JI. 2014. Characterisation of terrestrial acidophilic archaeal ammonia oxidisers and their inhibition and stimulation by organic compounds. FEMS Microbiol Ecol 89:542-552. http://dx.doi.org/10.1111/1574-6941.12353.

8. Suzuki I, Dular U, Kwok SC. 1974. Ammonia and ammonium ion as substrate for oxidation by *Nitrosomonas* cells and extracts. J Bacteriol 120: 556-558.

9. Burton SAQ, Prosser JI. 2001. Autotrophic ammonia oxidation at low pH through urea hydrolysis. Appl Environ Microbiol 67:2952-2957. http://dx.doi.org/10.1128/AEM.67.7.2952-2957.2001.

10. Allison SM, Prosser JI. 1993. Ammonia oxidation at low pH by attached populations of nitrifying bacteria. Soil Biol Biochem 25:935-941. http://dx.doi.org/10.1016/0038-0717(93)90096-T.

11. de Boer W, Gunnewiek PJ, Veenhuis M, Bock E, Laanbroek HJ. 1991. Nitrification at low pH by aggregated chemolithotrophic bacteria. Appl Environ Microbiol 57:3600-3604.

12. Arp DJ, Sayavedra-Soto LA, Hommes NG. 2002. Molecular biology and biochemistry of ammonia oxidation by *Nitrosomonas europaea*. Arch Microbiol 178:250-255. http://dx.doi.org/10.1007/s00203-002-0452-0.

13. Berg IA, Kockelkorn D, Ramos-Vera WH, Say RF, Zarzycki J, Hügler M, Alber BE, Fuchs G. 2010. Autotrophic carbon fixation in archaea. Nat Rev Microbiol 8:447-460. http://dx.doi.org/10.1038/nrmicro2365.

14. Slonczewski JL, Fujisawa M, Dopson M, Krulwich TA. 2009. Cytoplasmic pH measurement and homeostasis in bacteria and archaea. Adv Microb Physiol 55:1-79. http://dx.doi.org/10.1016/S0065-2911(09)05501-5.

15. Baker-Austin C, Dopson M. 2007. Life in acid: pH homeostasis in acidophiles. Trends Microbiol 15:165-171. http://dx.doi.org/10.1016/j.tim.2007.02.005.

16. Shimada H, Nemoto N, Shida Y, Oshima T, Yamagishi A. 2008. Effects of pH and temperature on the composition of polar lipids in *Thermoplasma acidophilum* HO-62. J Bacteriol 190:5404-5411. http://dx.doi.org/10.1128/JB.00415-08.

17. van de Vossenberg JLCM, Driessen AJM, Zillig W, Konings WN. 1998. Bioenergetics and cytoplasmic membrane stability of the extremely acidophilic, thermophilic archaeon *Picrophilus oshimae*. Extremophiles 2:67-74. http://dx.doi.org/10.1007/s007920050044.

18. Michels M, Bakker EP. 1985. Generation of a large protonophoresensitive proton motive force and pH difference in the acidophilic bacteria *Thermoplasma acidophilum* and *Bacillus acidocaldarius*. J Bacteriol 161: 231-237.

19. Sachs G, Weeks DL, Wen Y, Marcus EA, Scott DR. 2005. Acid acclimation by *Helicobacter pylori*. Physiology 20:429-438. http://dx.doi.org/10.1152/physiol.00032.2005.
20. Richard H, Foster JW. 2004. *Escherichia coli* glutamate- and arginine-dependent acid resistance systems increase internal pH and reverse trans-membrane potential. J Bacteriol **186**:6032–6041. http://dx.doi.org/10.1128/JB.186.18.6032-6041.2004.

21. Alexander B, Leach S, Ingledew WJ. 1987. The relationship between chemiosmotic parameters and sensitivity to anions and organic acids in the acidophile *Thiobacillus ferrooxidans*. J Gen Microbiol **133**:1171–1179.

22. Fütterer O, Angelov A, Liesegang B, Gottschalk G, Schleper C, Shepers B, Dock C, Antranikian G, Liebl W. 2004. Genome sequence of *Picrophilus torridus* and its implications for life around pH 0. Proc Natl Acad Sci U S A **101**:9091–9096. http://dx.doi.org/10.1073/pnas.0401356101.

23. Bramwell PA, Barallon RV, Rogers HJ, Bailey MJ. 1995. Extraction of microbial DNA from the rhizoplane, p 36–55. *In* Akkermans ADL, van Elsas JD, De Bruijn FJ (ed), Molecular microbial ecology manual. Kluwer Academic Publishers, Dordrecht, The Netherlands.

24. Myers EW, Sutton GG, Delcher AL, Dew IM, Fasulo DP, Flanigan MJ, Kravitz SA, Mobarry CM, Reinert KH, Remington KA, Anson EL, Bolanos RA, Chou HH, Jordan CM, Halpern AL, Lonardi S, Beasley EM, Brandon RC, Chen L, Dunn PJ, Lai Z, Liang Y, Nusskern DR, Zhan M, Zhang Q, Zheng X, Rubin GM, Adams MD, Venter JC. 2000. A whole-genome assembly of Drosophila. Science **287**:2196–2204. http://dx.doi.org/10.1126/science.287.5461.2196.

25. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahni N, Tesler G, Alekseyev MA, Pevzner PA. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol **19**:455–477. http://dx.doi.org/10.1089/cmb.2012.0021.

26. Darling ACE, Mau B, Blattner FR, Perna NT. 2004. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res **14**:1394–1403. http://dx.doi.org/10.1101/gr.2289704.

27. Vallenet D, Engelen S, Mornico D, Cruveiller S, Fleury L, Lajus A, Rouy Z, Roche D, Salvignol G, Scarpelli C, Médigue C. 2009. MicroScope: a platform for microbial genome annotation and comparative genomics. Database **2009**:bap021.

28. Schattner P, Brooks AN, Lowe TM. 2005. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res **33**:686–689.

29. Laslett D, Canback B. 2004. ARAGORN, a program for the detection of transfer RNA and transfer-messenger RNA genes in nucleotide sequences. Nucleic Acids Res **32**:11–16. http://dx.doi.org/10.1093/nar/gkh152.

30. Ren Q, Chen K, Paulsen IT. 2007. TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res **35**:D274–D279. http://dx.doi.org/10.1093/nar/gkl925.

31. Richter M, Rosselló-Móra R. 2009. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A **106**: 19126–19131. http://dx.doi.org/10.1073/pnas.0906412106.

32. Hammer Ø Harper DAT, Ryan PD. 2001. PAST: Paleontological statistics software package for education and data analysis. Palaeontol Electron **4**:9.

33. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. 2009. BLAST+: architecture and applications. BMC Bioinformatics **10**:421. http://dx.doi.org/10.1186/1471-2105-10-421.

34. Schouten S, Huguet C, Hopmans EC, Sinninghe Damsté JS. 2007. Improved analytical methodology of the TEX86 paleochemistry by high performance liquid chromatography/atmospheric pressure chemical ionization-mass spectrometry. Anal Chem **79**:2940–2944. http://dx.doi.org/10.1021/ac062339v.

35. Hall TA. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser (Oxf) **41**:95–98.

36. Lieberman RL, Rosenzweig AC. 2005. Crystal structure of a membrane-bound metalloenzyme that catalyses the biological oxidation of methane. Nature **434**:177–182. http://dx.doi.org/10.1038/nature03311.

37. Balasubramanian R, Smith SM, Rawat S, Yatsunyk LA, Stemmler TL, Rosenzweig AC. 2010. Oxidation of methane by a biological dicopper centre. Nature **465**:115–119. http://dx.doi.org/10.1038/nature08992.

38. Liew EF, Tong D, Coleman NV, Holmes AJ. 2014. Mutagenesis of the hydrocarbon monooxygenase indicates a metal centre in subunit-C, and not subunit-B, is essential for copper-containing membrane monooxygenase activity. Microbiology **160**:1267–1277. http://dx.doi.org/10.1099/mic.0.078584-0.

39. Lawton TJ, Ham J, Sun T, Rosenzweig AC. 2014. Structural conservation of the B subunit in the ammonia monooxygenase/particulate methane monooxygenase superfamily. Proteins **82**:2263–2267. http://dx.doi.org/10.1002/prot.24535.

40. Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. 2001. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol **305**:567–580. http://dx.doi.org/10.1006/jmb.2000.4315.

41. Hofmann K, Stoffel W. 1993. TMbase—a database of membrane spanning proteins segments. Biol Chem Hoppe-Seyler **374**:166.

42. von Heijne G. 1992. Membrane protein structure prediction, hydrophobicity analysis and the positive-inside rule. J Mol Biol **225**:487–494. http://dx.doi.org/10.1016/0022-2836(92)90934-C.

43. Bendtsen JD, Nielsen H, von Heijne G, Brunak S. 2004. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol **340**:783–795. http://dx.doi.org/10.1016/j.jmb.2004.05.028.

44. Käll L, Krogh A, Sonnhammer ELL. 2004. A combined transmembrane topology and signal peptide prediction method. J Mol Biol **338**:1027–1036. http://dx.doi.org/10.1016/j.jmb.2004.03.016.

45. Guindon S, Gascuel O. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol **52**:696–704. http://dx.doi.org/10.1080/10635150390235520.

46. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol **30**: 2725–2729. http://dx.doi.org/10.1093/molbev/mst197.

47. Nicol GW, Leininger S, Schleper C, Prosser JI. 2008. The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. Environ Microbiol **10**:2966–2978. http://dx.doi.org/10.1111/j.1462-2920.2008.01701.x.

48. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG. 2012. Primer3—new capabilities and interfaces. Nucleic Acids Res **40**:e115. http://dx.doi.org/10.1093/nar/gks596.

49. Spang A, Poehlein A, Offre P, Zumbrägl S, Haider S, Rychlik N, Nowka B, Schmeisser C, Lebedeva EV, Rattei T, Böhm C, Schmid M, Galushko A, Hatzenpichler R, Weinmaier T, Daniel R, Schleper C, Speick E, Streit W, Wagner M. 2012. The genome of the ammonia-oxidizing *Candidatus* Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations. Environ Microbiol **14**:3122–3145. http://dx.doi.org/10.1111/j.1462-2920.2012.02893.x.

50. Walker CB, de la Torre JR, Klotz MG, Urakawa H, Pinel N, Arp DJ, Brochier-Armanet C, Chain PS, Chan PP, Gollabgir A, Hemp J, Hügler M, Karr EA, Könneke M, Shin M, Lawton TJ, Lowe T, Martens-Habbena W, Sayavedra-Soto LA, Lang D, Sievert SM, Rosenzweig AC, Manning G, Stahl DA. 2010. Nitrosoopumilus maritimus genome reveals unique mechanisms for nitrification and autotrophy in globally distributed marine crenarchaea. Proc Natl Acad Sci U S A **107**:8818–8823. http://dx.doi.org/10.1073/pnas.0913533107.

51. Lewalter K, Muller V. 2006. Bioenergetics of archaea: ancient energy conserving mechanisms developed in the early history of life. Biochim Biophys Acta Bioenerg **1757**:437–445. http://dx.doi.org/10.1016/j.bbabio.2006.04.027.

52. Wilks JC, Kitko RD, Cleeton SH, Lee GE, Ugwu CS, Jones BD, Bon-Durant SS, Slonczewski JL. 2009. Acid and base stress and transcriptomic responses in *Bacillus subtilis*. Appl Environ Microbiol **75**:981–990. http://dx.doi.org/10.1128/AEM.01652-08.

53. Monnet C, Aymes F, Corrieu G. 2000. Diacetyl and α-acetolactate over-production by *Lactococcus lactis* subsp. *lactis* biovar diacetylactis mutants that are deficient in α-acetolactate decarboxylase and have a low lactate dehydrogenase activity. Appl Environ Microbiol **66**:5518–5520. http://dx.doi.org/10.1128/AEM.66.12.5518-5520.2000.

54. Albers SV, Meyer BH. 2011. The archaeal cell envelope. Nat Rev Microbiol **9**:414–426. http://dx.doi.org/10.1038/nrmicro2576.

55. Calo D, Kaminski L, Eichler J. 2010. Protein glycosylation in archaea: sweet and extreme. Glycobiology **20**:1065–1076. http://dx.doi.org/10.1093/glycob/cwq055.

56. Bliss JM, Silver RP. 1996. Coating the surface: a model for expression of capsular polysialic acid in *Escherichia coli* K1. Mol Microbiol **21**:221–231. http://dx.doi.org/10.1046/j.1365-2958.1996.6461357.x.

57. McNally DJ, Hui JPM, Aubry AJ, Mui KKK, Guerry P, Brisson JR, Logan SM, Soo EC. 2006. Functional characterization of the flagellar glycosylation locus in *Campylobacter jejuni* 81-176 using a focused metabolomics approach. J Biol Chem **281**:18489–18498. http://dx.doi.org/10.1074/jbc.M603777200.
58. Driks A. 1999. *Bacillus subtilis* spore coat. Microbiol Mol Biol Rev 63:1-20.

59. Elling FJ, Könneke M, Lipp JS, Becker KW, Gagen EJ, Hinrichsn KU. 2014. Effect of growth phase on the membrane lipid composition of the thaumarchaeon *Nitrosopumilus maritimus* and their implication for ar- chaeal lipid distribution in the marine environment. Geochim Cos- mochim Acta 141:579-597. http://dx.doi.org/10.1016/j.gca.2014.07.005.

60. Pitcher A, Hopmans EC, Mosier AC, Park SJ, Rhee SK, Francis CA, Schouten S, Sinninghe Damsté JS. 2011. Core and intact polar glycerol dibiphytanyl glycerol tetraether lipids of ammonia-oxidizing archaea en- riched from marine and estuarine sediments. Appl Environ Microbiol 77:3468-3477. http://dx.doi.org/10.1128/AEM.02758-10.

61. Sinninghe Damsté JS, Hopmans EC, Schouten S, Van Duin ACT, Geenevasen JAJ. 2002. Crenarchaeol: the characteristic core glycerol dibi- phytanyl glycerol tetraether membrane lipid of cosmopolitan pelagic cre- narchaeota. J Lip Res 43:1641-1651. http://dx.doi.org/10.1194/jlr .M200148-JLR200.

62. Boyd ES, Hamilton TL, Wang J, He L, Zhang CL. 2013. The role of tetraether lipid composition in the adaptation of thermophilic archaea to acidity. Front Microbiol 4:62. http://dx.doi.org/10.3389/fmicb.2013 .00062.

63. Boyd E, Pearson A, Pi Y, Li WJ, Zhang Y, He L, Zhang CL, Geesey G. 2011. Temperature and pH controls on glycerol dibiphytanyl glycerol tetraether lipid composition in the hyperthermophilic crenarchaeon *Acidi- lobus sulfurreducens*. Extremophiles 15:59-65. http://dx.doi.org/10.1007 /s00792-010-0339-y.

64. Gabriel JL, Chong PL. 2000. Molecular modelling of archaebacterial bipolar tetraether lipid membranes. Chem Phys Lipids 105:193-200. http: //dx.doi.org/10.1016/S0009-3084(00)00126-2.

65. Wang X, Lv B, Cai G, Fu L, Wu Y, Wang X, Ren B, Ma H. 2012. A proton shelter inspired by the sugar coating of acidophilic archaea. Sci Rep 2:892. http://dx.doi.org/10.1038/srep00892.

66. Chain P, Lamerdin J, Larimer F, Regala W, Lao V, Land M, Hauser L, Hooper A, Klotz M, Norton J, Sayavedra-Soto L, Arciero D, Hommes N, Whittaker M, Arp D. 2003. Complete genome sequence of the am- monia-oxidizing bacterium and obligate chemolithoautotroph *Nitro- somonas europaea*. J Bacteriol 185:2759-2773. http://dx.doi.org/10.1128 /JB.185.9.2759-2773.2003.

67. Offre P, Kerou M, Spang A, Schleper C. 2014. Variability of the trans- porter gene complement in ammonia-oxidizing archaea. Trends Micro- biol 22:665-675. http://dx.doi.org/10.1016/j.tim.2014.07.007.

68. Mayer M, Schaaf G, Mouro I, Lopez C, Colin Y, Neumann P, Cartron JP, Ludewig U. 2006. Different transport mechanisms in plant and hu- man AMT/Rh-type ammonium transporters. J Gen Physiol 127:133-144. http://dx.doi.org/10.1085/jgp.200509369.

69. Ludewig U. 2004. Electroneutral ammonium transport by basolateral rhesus B glycoprotein. J Physiol 559:751-759. http://dx.doi.org/10.1113 /jphysiol.2004.067728.

70. Weidinger K, Neuhäuser B, Gilch S, Ludewig U, Meyer O, Schmidt I. 2007. Functional and physiological evidence for a rhesus-type ammonia transporter in *Nitrosomonas europaea*. FEMS Microbiol Lett 273:260-267. http://dx.doi.org/10.1111/j.1574-6968.2007.00805.x.

71. Ullmann RT, Andrade SL, Ullmann GM. 2012. Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations. J Phys Chem B 116:9690-9703. http://dx.doi.org /10.1021/jp305440f.

72. Sogaard R, Alsterfjord M, MacAulay N, Zeuthen T. 2009. Ammonium ion transport by the AMT/Rh homolog TaAMT1;1 is stimulated by acidic pH. Pflugers Arch Eur J Physiol 458:733-743. http://dx.doi.org/10.1007 /s00424-009-0665-z.

73. Khademi S, O'Connell J, Remis J, Robles-Colmenares Y, Miercke LJ, Stroud RM. 2004. Mechanism of ammonia transport by Amt/MEP/Rh: structure of AmtB at 1.35 Å. Science 305:1587-1594. http://dx.doi.org/10 .1126/science.1101952.

74. Li X, Jayachandran S, Nguyen HT, Chan MK. 2007. Structure of the *Nitrosomonas europaea* Rh protein. Proc Natl Acad Sci U S A 104:19279- 19284. http://dx.doi.org/10.1073/pnas.0709710104.

75. Nakagawa T, Stahl DA. 2013. Transcriptional response of the archaeal ammonia oxidizer *Nitrosopumilus maritimus* to low and environmentally relevant ammonia concentrations. Appl Environ Microbiol 79:6911- 6916. http://dx.doi.org/10.1128/AEM.02028-13.

76. Arcondeguy T, Jack R, Merrick M. 2001. P<sub>II</sub> signal transduction proteins, pivotal players in microbial nitrogen control. Microbiol Mol Biol R 65: 80-105. http://dx.doi.org/10.1128/MMBR.65.1.80-105.2001.

77. Hollibaugh JT, Gifford S, Sharma S, Bano N, Moran MA. 2011. Meta- transcriptomic analysis of ammonia-oxidizing organisms in an estuarine bacterioplankton assemblage. The ISME J 5:866-878. http://dx.doi.org /10.1038/ismej.2010.172.

78. Violante A, Cozzolino V, Perelomov L, Caporale AG, Pigna M. 2010. Mobility and bioavailability of heavy metals and metalloids in soil envi- ronments. J Soil Sci Plant Nutr 10:268-292.

79. Vajrala N, Martens-Habbena W, Sayavedra-Soto LA, Schauer A, Bot- tomley PJ, Stahl DA, Arp DJ. 2013. Hydroxylamine as an intermediate in ammonia oxidation by globally abundant marine archaea. Proc Natl Acad Sci U S A 110:1006-1011. http://dx.doi.org/10.1073/pnas.1214272110.

80. Martens-Habbena W, Berube PM, Urakawa H, de la Torre JR, Stahl DA. 2009. Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461:976-979. http://dx.doi.org /10.1038/nature08465.

81. Kumar S, Nicholas DJ. 1983. Proton electrochemical gradients in washed cells of *Nitrosomonas europaea* and *Nitrobacter agilis*. J Bacteriol 154:65-71.

82. Sinninghe Damsté JS, Rijpstra WIC, Hopmans EC, Jung MY, Kim JG, Rhee S-K, Stieglmeier M, Schleper C. 2012. Intact polar and core glycerol dibiphytanyl glycerol tetraether lipids of group I.1a and I.1b Thaumar- chaeota in soil. Appl Environ Microbiol 78:6866-6874. http://dx.doi.org /10.1128/AEM.01681-12.

83. Frijlink MJ, Abee T, Laanbroek HJ, de Boer W, Konings WN. 1992. The bioenergetics of ammonia and hydroxylamine oxidation in *Nitrosomonas europaea* in acid and alkaline pH. Arch Microbiol 157:194-199. http://dx .doi.org/10.1007/BF00245290.