Skip to main content

n-Butanol production in Saccharomyces cerevisiae is limited by the availability of coenzyme A and cytosolic acetyl-CoA

Abstract

Background

Butanol isomers are regarded as more suitable fuel substitutes than bioethanol. n-Butanol is naturally produced by some Clostridia species, but due to inherent problems with clostridial fermentations, industrially more relevant organisms have been genetically engineered for n-butanol production. Although the yeast Saccharomyces cerevisiae holds significant advantages in terms of scalable industrial fermentation, n-butanol yields and titers obtained so far are only low.

Results

Here we report a thorough analysis and significant improvements of n-butanol production from glucose with yeast via the acetoacetyl-CoA-derived pathway. First, we established an improved n-butanol pathway by testing various isoenzymes of different pathway reactions. This resulted in n-butanol titers around 15 mg/L in synthetic medium after 74 h. As the initial substrate of the n-butanol pathway is acetyl-coenzyme A (acetyl-CoA) and most intermediates are bound to coenzyme A (CoA), we increased CoA synthesis by overexpression of the pantothenate kinase coaA gene from Escherichia coli. Supplementation with pantothenate increased n-butanol production up to 34 mg/L. Additional reduction of ethanol formation by deletion of alcohol dehydrogenase genes ADH1-5 led to n-butanol titers of 71 mg/L. Further expression of a mutant form of an ATP independent acetylating acetaldehyde dehydrogenase, adhEA267T/E568K, converting acetaldehyde into acetyl-CoA, resulted in 95 mg/L n-butanol. In the final strain, the n-butanol pathway genes, coaA and adhE A267T/E568K, were stably integrated into the yeast genome, thereby deleting another alcohol dehydrogenase gene, ADH6, and GPD2-encoding glycerol-3-phosphate dehydrogenase. This led to a further decrease in ethanol and glycerol by-product formation and elevated redox power in the form of NADH. With the addition of pantothenate, this strain produced n-butanol up to a titer of 130 ± 20 mg/L and a yield of 0.012 g/g glucose. These are the highest values reported so far for S. cerevisiae in synthetic medium via an acetoacetyl-CoA-derived n-butanol pathway.

Conclusions

By gradually increasing substrate supply and redox power in the form of CoA, acetyl-CoA, and NADH, and decreasing ethanol and glycerol formation, we could stepwise increase n-butanol production in S. cerevisiae. However, still further bottlenecks in the n-butanol pathway must be deciphered and improved for industrially relevant n-butanol production levels.

Background

Biofuels produced by microbial fermentations of renewable resources represent an important replacement for fossil fuels. Today, bioethanol is the most important renewable gasoline substitute. However, higher alcohols like butanol are regarded as more suitable fuel substitutes, due to their higher energy density, lower hygroscopicity, and therefore lower corrosiveness [1].

n-Butanol is produced by various Clostridium species via the acetone butanol ethanol (ABE) fermentation process with a ratio of 1:3:6 and titers up to 13 g/l [2, 3]. However, clostridial fermentations are associated with several problems, including sporulation, slow growth rates, bacteriophage infections, and strict anaerobic cultivations [4]. For this reason, other industrially more relevant organisms have been engineered for n-butanol production by introducing variants of the clostridial pathway [5] or using ketoacid degradation pathways [6]. Although n-butanol titers up to 15 g/l [7] were achieved in E. coli, large-scale industrial fermentations with this bacterium also suffer from the high risk of phage infections and bacterial contaminations, for instance [8, 9]. Compared to E. coli, S. cerevisiae holds significant advantages in terms of scalable industrial fermentations, due to long-lasting experiences with this yeast in fermentation processes, high robustness, and tolerance against inhibitory compounds [10].

S. cerevisiae is able to produce small amounts of n-butanol on rich medium via endogenous pathways dependent on threonine or glycine catabolism, the first one stimulated by deletion of ADH1 [11, 12]. By improving these pathways, n-butanol titers up to 242.8 and 92 mg/L, respectively, could be achieved. Combining the amino acid-dependent endogenous pathways with a synthetic ABE pathway and a mutant version of translation initiation factor eIF2B in an adh1 mutant strain resulted in up to 300 mg/L n-butanol, but only in YEPD medium after 15–20 days [13]. Nevertheless, as biofuel production from amino acids is economically not viable, for industrial use threonine or glycine overproduction strains would be required which, however, exhibit lower theoretical maximum yields [1].

Therefore, it is more promising to enhance n-butanol production by improvements of variants of the n-butanol pathway expressed in S. cerevisiae. In this pathway, two molecules of acetyl-CoA are condensed to acetoacetyl-CoA, which is reduced to 3-hydroxybutyryl-CoA and then dehydrated to crotonyl-CoA. This is further reduced to butyryl-CoA, butyraldehyde, and finally n-butanol (Fig. 1). Extending the pathway at the level of butyryl-CoA in a reversal of the β-oxidation cycle would enable even the production of longer molecules [14, 15]. Besides alcohols, also the intermediates of the pathway are valuable compounds or can be used to produce the corresponding aldehydes, fatty acids, and alkanes with different chain lengths. First attempts to establish n-butanol production in yeast using this acetoacetyl-CoA-derived pathway resulted in n-butanol titer as low as only 2.5 mg/L [16]. Further pathway optimization and engineering of substrate supply could enhance n-butanol titers up to 120 mg/L [9, 17]. However, these titers were obtained only in synthetic complete media in high-cell density fermentations and were critically reviewed in [1].

Fig. 1
figure 1

Metabolic pathway for n-butanol production via reverse β-oxidation in yeast. n-Butanol production and reverse β-oxidation are shown as a process of one or several rotations, respectively, of an acetoacetyl-CoA-derived synthesis pathway. This pathway can generate various compounds with different chain lengths. The different isoenzymes used in this study (Fig. 2, Table 2) are shown for each reaction. The reaction of Ss nphT7 (use of malonyl-CoA and acetyl-CoA for production of acetoacetyl-CoA) is shown in light blue

