|Home | About | Journals | Submit | Contact Us | Français|
Acidilobus saccharovorans is an anaerobic, organotrophic, thermoacidophilic crenarchaeon isolated from a terrestrial hot spring. We report the complete genome sequence of A. saccharovorans, which has permitted the prediction of genes for Embden-Meyerhof and Entner-Doudoroff pathways and genes associated with the oxidative tricarboxylic acid cycle. The electron transfer chain is branched with two sites of proton translocation and is linked to the reduction of elemental sulfur and thiosulfate. The genomic data suggest an important role of the order Acidilobales in thermoacidophilic ecosystems whereby its members can perform a complete oxidation of organic substrates, closing the anaerobic carbon cycle.
Acidophilic microorganisms are widely dispersed in natural acidic environments, including volcanic hot springs, and are, in the majority, aerobes (14). However, such anoxic, high-temperature, acidic environments are inhabited by metabolically versatile anaerobic thermoacidophiles of the archaeal phylum Crenarchaeota. Lithoautotrophic thermoacidophiles oxidize molecular hydrogen in the course of elemental sulfur (S0) respiration. Organotrophs couple the oxidation of organic substrates to the reduction of S0 or thiosulfate. They all belong to the genus Acidilobus in the family Acidilobaceae and to the genus Caldisphaera in the family Caldisphaeraceae (4, 13, 22, 24). Acidilobaceae and Caldisphaeraceae form the crenarchaeal order Acidilobales (24). Acidilobus saccharovorans was isolated from an acidic hot spring of Uzon Caldera, Kamchatka, Russia (24). It is an obligately anaerobic acidophile with a range of growth from pH 2.5 to 5.8 (optimum at pH 3.5 to 4) and a temperature range from 60 to 90°C (optimum at 80 to 85°C). It utilizes a wide range of proteinaceous and carbohydrate substrates and cannot grow lithoautotrophically on H2 and CO2 (24). S0 and thiosulfate stimulate growth and are reduced to H2S. Protons cannot serve as electron acceptors, since no H2 is produced during growth in the absence of S0 (24). Genomic sequences of aerobic, thermoacidophilic euryarchaea Thermoplasma acidophilum (26) and Picrophilus torridus (8) give an insight into the thermoacidophilic survival strategy. However, no genomes of obligately anaerobic, thermoacidophilic archaea were available until now. Here we present the genome of A. saccharovorans and show that it encodes numerous hydrolytic enzymes and metabolic pathways necessary for the utilization and complete mineralization of organic substrates in its natural habitat, acidic hot springs.
A. saccharovorans 345-15T (DSM 16705) was obtained from the culture collection of the Laboratory of Hyperthermophilic Microbial Communities, the Winogradsky Institute of Microbiology of the Russian Academy of Sciences, and grown as described previously (24). The A. saccharovorans genome was sequenced on a Roche GS FLX genome sequencer by the standard protocol for a shotgun genome library. The GS FLX run (1/4 plate) resulted in the generation of about 22 Mb of sequences with an average read length of 230 bp. The GS FLX reads were assembled into 11 contigs by a GS De Novo Assembler (Roche). The contigs were oriented into scaffolds, and the complete genome sequence was obtained upon the generation and sequencing of appropriate PCR fragments. The rRNA genes were identified with RNAmmer (17). tRNA genes were located with tRNAscan-SE (19). Protein-encoding genes were identified with the GLIMMER gene finder (5), followed by a round of manual curation. Genome annotation was performed with the AutoFACT tool (16). Protein similarity was analyzed by BLASTP search (1). Signal peptides were predicted with SignalP version 3.0 (3).
A. saccharovorans has a single circular chromosome of 1,496,453 bp with no extrachromosomal elements (see Fig. S1 in the supplemental material). There is a single copy of the 16S-23S rRNA operon and a single distantly located 5S rRNA gene. A total of 45 tRNA genes coding for all 20 amino acids are scattered over the genome. The cumulative GC skew profile (9, 18) displays two local minima that likely correspond to the DNA replication origins (see Fig. S1 in the supplemental material). Genome annotation reveals 1,499 potential protein-encoding genes, 972 (65%) of which can be functionally assigned. The functions of the remaining 527 genes cannot be predicted from the deduced amino acid sequences; 246 of them are unique to A. saccharovorans.
The numbers of predicted protein-encoding genes of A. saccharovorans that have homologs in other archaeal genomes are shown in Fig. Fig.11 A. The closest relative of A. saccharovorans is the aerobic, neutrophilic crenarchaeon Aeropyrum pernix; they share about half of their proteomes (~740 proteins). Fewer A. saccharovorans proteins have homologs in the genomes of other members of the crenarchaeal orders Desulfurococcales and Thermoproteales (Fig. (Fig.1B).1B). Notably, A. saccharovorans shares higher numbers of proteins (661 to 687, 16% of best BLASTP hits) with the phylogenetically distant aerobic thermoacidophiles of the order Sulfolobales (24). Predicted functions of most of the 199 proteins of A. saccharovorans showing the best BLASTP hits in Sulfolobales proteomes are related to transport (54 proteins) or are unknown (86 proteins). A. saccharovorans and Sulfolobales often coexist in acidic hot springs, and lateral gene transfer may be responsible for the relatively high similarity between their proteomes. Extensive lateral gene transfer was suggested to explain the high fraction of genes shared by Sulfolobales and aerobic euryarchaeal thermoacidophiles P. torridus and T. acidophilum (8, 26). However, it is unlikely that extensive gene transfer occurred between A. saccharovorans and euryarchaeal thermoacidophiles, because the numbers of proteins shared with P. torridus (428) and T. acidophilum (459), as well as with nonacidophilic euryarchaea Pyrococcus furiosus (440) and Archaeoglobus fulgidus (427), are similar.
Growth of A. saccharovorans on proteins and peptides may be enabled by the function of several extracellular proteolytic enzymes. The imported peptides can be fermented with concomitant substrate-level phosphorylation to generate ATP (Fig. (Fig.2;2; see Table S1 in the supplemental material). The presence of genes for extra- and intracellular glycoside hydrolases and sugar transporters (Fig. (Fig.22 and Table S1) correlates with the growth of A. saccharovorans on carbohydrates (24). Notably, the A. saccharovorans genome encodes at least two primary ABC-type sugar transport systems (Table S1) and more than 20 secondary transporters. A similar high ratio of secondary transport systems to ATP-consuming primary transport systems was found in the thermoacidophile P. torridus (8), suggesting that the high external proton concentration is extensively used by A. saccharovorans for transport processes.
Glucose in archaea is metabolized to pyruvate via either the Embden-Meyerhof (EM) pathway or the Entner-Doudoroff (ED) pathway. Both modified EM and ED pathways are encoded by A. saccharovorans genome (Fig. (Fig.2),2), suggesting their operation in parallel as found in Thermoproteus tenax (32). The ED pathway can be beneficial at the upper temperature range of growth, since it avoids the most heat-labile intermediates (31). The EM pathway may account for the synthesis of fructose-6-phosphate used in the ribulose monophosphate pathway (Fig. (Fig.22).
The genome contains 14 genes encoding esterases and genes encoding a complete β-oxidation pathway, suggesting the ability of A. saccharovorans to utilize triacylglycerides and long-chain fatty acids available from the environment (Fig. (Fig.22 and Table S1). To our knowledge, besides A. saccharovorans, A. fulgidus is the only archaeon to encode the enzymes of β-oxidation (15).
A. saccharovorans may conserve energy by substrate-level phosphorylation as well as by oxidative phosphorylation. During growth in the presence of S0 as a terminal electron acceptor, H2S is produced and acetate is formed (24). It is likely that in the presence of S0, acetyl coenzyme A (acetyl-CoA) generated from the oxidation of substrates is concomitantly converted to acetate by fermentation pathway as well as oxidized to CO2 by anaerobic respiration (Fig. (Fig.2).2). The A. saccharovorans genome encodes all eight enzymes of the oxidative tricarboxylic acid (TCA) cycle (Fig. (Fig.22 and Table S1). This cycle operates in T. tenax and Pyrobaculum islandicum and was proposed to function in both oxidative and reductive directions, enabling the complete oxidation of organic substrates to CO2 and H2S, as well as CO2 fixation during autotrophic growth (12, 30, 31, 32). In A. saccharovorans, which cannot grow autotrophically (24), the TCA cycle appears to operate only in the oxidative direction. Indeed, citrate lyase, the key enzyme of the reductive TCA cycle, is not encoded.
During growth of A. saccharovorans by anaerobic respiration in the presence of S0, reduced ferredoxin, NAD(P)H, and H produced in the TCA cycle can be oxidized by a concerted action of a set of membrane-bound protein complexes and cytoplasmic proteins, resulting in the generation of transmembrane proton gradient. The genome analysis reveals a putative membrane-bound protein complex (ASAC_0373 to ASAC_0383), which is partially homologous to bacterial NADH:quinone oxidoreductase (respiratory complex I). From the 14 subunits (NuoA to NuoN) present in Escherichia coli complex I (7), this putative A. saccharovorans complex contains 11 subunits (NuoA to NuoD and NuoH to NuoN), while three subunits (NuoEFG) known to form the dehydrogenase domain involved in NADH binding and oxidation (28) are missing. Thus, the A. saccharovorans complex most likely accepts electrons from a donor other than NADH. ASAC_0380 is also homologous to the MbxL subunit of the membrane-bound oxidoreductase (MBX) complex from P. furiosus (29). It has been shown that in the MBX complex, the oxidation of ferredoxin and the reduction of NAD(P)+ are coupled to the generation of a proton motive force in the presence of S0 (29). A similar function of the MBX complex in the crenarchaea D. kamchatkensis (25) and S. marinus (2), as well as in the euryarchaeon T. sibiricus (20), has been proposed. Based on these data, we propose that ASAC_0373-ASAC_0383 encodes a proton-translocating ferredoxin:NAD(P)+ oxidoreductase complex (FNOR) which accepts electrons from reduced ferredoxin and transfers them to NAD(P)+, thereby translocating protons and establishing a proton gradient (Fig. (Fig.2).2). NAD(P)H is then oxidized by the encoded NAD(P)H:elemental sulfur oxidoreductase (NSR) (ASAC_1028), forming H2S, or enters other processes linked to NAD(P)H reoxidation (e.g., thiosulfate reduction). In addition, the NuoB (Nqo6) subunit, which transfers electrons to menaquinone in Thermus thermophilus complex I (27), is present in the A. saccharovorans FNOR complex (ASAC_0382), suggesting the possibility of electron transfer not only to NAD(P)+ but also to quinones in the respiratory chain and further to membrane-bound sulfur reductase (SR; ASAC_1394 to ASAC_1397) (Fig. (Fig.2;2; see Table S1 in the supplemental material).
The transmembrane proton gradient could also be generated by the function of the predicted integral membrane protein H+-translocating pyrophosphatase (PPase; ASAC_1013) closely related to functionally characterized enzyme from Pyrobaculum aerophilum (6). The pyrophosphatase may work in concert with H+-ATP synthase to scavenge energy from biosynthetic waste pyrophosphate in order to maintain the proton gradient, especially under energy-limiting conditions.
Besides involving the proposed FNOR complex, the respiratory chain of A. saccharovorans presumably involves the succinate:quinone oxidoreductase complex (SQOR; ASAC_1440 to ASAC_1443), which is the membrane-bound component of the oxidative TCA cycle. Electrons from the oxidation of succinate in the TCA cycle would be delivered by SQOR to the quinones in the membrane.
The genome contains a gene cluster (corresponding to ASAC_0264 to ASAC_0267) encoding two subunits of electron transfer flavoprotein (EtfAB; ASAC_0264 and ASAC_0265, respectively), a ferredoxin-like protein (ASAC_0266), and an electron transfer flavoprotein-quinone oxidoreductase (Etf-QOR; ASAC_0267). In analogy to the mitochondrial enzymes, ASAC_0264 to ASAC_0267 might be involved in the oxidation of fatty acids and some amino acids (10). These substrates would be oxidized by acyl-CoA dehydrogenases (FadE) with the Etf as an acceptor, which would pass the electrons to Etf-QOR and further to quinone in the respiratory chain of A. saccharovorans.
The genome analysis also suggests a second site of S0 reduction in A. saccharovorans. It involves a predicted membrane-bound SR complex, which would function as the terminal oxidase of S0-dependent respiration. Based on the sequence analyses, we propose that the SR is composed of the hydrophilic catalytic subunit ASAC_1397, the electron transfer subunit ASAC_1396, and the membrane-bound subunits ASAC_1395 and ASAC_1394.
Summarizing the data on the electron transfer in A. saccharovorans, we propose that its respiratory chain functioning during the growth with S0 is branched (Fig. (Fig.2).2). First, the electrons from reduced ferredoxin could be utilized by FNOR and NSR. It would be linked to proton translocation and the first site of S0 reduction. Second, a part of the electrons from FNOR, as well as the electrons from SQOR and Etf-QOR, would be transferred to the quinone pool. Oxidation of quinol by the SR would also be linked to proton translocation and the second site of reduction of S0. During growth in the absence of S0, the TCA cycle is not functional and the respiratory chain involving the SQOR complex, quinones, and SR seems not to operate. However, a proton gradient still could be established by the oxidation of reduced ferredoxin in the reaction of the proposed FNOR complex linked to reoxidation of NAD(P)H in reactions yielding reduced fermentation products.
Some anaerobic Crenarchaeota were reported to use in addition to S0 other terminal electron acceptors, including sulfate, thiosulfate, and nitrate. The presence of thiosulfate stimulated the growth of A. saccharovorans and led to production of H2S (24). We have predicted the function of ASAC_1397 as a component of membrane-bound SR complex (see above). Nevertheless, in analogy with thiosulfate reductase of Salmonella enterica involved also in the reduction of S0 (11), the dual function of ASAC_1397 in the reduction of S0 as well as of thiosulfate can be suggested.
The established proton gradient can be used by ATP synthase for the synthesis of ATP. The genome of A. saccharovorans encodes a single A1Ao ATP synthase composed of nine subunits (Table S1). The subunit K (c) (ASAC_0393) of the membrane-embedded motor does not contain conserved residues involved in Na+ binding in Na+-ATP synthases (21), strongly suggesting that H+ is the coupling ion in the A1Ao ATP synthase. During growth of A. saccharovorans in the absence of S0, the required proton gradient could be established only partially by the function of the FNOR complex and H+-translocating PPase. An additional proton gradient required for maintenance of normal intracellular pH could be generated by the reverse activity of A1Ao ATP synthase.
Organic substances reaching terrestrial acidic hot springs from surrounding areas as well as produced by lithoautotrophic microorganisms serve as substrates for organotrophic thermoacidophiles. Metabolic activity of these organisms is supported by fast anaerobic oxidation of [14C]sucrose to CO2 in acidic hot springs (23). Catabolic pathways and energy generation mechanisms predicted from the genome of A. saccharovorans explain its growth on proteins, carbohydrates, and lipids (Fig. (Fig.2).2). Though phylogenetically Acidilobales are most closely related to Desulfurococcales, they metabolically resemble Thermoproteales. Similarly to T. tenax, A. saccharovorans possesses both the modified EM and ED pathways of glucose catabolism. Another feature common between A. saccharovorans and Thermoproteales is the presence of the TCA cycle, which enables a complete oxidation of organic substances in the presence of external anaerobic electron acceptors. Adaptation of A. saccharovorans to acidic conditions is reflected by a high ratio of secondary transporters to ATP-dependent primary transporters, the suggested reversibility of H+-ATP synthase, the function of the encoded H+-translocating PPase, and the proposed two proton translocation sites in the anaerobic respiratory chain.
Physiological functions predicted by A. saccharovorans genome analysis (Fig. (Fig.2)2) and the wide distribution of Acidilobus spp. in acidic hot springs (23, 24) suggest an important role of Acidilobales in microbial communities of such habitats, closing the anaerobic carbon cycle through complete mineralization of organic substrates via S0 respiration. Contrary to mesophilic microbial communities, in which the process of anaerobic destruction of organic substances consists of several stages performed by different groups of microorganisms (hydrolysis, fermentation, syntrophic reactions, and terminal oxidation), in thermal environments this could be achieved by a single group of organisms, Thermoproteales in neutral hot springs and Acidilobales in acidic ones.
The annotated genome sequence of A. saccharovorans 345-15T has been deposited in the GenBank database under accession number CP001742.
This work was supported by the Federal Agency for Science and Innovations of Russia (contract 02.512.11.2201).
The work of M.I.P. and E.A.B.-O. on the analysis of growth characteristics of A. saccharovorans was supported by the Molecular and Cellular Biology program of RAS.
Published ahead of print on 25 June 2010.
‡Supplemental material for this article may be found at http://aem.asm.org/.