Skip to main content

Anaerobic poly-3-d-hydroxybutyrate production from xylose in recombinant Saccharomyces cerevisiae using a NADH-dependent acetoacetyl-CoA reductase



Poly-3-d-hydroxybutyrate (PHB) that is a promising precursor for bioplastic with similar physical properties as polypropylene, is naturally produced by several bacterial species. The bacterial pathway is comprised of the three enzymes β-ketothiolase, acetoacetyl-CoA reductase (AAR) and PHB synthase, which all together convert acetyl-CoA into PHB. Heterologous expression of the pathway genes from Cupriavidus necator has enabled PHB production in the yeast Saccharomyces cerevisiae from glucose as well as from xylose, after introduction of the fungal xylose utilization pathway from Scheffersomyces stipitis including xylose reductase (XR) and xylitol dehydrogenase (XDH). However PHB titers are still low.


In this study the acetoacetyl-CoA reductase gene from C. necator (CnAAR), a NADPH-dependent enzyme, was replaced by the NADH-dependent AAR gene from Allochromatium vinosum (AvAAR) in recombinant xylose-utilizing S. cerevisiae and PHB production was compared. A. vinosum AAR was found to be active in S. cerevisiae and able to use both NADH and NADPH as cofactors. This resulted in improved PHB titers in S. cerevisiae when xylose was used as sole carbon source (5-fold in aerobic conditions and 8.4-fold under oxygen limited conditions) and PHB yields (4-fold in aerobic conditions and up to 5.6-fold under oxygen limited conditions). Moreover, the best strain was able to accumulate up to 14% of PHB per cell dry weight under fully anaerobic conditions.


This study reports a novel approach for boosting PHB accumulation in S. cerevisiae by replacement of the commonly used AAR from C. necator with the NADH-dependent alternative from A. vinosum. Additionally, to the best of our knowledge, it is the first demonstration of anaerobic PHB synthesis from xylose.


Most of today’s energy and materials are primarily coming from petroleum based resources that are non-renewable. Petroleum derivatives such as plastics are one example of an important commodity in contemporary society, with a world annual plastic production in 2014 of 311 million tons and that is estimated to double again in the next 20 years [1]. An alternative could be to use biopolymers that are produced in nature, for instance isoprenoids and polyhydroxyalkanoates (PHAs). Of the PHAs, poly-3-d-hydroxybutyrate (PHB) has been suggested to be one of the most promising biopolymers, as it shares many physical properties with polypropylene [2].

Several microorganisms are capable of accumulating intracellular granules of PHB through polymerisation of soluble molecules, hence preventing leakage of valuable compounds out of the cells. In both gram-positive and gram-negative bacteria, PHB is used as carbon storage when nutrient supplies are imbalanced [3]. Cupriavidus necator (formerly known as Ralstonia eutropha) is one of the most studied bacteria for its ability to accumulate considerable amounts of PHB [4]. The biosynthesis of PHB in this bacterium is catalyzed through three enzymatic reactions where the first reaction allow the condensation of two acetyl coenzyme A (acetyl-CoA) molecules into acetoacetyl-CoA by β-ketoacyl-CoA thiolase (encoded by phaA). The next step is the conversion of acetoacetyl-CoA to (R)-3-hydroxybutyryl-CoA catalyzed by the acetoacetyl-CoA reductase (AAR) (encoded by phaB). Lastly, the (R)-3-hydroxybutyryl-CoA monomers are polymerized into PHB by the polyhydroxyalkanoate synthase (encoded by phaC). All three enzymes required for PHB synthesis in bacteria are located in the cytoplasm of the cell where PHB granule accumulation takes place [5, 6]. Most characterized PHB-producers carry an NADPH-dependent AAR gene, whereas AAR with putative promiscuous cofactor utilization are only reported in few cases: for the halotolerant bacterium Allochromatium vinosum, formerly known as Chromatium vinosum [6], the anaerobic syntrophic bacterium Syntrophomonas wolfei [7] and the recently described AAR from the halophilic bacterium Halomonas boliviensis [8].

However the economic feasibility of PHB bio-production is dependent on using a cheap and highly available source of fermentable sugars, such as forestry and agriculture residues. During lignocellulosic biomass pretreatment, several inhibitory compounds are released [9], generating a harsh environment in which PHB-producing bacteria species are not well adapted. Therefore the industrial workhorse Saccharomyces cerevisiae that is outcompeting other cell factories for inhibitor tolerance [10, 11] has been explored as an alternative for the production of PHAs, possibly as a valuable side-product from the lignocellulosic ethanol anaerobic process. The first attempt to produce PHB in S. cerevisiae was from the available cytosolic 3-hydroxybutyryl-CoA derived from fatty acid β-oxidation pathway through the overexpression of phaC gene from C. necator [12]. Later, PHB production was detected from glucose by expressing the three bacterial genes phaA, phaB and phaC from C. necator [13, 14]. In addition the pool of cytosolic acetyl-CoA was increased in S. cerevisiae by overexpressing the alcohol dehydrogenase (ADH2), the acetaldehyde dehydrogenase (ALD6), the acetyl-CoA acetyltransferase (ERG10) and the Salmonella enterica acetyl-CoA synthetase variant (acs L641P) and it resulted in improved PHB productivity from glucose [15]. Finally since a major fraction of sugars from lignocellulosic biomass may consist of xylose, PHB pathway has been introduced in S. cerevisiae engineered for xylose utilization, leading to PHB synthesis from xylose [16]. In the present study we demonstrate that PHB can be produced anaerobically, and in combination with ethanol, from xylose by cofactor shift through the introduction of the NADH-dependent AAR alternative from A. vinosum in recombinant S. cerevisiae.