As n-butanol yield in yeast is still much lower than in E. coli, it was hypothesized that cytosolic acetyl-CoA as the initial substrate of the n-butanol pathway is a limiting factor [9]. In yeast, acetyl-CoA is present in four different compartments: cytosol, mitochondria, peroxisomes, and the nucleus. Cytosolic acetyl-CoA is produced via the decarboxylation of pyruvate into acetaldehyde, which is oxidized to acetate by aldehyde dehydrogenases and finally converted to acetyl-CoA, at the expense of two ATPs. This pathway is known as pyruvate dehydrogenase (Pdh) bypass. Aldehyde dehydrogenases compete with alcohol dehydrogenases (Adh) which due to the crabtree effect convert most of the sugar-derived acetaldehyde into ethanol [18]. On the other hand, mitochondrial acetyl-CoA obtained via pyruvate dehydrogenase obviously cannot be transported to the cytosol [19]. Different strategies were tested to increase the level of cytosolic acetyl-CoA. In one approach, alcohol dehydrogenases were deleted, the reactions of the Pdh bypass were overexpressed, and the withdrawal of acetyl-CoA into the glyoxylate cycle was prevented [18]. In another approach, an ATP independent pyruvate dehydrogenase complex from Enterococcus faecalis was expressed in the yeast cytosol replacing the Pdh bypass, but being dependent on supplementation with lipoic acid [20]. Another possibility of ATP independent cytosolic acetyl-CoA production is the use of acetylating acetaldehyde dehydrogenases or pyruvate formate lyases [21, 22].

Since the production of cytosolic acetyl-CoA proceeds from two substrates, acetate and coenzyme A (CoA), the level of CoA might limit acetyl-CoA production as well. Moreover, nearly all the intermediates of the n-butanol pathway exist as CoA esters which might additionally decrease CoA availability. Therefore, we hypothesized that engineering of CoA biosynthesis could further increase the concentration of acetyl-CoA in yeast, resulting in elevated n-butanol production. CoA biosynthesis starts from pantothenate. In yeast, pantothenate can be produced via an endogenous pathway starting from 2-ketoisovalerate, derived from the valine metabolic pathway, and spermine, derived from l-ornithine and methionine [23]. Besides, it can be taken up from the medium by the transporter Fen2 [24]. Pantothenate is then converted into 4-phosphopantothenate by pantothenate kinase Cab1, which is postulated to catalyze the rate-limiting step in CoA biosynthesis. Cysteine is incorporated by phosphopantothenate cysteine ligase Cab2, which is weakly glucose repressed like Cab1 [23]. The resulting product is finally converted into CoA in three further steps [25].

In the present study, we first established and improved n-butanol production in S. cerevisiae by testing different variants of the acetoacetyl-CoA-derived pathway. We then improved synthesis of CoA by overexpression of pantothenate kinase coaA from E. coli, which is not inhibited by acetyl-CoA [23], or by overexpression of the yeast Fen2 pantothenate transporter. The level of CoA could further be increased by adding pantothenate to the medium. By deletion of alcohol dehydrogenases and glycerol-3-phosphate dehydrogenase genes, we increased the availability of acetaldehyde and NADH as driving forces for the n-butanol pathway. Overexpression of a mutant form of an acetylating acetaldehyde dehydrogenase could further improve n-butanol production, resulting in n-butanol titers and yields of up to 130 Â± 20 mg/L and 0.012 g/g glucose, respectively. These are the highest values ever reported for S. cerevisiae in synthetic medium via the acetoacetyl-CoA-derived pathway.

Results and discussion

Comparison of different variants of an acetoacetyl-CoA-derived n-butanol pathway in yeast

In order to engineer S. cerevisiae for n-butanol production, we introduced different variants of an acetoacetyl-CoA-derived n-butanol pathway in CEN.PK113-5D (Figs. 1, 2, Additional file 1: Figure S1). In our first experiments, we chose multicopy vectors for high-level expression of heterologous genes but soon recognized that the use of centromeric vectors generally resulted in significantly higher n-butanol titers (data not shown) and also higher enzyme activities of most of the pathway enzymes in crude extracts (Tables 1 and 2). Indeed, such problems by using multicopy vectors have also been reported by others [9, 20, 26]. The problems might be related to instabilities, vector burden, or transcriptional obstacles due to the high copy numbers and the strong promoters.

Fig. 2
figure 2

n-Butanol production with different variants of acetoacetyl-CoA-derived n-butanol pathways in S. cerevisiae strain CEN.PK113-5D. a The n-butanol concentrations after 74 h of semi-anaerobic fermentations in SMD medium are shown. Error bars represent the standard deviation of three independent replicates. b The respective enzyme combinations of n-butanol pathway vectors are displayed. Genes from Saccharomyces cerevisiae (Sc), Clostridium acetobutylicum (Ca), E.coli (Ec), Euglena gracilis (Eg), Streptomyces collinus (Sco), Streptomyces sp. CL190 (Ss), and Treponema denticola (Td) are indicated by prefixes in superscript

Table 1 Yeast strains and plasmids used in this study
Table 2 In vitro activities of various enzymes involved in n-butanol production pathways expressed in S. cerevisiae

On the other hand, even without a heterologous pathway, adh1 mutants of S. cerevisiae are able to produce n-butanol in YEPD medium via an endogenous 2-keto acid-derived pathway from threonine [11]. As we observed comparable results in YEPD medium with Δadh1-5 (VSY0) and Δadh1,3,5 (JDY4) mutant strains (n-butanol production from 50 mg/L up to 120 mg/L; data not shown), all further fermentations were carried out in synthetic minimal medium with 20 g/L glucose (SMD), where no n-butanol production was observed for VSY0 (Additional file 1: Figure S2). For most reactions of the acetoacetyl-CoA-derived pathway, we tested several isoenzymes from different organisms. We performed fermentation experiments (Fig. 2, Additional file 1: Figure S1) and also tested enzyme activities in crude extracts (Table 2). We tested the different isoforms in various combinations, always combined on a single centromeric plasmid (Table 1). Several strong promoters and terminators were used for expression (Additional file 2: Table S1). Combinations were chosen due to the measured enzyme activities, and enzymes were also changed if they were assumed to catalyze critical steps.

