Characterizing acetogenic metabolism using a genome-scale metabolic reconstruction of Clostridium ljungdahlii

Background The metabolic capabilities of acetogens to ferment a wide range of sugars, to grow autotrophically on H2/CO2, and more importantly on synthesis gas (H2/CO/CO2) make them very attractive candidates as production hosts for biofuels and biocommodities. Acetogenic metabolism is considered one of the earliest modes of bacterial metabolism. A thorough understanding of various factors governing the metabolism, in particular energy conservation mechanisms, is critical for metabolic engineering of acetogens for targeted production of desired chemicals. Results Here, we present the genome-scale metabolic network of Clostridium ljungdahlii, the first such model for an acetogen. This genome-scale model (iHN637) consisting of 637 genes, 785 reactions, and 698 metabolites captures all the major central metabolic and biosynthetic pathways, in particular pathways involved in carbon fixation and energy conservation. A combination of metabolic modeling, with physiological and transcriptomic data provided insights into autotrophic metabolism as well as aided the characterization of a nitrate reduction pathway in C. ljungdahlii. Analysis of the iHN637 metabolic model revealed that flavin based electron bifurcation played a key role in energy conservation during autotrophic growth and helped identify genes for some of the critical steps in this mechanism. Conclusions iHN637 represents a predictive model that recapitulates experimental data, and provides valuable insights into the metabolic response of C. ljungdahlii to genetic perturbations under various growth conditions. Thus, the model will be instrumental in guiding metabolic engineering of C. ljungdahlii for the industrial production of biocommodities and biofuels.


Background
Acetogenic microorganisms have unique metabolic capabilities that, if understood, could be harnessed to greatly increase strain engineering design options for microbial production of biofuels and biocommodities. Acetogens were discovered for their ability to autotrophically reduce CO 2 to acetate and conserve energy simultaneously using the Wood-Ljungdahl pathway [1]. In addition to reducing CO 2 , acetogens can ferment a wide variety of sugars as well. Their ability for autptrophy allows them also to grow on synthesis gas (H 2 /CO/CO 2 ) by utilizing either H 2 /CO 2 or CO alone. More recently, acetogens like Sporomusa ovata and Clostridium ljungdahlii have been shown to be capable of another form of autotrophic metabolism, called microbial electrosynthesis [2]. Microbial electrosynthesis is a process in which microorganisms directly use electric current to reduce carbon dioxide to multi-carbon organic compounds that are excreted from the cells into extracellular medium [3]. These discoveries further expand the range of economically viable feedstocks that can be used for industrial production of biofuels and biochemicals [4].
However, to efficiently engineer acetogens into platform strains and production hosts for chemicals at an industrial scale, a thorough understanding of the metabolic capabilities and aspects of energy conservation is a necessary prerequisite. The Wood-Ljungdahl pathway of carbon fixation employed by acetogens is believed to be one of the most ancient metabolic pathways [5]. Physiological and biochemical aspects governing this metabolic capability have been poorly characterized with several fundamental discoveries being made only recently [6,7]. One of the fundamental characteristics of acetogenic metabolism discovered recently is the concept of flavin-based electron bifurcation. In this mechanism, there is a concomitant coupling of an endergonic redox reaction with the oxidation of the same electron donor with higher potential electron acceptors. This feature is believed to play a key role in the energy conservation mechanisms of acetogens [8].
Constraints based reconstruction and analysis (COBRA) has been a powerful technique for discovering and understanding new capabilities and content in microorganisms, as well as in guiding metabolic engineering efforts for targeted production [9]. The COBRA approach relies on a genome-scale metabolic network reconstruction, which enumerates the metabolic transformations and the genes encoding them in a mathematical format. This reconstructed network together with physiological data can then enable the prediction of the functionality of an organism under conditions of interest. Such a validated and accurate network can be utilized for prospective design and engineering of cellular networks [10]. Furthermore, genome-scale metabolic networks of bacteria and constraints-based modeling have been instrumental in guiding metabolic engineering at an industrial scale [11].
In this study, we reconstructed the first genome-scale metabolic network of an acetogen, Clostridium ljungdahlii. We characterized the metabolic phenotypes of this bacterium under heterotrophic and autotrophic growth conditions. We further utilized physiological and transcriptomic data to elucidate novel biological capabilities and aspects of energy conservation during autotrophic metabolism of C. ljungdahlii.

