This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Siebers, B. Articles by Hensel, R. Search for Related Content PubMed PubMed Citation Articles by Siebers, B. Articles by Hensel, R.

Reconstruction of the Central Carbohydrate Metabolism of Thermoproteus tenax by Use of Genomic and Biochemical Data

Department of Microbiology, Universität Duisburg-Essen, 45117 Essen,¹ Epidauros Biotechnologie AG, 82347 Bernried,³ Max-Planck Institut für Entwicklungsbiologie, 72076 Tübingen, Germany,⁴ Sun Center of Excellence for Visual Genomics, Faculty of Medicine, University of Calgary, Calgary, Alberta, Canada T2N 4N1²

Received 10 June 2003/ Accepted 13 November 2003

	ABSTRACT

Top Abstract Introduction Materials and Methods Results and Discussion References

 
The hyperthermophilic, facultatively heterotrophic crenarchaeum Thermoproteus tenax was analyzed using a low-coverage shotgun-sequencing approach. A total of 1.81 Mbp (representing 98.5% of the total genome), with an average gap size of 100 bp and 5.3-fold coverage, are reported, giving insights into the genome of T. tenax. Genome analysis and biochemical studies enabled us to reconstruct its central carbohydrate metabolism. T. tenax uses a variant of the reversible Embden-Meyerhof-Parnas (EMP) pathway and two different variants of the Entner-Doudoroff (ED) pathway (a nonphosphorylative variant and a semiphosphorylative variant) for carbohydrate catabolism. For the EMP pathway some new, unexpected enzymes were identified. The semiphosphorylative ED pathway, hitherto supposed to be active only in halophiles, is found in T. tenax. No evidence for a functional pentose phosphate pathway, which is essential for the generation of pentoses and NADPH for anabolic purposes in bacteria and eucarya, is found in T. tenax. Most genes involved in the reversible citric acid cycle were identified, suggesting the presence of a functional oxidative cycle under heterotrophic growth conditions and a reductive cycle for CO₂ fixation under autotrophic growth conditions. Almost all genes necessary for glycogen and trehalose metabolism were identified in the T. tenax genome.

	INTRODUCTION

Top Abstract Introduction Materials and Methods Results and Discussion References

 
Archaea were recognized as a distinct phylogenetic group more than 25 years ago (75). Their importance as the third major evolutionary line is well established, but our knowledge of their physiological capabilities remains limited. Central metabolic pathways within these organisms are far from fully understood (43, 47, 72).

Thermoproteus tenax was the first hyperthermophilic archaeum described (76). It is able to grow chemolithoautotrophically on H₂, CO₂, and S^o as well as chemoorganoheterotrophically in the presence of S^o and various organic substrates such as glucose, starch, amylose, glycerate, glycerol, ethanol, and malate (12, 76). Physiological and biochemical studies revealed T. tenax as a physiologically versatile organism with numerous archaeon-specific metabolic capabilities, regulation, and thermoadaptive traits.

Comparative studies of carbohydrate metabolism in hyperthermophilic archaea indicate that sugars are generally metabolized by variants of the Entner-Doudoroff (ED) and Embden-Meyerhof-Parnas (EMP) pathways. The so-called nonphosphorylative ED pathway (phosphorylation takes place only at the stage of glycerate) is the only pathway that has been identified for sugar degradation in the aerobes Sulfolobus solfataricus (9) and Thermoplasma acidophilum (7), whereas the anaerobes Pyrococcus furiosus (27, 28, 31, 66, 69), Thermococcus spp. (31, 44, 50), Desulfurococcus amylolyticus (20), and Archaeoglobus fulgidus (33) use modified versions of the EMP pathway. In contrast to findings for other hyperthermophilic archaea, T. tenax uses both variants—the ED and EMP pathways—for glucose metabolism, as shown by in vitro studies identifying specific intermediates and enzyme activities in cell extracts (49, 55). The utilization of both pathways was furthermore shown by in vivo ¹³C-labeling experiments with growing cells (57) and with cell suspensions (50). Several enzymes and their coding genes involved in the degradative pathways, predominantly in the EMP variant, have already been characterized (4, 5, 48, 54, 56, 57), but the brunt of core carbon biochemistry in T. tenax is still unknown.

The main difference between the T. tenax EMP variant and the classical version of this pathway concerns (i) the presence of a hexokinase (HK) with no notable regulatory properties (11), (ii) the replacement of the antagonistic enzyme couple ATP-dependent PFK and fructose bisphosphatase by a bidirectional PP_i-dependent enzyme (56), (iii) the presence of two different glyceraldehyde-3-phosphate (GAP) dehydrogenases (GAPDHs), differing in phosphate dependence, reversibility of the catalyzed reaction, and allosteric properties (4, 5, 21, 40) and (iv) a pyruvate kinase (PK) with reduced allosteric potential (48). Enzyme as well as transcript studies with T. tenax (late-log-phase cells) indicate that regulation takes place on the protein and gene levels. Under heterotrophic conditions, the catabolic flux seems to be forced by the preferred expression of the PP_i-dependent phosphofructokinase (PP_i-PFK) and fructose-bisphosphate aldolase (FBPA), which form an operon in T. tenax (fbp-pfp operon) (34, 54, 56), and of the PK (48). A strong influence on the catabolic flux is also exerted by the allosterically regulated nonphosphorylating GAPDH (GAPN), which catalyze the unidirectional oxidation of GAP to 3-phosphoglycerate. In the presence of activators such as glucose-1-phosphate, fructose-6-phosphate, AMP, and ADP, i.e., under conditions which are characterized by the availability of storage carbohydrates (e.g., glycogen) and/or a low-energy charge of the cell, the catabolic carbon flux is most effective. In contrast, the classical, phosphorylating NADP⁺-dependent GAPDH and phosphoglycerate kinase (PGK) are of predominant importance for anabolism rather than for catabolism, as indicated by the enzymatic features of the NADP⁺-dependent GAPDH, by higher transcript amounts of both genes (gap and pgk), and by higher-level enzyme activities in autotrophically grown cells (5). Phosphoenolpyruvate synthetase (PEPS) activity was detected only in extracts of cells grown under autotrophic conditions (53); thus, the PEPS seems to exert a key function in driving the carbon flux into the anabolic direction.

Enzymatic studies with cell extracts (49, 55) and analyses of characteristic intermediates (55) identified the nonphosphorylative version of the ED pathway, which had been described originally for S. solfataricus (9) and T. acidophilum (7). In contrast to findings for the haloarchaeal semiphosphorylative pathway (with phosphorylation at the level of 2-keto-3-deoxy-gluconate [KDG]), in this version of the ED pathway phosphorylation takes place only at the level of glycerate. Glyceraldehyde generated by the KDG aldolase (KDGA) is therefore the characteristic intermediate of the nonphosphorylative version. So far, only glucose dehydrogenase (GDH), the first enzyme of the nonphosphorylative ED pathway of T. tenax, has been purified and characterized (57).

Because of its essential function for providing pentoses and NADPH for anabolic purposes, the pentose phosphate pathway (PPP) has been assumed for all three domains of life. Genome analysis confirmed by biochemical studies indicates that glucose-6-phosphate dehydrogenase seems to be generally absent from methanoarchaea and hyperthermophilic archaea (51, 52); therefore, the conventional PPP does not seem to be operative in these organisms. Also, in halophiles—although activity of glucose-6-phosphate dehydrogenase and 6-phosphogluconate dehydrogenase was detected in Halococcus saccharolyticus (23) and gene homologs of these enzymes were found in the genome of Halobacterium sp. strain NRC1—PPP is not involved in sugar catabolism, as indicated by ¹³C nuclear magnetic resonance studies of H. saccharolyticus. Possibly, the fragmentary PPP serves just for the synthesis of pentoses in haloarchaea.