All fermentations were started with an OD600 of 0.3 and were performed for at least 100 h under semi-anaerobic conditions at 30 Â°C. The highest n-butanol titers were reached at about 74 h with close to 20 mg/L (Additional file 1: Figure S1). The gene combinations on plasmids pVS1, 6, 9, 10, and 11 turned out to perform best (Fig. 2). The bad performance of pVS7 and pVS8 seems to be due to the low activity of butanol dehydrogenase encoded by bdhB [27]. For the thiolase reaction, the yeast Erg10 enzyme and the C. acetobutylicum thlA resulted in comparable n-butanol production although the in vitro activity of the bacterial enzyme seems to be considerably higher. However, the enzyme assay is performed in the reverse reaction by conversion of acetoacetyl-CoA into two acetyl-CoA. NphT7 as a malonyl-CoA-dependent acyltransferase, which catalyzes an irreversible reaction [28], only marginally improved n-butanol production but works on the expense of one additional ATP for malonyl-CoA synthesis. For the conversion of crotonyl-CoA to butyryl-CoA, the trans-2-enoyl-CoA reductase ter enzymes [9] seem to be superior to crotonyl-CoA reductase ccr [29] (compare pVS10 and pVS9). Although the butyraldehyde dehydrogenase eutE [15] did not improve n-butanol production, assuming that the different ter enzymes did not have a great influence, we decided to include it in our further experiments due to its high enzyme activity and more reliable n-butanol production (see low standard deviation). Therefore, we continued with the enzyme combinations of pVS6, because of acetyl-CoA-dependent yeast’s own thiolase (Erg10) reaction, which is superior to ATP- and malonyl-CoA-consuming NphT7 of pVS9.

Overexpression of pantothenate kinase coaA and addition of pantothenate increase n-butanol production

Although we could demonstrate reasonable n-butanol production from glucose with the heterologous acetoacetyl-CoA-derived pathways, the yields and titers were still low. As the substrate acetyl-CoA and also most intermediates of the pathway are bound to coenzyme A (CoA), we speculated that a limiting availability of CoA might be a reason for the low n-butanol production. Therefore, we aimed to increase the CoA pool by optimization of its anabolic synthesis route. In yeast, CoA can either be produced via the endogenous pantothenate synthesis pathway [23] or from exogenously supplied pantothenate. We overexpressed pantothenate kinase as this step has been reported to limit CoA synthesis [25]. We chose the coaA enzyme from E. coli as this is not inhibited by acetyl-CoA [30]. On the other hand, we overexpressed the yeast pantothenate transporter Fen2. CoaA and FEN2 were expressed under control of the methionine-repressible MET25 promoter [31]. The MET25 promoter is most active in the absence of methionine and is gradually repressed up to 1 mM methionine, whereas 2 mM as a final concentration was chosen to guarantee repression. Therefore, fermentations were performed in SMD medium with strain CEN.PK113-5D transformed with pVS6 (n-butanol pathway) and either pVS5_4 (coaA) or pVS5_6 (FEN2), in the presence of different amounts of methionine. With decreasing methionine concentrations and therefore increasing expression of coaA or FEN2, in both cases the production of n-butanol increased up to 19 and 16 mg/L, respectively (Fig. 3, Additional file 1: Figure S3). Actually, the SMD medium contains 0.84 ÂµM pantothenate originating from the yeast nitrogen base (BD Difco YNB). Surprisingly, additional overexpression of phosphopantothenoylcysteine synthetase Cab2 or co-overexpression of coaA and Fen2 led to a decrease in n-butanol production (data not shown). For unknown reasons, it was not possible to get yeast transformants which overexpress all three genes together, even when expressed behind other promoters. Nevertheless, the results indicate that the availability of CoA is indeed a limiting factor for n-butanol production. However, it might be that also other CoA-derived products will increase. But, as we overexpress acetylating acetaldehyde dehydrogenase, which needs CoA for acetyl-CoA synthesis, most of the CoA will increase acetyl-CoA production. Furthermore, a lot of CoA is bound to the intermediates of n-butanol pathway and is not available for other reactions.

Fig. 3
figure 3

n-Butanol production with n-butanol pathway (pVS6) and coaA overexpression under control of the MET25 promoter (pVS5_4) in CEN.PK113-5D with addition of methionine. Five different concentrations of methionine were added to SMD medium: 0 mM (diamond), 0.125 mM (cross), 0.25 mM (circle), 0.5 mM (triangle), and 2 mM (square). Error bars represent the standard deviation of three independent replicates

To investigate whether higher concentrations of pantothenate can further increase n-butanol production, fermentations were performed with the coaA-overexpressing strain CEN.PK113-5D (pVS6, pVS5_4) after addition of 0–50 ÂµM pantothenate to SMD medium (without methionine). n-Butanol production further increased from 15 mg/L (0 ÂµM pantothenate) to 34 mg/L at a pantothenate concentration of 25 ÂµM but decreased slightly at 50 ÂµM pantothenate (Fig. 4). Nevertheless, for industrial applications pantothenate addition would be too expensive and not acceptable. Therefore, our results indicate that for industrial settings it might be promising to increase flux through the endogenous pantothenate synthesis pathway together with coaA overexpression.

Fig. 4
figure 4

n-Butanol production with n-butanol pathway (pVS6) and coaA overexpression (pVS5_4) in CEN.PK113-5D with addition of pantothenate. Five different concentrations of pantothenate were added to SMD medium: 0 ÂµM (square), 0.5 ÂµM (triangle), 5 ÂµM (circle), 25 ÂµM (cross), and 50 ÂµM (diamond). Error bars represent the standard deviation of three independent replicates

For all important fermentations in this study, n-butanol titers (mg/L) are compared in Fig. 5, n-butanol yields (g/g consumed glucose) in Additional file 1: Figure S4, glucose consumption/residual glucose in Additional file 1: Figure S5, ethanol and glycerol production in Additional file 1: Figure S6, and OD600 values in Additional file 1: Figure S7.

Fig. 5
figure 5

Comparison of n-butanol production with different strains. n-Butanol concentrations after 74 h of semi-anaerobic fermentations in SMD medium are shown. Compared are the fermentations of the wildtype CEN.PK113-5D, VSY0 (Δadh1-5), and VSY8 (Δadh1-5 sfa1 with adhE A267T/E568K) with plasmid pVS6 (n-butanol pathway), with or without coaA overexpression (pVS5_4), and strain VSY10 (Δadh1-6 sfa1 gpd2, with n-butanol pathway genes of pVS6, coaA, and adhE A267T/E568K), in the absence or presence of additional 25 ÂµM pantothenate, in SMD medium. Error bars represent the standard deviation of three independent replicates and a statistical analysis is shown in Additional file 2: Table S4