Results and discussion
Genome-scale reconstruction of the acetogen Clostridium ljungdahlii The genome-scale metabolic network for C. ljungdahlii was reconstructed using a four-step integrative reconciliatory workflow involving four published models of related clostridia species and two draft models ( Figure 1). The first draft metabolic model was generated based on the C. ljungdahlii genome annotation [12] using the AutoModel functionality of SimPheny (Genomatica, San Diego), while the second draft model was generated using the Model-SEED database [13]. In addition to these two draft models, homologs to C. ljungdahlii genes were identified in published genome-scale reconstructions of related clostridia species (C. acetobutylicum, C. thermocellum, and C. beijerenckii) [14][15][16][17] using the Smith-Waterman alignment. A 60% amino acid sequence identity cutoff was used to identify C. ljungdahlii homologs in the other clostridia genomes. The reactions corresponding to these genes in the respective clostridia models were compiled and reconciled with the two draft models. The list of discrepancies pertaining to nomenclature between the different databases and gene-reaction associations was manually curated with the aid of biochemical literature and databases such as KEGG [18] and SEED [13]. Manual evaluation of new content from the annotation and existing genome-scale reconstructions consisted of gathering genetic, biochemical, sequence, and physiological data and reconciling this information to determine the likelihood of each reaction being present in the organism. The curated reconstruction was evaluated for functional performance with the aid of a biomass objective function that was formulated using an existing template (see Methods). Using inference based on pathway function, as well as the SMILEY computational algorithm [19,20], which predicts reactions that fill gaps in a metabolic network, reaction content was added to the network so that it could produce the necessary biomass components. This resulted in a final network (iHN637) consisting 637 genes, 785 reactions, and 698 metabolites. This iHN637 reconstruction represents the first genome-scale metabolic model of an acetogen.

Functional testing of the model
The functional capabilities of the iHN637 reconstruction were evaluated by simulating growth under different conditions using the developed biomass objective function (Methods). iHN637 was successfully able to simulate heterotrophic growth on several of the known substrates [21] (fructose, glucose, gluconate, arabinose, ribose, xylose, ethanol, formate, serine, citrulline, pyruvate, arginine, aspartate, and glutamate). Moreover, autotrophic growth with H 2 and CO 2 as well as just CO alone was feasible using the iHN637 reconstruction. The predicted growth rates and acetate production rates for each of these various substrates are provided in Table 1. In addition to acetate and ethanol, the iHN637 metabolic model is also capable of producing lactate and 2,3-butanediol using native pathways as described elsewhere [22]. Furthermore, the model was able to successfully predict a growth rate consistent with experimental observations for heterotrophic growth on fructose. When the iHN637 model was constrained with experimentally determined uptake rates of fructose (1.8 mmol/gDWh), the growth rate predicted by the model (0.077 h −1 ) was in good agreement with the experimentally measured growth rates (0.072 h −1 ). The simulated acetate production rate (2.7 mmol/gDWh) was also consistent with the measured value of 2.85 mmol/gDWh (Additional file 1: Figure S1).