The catabolic fate of pyruvate was monitored in T. tenax by fermentation analyses and enzyme measurements in crude extracts, which showed that the final oxidation of carbohydrates occurs via an oxidative citric acid cycle (CAC) linked to the glycolytic pathways by a pyruvate synthase (POR) (pyruvate:ferredoxin oxidoreductase) (49).

No information is available yet about the reaction sequence of CO₂ fixation in T. tenax. However, there is some evidence indicating that the CO₂ fixation occurs via the reductive CAC, as demonstrated for T. neutrophilus, a close relative of T. tenax (1). Most of the enzymes of the CAC work reversibly and could therefore be engaged in both directions. Only three counteracting enzyme pairs are thought to determine the oxidative or reductive direction of the cycle: (i) citrate synthase/citrate lyase, (ii) 2-oxoglutarate dehydrogenase/2-oxoglutarate synthase (2-oxoglutarate-ferredoxin oxidoreductase [KOR]), and (iii) succinate dehydrogenase/fumarate reductase. Experimental data (49) give no support for the presence of alternative CO₂ fixation mechanisms such as the Calvin cycle via ribulose-1,5-bisphosphate carboxylase or the reductive acetyl coenzyme A (acetyl-CoA) pathway via carbon monoxide dehydrogenase (Wood-Ljungdahl pathway).

Two additional carbohydrates—glycogen and trehalose—were identified in T. tenax, thus raising questions about their metabolism and function. Glycogen is a branched polymer of -1,4-linked and -1,6-linked glucosyl moieties and represents an osmotic inactive carbon storage compound in all three domains of life. Although the presence of glycogen in T. tenax and other archaea (members of Thermoproteales and Sulfolobales) was demonstrated previously (32), not much is known about glycogen metabolism in archaea. Trehalose, a glucose disaccharide [-D-glucopyranosyl-(1,1)--D-glucopyranoside], is found in members of all three major phylogenetic domains. Trehalose plays an important role as a compatible solute, being involved in stress response to high-level osmolarity (Escherichia coli) and heat (Saccharomyces cerevisiae) (16, 39). Although trehalose has been identified in several archaea (e.g., T. acidophilum, Methanothermobacter thermautotrophicus, and T. tenax) (37, 39), the role of trehalose in archaea is still unknown. In T. tenax, trehalose was identified as the exclusive solute in concentrations of up to 0.3 µmol/mg of protein (37).

These insights into the physiological capabilities of T. tenax motivated us to study in more detail the biochemical and genetic mechanisms that select the various pathways, regulate the carbon flux through them, and govern the synthesis and breakdown of low- and high-molecular-mass carbohydrate compounds. For a better understanding of the central metabolism in T. tenax, as well as in archaea in general, we aimed at a representative genome sequence analysis of T. tenax to use the genetic information for reconstruction of the carbohydrate metabolism by assigning genes by sequence similarity and/or by the function of their recombinant gene products. So far, only four complete crenarchaeal genome sequences have been published: those for Aeropyrum pernix (26), Pyrobaculum aerophilum (13), and two species of Sulfolobus (25, 52). Here we report on the reconstruction of the T. tenax-specific EMP and ED pathways and the reversible CAC as well as on glycogen and trehalose metabolism.

	MATERIALS AND METHODS

Top Abstract Introduction Materials and Methods Results and Discussion References

 
Strains and growth conditions. Mass cultures of T. tenax Kra1 (DSM 2078) (12, 76) were grown under autotrophic conditions as described previously (5). E. coli strains DH5 (Invitrogen), BL21(DE3), and BL21-CodonPlus(DE3)-RIL (Stratagene) for use in cloning and expression studies were grown under standard conditions (46) following the instructions of the manufacturers.

Library construction, shotgun sequencing, and primer walking. Genomic DNA was prepared from T. tenax cells as described previously (48). For construction of a whole-genome shotgun library, genomic DNA from T. tenax was sheared with a HydroShear (GeneMachines, San Carlos, Calif.) to yield fragments with estimated sizes of 2.5 to 5.0 kbp. The fragmented DNA was separated in a low-melting-point agarose gel, and the dominant fraction of the desired size range was isolated. The fragments were purified with a QIAquick gel extraction kit (Qiagen, Hilden, Germany) and cloned using a TOPO shotgun subcloning kit (Invitrogen Corp., Carlsbad, Calif.). Positive clones were identified by their ampicillin resistance in combination with their white colony phenotype in a X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside) screening. The isolated vectors served as templates in sequencing reactions using the primers M13uni(-43) (5'-agggttttcccagtcacgacgtt-3') and M13rev(-49) (5'-gagcggataacaatttcacacagg-3'). The sequence data were collected using an ABI Prism DNA System 377 apparatus and an ABI Prism 3700 DNA analyzer and processed and assembled with Phred/Phrap/Consed software (University of Washington). Gaps were closed where needed by direct sequencing using chromosomal DNA as a template.

Automated analysis and annotation using the MAGPIE system. The genomic sequence was analyzed using Multipurpose Automated Genome Project Investigation Environment (MAGPIE) automated genome analysis and annotation software (15, 17). The function of open reading frames (ORFs) involved in metabolic pathways was verified manually, and gene ontology categories were attached to the functional assignments.

Cloning of genes and expression in E. coli. To prove the assignment of gene homologs in the T. tenax genome, the respective genes were cloned into a pET vector system via PCR mutagenesis. Expression of the gene products in E. coli BL21(DE3) and BL21-CodonPlus(DE3)-RIL was performed following the instructions of the manufacturer (Stratagene). For expression, the glucose-6-phosphate isomerase (PGI) gene (pgi) of T. tenax was cloned into vector pET-11c through the use of the two new restriction sites (NdeI and BamHI) introduced by PCR mutagenesis with the following primer set: PGI-NdeI (5'-TCCGAAGATAGCTGCATATGATAG-3' [sense]) (the mutations are shown in boldface characters and the resulting restriction sites are underlined) and PGI-BamHI (5'-TCGGCCAGGATAGGATCCTTGATC-3' [antisense]). The following primer set—devised on the basis of N- and C-terminal partial sequence information—was designed to amplify the glycogen phosphorylase (glgP) gene in pET-24a: GLGP-AseI-f (5'-GCCCCCCCGATTAATGATAGTGT-3' [sense]) and GLGP-HindIII-rev (5'-AGATCAAGCTTTTGGCCGAGCTCA-3' [antisense]). Both enzymes were expressed in BL21-CodonPlus(DE3)-RIL.

Isolation of recombinant hyperthermophilic enzymes. Recombinant E. coli cells were suspended in buffer (e.g., 100 mM HEPES-KOH [pH 7.5] or 100 mM TRIS-HCl [pH 7.0 at 70 to 90°C, depending on the temperature for heat precipitation]) containing 7.5 mM dithiothreitol (DTT) or 10 mM ß-mercaptoethanol and passed three times through a French pressure cell at 150 MPa. Cell debris and unbroken cells were removed by centrifugation (85,000 x g for 30 min at 4°C). The resulting crude extract was heat precipitated at 70 to 90°C for 20 to 30 min, depending on the stability of the respective enzyme. The fraction remaining after heat precipitation was centrifuged again (20,000 x g for 20 min at 4°C) and dialyzed overnight at 4°C. In general, after heat treatment only a few faint contaminant E. coli proteins were detectable by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and Coomassie blue staining.