Deletion of alcohol dehydrogenases creates a push to n-butanol production

The next step was to redirect the metabolic flux away from the predominant fermentation product ethanol to n-butanol, and to increase the redox power in the form of NADH. Therefore, alcohol dehydrogenase genes ADH1-5 were deleted in order to reduce NADH-dependent acetaldehyde reduction and the concomitant formation of ethanol. As expected, we observed less ethanol production, but therefore an increased glycerol formation as a result of an enhanced alternative NADH oxidation route of yeast (Additional file 1: Figure S6). Fermentations of VSY0 (Δadh1-5) transformed with pVS6 (n-butanol pathway) resulted in n-butanol production up to 30 mg/L (Additional file 1: Figure S2), which is a significant increase compared to the reference strain (13 mg/L) (Fig. 5). The effect of elevated CoA supply was also tested in this strain by introducing pVS5_4 (coaA) together with pVS6. n-Butanol production was significantly higher in the adh1-5 deletion strain with coaA overexpression (52 versus 30 mg/L) (Additional file 1: Figure S2, Fig. 5). After the addition of 25 ÂµM pantothenate, the final n-butanol titer even increased up to 71 mg/L in the Δadh1-5 (pVS6, pVS5_4) strain (Additional file 1: Figure S8, Fig. 5). These results show that providing more substrate in the form of acetaldehyde and more reducing power in the form of NADH can increase n-butanol production.

Improving n-butanol production with an acetylating acetaldehyde dehydrogenase

In yeast, cytosolic acetyl-CoA is mainly produced by acetyl-CoA synthetases Acs1 and 2 [18]. It has been shown before that acetyl-CoA-dependent product formation can be increased by increasing acetyl-CoA formation [8, 32]. However, acetyl-CoA synthetases use hydrolysis of ATP to AMP which is not favorable for anaerobic n-butanol production, since only glycolysis and not respiration can provide ATP. Therefore, we tested whether the expression of an ATP independent acetylating acetaldehyde dehydrogenase from E. coli, adhE, can further increase n-butanol production. Normally, the enzyme adhE catalyzes the reduction of acetyl-CoA into ethanol. Membrillo-Hernandez and coworkers [33] succeeded in engineering the enzyme to mainly perform the reverse reaction of converting acetaldehyde into acetyl-CoA. The specific acetaldehyde dehydrogenase activity could be increased by changing two amino acids: A267T and E568 K. For this reason, we chose to use adhEA267T/E568K to increase the cytosolic acetyl-CoA pool. The mutant adhE A267T/E568K allele could complement the growth defect of an acs2 mutant strain on glucose medium, indicating that it is functional in yeast (data not shown) [21, 22]. We integrated this allele, expressed under control of the strong glycolytic PFK1 promoter, into the SFA1 locus resulting in strain VSY8. As SFA1 encodes an additional but minor alcohol dehydrogenase [34], we chose this locus for integration to stabilize the strain and to prevent the occurrence of ethanol-producing suppressor mutants. Surprisingly, fermentations of VSY8 expressing the n-butanol pathway (pVS6) showed almost no n-butanol production, but increased ethanol production instead (Fig. 5, Additional file 1: Figure S6). One reason could be that adhEA267T/E568K mainly converts acetaldehyde or acetyl-CoA into ethanol if there is not enough free CoA, the second substrate for the acetylating acetaldehyde dehydrogenase activity, available. Indeed, concomitant overexpression of coaA reduced ethanol and increased n-butanol formation (Additional file 1: Figure S9, Additional file 1: Figure S6). On the other hand, it might be possible that the additional alcohol dehydrogenase gene ADH6 is upregulated under those conditions. Indeed, when ADH6 was deleted by integration of coaA, a further decrease in ethanol formation could be observed [calculated yield for ethanol: 0.1 g/g glucose in VSY10 and 0.13 g/g in VSY8 (pVS6, pVS5_4)] (Additional file 1: Figure S6). In fermentations of the Δadh1-5, adhE A267T/E568K-overexpressing strain VSY8 (pVS6, pVS5_4), n-butanol production was increased up to 75 mg/L, which is a significant improvement compared to the strain without adhE A267T/E568K overexpression (52 mg/L) (Fig. 5). To further increase CoA availability, we added 25 ÂµM pantothenate which resulted in up to 95 mg/L n-butanol production (Fig. 5, Additional file 1: Figure S9). These results show that increasing the cytosolic acetyl-CoA synthesis by adhEA267T/E568K does improve n-butanol production, especially when combined with elevated CoA synthesis. It should be noted that although we used butyraldehyde dehydrogenase eutE in pVS6, which was also shown to convert acetaldehyde into acetyl-CoA [22], we could not observe a clear improvement in n-butanol production due to this activity.

Disruption of glycerol synthesis via chromosomal integration of the n-butanol pathway genes into GPD2 further increases n-butanol production

GPD2 encodes glycerol-3-phosphate dehydrogenase II which is the major enzyme for glycerol production in yeast under anaerobic conditions [35]. Glycerol production is used to re-oxidize a surplus of NADH. Deletion of GPD2 would further increase NADH levels in the adh deletion mutants [36] and could therefore push n-butanol production. Moreover, in order to overcome the negative effects of using different plasmids we decided to integrate coaA into the ADH6 locus of VSY8 (resulting in strain VSY9), and the n-butanol pathway genes of pVS6 into the GPD2 locus of strain VSY9, finally resulting in strain VSY10. Fermentations of strain VSY10 resulted in n-butanol titers up to 108 mg/L, which is significantly higher compared to its plasmid-based parental strain VSY8 (pVS6, pVS5_4; 75 mg/L) (Fig. 5). The increase might be due to both the stable integration of n-butanol pathway genes into the genome as well as reduced glycerol production [calculated yield for glycerol: 0.17 g/g glucose in VSY10 and 0.27 g/g in VSY8 (pVS6, pVS5_4)] (Additional file 1: Figure S6). Nevertheless, VSY10 did not consume all of the glucose (Additional file 1: Figure S5). This might be due to deletions of adh1-6, sfa1, and gpd2, resulting in inefficient NADH oxidation. The capacity of the introduced n-butanol pathway is obviously not strong enough to compensate this deficiency. This can also be seen by measured OD600 values. Without deletions, the wildtype is able to reach an OD600 of about 6, but all deletion strains reached not more than an OD600 of about 2, except for VSY8 (pVS6), that produced ethanol instead of n-butanol (Additional file 1: Figure S7). As a result of incomplete glucose consumption the n-butanol yield of VSY10 (0.01 g/g glucose) is much higher than for all fermentations with VSY8 (up to 0.0057 g/g) (Additional file 1: Figure S4). The addition of 25 ÂµM pantothenate increased the n-butanol yield of strain VSY10 up to 0.012 g/g (Additional file 1: Figure S4), resulting in a titer of 130 Â± 20 mg/L (Fig. 5, Additional file 1: Figure S10). This is to our knowledge the highest titer reported so far for n-butanol production in synthetic media in yeast via an acetoacetyl-CoA-derived n-butanol pathway.