Plasmids used in the study are presented in Table 1. For the construction of YIpAGS3, the AAR gene from A. vinosum (GenBank Accession No. YP_003442070.1) was codon optimised for S. cerevisiae at Eurofins Genomics (Ebersberg, Germany) and synthesized by Genescript (Piscataway, NJ, USA). The coding sequence was designed to be under the control of the constitutive promoter-terminator pair of the gene TPI1 flanked by the restriction sites SacI and SpeI. The custom-synthesized coding sequence was cloned on the integrative plasmid YIpAGS2 [16] generating YIpAGS3 (Additional file 1).

Table 1 Plasmids used in the study

Strains, media and culture conditions

Yeast strains used in the study are listed in Table 2. Yeast strains were recovered, from 20% glycerol stocks stored at −80 °C, on solid YPD (10 g/L yeast extract, 20 g/L peptone, 20 g/L glucose, 15 g/L agar) for two days at 30 °C. Yeast cultures were grown in liquid YPD medium for 14–16 h, or less when required, at 30 °C and 180 rpm in an orbital shaker. Competent yeast cells were prepared and transformed according to the High efficiency transformation protocol [17]. Transformants were selected on solid Yeast Nitrogen Base (YNB) medium (6.7 g/L Yeast Nitrogen Base without amino acids (Becton, Dickinson and Company, USA) supplemented with 20 g/L glucose or xylose and 15 g/L agar). Leucine was added for complementation at a concentration of 30 mg/L when required.

Table 2 Saccharomyces cerevisiae strains used in the study

Escherichia coli strain NEB5-α (New England Biolabs) was used for sub-cloning of plasmid DNA. Heat shock competent E. coli cells were prepared according to the Inoue method [18] and transformed according to the supplier’s instructions. Transformants were selected on solid Luria–Bertani (LB) plates (5 g/L yeast extract, 10 g/L peptone, 5 g/L NaCl, 15 g/L agar, pH 7.0), supplemented with 100 mg/L of ampicillin, for 16 h at 37 °C. Cultures of transformed E. coli were recovered from 25% glycerol stocks stored at −80 °C and grown in liquid LB medium, supplemented with ampicillin 50 mg/L, for 14–16 h at 37 °C and 180 rpm in an orbital shaker.

Strain TMB3043 [19], an engineered strain overexpressing the non-oxidative pentose phosphate pathway for efficient pentose utilization, was used as background strain in this study. TMB3043 was transformed with the linearized vectors YIpDR7 or the YIpOB8 [20], generating the strains TMB4420 and TMB4440, respectively (Table 2). TMB4420 and TMB4440 were subsequently transformed with the linearized vector YIpAGS3 generating the strains TMB4425 and TMB4445, respectively. TMB4424 was obtained from transformation of TMB4420 with linearized YIplac128. The YIpDR7 and YIpOB8 plasmids were digested with Fast Digest EcoRV (Thermo Scientific), YIplac128 with Fast Digest Eco9I and YIpAGS3 with Fast Digest PpuMI (Thermo Scientific, USA). Transformed strains TMB4425 and TMB4445 were validated by amplifying stretches of the integrated cassette by PCR, using extracted genomic DNA as template and primers PhaC1_f: 5′CATATTACAATAATGGCCACTGGTAAAGG3′ and PhaC1_r: 5′CATTCATTCTTCAGACTTATGCCTTTGCTTTCACATAC 3′.

Enzymatic assays

Cells were cultivated in 50 mL conical centrifuge tubes containing 5 mL of YNB-glucose, at 30 °C 180 rpm overnight. At starting OD620 of 1.0 the cells were inoculated in the same media as described above and harvested during the late exponential phase. Whole-cell protein extract was obtained by treating the cells with Y-PER extraction solution (Pierce, Rockford, IL, USA) according to the manufacturer's instructions. The total protein concentrations were measured with a Bradford assay using Micro BCA™ Protein Assay (ThermoFisher Scientific, Waltham, MA, USA) and bovine albumin standard (ThermoFisher Scientific, Waltham, MA, USA). Absorbance was measured at 340 nm on a Multiskan Ascent (Thermo Electro Corporation, Finland) using 96-well plates (final volume of 250 µL, in three replicates). Acetoacetyl-CoA reductase kinetics was determined using 200 µM acetoacetyl-CoA as substrate, 10 µL of cell extract, 200 mM NADH or NADPH as cofactor, 100 mM MOPS pH 7.0 as buffering agent and 5 mM 2-mercaptoethanol as reducing additive agent. The measured initial rates were used to calculate specific activities (µmol/min mg total protein).

Aerobic cultivations in shake flasks

Single colonies were pre-inoculated in 5 mL YNB-xylose (13.4 g/L YNB without amino acids, 50 g/L xylose, 50 mM pH 5.5 phthalic acid) in 50 mL conical centrifuge tubes at 30 °C, 180 rpm overnight. Cells were harvested during the exponential phase and inoculated in 250 mL baffled shake flasks containing 50 mL buffered YNB-xylose, at starting OD620 of 0.05 and set at 180 rpm at 30 °C. The growth rate was followed by measuring optical density at 620 nm with an Ultrospec 2100 pro spectrophotometer (GE Healthcare Life Sciences, Sweden). The measurements were made in technical duplicates. Cultivations were performed in biological duplicates.

Oxygen-limited cultivations

Oxygen-limited cultivations were performed in 65 mL serum flasks containing 50 mL of buffered YNB-xylose. For oxygen limited growth, media was supplemented with Tween 80 (400 mg/L) and ergosterol (10 mg/L) and 7 mL of mineral oil on top. Temperature was maintained at 30 °C in a water bath and agitation with a magnetic stirrer was set to 180 rpm. Inoculated cells were pre-grown under aerobic conditions to allow biomass formation for starting at OD620 of 0.1. Batch experiments were performed in biological duplicates.

Anaerobic batch fermentations in bioreactors