For the recombinant PGI, 0.5 g of E. coli cells was resuspended in 3 ml of 100 mM TRIS-HCl (pH 7.0 at 90°C) containing 10 mM ß-mercaptoethanol. The heat precipitation was performed for 20 min at 90°C, and the protein was dialyzed overnight at 4°C with 1.6 liters of 100 mM TRIS-HCl (pH 7.0 at 55°C) containing 10 mM ß-mercaptoethanol. For the recombinant glycogen phosphorylase, 10 g of E. coli cells was resuspended in 30 ml of 100 mM HEPES-KOH (pH 7.5) containing 7.5 mM DTT and the heat precipitation procedure was performed for 30 min at 90°C. The extract was dialyzed overnight at 4°C with 2 liters of 50 mM HEPES-KOH (pH 7.5) containing 7.5 mM DTT, and the dialyzed fraction was applied on Q-Sepharose Fast-Flow (Amersham Pharmacia Biotech) equilibrated in the same buffer and eluted with a linear salt gradient of 0 to 600 mM KCl. Fractions containing the homogeneous enzyme solution were pooled, dialyzed, and concentrated.

Enzyme assays. PGI activity was determined in the anabolic direction (glucose 6-phosphate formation) at 50°C in a coupled assay with glucose-6-phosphate dehydrogenase (EC 1.1.1.49) from baker's yeast as an auxiliary enzyme. The assay (total volume, 0.5 ml) was performed at 50°C in 100 mM TRIS-HCl (pH 7.0) in the presence of 1 mM NADP⁺, 1 U of glucose-6-phosphate dehydrogenase, 0.05 to 7 mM fructose 6-phosphate, and 4.5 µg of PGI after heat precipitation. Enzyme activity in the catabolic direction (fructose 6-phosphate formation) was assayed at 50°C in a coupled assay using ATP-dependent phosphofructokinase (EC 2.7.1.11; Bacillus stearothermophilus), fructose-1,6-bisphosphate aldolase (EC 4.1.2.13; rabbit muscle), triosephosphate isomerase (TIM) (EC 5.3.1.1; rabbit muscle), and glycerol-3-phosphate dehydrogenase (EC 1.1.1.8; rabbit muscle) as auxiliary enzymes. The assay (total volume, 0.5 ml) was performed at 50°C in 100 mM TRIS-HCl (pH 7.0) in the presence of 1 mM NADH, 3 mM ATP, 5 mM MgCl₂, 2 U of phosphofructokinase, 1 U of fructose-1,6-bisphosphate aldolase, 50 U of TIM, 9 U of glycerol-3-phosphate dehydrogenase, 0.1 to 14 mM glucose 6-phosphate, and 0.9 µg of recombinant PGI after heat precipitation.