Conclusions

Compared to bacterial systems, n-butanol production in S. cerevisiae via heterologous acetoacetyl-CoA-derived pathways is very low. The reason for this is not fully understood. Our results confirm that several reasons contribute to the low efficiency of n-butanol production in yeast. By stepwise engineering various aspects of n-butanol production, we could gradually increase n-butanol titers and yields. Of course, a very important aspect was the composition of the heterologous pathway. By testing various isoenzymes for the different reactions, we could establish a functional pathway. n-Butanol production was then increased by stable expression of the heterologous genes. Expression from centromeric plasmids was clearly superior to expression from multicopy plasmids (Table 2). Other important aspects were the elevation of acetaldehyde levels and NADH by reduction of alcohol dehydrogenase activities and glycerol synthesis. Moreover, increasing acetyl-CoA availability by means of a mutant form of an acetylating acetaldehyde dehydrogenase from E. coli further improved n-butanol production. However, our results indicate that obviously free CoA levels in S. cerevisiae are limiting n-butanol production. By increasing synthesis of CoA via overexpression of pantothenate kinase coaA of E. coli and by feeding pantothenate to the medium, n-butanol production levels could be increased up to 130 Â± 20 mg/L, the highest level ever reached in yeast via an acetoacetyl-CoA-derived pathway in synthetic medium.

Methods

Strains and media

All yeast strains used in this work are listed in Table 1. The adh deletion strains are based on strain JDY4 (MATa; ura3-52; trp1-289; leu2-3_112; his3 Δ1; MAL2-8C; SUC2 adh1::loxP adh3::loxP adh5::kanMX) (lab stocks Boles group), which was constructed from strain CEN.PK484 (EUROSCARF, Frankfurt).

S. cerevisiae was cultured in YEPD medium (10 g/L yeast extract, 20 g/L bacteriological peptone, 20 g/L glucose) or synthetic minimal (SMD) medium (1.7 g/L yeast nitrogen base without amino acids, 5 g/L ammonium sulfate, 20 mM monopotassium phosphate, supplemented with uracil (0.171 mM), leucine (0.439 mM), tryptophan (0.093 mM), and histidine (0.124 mM), 20 g/L glucose) at 30°. Synthetic media were adjusted to pH 6.3 with potassium hydroxide. The carbon source was autoclaved separately and added afterwards. If needed, the medium was supplemented with G418 (200 Î¼g/mL), hygromycin B (200 Î¼g/mL), and/or clonNAT/nourseothricin (100 Î¼g/mL) for selection of kanMX, hphNT1, and/or natNT2 markers, respectively.

Escherichia coli DH10b was used for cloning procedures and was grown in lysogeny broth (LB) at 37 Â°C with 100 Î¼g/ml ampicillin for plasmid selection.

Plasmid and strain construction

Plasmids were assembled via homologous recombination in yeast with up to 20 overlapping fragments with 50 bp overlaps. All heterologous genes were codon optimized according to the yeast glycolytic codon usage [37]. The mitochondrial targeting sequence of Eg ter was excluded [38] and the start codon of Ca bdhB and Sco ccr was changed from GTG into ATG. In Additional file 2: Table S2, a list of used genes with identifiers is available. Restriction sites were included in order to change sequences. If sequences within a plasmid were changed, the assembly of the corresponding backbone and the new sequence was performed via Gibson Assembly method [39]. Yeast was transformed according to protocols by Gietz and Schiestl [40, 41] with single fragments generated by PCR using primers shown in Additional file 2: Table S3. The assembled plasmids were recovered by yeast DNA preparations and transformed in E. coli for amplification by standard procedures. Restriction digestion analysis and DNA sequencing verified the constructs. The integrative plasmids were generated by successive excision of the 2 Î¼ or CEN6ARS4 origin and the AmpColE1 region and the cassette was integrated via homologous recombination with 400 bp overhangs. Marker rescue was carried out by the loxP/cre recombinase system with kanMX, hphNT1, or natNT2 markers. Genome integrations and deletions were confirmed by PCR analysis; primers are listed in Additional file 2: Table S3.

Fermentations and n-butanol production analysis

To test the n-butanol production of the different strains, pre-cultures were inoculated in YEPD media and grown aerobically to exponential phase. After washing steps, 40 mL SMD medium was inoculated with a starting OD600 of 0.3 and grown for at least 100 h. These batch fermentations were carried out semi-anaerobically with a rubber plug and a fermentation lock-sealed 100-mL flasks on a magnetic stirrer with 120 rpm at 30°. Samples were collected at 0, 6, 25, 50, 74, and 100 h for determination of cell density and for HPLC analysis through an inserted, sterile needle and syringe.

To quantify n-butanol concentrations, the samples were centrifuged and 450 Âµl of the supernatant was mixed with 50 Âµl 50 % (w/v) 5-sulfosalicylic acid and centrifuged again (10 min, 4 Â°C, 16,000×g). The samples were analyzed in a UHPLC + system by Thermo Scientific (Dionex UltiMate 3000) equipped with a HyperREZ XP Carbohydrate H + 8 Î¼m column and a refractive index detector (Thermo Shodex RI-101). The HPLC was operated at 65 Â°C with 5 mM sulfuric acid and a flow rate of 0.6 ml/min. n-Butanol, glucose, ethanol, glycerol, and pyruvate were used as standards with concentrations of 0.5 to 20 g/L. Additional standard concentrations for n-butanol were from 0.0125 up to 0.25 g/L.

