The draft genome of the hyperthermophilic archaeon Pyrodictium delaneyi strain hulk, an iron and nitrate reducer, reveals the capacity for sulfate reduction

Pyrodictium delaneyi strain Hulk is a newly sequenced strain isolated from chimney samples collected from the Hulk sulfide mound on the main Endeavour Segment of the Juan de Fuca Ridge (47.9501 latitude, −129.0970 longitude, depth 2200 m) in the Northeast Pacific Ocean. The draft genome of strain Hulk shared 99.77% similarity with the complete genome of the type strain Su06T, which shares with strain Hulk the ability to reduce iron and nitrate for respiration. The annotation of the genome of strain Hulk identified genes for the reduction of several sulfur-containing electron acceptors, an unsuspected respiratory capability in this species that was experimentally confirmed for strain Hulk. This makes P. delaneyi strain Hulk the first hyperthermophilic archaeon known to gain energy for growth by reduction of iron, nitrate, and sulfur-containing electron acceptors. Here we present the most notable features of the genome of P. delaneyi strain Hulk and identify genes encoding proteins critical to its respiratory versatility at high temperatures. The description presented here corresponds to a draft genome sequence containing 2,042,801 bp in 9 contigs, 2019 protein-coding genes, 53 RNA genes, and 1365 hypothetical genes.


Introduction
The unifying metabolic feature of the first five species described in the family Pyrodictiaceae, in the archaeal order Desulfurococcales, was for long their ability to respire sulfur-containing electron acceptors, mainly elemental sulfur (S 0 ), thiosulfate (S 2 O 3 2− ) and sulfite (SO 3 2− ) [1]. Pyrolobus fumarii 1A T reduces nitrate in addition to thiosulfate [2]. Yet the recently described Pyrodictium delaneyi Su06 T reduces nitrate and iron but cannot use sulfur or thiosulfate [3].
Here we report on the isolation and genome sequencing and annotation of a novel strain of P. delaneyi, designated strain Hulk, capable of hyperthermophilic growth with iron, nitrate, and several sulfur-containing compounds. This makes the novel strain Hulk the first hyperthermophilic archaeon known to respire iron, nitrate, and sulfurcontaining electron acceptors. Furthermore, strain Hulk is the first member of the Pyrodictiaceae family able to use formate as an electron donor and carbon source. In addition, it oxidized peptides, an ability only reported for the two obligate peptide organotrophs in the family, Pyrodictium abyssi AV2 T [4] and Hyperthermus butylicus PLM1-5 T [5]. The ability to oxidize formate and peptides is of special environmental significance, as these are abundant electron donors in hydrothermal marine vent systems [6][7][8][9]. Analysis of the draft genome of strain Hulk reveals numerous pathways and genes that allow this archaeon to couple autotrophic and heterotrophic growth with these many electron donors and acceptors.

Classification and features
A novel strain of P. delaneyi designated strain Hulk (Fig. 1) was isolated from a hot sediment sample collected from the Hulk hydrothermal vent located on the Main Endeavour segment of the Juan de Fuca Ridge (47.9501 latitude −129.0970 longitude) in the Northeast Pacific Ocean approximately 300 miles west of Seattle, Washington, at the depth of 2200 m. Enrichment cultures used marine enrichment media [10] modified with the addition NaCl (20 g l −1 ) and yeast extract (25 μg l −1 ), and supplemented with pebble-milled cellulose (0.24%, wt/vol) as the electron donor and poorly crystalline Fe(III) oxide (100 mmol l −1 ) as electron acceptor. The headspace of the tubes was pressurized with N 2 :CO 2 (80:20%, v/v, 101 kPa). Isolation was in the same medium solidified with GELRITE gellan gum (Sigma-Aldrich), as previously described [10]. All incubations during enrichments and isolation were at 100°C.
The draft genome of P. delaneyi strain Hulk contained a single 16S rRNA gene, whose sequence was 100% identical to the 16S rDNA sequence in the complete genome of P. delaneyi strain Su06 T [3]. Pair-wise Fig. 1 Phylogenetic tree constructed with the maximum likelihood algorithm comparing the 16S rRNA gene sequence from P. delaneyi, strain Hulk to the species type strain Su06 T and other hyperthermophilic archaea. GenBank accession numbers are listed in parentheses. Bootstrap values displayed at branch points were determined from 100 replicates [72,73] genome comparisons between the two genomes [11] revealed an average nucleotide identity (ANI) of 99.77%, above the cutoff (92%) established for bacterial species definition [12]. The full-length sequence of the 16S rRNA gene of strain Hulk (1496 bp) was used to construct a phylogenetic tree in reference to 16S rRNA gene sequences from other hyperthermophilic archaea (Fig. 1). The Hyperthermus butylicus was the nearest neighbor outside of the P. delaneyi species group (99% identity).