Characterization of a nitrate reduction pathway in Clostridium ljungdahlii
A novel nitrate reduction pathway for clostridia was reconstructed in the iHN637 model based on the C. ljungdahlii genome annotation. In addition to the glutamine syn-thetase and glutamine:2-oxoglutarate aminotransferase or the fixation of dinitrogen using a molybdenum-dependent nitrogenase, this pathway could serve as a third possible nitrogen assimilation route in C. ljungdahlii. This pathway presents characteristics of both, assimilatory and respiratory nitrate reductases and the reduction of nitrate to ammonia via hydroxylamine resembles the nitrate respiratory system recently identified in Nautilia profundicola [23]. The nitrate reduction pathway consists of a soluble nitrate reductase (NTRARf), a nitrite reductase (NTRIR5), and an additional hydroxylamine reductase (HAMR) (Figure 2A). Transcriptomic profiling of C. ljungdahlii grown heterotrophically on fructose in medium containing ammonium or ammonium-free medium with nitrate as a nitrogen source, helped validate this pathway. The genes encoding for the nitrate reduction pathway were all significantly upregulated (>8-fold) when nitrate was provided in ammonium-free medium. Specifically, the genes encoding for nitrate reductase, nitrite reductase, and the nitrate transporter were upregulated greater than 400-fold ( Figure 2A). Growth experiments showed that nitrate was consumed with stoichiometric production of ammonium ( Figure 2B), further validating the proposed pathway. Furthermore, modeling simulations predict a 15% reduction in acetate flux when nitrate served as nitrogen source due to nitrate reduction acting as an additional electron sink (Additional file 2: Table S1). The diversion of electrons seems to impact the acetate production gained through fixation of CO 2 using electrons obtained from fermenting fructose. In addition, nitrate reduction has been thought of as a primitive mode of energy metabolism in certain clostridia with potential roles in balancing electrons and affecting growth yields during fermentation [24]. The identification of this pathway provides the opportunity for a detailed investigation on its potential functional effects on the fermentation of C. ljungdahlii.

Reconstructing the carbon fixation pathway in C. ljungdahlii
The C. ljungdahlii genome encodes for all the genes involved in the Wood-Ljungdahl pathway of carbon fixation [12]. Accordingly, this pathway was incorporated into the iHN637 reconstruction. The pathway consists of two different branches, one CO 2 contributing to the methyl group of acetyl CoA via the eastern branch and the other CO 2 contributing to the carbonyl group of the acetyl CoA via the western branch [25]. The key enzyme in this pathway is the CODH/ACS complex, which performs the dual activity of reducing CO 2 to CO via the carbon monoxide dehydrogenase (CODH) activity and the subsequent formation of acetyl CoA through the acetyl-CoA synthase (ACS). It has been suggested that the CODH/ACS complex functions in a manner that CO resulting from CODH activity is kept as a bound metabolite in the complex to diminish thermodynamic barriers involved in this energy intensive process [26]. Furthermore, the genome of C. ljungdahlii suggested that CO oxidation to CO 2 is likely to be catalyzed by a different carbon monoxide dehydrogenase gene (CLJU c09090-09110) other than the ones encoded by the CODH/ACS complex. Taking these into consideration, the iHN637 reconstruction incorporates the CODH/ACS reaction as the net reaction of carbon monoxide dehydrogenase and acetyl CoA synthase activity, and associates the gene cluster CLJU c09090-09110 to the reaction representing CO oxidation to CO 2 . The methylenetetrahydrofolate dehydrogenase reaction catalyzed by the bifunctional FolCD (CLJU_c37630) was assumed to be NADPH-dependent based on sequence similarity with the corresponding enzyme in the acetogen Moorella thermoacetica [25,27]. The energy conservation mechanisms associated with autotrophic growth using the Wood-Ljungdahl pathway is discussed in detail in the following section.