Enzyme assays

For enzyme assays, CEN.PK113-5D with the respective plasmid was inoculated in 50 ml YEPD liquid medium with corresponding antibiotics and grown aerobically to an OD600 of 0.8–1. The cells were centrifuged and washed with water. The cell pellet was resuspended in 500 µl of the respective enzyme buffer and then lysed by the addition of 2/3 volumes glass beads (∅ 0.25–0.5 mm) and 10 min of shaking on a VXR basic Vibrax (IKA) at 2000 rpm. Cell debris was pelleted by centrifugation (10 min, 16,000×g, 4 Â°C) and the supernatant was used as a cell extract. Protein concentrations were measured by the Bradford method with a Roti-Quant Kit (Carl-Roth GmbH and Co). All enzyme assays were performed using the Ultrospec 2100 pro spectrophotometer (GE Healthcare, USA) at 30 Â°C under aerobic conditions. For each reaction, a volume of 200 Âµl was used, whereas 10 Âµl of cell extract with protein concentrations of 0.05–5 Âµg/µl was used. The determination of enzyme activity was carried out by calculating the slope of linear sections of absorbance versus time, using GE Healthcare\Datrys\Resolution and the respective molar extinction coefficients.

Thiolase (Erg10, thlA) activity was measured with 0.2 mM acetoacetyl-CoA and 0.2 mM CoA as substrates and the decrease in acetoacetyl-CoA concentration was measured at 303 nm. The reaction buffer contained 100 mM Tris–HCl, 10 mM MgCl2, and 1 mM DTT with pH 8. The reaction was started by the addition of cell extract. For determination of enzyme activity, the molar extinction factor ε = 14 mM−1cm−1 was used (based on [5, 15]).

Hydroxybutyryl-CoA dehydrogenase (hbd) activity was measured with 0.1 mM acetoacetyl-CoA as a substrate and 0.23 mM NADH as a cofactor. The decrease in NADH concentration at 340 nm was monitored. The reaction buffer contained 50 mM MOPS, 15 mM MgCl2, 1 mM EDTA, and 5 mM DTT with pH 7.3. The reaction was started by the addition of acetoacetyl-CoA [42].

Crotonase (crt) activity was measured with 0.1 mM crotonyl-CoA as a substrate. The decrease of absorption at 263 nm (disruption of the α-β unsaturation of crotonyl-CoA) was monitored. The reaction buffer contained 100 mM Tris–HCl, pH 7.6. The reaction was started by the addition of cell extract. For determination of enzyme activity, the molar extinction factor ε = 6.7 mM−1cm−1 was used [7].

Trans-2-enoyl-CoA reductase (ter) activity was measured with 0.5 mM crotonyl-CoA as a substrate and 0.46 mM NADH as a cofactor and 2 ÂµM FAD. The decrease in NADH concentration at 340 nm was monitored. The reaction buffer contained 100 mM potassium phosphate, pH 6.2. The reaction was started by the addition of crotonyl-CoA [38].

Butyraldehyde dehydrogenase (adhE2, eutE) activity was measured with 1 mM butyryl-CoA as a substrate and 300 ÂµM NADH as a cofactor. The decrease in NADH concentration at 340 nm was monitored. The reaction buffer contained 100 mM Tris–HCl and 5 mM DTT at pH 7.5. The reaction was started by the addition of butyryl-CoA [7].

Butanol dehydrogenase (adhE2, bdhB) activity was measured with 50 mM butyraldehyde as a substrate and 300 ÂµM NADH as a cofactor. The decrease in NADH concentration at 340 nm was monitored. The reaction buffer contained 100 mM Tris–HCl and 5 mM DTT at pH 7.5. The reaction was started by the addition of cooled butyraldehyde [7].

Acetyl-CoA:malonyl-CoA acyltransferase (nphT7) activity was measured with 0.4 mM acetyl-CoA and 0.4 mM malonyl-CoA as substrates and the increase in acetoacetyl-CoA concentration was measured at 303 nm. The reaction buffer contained 100 mM Tris–HCl, 10 mM MgCl2, and 1 mM DTT with pH 8. The reaction was started by the addition of cell extract. For determination of enzyme activity, the molar extinction factor ε = 14 mM−1cm−1 was used (based on [5, 28]).

Crotonyl-CoA reductase (ccr) activity was measured with 0.5 mM crotonyl-CoA as a substrate and 0.46 mM NADPH as a cofactor. The decrease in NADPH concentration at 340 nm was monitored. The reaction buffer contained 100 mM potassium phosphate, pH 6.2. The reaction was started by the addition of crotonyl-CoA (based on [38]).

Abbreviations

OD600 :

optical density at 600 nm

SMD:

synthetic minimal medium containing glucose

ADH:

alcohol dehydrogenase