Genome sequencing information
Genome project history P. delaneyi strain Hulk was sequenced and annotated based on its phylogenetic position and its metabolic versatility. Genome comparisons [11] with the recently described species type strain Su06 T [3] and other members of the Pyrodictiaceae were performed to provide novel insights into this archaeal family and the metabolic potential of their members. The draft genome of P. delaneyi strain Hulk presented here is contained within 9 contigs with an average coverage of at least 380× (Table 2). This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession NCQP00000000 ( Table 2). The version described in this paper is version NCQP01000000. The genome project summary can be viewed at Genomes Online Database (Ga0169944) [8] (Table 2).

Growth conditions and genomic DNA preparation
Genomic DNA was extracted from exponentially-grown cells from formate-Fe(III) citrate cultures incubated at 90°C as previously described [13]. Cells were harvested , not directly observed for the living, isolated sample, but based on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from the Gene Ontology project [81] at room temperature by centrifugation for 20 min, washed with a marine wash buffer [14], and extracted using the MoBio Powersoil kit, except that an additional lysis buffer step was included, as reported elsewhere [14]. Quality of genomic DNA was verified by nanodrop and gel electrophoresis.

Genome sequencing and assembly
The assembly presented was generated from three paired-end Illumina libraries [15]. One of the libraries was sequenced at the Research Technology Support Facility center at Michigan State University and sequence data for the remaining libraries generated at Swift Biosciences. Genomic DNA was fragmented to an average size of 200 bp using a Covaris M220 (Covaris, Woburn, MA). Two genome libraries were generated using the Accel-NGS 2S DNA Library kit and Accel-NGS 1S DNA Library Kit (Swift Biosciences, Ann Arbor, MI). The third library was generated at the Michigan State University genomics core. These libraries were sequenced on an Illumina MiSeq (Illumina, San Diego, CA) using MiSeq Reagent kit v2. The read data from each of the three libraries were trimmed using fastq-mcf [16] and ConDeTri [17] and analyzed with the genome assembly program Velvet [17] to generate three independent draft de novo genome assemblies. A consensus . The reduction of DMSO supports increases in cell numbers compared to controls without the electron donor (c). All cultures were prepared at the optimum pH and salt concentration (7 and 2%, respectively) and incubated at the optimum temperature (90°C) assembly was then generated using Mix [18] with a minimum alignment length of 200 bp and a minimum contig length of 0 bp. Contigs with total length below 1000 bp were removed from the final assembly. This resulted in an assembly of 9 contigs with an average coverage depth exceeding 380×.

Genome annotation
The genome annotation used the Rapid Annotations via Subsystems Technology (RAST) server [19]. Coding regions were identified with RAST's GLIMMER tool [20]. This initial annotation of protein encoding regions in the genome was then manually refined for genes of interest. Selected genes were also analyzed with DELTA-BLAST and PSI-BLAST to identify conserved domains and homology to known proteins and to infer functions. Enzyme Commission numbers and Clusters of Orthologous Group categories were determined with a combination of DELTA-BLAST analysis of each annotated gene and the IMG-ER platform [21]. After COGs were identified the genome was again submitted for annotation via the Prokaryotic Genome Annotation Pipeline and verified the genes referenced here.
Putative c-type cytochromes were identified based on the presence of conserved heme-binding motifs (CXXCH), as previously described [13]. The presence of signal peptides or N-terminal helix membrane anchors in the proteins containing the conserved heme-binding motif was then assessed using PRED-TAT [22] and the TMHMM Server (v2.0) [23], respectively, to infer their cellular localization (exported or membrane-bound, respectively). DELTA-BLAST was then used to assess the homology of the putative c-type cytochrome proteins with known c-type cytochromes in the non-redundant NCBI database.

