|Home | About | Journals | Submit | Contact Us | Français|
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.
The unifying metabolic feature of the first five species described in the family 10.1601/nm.55 , in the archaeal order 10.1601/nm.30, was for long their ability to respire sulfur-containing electron acceptors, mainly elemental sulfur (S0), thiosulfate (S2O3 2−) and sulfite (SO3 2−) . 10.1601/nm.63 1AT reduces nitrate in addition to thiosulfate . Yet the recently described 10.1601/nm.28674 Su06T reduces nitrate and iron but cannot use sulfur or thiosulfate .
Here we report on the isolation and genome sequencing and annotation of a novel strain of 10.1601/nm.28674, 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 sulfur-containing electron acceptors. Furthermore, strain Hulk is the first member of the 10.1601/nm.55 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, 10.1601/nm.58 AV2T  and 10.1601/nm.61 PLM1–5T . The ability to oxidize formate and peptides is of special environmental significance, as these are abundant electron donors in hydrothermal marine vent systems [6–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.
A novel strain of 10.1601/nm.28674 designated strain Hulk (Fig. (Fig.1)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  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 N2:CO2 (80:20%, v/v, 101 kPa). Isolation was in the same medium solidified with GELRITE gellan gum (Sigma-Aldrich), as previously described . All incubations during enrichments and isolation were at 100 °C.
Strain Hulk was tested for electron donor/acceptor use in the modified marine medium described above, but incubating the cultures in the dark at the optimum growth temperature (90 °C). The following potential electron donors were tested: H2 (H2:CO2, 80:20%, v/v, 101 kPa), lactate (10 mM), pyruvate (10 mM), formate (10 mM), butyrate (2.5 mM), yeast extract (0.01%, wt/vol), methanol (5 mM), ethanol (5 mM), 1-propanol (5 mM), 3-methyl-1-butanol (5 mM), stearate (1 mM), starch (0.1%, wt/vol), D-glucose (5 mM), fructose (5 mM), maltose (2.5 mM), melibiose (2.5 mM), mannose (5 mM), galactose (5 mM), palmitate (1 mM), acetate (10 mM), malate (10 mM), succinate (10 mM), citrate (10 mM), fumarate (10 mM), peptone (0.1%, wt/vol), valeric acid (5 mM), propionate (5 mM), alanine (12 mM), histidine (6 mM), proline (8.7 mM), glycine (14 mM), isoleucine (7.6 mM), aspartic acid (7.6 mM), glutamic acid (6.8 mM), arginine (6 mM), L-cysteine (8 mM), serine (10 mM), asparagine (6 mM), cellulose (0.24%, wt/vol), chitin (0.4%, wt/vol), and cellobiose (0.2%, wt/vol). The electron acceptors tested were: uranium (1 mM), vanadium (1 mM), Fe(III) pyrophosphate (10 mM), sulfate (14 mM), thiosulfate (10 mM), sulfite (4 mM), selenate (10 mM), selenite (10 mM), arsenate (10 mM), nitrate (10 mM), nitrite (1 mM), poorly crystalline iron (Fe[III]) oxides (100 mmol l−1), manganese oxides (poorly crystalline Mn[IV], 20 mmol l−1), malate (10 mM), fumarate (50 mM), iron (Fe[III]) citrate (50 mM), dimethyl sulfoxide (DMSO, 1 mM), anthraquinone-2,6-disulfonate (AQDS, 5 mM), goethite (50 mmol l−1), and hematite (50 mmol l−1).
Table Table11 shows the general features of the novel isolate. The strain was a chemolithoautotroph that coupled the reduction of Fe(III) (provided as poorly crystalline Fe(III) oxides or soluble Fe(III) citrate) to the oxidation of H2, formate, peptone, cellobiose, and starch. This contrasts with the type strain Su06T, which cannot grow with formate and sugars . Cells of 10.1601/nm.28674 strain Hulk from cultures with formate and Fe(III) citrate or sulfate were cocci of ca. 1.0 μm in diameter, motile, and lophotrichously flagellated (up to 10 flagella per cell were apparent; Fig. Fig.2a).2a). Scanning electron micrographs of cells grown with formate (10 mM) as the electron donor and sulfate (14 mM) as the electron acceptor revealed extensive vesiculation, with some membrane vesicles appearing as filamentous aggregates attached to the cell surface (Fig. (Fig.2b).2b). Like strain Su06T , strain Hulk also gained energy for growth from the reduction of nitrate. But, unlike the type strain, strain Hulk was capable of growth with most of the sulfur-containing electron acceptors tested such as AQDS, DMSO, sulfate, and thiosulfate (Fig. (Fig.3).3). Also notable is the higher upper temperature limit for growth of strain Hulk (105 °C) compared to the type strain Su06T (97 °C) . Additionally, strain Hulk has a broader range of growth-supporting temperatures (70–105 °C, optimum at 90 °C) and pH (pH 3.3–7.7, optimum at pH 6.7) than strain Su06T (82–97 °C, optimum at 90–92 °C; and pH 0.9–3.6, optimum at pH 1.9) . The range of salt concentrations that supported growth of strain Hulk (0.5–4.5% (w/v) NaCl, with optima at 2% NaCl) was broader than those reported for the type strain (0.9–3.6% (w/v), optimum at 1.9% NaCl) .
The draft genome of 10.1601/nm.28674 strain Hulk contained a single 16S rRNA gene, whose sequence was 100% identical to the 16S rDNA sequence in the complete genome of 10.1601/nm.28674 strain Su06T . Pair-wise genome comparisons between the two genomes  revealed an average nucleotide identity (ANI) of 99.77%, above the cutoff (92%) established for bacterial species definition . 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. (Fig.1).1). The 10.1601/nm.61 was the nearest neighbor outside of the 10.1601/nm.28674 species group (99% identity).
10.1601/nm.28674 strain Hulk was sequenced and annotated based on its phylogenetic position and its metabolic versatility. Genome comparisons  with the recently described species type strain Su06T  and other members of the 10.1601/nm.55 were performed to provide novel insights into this archaeal family and the metabolic potential of their members. The draft genome of 10.1601/nm.28674 strain Hulk presented here is contained within 9 contigs with an average coverage of at least 380× (Table (Table2).2). This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession NCQP00000000 (Table (Table2).2). The version described in this paper is version NCQP01000000. The genome project summary can be viewed at Genomes Online Database (Ga0169944)  (Table (Table22).
Genomic DNA was extracted from exponentially-grown cells from formate-Fe(III) citrate cultures incubated at 90 °C as previously described . Cells were harvested at room temperature by centrifugation for 20 min, washed with a marine wash buffer , and extracted using the MoBio Powersoil kit, except that an additional lysis buffer step was included, as reported elsewhere . Quality of genomic DNA was verified by nanodrop and gel electrophoresis.
The assembly presented was generated from three paired-end Illumina libraries . 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  and ConDeTri  and analyzed with the genome assembly program Velvet  to generate three independent draft de novo genome assemblies. A consensus assembly was then generated using Mix  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×.
The genome annotation used the Rapid Annotations via Subsystems Technology (RAST) server . Coding regions were identified with RAST’s GLIMMER tool . 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 . 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 . 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  and the TMHMM Server (v2.0) , 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.
The draft genome sequence of 10.1601/nm.28674 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 (Table3).3). The draft genome size is close to that reported for the complete genome of the type strain of 10.1601/nm.28674 Su06T (2,023,836 bp) . 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 (Table3).3). Seventeen pseudogenes were identified by the Prokaryotic Genomes Automatic Annotation Pipeline . 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 Table44.
The preferred start codon in 10.1601/nm.28674 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 Su06T (ATG, 45.08%; GTG, 30.23%, and TTG, 24.69%). The closest relative, 10.1601/nm.61, also had ATG as the preferred codon (38%), followed closely by TTG (37%), then GTG (25%). Another close relative, 10.1601/nm.39 , 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 Table33 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 . We identified in the genome of 10.1601/nm.28674 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  assigned to 10.1601/nm.10785 (10.1601/strainfinder?urlappend=%3Fid%3DORI+10010050) and 10.1601/nm.66 islandicus (10.1601/strainfinder?urlappend=%3Fid%3DORI+10010117). 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 . Genes related to cell division and chromosome replication were identified in nearby locations (Pdsh_04035, Pdsh_05165, Pdsh_05230, and Pdsh_05215–05220).
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. (Fig.4a).4a). As is typical in many archaea , the modified EMP pathway replaces the glycolytic glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase enzymes with a glyceraldehyde-3-phosphate ferredoxin oxidoreductase (Pdsh_09265) for the conversion of glyceraldehyde-3-phosphate into 3-phosphoglyceric acid . Interestingly, the pathway is also present in the type strain P. delayeni Su06T, which, unlike strain Hulk, does not reportedly use sugars as electron donors and carbon sources . The genome of strain Hulk also contains homologs of glyceraldehyde-3-phosphate 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. (Fig.4b).4b). The genome encoded a pyruvate carboxylase (Pdsh_02345–02350), which fixes carbon as the bicarbonate anion (HCO3 −) to convert pyruvate into oxaloacetate. However, it lacked the citrate synthase enzyme that catalyzes the condensation of oxaloacetate and acetyl-CoA to generate citrate . 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 10.1601/nm.55 also lacked obvious citrate transporters (Table (Table5).5). Furthermore, the gene is also absent in the complete genomes of most other members of the order 10.1601/nm.30 (Table (Table55).
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 CO2 in a reaction that converts acetyl-CoA into pyruvate to feed into the TCA cycle (Fig. (Fig.4b).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. (Fig.4b).4b). This reverse mode bypasses the need for the citrate synthase step and also enables cells to fix CO2. This could explain why strain Hulk grows autotrophically with H2 and CO2 (Fig. (Fig.5a)5a) although its genome lacks key enzymes in other carbon fixation pathways such as the dicarboxylate/4-hydroxybutyrate 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 .
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 10.1601/nm.143 and 10.1601/nm.359 species have a complete non-oxidative PPP . 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 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 . 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/6-phospho-3-hexuloisomerase (HPF/PHI) and ribose 5-phosphate isomerase (Pdsh_00725) (Fig. (Fig.4c4c).
H2 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 H2 and formate [10, 13, 36–38]. Indeed, we demonstrated autotrophic growth with H2 for strain Hulk (Fig. (Fig.5a)5a) and growth with formate during the reduction of sulfur-containing electron acceptors (Fig. (Fig.3),3), Fe(III) citrate (Fig. (Fig.5b)5b) and nitrate (Fig. (Fig.6)6) .
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 H2 . 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 Su06T contained homologs of HypB, a maturation protein that scavenges nickel at low concentrations for its incorporation in the active site of the NiFe hydrogenase . The loss of HypB can be overcome if high levels of nickel are provided in the external environment . 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 H2 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” , 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 .
In contrast to the type strain Su06T, which cannot use sugars as sole electron donors or carbon sources for growth , 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 1-phosphate . Searches for homologs of the Pdsh_03720 cellobiose phosphorylase in the sequenced genomes of other members of the 10.1601/nm.30 only identified a protein with a low degree of homology (66% identity, 84% similarity) in 10.1601/nm.40 (BAN90220.1). By contrast, starch-degrading enzymes were widespread within the 10.1601/nm.30 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 . 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 . It is interesting to note that neither cellobiose nor starch are abundant in most deep-sea environments . Yet, the ability to scavenge these complex sugars has been proposed to allow cells to store carbon as glycogen . However, the genome of strain Hulk and other members of the order 10.1601/nm.30 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.
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 Su06T also identified all of these genes, suggesting that the ability to grow with peptones is not restricted to strain Hulk only.
The ability of strain Hulk to respire both insoluble (i.e. ferrihydrite) and chelated (i.e. Fe(III) citrate) forms of iron (Fig. (Fig.5)5) further confirms this is a widespread metabolic ability of hyperthermophilic organisms . Mechanistic studies in the model hyperthermophilic archaeon 10.1601/nm.413 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 . 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 (Table6).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 . 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 .
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 . Although most members of the family 10.1601/nm.55 are reportedly not motile , cells of strain Hulk were motile and assembled several flagella (Fig. (Fig.2a).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–57] in addition to genes encoding for hypothetical proteins [(Pdsh_02550–02570) (Fig. (Fig.5c)].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 10.1601/nm.20073 (AFH42914.1), a protein proposed to anchor the flagellar hook and filament . The protein encoded by Pdsh_02560 belongs to the DD superfamily that has been implicated in the selective signaling of several complexes . 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 10.1601/nm.2 archaella . However, the genome of strain Hulk did not contain any obvious chemotaxis genes. Furthermore, a search in the sequenced genomes within the order 10.1601/nm.30 only identified chemotaxis genes in two species (10.1601/nm.46 and 10.1601/nm.10019). 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  Indeed, fluid circulation at the Hulk sulfide chimney vent is high . Thus, the lack of chemotaxis genes within the genome of strain Hulk is consistent with lack of evolutionary pressure to sense chemical gradients.
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–64]. Archaea, including members of the 10.1601/nm.2, 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 . The two enzymes form the GS/GOGAT pathway for the cytoplasmic assimilation of nitrogen from ammonium . The genome also contained an NADPH-dependent glutamate synthase (Pdsh_04375), which converts glutamine into glutamate .
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. (Fig.6a).6a). Consistent with the genome prediction, 10.1601/nm.28674 strain Hulk also grew with nitrate as an electron acceptor (Fig. (Fig.6b).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 N2 reported for the type strain Su06T .
The annotation of the genome of strain Hulk revealed genes encoding proteins involved in assimilatory and dissimilatory metabolism of several sulfur-containing compounds (Fig. (Fig.3c).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 10.1601/nm.2931 (SO1429 and SO1430, respectively) . 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. (Fig.3a3a--bb).
Also annotated was a multifunctional sulfate/thiosulfate ATP-binding transport protein (Pdsh_03910) for the transport of sulfate (SO4 2−) and thiosulfate (S2O3 2−) anions inside the cell (Fig. (Fig.3c).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) . PAPS is subsequently converted to sulfite (SO3 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 non-specific 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. (Fig.3c).3c). The sulfite can also be reduced to sulfide by the membrane-bound sulfite reductase of the dissimilatory sulfate reduction pathway (Fig. (Fig.3c).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 Su06T annotated as hypothetical proteins, although Su06T reportedly cannot grow with elemental sulfur and thiosulfate . For example, the gene cluster encoding the dissimilatory sulfite reductase of strain Hulk is also present in the type strain Su06T (Pyrde_0485–0492). The first five genes in the Su06T 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 Su06T as well. If so, the metabolic hallmark of the family 10.1601/nm.55 is, as proposed earlier , the ability of these hyperthermophilic archaea to use sulfur-containing electron acceptors.
10.1601/nm.28674 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 10.1601/nm.55 [1, 3], strain Hulk grows autotrophically with H2. 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, 10.1601/nm.58 AV2T  and 10.1601/nm.61 PLM1–5T , 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.
We thank James Holden for the chimney samples, James Kremer for assistance with SEM.
This work was supported by Strategic Partnership Grant 08-SPG-Full-1400 from the Michigan State University Foundation to GR and KK. MM acknowledges support by a GAANN fellowship from the Department of Education.
KK conceived the project and, together with LD and GR, wrote the manuscript. LD annotated and analyzed the genome of strain Hulk, CM generated all of the physiologic data, RRS sequenced the genome of strain Hulk., SKS and MM assembling the genome presented here. MM extracted genomic DNA from P. delaneyi strain Hulk. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.