Glycogen phosphorylase activity was assayed in a coupled enzyme assay at 50°C with phosphoglucomutase (PGM) (EC 5.4.2.2; rabbit muscle) and glucose-6-phosphate dehydrogenase (EC 1.1.1.49; baker's yeast) as auxiliary enzymes. The standard assay (volume, 1 ml) was performed at 50°C in 100 mM TRIS-HCl (pH 7.5) containing 10 mM KH₂PO₄ (P_i), 1 mM NADP⁺, 10 mM MgCl₂, 3 U of PGM, 3 U of glucose-6-phosphate dehydrogenase, 0.1% (wt/vol) glycogen (bovine liver), and 5 to 15 µg of protein. For the different enzyme assays, the change in absorbance was monitored at 366 nm ( _{NADH 50°C} = 3.36 mM^-1 cm^-1; _{NADPH 50°C} = 3.43 mM^-1 cm^-1).

Nucleotide sequence accession numbers Gene homologs assigned in this paper have been submitted to the EMBL database, and accession numbers are given in the tables and text.

	RESULTS AND DISCUSSION

Top Abstract Introduction Materials and Methods Results and Discussion References

 
Genome sequence features. The genome assembly used for the annotation of the genes described in this paper consists of 19,850 sequence reads assembled in 53 contigs more than 2 kbp in length; the largest contig is 110 kbp in length. Contigs smaller than 2 kbp were excluded from the analysis. The total length of the contigs adds up to 1.81 Mbp, which represents about 98.5% of the estimated genome size of 1.84 Mbp. The average gap size is estimated to be about 100 bp. With sequence coverage at 5.3-fold, the sequence quality of the trimmed contigs is approximately 15 errors in 10,000 bp. The average G+C content of the genome is calculated to be 54%. The individual contigs were subjected to an automated analysis with MAGPIE software. The ORFs were identified on the basis of stop codons and their upstream start codons. An ORF size of 100 amino acids was chosen as the minimum cutoff level. The automated analysis tool identified and annotated 3,913 ORFs, which were then surveyed for the presence of genes related to central carbohydrate metabolism. Given an average G+C level of 54%, this leads to an overprediction of ORFs; however, at least 2,000 to 2,500 of these should be coding ORFs (given the number of coding regions in other organisms of the same size).

The central carbohydrate metabolism of T. tenax. (i) Variant of the EMP pathway. (a) Assignment of genes involved in the reversible EMP variant. The following enzymes (and their coding genes) essential for a functional reversible EMP pathway have been identified and characterized to date: PP_i-PFK (56), FBPA (34, 54), GAPDH (5), PGK (5), TIM (H. Walden, G. L. Taylor, E. Lorentzen, E. Pohl, H. Lilie, A. Schramm, T. Knura, K. Stubbe, B. Tjaden, and R. Hensel, submitted for publication), PK (48), and PEPS (63; B. Tjaden, C. Dörr, B. Siebers, and R. Hensel, unpublished data). As an additional constituent of the catabolic version of the pathway, the nonphosphorylating GAPN (4, 5, 21), and its coding gene were identified.

Genes still missing from the T. tenax genome for a functional EMP pathway were assigned (Table 1; Fig. 1). In the case of HK (11), enolase (ENO), PGI, and phosphoglycerate mutase (PGAM), the sequence similarity of the encoding genes allowed an unequivocal assignment. The HK was originally assigned as glucokinase, and only the characterization of the enzyme after expression of the respective gene in E. coli revealed that it is a true HK, catalyzing the phosphorylation of different sugars such as glucose, fructose, and mannose (11).

View larger version (22K): [in this window] [in a new window]   FIG. 1. Variant of the EMP pathway in T. tenax. The enzymes corresponding to homologs identified in this study are shown in grey-shaded boxes; enzymes identified and characterized previously are shown in open boxes. Open arrows indicate irreversible reactions. For abbreviations of enzyme names, see Table 1.

T. tenax does not contain genetic information for glucose-6-phosphatase, suggesting either that glucose-6-phosphatase is absent or that T. tenax possesses a novel enzyme type that has not yet been recognized as such. In addition, gene homologs for ADP-dependent glucokinase or ADP-dependent phosphofructokinases present in other anaerobic, hyperthermophilic archaea (e.g., P. furiosus, A. fulgidus, and Thermococcus species) and for ATP-dependent PFK (Desulfurococcus species) are absent from the T. tenax genome, underlining the central function of the reversible nonallosteric PP_i-PFK.

Surprisingly, three ORFs in the T. tenax genome showed convincing similarities to genes encoding enzymes with functions obviously equivalent to those of the assigned EMP enzymes fructose-1,6-bisphosphatase (FBP), ferredoxin-dependent GAP oxidoreductase (GAPOR), and pyruvate, phosphate dikinase (PPDK) (partial sequence information). The identified FBP homolog is a member of the recently described archaeal FBP (type IV) and related enzymes of the inositol-monophosphatase family (COG 0483) and shows 28% identity to the characterized FBP of P. furiosus (70) and 29% identity to the inositol-1-monophosphatase (IMP) of M. jannaschii (61). More recently a type V FBP (COG 1980) was detected in Thermococcus kodakaraensis, and homologs were identified in almost all archaeal genomes (42, 72). A type V FBP which is supposed to be the general archaeal gluconeogenetic FBP (72) is absent from T. tenax. Thus, biochemical studies of the T. tenax enzyme have to be awaited to elucidate a possible gluconeogenetic role in parallel to that of the bifunctional nonallosteric PP_i-PFK. Ferredoxin-dependent GAPOR—like GAPN of T. tenax—catalyzes the nonphosphorylating, irreversible conversion of glyceraldehyde 3-phosphate to 3-phosphoglycerate and thus also fulfills a catabolic function. The conditions under which the expression of GAPOR is favored or discriminated against are yet unknown. The complete sequence information for the PPDK gene was obtained by the identification (via Southern hybridization) and sequencing of a genomic clone. The assignment of the PPDK homolog was confirmed by the functional properties of the encoded enzyme (63; Tjaden et al., unpublished data). The recombinant enzyme expressed in E. coli catalyzes the interconversion of PEP and pyruvate according to the following equation: ATP + P_i + pyruvate AMP + PP_i + PEP. The reversibility of the PPDK reaction implies that the direction of conversion depends strongly on the ratio of the substrates and cosubstrates on both sides of the equation; therefore, the enzyme fulfills the functions of both the anabolic PEPS and the catabolic PK.

In summary, the genome data indicate that the variant of the EMP pathway in T. tenax is much more complex than previously suspected. In contrast to results seen with the classical pathway, no control point is present in the phosphorylating part of the EMP variant, as indicated by the presence of a HK with apparently reduced allosteric potential and a reversible nonallosteric PP_i-dependent PFK. The major control point is shifted to the middle of the pathway at the level of GAP, where three different enzymes (GAPN, GAPDH, and GAPOR) are involved in the regulation of the carbohydrate flux. In addition, a second control point seems to emerge at the level of pyruvate/PEP conversion, which is catalyzed again by three different enzymes (PK, PEPS, and PPDK) (48, 63).

Other modifications of the EMP pathway are found in the anaerobic species P. furiosus (27, 28, 31, 66, 69), Thermococcus species (31, 44, 50), D. amylolyticus (20), and A. fulgidus (33). The glycolysis of D. amylolyticus differs from that of the classical version of the EMP pathway by the presence of a ferredoxin-dependent GAPOR replacing the "conventional" phosphate- and pyridine nucleotide-dependent GAPDH. A more extended variation is seen with P. furiosus and Thermococcus species (as well as with A. fulgidus), which additionally replace the classical ATP-dependent HK and phosphofructokinase with nonallosteric ADP-dependent glucokinases and phosphofructokinases. In addition, ADP-dependent phosphofructokinase activity was demonstrated in crude extracts from members of the orders Methanococcales (e.g., Methanococcus maripaludis) and Methanosarcinales (e.g., Methanosarcina mazei) (73) and a novel bifunctional ADP-dependent glucokinase/phosphofructokinase is described for M. jannaschii (45, 73). The absence of allosteric sugar kinases and the presence of GAPN and GAPOR enzyme activities in several archaea, as well as the identification of GAPN and GAPOR gene homologs—in addition to GAPDH—in many archaeal genomes (e.g., those of Pyrococcus spp.) (72), indicates that among all the different modifications of the EMP pathway in archaea, the regulation at the level of GAP might represent a unique feature of central carbon metabolism.

(b) Organization of the EMP genes. Analyses of gene organization and transcription studies reveal that several EMP genes are organized in operons, suggesting specific functional relationships between the genes or their products (Table 1) (see Fig. 3). Thus, as previously shown (54), the pfp gene coding for the PP_i-dependent PFK, which catalyzes the reversible phosphorylation of fructose 6-phosphate, forms an operon with the fba gene encoding FBA (fba-pfp operon). The coordinated transcription of these genes in one operon seems meaningful for T. tenax, because both enzymes catalyze not only consecutive but also reversible reactions and are therefore of equivalent levels of importance for both directions of the pathway. Also, the phosphorylating NADP⁺-dependent GAPDH gene (gap) and the PGK gene (pgk) are clustered in an operon (pgk-gap operon) (5). In both operons, the two genes overlap by 4 bp (the coding regions overlap by 1 bp) and the formation of bicistronic as well as monocistronic mRNAs was demonstrated. Northern analyses reveal that the HK gene (hxk) is part of an operon (orfX-hxk operon) (11) and that, again, both genes overlap by 4 bp. The ORF located upstream (orfX) codes for a protein of unknown function, and only one homolog was identified in P. aerophilum (AE009933) (11).

View larger version (109K): [in this window] [in a new window]   FIG. 3. (A) Organization of functionally related genes of the central carbohydrate metabolism in T. tenax. (Gene names are given in alphabetical order at the end of the legend; gene data are given in Table 1 to Table 6.) (B) Putative promoter structures of genes organized in operons. The nucleotide sequences (5' end) and upstream regions of genes are shown. The putative crenarchaeal promoter sequences (3, 59, 60) with BRE sites [crenarchaeal consensus sequence (A/G)N(A/T)AA(A/T)] and the TATA box [crenarchaeal consensus sequence (C/T)TTTTAAA] are in light and dark grey boxes, respectively. Putative Shine-Dalgarno sequences (GAGG) (29) are underlined, the putative start codon is shown in boldface characters, and the stop codon of preceding genes is indicated with double underlining. The operon structures shown were confirmed for the EMP genes by Northern analysis. All other clusters shown represent putative operons. Gene names and corresponding enzyme names: acn, aconitase; adh, alcohol dehydrogenase; acs, acetyl-CoA synthetase; act, acetyl-CoA transferase; amyA, -amylase; cis1, citrate synthase 1; citE, citrate lyase ß chain; fba, FBPA; frdAB, fumarate reductase and ß chains; gaa, glucan-1,4--glucosidase; gad, GAD; gap, GAPDH; glgA, glycogen synthase; glgP, (glycogen) phosphorylase; gltX, glutamyl tRNA synthetase; gt, glycosyltransferase; hp, hypothetical protein, hxk, HK; kdgA, KDP(G)A; kdgK, KDGK; mthfs, 5-formyltetrahydrofolate cyclo-ligase; oorABCD, oxoacid:ferredoxin oxidoreductase , ß, , and chains; orfV, orfX, orfY, and orfZ, hypothetical proteins; pncB, nicotinate phosphoribosyltransferase; pfp, PPi-PFK; pgk, PGK; ptp, putative transport protein; rfbA, glucose-1-phosphate thymidylyltransferase; rfbB, dTDP-glucose-4,6-dehydratase; rfbC, dTDP-4-dehydrorhamnose-3,5-epimerase; rfbD, dTDP-4-dehydrorhamnose reductase; sdhABCD, succinate dehydrogenase , ß, , and chains; sucDC, succinyl-CoA synthetase and ß chains; snt, sugar phosphate nucleotidyltransferase; tpi, TIM; tpsp, TPSP; treS/treP, trehalose synthase or trehalose phosphorylase.

(ii) Variants of the Entner-Doudoroff (ED) pathway. (a) As-signment of genes involved in the non- and semiphosphorylative ED pathways. For reconstruction of the nonphosphorylative ED pathway of T. tenax, we screened for genes encoding the enzymes GDH, gluconolactonase, gluconate dehydratase (GAD), KDGA, aldehyde oxidoreductase/dehydrogenase, and glycerate kinase. Homologs of all these genes—with exception of a homolog of the gene encoding gluconolactonase—were found in the T. tenax genome sequence (Table 2; Fig. 2). The lack of the gluconolactonase gene homolog in T. tenax indicates that this enzyme (EC 3.1.1.17) is either absent or not related to known bacterial counterparts (e.g., Zymomonas mobilis). So far, archaeal gluconolactonase homologs were only identified in A. fulgidus and S. solfataricus (72). The assignment of the genes encoding GDH (gdh), GAD (gad), and KDGA (kdgA) was confirmed by identifying their encoded products. Thus, the gdh gene was unequivocally confirmed by identification of the N-terminal sequence of the GDH isolated from T. tenax cells (57) and of the activity of the recombinant protein expressed in E. coli. The GAD gene (gad) was recognized due to its neighborhood to the KDGA gene, and this assignment was confirmed by the enzymatic properties of the recombinant gene product (H. Ahmed, T. Ettema, B. Tjaden, J. Van der Oost, and B. Siebers, unpublished data). The gene showed significant similarity to those encoding members of the mandelate racemase/muconate-lactonizing enzyme family (ENO superfamily), which comprises numerous different enzyme activities (e.g., those of muconate cycloisomerase, galactonate dehydratase, putative isomerase, and the starvation-sensing protein RspA). The kdgA gene has been expressed in E. coli and biochemically characterized (Ahmed et al., unpublished data).

View larger version (26K): [in this window] [in a new window]   FIG. 2. Variants of the ED pathway in T. tenax. The enzymes identified in the T. tenax genome are shown in grey-shaded boxes, and the previously characterized GDH is shown in an open box. The central intermediate KDG of both variants is boxed. Closed arrows indicate the reactions of the nonphosphorylative version of the ED pathway; open arrows indicate the reactions of the semiphosphorylative version. Grey-shaded arrows mark reactions involved in both versions. For abbreviations of enzyme names, see Table 2.

Surprisingly, an ORF with high-level similarity to genes encoding KDG kinases (KDGKs) and fructokinases was identified in the T. tenax genome. The close neighborhood to the kdgA gene suggests an involvement in the ED pathway (see below). The KDGK gene was expressed heterologously in E. coli and the corresponding enzyme activity has been confirmed, indicating that the semiphosphorylative ED pathway is also active in T. tenax (Ahmed et al., unpublished data). This finding is rather unexpected, since the semiphosphorylative ED pathway has been assumed to be characteristic of haloarchaea. The localization of a gene homolog coding for glucan-1,4--glucosidase (GAA) (gaa; see Table 5) in the close neighborhood of kdgK suggests that the ED variants are involved in the hydrolytic degradation of polysaccharides (e.g., glycogen) in T. tenax (see below). A homolog of the glucose-6-phosphate dehydrogenase gene of the kind typical of the classical ED pathway was not found in the T. tenax genome.

(b) Organization of the ED genes. Analysis of the organization of ED genes (Table 2; Fig. 2) revealed the presence of a gene cluster which comprises the genes coding for GAD (gad), KD(P)GA (kdgA), KDGK (kdgK), and GAA (gaa). The kdgA and kdgK genes as well as the kdgK and gaa genes overlap by 4 bp, giving good evidence that the three genes form an operon (Fig. 3B). The gene encoding the GAD is localized 68 bp downstream of the putative kdgA-kdgK-gaa operon and is oriented in the opposite direction (Fig. 3A).

(iii) Pentose metabolism. Assignment of genes involved in pentose metabolism. T. tenax does not possess the complete gene set for oxidative degradation of carbohydrates via the PPP (Table 3). As seen with S. solfataricus (52), only the gene homologs for ribosephosphate isomerase and transketolase (partial sequence information for the N- and C-terminal sections) were found to date. Genes coding for the classic glucose-6-phosphate dehydrogenase, decarboxylating phosphogluconate dehydrogenase, ribulosephosphate-3-epimerase, and transaldolase seem to be absent. These findings might indicate that hyperthermophilic archaea either use different, nonhomologous enzymes (which show no sequence similarity to known counterparts) in the PPP or utilize a yet unknown, analogous pathway for carbohydrate catabolism. Thus, the biosynthesis of pentoses in hyperthermophilic archaea remains unclear. The presence of a nonoxidative PPP (as proposed for ribose biosynthesis in Methanococcus maripaludis) (67) seems to be improbable, since a gene homolog for the common ribulosephosphate-3-epimerase was not found in the T. tenax genome. However, two ORFs were detected in the T. tenax genome which showed similarity to genes coding for well-known enzymes of pentose metabolism such as ribokinase (sugar kinase of the ribokinase [pfkB] family) (ATP + D-ribose ADP + D-ribose 5-phosphate) and deoxyribose-phosphate aldolase (2-deoxy-D-ribose 5-phosphate D-glyceraldehyde 3-phosphate + acetaldehyde).

The identification of the genes encoding enzymes determining the direction of the cycle—with the exception the gene encoding citrate synthase—reveals some problems. (i) For the citrate lyase, only a gene homolog for the ß chain [citryl-CoA lyase subunit; (3S)-citryl-CoA acetyl-CoA + oxaloacetate] was found but homologs for subunits [citrate-CoA transferase subunit; acetyl-CoA + citrate acetate + (3S)-citryl-CoA] and (acyl carrier protein subunit) are missing [Leuconostoc mesenteroides (ß)₆ subunit composition] (2). To our knowledge nothing is known so far about the reaction mechanism of archaeal citrate lyases, although ATP citrate lyase activity was demonstrated in cell extracts of T. neutrophilus (1). The gene coding for the ß subunit seems to form an operon with three gene homologs coding for a putative transport protein (ptp; AJ621305), putative acetyl-CoA synthetases (acetate-CoA ligase gene and acs; AJ621304), and acetyl-CoA transferases or carnitine dehydratases (act; AJ621303). Thus, acetate-CoA ligase (EC 6.2.1.1; ATP + acetate + CoA AMP + P_i + acetyl-CoA) might be responsible for the generation of acetyl-CoA, which is considered to be a free analog of the acetylated acyl carrier protein ( subunit) (2) and the putative acetyl-CoA transferase might represent the citrate-CoA transferase subunit . However, for a confidential assignment of citrate lyase experimental analyses have to be awaited. (ii) For the catabolic enzyme 2-oxoglutarate dehydrogenase, only the E3 component (dihydrolipoamide dehydrogenase; two candidates) of the multienzyme complex consisting of E1, dihydrolipoamide-S-succinyltransferase (E2), and E3 was found, which can be shared by different 2-oxo acid dehydrogenases and the glycine cleavage system (22, 74). (iii) KOR, which is thought to catalyze the reaction in the anabolic direction, could not be assigned unequivocally because of the high level of similarity between the different ferredoxin-dependent oxidoreductases. Two clusters of gene homologs were identified with similarity to genes encoding POR (pyruvate:ferredoxin oxidoreductase) and related 2-oxoacid synthases (2-oxoacid:ferredoxin oxidoreductases [OOR]). Both are organized in an operon which comprises genes encoding the , ß, , and subunits and thus—because of their four-subunit composition—might represent either POR or KOR (22, 35, 62). In addition, four additional gene clusters were identified which encode homologs for 2-oxoacid synthases (accession no. AJ621306, AJ621307, AJ621313, AJ621314, AJ621318, AJ621319, AJ621320, AJ621352, and AJ621353). (iv) For the succinate dehydrogenase (catabolic) and fumarate reductase (anabolic), two gene clusters were identified which allow no differentiation between the two enzymes due to their similarity. However, the genetic organization might give some evidence for the respective enzyme activities. We assume that the putative succinate dehydrogenase (candidate 1) is organized in an operon which encodes all four subunits (, ß, , and ), whereas the fumarate reductase (candidate 2) is organized in a gene cluster with homologs coding for the flavoprotein subunit (chain ), the iron-sulfur protein (chain ß; homology only to the N-terminal part), and a zinc-dependent class III alcohol dehydrogenase (AJ621280). The gene displacement in noncatalytic subunits of fumarate reductases has been previously described (22); thus, whereas the membrane-bound enzymes are composed of four or three subunits, the cytoplasmatic enzyme exhibits only one type of subunit.

In summary, all these findings suggest that the CAC in T. tenax is operative in the reductive and oxidative direction. The only missing enzyme is 2-oxoglutarate dehydrogenase. However, measurements in crude extracts demonstrate 2-oxoglutarate:benzylviologen oxidoreductase activity (49) and more recent reports about the POR in Clostridium thermoaceticum (14) demonstrate that the enzyme is active with ferredoxin as a cosubstrate in the oxidative as well as in the reductive reaction. Thus, 2-oxoglutarate oxidoreductase might be involved in both the oxidative and the reductive CAC in T. tenax.

No convincing candidates for the two key enzymes of the glyoxylate cycle—isocitrate lyase and malate synthase—were identified, suggesting that a functional glyoxylate cycle is absent from T. tenax. Surprisingly, putative gene homologs coding for the three subunits of the carbon monoxide dehydrogenase (medium, small, and large chains; AJ621345, AJ621344, and AJ621343, respectively) were identified in the T. tenax genome. All three subunits overlap by 4 bp and thus obviously form an operon. However, none of the other genes for a functional Wood-Ljungdahl pathway (e.g., from formylmethanofuran to methyl-tetrahydromethanopterin) were found, suggesting that a reductive acetyl-CoA pathway, which is described for several archaea (e.g., M. thermautotrophicus, Methanopyrus kandleri and in autotrophic species of Archaeoglobus [47]), is absent from T. tenax. However, the carbon monoxide dehydrogenase might be involved in coupling the POR reaction to CO oxidation via ferredoxin, as described for C. thermoaceticum (14).

(b) Organization of the CAC genes. Obviously, several CAC genes are organized in operons or operon-like structures (Fig. 3; Table 4). Seven operons have been identified so far. The gene coding for the ß subunit of citrate lyase (citE) seems to form an operon with the upstream-located gene homologs coding for a putative transport protein (ptp), acetyl-CoA synthetase (acs), and acetyl-CoA transferases (act) or carnitine dehydratases. The ptp and acs as well as the act and citE genes overlap by 4 bp and the acs and act genes overlap by 17 bp, forming the ptp-acs-act-citE operon and suggesting a functional connection of the genes. Further, the genes coding for the reversible enzymes TIM (tpi) and aconitase (acn) (see above) are linked in an operon (acn-tpi) (as confirmed by Northern analysis) (63), indicating a strong correlation between the reversible CAC and EMP pathways. The genes coding for the two putative candidates for KOR, e.g., oxoacid synthase (oxoacid:ferredoxin oxidoreductase; oor gene), are organized in operons which comprise all four subunits (, ß, , and ). Whereas for candidate 1 each subunit is encoded by one gene (operon oorC-D-A-B), the second candidate comprises subunit and in one ORF (operon oorCD-A-B). All genes overlap by 4 bp. The fifth operon comprises two genes coding for and ß subunits of succinyl-CoA synthetase (sucC-D operon) which overlap by 4 bp. The sixth operon harbors four genes coding for the putative succinate dehydrogenase (sdh gene; candidate 1 Sdh/Frd). Whereas the downstream sdhA-located genes coding for the ß, , and subunits overlap by 4 bp, the upstream-located sdhA gene coding for the subunit is separated by 2 bp. The presence of putative promoter structures in front of sdhA only (Fig. 3B) indicates that all four genes form the sdhA-B-C-D operon. A second candidate, which encodes a putative fumarate reductase and comprises two genes with homology to sdhA/frdA and sdhB/frdB, was identified in the T. tenax genome. Here both ORFs overlap by 56 bp but are oriented in opposite directions, and the sdhB/frdB gene homolog overlaps by 4 bp with the adh gene which codes for a putative Zn-dependent alcohol dehydrogenase (adh-frdB operon).

(v) Glycogen metabolism. (a) Assignment of genes involved in glycogen metabolism. By the use of the T. tenax genome, gene homologs for almost all enzymes necessary for catalyzing the synthesis and degradation of glycogen were found (Table 5). As determined on the basis of sequence similarity, seven genes are thought to encode sugar phosphate nucleotidyl transferases (snt) catalyzing the synthesis of the nucleotide-activated glucose. A possible role in glycogen synthesis is proposed for candidate 3 (snt3), which is found closely neighboring genes involved in glycogen metabolism (Fig. 3A). Because of the relatedness and numerous different annotations of gene homologs in the databases, a prediction of substrate specificity was not possible. Only for candidate 7, glucose-1-phosphate thymidylyltransferase (dTDP-glucose pyrophosphorylase [RfbA]; EC 2.7.7.24), can activity be suggested due to the closely neighboring genes coding for dTDP-glucose-4,6-dehydratase (RfbB; EC 4.2.1.46), dTDP-4-dehydrorhamnose reductase (RfbD; EC 1.1.1.133), and dTDP-4-dehydrorhamnose-3,5-epimerase (RfbC; EC 5.1.3.13). All four genes are involved in the biosynthesis of dTDP-rhamnose, which represents a precursor for the O-specific polysaccharide chain in the lipopolysaccharide of the outer membrane in most gram-negative bacteria (36), thus raising questions about their functions in archaea.

Only one gene was identified which shows convincing similarity to genes encoding glycogen synthase, which catalyzes the succeeding polymerization reaction. No gene homolog was found for the 1,4--glucan branching enzyme, which catalyzes the formation of -1-6 bonds. For glycogen degradation, four different genes were identified by sequence similarity: the genes coding for glycogen phosphorylase, -amylase, GAA, and -glucosidase (maltase).

The glgP gene coding for glycogen phosphorylase was identified in T. tenax due to its similarity to genes encoding archaeal and bacterial glucan, glycogen, and maltodextrin phosphorylases. The glgP gene (1,458 bp), which codes for a protein of 485 amino acids, was expressed in E. coli. The recombinant protein was purified by heat precipitation and ion exchange chromatography and revealed a subunit molecular mass of 55 kDa (as determined by SDS-PAGE), a result in good agreement with the calculated molecular mass of 55.56 kDa. The enzyme catalyzes the P_i-dependent degradation of glycogen, forming glucose 1-phosphate. Glycogen phosphorylase follows Michaelis-Menten kinetics, with K_m and V_max values of 0.09 mg/ml and 0.9 U/mg of protein for glycogen and 0.09 mM and 0.9 U/mg of protein for P_i at 50°C. Thus, the glycogen phosphorylase seems to play a major role in the phosphorolytic degradation of glycogen.

Genes encoding debranching enzymes (such as pullulanases) seem to be missing from the T. tenax genome. Since most of the GAAs described are able to hydrolyze 1,6--bonds also, however, this enzyme might compensate for the missing activity of debranching enzyme. The gaa gene is found closely neighboring the genes coding for the key enzymes of the modified ED pathways, indicating the close functional relationship between hydrolytic degradation of polysaccharides and the ED variants. Finally, one gene homolog was found with similarity to genes encoding phosphomannomutase or PGM. Whether the respective gene product exhibits PGM activity and thus is engaged in the catabolic and anabolic directions remains to be analyzed.

In summary, most gene homologs necessary for glycogen metabolism were identified in the T. tenax genome, and the glycogen phosphorylase, which is involved in the phosphorolytic degradation of glycogen, is characterized. The only missing counterparts are enzymes that are involved in the formation (branching enzyme) and hydrolysis of 1,6- bonds. Whether these enzymes show no similarity to known enzymes or are generally absent and their activity is taken over by glycogen synthase and GAA remains to be proven.

(b) Organization of genes involved in glycogen metabolism. Three of the seven gene homologs coding for sugar phosphate nucleotidyl transferases are organized in operons (Fig. 3; Table 5). The gene of candidate 1 (snt1) forms an operon with the gltX gene (AJ621274 and AJ621275; coding for glutamyl tRNA synthetase) and orfV (AJ621276; coding for a hypothetical protein). snt1 and gltX overlap by 14 bp, and gltX and orfV overlap by 4 bp, forming the snt1-gltX-orfV operon. The gene of candidate 2 (snt2) overlaps by 4 bp with orfZ (AJ621323), which encodes a hypothetical protein in P. aerophilum (orfZ-snt2 operon). A functional connection is proposed for the putative glucose-1-phosphate thymidylyltransferase gene (rfbA; candidate 7) which forms an operon with the rfbB (4-bp overlap; dTDP-glucose-4,6-dehydratase; AJ621291), rfbD (10-bp overlap; dTDP-4-dehydrorhamnose reductase; AJ621290) and rfbC (3-bp gap; dTDP-4-dehydrorhamnose-3,5-epimerase; AJ621289) genes (rfbA-B-D-C operon). Four of the genes involved in the glycogen metabolism, the genes encoding glycogen phosphorylase (glgP), -amylase (amyA), glycogen synthase (glgA), and sugar phosphate nucleotidyl transferases (candidate 3; snt3), are clustered in the genome of T. tenax. The glgA gene is located upstream of a gene homolog (mthfs; AJ621295) which encodes a putative 5-formyltetrahydrofolate cyclo-ligase and is oriented in the opposite direction. Both ORFs overlap by 22 bp. The amyA, glgA, and (separated by 44 bp) snt3 genes are localized downstream of the mthfs gene (75 bp) and in the same orientation as glgP. Since both amyA and glgA overlap by 1 bp (amyA-glgA operon), we assume that a functional linkage existed between these two genes at least, suggesting that the -amylase is involved in glycogen metabolism. A similar combination of genes involved in glycogen degradation and synthesis has been previously reported for bacterial species (41). As described above, the gaa gene coding for GAA is organized in an operon with genes coding for the key enzymes of the ED pathway (see "Organization of the ED genes," above).

(vi) Trehalose metabolism. (a) Assignment of genes involved in trehalose metabolism. Sequence searches clearly indicated that trehalose is synthesized in T. tenax via the OtsA-OtsB pathway (Table 6). In this pathway, trehalose 6-phosphate is formed from glucose 6-phosphate and UDP-glucose by trehalose-6-phosphate synthase (TPS, or OtsA) and trehalose 6-phosphate is dephosphorylated by trehalose-6-phosphate phosphatase (TPP, or OtsB), as shown for E. coli and S. cerevisiae (16, 24). Strikingly, only one gene that was obviously a chimera composed of two domains, one coding for TPS and the other coding for TPP, was identified in T. tenax. Thus, the protein represents a TPS-phosphatase (TPSP). First studies of the recombinant TPSP demonstrate that the encoded enzyme exhibits both activities (M. Zaparty, unpublished data). Besides this key enzyme of trehalose synthesis, further enzymes involved in the trehalose metabolism were identified: PGM and (NTP) sugar phosphate nucleotidyltransferases, which are engaged in the glycogen metabolism, are also constituents of trehalose metabolism. Additionally, a gene homolog with similarity to either the gene encoding trehalose synthase (TreS; treS gene) or the gene encoding ,-trehalose phosphorylase (TreP; treP gene) seems to be present in T. tenax. Whereas TreS catalyzes the conversion of the -1,4 linkage in maltose to -1,1 generating trehalose (specific for some bacteria, e.g., Thermus aquaticus and Pimelobacter) (65), TreP catalyzes the phosphorolytic degradation of trehalose, forming glucose and glucose 1-phosphate. So far, no biochemical information is available concerning whether the gene product represents an alternative pathway for trehalose synthesis or is involved in trehalose degradation.

(b) Organization of genes involved in trehalose metabolism. The organization of genes also involved in glycogen metabolism is discussed above (see also Table 5). The tpsp gene coding for the unusual TPSP is localized downstream of an ORF coding for a putative glycosyltransferase (gt); the latter overlaps by 4 bp with an upstream-located ORF (hp; AJ621288) which shows homology only to a gene coding for a hypothetical protein in P. aerophilum (Fig. 3; Table 6). Thus, at least the genes coding for the putative glycosyltransferase and the hypothetical protein seem to form an operon (hp-gt operon). Also, an ORF (orfY; AJ621342) (upstream of the treS/treP gene coding for the trehalose synthase or trehalose phosphorylase) was found with similarity to a gene encoding a hypothetical protein in P. aerophilum, S. solfataricus, and S. tokodaii. Both genes overlap by 4 bp and form the orfY-treS/treP operon (Fig. 3), but so far we have no indication about the function of the two genes.

Conclusions. T. tenax genome data reveal the presence of pathways involved in central carbohydrate metabolism as well as in glycogen and trehalose turnover (Fig. 4). So far, T. tenax and S. solfataricus are the only facultative heterotrophic archaea sequenced. Thus, the presented analysis not only contributes to the diversity of carbon metabolism in archaea but also represents an important step towards understanding of archaeal metabolic networks, which allow adaptation to alternative growth conditions.

View larger version (21K): [in this window] [in a new window]   FIG. 4. The central carbohydrate metabolism in T. tenax. A selection of growth substrates is given at the top; the labels of the different pathways are marked by grey-shaded boxes.

The genome analysis strongly suggests the presence of three pathways for carbohydrate metabolism: (i) the EMP variant and (ii) the semiphosphorylative and (iii) the nonphosphorylative versions of the ED pathway. At the moment, we can only speculate about the physiological significance of the three different pathways. The organization of ED genes suggests that the ED variants are involved in hydrolytic degradation of glycogen, whereas the EMP variant is responsible for the phosphorolytic breakdown. The energy demand of the cell seems to have a strong influence on the selection of the different pathways. Whereas no ATP is generated by the ED variants, one ATP is gained using the EMP variant (taking into account that PP_i, the phosphoryl donor of the phosphofructokinase, is a waste product of the cell). In T. tenax, GAP, the common intermediate of EMP pathway and semiphosphorylative ED version, can be converted into 3-phosphoglycerate in a one-step reaction—via GAPN or GAPOR—thus circumventing the formation of the extremely heat-labile intermediate 1,3-bisphosphoglycerate. Since the nonphosphorylative ED pathway avoids the most heat-labile intermediates (1,3-bisphosphoglycerate, GAP, and dihydroxyacetone phosphate), this pathway could be preferred at the upper temperature range of growth.

Further, the sequence data support the idea of the existence of a reversible CAC in T. tenax, indicating that—as already suggested—the organism uses the reaction cycle not only for catabolic purposes but also for CO₂ fixation (1). Almost all gene homologs necessary for the glycogen metabolism (with the exception of the branching and debranching enzyme) were identified in the T. tenax genome. In addition, trehalose synthesis was shown to proceed via the OtsA-OtsB pathway in T. tenax. No evidence was found for a functional PPP, which has been thought to be essential for generation of pentoses and NADPH for anabolic purposes. Together with previous genome analyses of methanoarchaea and hyperthermophilic archaea, this finding suggests that the oxidative PPP is generally absent from archaea or that an analogous pathway exists with enzymes showing no significant similarity to known bacterial and eucaryal counterparts.

.

	ACKNOWLEDGMENTS

 
We thank Bettina Haberl (EPIDAUROS Biotechnologie AG) for assistance in the setup of sequencing reactions, Klaus Gellner (EPIDAUROS Biotechnologie AG) and Günter Raddatz (MPI für Entwicklungsbiologie, Tübingen) for bioinformatic support, and William Martin (Heinrich-Heine Universität Düsseldorf) for critical review of the manuscript.

This work was supported by the Deutsche Forschungsgemeinschaft. The MAGPIE project is supported by Genome Prairie, Canada.

	FOOTNOTES

 
* Corresponding author. Mailing address: FB 9, Mikrobiologie, Universität Duisburg-Essen, Universitätsstr. 5, 45117 Essen, Germany. Phone: 0049-201-1833442. Fax: 0049-201-1833990. E-mail: bettina.siebers{at}uni-essen.de.

	REFERENCES

Top Abstract Introduction Materials and Methods Results and Discussion References

 

1. Beh, M., G. Strauss, R. Huber, K.-O. Stetter, and G. Fuchs. 1993. Enzymes of the reductive citric acid cycle in the autotrophic eubacterium Aquifex pyrophilus and in the archaebacterium Thermoproteus neutrophilus. Arch. Microbiol. 160:306-311.[CrossRef]
2. Bekal, S., J. van Beeumen, B. Samyn, D. Garmyn, S. Henini, C. Diviès, and H. Prévost. 1998. Purification of Leuconostoc mesenteroides citrate lyase and cloning and characterization of the citCDEFG gene cluster. J. Bacteriol. 180:647-654.[Abstract/Free Full Text]
3. Bell, S. D., and S. P. Jackson. 2001. Mechanism and regulation of transcription in archaea. Curr. Opin. Microbiol. 4:208-213.[CrossRef][Medline]
4. Brunner, N. A., H. Brinkmann, B. Siebers, and R. Hensel. 1998. NAD⁺-dependent glyceraldehyde-3-phosphate dehydrogenase from Thermoproteus tenax. The first identified archaeal member of the aldehyde dehydrogenase superfamily is a glycolytic enzyme with unusual regulatory properties. J. Biol. Chem. 273:6149-6156.[Abstract/Free Full Text]
5. Brunner, N. A., B. Siebers, and R. Hensel. 2001. Role of two different glyceraldehyde-3-phosphate dehydrogenases in controlling the reversible Embden-Meyerhof-Parnas pathway in Thermoproteus tenax: regulation on protein and transcript level. Extremophiles 5:101-109.[CrossRef][Medline]
6. Buchanan, C. L., H. Connaris, M. J. Danson, C. D. Reeve, and D. W. Hough. 1999. An extremely thermostable aldolase from Sulfolobus solfataricus with specificity for non-phosphorylated substrates. Biochem. J. 343:563-570.
7. Budgen, N., and M. J. Danson. 1986. Metabolism of glucose via a modified Entner-Doudoroff pathway in the thermoacidophilic archaebacterium Thermoplasma acidophilum. FEBS Lett. 196:207-210.[CrossRef]
8. Colombo, S., M. Grisa, P. Tortora, and M. Vanoni. 1994. Molecular cloning, nucleotide sequence and expression of a Sulfolobus solfataricus gene encoding a class II fumarase. FEBS Lett. 340:93-98.[CrossRef][Medline]
9. De Rosa, M., A. Gambacorta, B. Nicolaus, P. Giardina, E. Poerio, and V. Buonocore. 1984. Glucose metabolism in the extreme thermoacidophilic archaebacterium Sulfolobus solfataricus. Biochem. J. 224:407-414.[Medline]
10. Di Lernia, I., A. Morana, A. Ottombrino, S. Fusco, M. Rossi, and M. De Rosa. 1998. Enzymes from Sulfolobus shibatae for the production of trehalose and glucose from starch. Extremophiles 2:409-416.[CrossRef][Medline]
11. Dörr, Ch., M. Zaparty, B. Tjaden, H. Brinkmann, and B. Siebers. 2003. The hexokinase of the hyperthermophile Thermoproteus tenax: ATP-dependent hexokinases and ADP-dependent glucokinases, two alternatives for glucose phosphorylation in archaea. J. Biol. Chem. 278:18744-18753.[Abstract/Free Full Text]
12. Fischer, F., W. Zillig, K. O. Stetter, and G. Schreiber. 1983. Chemolithoautotrophic metabolism of anaerobic extremely thermophilic archaebacteria. Nature 301:511-513.[CrossRef][Medline]
13. Fitz-Gibbon, S. T., H. Ladner, U. J. Kim, K. O. Stetter, M. I. Simon, and J. H. Miller. 2002. Genome sequence of the hyperthermophilic crenarchaeon Pyrobaculum aerophilum. Proc. Natl. Acad. Sci. USA 99:984-989.[Abstract/Free Full Text]
14. Furdui, C., and S. W. Ragsdale. 2000. The role of pyruvate ferredoxin oxidoreductase in pyruvate synthesis during autotrophic growth by the Wood-Ljungdahl pathway. J. Biol. Chem. 275:28494-28499.[Abstract/Free Full Text]
15. Gasterland, T., and C. W. Sensen. 1996. Fully automated genome analysis that reflects user needs and preferences. A detailed introduction to the MAGPIE system architecture. Biochimie 78:302-310.[Medline]
16. Giver, H. M., O. B. Styrvold, I. Kaasen, and A. R. Strom. 1988. Biochemical and genetic characterization of osmoregulatory trehalose synthesis in Escherichia coli. J. Bacteriol. 170:2841-2849.[Abstract/Free Full Text]
17. Gordon, P., T. Gaasterland, and C. W. Sensen. 2002. Genomic data representation through images: MAGPIE as an example, pp. 345-363. In C. W. Sensen (ed.), Essentials of genomics and bioinformatics. Wiley-VCH, Weinheim, Germany.
18. Graham, D. E., H. Xu, and R. H. White. 2002. A divergent archaeal member of the alkaline phosphatase binuclear metalloenzyme superfamily has phosphoglycerate mutase activity. FEBS Lett. 517:190-194.[CrossRef][Medline]
19. Hansen, T., M. Oehlmann, and P. Schönheit. 2001. Novel type of glucose-6-phosphate isomerase in the hyperthermophilic archaeon Pyrococcus furiosus. J. Bacteriol. 183:3428-3435.[Abstract/Free Full Text]
20. Hansen, T., and P. Schönheit. 2000. Purification and properties of the first-identified, archaeal, ATP-dependent 6-phosphofructokinase, an extremely thermophilic non-allosteric enzyme, from the hyperthermophile Desulfurococcus amylolyticus. Arch. Microbiol. 173:103-109.[CrossRef][Medline]
21. Hensel, R., S. Laumann, J. Lang, H. Heumann, and F. Lottspeich. 1987. Characterization of two D-glyceraldehyde-3-phosphate dehydrogenases from the extremely thermophilic archaebacterium Thermoproteus tenax. Eur. J. Biochem. 170:325-333.[Medline]
22. Huynen, M. A., T. Dandekar, and P. Bork. 1999. Variation and evolution of the citric-acid cycle: a genomic perspective. Trends Microbiol. 7:281-291.[CrossRef][Medline]