Genome properties
The draft genome sequence of P. delaneyi strain Hulk was assembled into 9 contigs (N50 557,338 bp, total length 2,042,801 bp) with a GC content of 53.88% (Table 3). The draft genome size is close to that reported for the complete genome of the type strain of P. delaneyi Su06 T (2,023,836 bp) [3]. Out of the total of 2089 genes identified in the genome sequence of strain Hulk, 53 were predicted to encode RNAs and 2019 proteins (Table 3). Seventeen pseudogenes were identified by the Prokaryotic Genomes Automatic Annotation Pipeline [24]. Furthermore, 67.73% of the predicted genes (2089) are represented by COG functional categories. Distribution of these genes and their percentage representation are listed in Table 4.  The preferred start codon in P. delaneyi strain Hulk is ATG (45.07%), but the start codons GTG (30.23%) and TTG (24.70%) remain significant. Such start codon preference is similar to that of the type strain Su06 T (ATG, 45.08%; GTG, 30.23%, and TTG, 24.69%). The closest relative, Hyperthermus butylicus, also had ATG as the preferred codon (38%), followed closely by TTG (37%), then GTG (25%). Another close relative, Aeropyrum pernix, had a preference for the TTG codon (52%), followed by ATG (28%) and GTG (20%) [25,26]. The draft genome of strain Hulk also contained one SSU RNA (Pdsh_09290), one LSU RNA (Pdsh_09285) and two copies of the 5S rRNA gene (Pdsh_05305 and Pdsh_07510). Table 3 shows additional genome statistics.
Archaeal origins of replication are often AT rich, contain one or more DNA unwinding elements, are in intergenic regions, and possess binding sites for origin binding proteins such as Cdc8 and Orc1 [27]. We identified in the genome of P. delaneyi strain Hulk two regions (contig 2, 205,575-205,655; and contig 1, 179,501-179,576) homologous to OR sequences in the DoriC database [28] assigned to Ignicoccus hospitalis (ORI10010050) and Sulfolobus islandicus (ORI10010117). The draft genome of strain Hulk also contained homologs of the DNA replication initiator proteins Cdc6 (Pdsh_04070) and Orc1 (Pdsh_05200), though none were located downstream of either of the predicted origins of replication, as is typical in other archaeal origins of replication [27]. Genes related to cell division and chromosome replication were identified in nearby locations (Pdsh_04035, Pdsh_05165, Pdsh_05230, and Pdsh_05215-05220).