Analysis of energy conservation steps in C. ljungdahlii
Acetogens have long been thought to be living at the thermodynamic limit due to the energy requirements of the Wood-Ljungdahl pathway [7]. Given that the ATP generated from acetate production is required for the activation of formate, it has been proposed that acetogens should have additional energy conservation mechanisms [7]. Acetogens have generally been classified into those containing respiratory cytochromes that establish a proton gradient (e.g., M. thermoacetica) and those that do not. Organisms such as Acetobacterium woodii have been shown to employ a sodium gradient as an energy conservation mechanism [7]. The C. ljungdahlii genome sequence has revealed that it falls into a third class of acetogens that neither uses respiratory cytochromes nor a sodium gradient for energy conservation [12]. As reviewed recently, flavin-based electron bifurcation is expected to serve as an alternative mechanism of energy conservation in acetogens [8]. Analysis of the iHN637 reconstruction of C. ljungdahlii reveals that electron bifurcation and proton translocating ferredoxin oxidation are critical mechanisms for energy conservation during autotrophic growth in C. ljungdahlii (Figure 3). The iHN637 model predicted proton translocation by the membranebound Rnf complex while oxidizing ferredoxin and reducing NAD was essential for autotrophic growth. Specifically, the proton gradient established by the Rnf complex is essential for driving ATP generation through ATP synthase ( Figure 3B). This prediction is consistent with recent Figure 2 Characterizing the nitrate reduction pathway in C. ljungdahlii. (A) Pathway map of the reconstructed nitrate reduction pathway with the respective genes associated to the reactions. Inset shows the differential expression of the genes in the nitrate reduction pathway under ammonium-free medium with nitrate compared to nitrate-free medium with ammonium as the nitrogen source. Reaction abbreviations: NO3t2 (nitrate transporter), NO2t2r (nitrite transporter), NTRARf (nitrate reductase), NTRIR5 (nitrite reductase), HAMR (hydroxylamine reductase), NH4t (ammonium transporter). (B) Physiological data when C. ljungdahlii is grown on fructose in ammonium-free media with nitrate as the nitrogen source. Purple line represents fructose, green line represents OD 600 , brown line represents acetate, red line represents nitrate and orange line represents ammonia measurements respectively.
experimental observations of the lack of autotrophic growth in a strain where the Rnf complex was deleted [28]. An electron-bifurcating transhydrogenase (Nfn complex) that couples exergonic reduction of NADP with ferredoxin to drive the endergonic reduction of NADP with NADH had been characterized in C. kluyveri [29]. Recently, a similar electron bifurcating NADPH transhydrogenase activity was also reported in the acetogen M.
thermoacetica [30]. A homolog to the nfnA and nfnB genes in C. kluyveri was identified in the C. ljungdahlii genome. This gene (CLJU_c37240) was previously annotated as glutamate synthase. Another recent report describes a similar electron-bifurcating NADPH transhydrogenase activity in C. autoethanogenum and further indicates the presence of a homolog in C. ljungdahlii [31]. As part of iHN637, this gene (CLJU_c37240) has  been assigned the electron bifurcating transhydrogenase reaction ( Figure 3C). Simulatons further reveal that this reaction is essential for the interconversion between NADPH and NADH during autotrophic growth.
It was speculated that the highly exergonic reduction of methylene-tetrahydrofolate by NADH could be a site for electron bifurcation to generate additional reduced ferredoxin from NADH [12]. Investigation of the A. woodii genome revealed that the small subunit of the methylenetetrahydrofolate reductase (MetV) can act as a flavoprotein that could aid in this proposed electron bifurcation, thereby indirectly establishing a proton gradient through the Rnf complex [7]. A homolog to the MetV gene from A. woodii was identified next to the MetF gene in C. ljungdahlii and investigation of the intergenic spacing in the genome revealed a possible organization of MetV and MetF in the same operon (CLJU c37610 and CLJU c37620). Hence, an electron bifurcating methylenetetrahydrofolate reductase reaction was included in the iHN637 reconstruction ( Figure 3A).