References

  1. Generoso WC, Schadeweg V, Oreb M, Boles E. Metabolic engineering of Saccharomyces cerevisiae for production of butanol isomers. Curr Opin in Biotechnol. 2015;33:1–7.

    Article  CAS  Google Scholar 

  2. Huang H, Liu H, Gan Y. Genetic modification of critical enzymes and involved genes in butanol biosynthesis from biomass. Biotechnol Adv. 2010;28:651–7.

    Article  CAS  Google Scholar 

  3. Lee SY, Park JH, Jang SH, Nielsen LK, Kim J, Jung KS. Fermentative butanol production by Clostridia. Biotechnol Bioeng. 2008;101:209–28.

    Article  CAS  Google Scholar 

  4. Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC. Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008;10:305–11.

    Article  CAS  Google Scholar 

  5. Inui M, Suda M, Kimura S, Yasuda K, Suzuki H, Toda H, Yamamoto S, Okino S, Suzuki N, Yukawa H. Expression of Clostridium acetobutylicum butanol synthetic genes in Escherichia coli. Appl Microbiol Biotechnol. 2008;77:1305–16.

    Article  CAS  Google Scholar 

  6. Shen CR, Liao JC. Metabolic engineering of Escherichia coli for 1-butanol and 1-propanol production via the keto-acid pathways. Metab Eng. 2008;10:312–20.

    Article  CAS  Google Scholar 

  7. Shen CR, Lan EI, Dekishima Y, Baez A, Cho KM, Liao JC. Driving forces enable high-titer anaerobic 1-butanol synthesis in Escherichia coli. Appl Environ Microb. 2011;77:2905–15.

    Article  CAS  Google Scholar 

  8. Chen Y, Daviet L, Schalk M, Siewers V, Nielsen J. Establishing a platform cell factory through engineering of yeast acetyl-CoA metabolism. Metab Eng. 2013;15:48–54.

    Article  CAS  Google Scholar 

  9. Krivoruchko A, Serrano-Amatriain C, Chen Y, Siewers V, Nielsen J. Improving biobutanol production in engineered Saccharomyces cerevisiae by manipulation of acetyl-CoA metabolism. J Ind Microbiol Biotechnol. 2013;40:1051–6.

    Article  CAS  Google Scholar 

  10. Weber C, Farwick A, Benisch F, Brat D, Dietz H, Subtil T, Boles E. Trends and challenges in the microbial production of lignocellulosic bioalcohol fuels. Appl Microbiol Biotechnol. 2010;87:1303–15.

    Article  CAS  Google Scholar 

  11. Si T, Luo Y, Xiao H, Zhao H. Utilizing an endogenous pathway for 1-butanol production in Saccharomyces cerevisiae. Metab Eng. 2014;22:60–8.

    Article  CAS  Google Scholar 

  12. Branduardi P, Longo V, Berterame NM, Rossi G, Porro D. A novel pathway to produce butanol and isobutanol in Saccharomyces cerevisiae. Biotechnol Biofuels. 2013;6:68.

    Article  CAS  Google Scholar 

  13. Swidah R, Wang H, Reid PJ, Ahmed HZ, Pisanelli AM, Persaud KC, Grant CM, Ashe MP. Butanol production in S. cerevisiae via a synthetic ABE pathway is enhanced by specific metabolic engineering and butanol resistance. Biotechnol Biofuels. 2015;8:188.

    Article  Google Scholar 

  14. Dellomonaco C, Clomburg JM, Miller EN, Gonzalez R. Engineered reversal of the β-oxidation cycle for the synthesis of fuels and chemicals. Nature. 2011;476:355–9.

    Article  CAS  Google Scholar 

  15. Lian J, Zhao H. Reversal of the β-oxidation cycle in Saccharomyces cerevisiae for production of fuels and chemicals. ACS Synth Biol. 2015;4:332–41.

    Article  CAS  Google Scholar 

  16. Steen EJ, Chan R, Prasad N, Myers S, Petzold CJ, Redding A, Ouellet M, Keasling JD. Metabolic engineering of Saccharomyces cerevisiae for the production of n-butanol. Microb Cell Fact. 2008;7:36.

    Article  Google Scholar 

  17. Lian J, Si T, Nair NU, Zhao H. Design and construction of acetyl-CoA overproducing Saccharomyces cerevisiae strains. Metab Eng. 2014;24:139–49.

    Article  CAS  Google Scholar 

  18. Nielsen J. Synthetic biology for engineering acetyl coenzyme A metabolism in yeast. MBio. 2014;5:e02153.

    Article  CAS  Google Scholar 

  19. Strijbis K, Distel B. Intracellular acetyl unit transport in fungal carbon metabolism. Eukaryot Cell. 2010;9:1809–15.

    Article  CAS  Google Scholar 

  20. Kozak BU, van Rossum HM, Luttik MAH, Akeroyd M, Benjamin KR, Wu L, de Vries S, Daran J, Pronk JT, van Maris AJA. Engineering acetyl coenzyme A supply: functional expression of a bacterial pyruvate dehydrogenase complex in the cytosol of Saccharomyces cerevisiae. MBio. 2014;5:e01696.

    Article  CAS  Google Scholar 

  21. Boles E, Duvnjak JS, Festel G, Weber C. Gesteigerte Produktion von Acetyl-CoA. Patent. 2008; DE102008004253.

  22. Kozak BU, van Rossum HM, Benjamin KR, Wu L, Daran JG, Pronk JT, van Maris AJA. Replacement of the Saccharomyces cerevisiae acetyl-CoA synthetases by alternative pathways for cytosolic acetyl-CoA synthesis. Metab Eng. 2013;21:1051–6.

    Google Scholar 

  23. Olzhausen J, Schübbe S, Schüller H. Genetic analysis of coenzyme A biosynthesis in the yeast Saccharomyces cerevisiae: identification of a conditional mutation in the pantothenate kinase gene CAB1. Curr Genet. 2009;55:163–73.

    Article  CAS  Google Scholar 

  24. Stolz J, Sauer N. The fenpropimorph resistance gene FEN2 from Saccharomyces cerevisiae encodes a plasma membrane H + -pantothenate symporter. J Biol Chem. 1999;274:18747–52.

    Article  CAS  Google Scholar 

  25. Vadali RV, Bennett GN, San K. Cofactor engineering of intracellular CoA/acetyl-CoA and its effect on metabolic flux redistribution in Escherichia coli. Metab Eng. 2004;6:133–9.

    Article  CAS  Google Scholar 

  26. de Jong BW, Shi S, Valle-Rodríguez JO, Siewers V, Nielsen J. Metabolic pathway engineering for fatty acid ethyl ester production in Saccharomyces cerevisiae using stable chromosomal integration. J Ind Microbiol Biotechnol. 2015;42:477–86.

    Article  Google Scholar 

  27. Nielsen DR, Leonard E, Yoon S, Tseng H, Yuan C, Prather K. Engineering alternative butanol production platforms in heterologous bacteria. Metab Eng. 2009;11:262–73.

    Article  CAS  Google Scholar 

  28. Lan EI, Liao JC. ATP drives direct photosynthetic production of 1-butanol in cyanobacteria. Proc Natl Acad Sci USA. 2012;109:6018–23.

    Article  CAS  Google Scholar 

  29. Wallace KK, Bao Z, Dai H, Digate R, Schuler G, Speedie MK, Reynolds KA. Purification of crotonyl-CoA reductase from Streptomyces collinus and cloning, sequencing and expression of the corresponding gene in Escherichia coli. Eur J Biochem. 1995;233:954–62.

    Article  CAS  Google Scholar 

  30. Rock CO, Calder RB, Karim MA, Jackowski S. Pantothenate kinase regulation of the intracellular concentration of coenzyme A. J Biol Chem. 2000;275:1377–83.

    Article  CAS  Google Scholar 

  31. Mumberg D, Müller R, Funk M. Regulatable promoters of Saccharomyces cerevisiae: comparison of transcriptional activity and their use for heterologous expression. Nucleic Acids Res. 1994;22:5767–8.

    Article  CAS  Google Scholar 

  32. Shiba Y, Paradise EM, Kirby J, Ro D, Keasling JD. Engineering of the pyruvate dehydrogenase bypass in Saccharomyces cerevisiae for high-level production of isoprenoids. Metab Eng. 2007;9:160–8.

    Article  CAS  Google Scholar 

  33. Membrillo-Hernandez J. Evolution of the adhE gene product of Escherichia coli from a functional reductase to a dehydrogenase. Genetic and biochemical studies of the mutant proteins. J Biol Chem. 2000;275:33869–75.

    Article  CAS  Google Scholar 

  34. Weimer E, Rao E, Brendel M. Molecular structure and genetic regulation of SFA, a gene responsible for resistance to formaldehyde in Saccharomyces cerevisiae, and characterization of its protein product. Mol Gen Genet. 1993;237:351–8.

    Google Scholar 

  35. Ansell R, Granath K, Hohmann S, Thevelein JM, Adler L. The two isoenzymes for yeast NAD + -dependent glycerol 3-phosphate dehydrogenase encoded by GPD1 and GPD2 have distinct roles in osmoadaptation and redox regulation. EMBO J. 1997;16:2179–87.

    Article  CAS  Google Scholar 

  36. Kim S, Hahn J. Efficient production of 2,3-butanediol in Saccharomyces cerevisiae by eliminating ethanol and glycerol production and redox rebalancing. Metab Eng. 2015;31:94–101.

    Article  CAS  Google Scholar 

  37. Wiedemann B, Boles E. Codon-optimized bacterial genes improve L-arabinose fermentation in recombinant Saccharomyces cerevisiae. Appl Environ Microbiol. 2008;74:2043–50.

    Article  CAS  Google Scholar 

  38. Hoffmeister M. Mitochondrial trans-2-enoyl-CoA reductase of wax ester fermentation from Euglena gracilis defines a new family of enzymes involved in lipid synthesis. J Biol Chem. 2004;280:4329–38.

    Article  Google Scholar 

  39. Gibson DG, Young L, Chuang R, Venter JC, Hutchison CA, Smith HO. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6:343–5.

    Article  CAS  Google Scholar 

  40. Gietz RD, Schiestl RH. Frozen competent yeast cells that can be transformed with high efficiency using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007;2:1–4.

    Article  CAS  Google Scholar 

  41. Gietz RD, Schiestl RH. High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007;2:31–4.

    Article  CAS  Google Scholar 

  42. Han L, Reynolds KA. A novel alternate anaplerotic pathway to the glyoxylate cycle in streptomycetes. J Bacteriol. 1997;179:5157–64.

    CAS  Google Scholar 

  43. Taxis C, Knop M. System of centromeric, episomal, and integrative vectors based on drug resistance markers for Saccharomyces cerevisiae. Biotechniques. 2006;40:73–8.

    Article  CAS  Google Scholar 

  44. Farwick A, Bruder S, Schadeweg V, Oreb M, Boles E. Engineering of yeast hexose transporters to transport D-xylose without inhibition by D-glucose. Proc Natl Acad Sci USA. 2014;111:5159–64.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

