Substrates and oxygen dependent citric acid production by Yarrowia lipolytica: insights through transcriptome and fluxome analyses
- Wael Sabra†1View ORCID ID profile,
- Rajesh Reddy Bommareddy†1, 3,
- Garima Maheshwari1,
- Seraphim Papanikolaou2 and
- An-Ping Zeng1Email author
© The Author(s) 2017
Received: 21 December 2016
Accepted: 23 April 2017
Published: 8 May 2017
Unlike the well-studied backer yeast where catabolite repression represents a burden for mixed substrate fermentation, Yarrowia lipolytica, an oleaginous yeast, is recognized for its potential to produce single cell oils and citric acid from different feedstocks. These versatilities of Y. lipolytica with regards to substrate utilization make it an attractive host for biorefinery application. However, to develop a commercial process for the production of citric acid by Y. lipolytica, it is necessary to better understand the primary metabolism and its regulation, especially for growth on mixed substrate.
Controlling the dissolved oxygen concentration (pO2) in Y. lipolytica cultures enhanced citric acid production significantly in cultures grown on glucose in mono- or dual substrate fermentations, whereas with glycerol as mono-substrate no significant effect of pO2 was found on citrate production. Growth on mixed substrate with glucose and glycerol revealed a relative preference of glycerol utilization by Y. lipolytica. Under optimized conditions with pO2 control, the citric acid titer on glucose in mono- or in dual substrate cultures was 55 and 50 g/L (with productivity of 0.6 g/L*h in both cultures), respectively, compared to a maximum of 18 g/L (0.2 g/L*h) with glycerol in monosubstrate culture. Additionally, in dual substrate fermentation, glycerol limitation was found to trigger citrate consumption despite the presence of enough glucose in pO2-limited culture. The metabolic behavior of this yeast on different substrates was investigated at transcriptomic and 13C-based fluxomics levels.
Upregulation of most of the genes of the pentose phosphate pathway was found in cultures with highest citrate production with glucose in mono- or in dual substrate fermentation with pO2 control. The activation of the glyoxylate cycle in the oxygen limited cultures and the imbalance caused by glycerol limitation might be the reason for the re-consumption of citrate in dual substrate fermentations. This study provides interesting targets for metabolic engineering of this industrial yeast.
With increasing focus on the development of sustainable technologies, the necessity for converting alternative, cheaper and waste carbon sources is becoming obvious. A further and equally important consideration is the development of sustainable microbial platforms to produce a variety of products from a vast array of biomass. Raw glycerol, which has become available in large quantities, and glucose, the most abundant carbon source in nature are applied in different industrial scale bioprocesses. Indeed, and with respect to price fluctuation and availability, finding efficient cell factories for the simultaneous conversion of these substrates to value added products will offer greater flexibility. Yarrowia lipolytica, an oleaginous yeast, is one of the most extensively studied ‘‘non-conventional’’ yeasts. Being a strictly aerobic microorganism it is capable of producing important metabolites and accumulates large amount of lipids from variety of substrates [1–4]. Y. lipolytica is known for its ability to secrete organic acids including tricarboxylic acid (TCA) cycle intermediates or precursors such as succinic, acetic, citric and isocitric acids [5–8]. Availability of genome sequence and the existence of suitable tools for genetic manipulation have made it possible to use the metabolic function of this species for biotechnological applications [9, 10]. One of the most striking features of this yeast is the presence of several multigene families involved in these metabolic pathways. The complexity and multiplicity of these genes enable Y. lipolytica to use and valorize a wide range of raw substrates and wastes [11–16].
It is well documented that in various yeast strains including Y. lipolytica, nitrogen limitation gradually change the metabolism from cell growth phase to citric acid production and/or storage lipid accumulation phase . Intracellular citrate is known to be the prime carbon source for fatty acid synthesis, and can also regulates glucose metabolism via its allosteric inhibition of phosphofructokinase [7, 18, 19]. Citrate is a key intermediate in both catabolism and anabolism, and it occupies a prominent position in the yeast energy metabolism, and therefore, its production is affected by the cultivation conditions and the energy and C-source used. Despite an increased number of biotechnological applications performed by several wild or genetically engineered Y. lipolytica strains [13–15, 20–22] relatively few studies have focused on the metabolism, growth kinetics and product formation in mixed substrate fermentation [7, 23]. In general, dual substrate fermentation by oleaginous yeasts has been little investigated [24, 25]. Hence, the starting point of this work is to compare the growth of Y. lipolytica and products formation in dual substrate of glucose and glycerol versus mono-substrate fermentation. Since growth and citric acid production are all known to be affected by oxygen availability, mono- and dual substrate fermentations were carried out in cultures with or without control of dissolved oxygen tension (pO2). Here, we report that controlling the pO2 at 50% air saturation enhance citric acid production when glucose was present in mono- or in dual substrate fermentation, whereas, no significant effect of the pO2 on citric acid production on glycerol was observed. Moreover, utilization of citric acid despite the presence of excess glucose in the medium was observed after glycerol limitation in the pO2 uncontrolled cultures. In order to understand such metabolic behavior, RNA-seq based transcriptome and 13C-based metabolic flux analyses were performed.
Results and discussion
Glycerol as a favored substrate for growth in the dual substrate fermentation in shake flask cultures
Biomass (X), lipid accumulation and citric acid production (Cit) of Y. lipolytica grown in nitrogen limited medium on glucose (Glc), glycerol (Glol) or blend of both
Lipid (%, w/w)
Oxygen requirements for citric acid production depends on the nature of the carbon source
Interestingly, with glucose as the sole carbon source, four different phases can be seen. The first phase in which the cells grew in the lag and the early logarithmic phase with excess dissolved oxygen concentration. The second phase was characterized by increased oxygen consumption and CO2 production rate parallel to exponential growth and hence a pO2 limited phase. An initiation of citric acid production characterizes the third phase which was parallel to the abrupt increase in pO2 values, indicating a nutrient limitation, and finally the fourth phase which is a stationary growth phase with continued citrate production under pO2 limited conditions (Fig. 1). On the other hand, with glycerol as the sole carbon source, citric acid production was initiated at the onset of oxygen limitation. In both mono-substrate fermentations with the individual substrate, the citrate production was initiated at the mid logarithmic phase and continued with the onset of the stationary phase (Fig. 1). Previously, Morgunov et al.  reported that after nitrogen limitation and at the onset of stationary phase in Y. lipolytica, the activity of citrate synthase remained high, whereas the activities of citric acid degrading enzymes such as aconitate hydratase and NAD-dependant isocitrate dehydrogenase were considerably decreased. Since oxygen availability will affect most of the TCA cycle enzymes, the aim of the following part was to compare the growth and citric acid production using mono-substrate fermentation in pO2 controlled bioreactors.
The effect of controlling the dissolved oxygen tension (pO2) on biomass, lipid and citric acid production by Y. lipolytica grown in nitrogen limited medium with glucose, glycerol or blend of both
Citric acid yield in g/g substrate
(Maximum) in stationary phase
Blend with BH
Blend with BH
During growth of Y. lipolytica strains under nitrogen-limited condition, three types of lipid metabolism have been presented in literatures. The typical “oleaginous” metabolism; in which lipid accumulation is triggered after nitrogen limitation in significant quantities, while low or negligible quantities of extra-cellular metabolites are simultaneously produced (mostly citric acid and to lesser extent polyols ). An atypical “oleaginous” metabolism; in which lipid is initially stored after nitrogen limitation, after which lipid accumulation decreases while citric acid production occurs uninterruptedly . A third atypical “oleaginous” metabolism was also described in which lipid accumulation starts at slow rates inside the yeast cells without degradation and with citric acid being constantly secreted extracellularly . By taking into consideration the results achieved with the strain used in the current investigation, it can be considered that ACA-DC 50109 follows the second or the third category of metabolism.
In fact, the most significant changes were recorded in the citrate production by Y. lipolytica on both substrates. Oxygen stimulated citrate production only if glucose was used as the C-source (average citric acid titers were 55 and 16 g/L at pO2 controlled and uncontrolled, respectively) whereas no significant change of citrate production was detected on glycerol (18 and 23 g/L at pO2 controlled and uncontrolled, respectively, Table 2). Previously, Kamzolova et al.  reported that oxygen requirements for growth and citric acid synthesis were found to depend on the iron concentration of the medium. Only at lower iron concentration (0.7 mg/L) higher oxygenation was reported to enhance citric acid production in Y. lipolytica with glycerol . This in fact can explain the results obtained with glycerol, where no significant differences in citric acid production was obtained with and without pO2 control in our iron rich medium (150 mg/L) (Fig. 1). In general, the difference in the nature of metabolism of both glycerol, an energy-poor carbon source and glucose and their effect on acid production are still not understood and deserve more in depth studies.
Higher fluxes to pentose phosphate pathway (PPP) and lower to TCA cycle were recorded with cells grown on glucose compared to glycerol at pO2-controlled cultures
Citrate re-consumption after glycerol limitation in the oxygen limited dual substrate culture
As shown in Fig. 4, after glycerol limitation, and despite the presence of excess glucose, utilization of citrate was observed in the pO2 uncontrolled culture. In the same culture, growth arrest was also observed after glycerol limitation (Fig. 4b). The same trend was also noticed using blend with spruce biomass hydrolysate in pO2 uncontrolled batch cultures (Additional file 1: Figure S1). Oxygen limitation in mono-substrate cultivation on either glucose or glycerol as shown in Fig. 1 didn’t trigger citrate re-consumption (with the current strain, or with Y. lipolytica A101 as reported recently ). This physiological behavior in dual-substrate fermentation is reported here for the first time. Y. lipolytica was reported to degrade citric acid after C-source limitation in mono-substrate fermentation , but was never investigated in dual substrate fermentation and especially when a carbon source like glucose was in excess. On the other hand, at controlled pO2 conditions, the yeast continued to produce citric acid as the major product, even after the assimilation of the favorable carbon source (glycerol). The citric acid yield was very similar to the mono-substrate fermentation in presence of glucose at pO2 controlled conditions (0.6 g/g substrate, Table 1). To understand the behavior difference in Y. lipolytica metabolism in dual substrate fermentation at oxygen limited and excess conditions, detailed transcriptomic analysis were done.
General stress metabolism after glycerol limitation in dual substrate fermentation with Y.lipolytica
Yeast cells are known to initiate a common gene expression program that generally protects the cell under stress . This response, referred to as the environmental stress response, includes ~900 genes whose expression is stereotypically altered when yeast cells are shifted to stressful environments . Irrespective of the pO2 control, several stress genes were found to be highly regulated after glycerol consumption. Interestingly, glycerol assimilating enzymes, namely glycerol kinase encoded by GUT1 (YALI0_F00484 g) and glycerol 3P dehydrogenase encoded by GUT2 (YALI0_B13970 g), were not significantly expressed (NSE) after glycerol limitation at both controlled and uncontrolled pO2 conditions (<2-fold change). Stress responses initiated by nitrogen limitation can be neglected as all samples for transcriptomic analysis were collected after nitrogen limitation. Examples are genes normally upregulated in response to osmotic stress (YALI0_F01210 g, 6.6-fold and YALI0_C00759 g, 2.3-fold at pO2-uncontrolled; YALI0_F01210 g 2.9 at pO2-controlled); starvation (YALI0_E31757 g, 5.2-fold at pO2-uncontrolled; YALI0_F14377 g, 2.9, YALI0_E20163 g, 2.6-fold at pO2-controlled), cell aging (YALI0_F31977 g, 2.4-fold at pO2-uncontrolled), oxidative stress genes (YALI0_E34265 g, 2.3-fold at pO2-uncontrolled), or chaperonin and other cellular stimuli (YALI0_C04884 g, 4.3-fold; YALI0_F25905 g, 2.4-fold at pO2-uncontrolled, YALI0_F00880 g, 2.4-fold at pO2-controlled) (Additional file 1: Table 1).
Growth arrest after glycerol limitation and irrespective to the pO2 control was also seen form the down-regulation of several growth related genes. Examples are some cell division control proteins that were downregulated in both cultures after glycerol limitation (CDC15:YALI0_F08165 g, −2.9 and −3.4-fold; CDC6:YALI0C00671 g, −4.6-fold and NSE; cell division cycle 20: YALI0_C03377 g, −2.1 and −2.4-fold; cell division cycle 14: YALI0_E16038 g, NSE and -2.4-fold at pO2-controlled and uncontrolled conditions, respectively, Additional file 1: Table 1). Growth cessation triggered by the limitation of the favorable substrate was previously investigated in Clostridium pasteurianum in a blend of glucose and glycerol .
Moreover, a large number of transport proteins, including multiple oligopeptide transporters were highly induced after glycerol limitation especially with pO2-uncontrolled conditions (YALI0_A00110 g, 11.9-fold and YALI0_D23485 g, 8.4-fold; YALI0_E19294 g, 6.9-fold; YALI0_E04510 g, 4.7-fold; YALI0_F30987 g, 4.3-fold; YALI0_B02398 g, 4.2-fold; YALI0_C18491 g, 4.0-fold; YALI0_C22616 g, 3.2-fold; at pO2-uncontrolled compared to, YALI0_D23485 g, 3.1-fold; YALI0_F18964 g, −5.0-fold; YALI0_C18491 g, −3.1-fold and YALI0_D19558 g, −2.4-fold; at pO2-controlled conditions). This may suggest that Y. lipolytica may use small peptides as nutrients, and such a behavior could support the growth of this yeast in natural environments where Yarrowia is naturally found (cheeses, yoghurts and meat). Similar behavior was recorded in Candida albicans when cells shifted from glucose, the preferred C-source to glycerol containing medium .
Allantoate permease was also highly regulated after glycerol limitation (relatively more at pO2-uncontrolled conditions, YALI0_A03091 g, 11.4-fold; YALI0_F28193 g, 10.8-fold; YALI0_F16423 g, 7.2-fold; YALI0_D13970 g, 4.5-fold; YALI0_C08569 g, 4.4-fold; YALI0_A00869 g, 3.8-fold; YALI0_C13574 g, 3.8-fold; YALI0_D05819 g, 2.7-fold; YALI0_C16907 g, 2.2-fold, compared to YALI0_A03091 g, 4.7-fold; YALI0_D13970 g, 4.2-fold; YALI0D13970 g, 3.8-fold; YALI0_C08569 g, 3.2-fold, YALI0_D13970 g, 2.4-fold at pO2-controlled experiment). Similarly, the gene YALI0_C07590 g for organic acid: sodium symporter activity was highly regulated after glycerol limitation especially in pO2-uncontrolled conditions (6.3-fold compared to 2.9 at pO2 controlled conditions). All these genes are related to RTG (retrograde signaling pathways) which are activated during a deficient mitochondrial function or in response to deficiencies in respiration.
Upregulation of the glyoxylate cycle after glycerol limitation in the pO2-uncontrolled culture
Important genes associated with citric acid and central carbon metabolism in Y. lipolytica were studied in culture without pO2 control before and after glycerol limitation (Fig. 5). Although the glyoxylate and TCA cycles share common reactions, it is only the genes for the glyoxylate cycle that were upregulated after glycerol consumption in the pO2 un-controlled cultures (malate synthase: YALI0_D19140 g, 4.5-fold and YALI0_E15708 g, 5.4-fold, Fig. 5). Moreover, enzymes for precursors synthesis needed for gyloxylate cycle were found to be up-regulated. In addition to the glycolate oxidase for the synthesis of glyoxylate (YALI0_D12661 g, 2.2-fold), several genes for the replenishment of acetyl CoA were found to be up-regulated after glycerol consumption. Examples are several genes involved in fatty acid degradation (YALI0_C23859 g, 3.9-fold; YALI0_D24750 g, 5.6-fold; YALI0_F10857 g, 4.5-fold; YALI0_B10406 g, 3.5; YALI0_E11099 g, 3.8-fold), and amino acid degradation (especially leucine, valine and isoleucine, YALI0_D08690 g, 6.6-fold; YALI0_D23815 g, 5.2-fold, YALI0_B10406 g, 3.5-fold; YALI0_B22550 g, 2.3-fold; YALI0_E11099 g, 3.8-fold). Also several genes involved in peroxisome biogenesis, a compartment of central role in glyoxylate cycle function, were found to be up-regulated in the glycerol limited phase (YALI0_B22660 g, 4.2-fold; YALI0_A20944 g, 3.5-fold; YALI0C18689 g, 2.8-fold; YALI0_F28457 g, 2.8-fold; YALI0_E09405 g, 2.8-fold; YALI0_C05775 g, 2.2-fold; YALI0D26642 g, 2.1; YALI0_C09504 g, 2.1-fold). In control experiments without the metabolic shift with pO2 controlled condition, neither genes expression of glyoxylate cycle nor peroxisome biogenesis were changed significantly in response to glycerol consumption (Fig. 4a, Additional file 1: Figure S4).
The glyoxylate cycle is known to serve as a link between catabolic activities and biosynthetic capacities and enables cells to utilize fatty acids or C2-units as sole carbon source . Unlike S. cerevisiae, Y. lipolytica prefers glycerol as C-source and quantitative studies on dual substrate fermentation with glucose are still very few [7, 23]. Therefore, known metabolic models for the traditionally applied yeast cell factories can hardly be used to explain the metabolic shift observed in Y. lipolytica in response to glycerol and oxygen limitation (Figs. 4b, 5). Definitely, oxygen limitation certainly affects the main function of mitochondria and energy generation by oxidative phosphorylation, and as shown in transcriptomic data (Additional file 1: Table 1), the four isozymes of the NADH dehydrogenases were upregulated after glycerol limitation at the pO2-controlled condition (YALI0_C17853 g, 4.3-fold and YALI0_C22319 g, 2.6, YALI0_C22319 g, 2.59-fold at pO2 controlled conditions) while all of these genes were found to be non-significantly expressed (NSE) under oxygen limited conditions. It seems likely that, under reduced function of mitochondria as a result of oxygen limitation, citrate which is produced in the mitochondria serves predominantly as a key substrate for energy generation (Fig. 5). Citrate exhaustion in the cytoplasm may ultimately trigger its re-consumption from the medium through organic acid transporters. Recently, Guo et al.  identified 6 putative genes responsible for the transport of organic acids in Y. lipolytica, from which only 5 genes encode proteins that can transport citrate. Form these 5 genes, only 3 genes were differentially expressed: YALI0_B19470 g (−4.5 and 2.3-folds), YLI0_C15488 g (3.3 and 4.3-folds), and YALI0E32901 g (NSE and 7.73-folds at pO2 controlled and uncontrolled conditions; respectively). Still it is not known whether efflux and influx of citrate are shared by the same proteins.
Upregulation of the non-oxidative PPP is correlated to high citric acid production in Y. lipolytica
Higher citric acid production by Y. lipolytica was found to be dependent on the carbon source and oxygen availability (Figs. 2, 3). The high citrate production with glucose in monosubstrate cultivation (Fig. 2) coincide with the relatively higher fluxes toward PPP than with glycerol as sole C-source (Fig. 3). Interestingly, when the glycerol limited phases at both pO2 controlled and uncontrolled conditions were compared, up-regulation of almost all the genes of the non-oxidative PPP were observed in the pO2 controlled culture compared to the pO2 uncontrolled one (Fig. 6). In fact, in comparison to glycerol and similar to the glucose as monosubstrates, the PPP fluxes were also significantly higher on dual substrate fermentation with pO2 control (Additional file 1: Figure S5). The role of PPP on the metabolism and citric acid production in Y. lipolytica is still not clear and deserve more studies. Whereas an oxidative PPP was shown to be essential for lipid production from glucose in Y. lipolytica , our results showed that the upregulation of the non-oxidative PPP genes could lead to high citrate production from glucose (Fig. 6). Essentially the glycolytic pathway is common to all yeast species and the carbon flux regulation is being done at the level of the pentose phosphate pathway. Whether glucose 6-phosphate enters glycolysis or the pentose phosphate pathway depends on the current needs of the cell and on the concentration of NADP+ in the cytosol. Previously, it was reported that the marked effect of the NADP+ level on the rate of the oxidative PPP ensures that NADPH generation is tightly coupled to its utilization in reductive biosynthesis such as lipid accumulation, whereas, the nonoxidative phase of the PPP is controlled primarily by the availability of substrates .
Yarrowia lipolytica is known for both efficient citrate excretion and high lipid productivity under stress conditions such as nitrogen limitation. Under N-limitation, the prevalence of oxygen limitation was shown to enhance lipid production while decrease citrate production on glucose significantly (Table 2 and ). Besides the complex interplay between the oxygen effect and the nature of substrate in regulating citrate production, the availability of the C-source and the efficiency of the glycolytic flux were recently found to have a significant impact on Y. lipolytica metabolism . Decreasing the glycolytic flux through glucose uptake rate in controlled fed batch culture was reported to decrease citrate production and enhanced lipid production by this yeast . In addition to the double limitations exerted by N- and O2 in Y. lipolytica cells grown in pO2 uncontrolled condition, the onset of an additional glycerol limitation shifted the metabolism significantly (Fig. 4b), and the resulted low glycolytic fluxes triggered a decreased citric acid production and consequently re-consumption of the medium. Indeed, the different in the fluxes toward PPP with the two substrates (Fig. 3) also played here a significant role. As Y. lipolytica consumes both these substrates simultaneously with preference to glycerol, it looks like that initially glucose is only used to generate the PP flux for NADPH, whereas glycerol mainly is driven towards lipid and lower glycolysis. Hence, at the moment of glycerol limitation, the cells fail to redirect glucose through the lower glycolysis and TCA cycle rather, they use citric acid as the preferred carbon source. For the development of Y. lipolytica as a microbial cell factory for citric acid production and other biotechnological applications , the differences between the fluxes with respect to the activity of the PPP, TCA cycle and glyoxylate cycle had to be optimized and understood especially with dual substrates.
Citric acid production by Y. lipolytica grown on either glucose, glycerol, or mixture of them was investigated in cultures with controlled or limited pO2. The different nature of the two substrates, especially their reduction degree was reflected by the different fluxes toward PPP and TCA cycle and finally citrate as evidenced by 13C analysis for mono-substrate fermentations. Relatively higher PPP and lower TCA cycle fluxes were recorded in Y. lipolytica cells grown under sufficient oxygen supply on glucose or on blend with glycerol but not in cells grown with glycerol as the sole C-source. Down-regulation of the TCA cycle fluxes triggered by oxygen limitation led to an increased flux to lipid formation, and less to citrate production. The TCA cycle fluxes were further decreased in oxygen limited—dual substrate fermentation after glycerol limitation. The imbalance exerted in the fluxes of PPP and TCA cycle as a result of glycerol limitation resulted in up-regulation of the glyoxylate pathway and the re-utilization of citrate despite the presence of excess glucose. This study highlights the importance of controlling the oxygen supply for citrate production by Y. lipolytica on mixed substrates.
Microorganism and media
Yarrowia lipolytica ACA DC 50109 was used in the present. The strain was maintained at −80 °C on potato dextrose medium with 20% glycerol (w/v). The medium for seed culture and batch fermentations are similar to that reported previously . For isotope labeled experiment, yeast extract was replaced by vitamins solution which contained (per liter): biotin 0.05 g; p-amino benzoic acid 0.2 g; nicotinic acid 1 g; Ca-pantothenate 1 g; 19 pyridoxine–HCl 1 g and thiamine–HCl 1 g. The medium was supplemented with 30 g/L glucose or glycerol and for mixed substrate cultivations 1:1 mixture of glucose and glycerol were added. For the labelled experiments, 100% 1-13C glucose, 99% atom (Sigma, Germany) or 20% U-13C6 glucose 99% atom (EURISO-TOP GmbH, Germany) or 20% U-13C3 glycerol (99% atom, Cortecnet, France) were used. For mixed substrate labelled experiments, 50% U-13C3 glycerol, 20% 1-13C glucose and 30% naturally labelled glucose were used.
Flask experiments were conducted in 250-mL non-baffled flasks, containing 50 mL of nitrogen limited growth medium. Flasks were shacked at 200 rpm in an orbital shaker (Lab-Line, USA) at 28 ± 1 °C. To minimize the drop of pH KH2PO4 and Na2HPO4 at initial concentrations of 12 g/L of each were added. The initial pH in these media was 6.4 ± 0.1, and effort was done to maintain the pH of the culture medium to values >5.0, by aseptically adding few drops of concentrated KOH solution (≈5 M). The initial carbon source concentrations used were 50 g/L in all flask experiments. Mono-substrate and dual substrate fermentation with different variations of glucose to glycerol concentration were also compared.
Batch fermentations were performed in Dasgip fermenters (Dasgip, Jülich, Germany). 1.5 L capacity vessels with an initial working volume of 500 mL were used. For 13C labelled experiments 200 mL capacity fermenter vessels with an initial working volume of 120 mL were used. One replicate with labeled glycerol (20% U-13C3) and two replicates with un-labeled glycerol were performed. In case of glucose as the substrate, one replicate with 100% 1-13C glucose and second replicate with 20% U-13C6 glucose and third replicate with unlabeled glucose were performed. In case of mixed substrates one set of experiment was performed with labeled substrates and the remaining two sets were performed with naturally labeled substrates. Agitation and temperature were maintained at 400 rpm and 30 °C respectively. pH was maintained at 6.0 by automatic addition of 5 M NaOH. Dissolved oxygen was maintained at 50% air saturation throughout the experiment by adjusting the oxygen concentration in the inlet gas. Aeration was set to 0.5 vvm. The off-gas from the fermenters was analyzed by the Dasgip off-gas analyzing system to determine the real time oxygen uptake and CO2 evolution rates.
13C labeled substrate experiments and metabolic flux estimation
For 13C-labeled experiments, cells from the precultures were washed prior to inoculation in order to reduce the effect of unlabeled biomass. The washed cells were used to inoculate the fermenters to an initial OD600 of 0.03. Media and initial substrate concentrations were similar to what was described before. Samples collected during the mid-exponential phase were treated and analyzed using GC–MS (Agilent GC 7890B coupled to a 5977 MSD) as described previously .
The compartmentalized metabolic network model of Y. lipolytica metabolism was constructed based on a previous model of R. toruloides . Briefly, the metabolic network is comprised of glycolysis, PP pathway, TCA cycle, fatty acid and biomass synthesis. Citrate from mitochondria is transported to the cytosol through the citrate-malate shuttle. Acetyl-CoA in the cytosol is assumed to be formed only through the ATP: citrate lyase (ACL) reaction. Oxaloacetate (OAA) in the cytosol was assumed to be formed by the anaplerotic reaction, pyruvate carboxylase and the ACL. A cytosolic malate dehydrogenase (NAD+ forming) converting OAA to malate in the cytosol was included. Malate transport from cytosol to mitochondrial membrane via a dicarboxylate carrier was included. The software toolbox INCA  was used for parameter estimation and evaluation in Matlab 2013a.
Extracellular metabolites and lipid analysis
Cells were harvested during the cultivations after centrifugation (5000 rpm, 10 min at 4 °C). Cell dry weight was determined gravimetrically after drying the harvested cells in an oven at 80 °C to a constant weight. Total lipid extraction with chloroform/methanol mixture (2:1 v/v) was performed as described in previous work . Fatty acid composition was done as described previously by Papanikolaou et al. . Lipid free biomass was calculated after subtracting the lipid content from the total biomass. Quantification of glucose, glycerol and organic acids was carried out using high-performance liquid chromatography (HPLC; Kontron Instruments, United Kingdom) with separation on an Aminex HPX-87H column at 60 °C with 0.005 M H2SO4 and detection via refractive index or by UV absorption at 210 nm. Ammonia concentration in the supernatant was determined by photometric measurements using a kit from Macherey–Nagel, Germany.
RNA extraction and sequencing
For whole transcriptome sequencing, samples from the mixed substrate (1:1 mixture of glucose and glycerol) cultivations were withdrawn before and after glycerol limitation. Approximately 3 × 108 cells were collected by centrifugation (5000 rpm, 2 °C for 5 min) and the biomass was immediately re-suspended in RNAlater ® Stabilization Solution (Life technologies, Ambion). The total RNA was extracted using the RiboPure™ RNA purification kit for yeast (Life technologies, Ambion). The quality and quantity of the RNA was assayed using a BioAnalyzer (Agilent, USA). About 2 ng of total RNA was used for sequencing based on Illumina HiSeq 2000 platform. Library construction was performed using the Illumina Trueseq kit at BGI Tech, Hong Kong.
RNA-seq data analysis
The generated pair-end reads (FASTQ format) were checked for their quality using the FastQC tool. The files were groomed to meet the Sanger-scaled quality values with ASCII offset 33. The groomed paired-end reads were then mapped to the reference genome of Y. lipolytica CLIB122 GCA_000002525  using Tophat2 . The gene and transcript expression levels were estimated and normalized based on FPKM (fragments per kilobase of exon per million fragments mapped) which were estimated using Cufflinks  with standard parameters and a p value cut-off set to 0.05. Gene ontology (GO) analysis was performed by using Blast2GO programs [54–56]. All raw sequence reads were deposited into the NCBI Sequence Read under the accession number of GSE83109.
WS and RRB designed and performed the bioreactor experiments with GM. WS interpreted the transcriptome results and wrote the manuscript. RRB carried out RNA-seq and 13C flux analysis, interpreted the results and drafted the manuscript. SP designed and performed the flask experiment, drafted the manuscript. APZ supervised the research. All authors read and approved the final manuscript.
This publication was supported by the German Research Foundation (DFG) and the Hamburg University of Technology (TUHH) in the funding programme ‘Open Access Publishing’.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this article and additional file. The transcriptome datasets used and/or analysed during the current study available from the corresponding author on request.
Ethics approval and consent to participate
This manuscript dose not report data collected from humans or animals.
This work was funded by from the BMBF under the project ‘Bio4Oil project number: 03SF0467’. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis 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.
- Beopoulos A, Mrozova Z, Thevenieau F, Le Dall MT, Hapala I, Papanikolaou S, Chardot T, Nicaud JM. Control of lipid accumulation in the yeast Yarrowia lipolytica. Appl Environ Microbiol. 2008;74:7779–89.View ArticleGoogle Scholar
- Papanikolaou S, Chevalot I, Komaitis M, Aggelis G, Marc I. Kinetic profile of the cellular lipid composition in an oleaginous Yarrowia lipolytica capable of producing a cocoa-butter substitute from industrial fats. Antonie Van Leeuwenhoek. 2001;80:215–24.View ArticleGoogle Scholar
- Ledesma-Amaro R, Nicaud JM. Yarrowia lipolytica as a biotechnological chassis to produce usual and unusual fatty acids. Prog Lipid Res. 2016;61:40–50.View ArticleGoogle Scholar
- Ledesma-Amaro R, Dulermo T, Nicaud JM. Engineering Yarrowia lipolytica to produce biodiesel from raw starch. Biotechnol Biofuels. 2015;8:148.View ArticleGoogle Scholar
- Kamzolova SV, Shishkanova NV, Morgunov IG, Finogenova TV. Oxygen requirements for growth and citric acid production of Yarrowia lipolytica. FEMS Yeast Res. 2003;3:217–22.View ArticleGoogle Scholar
- Kamzolova SV, Vinokurova NG, Shemshura ON, Bekmakhanova NE, Lunina JN, Samoilenko VA, Morgunov IG. The production of succinic acid by yeast Yarrowia lipolytica through a two-step process. Appl Microbiol Biotechnol. 2014;98:7959–69.View ArticleGoogle Scholar
- Papanikolaou S, Muniglia L, Chevalot I, Aggelis G, Marc I. Yarrowia lipolytica as a potential producer of citric acid from raw glycerol. J Appl Microbiol. 2002;92:737–44.View ArticleGoogle Scholar
- Chatzifragkou A, Makri A, Belka A, Bellou S, Mavrou M, Mastoridou M, Mystrioti P, Onjaro G, Aggelis G, Papanikolaou S. Biotechnological conversions of biodiesel derived waste glycerol by yeast and fungal species. Energy. 2011;36:1097–108.View ArticleGoogle Scholar
- Qiao K, Wasylenko TM, Zhou K, Xu P, Stephanopoulos G. Lipid production in Yarrowia lipolytica is maximized by engineering cytosolic redox metabolism. Nat Biotechnol. 2017;35:173–7.View ArticleGoogle Scholar
- Xu P, Qiao K, Ahn WS, Stephanopoulos G. Engineering Yarrowia lipolytica as a platform for synthesis of drop-in transportation fuels and oleochemicals. Proc Natl Acad Sci USA. 2016;113:10848–53.View ArticleGoogle Scholar
- Papanikolaou S, Chevalot I, Komaitis M, Marc I, Aggelis G. Single cell oil production by Yarrowia lipolytica growing on an industrial derivative of animal fat in batch cultures. Appl Microbiol Biotechnol. 2002;58:308–12.View ArticleGoogle Scholar
- Papanikolaou S, Aggelis G. Biotechnological valorization of biodiesel derived glycerol waste through production of single cell oil and citric acid by Yarrowia lipolytica. Lipid Technol. 2009;21:83–7.View ArticleGoogle Scholar
- Papanikolaou S, Galiotou-Panayotou M, Chevalot I, Komaitis M, Marc I, Aggelis G. Influence of glucose and saturated free-fatty acid mixtures on citric acid and lipid production by Yarrowia lipolytica. Curr Microbiol. 2006;52:134–42.View ArticleGoogle Scholar
- Papanikolaou S, Galiotou-Panayotou M, Fakas S, Komaitis M, Aggelis G. Citric acid production by Yarrowia lipolytica cultivated on olive-mill wastewater-based media. Bioresour Technol. 2008;99:2419–28.View ArticleGoogle Scholar
- Papanikolaou S, Beopoulos A, Koletti A, Thevenieau F, Koutinas AA, Nicaud JM, Aggelis G. Importance of the methyl-citrate cycle on glycerol metabolism in the yeast Yarrowia lipolytica. J Biotechnol. 2013;168:303–14.View ArticleGoogle Scholar
- Mironczuk AM, Rzechonek DA, Biegalska A, Rakicka M, Dobrowolski A. A novel strain of Yarrowia lipolytica as a platform for value-added product synthesis from glycerol. Biotechnol Biofuels. 2016;9:180.View ArticleGoogle Scholar
- Sarris D, Galiotou-Panayotou M, Koutinas AA, Komaitis M, Papanikolaou S. Citric acid, biomass and cellular lipid production by Yarrowia lipolytica strains cultivated on olive mill wastewater-based media. J Chem Technol Biotechnol. 2011;86:1439–48.View ArticleGoogle Scholar
- Sun J, Aluvila S, Kotaria R, Mayor JA, Walters DE, Kaplan RS. Mitochondrial and plasma membrane citrate transporters: discovery of selective inhibitors and application to structure/function analysis. Mol Cell Pharmacol. 2010;2:101–10.Google Scholar
- Morgunov IG, Kamzolova SV, Lunina JN. The citric acid production from raw glycerol by Yarrowia lipolytica yeast and its regulation. Appl Microbiol Biotechnol. 2013;97:7387–97.View ArticleGoogle Scholar
- Bellou S, Triantaphyllidou IE, Mizerakis P, Aggelis G. High lipid accumulation in Yarrowia lipolytica cultivated under double limitation of nitrogen and magnesium. J Biotechnol. 2016;234:116–26.View ArticleGoogle Scholar
- Tomaszewska L, Rakicka M, Rymowicz W, Rywinska A. A comparative study on glycerol metabolism to erythritol and citric acid in Yarrowia lipolytica yeast cells. FEMS Yeast Res. 2014;14:966–76.View ArticleGoogle Scholar
- Rywinska A, Rymowicz W, Zarowska B, Skrzypinski A. Comparison of citric acid production from glycerol and glucose by different strains of Yarrowia lipolytica. World J Microbiol Biotechnol. 2010;26:1217–24.View ArticleGoogle Scholar
- Workman M, Holt P, Thykaer J. Comparing cellular performance of Yarrowia lipolytica during growth on glucose and glycerol in submerged cultivations. AMB Express. 2013;3:58.View ArticleGoogle Scholar
- Ledesma-Amaro R, Nicaud JM. Metabolic engineering for expanding the substrate range of Yarrowia lipolytica. Trends Biotechnol. 2016;34:798–809.View ArticleGoogle Scholar
- Lane S, Zhang S, Wei N, Rao C, Jin YS. Development and physiological characterization of cellobiose-consuming Yarrowia lipolytica. Biotechnol Bioeng. 2015;112:1012–22.View ArticleGoogle Scholar
- Rymowicz W, Fatykhova AR, Kamzolova SV, Rywinska A, Morgunov IG. Citric acid production from glycerol-containing waste of biodiesel industry by Yarrowia lipolytica in batch, repeated batch, and cell recycle regimes. Appl Microbiol Biotechnol. 2010;87:971–9.View ArticleGoogle Scholar
- Kamzolova SV, Vinokurova NG, Lunina JN, Zelenkova NF, Morgunov IG. Production of technical-grade sodium citrate from glycerol-containing biodiesel waste by Yarrowia lipolytica. Bioresour Technol. 2015;193:250–5.View ArticleGoogle Scholar
- Tchakouteu SS, Chatzifragkou A, Kalantzi O, Papanikolaou S. Oleaginous yeast Cryptococcus curvatus exhibits interplay between biosynthesis of intracellular sugars and lipids: biochemical features of Cryptococcus curvatus. Eur J Lipid Sci Technol. 2014;117:657–72.View ArticleGoogle Scholar
- Bommareddy RR, Sabra W, Zeng AP. Glucose-mediated regulation of glycerol uptake in Rhodosporidium toruloides: insights through transcriptomic analysis on dual substrate fermentation. Eng Life Sci. 2017;17(3):282–91.View ArticleGoogle Scholar
- Andre A, Chatzifragkou A, Diamantopoulou P, Sarris D, Philippoussis A, Galiotou-Panayotou M, Komaitis M, Papanikolaou S. Biotechnological conversions of bio-diesel-derived crude glycerol by Yarrowia lipolytica strains. Eng Life Sci. 2009;9:468–78.View ArticleGoogle Scholar
- Makri A, Fakas S, Aggelis G. Metabolic activities of biotechnological interest in Yarrowia lipolytica grown on glycerol in repeated batch cultures. Bioresour Technol. 2010;101:2351–8.View ArticleGoogle Scholar
- Tomaszewska L, Rywinska A, Gladkowski W. Production of erythritol and mannitol by Yarrowia lipolytica yeast in media containing glycerol. J Ind Microbiol Biotechnol. 2012;39:1333–43.View ArticleGoogle Scholar
- Tomaszewska L, Rymowicz W, Rywinska A. Mineral supplementation increases erythrose reductase activity in erythritol biosynthesis from glycerol by Yarrowia lipolytica. Appl Biochem Biotechnol. 2014;172:3069–78.View ArticleGoogle Scholar
- Papanikolaou S, Muniglia L, Chevalot I, Aggelis G, Marc I. Accumulation of a cocoa-butter-like lipid by Yarrowia lipolytica cultivated on agro-industrial residues. Curr Microbiol. 2003;46:124–30.View ArticleGoogle Scholar
- Bommareddy RR, Sabra W, Maheshwari G, Zeng AP. Metabolic network analysis and experimental study of lipid production in Rhodosporidium toruloides grown on single and mixed substrates. Microb Cell Fact. 2015;14:36.View ArticleGoogle Scholar
- Papanikolaou S, Aggelis G. Lipid of oleaginous yeast. Part I: biochemistry of single cell oil production. Eur J Lipid Sci Technol. 2011;113:1031–51.View ArticleGoogle Scholar
- Tsigie YA, Wang CY, Truong CT, Ju YH. Lipid production from Yarrowia lipolytica Po1 g grown in sugarcane bagasse hydrolysate. Bioresour Technol. 2011;102:9216–22.View ArticleGoogle Scholar
- Li X, Wang P, Ge Y, Wang W, Abbas A, Zhu G. NADP(+)-specific isocitrate dehydrogenase from oleaginous yeast Yarrowia lipolytica CLIB122: biochemical characterization and coenzyme sites evaluation. Appl Biochem Biotechnol. 2013;171:403–16.View ArticleGoogle Scholar
- Yeh JI, Chinte U, Du S. Structure of glycerol-3-phosphate dehydrogenase, an essential monotopic membrane enzyme involved in respiration and metabolism. Proc Natl Acad Sci USA. 2008;105:3280–5.View ArticleGoogle Scholar
- Gasch AP, Werner-Washburne M. The genomics of yeast responses to environmental stress and starvation. Funct Integr Genom. 2002;2:181–92.View ArticleGoogle Scholar
- Sabra W, Wang W, Surandram S, Groeger C, Zeng AP. Fermentation of mixed substrates by Clostridium pasteurianum and its physiological, metabolic and proteomic characterizations. Microb Cell Fact. 2016;15:114.View ArticleGoogle Scholar
- Lorenz MC, Bender JA, Fink GR. Transcriptional response of Candida albicans upon internalization by macrophages. Eukaryot Cell. 2004;3:1076–87.View ArticleGoogle Scholar
- Kunze M, Pracharoenwattana I, Smith SM, Hartig A. A central role for the peroxisomal membrane in glyoxylate cycle function. Biochim Biophys Acta. 2006;1763:1441–52.View ArticleGoogle Scholar
- Guo H, Liu P, Madzak C, Du G, Zhou J, Chen J. Identification and application of keto acids transporters in Yarrowia lipolytica. Sci Rep. 2015;5:8138.View ArticleGoogle Scholar
- Wasylenko TM, Ahn WS, Stephanopoulos G. The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica. Metab Eng. 2015;30:27–39.View ArticleGoogle Scholar
- Berg JM, Tymoczko JL, Stryer L. The metabolism of glucose 6-phosphate by the pentose phosphate pathway is coordinated with glycolysis. In: Biochemistry. 5th ed. New York: Freeman WH and Co Ltd; 2002.Google Scholar
- Kavscek M, Bhutada G, Madl T, Natter K. Optimization of lipid production with a genome-scale model of Yarrowia lipolytica. BMC Syst Biol. 2015;9:72.View ArticleGoogle Scholar
- Groenewald M, Boekhout T, Neuveglise C, Gaillardin C, van Dijck PW, Wyss M. Yarrowia lipolytica: safety assessment of an oleaginous yeast with a great industrial potential. Crit Rev Microbiol. 2014;40:187–206.View ArticleGoogle Scholar
- Zamboni N, Fendt SM, Ruhl M, Sauer U. (13)C-based metabolic flux analysis. Nat Protoc. 2009;4:878–92.View ArticleGoogle Scholar
- Young JD. INCA: a computational platform for isotopically non-stationary metabolic flux analysis. Bioinformatics. 2014;30:1333–5.View ArticleGoogle Scholar
- Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I, De Montigny J, Marck C, Neuveglise C, Talla E, et al. Genome evolution in yeasts. Nature. 2004;430:35–44.View ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.View ArticleGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25:25–9.View ArticleGoogle Scholar
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S. AmiGO: online access to ontology and annotation data. Bioinformatics. 2009;25:288–9.View ArticleGoogle Scholar