Modeling of heterotrophic and autotrophic metabolism
Transcriptome profiling of C. ljungdahlii when grown on fructose and autotrophically on H 2 /CO 2 identified 114 differentially expressed genes (50 upregulated and 64 downregulated in autotrophic conditions) (Additional file 3: Table S2). Among the key genes that were downregulated during autotrophic growth were genes involved in fructose metabolism, pyrimidine biosynthesis, and amino acid biosynthesis. The downregulation of biosynthetic genes in autotrophic conditions is likely due to the slower growth rate during autotrophic growth compared to heterotrophic growth. Among the key genes upregulated was the cooS1 (CODH) gene. The cooS1 gene in Carboxydothermus hydrogenoformans is hypothesized to play a role in facilitating electron transfer through the hydrogenase [32]. Hence it is possible that this gene in C. ljungdahlii has a similar role while growing on hydrogen. Importantly, transcriptome profiling during autotrophic growth revealed that the gene annotated as metV was co-transcribed along with metF, further validating the annotation and assignment of its function for electron bifurcation with metF (Additional file 4: Figure S2). It must also be noted that genes encoding the Wood-Ljungdahl pathway are not differentially expressed between the conditions examined in this study (heterotrophic growth on fructose versus autotrophic growth on H 2 /CO 2 ). This is probably due to the fact that the Wood-Ljungdahl pathway is employed for fixation of CO 2 using electrons obtained from fermenting sugars [1]. However, a recent transcriptomic study on syngas utilization in C. ljungdahlii showed upregulation of the genes encoding the Wood-Ljungdahl pathway under autotrophic growth (4:1 CO:CO 2 ) [33]. This difference can be attributed to the differences in the autotrophic conditions investigated in both these studies with one being CO 2 /H 2 and the other being a condition rich in CO.
The genome-scale model of C. ljungdahlii has revealed that flavin-based electron bifurcation is found to play a critical role in the essential reactions for autotrophic growth, thereby providing an explicit account of the requirement of reducing equivalents for the fixation of CO 2 to acetyl CoA ( Figure 4A). Another instance of electron bifurcation in the autotrophic growth of C. ljungdahlii was observed in the first step of H 2 activation to reduce ferredoxin. A bifurcating hydrogenase (HydABC), which couples the exergonic reduction of NADH from H 2 to drive the endergonic reduction of ferredoxin from H 2 has been identified to perform this process in A. woodii [7,34]. A homolog to this gene cluster has been identified in C. ljungdahlii ( Figure 3D). However, in another recent report, a similar electron-bifurcating hydrogenase that is specific to NADPH has been identified and characterized in C. autoethanogenum, an acetogen closely related to C. ljungdahlii [31]. This report also mentions that the hydrogenase of C. ljungdahlii is homologous to the NADPspecific hydrogenase in C. autoethanogenum. Since there is no definitive biochemical data available for C. ljungdahlii the metabolic reconstruction (iHN637) has both the NAD-specific and NADP-specific electron bifurcating hydrogenases associated to this gene cluster, with the NADP-specific one being set as the default.
The impact of this cofactor specificity on the metabolic network and the importance of detailing these energy conservation mechanisms was further highlighted when genetic perturbations were modeled under heterotrophic and autotrophic conditions. Specifically, the effect of knocking out acetate kinase (ackA) was modeled under heterotrophic growth on fructose and autotrophic growth under the major constituents of syngas (CO alone as well as H 2 /CO 2 ) ( Figure 5). The deletion of ackA is one of the primary targets for optimizing ethanol production in acetogens.
Modeling simulations revealed that the activity of acetate kinase was not essential for growth on fructose. As expected, the lost ATP from the acetate kinase reaction adversely impacted the growth resulting in a 20% reduction in growth yield with the model predicting a diversion of flux to ethanol ( Figure 5A).
When ethanol is the desired product from carbon fixation through the Wood-Ljungdahl pathway, 2 reducing equivalents in the form of NADH are required in addition to the 4 reducing equivalents (1 NADPH, 1 ferredoxin, and 2 NADH) and 1 ATP ( Figure 4B). Since ethanol production from syngas is an industrially relevant process, iHN637 was used to simulate the effect of this important deletion (ackA) under two different autotrophic growth conditions (CO alone as well as H 2 /CO 2 ) that reflect the composition of syngas. Simulations of the acetate kinase knock out under these conditions showed that the essentiality of ackA gene was dependent on the source of electrons (CO or H 2 ) and the cofactor specificity of the bifurcating hydrogenase. When CO is the sole electron donor, all 6 required reducing equivalents are obtained solely in the form of low-potential reduced ferredoxin. Hence, flux is driven through the Rnf complex to get the necessary NADH for the Wood-Ljungdahl pathway and ethanol formation. Since the Rnf complex couples proton translocation with the oxidation of ferredoxin, this results in the formation of the required ATP to fix CO 2 using the Wood-Ljungdahl pathway ( Figure 5B).
When H 2 is the electron donor, and the bifurcating hydrogenase is NADP dependent, the 6 reducing equivalents are obtained in an equimolar split of NADPH and ferredoxin ( Figure 5C). The Nfn complex ensures conversion of 2 NADPH to 1 NADH and 1 reduced ferredoxin. However, to get the required balance of ferredoxin and NADH, there is a requirement of flux through the proton translocating Rnf complex that interconverts ferredoxin and NADH. While converting 3 ferredoxin to 3 NADH, the proton motive force at the Rnf complex enables the generation of 1.5 ATP. This ensures that carbon fixation and ethanol production is feasible in the absence of ackA. However, if the bifurcating hydrogenase is NADH dependent, the required 6 reducing equivalents from this reaction are obtained in the form of 3 NADH and 3 ferredoxin ( Figure 5D). The Nfn complex and Rnf complex ensure the interconversion of these reducing equivalents to obtain the necessary NADPH and NADH. However, due to the lower demand for NADH in this scenario, the flux through the proton-translocating Rnf complex is reduced, resulting in the generation of only 0.75 ATP.
This ATP generation due to the lower flux through Rnf complex explains the infeasibility of autotrophic growth of the ackA mutant on H 2 /CO 2 if the cofactor specificity of the bifurcating hydrogenase is NADH. Similarly, since H 2 is a major component of syngas along with CO, this impact of the cofactor specificity of the bifurcating hydrogenase on energy conservation during autotrophic growth is fundamental to engineering acetogens for the production of desired compounds from syngas.