VS and EB conceived the study. VS conducted the experiments. VS and EB analyzed the data. VS and EB wrote the paper. Both authors read and approved the manuscript.

Acknowledgements

We thank Dr. Mislav Oreb (Frankfurt, Germany) for helpful advice. This work has been financially supported by BASF SE.

Competing interests

The authors declare they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Eckhard Boles.

Additional files

13068_2016_456_MOESM1_ESM.pdf

Additional file 1: Figure S1. n-Butanol production with different variants of acetoacetyl-CoA derived n-butanol pathways in S. cerevisiae strain CEN.PK113-5D. Figure S2. Comparison of n-butanol production with and without overexpression of coaA (pVS5_4) and adhE A267T/E568K (pVS4) in VSY0 (Δadh1-5). Figure S3. n-Butanol production with n-butanol pathway (pVS6) and FEN2 overexpression under control of the MET25 promoter (pVS5_6) in CEN.PK113-5D with addition of methionine. Figure S4. Comparison of n-butanol yields of different strains. Figure S5. Consumption of glucose in semi-anaerobic fermentations started with 20 g/L glucose. Figure S6. Ethanol and glycerol production in semi-anaerobic fermentations. Figure S7. OD600 after 74 h in semi-anaerobic fermentations. Figure S8. n-Butanol production with n-butanol pathway (pVS6) and coaA overexpression (pVS5_4) in VSY0 (Δadh1-5) with addition of pantothenate. Figure S9. n-Butanol production with n-butanol pathway (pVS6) and coaA overexpression (pVS5_4) in VSY8 (Δadh1-5 sfa1, integrated adhE A267T/E568K) with addition of pantothenate. Figure S10. n-Butanol production of strain VSY10 (Δadh1-5 sfa1 adh6 gpd2 with integrated adhE A267T/E568K, coaA and n-butanol pathway) with addition of pantothenate.

13068_2016_456_MOESM2_ESM.pdf

Additional file 2: Table S1. Yeast promoters and terminators used for the expression of n-butanol pathway genes in this study. Table S2. Genes used in this study and their corresponding identifiers. Table S3. Relevant primers for this study. Table S4. Statistical analysis of n-butanol production.

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schadeweg, V., Boles, E. n-Butanol production in Saccharomyces cerevisiae is limited by the availability of coenzyme A and cytosolic acetyl-CoA. Biotechnol Biofuels 9, 44 (2016). https://doi.org/10.1186/s13068-016-0456-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-016-0456-7

Keywords