Central metabolism
Consistent with the ability of strain Hulk to use with glucose-containing sugars such as cellobiose, (a β(1,4) glucose disaccharide) and starch (an α-glucan) as electron donors, we identified in the annotated draft genome genes encoding proteins for the oxidation and assimilation of glucose through a modified Embden-Meyerhof-Parnas glycolytic pathway (Fig. 4a). As is typical in many archaea [29], the modified EMP pathway replaces the glycolytic glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase enzymes with a glyceraldehyde-3phosphate ferredoxin oxidoreductase (Pdsh_09265) for the conversion of glyceraldehyde-3-phosphate into 3- The total is based on the total number of protein coding genes in the genome phosphoglyceric acid [29]. Interestingly, the pathway is also present in the type strain P. delayeni Su06 T , which, unlike strain Hulk, does not reportedly use sugars as electron donors and carbon sources [3]. The genome of strain Hulk also contains homologs of glyceraldehyde-3phosphate dehydrogenase (Pdsh_05050) and phosphoglycerate kinase (Pdsh_04830) enzymes, which could function during gluconeogenesis. Genome analyses also identified most of the enzymes of the tricarboxylic acid cycle (Fig. 4b). The genome encoded a pyruvate carboxylase (Pdsh_02345-02350), which fixes carbon as the bicarbonate anion (HCO 3 − ) to convert pyruvate into oxaloacetate. However, it lacked the citrate synthase enzyme that catalyzes the condensation of oxaloacetate and acetyl-CoA to generate citrate [30]. It is unlikely that the absence of a citrate synthase was the result of the incomplete genome coverage because genomes from other members of the family Pyrodictiaceae also lacked obvious citrate transporters (Table 5). Furthermore, the gene is also absent in the complete genomes of most other members of the order Desulfurococcales (Table 5).
Interestingly, we identified in the genome of strain Hulk transporters for acetate (ActP acetate permease, Pdsh_04815), a carbon source that, though common in hydrothermal environments [9,31,32], does not serve as electron donor for strain Hulk. Furthermore, the genome contained an enzyme for the conversion of acetate into acetyl-CoA (acetyl-CoA synthetase, Pdsh_04765) and several proteins involved in coenzyme A biosynthesis (Pdsh_01535, Pdsh_01545-01560, Pdsh_03865, Pdsh_03830-03880, Pdsh_08205, and Pdsh_08280). We also identified a pyruvate carboxylase enzyme (Pdsh_06115-06130), which fixes CO 2 in a reaction that converts acetyl-CoA into pyruvate to feed into the TCA cycle (Fig. 4b). This suggests that acetate carbon is assimilated in the TCA cycle yet it is not oxidized because it lacks a citrate synthase. Instead, the TCA cycle can run in reverse, that is, through reductive reactions from oxaloacetate to malate and so on (Fig. 4b). This reverse mode bypasses the need for the citrate synthase step and also enables cells to fix CO 2 . This could explain why strain Hulk grows autotrophically with H 2 and CO 2 (Fig. 5a) although its genome lacks key enzymes in other carbon fixation pathways such as the dicarboxylate/4hydroxybutyrate cycle (DC/HB), the 3-hydroxypropionate/ 4-hydroxybutyrate cycle (HP/HB), the fuchs-holo bi-cycle, the Calvin-Benson-Bassham cycle, and the Wood-Ljungdahl pathway. Indeed, the reverse TCA cycle has been proposed to be the major autotrophic pathway in hydrothermal vent environments [33].