Conclusions
This study presents the first genome-scale metabolic network for an acetogen. The model is predictive of the metabolic capabilities of C. ljungdahlii and together with the recently developed genetic system [35] can aid in the characterization of its metabolic phenotype and guide engineering strategies. For instance, C. ljungdahlii is an important organism in the conversion of synthesis gas. Eliminating acetate production is typically the first choice when engineering acetogens like C. ljungdahlii for the production of chemicals like ethanol. Our results highlight the genetic and energetic constraints while rerouting production of target molecules away from acetate. Therefore, engineering strategies have to account for the appropriate cofactor dependencies and the specific energy conservation mechanisms employed depending on the electron source. Furthermore, C. ljungdahlii is also capable of microbial electrosynthesis, thereby fixing CO 2 using electrons directly provided by an electrode [2]. This metabolic model can be extended to obtain insights into potential electron transfer components involved in microbial electrosynthesis. Modeling microbial electrosynthesis would require a systematic evaluation of the different electron transfer pathways through potential carriers such as cytochromes, quinones, or flavins. So far none of these redox active molecules have been identified in C. ljungdahlii. Preliminary studies using a quinone extraction procedure and LC-MS have hinted at the presence of a menaquinone-like derivative in C. ljungdahlii. While beyond the scope of this study, a thorough investigation of the structure of this molecule and its potential role in electron transfer could significantly aid modeling and characterization of microbial electrosynthesis. The model could also serve as the basis for gaining further insight into the metabolism of other acetogens. Given the different energy conservation mechanisms adopted by acetogens like M. thermoacetica, the model enables a systematic comparison of the functional impact of these different mechanisms on acetogenic metabolism in general. Overall, these representative examples emphasize how the detailed metabolic reconstruction of an acetogen can aid in the discovery of new function and guiding strain-design strategies for advanced bioproduction from waste streams.

Reconstruction process
The genome-scale metabolic network for C. ljungdahlii was reconstructed using a four-step integrative reconciliatory workflow involving four published models of related Clostridia species and two draft models. The first draft metabolic model was generated using the ModelSEED database [13]. This draft model is generated automatically in a high-throughput manner by comparing the existing genome annotations in the SEED database. To complement this, we generated another automated draft model that relies on a database of curated genome-scale models. This was done using the AutoModel functionality of SimPheny (Genomatica, San Diego). In addition to these two draft models, homologs to C. ljungdahlii genes were identified in published genome-scale reconstructions of related Clostridia species (C. acetobutylicum [14,15], C. thermocellum [16], and C. beijerenckii [17]) using the Smith-Waterman alignment. A 60% amino acid sequence identity was employed to identify C. ljungdahlii homologs in the other Clostridia genomes. The final curated reconstruction was built in the SimPheny platform (Genomatica, San Diego) to ensure that the reactions in the network were all elementally and charge balanced.