Fully anaerobic fermentations were performed in Multifors 1.4 L vessels (Infors, Switzerland) containing 1000 mL unbuffered YNB-xylose supplemented with Tween 80 (400 mg/L) and ergosterol (10 mg/L). Bioreactors were initially inoculated with freshly grown cells to a starting OD620nm of 0.5. The cells used for the inoculum were grown aerobically as described above in 50 mL conical centrifuge tubes and then transferred to 250 mL baffled shake-flasks to get higher levels of biomass. The bioreactors were sparged with 0.1 vvm nitrogen at 30 °C and 200 rpm. Silicon antifoam RD emulsion (Dow Corning, USA) was added to avoid excess foam formation. The pH was kept at 5.5 with the addition of 3M KOH and 3M H2SO4 solutions. The experiments were run independently twice.

Metabolite determination and cell dry weight analysis

Cultures were regularly sampled for metabolite analysis. Extracellular metabolite levels were measured from the sample supernatants with HPLC using a Waters (Milford, USA) system equipped with an Aminex HPX-87H column (Bio-Rad, Richmond, USA) that operates at 60 °C with a mobile phase of 5 mM H2SO4 and 0.6 mL/min flow rate. Concentrations of glucose, xylose, xylitol, glycerol, acetate and ethanol were calculated from an external standard calibration curve. Samples were analysed in technical duplicates. Cell dry weight was determined by filtering samples through a pre-weighed 0.45 µm membrane filter (Pall Corporation, New York, USA) The membranes were washed with ultrapure H2O and dried in a microwave oven at 350 W for 8 min. The final weight was measured after equilibration to room temperature in a desiccator. The measurements were made in technical duplicates. In the case of the bioreactor fermentation, biomass was measured at each time point.

PHB quantification

