Proteome Analyses of Hydrogen-producing Hyperthermophilic Archaeon Thermococcus onnurineus NA1 in Different One-carbon Substrate Culture Conditions*

Thermococcus onnurineus NA1, a sulfur-reducing hyperthermophilic archaeon, is capable of H2-producing growth, considered to be hydrogenogenic carboxydotrophy. Utilization of formate as a sole energy source has been well studied in T. onnurineus NA1. However, whether formate can be used as its carbon source remains unknown. To obtain a global view of the metabolic characteristics of H2-producing growth, a quantitative proteome analysis of T. onnurineus NA1 grown on formate, CO, and starch was performed by combining one-dimensional SDS-PAGE with nano UPLC-MSE. A total of 587 proteins corresponding to 29.7% of the encoding genes were identified, and the major metabolic pathways (especially energy metabolism) were characterized at the protein level. Expression of glycolytic enzymes was common but more highly induced in starch-grown cells. In contrast, enzymes involved in key steps of the gluconeogenesis and pentose phosphate pathways were strongly up-regulated in formate-grown cells, suggesting that formate could be utilized as a carbon source by T. onnurineus NA1. In accordance with the genomic analysis, comprehensive proteomic analysis also revealed a number of hydrogenase clusters apparently associated with formate metabolism. On the other hand, CODH and CO-induced hydrogenases belonging to the Hyg4-II cluster, as well as sulfhydrogenase-I and Mbx, were prominently expressed during CO culture. Our data suggest that CO can be utilized as a sole energy source for H2 production via an electron transport mechanism and that CO2 produced from catabolism or CO oxidation by CODH and CO-induced hydrogenases may subsequently be assimilated into the organic carbon. Overall, proteomic comparison of formate- and CO-grown cells with starch-grown cells revealed that a single carbon compound, such as formate and CO, can be utilized as an efficient substrate to provide cellular carbon and/or energy by T. onnurineus NA1.

Hyperthermophilic archaea can use a wide variety of carbon and energy sources. Hyperthermophiles are widely distributed in extreme habitats such as deep-sea thermal vents, hot springs, and deep oil reservoirs (1)(2)(3). So far, the most frequently studied hyperthermophiles are from the genera Thermococcus and Pyrococcus, which belong to the order Thermococcales (4). These are ecologically important hyperthermophilic archaea for understanding the physiology and metabolic activity of microbial consortia within marine hotwater ecosystems. Members of the order Thermococcales are anaerobic heterotrophs that utilize various complex substrates with elemental sulfur (S 0 ) or protons as electron acceptors (4 -6). Unlike other Thermococcales, Thermococcus strain AM4 (7) and Thermococcus onnurineus NA1 (8) are capable of lithotrophic carbon monoxide-dependent hydrogenogenic growth. These Thermococcus strains use CO as a carbon and energy source by converting it into carbon dioxide (CO 2 ). In addition, several hyperthermophilic archaea of the genus Thermococcus can grow with formate as an electron donor, producing hydrogen gas (9). T. onnurineus NA1 is a sulfur-reducing hyperthermophilic archaeon isolated from a deep sea hydrothermal vent area in the Eastern Manus Basin of Papua New Guinea (10). It is one of the more metabolically versatile hyperthermophiles in that it can use one-carbon (C1) substrates such as formate and CO, as well as multi-carbon substrates such as starch as the sole source of carbon and energy. Despite new insights into the mechanisms involved in the coupling of formate oxidation to energy conservation in T. onnurineus NA1 (9), the metabolic processes unique to growth on formate remain to be characterized. Typically, Most enzymes in the catabolic pathway of hyperthermophilic archaea have been identified and characterized. However, the anabolic pathway used by T. onnurineus NA1 to anaerobically convert these one-carbon substrates into cellular carbon is poorly understood. Therefore, study of one-carbon metabolism involved in H 2 -producing growth in T. onnurineus NA1 is of considerable interest because one-carbon substrates (i.e. methanol, carbon dioxide, and formate) are important for energy metabolism and carbon fixation pathways in some archaea (e.g. methanogenic archaea) and thus are expected to have many potential applications to industrial processes.
Recently, complete genome sequences have been determined for a number of Thermococcales, including three representative Pyrococcus species (11)(12)(13) and four Thermococcus strains: Thermococcus sibiricus (6), T. onnurineus NA1 (8), Thermococcus kodakaraensis KOD1 (14), and Thermococcus gammatolerans (15). Comparative genomic studies on these Thermococcales have paved the way for exploring the features and functions of genes involved in major metabolic pathways. Using a whole genome shotgun approach, the structure of the T. onnurineus NA1 genome was accurately determined, indicating that it contains 1,847,607 base pairs and 1976 predicted open reading frames (8). Moreover, genetic predictions based on these genomic data have provided clues for identifying genes correlated with biological functions (6,8,14). In particular, genomic analysis of T. onnurineus NA1 revealed several distinct hydrogenase gene clusters involved in H 2 metabolism and carboxydotrophic pathways (16). Compared with genome analysis, only a handful of proteomic studies have been performed on Thermococcus strains. In T. gammatolerans, the first archaeal genome-wide proteome investigation was performed using an LC-MS/MS shotgun approach at the primary genome annotation stage (15). This study demonstrated that T. gammatolerans utilizes various metabolic pathways during growth in nutrient-rich medium (15). However, few studies have examined changes in the protein expression profile of T. onnurineus NA1 (17,18). A previous 2-DE/MS-MS proteome analysis provided the first global view of the metabolic pathways of T. onnurineus NA1 during heterotrophic growth (17). Using an SDS-PAGE/LC-MS/MS shotgun proteomic approach, we recently made significant progress in identifying the metabolic enzymes specific for hydrogen production under carboxydotrophic culture conditions (18). Although recent genome analyses and proteomic characterizations have revealed general genomic features and metabolic pathways, a comprehensive understanding of the overall metabolism of T. onnurineus NA1 is lacking. Furthermore, a quantitative comparative analysis of all enzymes needed for central carbon flow, electron transfer, and/or energy conservation under different growth conditions has not been published.
Here we describe a quantitative proteome analysis of T. onnurineus NA1 cells using one-dimensional SDS-PAGE cou-pled with nano-UPLC-MS E . 1 We examined the metabolism of T. onnurineus NA1 during H 2 -producing growth on different substrates. The major metabolic pathways proposed from a genomic analysis were characterized by comparing the protein expression profiles of cultures grown in the presence of formate, CO, and starch as sole carbon and/or energy sources. For the first time, we can draw a more complete metabolic picture of T. onnurineus NA1 at the protein level, which may facilitate a greater understanding of the metabolic adaptation of hyperthermophiles to extreme environments.

