Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2008, p. 6720-6729, Vol. 74, No. 21
0099-2240/08/$08.00+0 doi:10.1128/AEM.00968-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Marcel R. A. Verhaart,1,
Amy L. VanFossen,2,
Karin Willquist,3
Derrick L. Lewis,2
Jason D. Nichols,2
Heleen P. Goorissen,1
Emmanuel F. Mongodin,4
Karen E. Nelson,5
Ed W. J. van Niel,3
Alfons J. M. Stams,1
Donald E. Ward,6
Willem M. de Vos,1
John van der Oost,1
Robert M. Kelly,2 and
Servé W. M. Kengen1*
Laboratory of Microbiology, Wageningen University, Wageningen, The Netherlands,1 Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695-7905,2 Applied Microbiology, Lund Institute of Technology, Lund University, SE-221 00 Lund, Sweden,3 Institute for Genome Sciences, University of Maryland School of Medicine, Baltimore, Maryland 21201,4 The J. Craig Venter Institute, 9712 Medical Center Drive, Rockville, Maryland 20850,5 Genencor International, a Danisco Company, 925 Page Mill Road, Palo Alto, California 943046
Received 29 April 2008/ Accepted 27 August 2008
|
|
|---|
|
|
|---|
|
|
|---|
Genome sequencing and assembly.
High-molecular-weight genomic DNA was provided to the U.S. Department of Energy Joint Genome Institute (http://www.jgi.doe.gov/) for cloning and shotgun sequencing. A combination of small (average insert sizes, 3 and 8 kb) and large (40 kb, fosmid) insert libraries were prepared and used for analysis, as indicated at http://www.jgi.doe.gov/.
Genome annotation and comparative analysis.
Critica (2) and Glimmer (13) software programs were used for coding region detection and gene identification. TMHMM 2.0 (23) was used to predict transmembrane helices in translated sequences. SignalP v2.0b2 (33) was used to predict the presence and location of N-terminal signal peptides. All automatic gene and function predictions were manually checked using BLAST programs (1), InterPro (31), and the Integrated Microbial Genomes system (28) and corrected if necessary. Protein functions were checked by using the Carbohydrate-Active enzymes (CAZy) (12; http://www.cazy.org) classification. A comparative analysis was performed based on assignment and classification of Clusters of Orthologous Groups (COG) of proteins (44) by using the Integrated Microbial Genomes system.
Growth experiments and RNA isolation.
C. saccharolyticus was subcultured overnight three times using the substrate of interest in modified DSMZ 640 medium before it was inoculated into a pH-controlled (pH 7) 1-liter fermentor containing 4 g substrate per liter. Cells were grown at 70°C until mid-logarithmic phase (optical density at 660 nm, 0.3 to 0.4) and were harvested by centrifugation and rapid cooling to 4°C and then stored at –80°C. Total RNA was isolated using a modified Trizol (Invitrogen) protocol in combination with an RNAeasy kit (Qiagen). Quality was determined with an Experion bioanalyzer (Bio-Rad), and cDNA was constructed with Superscript III reverse transcriptase (Invitrogen).
Whole-genome oligonucleotide DNA microarray design and construction.
A DNA microarray was designed and constructed based on 2,695 protein-encoding sequences in the C. saccharolyticus genome obtained from the Department of Energy's Joint Genome Institute (http://genome.ornl.gov/microbial/csac). OligoArray 2.0 (40) was used to generate one 60-mer oligonucleotide probe sequence for each open reading frame. The probes were synthesized (Integrated DNA Technologies, Iowa), resuspended in 50% dimethyl sulfoxide, and printed on Ultragap microarray slides (Corning, New York) using a QArrayMini arrayer (Genetix, United Kingdom). Each probe was spotted five times on each array to fortify the statistical analysis.
Microarray hybridization.
The cDNA samples were processed using a QIAquick purification kit (Qiagen, California), and they were eluted using phosphate buffer. The quantity and quality of the recovered cDNA samples were subsequently analyzed by using the ratio of absorbance at 260 to absorbance at 280 nm. Cyanine-3 and Cyanine-5 (Amersham, United Kingdom) dye labeling and sample hybridization were performed by following the instructions provided by TIGR (http://www.tigr.org/tdb/microarray/protocolsTIGR.shtml), with minor adjustments to accommodate long-oligonucleotide platforms. Samples were hybridized using a four-slide loop (see Fig. S1 in the supplemental material).
Data collection and analysis.
After incubation, slides were washed to remove nonspecifically bound material and scanned with a ScanArray Lite microarray scanner (Perkin Elmer, Massachusetts). Data acquisition and spot quantitation were performed with the ScanArray Express software. Once all the slides were quantitated, data from the loop were analyzed with JMP Genomics 3.0 (SAS, North Carolina), as described previously (36), using a mixed-effects analysis of variance model (49).
Nucleotide sequence and gene expression data accession numbers.
The complete final genome assembly was released on 8 May 2007 and was deposited in the GenBank database under accession number CP000679 (http://genome.jgi-psf.org/finished_microbes/calsa/calsa.home.html). The gene expression data have been deposited in the Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/projects/geo/) under accession number GSE11153.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. General features of the C. saccharolyticus genome
|
![]() View larger version (60K): [in a new window] |
FIG. 1. Diagram of the C. saccharolyticus chromosome. From outside to inside the circles show (i) the genomic positions, (ii) the coding sequences on the positive negative strand, (iii) the coding sequences on the negative strand (the colors used for the coding sequences indicate the COG functional categories), (iv) the tRNA genes, (v) the G+C content, and (vi) the GC skew. The Microbial Genome Viewer was used to construct the circular chromosome wheel (21).
|
|
View this table: [in a new window] |
TABLE 2. Functional categories of predicted open reading frames in the genomes of hydrogen-producing organisms
|
Central carbon metabolism.
C. saccharolyticus is able to metabolize a wide variety of carbohydrates, including the monosaccharides D-glucose, D-xylose, D-fructose, D-galactose, DL-arabinose, D-mannose, L-rhamnose, and L-fucose, as well as
- and β-linked di- and polysaccharides, including maltose, starch, pullulan, sucrose, trehalose, amorphous and microcrystalline cellulose, xylan, locust bean gum, and pectin (37). Once they are hydrolyzed, sugars are channeled to the central catabolic pathways (Fig. 2). Analysis of the genome sequence revealed components of a complete Embden-Meyerhof (EM) pathway, including a ROK family glucokinase (Csac_0778), 6-phosphofructokinase, (Csac_2366 and Csac_1830), a bifunctional phosphoglucose/phosphomannose isomerase (Csac_1187), fructose-1,6-bisphosphate aldolase (Csac_1189), and pyruvate kinase (Csac_1831), as well as pyruvate-phosphate dikinase (PPDK) (Csac_1955). Also, a gapA operon was found, consisting of genes encoding glyceraldehyde-3-phosphate dehydrogenase (Csac_1953), the fusion protein phosphoglycerate kinase/triose-phosphate isomerase (Csac_1952), phosphoglycerate mutase (Csac_1951), and enolase (Csac_1950) (Fig. 2). However, the oxidative branch of the pentose phosphate pathway (PPP) and the Entner-Doudoroff pathway were not detected, which is consistent with previous findings obtained using 13C nuclear magnetic resonance (14). The absence of the oxidative branch of the PPP, however, raises questions about how NADPH is produced for biosynthesis. The only other obvious NADPH-producing reaction is the isocitrate dehydrogenase (Csac_0751) reaction. However, based on sequence homology, the isocitrate dehydrogenase is likely to produce NADH instead of NADPH. Also, no obvious homolog of an NADPH-producing glyceraldehyde-3-phosphate dehydrogenase could be identified, as has been reported previously for Streptococcus species and some clostridia (8). NADPH can also be synthesized from NADH by a transhydrogenase, which is either membrane bound or soluble. In the genome there are no genes encoding orthologs of subunits of the membrane-bound type, but several homologous genes coding for soluble flavin adenine dinucleotide-dependent pyridine nucleotide-disulfide oxidoreductases are present (Csac_0759, Csac_2199, and Csac_0402). The exact physiological roles of these enzymes have not been determined. Furthermore, no ferredoxin:NADPH reductase homolog is present, although such activity has been found in some Thermoanaerobacter spp. (19).
![]() View larger version (63K): [in a new window] |
FIG. 2. Overview of the carbon metabolism and transport systems in C. saccharolyticus. The identities of the various ABC-type sugar transporters are not known. Secondary transport systems may be involved as well. Abbreviations: Rhu-1-P, rhamnulose-1-phosphate; F-1-P, fructose-1-phosphate; G-1-P, glucose-1-phosphate; G-6-P, glucose-6-phosphate; F-6-P, fructose-6-phosphate; F-1,6-bisP, fructose-1,6-bisphosphate; L-Ru-5-P, L-ribulose-5-phosphate; D-Ru-5-P, D-ribulose-5-phosphate; D-Xu-5-P, D-xylulose-5-phosphate; D-Ri-5-P, D-ribose-5-phosphate; DKI, 5-keto-4-deoxyuronate; DKII, 2,5-diketo-3-deoxygluconate; DHAP, dihydroxyacetone phosphate; GAP, glyceraldehyde-3-phosphate; KDPG, KDG phosphate; Fdred, reduced ferredoxin; 3-PG, 3-phosphoglycerate; 2-PG, 2-phosphoglycerate; FADH, reduced flavin adenine dinucleotide; Ac-CoA, acetyl coenzyme A; Ac-P, acetyl phosphate; br., branched.
|
Pyruvate, the end product of the EM pathway, is decarboxylated to acetyl coenzyme A (acetyl-CoA) by pyruvate:ferredoxin oxidoreductase (POR). C. saccharolyticus contains three 2-oxoacid:ferredoxin oxidoreductase enzyme complexes (Csac_2248 and Csac_2249, Csac_1458 to Csac_1461, and Csac_1548 to Csac_1551). According to transcriptional response data (see below), the true POR is probably encoded by the Csac_1458 to Csac_1461 genes. Acetyl-CoA is used to generate acetate and ATP (Csac_2040 and Csac_2041), or it enters the tricarboxylic acid (TCA) cycle for biosynthetic purposes. The TCA cycle in C. saccharolyticus is incomplete; it has an oxidative branch to succinyl-CoA catalyzed by a citrate (Re) synthase (Csac_0746), aconitate hydratase (Csac_0750), isocitrate dehydrogenase (Csac_0751), and the 2-oxoglutarate:ferredoxin oxidoreductase complex (Csac_1548 to Csac_1551). In the reductive direction, only orthologs of the subunits of fumarate hydratase were detected with a high level of confidence (Csac_2759/Csac_0738). Malate dehydrogenase (oxaloacetate decarboxylating) (Csac_2059) may be used to generate malate directly from pyruvate instead of from oxaloacetate. Fumarate reductase, however, could not be identified, which is in agreement with the absence of this enzyme in related clostridia. In addition to the activity of malate dehydrogenase, TCA metabolites could be replenished by a putative sodium pump oxaloacetate decarboxylase enzyme complex (Csac_2482 to Csac_2485).
Polysaccharide-degrading enzymes.
The capacity of C. saccharolyticus to hydrolyze a broad range of polysaccharides prior to fermentation differentiates this bacterium from many thermophilic anaerobes. Indeed, the genome of C. saccharolyticus encodes a wide range of enzymes active with carbohydrates (see Table S3 in the supplemental material). The genes encoding these carbohydrate-utilizing enzymes are often clustered on the chromosome, and the enzymes can be assigned to substrate-specific catabolic pathways for cellulose, hemicellulose, and, to a lesser extent, starch and pectin.
-1,4-Glucan polymers, for instance, can be transported into the cell using the maltodextrin ABC transport system proteins (Csac_0427 to Csac_0428/Csac_0431). An intracellular
-amylase (Csac_0426) and a 1,4-
-glucan phosphorylase (Csac_0429) further degrade the intracellular maltodextrins, releasing glucose-1-phosphate. Remarkably, a transcriptional regulator of the LacI family (Csac_0430) is also in this maltodextrin cluster and therefore is a good candidate for controlling expression of this maltodextrin-degrading pathway at the transcriptional level. In addition, a GCAAACGTTTGC consensus sequence was found in sequences upstream of genes in this transport cluster and genes encoding several starch-degrading enzymes, such as an
-amylase precursor (Csac_0408), an oligo-1,6-glucosidase (Csac_2428), a pullulanase (Csac_0689), a 4-
-glucanotransferase (Csac_0203), and a putative glucan 1,4-
-glucosidase (Csac_0130). The consensus sequence resembles the binding site (CGCAAACGTTTGCGT) of the maltose/maltodextrin transcriptional repressor MalR from the gram-positive organism Streptococcus pneumoniae (34). Besides this putative starch-degrading regulon, C. saccharolyticus has genes encoding a glycogen metabolic cluster (Csac_0780 to Csac_0784), a maltose ABC transport system (Csac_2491 to Csac_2493), and a second pullulanase (Csac_0671). Thus, C. saccharolyticus is well equipped for starch utilization.
An important feature of C. saccharolyticus is its ability to produce H2 not only from
-linked polymers but also from complex β-linked glycans, such as cellulose, hemicellulose, laminarin, and galactomannan. Growth on cellulosic substrates is rare among (hyper)thermophilic microorganisms. C. saccharolyticus does not metabolize cellulose by means of a cellulosome (45). For example, genes encoding the typical molecular components of a cellulosome (i.e., dockerin domains and scaffolding proteins) were not identified in the genome. Nevertheless, a gene cluster (Csac_1076 to Csac_1081) containing genes encoding cellulase precursors is present. The cellulases are potentially capable of degrading cellulose, a plant polysaccharide (6) (see Table S3 in the supplemental material). Moreover, the product of another gene cluster (Csac_1089 to Csac_1091) and an extracellular cellulase (Csac_0678) may assist in completely hydrolyzing cellulose to glucose. A remarkable aspect of the β-glycanases in Caldicellulosiruptor species is their bifunctional domain architecture consisting of central cellulose binding domains which are bordered by distinct catalytic domains (6). The C. saccharolyticus genome sequence confirmed the presence of genes encoding these multidomain hydrolases (Csac_1076 to Csac_1079 and Csac_2411), which comprise cellulase, mannanase, xylanase and arabinofuranosidase domains (6, 27, 30, 45). No genes encoding additional bifunctional β-glycanases were found in the genome. The bifunctional Csac_2411 gene is part of another large gene cluster (Csac_2404 to Csac_2411) coding for enzymes involved in hydrolysis of the plant polysaccharide xylan (hemicellulose). The genes encoding these mostly extracellular enzymes might be coexpressed with genes in a smaller putative xylan-utilizing cluster (Csac_0203 to Csac_0205). The genes in the latter cluster were not significantly upregulated during growth on xylose, in contrast to the genes in the former cluster. Furthermore, putative genes that encode enzymes that degrade galactomannan (Csac_0663 and Csac_0664), galactoarabinan (Csac_1560 to Csac_1562), and laminarin (Csac_2548) were identified.
The plant cell wall component pectin consists of an
-1,4-linked D-galacturonic acid backbone, sometimes with L-rhamnose interspersed, and side chains consisting of monosaccharides, such as D-galactose, D-xylose, and L-arabinose (38). Degradation of the main pectin component, D-galacturonate, requires a galacturonate isomerase, a tagaturonate reductase, and an altronate dehydratase to form 2-keto-3-deoxygluconate (KDG). Galacturonate isomerization may be catalyzed by glucuronate isomerase (Csac_1949). However, genes encoding tagaturonate reductase and altronate dehydratase were not detected in the genome of C. saccharolyticus. Apparently, novel enzymes or a novel pathway is responsible for the degradation of galacturonate. In contrast, a gene cluster encoding enzymes involved in the conversion of glucuronic acid to KDG (Csac_2686 to Csac_2689) was identified, and it includes genes encoding fructuronate reductase, mannonate dehydratase, a putative β-galactosidase/β-glucuronidase, and an
-glucuronidase. Glucuronic acid is a common substituent of xylan. Enzymes involved in the conversion of KDG to pyruvate and glyceraldehyde-3-phosphate, including KDG kinase (Csac_0355 or Csac_2720) and KDG-6-phosphate aldolase (Csac_0354), are present as well. The genes encoding enzymes involved in the last two steps are clustered with genes (Csac_0356 and Csac_0357, as well as Csac_2718 and Csac_2719) that encode enzymes involved in the conversion of 5-keto-4-deoxyuronate, an unsaturated cleavage product of pectate, to KDG. The genes encoding enzymes that are able to hydrolyze the pectate backbone and the side chains (e.g., unsaturated rhamnogalacturonyl hydrolase [Csac_0360], galacturan 1,4-
-galacturonidase [Csac_0361], β-galactosidase [Csac_0362], and a glycoside hydrolase with unknown substrate specificity [Csac_0363]) are close to the genes encoding the KDG metabolic enzymes as well. However, neither a pectate lyase gene nor a methylesterase gene could be definitively identified in the genome, although Csac_2721 and Csac_2728 might be candidates for a pectate lyase gene based on distant homology to known lyase genes.
C. saccharolyticus is also able to grow on L-rhamnose and on L-fucose, producing 1,2-propanediol as the end product (unpublished data). A putative rhamnose catabolic pathway was identified, in which generation of dihydroxyacetone phosphate and 1,2-propanediol is catalyzed by an L-rhamnose isomerase (Csac_0876), a putative L-rhamnulokinase (Csac_0989), an L-rhamnulose-1-phosphate aldolase (Csac_0865), and a putative lactaldehyde reductase (Csac_0407). Fucose can be processed by a similar pathway, using the lactaldehyde reductase mentioned above and as-yet-unidentified versions of L-fuculokinase, a bifunctional L-fucose isomerase/D-arabinose isomerase (Csac_1339), and fuculose-1-phosphate aldolase (Csac_0425).
Fermentation products.
Reducing equivalents are produced at the level of NAD and ferredoxin (Csac_0737). Since C. saccharolyticus can produce almost 4 mol of H2 per mol of glucose (14), both NADH and reduced ferredoxin should ultimately be able to transfer their reducing equivalents to protons to form hydrogen. In the genome, two hydrogenase gene clusters were identified, which are very similar to the two related clusters in T. tengcongensis (42). The first cluster includes genes encoding subunits of an Ni-Fe hydrogenase (EchA to EchF) (Csac_1534 to Csac_1539) and various genes required for maturation of the hydrogenase complex (HypA to HypF) (Csac_1540 to Csac_1545). In T. tengcongensis, the Ni-Fe hydrogenase is ferredoxin dependent and membrane bound and may act as a proton pump to generate a proton motive force. The second cluster (Csac_1860 to Csac_1864) codes for an Fe-only hydrogenase (HydA to HydD), which is NAD dependent and located in the cytoplasm, similar to the situation in T. tengcongensis (42). Hydrogenases that form H2 directly from NADH are unusual and make an NAD:ferredoxin oxidoreductase (Nfo) redundant. Nfos (also known as Rnfs) are membrane-bound multisubunit complexes that use or create an Na+ gradient coupled to the transfer of reducing equivalents between NADH and ferredoxin (7). An Nfo gene cluster has been identified in the genomes of C. thermocellum, T. maritima, and Thermotoga ethanolicus, but not in the genomes of T. tengcongensis and C. saccharolyticus. The absence of an Nfo in C. saccharolyticus also implies either that in the presence of elevated levels of H2 reduced ferredoxin may not be used to produce NADH or that a novel type of enzyme (complex) performs this reaction. Altogether, the information available suggests that C. saccharolyticus is able to produce hydrogen from ferredoxin but can also produce hydrogen directly from NADH. Production of hydrogen seems to be preferable, because under these conditions all pyruvate is converted to acetate (and CO2), which is coupled to the synthesis of ATP.
When the hydrogen partial pressure is too high, formation of hydrogen from NADH is no longer thermodynamically favorable. In such a case, NADH is oxidized through the formation of lactate or ethanol. A gene for a lactate dehydrogenase was identified (Csac_1027), but genes for acetaldehyde dehydrogenase and alcohol dehydrogenase were not obvious. In T. tengcongensis and T. ethanolicus, ethanol formation is NADPH dependent and catalyzed by a bifunctional alcohol dehydrogenase/acetyl-CoA thioesterase; this enzyme also has a homolog in C. saccharolyticus (Csac_0395).
A third small hydrogenase-like cluster was detected in the C. saccharolyticus genome and was composed of four genes encoding two NADH-binding proteins (Csac_0619 and Csac_0620), a molybdopterin oxidoreductase containing NAD and 4Fe-4S binding regions (Csac_0621), and an iron-containing alcohol dehydrogenase (Csac_0622). The function of this cluster is not known.
Transport systems.
As mentioned above, a number of genes encoding ABC transporters have been found in the C. saccharolyticus genome, including genes encoding the previously described carbohydrate-specific maltodextrin ABC transport system (Csac_0427, Csac_0428, and Csac_0431) and the maltose ABC transport system (Csac_2491 to Csac_2493). As observed for certain T. maritima maltose transporters, both sets of transport proteins lack ATP-binding subunits. In many bacteria, the intracellular ATPase used in the system is not encoded in the same operon. The genes encoding both of these sets of ABC transporters are located downstream from the genes encoding a two-component system consisting of a sensor histidine protein kinase and a response regulator. In C. saccharolyticus
50% of the genes encoding ABC carbohydrate transport systems are located near the genes encoding the two-component system on the chromosome.
Comparative analysis of C. saccharolyticus sugar binding proteins (SBPs) revealed that about two-thirds of them belong to the COG1653 category. This category includes the CUT1 subfamily members TM0432, TM0595, and TM1855, which transport a variety of di- and oligosaccharides, such as maltose. In genomes, the genes encoding more than one-half of COG1653 members are close to genes encoding glycoside hydrolases, supporting the identification of these proteins as ABC transporters involved in carbohydrate utilization. Putative SBPs encoded by Csac_0242, Csac_0391, Csac_2326, and Csac_2507 belong to COG1879. Csac_2506 and Csac_2510 are associated with the xylose transport-specific COG4213 category. As in T. maritima, a few putative SBPs (Csac_0261 and Csac_4166) are annotated as peptide transporters, although their actual functions are unknown. Components of phosphotransferase systems have been identified in C. saccharolyticus (although only one set of carbohydrate-specific EII has been identified), as have a few putative members of the major facilitator superfamily (Csac_0685, Csac_0786, Csac_1100, Csac1170, and Csac_2298). However, it is likely that carbohydrate utilization proceeds mainly through ABC transporters.
Transcriptional regulation.
The ability of C. saccharolyticus to utilize many different carbohydrates suggests that there is tight regulation in the pathways. The genes encoding many enzymes involved in carbohydrate utilization pathways appear to be regulated at the transcriptional level. Apart from the RNA polymerase core enzyme subunits (Csac_2259, Csac_2085, Csac_0951, and Csac_0952), this gram-positive species has 12 different
factors for construction of the RNA polymerase holoenzyme. In addition, multiple copies of genes encoding many sugar transcriptional regulators are present in the genome; for example, there are nine proteins belonging to the LacI family, six proteins belonging to the DeoR family, and eight proteins belonging to the GntR family, as well as 19 receiver proteins involved in a two-component system with a helix-turn-helix AraC domain. The genes encoding the latter proteins are always clustered with genes encoding sugar transporters and sugar hydrolytic enzymes.
Carbon catabolite repression (CCR) by glucose was not observed in C. saccharolyticus (Fig. 3). Nevertheless, genes encoding some indicators of the presence of carbon control protein A (CcpA)-dependent CCR of gram-positive bacteria are present in the genome, including (i) a gene encoding a histidine-containing phosphocarrier (HPr) (Csac_2438) that is close to the gene encoding the only PEP-dependent phosphotransferase system, which is fructose specific; (ii) a gene encoding an HPr(Ser) kinase (Csac_1186); and (iii) genes encoding a catabolite repression HPr (CrH) (Csac_1163) and nine members of the CcpA-containing LacI family. Binding sites for a putative CcpA, the catabolite-responsive element (cre), could not be identified. The Bacillus subtilis consensus sequence WWTGNAARCGNWWWCAWW (29), for instance, was detected only twice, once in the upstream region of the gene encoding the
-amylase precursor mentioned above (Csac_0408), where it overlaps the putative MalR binding site, and once in the middle of the gene encoding fumarate hydratase subunit
(Csac_2759). Nevertheless, CCR probably occurs in C. saccharolyticus, although the metabolite that induces this repression is not known.
![]() View larger version (19K): [in a new window] |
FIG. 3. Growth of C. saccharolyticus on a xylose-glucose mixture (1:1, wt/wt). , optical density at 660 nm (OD660); , hydrogen; , acetate; , glucose; , xylose; , lactate.
|
-linked glucan transcriptional regulator (Csac_0430), a regulator of the oxidative branch of the TCA cycle to oxoglutarate (Csac_0752), a repressor of L-arabinose metabolism (Csac_0722), and a putative response regulator receiver protein involved in glucuronate degradation (Csac_2690).
Transcriptome analysis of C. saccharolyticus.
One of the features of C. saccharolyticus that is beneficial for hydrogen production is its ability to degrade cellulosic substrates as well as hemicellulose. Moreover, mixtures of glucose and xylose can be fermented simultaneously (Fig. 3), suggesting that classical CCR by glucose does not occur. To elucidate the central carbon metabolic pathways and their regulation, a transcriptome analysis was performed after growth on glucose, xylose, and a 1:1 mixture of these two substrates. L-Rhamnose, which was likely to be metabolized by another pathway, was used as a reference substrate. The transcriptional data clearly show that glucose, xylose, and the glucose-xylose mixture all triggered upregulation of genes in the EM pathway, compared to the results for rhamnose (Fig. 4; see Table S4 in the supplemental material). In particular, fructose bisphosphate aldolase, glyceraldehyde-3-phosphate dehydrogenase, PPDK, and POR were significantly stimulated. The ultimate acetate-forming acetate kinase was also highly upregulated. A catabolic role for PPDK is intriguing, since this enzyme normally is associated with gluconeogenesis (as it is in propionic acid bacteria and plants) and PEP is usually converted by pyruvate kinase. However, homologs of PPDK are also present in related clostridia and Thermoanaerobacter species.
![]() View larger version (41K): [in a new window] |
FIG. 4. Intensity ratios for transcript levels of selected genes that responded to growth on glucose (open bars), xylose (gray bars), or a mixture of glucose and xylose (black bars), compared to the results obtained with rhamnose. Ratios are expressed as log2 values. Abbreviations: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; NdH, NAD-dependent hydrogenase; FdH, ferredoxin-dependent hydrogenase; dep., dependent; G6P, glucose-6-phosphate; FbisP, fructose bisphosphate; DH, dehydrogenase.
|
The transcriptional response to growth on monosaccharides enabled identification of genes and groups of adjacent genes (gene clusters) that were specifically upregulated in response to either glucose or xylose. On glucose, several genes coding for
-glucan hydrolases responded. The most striking observation, however, was the upregulation of an entire gene cluster (Csac_1991 to Csac_2000) involved in purine synthesis, which was not observed with xylose. Upregulation of purine biosynthesis genes was also detected in the transcriptome of Escherichia coli growing on glucose compared to the results obtained with xylose (16). On xylose, several gene clusters required for xylan or xylose conversion were upregulated (Csac0692 to Csac_0696, Csac0240 to Csac_0242, and Csac2416 to Csac_2419). These clusters encode ABC transport systems, transcriptional regulators, and endoxylanases. In addition, genes specifically required for growth on rhamnose were highly upregulated during growth on rhamnose, indicating the pathway for utilization of this sugar.
|
|
|---|
The absence of catabolite repression by glucose is an important characteristic for biohydrogen producers since it allows them to process an array of biomass-derived substrates simultaneously. Glucose is generally known to repress the use of xylose by CCR (18), but this was not observed in C. saccharolyticus. Moreover, the transcriptome showed that the various components of a CCR system encoded in the genome (CcpA homologs, HPr, HPr kinase) were not differentially transcribed under the conditions examined, suggesting that this type of regulation is not triggered by glucose or xylose. There were no obvious differences in the transcriptome for central metabolic pathways during growth on glucose, xylose, or a mixture of glucose and xylose. The EM pathway was not affected by the hexose or pentose substrate, which is in contrast to the results of a transcriptome analysis of E. coli for growth on the same substrates (16). Remarkably, the same specific ABC transporter is upregulated with both substrates, which is also consistent with the nonpreferential behavior of C. saccharolyticus with these two monomeric sugars.
Detailed knowledge of the metabolic pathways leading to hydrogen production enables workers to identify key enzymes that may be targets for improving the H2 yield by metabolic engineering. Currently, a genetic system for C. saccharolyticus is being developed, which initially will target the dehydrogenases involved in lactate and ethanol formation. Alternatively, genes could be introduced to form an oxidative branch of the PPP in order to obtain H2 yields greater than 4 mol per mol of glucose (50). In any case, the C. saccharolyticus genome provides new insights into the metabolic features of a versatile biohydrogen producer, which could inspire efforts to optimize microbial bioenergy systems.
This research was supported by Earth and Life Sciences Foundation and Chemische Wetenschappen project 050.50.206 (subsidized by The Netherlands Organization for Scientific Research) and by the European Commission (HYVOLUTION project reference no. 19825). Microarray fabrication was financially supported in part by Danisco (Genencor International, Inc., Palo Alto, CA). R.M.K. acknowledges support from the U.S. NSF Biotechnology Program. A.L.V. and D.L.L. acknowledge support from a U.S. Department of Education GAANN Fellowship.
Published ahead of print on 5 September 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
H.J.G.V.D.W., M.R.A.V., and A.L.V. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»