The method used relied on the quantitative conversion of PHB to crotonic acid catalysed by hot concentrated sulfuric acid [21]. PHB content was analysed by harvesting 2 mL of culture by centrifugation for 5 min at 6000×g. Supernatant was carefully removed and the pellets were resuspended in ultrapure H2O, followed by centrifugation for 5 min at 6000×g. This step was repeated twice, then 0.5 mL of 99% sulfuric acid was added onto the pellets that were incubated at 95 °C for 1 h with open lids in a heat block (Grant QBD1, Grant Instruments, UK). The resulted solution was diluted 20 times by serial dilution and further analysed by HPLC using the same conditions as for the metabolite analysis. Commercially available PHB (#363502, Sigma-Aldrich) processed in parallel with the samples was used as control. Concentrations of crotonic acid were calculated from an external ten-point standard calibration curve. Samples were analysed in technical duplicates.

Enzyme structure homology modelling

Protein structure models were generated with automated homology modeling using SWISS-MODEL [22] by submitting translated gene sequences. The translated AAR gene from A. vinosum (GenBank Accession No. YP_003442070.1) was submitted, using the C. necator AAR structure as designated template (PDB: 3VZP) [23].


Selection of putative NADH-dependent acetoacetyl-CoA reductase

The AAR from C. necator uses NADPH as a cofactor [24], which competes with biomass synthesis in recombinant S. cerevisiae and could explain the low PHB titers observed on glucose and xylose [15, 16]. This led us to look for NADH-dependent alternatives in reported studies. Three AARs appeared to be promising candidates: the AAR from A. vinosum (GenBank: YP 003442070.1) that is a thiosulfate oxidizing, gram-negative bacterium accumulating PHB up to 58% of its CDW [25] and two AARs from S. wolfei species adapted for anaerobic synthropic growth [GenBank: ABI67978.1 (Swol_0651) and ABI69207.1 (Swol_1910)] [7, 26]. Since the crystal structures of these AARs were not available, protein structures were made with automated homology modeling using SWISS-MODEL [22]. The AvAAR displayed 56% protein sequence identity to the C. necator AAR whose crystal structure was used as template (PDB: 3VZP) [23]. The S. wolfei AARs had a protein sequence identity of 30% (GenBank ABI67978.1) and 36% (GenBank ABI69207.1) respectively, compared to the CnAAR. Both S. wolfei homology models (Swol_0651 modelled on PDB: 4MOW [27] and Swol_1910 modelled on PDB: 4DMM [28]) indicated a NADPH-binding arginine motif [23], typical of the CnAAR and related AARs. In contrast, the cofactor binding site of the AvAAR contained a dominant glutamate instead of an arginine (compared to the C. necator AAR) (Fig. 1). Acidic residues, such as glutamate in the case of the AvAAR, can create hydrogen bonds to the 2′- and 3′-hydroxyl groups on the adenine ribose of NADH [29], in comparison to arginine, which are ideally forming contacts to the 2′-phosphate in NADPH [23]. The AvAAR was also assumed to generate the correct enantiomer, (R)-3-hydroxybutyryl-CoA, preferred by the C. necator polyhydroxyalkanoate synthase for PHB accumulation as the binding site for the prochiral keto functionality on the acetoacetyl-CoA substrate was almost identical in the AvAAR homology model as in the CnAAR crystal structure. From the results above, the AvAAR was further used and evaluated in a recombinant PHB-producing S. cerevisiae host.

Fig. 1
figure 1

Cofactor binding sites of the AAR from A. vinosum, made through homology modeling (left) and the AAR from C. necator (PDB: 3VZP) (right). The cofactor NADPH is modeled in the binding site and Glu-37 is marked in the A. vinosum AAR, and the Arg-40, which forms several polar contacts to NADPH, is marked in the C. necator AAR

Construction of NADH-dependent PHB pathways in recombinant xylose-utilising strains

The previously described PHB-producing xylose-utilizing S. cerevisiae strains carried the integrative yeast plasmid YIpAGS2 that contains codon-optimised genes encoding the β-ketothiolase, acetoacetyl-CoA reductase (AAR) and PHB synthase from C. necator under the control of constitutive promoters [16]. In this study YIpAGS3 (Table 1; Additional file 1) was constructed by replacing the AAR from C. necator (CnAAR) in YIpAGS2 with the codon-optimised homolog from A. vinosum (AvAAR). YIpAGS3 was then introduced in two xylose-consuming strain backgrounds: TMB4420 that carries a Asn272Asp substitution in the xylose reductase (XRmut) from Scheffersomyces stipitis [20] and TMB4440 that carries a wild type xylose reductase (XRwt) from S. stipitis [30]. These two alternatives were chosen because XRwt is known to have a nearly exclusive NADPH-selectivity, consequently depletion of NADPH in the xylose reduction step is expected to impact both growth and NADPH-dependent PHB production in xylose-rich media [16, 20]. In contrast, XRmut contains the N272D mutation in the cofactor binding that increases the enzyme specificity for NADH [20]. This is known to result in an improvement in the cofactor balance for xylose uptake (Fig. 2) that translates into enhanced anaerobic growth [20]. But the effect on PHB production is unknown.

Fig. 2
figure 2

Metabolic diagram of the main central carbon metabolism including the XR and PHB pathway steps. XR xylose reductase, NADH nicotinamide adenine dinucleotide, NADPH nicotinamide adenine dinucleotide phosphate, ATP adenosine triphosphate, AMP adenine monophosphate, TCA tricarboxylic acid cycle, PHB poly-3-D-hydroxybutyrate, AAR acetoacetyl-CoA reductase, Ac-CoA acetyl coenzyme A, G3P glyceraldehyde 3-phosphate, PPP pentose phosphate pathway

The integrative empty vector YIplac128 and the newly constructed YIpAGS3, were used to transform TMB4420 (XRmut), generating strains TMB4424 (XRmut) and TMB4425 (XRmut, AvAAR), respectively. YIpAGS3 was also used to transform TMB4440 (XRwt), generating strain TMB4445 (XRwt, AvAAR). The strains were then compared with the control strain TMB4444 (XRwt) and with TMB4443 (XRwt, CnAAR) [16].

In vitro activities of acetoacetyl-CoA reductases expressed in S. cerevisiae

In order to check whether the introduced A. vinosum AAR gene was functional in S. cerevisiae, enzymatic assays were performed in strains TMB4443 (XRwt, CnAAR) and TMB4445 (XRwt, AvAAR) carrying C. necator and A. vinosum AAR, respectively (Table 3).

Table 3 Acetoacetyl-CoA reductase (AAR) specific activity

The expressed AAR from A. vinosum showed specific activity for both NADH and NADPH cofactors but it displayed 4.5-fold higher specific activity for NADH than for NADPH. In contrast the C. necator AAR showed significant activity for NADPH only. In addition, the NADPH-dependent activity was 4.3-fold higher for A. vinosum than for C. necator AAR.

Altogether the results demonstrated that the newly expressed AAR from A. vinosum was functionally expressed in S. cerevisiae and had a preference for NADH as cofactor.

Aerobic PHB synthesis from xylose in shake flasks

The strains carrying the NADH-dependent AAR, i.e. TMB4425 (XRmut, AvAAR) and TMB4445 (XRwt, AvAAR), were studied under aerobic conditions in xylose using shake flasks, i.e. under the exact same conditions as the previously described strain TMB4443 (XRwt, CnAAR) [16]. Figure 3 shows a representative pattern for growth, xylose consumption, ethanol, acetate, glycerol formation and PHB production for the three strains. For comparison purposes, the control strains TMB4424 (XRmut) and TMB4444 (XRwt) that were previously cultivated under the same conditions [16], were also included in the result summary that is displayed in Additional file 2.

Fig. 3
figure 3

Aerobic growth and metabolite profiles from xylose for recombinant S. cerevisiae strains carrying a NAD(P)H-dependent PHB pathway. a TMB4443 (XRwt, CnAAR) [16] (rights and permission from AMB Express), b TMB4445 (XRwt, AvAAR) and c TMB4425 (XRmut, AvAAR). Strains were cultivated in biological duplicates on xylose in defined buffered medium using baffled shake flask. Results display a representative growth from two biological duplicates. Error bars represent standard deviation

PHB production was detected during the exponential growth in all strains carrying a PHB pathway (Fig. 3). The results of the previously described strain that carries the XRwt and the NADPH-dependent PHB pathway (TMB4443) are shown in Fig. 3a. In comparison, the newly generated strain harboring XRwt and the NADH-dependent PHB pathway (TMB4445, Fig. 3b) displayed significantly higher PHB yield (4.3-fold), PHB titer (5.1-fold) and final PHB content per CDW (3.4-fold). When comparing strains carrying the NADH-dependent PHB pathway, PHB titers were 1.3-fold higher in TMB4445 (XRwt) than in TMB4425 (XRmut) but the final PHB content per CDW was 1.4-fold higher in TMB4425 than in TMB4445, as result of less xylose being used in the XRmut strain (Fig. 3c; Additional file 2).

More generally strains harboring XRwt (TMB4443, TMB4444, and TMB4445) had a higher growth rate than the ones carrying XRmut (TMB4424, TMB4425). Xylitol was the main by-product from xylose fermentation and it was 2.5-fold higher in TMB4445 (XRwt) than in TMB4425 (XRmut). In both strains the glycerol yield remained low, while ethanol yields were lower in TMB4425 (XRmut) than in TMB4445 (XRwt) which correlates to a higher yield in acetate and PHB for TMB4425.

Oxygen-limited cultivations

To investigate the impact of the redox balance on PHB synthesis of the cell, cultivations were performed for the five strains under oxygen-limiting conditions using serum flasks (Additional file 2). Figure 4 shows the representative growth, xylose consumption, ethanol, acetate, glycerol formation and PHB during growth in oxygen-limited serum flasks for the three strains expressing a PHB pathway (TMB4443, TMB4445 and TMB4425).

Fig. 4
figure 4

Oxygen-limited growth and metabolite profiles on xylose for recombinant S. cerevisiae strains carrying a PHB-pathway. a TMB4443 (XRwt, CnAAR). b TMB4445 (XRwt, AvAAR). c TMB4425 (XRmut, AvAAR). Strains were cultivated in biological duplicates on xylose in defined buffered medium using serum flasks. Results display a representative graph from two biological duplicates. Error bars represent standard deviation

TMB4425 (XRmut, AvAAR) reached the highest PHB titer, 252 mg/L, followed by TMB4445 (XRwt, AvAAR), 131 mg PHB/L, and TMB4443 (XRwt, CnAAR), 29 mg PHB/L. However, since the strains carrying XRwt grew less efficiently on xylose under oxygen-limited conditions, the final PHB content per biomass was higher for TMB4445 (XRwt, AvAAR) (14.7%) than for TMB4425 (XRmut, AvAAR) (9.4%). In all cases, PHB was detected after 48 h of cultivation; also the strains showed higher biomass yields than their respective controls lacking the PHB pathway, which might be attributed to an increase of the cellular weight as a result of internal accumulation of PHB granules (Additional file 2).

At the metabolite level under oxygen-limiting conditions, the xylitol titer was higher for the strains expressing the XRwt, TMB4443 (7.63 ± 1.40 g) and TMB44445 (6.81 ± 0.09 g) than for the strain expressing the XRmut, TMB4425 (5.34 ± 0.15 g), which also correlates with their respective xylitol yields (0.28 ± 0.02; 0.25 ± 0.01 and 0.15 ± 0.01 g/g). This results confirmed what previous works showed about the importance of a more balanced co-factor ratio in the xylose metabolisation steps (see e.g. [3133]). In particular, strains expressing XRmut have higher affinity for NADH than XRwt and therefore demand less usage of NADPH for xylose reduction thereby increasing the availability of NADPH for biomass formation and a cofactor balanced xylose conversion to xylulose [34].

Ethanol was produced and accumulated during exponential growth, but it was not re-assimilated. Final ethanol yields were slightly lower in the strain carrying XRmut and the NADH-PHB pathway (TMB4425) than its respective control strain (TMB4424). Finally, glycerol and acetate levels remained low during all the cultivation time (Additional file 2).

Anaerobic batch fermentations

The strain that displayed the most promising results in terms of PHB production, i.e. TMB4425 carrying a NADH-dependent PHB pathway and XRmut, was further evaluated in bioreactor under fully anaerobic conditions (Fig. 5).

Fig. 5
figure 5

Representative anaerobic growth and metabolite profiles for recombinant S. cerevisiae TMB4425 (XRmut, AvAAR) strain on xylose defined medium at pH 5.5 using controlled bioreactors. Error bars represent standard deviation from two independent biological duplicates

TMB4425 grew at a maximum growth rate of 0.07/h, i.e. at similar values that those were previously obtained by Runquist and coworkers [20]. High xylose conversion efficiency was achieved, with only 2.1% of the total xylose remaining at the end of the fermentation. This resulted in a higher biomass level than under oxygen-limited conditions (Additional file 2). Ethanol was the most abundant metabolite, reaching a maximum titer of 16.4 g/L. After 150 h, 360 mg PHB/L was obtained, corresponding to a 1.4- and 2.0-fold increase in PHB production as compared to the oxygen-limited and aerobic cultivations, respectively. It also corresponded to the highest PHB content per biomass (14.2% of PHB of the total CDW) and the highest PHB yield per xylose consumed for this strain. The acetate, xylitol and glycerol yields were lower than under oxygen-limited conditions, suggesting that the carbon flux increased towards PHB production.


The present study proposes a metabolic engineering strategy to enhance the production of PHB from xylose in recombinant S. cerevisiae. The bacterial pathway for PHB production from C. necator has previously been integrated into xylose-utilising S. cerevisiae strains as part of a strategy to convert xylose into PHB [16]. Here we show that replacing the NADPH-dependent acetoacetyl-CoA reductase (AAR) from C. necator with the NADH-dependent AAR from A. vinosum is essential for boosting PHB production from xylose under aerobic and oxygen-limiting conditions. We also obtain anaerobic PHB production from xylose and the highest ever reported PHB content per cell in S. cerevisiae, in a strain combining the new AAR with a co-factor balanced xylose pathway. These results offer interesting insights for the potential of industrial anaerobic ethanol-PHB co-production from lignocellulosic feedstocks and other agricultural/industrial residues because during the anaerobic processes the metabolic flux for synthesis of ethanol and PHB is favored. In addition the need for sparging the reactor with air is no longer necessary, leading to a reduction of investments in equipment and energy supply that are necessary to increase oxygen transfer and cooling [35].

In native organisms PHB is a storage carbohydrate generated under starvation conditions, which implies that its production is growth-uncoupled [36]. This may explain why most PHB pathways can use NADPH-dependent AAR as it does not have to compete with the NADPH-demanding biosynthetic reactions. In S. cerevisiae, however, this is a drawback because the deregulated PHB pathway competes with biosynthesis for NADPH. Also, cytosolic availability of NADPH is known to be lower than for NADH [37]. Altogether this may contribute to the limited PHB yields that have been observed in this organism [13, 15, 16]. Using instead a NADH-dependent AAR should enable PHB production that is less dependent from biomass formation than with the regular C. necator pathway. Indeed, we were able to show that PHB volumetric titers, yields and PHB content per cell could be improved by using the NADH-dependent AAR from A. vinosum. It cannot be ruled out that the improvement may result from the overall higher NAD(P)H acetoacetyl-CoA conversion rate in the novel AAR. However the increasing PHB content per cell under oxygen-limiting conditions suggests that the availability of intracellular NADH plays a key role in the overall improvement.

A growth-coupled PHB production also implies that part of the cytosolic acetyl-CoA that is needed for biosynthesis has to be directed towards PHB formation. This can explain why PHB production is found to be inversely proportional to biomass formation and growth rate under all tested conditions. For example, TMB4425 carrying the NADH-favoring XRmut has a lower growth rate than TMB4445 (carrying the NADPH-favoring XRwt) under aerobic conditions, as a result of limiting NADH availability under respiratory conditions, and it accumulates 1.4 times more PHB per cell; instead, when the oxygen availability is limiting, TMB4425 grows faster than TMB4445 (XRwt) due higher availability of NADH but the PHB content per cell is one-third lower than the one observed in TMB4445 (XRwt). More generally a reduction of the carbon flux towards biomass synthesis via the TCA cycle is expected to result in an increase in the flow towards ethanol and acetate, thereby increasing the levels of PHB precursors.

The choice of XR is also critical for anaerobic PHB production from xylose. We show that a combination of XRmut and NADH-dependent AAR, i.e. TMB4425, should be used under these conditions to produce significant amounts of PHB because they favour NADH excess, which benefits both XRmut and AvAAR activity, as well as a shift towards the fermentative pathway, thereby the PHB precursors, as a result of the absence of mitochondrial respiration. Still ethanol remains the major by-product of xylose fermentation, therefore it will be necessary to further increase the carbon flow towards the acetate and acetyl-CoA precursors, for instance by combining the current changes with the reported acetyl-CoA pool engineering strategies consisting in expressing heterologous pathways [3841] or modifying the native pathways [39, 4244].

Under oxygen-limiting conditions the strains carrying the PHB pathway (TMB4425 and TMB4445) were found to have higher glycerol yields and lower ethanol yields than the corresponding reference strains (TMB4424 and TMB4444). This can be expected when a part of acetaldehyde, i.e. the precursor of that contributes to the re-oxidation of NADH from glycolysis, is redirected towards acetate formation for PHB production. Indeed we observed that the difference in ethanol molar yield corresponded to the increase in glycerol molar yield (Additional file 3). This would mean that PHB synthesis cannot per se act as an alternative redox sink to ethanol, since its synthesis only regenerates 0.5 NAD+ per pyruvate, i.e. half of the NAD+ regenerated via ethanol formation. Therefore glycerol formation is preferred. We also hypothesize that the redirected flux from ethanol was used for PHB production; however the produced PHB level per g xylose was much lower than the potentially redirected acetyl-CoA precursor (Additional file 3); so a part of the redirected flux might be used for other purposes, which could for instance also explain the increase in biomass yield in the PHB-producing strains as compared to their controls.


This work demonstrated the advantages of using a NADH-dependent AAR for PHB synthesis from xylose in S. cerevisiae compared to the regularly used NADPH-dependent AAR of C. necator. The NADH-dependent AAR from A. vinosum not only improved the conversion of xylose into PHB under aerobic conditions but it also enabled anaerobic PHB production from xylose when combined with redox engineered XR enzyme in recombinant S. cerevisiae. Further optimisation for higher titers, volumetric yields and productivities will notably require to engineer the central carbon metabolism in order to boost the production of the acetyl-CoA pathway precursor.



acetyl coenzyme A


acetyl-coA reductase


nicotinamide adenine dinucleotide


nicotinamide adenine dinucleotide phosphate


polymerase chain reaction


xylitol dehydrogenase

XRmut :

mutated xylose reductase

XRwt :

wild type xylose reductase

YsX :

yield of biomass on xylose

YsEtOH :

yield of ethanol on xylose

YsAc :

yield of acetate on xylose

YsXylitol :

yield of xylitol on xylose

YsGly :

yield of glycerol on xylose


yield of PHB on xylose


  1. Plastics—the facts 2015. An analysis of European plastics production, demand and waste data.—the-facts-2015.aspx.

  2. Griffin GJL. Chemistry and technology of biodegradable polymers. J Chem Educ. 1995;72:A73.

    CAS  Google Scholar 

  3. Madison LL, Huisman GW. Metabolic engineering of poly(3-hydroxyalkanoates): from DNA to plastic. Microbiol Mol Biol Rev. 1999;63:21–53.

    CAS  Google Scholar 

  4. Brigham CJ, Reimer EN, Rha CK, Sinskey AJ. Examination of PHB depolymerases in Ralstonia eutropha: further elucidation of the roles of enzymes in PHB homeostasis. AMB Express. 2012;2:1.

    Article  Google Scholar 

  5. Suriyamongkol P, Weselake R, Narine S, Moloney M, Shah S. Biotechnological approaches for the production of polyhydroxyalkanoates in microorganisms and plants—a review. Biotechnol Adv. 2007;25:148–75.

    Article  CAS  Google Scholar 

  6. Liebergesell M, Steinbüchel A. Cloning and nucleotide sequences of genes relevant for biosynthesis of poly(3-hydroxybutyric acid) in Chromatium vinosum strain D. Eur J Biochem. 1992;209:135–50.

    Article  CAS  Google Scholar 

  7. Amos DA, McInerney MJ. Formation of d-3-hydroxybutyryl-coenzyme A by an acetoacetyl-coenzyme A reductase in Syntrophomonas wolfei subsp. wolfei. Arch Microbiol. 1993;159:16–20.

    Article  CAS  Google Scholar 

  8. Perez-Zabaleta M, Sjöberg G, Guevara-Martínez M, Jarmander J, Gustavsson M, Quillaguamán J, Larsson G. Increasing the production of (R)-3-hydroxybutyrate in recombinant Escherichia coli by improved cofactor supply. Microb Cell Fact. 2016;15:1–10.

    Article  Google Scholar 

  9. Larsson S, Palmqvist E, Hahn-Hägerdal B, Tengborg C, Stenberg K, Zacchi G, Nilvebrant N-O. The generation of fermentation inhibitors during dilute acid hydrolysis of softwood. Enzym Microbial Technol. 1999;24:151–9.

    Article  CAS  Google Scholar 

  10. Lau MW, Gunawan C, Balan V, Dale BE. Comparing the fermentation performance of Escherichia coli KO11, Saccharomyces cerevisiae 424A (LNH-ST) and Zymomonas mobilis AX101 for cellulosic ethanol production. Biotechnol Biofuels. 2010;3:1.

    Article  Google Scholar 

  11. Hahn-Hägerdal B, Pamment N. Special session A microbial pentose metabolism. Appl Biochem Biotechnol. 2004;116:1207–9.

    Article  Google Scholar 

  12. Leaf T, Peterson M, Stoup S, Somers D, Srienc F. Saccharomyces cerevisiae expressing bacterial polyhydroxybutyrate synthase produces poly-3-hydroxybutyrate. Microbiology. 1996;142(Pt 5):1169–80.

    Article  CAS  Google Scholar 

  13. Carlson R, Srienc F. Effects of recombinant precursor pathway variations on poly[(R)-3-hydroxybutyrate] synthesis in Saccharomyces cerevisiae. J Biotechnol. 2006;124:561–73.

    Article  CAS  Google Scholar 

  14. Breuer U, Terentiev Y, Kunze G, Babel W. Yeasts as producers of polyhydroxyalkanoates: genetic engineering of Saccharomyces cerevisiae. Macromol Biosci. 2002;2:380–6.

    Article  CAS  Google Scholar 

  15. Kocharin K, Chen Y, Siewers V, Nielsen J. Engineering of acetyl-CoA metabolism for the improved production of polyhydroxybutyrate in Saccharomyces cerevisiae. AMB Express. 2012;2:52.

    Article  Google Scholar 

  16. Sandström A, Muñoz de las Heras A, Portugal-Nunes D, Gorwa-Grauslund MF. Engineering of Saccharomyces cerevisiae for the production of poly-3-d-hydroxybutyrate from xylose. AMB Express. 2015;5:14.

    Article  Google Scholar 

  17. Gietz R, Schiestl R. Quick and easy yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007;2:35–7.

    Article  CAS  Google Scholar 

  18. Sambrook J, Russell D. The inoue method for preparation and transformation of competent E. coli: “ultra-competent” cells. CSH Protoc. 2006;2:3944.

    Google Scholar 

  19. Karhumaa K, Hahn-Hägerdal B, Gorwa-Grauslund M. Investigation of limiting metabolic steps in the utilization of xylose by recombinant Saccharomyces cerevisiae using metabolic engineering. Yeast. 2005;22:359–68.

    Article  CAS  Google Scholar 

  20. Runquist D, Hahn-Hägerdal B, Bettiga M. Increased ethanol productivity in xylose-utilizing Saccharomyces cerevisiae via a randomly mutagenized xylose reductase. Appl Environ Microbiol. 2010;76:7796–802.

    Article  CAS  Google Scholar 

  21. Law JH, Slepecky RA. Assay of poly-β-hydroxybutyric acid. J Bacteriol. 1961;82:33–6.

    CAS  Google Scholar 

  22. Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T. Protein structure homology modeling using SWISS-MODEL workspace. Nat Protoc. 2009;4:1–13.

    Article  CAS  Google Scholar 

  23. Matsumoto K, Tanaka Y, Watanabe T, Motohashi R, Ikeda K, Tobitani K, Yao M, Tanaka I, Taguchi S. Directed evolution and structural analysis of NADPH-dependent acetoacetyl-CoA reductase from Ralstonia eutropha reveals two mutations responsible for enhanced kinetics. Appl Environ Microbiol. 2013;79:6134–9.

    Article  CAS  Google Scholar 

  24. Kim J, Chang JH. Kim E-JJ, Kim K-JJ: Crystal structure of (R)-3-hydroxybutyryl-CoA dehydrogenase PhaB from Ralstonia eutropha. Biochem Biophys Res Comm. 2014;443:783–8.

    Article  CAS  Google Scholar 

  25. Lierbergesell M, Steinbüchel A. Cloning and nucleotide sequences of genes relevant for biosynthesis of poly(3-hydroxybutyric acid) in Chromatium vinosum strain D. Eur J Biochem. 1992;209:135–50.

    Article  Google Scholar 

  26. Sieber JR, Sims DR, Han C, Kim E, Lykidis A, Lapidus AL, McDonnald E, Rohlin L, Culley DE, Gunsalus R, McInerney MJ. The genome of Syntrophomonas wolfei: new insights into syntrophic metabolism and biohydrogen production. Environ Microbiol. 2010;12:2289–301.

    CAS  Google Scholar 

  27. Lukacs CM, FJ, Edwards TE, Lorimer D. Crystal structure of a putative glucose 1-dehydrogenase from Burkholderia cenocepacia J2315. 2013. Accessed 15 Nov 2016.

  28. Chen C, Zhuang NN, Lee KH. 3-oxoacyl-[acyl-carrier-protein] reductase from Synechococcus elongatus PCC 7942; 2012.

  29. Suhadolnik RJ, Lennon MB, Uematsu T. Role of adenine ring and adenine ribose of nicotinamide adenine dinucleotide in binding and catalysis with alcohol, lactate, and glyceraldehyde-3-phospate dehydrogenases. J Biol Chem. 1977;2252:4125–33.

    Google Scholar 

  30. Bengtsson O, Hahn-Hägerdal B, Gorwa-Grauslund MF. Xylose reductase from Pichia stipitis with altered coenzyme preference improves ethanolic xylose fermentation by recombinant Saccharomyces cerevisiae. Biotechnol Biofuels. 2009;2:1–10.

    Article  Google Scholar 

  31. Bruinenberg PM, Bot PHM, Dijken JP, Scheffers AW. The role of redox balances in the anaerobic fermentation of xylose by yeasts. Eur J Appl Microbiol Biotechnol. 1983;18:287–92.

    Article  CAS  Google Scholar 

  32. Klimacek M, Kirl E, Krahulec S, Longus K, Novy V, Nidetzky B. Stepwise metabolic adaption from pure metabolization to balanced anaerobic growth on xylose explored for recombinant Saccharomyces cerevisiae. Microb Cell Fact. 2014;13:1–12.

    Article  Google Scholar 

  33. Cadete RM, de las Heras AM, Sandström AG, Ferreira C, Gírio F, Gorwa-Grauslund M-F, Rosa CA, Fonseca C. Exploring xylose metabolism in Spathaspora species: XYL1.2 from Spathaspora passalidarum as the key for efficient anaerobic xylose fermentation in metabolic engineered Saccharomyces cerevisiae. Biotechnol Biofuels. 2016;9:167.

    Article  Google Scholar 

  34. Bakker BM, Overkamp KM, van Maris AJ, Kötter P, Luttik MA, van Dijken JP, Pronk JT. Stoichiometry and compartmentation of NADH metabolism in Saccharomyces cerevisiae. FEMS Microbiol Rev. 2001;25:15–37.

    Article  CAS  Google Scholar 

  35. Weusthuis RA, Lamot I, Oost J, Sanders JPM. Microbial production of bulk chemicals: development of anaerobic processes. Trends Biotechnol. 2011;29:153–8.

    Article  CAS  Google Scholar 

  36. Reinecke F, Steinbüchel A. Ralstonia eutropha strain H16 as model organism for PHA metabolism and for biotechnological production of technically interesting biopolymers. J Mol Microbiol Biotechnol. 2008;16:91–108.

    Article  Google Scholar 

  37. Nissen TL, Anderlund M, Nielsen J, Villadsen J, Kielland-Brandt MC. Expression of a cytoplasmic transhydrogenase in Saccharomyces cerevisiae results in formation of 2-oxoglutarate due to depletion of the NADPH pool. Yeast. 2001;18:19–32.

    Article  CAS  Google Scholar 

  38. Medina VG, Almering MJH. Elimination of glycerol production in anaerobic cultures of a Saccharomyces cerevisiae strain engineered to use acetic acid as an electron acceptor. Appl Environ Microbiol. 2010;76:190–5.

    Article  CAS  Google Scholar 

  39. Shiba Y, Paradise EM, Kirby J, Ro D-K, 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 

  40. Kocharin K, Siewers V, Nielsen J. Improved polyhydroxybutyrate production by Saccharomyces cerevisiae through the use of the phosphoketolase pathway. Biotechnol Bioeng. 2013;110:2216–24.

    Article  CAS  Google Scholar 

  41. Kozak BU, van Rossum HM, Luttik MAH, Akeroyd M, Benjamin KR, Wu L, de Vries S, Daran J-M, 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-14.

    Article  Google Scholar 

  42. Tang X, Feng H, Chen WN. Metabolic engineering for enhanced fatty acids synthesis in Saccharomyces cerevisiae. Metab Eng. 2013;16C:95–102.

    Article  Google Scholar 

  43. 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 

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

    Article  Google Scholar 

  45. Gietz DR, Akio S. New yeast-Escherichia coli shuttle vectors constructed with in vitro mutagenized yeast genes lacking six-base pair restriction sites. Gene. 1988;74:527–34.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

AMH drafted the manuscript, carried out molecular genetic engineering, enzymatic assay, growth assays and HPLC metabolite analysis. DPN carried out bioreactor experiments and HPLC metabolite analysis and helped drafting the manuscript. NR carried out strain construction, growth assays and enzymatic assays. AS participated in the design of the study and its coordination, carried out molecular genetic engineering, structural biology analysis, enzymatic assays and helped drafting the manuscript. MFGG conceived the study, participated in the design of the study and its coordination and helped drafting and revising the manuscript. All authors read and approved the final manuscript.


We would like to thank Lisa Wasserstrom and Magnus Carlquist for their input in the experimental work and in the revision of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets analysed during the current study are available from the corresponding author on reasonable request.


EU Commission (FP7 project BRIGIT, Contract No. 311935) is acknowledged for their financial support.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Marie F. Gorwa-Grauslund.

Additional files


Additional file 1. Schematic representation of genes and regulatory sequences in the integrative vector YIpAGS3.


Additional file 2. The values reported are calculated from cultivations on defined media with xylose (50 g/L) as carbon source. The specific growth rate is calculated from the exponential phase. Yields and titers are calculated from a single time point in the end of the cultivation. Reported values represent the mean ± SD of at least two independent cultivations performed. YsX: yield of biomass on xylose, YsEtOH: yield of ethanol on xylose, YsAc: yield of acetate on xylose, YsXylitol: yield of xylitol on xylose, YsGly: yield of glycerol on xylose, YsPHB: yield of PHB on xylose. PHB titer: the volumetric PHB titer. PHB/CDW: The PHB component as percentage of total cell dry weight.


Additional file 3. Theoretical analysis of flux redirection in the case of TMB4425 under oxygen-limited conditions.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( 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

de las Heras, A.M., Portugal-Nunes, D.J., Rizza, N. et al. Anaerobic poly-3-d-hydroxybutyrate production from xylose in recombinant Saccharomyces cerevisiae using a NADH-dependent acetoacetyl-CoA reductase. Microb Cell Fact 15, 197 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Saccharomyces cerevisiae
  • Poly-3-d-hydroxybutyrate (PHB)
  • Xylose
  • NADH
  • Acetoacetyl-CoA reductase