Generation of the biomass objective function
The Biomass Objective Functions for C. ljungdahlii was formulated using a previous template [36]. The protein content was determined to be 43% and for the rest of the macromolecular breakdown, the biomass content previously determined for other Gram-positive bacteria, B. subtilis and C. beijerenckii, was used [17,37]. However, for modeling teichoic acid composition, the distribution based on S. aureus was used due to similarity in terms of low G + C content [38]. It should be noted that prediction of growth rate and unmeasured uptake rates are relatively insensitive to realistic variations in biomass macromolecular weight fractions [39].

Flux balance analysis simulations
The reconstructed metabolic network was represented in a mathematical format in a stoichiometric matrix S, where the rows correspond to the metabolites and columns correspond to the reactions in the network. Flux Balance Analysis simulations were carried out as described previously [40] using the COBRA Toolbox [41] and the SimPheny framework (Genomatica, Inc., San Diego, CA) was used for simulations. The objective used in the simulations was maximizing growth through the biomass objective function. In silico simulations of gene knockout was carried out using the "singleGeneDeletion" function of the COBRA toolbox [41]. The metabolic models for simulating heterotrophic and autotrophic growth are available for download (Additional file 5: Model 1 and Additional file 6: Model 2).

Bacterial growth conditions and chemical analysis
Clostridium ljungdahlii (ATCC 55383) was grown in 125 mL serum bottles under anaerobic conditions containing 100 mL of PETC medium (ATCC medium 1754) at 37°C. For growth experiments with nitrate, 1 mL of an anoxic NaNO 3 stock solution (1 M) was added to the medium. Ammonium was omitted in some of the experiments from the medium containing nitrate. Growth was routinely determined by measurement of the OD 600 . Hydrogen and other gases in the cultures were quantified using gas chromatography as previously described [42]. Concentrations of fructose, acetate, and ethanol were determined by high-performance liquid chromatography (Waters) as previously described [43]. Detection was performed by UV absorption at 410 nm. Nitrate, nitrite, and ammonia concentrations were analyzed using an automated segmented flow analyzer (AA3 HR, SEAL Analytical). Ammonium concentrations were determined by flow injection analysis modified for small sample volumes [44]. The sum of nitrate and nitrite was determined spectrophotometrically after reduction of samples with cadmium [45]. The procedure was identical for determination of nitrite but without the cadmium reduction step.
Transcriptomic profiling of C. ljungdahlii The transcriptome of C. ljungdahlii was profiled under three different growth conditions: heterotrophically on fructose, heterotrophically on fructose with nitrate as the sole nitrogen source, and autotrophically on H 2 /CO 2 . For each of these conditions, total RNA was isolated from C. ljungdahlii cells growing in mid-log phase (OD 600 of~0.4 for both the heterotrophic conditions and an OD 600 of~0.09 for the autotrophic condition on H 2 /CO 2 ) using the QIAGEN RNeasy Mini kit with on-column DNase I (QIAGEN) treatment. rRNA was subtracted from 2.5 μg of total RNA using the Gram-Positive Ribo-Zero rRNA removal kit (Epicentre). Then, paired-end, strand specific RNA sequencing (RNA-seq) was performed using the dUTP method [46,47] with the following changes. The rRNA subtracted RNA was fragmented with RNA fragmentation reagents (Ambion) for 3 minutes at 70°C. To synthesize first strand cDNA, random hexamer primers were used (Invitrogen). Final libraries were assayed for quality using High Sensitivity DNA chips on a Bioanalyzer (Agilent) and quantified using a Qubit (Invitrogen). A total of 12.5 pM was loaded onto a MiSeq (Illumina) to yield 31 bp paired-end reads.
The RNA-seq reads were aligned to the genome sequence of C. ljungdahlii (RefSeq: NC_014328.) using Bowtie [48] with two mismatches allowed per read alignment.