SDS-PAGE and Trypsin
Digestion-Protein content was quantified using a Bradford assay kit (Bio-Rad). Crude cell extracts were separated by 12% SDS-PAGE (Mini-Protean; Bio-Rad) according to a method described previously (20). The gels were stained with Coomassie Brilliant Blue R250, and 100 g of proteins were loaded on the gel. The lanes were divided into eight fractions according to molecular weight, and slices were prepared for trypsin digestion (Promega, Madison, WI). The sliced gels were destained in a solution containing 50% ACN and 10 mM ammonium bicarbonate (ABC), and then rinsed with 100% ACN to stop the reaction. After drying in a speed vacuum concentrator (Mivac Quattro; Genevac Ltd.), the gels were incubated in a solution containing 10 mM DTT in 100 mM ABC at 56°C for 1 h to reduce protein disulfide bonds and then in the same volume of 55 mM iodoacetamide in 100 mM ABC added to alkylate cysteines in the dark for 45 min. The gels were washed with three volumes of distilled water and dried using a speed vacuum concentrator. After swelling with 300 l of 50 mM ABC, the proteins were digested with 15 l of trypsin (0.1 g/l) at 37°C for 18 h. The digested peptides were recovered by extracting twice with a solution containing 50 mM ammonium bicarbonate, 50% ACN, and 5% TFA. Peptide extracts were pooled at each gel line, lyophilized, and stored at Ϫ80°C until required.
Analysis by Nano-UPLC-MS E Tandem Mass Spectrometry and Quantitative Analysis-Separations were performed on a nano-UPLC C18 RP column (75 m ϫ 250 mm; particle size, 1.7 m) and an enrichment Symmetry C18 RP column (180 m ϫ 20 mm; particle size, 5 m) using a nano-ACQUITY Ultra Performance TM chromatography system (Waters Corporation).
Tryptic-digested peptides (5 l) were loaded onto the enrichment column with mobile phase A (3% acetonitrile in water with 0.1% formic acid). A step gradient was used at a flow rate of 300 nl/min. This included a 3-40% mobile phase B (97% acetonitrile in water with 0.1% formic acid) over 95 min and a 40 -70% mobile phase B over 20 min, followed by a sharp increase to 80% B within 10 min. Sodium formate (1 mol/min) was used to calibrate the TOF analyzer in the range of m/z 50 -2000, and [Glu 1 ]-fibrinopeptide (m/z 785.8426) was used for lock mass correction. During data acquisition, the collision energy of low energy MS mode and elevated energy mode (MS E ) were set to 4 and 15-40 eV energy ramping mode, respectively. One cycle of MS and MS E was performed every 3.2 s. In each cycle, MS spectra were acquired for 1.5 s with a 0.1-s interscan delay (m/z 300 -1990), and the ions exceeding 50 counts were selected for MS E fragmentation in the collision cell (m/z 50 -2000).
LC-MS E data were processed and searched using ProteinLynx GlobalServer version 2.3.3 (Waters Corporation). Proteins were identified by searching the T. onnurineus NA1 database on the NCBI website (http://www.ncbi.nlm.nih.gov), including 4492 entries. For the ProteinLynx GlobalServer search, the parent ion tolerance was set at 100 ppm, and the fragment ion tolerance was set at 0.2 Da. Carbamidomethylation (ϩ57 Da) of cysteine and methionine oxidation (ϩ16 Da) were chosen as the fixed and variable modifications, respectively. Analysis of quantitative changes in protein abundance, based on measurements of peptide ion peak intensities observed in low collision energy mode (MS) in a triplicate set, was carried out using Expression TM software (version 2). To normalize comparative proteomic data, the "auto normalization" function of ProteinLynx GlobalServer was used. Protein identification was allowed only if the confidence was greater than 95% on the basis of the IDENTITY E algorithm (21). The Expression software generated results as an exact mass retention time table with quantitative protein and peptide values.

RESULTS AND DISCUSSION
H2-producing Growth on CO, Formate, and Starch-A previous study showed that T. onnurineus NA1 produced H 2 when cultured with formate, CO, and starch as growth substrates (22). The highest level of growth and H 2 production was observed in cultures grown on formate. The maximal yield of H 2 was also achieved on formate or CO (22). These results indicate that formate and CO can serve as major energy sources for T. onnurineus NA1. We expected that T. onnurineus NA1 would exhibit specific responses to these single-carbon substrates, giving rise to changes in the proteome that would subsequently permit us to identify specific metabolic proteins involved in H 2 -producing growth.
Identification of Proteins Expressed during Growth on Different Substrates-To obtain a more complete picture of the metabolism involved in growth, a genome-wide proteome analysis of T. onnurineus NA1 grown on formate, CO, and starch was performed by one-dimensional SDS-PAGE sample fractionation and LC-MS E -based protein identification. Compared with previous approaches, the nano-UPLC-MS E method based on label-free quantitative analysis enables large scale comparison of all detectable proteins in complex samples separated by one-dimensional PAGE (21,23).
Nano-UPLC-coupled Q-TOF tandem mass spectrometry analysis was used to identify, quantify, and compare individual proteins among CO-grown, formate-grown, and starchgrown cells. In total, 587 distinct proteins with more than two peptide fragmentations were identified at least twice across the experiment, which accounts for 29.7% of the 1976 genes predicted in the T. onnurineus NA1 genome. A list of identified proteins with quantitative data is provided in supplemental Table S1. Of the 587 proteins, we defined those whose expression was up-regulated Ն1.5 times or down-regulated Յ0.67 times with a p value Ͻ 0.05 as differentially regulated under the particular growth condition (Table I). A total of 447 proteins for starch versus CO-grown cells, 453 for starch versus formate-grown cells, and 422 for CO versus formategrown cells were selected. Using our criterion (Ն1.5-fold change), 91, 97, and 112 proteins were found to be differentially abundant in starch versus CO-grown, starch versus formate-grown, and CO versus formate-grown cells, respectively (Table I). These differentially abundant proteins were consid- ered candidates for executing the predominant roles in the energy and carbon metabolisms induced by formate, CO, and starch.
The identified proteins were categorized by their predicted metabolic function (supplemental Table S1). The proteins were classified into 16 different categories based on the Kyoto Encyclopedia of Genes and Genomes database. Those that were significantly homologous with proteins of T. kodakaraensis KOD1 were classified into the same functional categories. However, the T. onnurineus NA1 proteome is not completely annotated; thus, 119 of the 587 identified proteins were assigned as "hypothetical." Notably, large proportions (28, 21, and 27%) of the identified proteins in starch versus COgrown, starch versus formate-grown, and CO versus formategrown cells, respectively, were involved in energy metabolism. This result is consistent with the idea that T. onnurineus NA1 can grow under chemolithotrophic culture conditions by utilizing formate and/or CO as energy sources, because this organism possesses a large variety of metabolic pathways and enzymes involved in energy metabolism and H 2 production.
Overview of General Metabolism in T. onnurineus NA1-Comparative proteomics enabled us to unambiguously identify the major metabolic pathways utilized by T. onnurineus NA1 during growth on formate, CO, and starch as sources of energy and/or carbon. The proteome induced in starch was used as a control. An overview of the putative metabolic pathways is shown in Fig. 1. Our proteomic analysis identified most of the enzymes involved in the modified Embden-Meyerhof-Parnas pathway, the nonoxidative branch of the pentose phosphate pathway (PPP), amino acid degradation, aldehyde metabolism, and the pseudo-TCA cycle ( Fig. 1 and Table II). These central metabolic pathways are similar to those proposed in genomic analyses of other Thermococcales (14, 15) but provide more information regarding which metabolic enzymes are specifically induced during growth on different substrates. In particular, the Embden-Meyerhof-Parnas pathway includes glyceraldehyde-3-phosphate dehydrogenase (TON_0639), 2-phosphoglycerate kinase (TON_0742), and FBPase (TON_1497), key enzymes in gluconeogenesis (24,25); these enzymes were strongly up-regulated during growth on formate. Intriguingly, unlike the FBPase (Class V) of thermophilic Archaea including T. kodakaraensis KOD1 (26,27), the recently characterized FBPase (TON_1497) from T. onnurineus NA1 has a different range of substrate specificities, although it has significant sequence similarity to T. kodakaraensis KOD1 (28). T. onnurineus NA1 also possesses key enzymes of the methylmalonyl pathway for propionate utilization, recently revealed in the thermophilic sulfur-reducing bacterium, Deferribacter desulfuricans SSM1 (29).
In the present study, we found that most hydrogenases, as well as enzymes involved in the stress response, were prominently expressed during growth on formate or CO (Table II  and supplemental Table S1). Moreover, as recently described in other Thermococcales (14,15), various ABC-type trans-porters, monovalent cation transporters that were differentially expressed during growth on formate, CO, and starch were identified. T. onnurineus NA1 harbors two TRAP-type transporter components (TON_0172 and TON_0174) in its genome, as does T. kodakaraensis KOD1 (14). In our proteomic analysis, only one TRAP-type transporter (TON_0174) was strongly up-regulated in cells grown on starch. The use of TRAP-type transporters might enable these bacteria to adapt to Na ϩ -rich and nutritionally poor environments by lowering the energetic cost of catalyzing the uptake of specific substrates (30). In addition, a branched chain amino acid ABC transporter (braC-4, TON_0162) was expressed in T. onnurineus NA1 cells during growth on starch. As in T. gammatolerans (15), an ABC-type dipeptide/oligopeptide transporter system (TON_1764 and TON_1768) was also found in the T. onnurineus NA1 proteome.
In contrast, most of the proteins involved in DNA metabolism were up-regulated during growth on formate compared with starch (Table II and supplemental Table S1). Currently, we do not know why formate induces DNA repair proteins in T. onnurineus NA1. However, formic acid depurinates DNA in vitro, leading to the loss of a purine from the DNA backbone (31). In vivo, depurination causes the formation of apurinic or apyrimidinic sites in DNA, which can result in strand breaks or point mutations by adenine mispairing during DNA synthesis (32,33). This depurination rate is greatly accelerated at low pH and elevated temperatures (34,35). Thus, the addition of formate at a high temperature may lead to the inhibition of DNA synthesis by damaging the DNA molecule. However, the spontaneous mutation rates of thermophilic archaea are reported to be lower than those of other mesophilic bacteria, despite the higher chemical instability of their DNA (36,37). Therefore, hyperthermophiles, including T. onnurineus NA1, seem to possess molecular strategies for DNA repair, to compensate for the intrinsic instability of DNA at their preferred growth temperatures and under acidic conditions (36,38).
Gluconeogenesis and Pentose Phosphate Pathway-T. onnurineus NA1 can grow using formate as its sole source of energy, producing hydrogen gas as a by-product (9). However, it is not clear whether formate can be used as a carbon source, and little is known about the metabolic pathways of formate assimilation. In T. onnurineus NA1, a potential reaction for formate assimilation is direct or indirect conversion of formate and acetyl-CoA into pyruvate (Fig. 1), the reversible reaction catalyzed by pyruvate-formate lyase (PFL). PFL orthologs are ubiquitously distributed in bacteria, but PFL has only been found in archaea such as T. sibiricus (6), T. kodakaraensis KOD1 (14), D. desulfuricans SSM1 (29), Archaeoglobus fulgidus (39), and Methanobacterium thermoautotrophicum (40). Therefore, it is conceivable that assimilation of pyruvate may commence with formate and acetyl-CoA. This reaction is the first committed step of anaerobic glucose metabolism and is important in Escherichia coli and other facultative anaerobes (41,42). Genomic analysis showed that no gene was annotated as PFL, suggesting that protein(s) with PFL activity in T. onnurineus NA1 may not be homologous to the PFLs of other archaea. Further study is necessary to identify PFL-like proteins in T. onnurineus NA1.
PFL is activated by pyruvate-formate lyase-activating enzyme (PFL-AE), which generates the stable glycyl radical at Gly-734 of PFL (43,44). PFL-AE is a member of the novel enzyme family that utilizes S-adenosylmethionine to initiate FIG. 1. A putative scheme of the cellular metabolism in the hydrogen producing arachaeon T. onnurineus NA1. A, glycolysis and gluconeogenesis. B, pyruvate degradation. C, pentose phosphate synthesis and carbon dioxide fixation. D, pseudo-TCA cycle. E, amino acid degradation and aldehyde metabolism. F, propionate utilization. Dashed lines or faint italic letters represent pathways or enzymes not yet experimentally validated in Thermococcales species. Each protein with its predicted function in ion or solute transport, or H 2 evolution or consumption, is illustrated on the membrane. Blue, red, and gray indicate proteins expressed prominently during growth on formate, CO, and starch, respectively. Major metabolic proteins deduced from the annotatable coding sequences (CDSs) were categorized by their induction on different substrates: formate (blue), CO (red), and starch (black). Other enzymes of the Embden-Meyerhof-Parnas pathway are designated by purple text, although they were mainly induced during growth on starch. Enzymes involved in sulfur reduction or ATP synthesis are shown in green, irrespective of their substrate specificity. Most proteins identified in T. onnurineus NA1 are shared by other Thermococcales: NPSOR (TON_0129), NADH:polysulfide oxidoreductase; FNR (TON_0056), ferredoxin-NADP(ϩ) reductase; braC-4 (TON_0162), branched chain amino acid transporter; CbiOQ (TON_0243), ABC-type cobalt transport system; FepBCD (TON_0296), ABC-type Fe 3ϩ -siderophore transporter; FbpABC (TON_0983 and TON_0985), ABC-type Fe 3ϩ transporter; Pst (TON_1464 and TON_1551), hypothetical phosphate transporter; NhaC (TON_1517), Na ϩ /H ϩ antiporter; MnhB (TON_1577), multisubunit Na ϩ /H ϩ antiporter; AppABC/OppBCDF (TON_1764 and TON_1768), ABCtype dipeptide/oligopeptide transporter; MalEFGK (TON_1791 and TON_1795), ABC-type maltodextrin transporter; AAT, additional amino acid ABC transporter; Fdred, reduced ferredoxin; Fdox, oxidized ferredoxin. The abbreviations of all of the identified proteins are presented separately in Table II and supplemental Table S1.

Analysis of One-carbon Metabolism in T. onnurineus NA1
radical catalysis (43,44). TON_0415, which is a PFL-AE homolog belonging to the radical S-adenosylmethionine superfamily (EC 1.97.1.4), was strongly up-regulated during growth on starch, suggesting that TON_0415 contributes to glycolytic flux in the Embden-Meyerhof-Parnas pathway ( Fig. 1 and Table II). This prompted us to inquire whether other isozymes of PFL-AE are present in T. onnurineus NA1. We pooled all of the identified proteins, investigated the domain assignments of the radical S-adenosylmethionine enzyme family in T. onnurineus NA1, and identified several possible PFL-AE candidates (TON_0789, TON_1269, TON_1353, TON_1410, and TON_1800). Of these, TON_1269, TON_1353, and TON_1410 were strongly up-regulated during growth on formate (Table II). Pyruvate formed by the reverse reaction of PFL can then be converted into various biosynthetic intermediates by the action of enzymes involved in gluconeogenesis and the PPP. As shown in Fig. 1 and Table II, four enzymes involved in glycolysis were strongly up-regulated during growth on starch (GLK, TON_0246; PGI, TON_0247; PFK, TON_1012; GAPOR, TON_1498), whereas three enzymes involved in key steps of gluconeogenesis were robustly up-regulated in formategrown cultures (glyceraldehyde-3-phosphate dehydrogenase, TON_0639; 2-phosphoglycerate kinase, TON_0742; FBPase, TON_1497). These data suggest that three gluconeogenic enzymes generate biosynthetic intermediates during growth on formate. More than seven enzymes of the PPP were identified (Fig. 1). Four related enzymes, ribose-5-phosphate isomerase (TON_0168), adenine phosphoribosyltransferase (TON_0120), AMP phosphorylase homolog DeoA (TON_1062), and HPS/PHI (TON_0336), were induced in response to formate ( Fig. 1 and Table II), the exceptions being PRPPS (TON_1403) and an RBPI homolog (TON_1296). The relative levels of these enzymes among formate-, CO-, and starch-grown cells suggest that biosynthesis of PPP products is favored during growth on formate. These data collectively show that formate can be utilized as a source of organic carbon via gluconeogenesis and the PPP. However, HPS/PHI (TON_0336) was detected at a 2-fold lower abundance in formate-grown cells compared with starch-grown cells (Table  II), indicating that HPS/PHI (TON_0336) expression during growth on formate can be suppressed when an energetically superior carbon source is present. Starch is a preferred substrate for growth so long as low levels of formate are provided. This process is reversible under certain conditions (e.g. glucose depletion) and may thus facilitate formate utilization for anabolic flux. Therefore, T. onnurineus NA1 can survive heterotrophically in environments with limited organic carbon, such as around hydrothermal vents. In addition, our observation that specific metabolic enzymes were elevated in response to formate provides the first evidence of the presence of a substrate-induced mechanism governing cellular metabolic flux in T. onnurineus NA1.
Pseudo Tricarboxylic Acid Cycle and Propionate Utilization-The T. onnurineus NA1 proteome indicated the pres-ence of another metabolic pathway, the reductive tricarboxylic acid cycle. Two key enzymes are required to reverse the TCA cycle (reductive TCA cycle) in archaea: KGOR and ATPdependent citrate lyase. Although ATP-dependent citrate lyase has not yet been detected in the T. onnurineus NA1 proteome, another key enzyme, KGOR (TON_0583-4 and TON_0586 -7) was detected in this study. KGOR catalyzes the reductive carboxylation of succinyl-CoA to form 2-oxoglutarate using reduced ferredoxin. However, the reductive TCA cycle of T. onnurineus NA1 is incomplete because of the absence of malate dehydrogenase, fumarate hydratase, and fumarate reductase. Similar results were reported in genomic analyses of other Thermococcales (14,15). Therefore, these cycles are called pseudo-TCA cycles. Despite the incompleteness of the reductive TCA cycle, our proteomic analysis showed significant up-regulation of proteins involved in COdependent anabolic flux for assimilation of CO 2 via a reductive TCA cycle (29,(45)(46)(47), including pyruvate ferredoxin oxidoreductase (POR, TON_1480 -82), phosphoenolpyruvae synthase (TON_0311), phosphoenolpyruvate carboxykinase (TON_0192), and KGOR (TON_0583-4 and TON_0586 -7) as shown in Fig. 1. This CO-dependent anabolic flux appears to be important for T. onnurineus NA1, because it could assimilate CO 2 into cellular carbon via the reductive TCA cycle. In contrast, SCS (TON_1665) and PC (TON_0904), enzymes of the pseudo-TCA cycle, were abundantly expressed in cells grown on starch ( Fig. 1 and Table II).
In addition to key enzymes in the pseudo-TCA cycle, we found up-regulation of type III ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO, TON_1234) in CO-grown cells, as shown in Table II. RuBisCO activities were detected in T. kodakaraensis KOD1 (48,49), A. fulgidus (50), and several other methanogenic archaea (51,52). Intriguingly, type III RuBisCO from T. kodakaraensis KOD1 has been shown to complement a RuBisCO-deficient strain of the purple nonsulfur bacterium Rhodopseudomonas palustris and therefore may allow the Calvin-Benson-Bassham pathway mode (49). As recently described in other Thermococcales (14,15), the CO 2 produced from catabolism or CO oxidation may be a suitable substrate for the type III RuBisCO (TON_1234), which together with DeoA (TON_1062) and RBPI (TON_1296) replenish the 3-phosphoglycerate used in glycolysis and gluconeogenesis (Fig. 1). This CO-dependent functional plasticity of RuBisCO may be important for the survival and persistence of T. onnurineus NA1 living heterotrophically around hydrothermal vents where organic carbon is limited. Thus, this finding provides additional evidence that T. onnurineus NA1 can grow on CO as the sole carbon source.
Amino Acid Degradation and Aldehyde Metabolism-In T. onnurineus NA1, the metabolic fate of aldehydes produced during anaerobic fermentation is unclear. Interestingly, one unique alcohol dehydrogenase (ADH) and other key enzymes required to metabolize aldehydes into ethanol or carboxylic acids were detected in our proteomic analysis ( Fig. 1 and Table II). In most anaerobic archaea (14,15,53), pyruvate is oxidatively decarboxylated into acetyl-CoA and CO 2 via POR (TON_1480 -82), which is then converted into acetate as an end product by the reverse reaction of acetyl-CoA synthetase (acetyl-CoA synthetase I, TON_1097 and TON_0327; and acetyl-CoA synthetase II, TON_1002 and TON_1313). Hence, a key question is how 2-oxoacids such as pyruvate are converted into ethanol in this bacterium. For hyperthermophilic archaea, an understanding of anaerobic fermentation will facilitate the development of a suitable enzyme system for producing ethanol at high temperatures with high efficiency. In addition to POR, three other types of 2-oxoacid ferredoxin oxidoreductase (KGOR, TON_0583-84 and TON_0586 -87; IOR, TON_1311-12; VOR, TON_1477-79) metabolize the 2-oxoacids generated by transamination of amino acids to their corresponding acyl-CoA derivatives (Fig. 1). Alternatively, 2-oxoacids may also be metabolized by 2-oxoacid: ferredoxin oxidoreductase (TON_0860 -61) into their corresponding aldehydes (53,54). Thus, aldehydes would be subsequently oxidized into their corresponding carboxylic acids by two different tungsten-containing aldehyde ferredoxin oxidoreductases, aor-2-tungsten-containing aldehyde ferredoxin oxidoreductase (TON_1354), and formaldehyde:ferredoxin oxidoreductase (TON_0749). Therefore, we propose that 2-oxoacid:ferredoxin oxidoreductase (TON_0860 -61) is the bridge linking the conversion of 2-oxoacid to its corresponding acyl-CoA or aldehyde derivatives in T. onnurineus NA1. The known steps of 2-oxoacid degradation are shown in Fig. 1 (solid arrows). Notably, these enzymes were robustly up-regulated in CO-grown cells (Table II), indicating that they could function in aldehyde metabolism during growth on CO. From these data, we infer that under CO-added conditions, the aldehyde produced from 2-oxoacids is oxidized into carboxylic acids by aldehyde ferredoxin oxidoreductase or formaldehyde:ferredoxin oxidoreductase, wherein the excess reductant (reduced ferredoxin) is utilized for H 2 S production, as proposed previously for other Thermococcales (54 -56). ADH (TON_0544), which catalyzes aldehydes into ethanol, was also detected. However, this was more abundantly expressed in starch-grown cells than in CO-grown cells (Table II). In the case of T. kodakaraensis KOD1, when pyruvate or starch was added to the medium, the cells predominantly produced H 2 rather than H 2 S (57). Under S 0 -limiting conditions, H 2 and alcohol production increase dramatically, as does activation of ADH (54). Therefore, we propose that a fermentation process shift may occur under S 0 -limited and starch excess conditions, wherein both H 2 and alcohol production increase but the amount of acids and H 2 S decrease, thereby serving to dispose of excess reductant.
Formate-dependent Hydrogenase and Hydrogenase Maturation System-Our proteomic analysis revealed a number of hydrogen clusters apparently associated with formate metabolism, consistent with previous results (9). These include two tripartite gene clusters containing formate dehydrogenases (Fdh): fdh1-mfh1-mnh1 (Hyg4-I) and fdh2-mfh2-mnh2 (Hyg4-III). In the comparative proteomic analysis, TON_1559 -1577, belonging to the Hyg4-III cluster, was substantially abundant in formate-grown cells relative to other culture conditions (Table II and Fig. 2). In particular, Fdh2 (TON_1563) was highly induced on formate compared with both starch (Ͼ12.50-fold) and CO (Ͼ17.46-fold). Of these Hyg4 III components, the most strongly induced member was mfh2 (TON_1569 -1572), which was also abundantly expressed in gene expression analysis (9). In a recent study, the hydrogenase encoded by mfh2 was proven to be responsible for growth and H 2 production when formate was supplied as a substrate (9). Our proteomic results support the role of mfh2 in formate metabolism. Furthermore, the up-regulation of an Na ϩ /H ϩ antiporter (TON_1577) in cells grown on formate is consistent with the previous report (9). Other Hyg4-III hydrogenases that were up-regulated under formate culture conditions were the F420 hydrogenases Frh TON_1559 and Frh TON_1561 and formate dehydrogenase FdhD (TON_1562). These protein complexes (FdhD and Frh), known as formate hydrogen lyase systems, convert formate into CO 2 and H 2 . Therefore, during growth on formate, H 2 can be produced by either the Fdh/Frh or fdh2-mfh2 Hyg4-III systems, probably in equilibrium with formate concentration.
F420-reducing hydrogenase is found solely in methanogenic archaea (58 -60), where it oxidizes H 2 to provide reduced coenzyme F420 (F420H 2 ), which acts as the electron donor for CO 2 reduction via the methane metabolic pathway. Recently, this protein was also discovered in Thermococcus strains such as T. onnurineus NA1 and T. gammatolerans (8,15). Although further examination of its role in T. onnurineus NA1 is necessary, generation of H 2 by F420 hydrogenase may be required to maintain redox balance during formatedependent growth. Coenzyme F420 may act as a low redox potential electron carrier, as proposed for methanogenic metabolism (58 -61). In methane-producing cells, the intracellular concentration ratios of reduced and oxidized coenzyme F420 are in equilibrium with the H 2 concentration in the medium (61). In view of this, formate metabolism may be modulated by the cellular concentration of H 2 , F420, and formate.
On the other hand, most Hyg4-I cluster proteins (TON_ 0266 -0282) were differentially expressed during growth on CO and starch as well as formate, indicating that this cluster has a minor (or no) relationship with formate metabolism. This result is consistent with a previous gene expression analysis in which genes encoding the mfh1 cluster (TON_0273-78) were slightly down-regulated during growth on formate (9). Overall, the proteomic and gene expression data for the formate metabolic pathway studied correspond relatively well, indicating that this pathway is not subject to additional translational control. Along with the hydrogenase (Hyg4-I) cluster, we identified the known auxiliary proteins Hyc I (TON_0263), Hyp F (TON_0286), and Hyp E (TON_0287), which are involved in hydrogenase maturation (62,63). Of these, Hyc I (TON_0263) was prominently up-regulated in cells grown on formate (Table II and Fig. 2). This finding is closely correlated with the fact that, in E. coli, expression of both the fdh gene and the hyc operon depends on formate and acidification of the growth medium (64). Also, the non-energy-conserving formate hydrogen lyase complex (Fdh-Mhy) is generally involved in formate dissimilation, thereby preventing acidification of the cytoplasm. Although formate dehydrogenase (FdhA, TON_0281) was not detected in cells grown on any substrates, the expression of Mhy1 (TON_0276) was up-regulated in cells grown on CO (Table II and Fig. 2). Presumably, the physiological functions of the FdhA-Mhy1 complex are linked to formate removal likely generated during amino acid catabolism, as proposed in Thermococcus litoralis (65).
CO-responsive Hydrogenase and Other Membrane-bound Hydrogenases-Previous studies revealed the presence of CO-dependent hydrogenase gene clusters in T. onnurineus NA1 (8,18). In accordance with these studies, our proteomic analysis showed marked expression of hydrogenase gene clusters belonging to Hyg4_II during growth on CO (Table II and Fig. 2). In particular, both the CO dehydrogenase (CODH, TON_1017-19) and Mch hydrogenase (TON_1023-24) of the Hyg4_II cluster were substantially elevated in CO-grown cells, suggesting that both hydrogenases function during carboxydotrophic growth. This result is consistent with our previous microarray studies that showed that CODH and a large Mch subunit (TON_1023) were significantly up-regulated in COgrown cells (data not shown). According to a previous report (66), most carboxydotrophic hydrogenogenic bacteria share a similar gene organization composed of CODH and an energyconverting hydrogenase. Although we detected no evidence of a functional energy-converting hydrogenase-like protein in T. onnurineus NA1, it is possible that CO-induced Mch hydrogenase clusters putatively encode a similar protein. Indeed, T. onnurineus NA1 encodes a putative membranebound hydrogenase (TON_1021) with NADH:quinone oxidoreductase activity, in addition to other Mch hydrogenase components (TON_1022-24). Even though it is membranebound, TON_1021 was not detected in our proteomic analysis; it is possible that this Mch hydrogenase may play a key role in carboxydotrophic hydrogenogenic metabolism (66,67).
Also, in CO-grown cells, the marked expression of monofunctional CODH (TON_1018), homologous to CooS and CooF (TON_1017), indicates the possibility that T. onnurineus NA1 can evolve H 2 during growth on CO, as displayed by Rhodospirillum rubrum (68) and Carboxydothermus hydrogenoformans (69,70). In addition, it has been reported that energy-converting hydrogenase-type hydrogenases couple the formation of H 2 to the generation of the proton motive force that drives ATP synthesis (66,71). Likewise, coupling of CO oxidation to proton translocation in T. onnurineus NA1 may enable it to use CO as a sole energy source and generate H 2 .
Sulfhydrogenase-I was also up-regulated in cells grown on CO (Table II and Fig. 2). The CO-induced sulfhydrogenase-I cluster contains a number of genes that encode enzymes putatively involved in sulfur reduction and H 2 recycling to provide NADPH for biosynthesis (72,73). As shown in Table II and Fig. 2, cytoplasmic (Ni-Fe) hydrogenases (encoded by TON_0534 to TON_0536) homologous to the Hyh-I of T. kodakaraensis KOD1 were identified in the sulfhydrogenase-I cluster of T. onnurineus NA1. Like Hyh-I of T. kodakaraensis KOD1 (74), these enzymes may function as an H 2 uptake hydrogenase during H 2 production. In addition to sulfhydrogenase (TON_0537), a glutamate synthase ␤ subunit (TON_0542) similar to that coded by KOD1-gltA was also present in T. onnurineus NA1. Glutamate synthase shows homology to the sulfide dehydrogenase that is involved in S 0 reduction using NAD(P)H as the electron donor (75,76). Up-regulation of formate dehydrogenase (TON_0539) was also detected in CO-grown cells but may be attributable to re-oxidation of formate rather than to formate production, as described in A. fulgidus (77). Although ADH (TON_0544) most likely functions in aldehyde reduction, it may also provide reductants for S 0 reduction, as proposed in P. furiosus (78).
Furthermore, a membrane-bound oxidoreductase (Mbx, TON_0488) was also strongly up-regulated in CO-grown cells (Fig. 2). Mbx appears to be directly involved in oxidizing reduced ferredoxin to transfer electrons to elemental sulfur (79). However, expression of a membrane-bound hydrogenase (Mbh: TON_1593) was higher during growth on starch (Fig. 2), which is consistent with its attributed function in sugar metabolism as seen in other Thermococcales (5,57,63).
Conclusion-To the best of our knowledge, this is the most comprehensive quantitative analysis of T. onnurineus NA1, the first global survey of metabolic activities at the proteome level, and the first detailed proteomic investigation of anabolic reactions in hyperthermophilic archaea. We analyzed and compared the soluble proteomes of T. onnurineus NA1 cells grown on formate, CO, and starch using one-dimensional SDS-PAGE in combination with nano UPLC-MS E . Our data provide a solid basis for further studies of the metabolic characteristics of T. onnurineus NA1. One-carbon metabolism in the H 2 -producing hyperthermophilic archaeon T. onnurineus NA1 is an essential process that relies on at least two different one-carbon donor molecules: CO and formate. In particular, formate was metabolized via gluconeogenesis and the PPP into cellular carbon. On the other hand, CO 2 formed from CO oxidation by CODH and CO-induced hydrogenases is a suitable substrate for carbon assimilation by type III RuBisCO. In addition, our data suggest that pyruvate formed from acetyl-CoA and CO 2 by the reverse reaction of POR can enter the CO-dependent anabolic pathway, generating cellular carbon via the incomplete TCA cycle. Also, CO was utilized as an energy source for a CO oxidation-coupled electron transport mechanism that resulted in H 2 production. Our data provide novel insights into the poorly understood metabolism of one-carbon compounds in hyperthermophilic archaea. Furthermore, identified proteins with known sequence or motif/ domain homologies were classified into groups according to their known and putative biological functions. It is also intriguing to note that enzymes involved in energy metabolism were found as major proteins in our proteomic analysis. In addition to these key metabolic enzymes, we identified and quantified proteins involved in other cellular processes, including components of the translational machinery, amino acid metabolism, DNA metabolism, nucleotide metabolism, transcription, the stress response, and transporters. This study also provides a pool of potentially important proteins with unknown functions that are specific to hyperthermophilic archaea and adds to our understanding of their biological functions and metabolic adaptation in extreme environments. These proteins represent an attractive target for further directed genetic and biochemical analyses to shed light on the physiology of T. onnurineus NA1.