Clostridium acetobutylicum possesses two homologous adhE genes, adhE1 and adhE2, which have been proposed to be responsible for butanol production in solventogenic and alcohologenic cultures, respectively. To investigate their contributions in detail, in-frame deletion mutants of each gene were constructed and subjected to quantitative transcriptomic (mRNA molecules/cell) and fluxomic analyses in acidogenic, solventogenic, and alcohologenic chemostat cultures.
Under solventogenesis, compared to the control strain, only ΔadhE1 mutant exhibited significant changes showing decreased butanol production and transcriptional expression changes in numerous genes. In particular, adhE2 was over expressed (126-fold); thus, AdhE2 can partially replace AdhE1 for butanol production (more than 30 % of the in vivo butanol flux) under solventogenesis. Under alcohologenesis, only ΔadhE2 mutant exhibited striking changes in gene expression and metabolic fluxes, and butanol production was completely lost. Therefore, it was demonstrated that AdhE2 is essential for butanol production and thus metabolic fluxes were redirected toward butyrate formation. Under acidogenesis, metabolic fluxes were not significantly changed in both mutants except the complete loss of butanol formation in ΔadhE2, but numerous changes in gene expression were observed. Furthermore, most of the significantly up- or down-regulated genes under this condition showed the same pattern of change in both mutants.
This quantitative system-scale analysis confirms the proposed roles of AdhE1 and AdhE2 in butanol formation that AdhE1 is the key enzyme under solventogenesis, whereas AdhE2 is the key enzyme for butanol formation under acidogenesis and alcohologenesis. Our study also highlights the metabolic flexibility of C. acetobutylicum to genetic alterations of its primary metabolism.
Clostridium acetobutylicum is now considered as the model organism for the study of solventogenic Clostridia [1, 2]. The superiority of butanol over ethanol as an alternative biofuel has attracted research interest into C. acetobutylicum and other recombinant bacteria producing butanol as major products .
In phosphate-limited chemostat cultures, C. acetobutylicum can be maintained in three different stable metabolic states [4–8] without cellular differentiation : acidogenic (producing acetate and butyrate) when grown at neutral pH with glucose; solventogenic (producing acetone, butanol, and ethanol) when grown at low pH with glucose; and alcohologenic (forming butanol and ethanol but not acetone) when grown at neutral pH under conditions of high NAD(P)H availability [5, 6, 10].
AdhE1 (CA_P0162 gene product, also referred to as Aad) has long been considered as an NADH-dependent bifunctional alcohol/aldehyde dehydrogenase responsible for alcohol formation in solventogenic C. acetobutylicum cultures [1, 2, 11]. Recently, however, AdhE1 was purified and shown to have lost most of its alcohol dehydrogenase activity despite its NADH-dependent aldehyde dehydrogenase activity .
Prior to the identification of adhE2 (CA_P0035), the existence of alcohologenesis-specific gene(s) responsible for alcohol formation was predicted because (i) there was high NADH-dependent butanol dehydrogenase activity in alcohologenesis versus high NADPH-dependent butanol dehydrogenase activity in solventogenesis [5, 7] and (ii) previously identified genes related to butanol production (bdhA, bdhB, adhE1) were not induced in alcohologenic cultures . The adhE2 gene is the second aldehyde/alcohol dehydrogenase-encoding gene and is carried by the pSol1 megaplasmid, as is adhE1 . The two genes are not clustered, in contrast to the observations for C. ljungdahlii  and their expression patterns differ [9, 12]. adhE1, ctfA, and ctfB (CA_P0163 and CA_P0164) form the sol operon [1, 11]; ctfA and ctfB encode the CoA-transferase responsible for the first step of acetone formation, while the second step, catalyzed by acetoacetate decarboxylase, is encoded by adc (CA_P0165), located downstream of the sol operon. However, adc is transcribed under the control of its own promoter, which is oriented in the opposite direction of the sol operon .
In the three metabolic states, the contributions of the different enzymes responsible for the butyraldehyde dehydrogenase and butanol dehydrogenase activities to butanol flux have recently been characterized . Under acidogenesis, the low butanol flux is catalyzed by AdhE2 (100 %) for butyraldehyde dehydrogenase activity, while BdhB and BdhA are responsible for butanol dehydrogenase activity. Under solventogenesis, AdhE1 (95 %; the other 5 % is contributed by AdhE2) is the key player responsible for butyraldehyde dehydrogenase activity, while BdhB, BdhA, and BdhC are responsible for butanol dehydrogenase activity. Under alcohologenesis, AdhE2 plays a major role in both butyraldehyde dehydrogenase (100 %) and butanol dehydrogenase activities. In the study of Cooksley et al. , adhE1 and adhE2 knockout mutants were (i) constructed using the ClosTron method  and (ii) phenotypically characterized in batch culture using Clostridium basal medium (CBMS) without pH adjustment. The adhE1 knockout mutant obtained in their study exhibited low ethanol and no butanol formation along with scant acetone production; these findings were consistent with the polar effect of the intron on ctfAB transcription . Using the adhE2 knockout mutant, no alteration of solvent production was observed; however, the adhE2 knockout mutant has not been evaluated under alcohologenic conditions, under which it is normally thought to play a major role .
The aim of this study was to perform clean individual in-frame deletions of adhE1 and adhE2 to characterize their roles in butanol formation in the three different metabolic states in more detail. Furthermore, to study the metabolic flexibility of C. acetobutylicum in response to each of these gene deletions, a complete fluxomic and quantitative transcriptomic analysis was also performed in the three conditions known for the wild-type strains: acidogenic, solventogenic, and alcohologenic states. The results presented here not only support our previous studies [12, 14] on the roles of AdhE1 and AdhE2 in butanol formation in different metabolic states but also highlight the metabolic flexibility of C. acetobutylicum to genetically alter its primary metabolism.
Results and discussion
Construction of ΔadhE1 and ΔadhE2 mutant strains
Construction of the ΔadhE2 mutant was relatively straightforward, as adhE2 is expressed in a monocistronic operon  (Fig. 1a). However, the position of adhE1 as the first gene of the sol operon made the construction of ΔadhE1 more complicated because the transcription of downstream ctfAB genes could be affected. Figure 1b–d shows different configurations of the sol operon promoter, ctfAB genes, and either catP cassette with two FRT (Flippase Recognition Target) sites or a single FRT site remaining after Flippase (Flp)-FRT recombination of the three different types of ΔadhE1 mutants generated in this study. The first constructed ΔadhE1 mutant, ΔCA_C1502ΔuppΔadhE1::catP (Fig. 1b), was unable to form acetone as predicted because a transcriptional terminator was included in the catP cassette, which is located upstream of ctfAB encoding the acetoacetyl coenzyme A:acetate/butyrate:coenzyme A transferase that is responsible for the first specific step of acetone formation . However, after removing the catP cassette from ΔCA_C1502ΔuppΔadhE1::catP, acetone production was unexpectedly not recovered in ΔCA_C1502ΔuppΔadhE1 (Fig. 1c). The presence of the megaplasmid pSOL1 was confirmed by the production of ethanol and butanol under alcohologenic conditions and was attributed to adhE2 expression. By sequencing the pSOL1 region around the adhE1 deletion, we confirmed that there was no mutation in the sol promoter, ctfAB and adc (encoding acetoacetate decarboxylase, which is responsible for the last step of acetone production). Based on these results, the possibility of unsuspected early transcriptional termination by the FRT site remaining after catP removal was deduced. To confirm the early termination of transcription by an FRT site and to eliminate this polar effect on acetone production, a new plasmid was constructed to position both of the FRT sites carried by the catP cassette upstream of the sol operon promoter and was used to construct the ΔadhE1 mutant ΔCA_C1502ΔuppΔadhE1::catP-A1A4 mutant (Fig. 1d). Consistent with our hypothesis, this last ΔadhE1 mutant recovered acetone production (Fig. 2, Additional file 1: Fig. S3). To the best of our knowledge, the potential role of an FRT site as a transcriptional terminator was reported once in Salmonella  and twice in yeast [19, 20], although the FRT site is not generally recognized as possessing this additional activity. However, the high score of the FRT site hit from the “Dimers and Hairpin Loops analysis” in Vector NTI software (Invitrogen) and the detection of this activity upon deleting adhE1 in C. acetobutylicum unambiguously demonstrate that the FRT site can function as a transcriptional terminator.
Hereafter, C. acetobutylicum ΔCA_C1502ΔuppΔadhE1::catP-A1A4 (Fig. 1d) is referred to as ΔadhE1 in all the chemostat culture experiments.
Carbon and electron fluxes of ΔadhE1 and ΔadhE2 mutants under different physiological conditions
The ΔadhE1 and ΔadhE2 mutants were first evaluated under acidogenic conditions and compared to previously published data for the control strain . All the strains behaved the same, and no significant changes in the metabolic fluxes were recorded (Additional file 1: Fig. S3), except that butanol production was completely abolished in the ΔadhE2 mutant strain (Fig. 2, Additional file 1: Fig. S3).
The two mutant strains were then evaluated under solventogenic conditions and compared to previously published data for the control strain . The control and ΔadhE2 strains behaved the same, with no significant change in metabolic fluxes (Additional file 1: Fig. S3). However, the ΔadhE1 mutant exhibited a completely different behavior. In the first phase, before the “pseudo steady state” was reached, this mutant exhibited considerable fluctuations in growth, glucose consumption, and metabolite profiles. Under “pseudo steady state conditions,” the butanol and acetone fluxes were stable, while the butyrate flux showed fluctuations between 2.2 and 2.9 mmol g−1 h−1. In ΔadhE1, the butanol, ethanol, and acetone fluxes decreased by 60, 49, and 46 %, respectively (Additional file 1: Fig. S3), compared to the control strain; thus, the acetone and ethanol fluxes were not reduced as greatly as the butanol fluxes. These results support the previously proposed [1, 11, 12, 14] key role of AdhE1 in butanol production under solventogenic conditions and demonstrate that an adhE1 knockout strain with no polar effect on ctfAB transcription can still produce acetone. The level of ctfAB expression was 3-fold higher in the adhE1 knockout compared to the control strain. This indicates that the lower flux of acetone production is the result of a control at the enzyme level due to a lower acetoacetyl-CoA concentration and/or higher acetyl-CoA/butyryl-CoA concentrations. The remaining ability of the ΔadhE1 strain to produce butanol under solventogenesis is explained by the higher adhE2 expression (~127-fold higher than the control strain, but only 25 mRNA molecules/cell) (Table 1, Additional file 2: Dataset S1). For the ΔadhE1 mutant, the butyrate flux increased by 5-fold compared to the control strain (Additional file 1: Fig. S3), although neither ptb-buk (CA_C3076–CA_C3075) nor buk2 (CA_C1660) experienced a significant transcriptional increase (Additional file 2: Dataset S1). Thus, flux is controlled at the enzyme level via an increase in the butyryl-CoA pool due to the lower flux in the butanol pathway. However, as the AdhE2 level in the mutant is the same as the AdhE1 level in the control (6.31 × 104 versus 5.99 × 104 protein molecules/cell), the lower flux of butanol production can be explained by (i) a lower catalytic efficiency of AdhE2 for butyryl-CoA and/or NADH or (ii) a lower intracellular pH under solventogenic conditions that would be less optimal for AdhE2 that is normally expressed under alcohologenic conditions at neutral pH. The second hypothesis can be eliminated as the previously measured intracellular pH [4, 21] in solventogenic and alcohologenic cells are relatively close (5.5 and 5.95, respectively) as the ΔpH is inverted (more acidic inside) under alcohologenic conditions . Finally, as we will see below, the fact that ethanol flux is less affected than the butanol flux might be explained by the existence of an ethanol flux through the Pdc (pyruvate decarboxylase, encoded by CA_P0025) and bdhA/BdhB.
The two mutant strains were also evaluated under alcohologenic conditions and compared to previously published data for the control strain . The control and ΔadhE1 strains behaved the same, with no significant changes in metabolic fluxes (Additional file 1: Fig. S3). However, the ΔadhE2 mutant exhibited a completely different behavior; no flux toward butanol was detected, whereas fluxes toward butyrate became the primary fluxes, as opposed to butanol in the control strain (Additional file 1: Fig. S3). In addition, acetate levels increased by ~3-fold, and such changes were accompanied by changes in electron fluxes (Fig. 3), which are described in detail below. These phenomena were not observed by Cooksley et al.  with their adhE2 knockout mutant, as they performed batch fermentation without promoting alcohologenic conditions. As adhE1 was not expressed under the “alcohologenic conditions” of the ΔadhE2 mutant, the physiological function of adhE2 does not appear to be compensated by adhE1 (Table 1). To verify that loss of the butanol-producing ability under alcohologenesis did not result from loss of the pSOL1 megaplasmid [22, 23] but rather from the deletion of adhE2, the culture was switched to solventogenic conditions before the experiment was ended; under solventogenic conditions, high butanol and acetone production fluxes were recovered (data not shown).
The butanol pathway was analyzed for three different conditions in the respective mutants (Additional file 1: Fig. S2) by calculating the contribution of each of the five enzymes potentially involved in each of the two steps to the fluxes (see methods for the calculation).
Under acidogenesis, adhE1 was not expressed, and thus AdhE1 could not replace AdhE2 for the conversion of butyryl-CoA to butyraldehyde in the ΔadhE2 mutant (Additional file 1: Fig. S2). This failure of AdhE1 to replace AdhE2 led to the absence of butanol production in the ΔadhE1 mutant, which behaved the same as the control strain, leaving AdhE2 responsible for all the conversion. The ΔadhE1 mutant behaved the same as the control strain with respect to the conversion of butyraldehyde to butanol under these conditions, and AdhE2 (45 % of the flux), BdhB (34 % of the flux), and BdhA (14 % of the flux) were the main contributors (Additional file 1: Fig. S2). The ΔadhE2 mutant was not analyzed because it does not produce butanol.
Under solventogenesis, AdhE2 replaced AdhE1 for the conversion of butyryl-CoA to butyraldehyde in the ΔadhE1 mutant, while in the ΔadhE2 mutant, which behaved the same as the control strain, AdhE1 was responsible for all the conversion. The two main contributors to the conversion of butyraldehyde to butanol under these conditions were AdhE2 (67 % of the flux) and BdhB (30 % of the flux) in the ΔadhE1 mutant, while in the ΔadhE2 mutant, which behaved the same as the control strain, BdhB (75 % of the flux) and BdhA (16 % of the flux) were the main contributors (Additional file 1: Fig. S2).
Under alcohologenesis, adhE1 was not expressed (Table 1, Additional file 2: Dataset S1), and thus, AdhE1 could not replace AdhE2 for the conversion of butyryl-CoA to butyraldehyde in the ΔadhE2 mutant. This failure of AdhE1 to replace AdhE2 led to the absence of butanol production, while in the ΔadhE1 mutant, which behaved the same as the control strain, AdhE2 was responsible for all the conversion. The ΔadhE1 mutant behaved the same as the control strain with respect to the conversion of butyraldehyde to butanol under these conditions, and AdhE2 was the main contributor (Additional file 1: Fig. S2). The ΔadhE2 mutant was not analyzed because it does not produce butanol.
Two possible routes are known for the conversion of pyruvate to acetaldehyde in C. acetobutylicum: (i) a two-step reaction by pyruvate:ferredoxin oxidoreductase (PFOR) and acetaldehyde dehydrogenase via acetyl-CoA production or (ii) a one-step reaction by pyruvate decarboxylase (Pdc, encoded by CA_P0025) . In the wild-type strain, the former route is considered as the primary pathway [2, 25]. Under acidogenic and alcohologenic conditions of the ΔadhE2 mutant, ethanol production was observed, but no butanol production was detected (Fig. 2, Additional file 1: Fig. S3). As previously reported , AdhE1 retains only aldehyde dehydrogenase activity, whereas AdhE2 possesses both aldehyde and alcohol dehydrogenases activities. Thus, the ethanol production of the ΔadhE2 mutant suggests that the latter route is active. In other words, Pdc could be functional, and the ethanol dehydrogenase activity in acidogenesis could be due to BdhA, BdhB, or BdhC (Table 1). The same pathway might also be functional in solventogenesis and explains why in the ΔadhE1 mutant the ethanol flux was less affected than the butanol flux.
Because the predominant use of reduced ferredoxin is for hydrogen production , no significant effects were observed under acidogenesis in both the ΔadhE1 and ΔadhE2 mutants with respect to electron flux (Fig. 3). In addition, solventogenesis of the ΔadhE2 mutant exhibited similar flux levels to the control strain due to the small contribution of AdhE2 (5 % for butyraldehyde dehydrogenase function and 9 % for butanol dehydrogenase function) under these conditions in the control strain. However, under the same conditions as for ΔadhE1, both the fluxes for NADH, known as the partner of AdhE1 and AdhE2, and for NADPH, known as the partner of BdhA, BdhB, and BdhC, were reduced (by ~2.7-fold and 1.8-fold, respectively) due to decreased carbon fluxes toward alcohols (Fig. 3, Additional file 1: Fig. S3). The most striking changes were observed in the ΔadhE2 mutant under alcohologenesis, in which the primary use of reduced ferredoxin was switched from NADH to hydrogen production. The absence of butanol formation resulted in a ~3.6-fold decreased flux toward NADH production and a 1.7-fold increased flux toward hydrogen production (Fig. 3).
Common criteria used for quantitative transcriptomic analysis
To filter the data from only significant results, the same criteria used to compare the wild-type strain under different physiological conditions  were used to compare the mutants to the control strain. The first criterion was >4.0-fold higher expression or >4.0-fold lower expression in ΔadhE1 or ΔadhE2 than in the control strain under the same physiological condition, and the second criterion was >0.2 mRNA molecules per cell in at least one of the two strains being compared.
Genes affected by adhE1 or adhE2 deletion under acidogenesis
As alcohols are minor products under acidogenesis, the deletion of adhE1 or adhE2 did not significantly alter the metabolic flux map (Additional file 1: Fig. S3). However, a surprisingly large number of genes (100 genes increased in ΔadhE1, 108 genes decreased in ΔadhE1, 119 genes increased in ΔadhE2, 170 genes decreased in ΔadhE2) showed significant changes in mRNA molecules/cell in response to the deletion of each gene (Table 2). Furthermore, 50 genes (>4-fold increase) and 87 genes (>4-fold decrease) revealed the same patterns of change in both the ΔadhE1 and ΔadhE2 mutants (Table 2). The primary metabolism-related genes that influence metabolic fluxes did not exhibit significant changes, whereas mostly subordinate metabolism-related genes were affected (Additional file 1: Table S2, Additional file 1: S3, and Fig. 4).
Interestingly, a large portion (18 genes showed >a 4-fold increase, and 2 genes showed a >2.8-fold increase out of 30 genes proposed by Wang et al. ) of the cysteine metabolism regulator (CymR) regulon showed significantly increased expression in both mutants under acidogenesis (CymR regulons are indicated in Table 3). In particular, an operon involved in cysteine and sulfur metabolism (CA_C0102–CA_C0110) showed a >10-fold increase in both mutants. This operon was reported to respond to butyrate/butanol stresses and to be up-regulated under alcohologenesis in wild-type strains [12, 26, 27] and under solventogenesis in the Δptb mutant . In addition, the expression of two putative cysteine ABC transporter operons belonging to the CymR regulon [26, 27], namely CA_C0878–CA_C0880 and CA_C3325–CA_C3327), was also up-regulated.
A long gene cluster linked to iron/sulfur/molybdenum metabolism (CA_C1988–CA_C2019) exhibited significantly decreased expression (except for CA_C1988, CA_C1990, CA_C1992 and CA_C1995, for which some values were below the significance criterion of 4-fold but were higher than 3-fold) (Table 3, Additional file 2: Dataset S1). A part of this cluster, CA_C1988–CA_C1996, was previously reported to be down-regulated under oxygen-exposed conditions . Moreover, this cluster was shown by Schwarz et al.  to be repressed by butanol stress in an acidogenic chemostat.
Transcriptional changes due to adhE1 or adhE2 deletion under solventogenesis
Under solventogenesis, a drastic change in fluxes was observed in the ΔadhE1 mutant, while the fluxes remained unchanged in the ΔadhE2 mutant; additionally, as expected, more genes showed significant changes in ΔadhE1 than in ΔadhE2 (Table 2, Additional file 1: Table S4, Additional file 1: S5). Specifically, in ΔadhE1, 55 genes were up-regulated, and 127 genes were down-regulated (Table 2). In ΔadhE2, 22 genes were up-regulated, and 17 genes were down-regulated (Table 2). In contrast to the observations previously made under acidogenesis, no gene was significantly increased in both the ΔadhE1 and ΔadhE2 mutants, and only 1 gene (CA_C3612, encoding a hypothetical protein) was significantly decreased in both mutants.
In ΔadhE1, the CA_C0102–CA_C0110 operon which was shown to be up-regulated in acidogenesis and belongs to the CymR regulon, was also up-regulated by >18-fold under solventogenesis (Additional file 1: Table S4). However, the up-regulation of this operon (under alcohologenesis in the control strain, acidogenesis and solventogenesis in ΔadhE1, or acidogenesis in ΔadhE2) did not have striking shared features with the main product profile.
Interestingly, expression of the natAB operon (CA_C3551–CA_C3550) (>10-fold), encoding a potential Na+-ABC transporter, and the kdp gene cluster (CA_C3678–CA_C3682), encoding a potential K+ transporter (>20-fold), was highly up-regulated under solventogenesis (Additional file 1: Table S4, Additional file 2: Dataset S1) in ΔadhE1. The natAB operon and the kdp gene cluster have previously been reported to be up-regulated by both acetate and butyrate stress . As the ability of the ΔadhE1 mutant to produce butanol was highly affected and as butyrate and acetate were the primary fermentation products (Fig. 2), this strain struggled to survive under acidic conditions (i.e., under the pH of 4.4 for solventogenesis); consequently, genes involved in ion transport were up-regulated.
The operon CA_P0029–CA_P0030, which potentially encodes a transporter and an isochorismatase, was up-regulated under acidogenesis in both mutants as well as under solventogenesis in ΔadhE2 (>20-fold) (Table 2, Additional file 1: Table S5). Two neighboring genes, CA_C3604 (ilvD), encoding dihydroxyacid dehydratase linked to valine/leucine/isoleucine biosynthesis, and CA_C3605 (gntP), encoding high affinity gluconate/L-idonate permease, exhibited striking increases (>120-fold) (Additional file 1: Table S5) in ΔadhE2.
As described above, the solventogenic culture of ΔadhE1 has a lower glucose consumption rate than the control strain (Fig. 2) and consequently more glucose remained unconsumed in the medium. Accordingly, numerous genes related to sugar metabolism were down-regulated under this metabolic state. For instance, all the structural genes on the mannitol phosphotransferase system (PTS)-related operon mtlARFD (CA_C0154–CA_C0157) and the mannose PTS-related operon (CA_P0066–CA_P0068) were decreased by >10-fold (Additional file 1: Table S4).
Interestingly, one of two operons encoding a quorum-sensing system and putatively involved in sporulation, CA_C0078–CA_C0079 (agrBD) , was strongly down-regulated (infinity-fold for CA_C0078 and 667-fold for CA_C0078) in ΔadhE2 relative to the control strain (Additional file 1: Table S5). However, the other operon, CA_C0080–CA_C0081 (agrCA), did not significantly change (<3-fold decreases) (Additional file 2: Dataset S1). Quantitatively, less than 1 agrCA mRNA molecule was found per cell, whereas more than 1 agrBD mRNA molecule was found per cell under all conditions in the control strain . These different expression levels are not surprising because agrBD and agrCA are independently transcribed [31–33]. In addition, agrBD was repressed under all conditions in ΔadhE2, although the sporulation of this mutant was not affected (Additional file 2: Dataset S1).
Transcriptional changes due to adhE1 or adhE2 deletion under alcohologenesis
Under alcohologenesis, a drastic change in fluxes was observed in the ΔadhE2 mutant, while in the ΔadhE1 mutant, the fluxes remained unchanged. As expected, more genes showed significant changes in the ΔadhE2 mutant than in the ΔadhE1 mutant (Table 2). Specifically, in ΔadhE1, only 1 gene was up-regulated (agrB), and 14 genes were down-regulated, while in ΔadhE2, 35 genes were up-regulated, and 38 genes were down-regulated.
The most dynamic changes in the ΔadhE2 mutant were observed in CA_C3604 (ilvD, 297-fold) and CA_C3605 (gntP, 301-fold) (Additional file 1: Table S7). As mentioned previously, these genes were highly up-regulated (>84-fold) under all the conditions in the ΔadhE2 mutant (Additional file 2: Dataset S1). Interestingly, two genes located immediately downstream of adhE2, CA_P0036, which encodes a cytosolic protein of unknown function, and CA_P0037, which encodes a potential transcriptional regulator, exhibited a ~ 9-fold increase under alcohologenesis (Additional file 1: Table S7) in ΔadhE2.
A sucrose metabolism operon comprising scrAKB (CA_C0423–CA_C0425), encoding a PTS IIBCA domain on a single gene, fructokinase and sucrose-6-P hydrolase [35, 36], was strikingly down-regulated (>47-fold) (Additional file 1: Table S6). Moreover, the gene immediately upstream, scrT (CA_C0422) (encoding a putative transcriptional antiterminator), and the gene downstream, CA_C0426, encoding a putative AraC-type of regulator, were also decreased, by 9.3-fold and 8-fold, respectively (Additional file 1: Table S6). The similar expression patterns of CA_C0422, CA_C0426, and scrAKB support the hypotheses of previous studies regarding their roles in regulating scrAKB [35, 36].
As expected based on the reduced consumption of glycerol (approximately one-fourth of the control strain) (Fig. 2) in ΔadhE2, the gene cluster for glycerol transport and utilization (CA_C1319-CA_C1322) was down-regulated (>4.3-fold) under these conditions (Additional file 1: Table S7).
Most arginine biosynthesis-related genes known to respond negatively to butanol and butyrate stress  (i.e., CA_C0316 (argF/I), CA_C0973–CA_C0974 (argGH), CA_C2389–CA_C2388 (argBD), CA_C2390–CA_C2391 (argCJ), CA_C2644 (carB), and CA_C2645 (carA)) were significantly down-regulated (>4-fold decrease) (Additional file 1: Table S7) in ΔadhE2. As “alcohologenic cultures” of ΔadhE2 produced 70 mM of butyrate and no butanol (Fig. 2), this down-regulation is consistent with the high butyrate stress (50 mM) response .
CA_C3486, which encodes a multimeric flavodoxin, was decreased by 4.4-fold in ΔadhE2 (Additional file 1: Table S7), resulting in a loss of butanol production under alcohologenesis. This finding is consistent with the proposed hypothesis  that under alcohologenesis, the gene product of CA_C3486 may function as a redox partner between the hydrogenase and ferredoxin-NAD+ reductase and may participate in the redistribution of electron fluxes in favor of butanol formation.
The results presented here support the hypothesis of the roles of AdhE1 and AdhE2 in butanol formation, namely that AdhE1 is the key enzyme for butanol formation in solventogenesis and that AdhE2 is the key enzyme for butanol formation in alcohologenesis. Furthermore, this study also demonstrates the metabolic flexibility of C. acetobutylicum in response to genetic alteration of its primary metabolism.
Bacterial strains and plasmid construction
All C. acetobutylicum strains used in this study and in the control study were constructed from the C. acetobutylicum ATCC 824 ΔCA_C1502 Δupp mutant strain, which was constructed for rapid gene knockout and gene knockin . Detailed procedures, including all strains and primers used, are described in the online supporting information (Supplementary experimental procedures).
All batch cultures were performed under strict anaerobic conditions in synthetic medium (MS), as previously described . C. acetobutylicum was stored in spore form at −20 °C after sporulation in MS medium. Heat shock was performed for spore germination by immersing the 30- or 60-mL bottle into a water bath at 80 °C for 15 min.
All the phosphate-limited continuous cultivations were performed as previously described by Vasconcelos et al.  and Girbal et al.  like in the control strain study . The chemostat was fed a constant total of 995 mM of carbon and maintained at a dilution rate of 0.05 h−1. The maintained pH of the bioreactor and the supplied carbon sources of each metabolic state were as follows: for acidogenesis, pH 6.3, with 995 mM of carbon from glucose; for solventogenesis, pH 4.4, with 995 mM of carbon from glucose; and for alcohologenesis, pH 6.3, with 498 mM of carbon from glucose and 498 mM of carbon from glycerol.
RNA extraction and microarray
Total RNA isolation and microarray experiments were performed as previously described . Briefly, 3 mL of chemostat cultures was sampled, immediately frozen in liquid nitrogen and ground with 2-mercaptoethanol. RNA was extracted by using an RNeasy Midi kit (Qiagen, Courtaboeuf, France) and RNase-Free DNase Set (Qiagen) per the manufacturer’s protocol. The RNA quantity and integrity were monitored using an Agilent 2100 Bioanalyzer (Agilent Technologies, Massy, France) and a NanoDrop ND-1000 spectrophotometer (Labtech France, Paris, France) at 260 and 280 nm. All microarray procedures were performed per the manufacturer’s protocol (Agilent One-Color Microarray-Based Exon Analysis).
The optical density at 620 nm (OD620 nm) was monitored and used to calculate the biomass concentration with the correlation factor between dry cell weight and OD620 nm (path length 1 cm) of 0.28, which was experimentally determined from continuous cultures and was used in a control strain study . The glucose, glycerol, acetate, butyrate, lactate, pyruvate, acetoin, acetone, ethanol, and butanol concentrations were determined using high-performance liquid chromatography (HPLC), as described by Dusséaux et al. . The concentration of the eluent H2SO4 was changed to 0.5 mM, as this concentration was required to optimize the mobile phase for the control strain study .
Calculation of the cytosolic proteins concentration (protein molecules per cell)
In a previously published work , we quantified the amount of (i) mRNA molecules per cell for all genes and (ii) protein molecules per cell (for approximately 700 cytosolic proteins) for steady-state chemostat cultures (at a specific growth rate of 0.05 h−1) of C. acetobutylicum under different physiological conditions. For 96 % of the cytosolic proteins that could be quantified, a linear relationship was obtained, with an R2 > 0.9, when the numbers of protein molecules per cell were plotted against the numbers of mRNA molecules per cell. This result indicated that for steady-state continuous cultures run at the same specific growth rate and with the same total amount of carbon supplied, the rate of protein turnover is proportional to the mRNA content for 96 % of the genes. As the mutants were cultivated in chemostat culture at the same growth rate (0.05 h−1), we used the absolute protein synthesis rates previously calculated  for each of the 700 genes to calculate the amount of protein molecule per cell for each of these 700 genes in the different mutants. (Additional file 2: Dataset S1).
Calculation of the contribution of different enzymes on the butanol flux
The contribution of the 5 proteins potentially involved in the butanol pathway, namely AdhE1, AdhE2, BdhA, BdhB, and BdhC, was made as previously described  by assuming that all five enzymes function at their Vmax and using the calculated amount of each protein per cell (Additional file 2: Dataset S1).
flippase recognition target
Nair RV, Bennett GN, Papoutsakis ET. Molecular characterization of an aldehyde/alcohol dehydrogenase gene from Clostridium acetobutylicum ATCC 824. J Bacteriol. 1994;176(3):871–85.
Vasconcelos I, Girbal L, Soucaille P. Regulation of carbon and electron flow in Clostridium acetobutylicum grown in chemostat culture at neutral pH on mixtures of glucose and glycerol. J Bacteriol. 1994;176(5):1443–50.
Girbal L, Vasconcelos I, Saint-Amans S, Soucaille P. How neutral red modified carbon and electron flow in Clostridium acetobutylicum grown in chemostat culture at neutral pH. FEMS Microbiol Rev. 1995;16(2):151–62.
Girbal L, Soucaille P. Regulation of Clostridium acetobutylicum metabolism as revealed by mixed-substrate steady-state continuous cultures: role of NADH/NAD ratio and ATP pool. J Bacteriol. 1994;176(21):6433–8.
Bahl H, Andersch W, Gottschalk G. Continuous production of acetone and butanol by Clostridium acetobutylicum in a two-stage phosphate limited chemostat. Eur J Appl Microbiol Biotechnol. 1982;15(4):201–5.
Fischer RJ, Helms J, Durre P. Cloning, sequencing, and molecular analysis of the sol operon of Clostridium acetobutylicum, a chromosomal locus involved in solventogenesis. J Bacteriol. 1993;175(21):6959–69.
Yoo M, Bestel-Corre G, Croux C, Riviere A, Meynial-Salles I, Soucaille P. A Quantitative System-Scale Characterization of the Metabolism of Clostridium acetobutylicum. MBio 2015, 6(6):e01808–01815.
Sauer U, Dürre P. Differential induction of genes related to solvent formation during the shift from acidogenesis to solventogenesis in continuous culture of Clostridium acetobutylicum. FEMS Microbiol Lett. 1995;125(1):115–20.
Fontaine L, Meynial-Salles I, Girbal L, Yang X, Croux C, Soucaille P. Molecular characterization and transcriptional analysis of adhE2, the gene encoding the NADH-dependent aldehyde/alcohol dehydrogenase responsible for butanol production in alcohologenic cultures of Clostridium acetobutylicum ATCC 824. J Bacteriol. 2002;184(3):821–30.
Leang C, Ueki T, Nevin KP, Lovley DR. A genetic system for Clostridium ljungdahlii: a chassis for autotrophic production of biocommodities and a model homoacetogen. Appl Environ Microbiol. 2013;79(4):1102–9.
Waghmare SK, Caputo V, Radovic S, Bruschi CV. Specific targeted integration of kanamycin resistance-associated nonselectable DNA in the genome of the yeast Saccharomyces cerevisiae. Biotechniques 2003, 34(5):1024–1028, 1033.
Storici F, Bruschi CV. Involvement of the inverted repeat of the yeast 2-micron plasmid in Flp site-specific and RAD52-dependent homologous recombination. Mol Gen Genet. 2000;263(1):81–9.
Cornillot E, Nair RV, Papoutsakis ET, Soucaille P. The genes for butanol and acetone formation in Clostridium acetobutylicum ATCC 824 reside on a large plasmid whose loss leads to degeneration of the strain. J Bacteriol. 1997;179(17):5442–7.
Wang Q, Venkataramanan KP, Huang H, Papoutsakis ET, Wu CH. Transcription factors and genetic circuits orchestrating the complex, multilayered response of Clostridium acetobutylicum to butanol and butyrate stress. BMC Syst Biol. 2013;7:120.
Alsaker KV, Paredes C, Papoutsakis ET. Metabolite stress and tolerance in the production of biofuels and chemicals: gene-expression-based systems analysis of butanol, butyrate, and acetate stresses in the anaerobe Clostridium acetobutylicum. Biotechnol Bioeng. 2010;105(6):1131–47.
Honicke D, Lutke-Eversloh T, Liu Z, Lehmann D, Liebl W, Ehrenreich A. Chemostat cultivation and transcriptional analyses of Clostridium acetobutylicum mutants with defects in the acid and acetone biosynthetic pathways. Appl Microbiol Biotechnol. 2014;98(23):9777–94.
Schwarz KM, Kuit W, Grimmler C, Ehrenreich A, Kengen SW. A transcriptional study of acidogenic chemostat cells of Clostridium acetobutylicum–Cellular behavior in adaptation to n-butanol. J Biotechnol. 2012;161(3):366–77.
Liu D, Chen Y, Ding F, Guo T, Xie J, Zhuang W, Niu H, Shi X, Zhu C, Ying H. Simultaneous production of butanol and acetoin by metabolically engineered Clostridium acetobutylicum. Metab Eng. 2015;27:107–14.
Fischer RJ, Oehmcke S, Meyer U, Mix M, Schwarz K, Fiedler T, Bahl H. Transcription of the pst operon of Clostridium acetobutylicum is dependent on phosphate concentration and pH. J Bacteriol. 2006;188(15):5469–78.
Croux C, Nguyen NPT, Lee J, Raynaud C, Saint-Prix F, Gonzalez-Pajuelo M, Meynial-Salles I, Soucaille P. Construction of a restriction-less, marker-less mutant useful for functional genomic and metabolic engineering of the biofuel producer Clostridium acetobutylicum. Biotechnol Biofuel. 2016;9:21.
Dusseaux S, Croux C, Soucaille P, Meynial-Salles I. Metabolic engineering of Clostridium acetobutylicum ATCC 824 for the high-yield production of a biofuel composed of an isopropanol/butanol/ethanol mixture. Metab Eng. 2013;18:1–8.
CC, IMS, and PS conceived the study; MY performed all the experimental work. MY and PS performed the data analysis and drafted the manuscript. PS supervised the work. All authors read and approved the final manuscript.
We thank Sophie Lamarre and Lidwine Trouilh for help with the data analysis.
Availability of supporting data
Microarray data can be accessed at GEO through accession number GSE69973.
The authors declare that they have no competing interests.
Consent for publication
Ethical Approval and Consent to participate
This work was financially supported by the European Community’s Seventh Framework Program “CLOSTNET” (PEOPLE-ITN-2008-237,942) to Minyeong Yoo.
Authors and Affiliations
INSA, UPS, INP, LISBP, Université de Toulouse, Toulouse, France
Minyeong Yoo, Christian Croux, Isabelle Meynial-Salles & Philippe Soucaille
INRA, UMR792, Toulouse, France
Minyeong Yoo, Christian Croux, Isabelle Meynial-Salles & Philippe Soucaille
CNRS, UMR5504, Toulouse, France
Minyeong Yoo, Christian Croux, Isabelle Meynial-Salles & Philippe Soucaille
Metabolic Explorer, Biopôle Clermont-Limagne, Saint Beauzire, France
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Yoo, M., Croux, C., Meynial-Salles, I. et al. Elucidation of the roles of adhE1 and adhE2 in the primary metabolism of Clostridium acetobutylicum by combining in-frame gene deletion and a quantitative system-scale approach.
Biotechnol Biofuels9, 92 (2016). https://doi.org/10.1186/s13068-016-0507-0