Ribulose monophosphate pathway
The pentose phosphate pathway allows many bacteria and eukaryotes to generate reducing power and precursor metabolites needed for the synthesis of nucleotides and aromatic amino acids [29,34]. However, archaea are not known to have a complete oxidative PPP, and only Methanococci and Thermoplasmatales species have a complete non-oxidative PPP [29]. We identified an incomplete non-oxidative PPP in the genome of strain Hulk, but were unable to find homologs of any of the enzymes involved in the oxidative PPP. As is common  Citrate transporter --- within the Archaea, a combination of the incomplete non-oxidative PPP and ribulose monophosphate pathway replaces the oxidative PPP to provide reactions for the synthesis of nucleotide and aromatic amino acid precursors such as ribose-5-phosphate and erythrose-4-phosphate respectively [34]. Indeed, the genome of strain Hulk contains all of the genes of the ribulose monophosphate pathway, including a bifunctional D-arabino-3-hexulose 6-phosphate formaldehyde lyase/6phospho-3-hexuloisomerase (HPF/PHI) and ribose 5phosphate isomerase (Pdsh_00725) (Fig. 4c).
Hydrogen and formate as an electron donor during autotrophic growth H 2 and formate are abundant electron donors in hydrothermal vent environments, due to their continuous replenishment through intense serpentinization [6,9,35]. These environments are also rich in Fe(III) minerals, creating conditions optimal for the growth of hyperthermophilic iron reducers with H 2 and formate [10,13,[36][37][38]. Indeed, we demonstrated autotrophic growth with H 2 for strain Hulk (Fig. 5a) and growth with formate during the reduction of sulfur-containing electron acceptors (Fig. 3), Fe(III) citrate (Fig. 5b) and nitrate (Fig. 6) [39]. Consistent with the experimental results, we identified in the genome of strain Hulk a formate dehydrogenase gene cluster (Pdsh_ 05790-05805) encoding two formate dehydrogenase formation proteins (FdhE and FdhD), a putative formate transporter (FdhC), and the formate dehydrogenase α (FdhA) and β (FdhB) subunits. We also identified in the genome of strain Hulk genes encoding the small (Pdsh_08510) and large (Pdsh_09845) subunits of a NiFe hydrogenase, which catalyzes the reversible oxidation of H 2 [40]. Interestingly, the small NiFe hydrogenase subunit of strain Hulk (Pdsh_08510) did not contain the FeS domains that are needed to transfer electrons from the enzyme center to the electron-accepting substrate [40,41]. Yet an FeS-containing protein encoded a downstream gene (Pdsh_08500) could mediate this reaction.
The genome of strain Hulk also contained several of the proteins (HypABCDEF) required for the maturation of NiFe hydrogenases [42,43]. We identified, for example, the hydrogenase maturation protein HypA (Pdsh_01645 and Pdsh_04600) as well as a cluster containing the HypCDEF proteins (Pdsh_08530-08545). Interestingly, neither strain Hulk nor the type strain Su06 T contained homologs of HypB, a maturation protein that scavenges nickel at low concentrations for its incorporation in the active site of the NiFe hydrogenase [44]. The loss of HypB can be overcome if high levels of nickel are provided in the external environment [45]. This suggests that the nickel concentrations may be high in hydrothermal environments. Alternatively, yet to be identified nickel-scavenging proteins may be used for the maturation of NiFe hydrogenases in these archaea.
The oxidation of H 2 or formate by membrane-bound dehydrogenase enzymes contributes to the formation of a proton gradient, which is harnessed to generate energy for growth. Two NADH ubiquinone oxidoreductase (Nuo) proteins could boost the proton motive force. These proteins are encoded in two gene clusters, the first encoding subunits A-D and H-N (Pdsh_06005-05955) and the second coding subunits B-D, L-N and I (Pdsh_06660-06630). NADH ubiquinone oxidoreductase, also known as "respiratory complex 1" [40], is a proton pump able to transfer electrons from NADH to ubiquinone. As a result, the enzyme generates a proton gradient that can be harnessed to synthesize ATP [40].

Cellobiose and starch as electron donors and carbon sources
In contrast to the type strain Su06 T , which cannot use sugars as sole electron donors or carbon sources for growth [3], strain Hulk grew with cellobiose and starch as sole electron donors. Consistent with this, the genome of strain Hulk codes for a cellobiose phosphorylase (Pdsh_03720). This enzyme catalyzes the phosphorolysis of the cellulose disaccharide, thus conserving the energy in the β(1,4) glycosidic bond to generate glucose 1phosphate [46]. Searches for homologs of the Pdsh_03720 cellobiose phosphorylase in the sequenced genomes of other members of the Desulfurococcales only identified a protein with a low degree of homology (66% identity, 84% similarity) in Aeropyrum camini (BAN90220.1). By contrast, starch-degrading enzymes were widespread within the Desulfurococcales order. The genome of strain Hulk contained, for example, homologs of a glycosidase (Pdsh_09725) and an amylopullulanase (Pdsh_09650) enzyme. Glycosidase enzymes break the α(1,4) glycosidic bonds in the linear glucose strands of the polymer, whereas amylopullulanases also cleave the branched α(1,6) glycosidic linkages [47]. The genome of strain Hulk also encoded for several glycosyltransferases (Pdsh_02755 and Pdsh_03780), which are enzymes that catalyze the transfer of sugar moieties from activated donor molecules to specific acceptor molecules [48]. It is interesting to note that neither cellobiose nor starch are abundant in most deepsea environments [49]. Yet, the ability to scavenge these complex sugars has been proposed to allow cells to store carbon as glycogen [49]. However, the genome of strain Hulk and other members of the order Desulfurococcales lack obvious glycogen synthase genes. This could indicate that these organisms carry highly divergent glycogen synthases or that carbon is stored as polymers other than glycogen. Alternatively, yet to be identified sources of complex carbon compounds may be available in these hydrothermal vent environments that select for genomes containing genes involved in the oxidation and assimilation of complex sugars.

Metabolism of peptones
The availability of free peptides in hydrothermal marine vent systems [7,8] provides an abundant source of carbon and energy for microorganisms inhabiting these ecosystems. Indeed, peptone served as electron donor to support the growth of strain Hulk. However, growth was not supported with individual amino acids such as histidine, cysteine, and leucine as electron donors. The genome of strain Hulk contains three gene clusters that encode ABC branched amino acid transporters (Pdsh_10390-10,365, Pdsh_01510-01495, and Pdsh_07600-07575). Two gene clusters were identified that encode ABC oligopeptide transporters (Pdsh_05670-05665 and Pdsh_07635-07620), and each oligopeptide transport cluster had dipeptide transport subunits directly upstream (Pdsh_05675 and Pdsh_07640-07645). We also identified in the genome 13 putative peptidases, including a metallocarboxypeptidase (M32) (Pdsh_06240), three aminopeptidases (Pdsh_04150, Pdsh_05315, and Pdsh_07085), a metalloprotease (M50) (Pdsh_05275), and several omega peptidases (Pdsh_06520, Pdsh_06735, Pdsh_07355, Pdsh_07465, and Pdsh_07745). We also identified prolyl oligopeptidases (Pdsh_07465 and Pdsh_07695), an isoaspartyl aminopeptidase (Pdsh_07355), pyrrolidone-carboxylate peptidases (Pdsh_06520 and Pdsh_06735), an exoaminopeptidase (M42) (Pdsh_02360), and a D-aminopeptidase (M55) (Pdsh_08590). The arsenal of peptidase enzymes encoded in the genome of strain Hulk likely maximizes the efficient utilization of the peptides available at hydrothermal sites. Analysis of the complete genome of the type strain Su06 T also identified all of these genes, suggesting that the ability to grow with peptones is not restricted to strain Hulk only.  (7), and osmolarity (2% NaCl) in reference to controls without the electron donor (formate)

Iron respiration
The ability of strain Hulk to respire both insoluble (i.e. ferrihydrite) and chelated (i.e. Fe(III) citrate) forms of iron (Fig. 5) further confirms this is a widespread metabolic ability of hyperthermophilic organisms [1]. Mechanistic studies in the model hyperthermophilic archaeon Geoglobus ahangari demonstrated the need for direct contact between the cell and the iron mineral to transfer respiratory electrons, a reaction that is mediated by c-type cytochromes on the outer surface of the cell [13]. Supporting a similar mechanism, we identified in the genome of strain Hulk five genes carrying the conserved heme-binding motif (CXXCH) of c-type cytochromes (Pdsh_01720, Pdsh_01730, Pdsh_08690, Pdsh_09040, and Pdsh_10070) ( Table 6). All of these genes encoded proteins with signal peptides, as expected of proteins that are exported across the membrane, and none had homology to any known c-type cytochromes in the non-redundant NCBI database. The genome also contains 9 genes encoding homologues of the Ccm cytochrome maturation pathway (CcmABCDEFGHI). We identified, for example, ccmA (Pdsh_08980), ccmB (Pdsh_09025), ccmC (Pdsh_09010), and ccmF (Pdsh_01745 and Pdsh_09015). Although this cytochrome-c maturation system is found in many archaea [50,51], two of its components (CcmE and CcmH) are often missing or are highly divergent in archaeal genomes [52]. Indeed, the annotation of the genome strain Hulk did not identify a CcmE homolog but included one gene (Pdsh_09020) encoding a protein of the CcdA superfamily that is predicted to be a functional homolog of CcmH [50].
To maximize access to the iron minerals, microorganisms with a direct contact mechanism of electron transfer often rely on flagellar motility and chemotaxis to access the minerals [53]. Although most members of the family Pyrodictiaceae are reportedly not motile [54], cells of strain Hulk were motile and assembled several flagella (Fig. 2a). A cluster of genes (Pdsh_02535-02575) was annotated in the genome that comprised the minimal gene set (flaB, flaH, flaI, and flaJ) needed to assemble the archaeal flagellum or archaellum [27,[55][56][57] in addition to genes encoding for hypothetical proteins [(Pdsh_02550-02570) (Fig. 5c)]. The protein encoded by Pdsh_02555 had conserved domains described in the archaeal-flagellum superfamily and also had weak homology (27% identify, 44% similarity) to FlaG in Fervidicoccus fontis (AFH42914.1), a protein proposed to anchor the flagellar hook and filament [58]. The protein encoded by Pdsh_02560 belongs to the DD superfamily that has been implicated in the selective signaling of several complexes [59]. The presence of flaB, flaG, flaH, flaI, and flaJ homologs and the absence of strong hits to flaC, flaD and flaF is consistent with the type Fla2 gene cluster that is typical of most Crenarchaeota archaella [60]. However, the genome of strain Hulk did not contain any obvious chemotaxis genes. Furthermore, a search in the sequenced genomes within the order Desulfurococcales only identified chemotaxis genes in two species (Staphylothermus hellenicus and Ignisphaera aggregans). The lack of chemotaxis genes has been proposed to reflect an adaptive response of organisms to hydrothermal environments where intense fluid circulation replenishes electron acceptors and donors [13] Indeed, fluid circulation at the Hulk sulfide chimney vent is high [61]. Thus, the lack of chemotaxis genes within the genome of strain Hulk is consistent with lack of evolutionary pressure to sense chemical gradients.

Assimilatory and dissimilatory nitrogen metabolism
Nitrogen-containing compounds such as ammonium, nitrate, and nitrite are key intermediates in the marine nitrogen cycle and available for assimilation and respiration in anoxic marine environments [62][63][64]. Archaea, including members of the Crenarchaeota, also contribute to the cycling of nitrogen in deep submarine ecosystems [65,66]. Consistent with the availability of ammonium in these environments, the genome of strain Hulk contained two ammonium transporters clustered next to a nitrogen regulatory P-II protein (Pdsh_06180-06175 and Pdsh_00700-00695). Also present are a glutamine synthase (Pdsh_02890) for the cytoplasmic incorporation of ammonium into glutamine using α-ketoglutarate as a substrate and a glutamate dehydrogenase (Pdsh_05430), which uses NADPH with α-ketoglutarate and ammonium to produce glutamate [67]. The two enzymes form the GS/GOGAT pathway for the cytoplasmic assimilation of nitrogen from ammonium [67]. The genome also contained an NADPHdependent glutamate synthase (Pdsh_04375), which converts glutamine into glutamate [67].
By contrast, the genome of strain Hulk lacked nitrate or nitrite transporters for the assimilation of these oxidized nitrogen species. We identified instead homologs of the genes encoding the periplasmic bacterial, dissimilatory nitrate reductase Nap complex (NapADFGH) (Pdsh_02795 and Pdsh_02785-02770) (Fig. 6a). Consistent with the genome prediction, P. delaneyi strain Hulk also grew with nitrate as an electron acceptor (Fig. 6b). Also absent in the genome was an obvious homolog of a nitrite reductase. The closest gene (Pdsh_01730) encodes a putative cytochrome-c 7 protein containing an NrfH domain. In agreement with the genome prediction, strain Hulk did not use nitrite as an electron acceptor. This suggests that the strain reduces nitrate to nitrite but cannot carry out the complete denitrification of nitrate to N 2 reported for the type strain Su06 T [3].

Respiration of sulfur-containing electron acceptors
The annotation of the genome of strain Hulk revealed genes encoding proteins involved in assimilatory and dissimilatory metabolism of several sulfur-containing compounds (Fig. 3c). A manual search identified, for example, two genes (Pdsh_01045 and 01050) encoding homologs of the DmsA (25% identity, 41% similarity) and DmsB (29% identity, 42% similarity) subunits of the DMSO reductase complex of the iron-reducing bacterium Shewanella oneidensis (SO1429 and SO1430, respectively) [68]. However, the genome of strain Hulk did not contain a homolog of the membrane-bound DmsC protein that anchors the DMSO reductase complex to the periplasmic side of the inner membrane in this bacterium, consistent with the different biochemical composition of the bacterial and archaeal membranes. Yet strain Hulk was able to reduce AQDS and DMSO ( Fig. 3a-b). Also annotated was a multifunctional sulfate/thiosulfate ATP-binding transport protein (Pdsh_03910) for the transport of sulfate (SO 4 2− ) and thiosulfate (S 2 O 3 2− ) anions inside the cell (Fig. 3c). The genome also contained a sulfate adenylyltransferase (Pdsh_03730) for the ATP-adenylation of sulfate into adenosine-5′-phosphosulfate (APS), a key step in both the assimilatory and dissimilatory sulfate reduction pathways [69,70]. Once sulfate is activated, the assimilatory pathway converts APS to 3′-phosphoadenosine 5′-phosphosulfate (PAPS) via an adenylyl sulfate kinase (Pdsh_03715) [69]. PAPS is subsequently converted to sulfite (SO 3 2− ) by a phosphoadenosine phosphosulfate reductase (PAPS-R) (Pdsh_07330) [69,71]. Sulfite can also be produced in energy-generating reactions of the sulfate dissimilatory pathway. The canonical dissimilatory sulfate reduction pathway uses an adenylsulfate reductase enzyme to convert APS to sulfite. However, the genome of strain Hulk lacked a clear homolog. Instead, we identified a phosphoadenosine phosphosulfate reductase (PUA_PAPS) (Pdsh_01200 and Pdsh_01605), a bifunctional enzyme that contains nonspecific APS reductase and PAPS reductase domains. Sulfite is then reduced to sulfide by a membrane-bound sulfite reductase complex (Pdsh_08990-09020). The sulfide gas can either be expelled (dissimilatory pathway) or assimilated into sulfur-containing compounds such as cysteine via a cysteine synthase (Pdsh_03275).
We also identified a gene cluster (Pdsh_01045-01055) encoding the subunits of a thiosulfate reductase, which catalyzes the reduction of thiosulfate to sulfite in a reaction that generates sulfide gas (Fig. 3c). The sulfite can also be reduced to sulfide by the membrane-bound sulfite reductase of the dissimilatory sulfate reduction pathway (Fig. 3c). Consistent with the genome predictions, strain Hulk grew in media with formate as the electron donor, utilizing sulfate, thiosulfate and sulfite as terminal electron acceptors.
Interestingly, all sulfate/thiosulfate reduction genes of strain Hulk had clear homologs in genes of the type strain Su06 T annotated as hypothetical proteins, although Su06 T reportedly cannot grow with elemental sulfur and thiosulfate [3]. For example, the gene cluster encoding the dissimilatory sulfite reductase of strain Hulk is also present in the type strain Su06 T (Pyrde_0485-0492). The first five genes in the Su06 T cluster (Pyrde_486-490) are 100% identical to the homologues in strain Hulk (Pdsh_08990-09010), whereas the last two genes (Pyrde_490 and 491) are 99 and 98% identical to Pdsh_09015 and Pdsh_09020 in the genome of strain Hulk. This suggests that the reduction of sulfite, and perhaps of other sulfur-containing compounds, may be an unsuspected respiratory capability of the type strain Su06 T as well. If so, the metabolic hallmark of the family Pyrodictiaceae is, as proposed earlier [1], the ability of these hyperthermophilic archaea to use sulfur-containing electron acceptors.

Conclusions
P. delaneyi strain Hulk is the first hyperthermophilic archaeon able to support growth through the reduction of iron, nitrate, and sulfur-containing compounds. Like most other members of the family Pyrodictiaceae [1,3], strain Hulk grows autotrophically with H 2 . However, it is the first member of the family reported to use formate and sugars as electron donors. It also shares with two species in the family, P. abyssi AV2 T [4] and H. butylicus PLM1-5 T [5], the ability to use peptides as electron donors. The annotated genome revealed the molecular basis for such remarkable metabolic versatility, with numerous pathways for autotrophic and heterotrophic growth using electron acceptors that are abundant in hydrothermal vent environments. This suggests that hyperthermophilic genomes are under selective pressure to maximize the use of the available resources.