Skip to main content
  • Research article
  • Open access
  • Published:

Development of cellobiose-degrading ability in Yarrowia lipolytica strain by overexpression of endogenous genes



Yarrowia lipolytica, one of the most widely studied “nonconventional” oleaginous yeast species, is unable to grow on cellobiose. Engineering cellobiose-degrading ability into this yeast is a vital step towards the development of cellulolytic biocatalysts suitable for consolidated bioprocessing.


In the present work, we identified six genes encoding putative β-glucosidases in the Y. lipolytica genome. To study these, homologous expression was attempted in Y. lipolytica JMY1212 Zeta. Two strains overexpressing BGL1 (YALI0F16027g) and BGL2 (YALI0B14289g) produced β-glucosidase activity and were able to degrade cellobiose, while the other four did not display any detectable activity. The two active β-glucosidases, one of which was mainly cell-associated while the other was present in the extracellular medium, were purified and characterized. The two Bgls were most active at 40–45°C and pH 4.0–4.5, and exhibited hydrolytic activity on various β-glycoside substrates. Specifically, Bgl1 displayed 12.5-fold higher catalytic efficiency on cellobiose than Bgl2. Significantly, in experiments where cellobiose or cellulose (performed in the presence of a β-glucosidase-deficient commercial cellulase cocktail produced by Trichoderma reseei) was used as carbon source for aerobic cultivation, Y. lipolyticapox co-expressing BGL1 and BGL2 grew better than the Y. lipolytica strains expressing single BGLs. The specific growth rate and biomass yield of Y. lipolytica JMY1212 co-expressing BGL1 and BGL2 were 0.15 h−1 and 0.50 g-DCW/g-cellobiose, respectively, similar to that of the control grown on glucose.


We conclude that the bi-functional Y. lipolytica developed in the current study represents a vital step towards the creation of a cellulolytic yeast strain that can be used for lipid production from lignocellulosic biomass. When used in combination with commercial cellulolytic cocktails, this strain will no doubt reduce enzyme requirements and thus costs.


It is widely recognized that lignocellulosic biomass (or LC biomass) will form an important part of the future bio-economy. However, the use of this renewable resource as feedstock for industrial activities poses a major challenge, because its deconstruction to sugars and lignin is complex, requiring a series of unit operations. These include costly pretreatment and enzymatic hydrolysis steps, the latter requiring the action of several types of enzymes [1, 2]. Indeed, the hydrolysis of cellulose alone requires the synergistic action of endoglucanases (EC, cellobiohydrolases (EC and β-glucosidases (EC [3]. Endoglucanases are active on the internal bonds in cellulose and release free reducing and non-reducing extremities, which are used by cellobiohydrolases as starting points for exo-processive hydrolysis that yields cellodextrins as products. Finally, β-glucosidases convert cellodextrins into glucose [4].

One strategy to reduce investment and operational costs in LC biomass processing is to internalize enzyme production and combine enzymatic hydrolysis with fermentation. This is known as consolidated bioprocessing (CBP) and can be achieved using a microorganism that possesses the dual ability to produce biomass-hydrolyzing enzymes and ferment sugars to products of commercial interest, thus allowing a one-pot type bioconversion process in which process integration is maximized [5]. While CBP is considered to be an ultimate aim for biorefining, the ways to achieve this goal are not simple. Although the number of naturally occurring, biomass-degrading microorganisms is no doubt large, those that possess the ability to hydrolyze LC biomass and ferment free sugars into target products, such as ethanol, butanol, hydrogen, fatty acid ethyl esters (FAEE) or isopropanol, at industrially compatible rates and titers, are probably very rare and so far undiscovered [6]. Additionally, many of the best known biomass-degrading microorganisms display low β-glucosidase (cellobiase) activity, meaning that the hydrolysis of cellobiose constitutes a rate-limiting step during the enzymatic processing of cellulose [79]. Therefore, engineering cellobiose-degrading ability into microorganisms is a vital step towards the development of cellulolytic biocatalysts suitable for CBP. In this respect, examples of recent work performed on Saccharomyces cerevisiae, the current workhorse of biotechnological processes, are noteworthy [1012]. In these studies, even though the engineered S. cerevisiae strains exhibited poor cellulose-degrading ability, the fact that they both produce significant cellobiase activity means that their incorporation into a simultaneous saccharification and fermentation (SSF) process is likely to reduce the loading of external cellulases and thus overall process cost [10].

Although ethanol is the target molecule in many biorefinery concepts, Fatty Acid Esters (FAEs) such as those used in biodiesel, are also attractive targets. This is because FAEs display high energy density and are well-tolerated by production strains [13]. Currently, FAEs are mainly produced by transesterification of plant oils using an alcohol (methanol or ethanol) and base, acid or enzyme catalysts [14]. However, the high cost of this process and various issues surrounding the production of plant oils for non-food purposes make the search for alternative routes both attractive and strategically pertinent. In this respect, microbial production of biofuels (so-called microdiesel and microkerosene) represents a sustainable and quite economical way to produce FAEs. For this purpose, both Escherichia coli and S. cerevisiae have been engineered to produce structurally tailored fatty esters [1517]. However, neither of these microorganisms is naturally able to accumulate high amounts of lipids, nor able to degrade cellulose. Moreover, in these microorganisms the biosynthesis of fatty acid is highly regulated [18], thus limiting the possibility to improve lipid production [16, 17, 19].

So-called oleaginous microorganisms, which naturally accumulate lipids to more than 20% of their dry cell weight (DCW) [20, 21], have already been exploited for the production of commercially useful lipids, such as substitutes for cocoa butter and polyunsaturated fatty acids [22]. Therefore, it is unsurprising that microbial lipid or single cell oil is also being considered for biodiesel production, especially because this route implies shorter production times, reduced labor costs and simpler scale-up [23]. Prominent among the oleaginous microorganisms, Yarrowia lipolytica has been extensively studied and is known to accumulate lipids up to 50% of its dry weight depending on culture conditions [20, 21, 24]. Advantageously, since Y. lipolytica is already widely used in the detergent, food, pharmaceutical and environmental industries it has been classified by the FDA (Food and Drug Administration) as “Generally Recognized as Safe” (GRAS) for numerous processes [25]. Nevertheless, despite these advantages, Y. lipolytica displays limited ability for sugar use and is unable to use cellulose as carbon source [26].

In a recent paper, the use of cellobiose by Y. lipolytica was tackled for the first time, thus opening the way towards the development of an efficient yeast-based CBP microorganism capable of consuming cellulose-derived glucose and converting it into lipids and derivatives thereof [27]. Herein, we present work that shares this aim, but which has employed a different strategy that relies upon the activation of endogenous β-glucosidase activity (Fig. 1).

Fig. 1
figure 1

The strategies used in the current study to develop the cellobiose-degrading ability in Y. lipolytica.


Identification of genes encoding active β-glucosidases in Y. lipolytica

Analysis of the Y. lipolytica genome using BLAST revealed the presence of six sequences that were identified as putative family GH3 β-glucosidases (see Additional file 1: Table S1; Additional file 2: Fig. S1) on the basis of high amino acid sequence identity with other yeast β-glucosidases (Fig. 2, Additional file 2: Fig. S1). However, in the absence of biochemical data it was impossible to assert at this stage that these sequences actually encode β-glucosidases, since family GH3 contains glycoside hydrolases that display other specificities. Moreover, Y. lipolytica does not grow on cellobiose and has not been found to express a detectable level of β-glucosidase activity (Additional file 2: Fig. S2), despite the fact that preliminary transcriptional analysis revealed that the six genes are weakly transcribed when Y. lipolytica is grown on glucose, although no further induction was observed in the presence of cellobiose (Additional file 2: Fig. S3). In this respect, it was particularly gratifying to observe that overexpression of BGL1 (YALI0F16027g) or BGL2 (YALI0B14289g) in Y. lipolytica (strains ZetaB 1 and ZetaB 2, respectively) enhanced the transcription of the genes and conferred the ability to grow on solid medium containing cellobiose as the sole carbon source (Additional file 2: Figs. S2, 3). Additionally, when these recombinant strains were grown on YNB-pNP-βGlc (p-nitrophenyl-β-d-glucoside) plates, yellow halos surrounding the colonies were clearly visualized, indicating β-glucosidase activity (Additional file 2: Fig. S3). Finally, after growth in liquid YTD medium, β-glucosidase activity could be measured in the cell extract of ZetaB 1 (3.2 ± 0.2 IU/mg) and in the culture supernatant of ZetaB 2 (2.6 ± 0.1 U/mL), while much lower activities were measured in the culture supernatant of ZetaB 1 (0.33 ± 0.02 U/mL) and in the cell extract of ZetaB 2 (0.42 ± 0.01 IU/mg). Regarding the remaining four putative β-glucosidases, expression of their encoding sequences (YALI0F01672g, YALI0D18381g, YALI0B14333g and YALI0E20185g) in Y. lipolytica failed to produce any detectable β-glucosidase activity or sustain yeast growth on solid medium containing cellobiose as the sole carbon source.

Fig. 2
figure 2

Phylogenetic tree of the structurally characterized family three glycosyl hydrolases. Multiple-sequence alignment was performed using CLUSTALW (, and a phylogenetic tree was constructed using MEGA5 ( Sequences are: Bacillus subtilus Bgl (PDB accession number 4GYJ_A); Streptomyces venezuelae Bgl (PDB accession number 4I3G_A); Kluyveromyces marxianus Bgl (PDB accession number 3AC0_A); Aspergillus aculeatus Bgl (PDB accession number 4IIB_A); Thermotoga neapolitana Bgl (PDB accession number 2X42_A); Trichoderma reesei Bgl (PDB accession number 4I8D_A); and include the 6 putative GH3 sequences from Yarrowia lipolytica genome.

To further investigate the production of the six β-glucosidases, western blot analysis was carried out using anti-His6 antibodies. This revealed the presence of Bgl1-His6 in both the culture supernatant and cell extract (Fig. 3a, b). Bgl2-His6 could only be detected in concentrated (tenfold) cell extract (Fig. 3c). However, this method failed to detect Bgl2-His6 in the culture media (Fig. 3d), consistent with the fact that no β-glucosidase activity was detected in this fraction. Significantly, the expression of the native BGL2 sequence (i.e. without the His6-tag) provided much more satisfactory expression, implying that the presence of the His6-tag on Bgl2 somehow impairs the expression and/or secretion of this protein.

Fig. 3
figure 3

Western blot detection of the expressed β-glucosidases a M, molecular weight standards; lane 1, intracellular Bgl1, b lane 1, extracellular Bgl1-His6, c lane 1, intracellular Bgl2-His6, d lane 1, extracellular Bgl2-His6, and SDS-PAGE analysis of the purified β-glucosidases from Y. lipolytica JMY1212 transformants e lane 1, purified Bgl1-His6, and f lane 1, purified Bgl2; lane 2, endo-H treated Bgl2 (The lower band in lane 2 represents the expected size of Endo-H).

Localization of β-glucosidases in ZetaB 1 and ZetaB 2

To determine the localization of Bgl1 and Bgl2, yeast cells producing these enzymes (without His6 tag) were fractionated, generating on the one hand extracellular samples (culture supernatant), and on the other cell-associated periplasmic, cytoplasmic and membrane fractions. Measurement of the β-glucosidase activities in each of these fractions revealed that Bgl1 was primarily localized in the periplasm (61%), but was also present in the cytoplasm (30%), while Bgl2 was mainly in the supernatant (80%), and to a lesser extent (26%) in the periplasm (Table 1). The apparent ambiguity of these results is undoubtedly due to the separation method that was employed to isolate the different cellular fractions, which inevitably led to a low level of cross-contamination between the samples. Nevertheless, taking this into account, it is reasonable to deduce that Bgl1 is primarily localized in the periplasmic space, while Bgl2 is secreted into the culture medium.

Table 1 Distribution of β-glucosidase activity in recombinant strains ZetaB 1 and ZetaB 2

Production, purification and characterization of Bgl1 and Bgl2

Production of Bgl1-His6 and native Bgl2 was achieved by growing the appropriate Y. lipolytica strains on YTD, with expression of both enzymes increasing until complete depletion of glucose was reached (36 h).

Regarding purification of Bgl1-His6, yeast cells from a 200-mL culture volume yielded approximately 550 U of enzyme in the crude cell extract. However, after one step of affinity-purification, only 17% of Bgl1 was recovered (Table 2). In the case of Bgl2, a two-step protocol using anion exchange chromatography and gel filtration allowed its purification to near homogeneity, but led to significant loss of protein (8.8% recovery). SDS-PAGE analysis of the two purified protein samples revealed that while the M r of Bgl1 was consistent with the expected value (i.e. M r of 90.4 kDa for the protein lacking the putative signal peptide), that of Bgl2 was significantly higher (>250 kDa) (Fig. 3e, f). To investigate whether this discrepancy was due to glycosylation, the amino acid sequence of Bgl2 was analyzed using the glycosylation predictor GlycoEP ( [28]. This revealed that Bgl2 harbors 18 potential N-glycosylation sites (Additional file 1: Table S2). Glycosylation was finally confirmed by treating the purified Bgl2 protein with endoglycosidase H (EndoH) and migrating it on a SDS-PAGE. This analysis revealed that the M r of the recombinant EndoH-treated Bgl2 was approximately 95 kDa, quite consistent with the theoretical Mr of 92.9 kDa (Fig. 3f). Finally, N-terminal amino acid sequence analysis of Bgl1 and Bgl2 confirmed that the signal sequences of both proteins had been cleaved and allowed the accurate localization of the cleavage sites (Additional file 2: Fig. S4).

Table 2 Purification of intracellular Bgl1-His6 and extracellular Bgl2 produced by Y. lipolytica overexpressing strains

Preliminary characterization of Bgl1 and Bgl2 using pNP-βGlc as substrate revealed that Bgl1 was 5-fold more active (102.8 U/mg) on this substrate than Bgl2 (25.8 U/mg). The activity of Bgl1 was highest at approximately pH 4.5 and 45°C, and was stable in the pH range of 4.0–5.0 and below 40°C. Regarding Bgl2, it was found to display highest activity at pH 4.0 and 50°C and was stable in the pH range of 3.5–7.0 and below 50°C (Additional file 2: Figs. S5, S6). It is noteworthy that deglycosylation of Bgl2 led to a 60% decrease in specific activity, which was probably due to its instability at 40°C (Additional file 2: Fig. S7).

Substrate specificity and kinetic parameters of Bgl1 and Bgl2

The substrate specificity of the purified β-glucosidases was examined using different substrates displaying α and β configurations. The results showed that both β-glucosidases were maximally active against pNP-βGlc (Fig. 4). However, using activity on pNP-βGlc as benchmark, it is noteworthy that both enzymes were active on pNP-β-d-cellobioside (Bgl1, 24% and Bgl2, 27%), but only Bgl1 displayed significant activity (10%) on pNP-β-d-xylopyranoside. Neither enzyme displayed activity on pNP-β-d-galactopyranoside and pNP-α-d-glucopyranoside.

Fig. 4
figure 4

Comparison of the hydrolytic activity of β-glucosidases from Y. lipolytica JMY1212. Bgl1-His and Bgl2 on a pNP-derived substrates, and b natural glycosyl substrates with different β-configurations.

When the activity of Bgl1 and Bgl2 on cellobiose was compared with that on other oligosaccharides, it was found that both enzymes displayed highest activity on laminaribiose (β-1, 3-linkage), followed by gentiobiose (β-1, 6-linkage), octylglucoside, sophorose (β-1, 2-linkage), cello-oligosaccharides (β-1, 4-linkage, trisaccharide and higher)) and cellobiose (β-1, 4-linkage). It is noteworthy that the hydrolytic activity of Bgl1 was less dependent on the chain length of cello-oligosaccharides, while hydrolytic activity of Bgl2 increased as the length of cello-oligosaccharides increased. Both enzymes recognized methylglucoside as substrate, but the hydrolytic activities were low compared with the other substrates (Fig. 4), indicating that correct occupation of subsite +1 is important for catalysis.

The determination of the apparent kinetic parameters of reactions catalyzed by Bgl1 and Bgl2 and containing various glucosyl disaccharides and cello-oligosaccharides revealed that the values of K M(app) and k cat /K M for Bgl2-catalyzed reactions increased as a function of degree of polymerization (DP) of the cello-oligosaccharides (Table 3). In the case of Bgl1, increased DP was associated with increased K M(app) values, but not k cat /K M values. Overall, considering the performance constant (k cat /K M), cellobiose and cellohexaose were the best substrates for Bgl1 and Bgl2, respectively. Additionally, the performance constant of Bgl1 measured on cellobiose was 12.5-fold higher than that describing Bgl2. Regarding other glucosyl substrates (i.e. those containing linkages other than β-1, 4), both Bgls displayed the highest performance constants on laminaribiose. Nevertheless, comparison of the performance constants on each of the substrates revealed that Bgl2 is less regioselective, since the k cat /K M values were always lower in reactions catalyzed by Bgl1 (86% for sophorose, 47% for laminaribiose, 37% for gentiobiose, 18% for methylglucoside and 45% for octylglucoside). Finally, the lowest performance constants for both Bgls were measured for reactions containing methylglucoside.

Table 3 Kinetic parameters of Y. lipolytica Bgls for various glycoside-substrates

Cellobiose and cello-oligosaccharide fermentation with Y. lipolytica recombinant strains

Yeast strains ZetaB 1 and ZetaB 2, expressing BGL1 and BGL2, respectively, were grown in micro cultivation plates under aerobic conditions in the presence of cellobiose or cellodextrins (until Glc × 6) as sole carbon sources, using wild-type Y. lipolytica ZetaW as the control. The maximum specific growth rates (μmax) of the transformants on cellobiose were essentially the same as that of the control grown on glucose (Fig. 5a, b). ZetaB 1 grew faster than ZetaB 2 on cellobiose and cellodextrins (Fig. 5b–f), while the control was unable to grow on either of these substrates. Surprisingly, despite indications that the wild-type strain cannot grow on cellobiose or cellodextrins, the control culture (ZetaW) reached an OD value of 2.0. However, further investigation using a defined medium revealed that this unexpected growth could be attributed to the presence in the medium of 0.2% w/v casamino acids, which acted as a suitable carbon source. In defined medium, ZetaB 1 consumed 80% cellobiose over 48 h. However, upon further incubation, the remaining 20% cellobiose (at the concentration of 2 g/L) was not consumed (Fig. 6). In contrast, ZetaB 2 consumed all of the cellobiose over 60 h. Furthermore, ZetaB 1 sustained a specific growth rate (μmax) of 0.10 h−1, whereas ZetaB 2 exhibited a long lag phase on cellobiose after which two subsequent growth phases (μmax values of 0.08 h−1 and then 0.16 h−1) were observed (Fig. 6; Table 4). In order to combine the advantages procured by the overexpression of BGL1 and BGL2 (i.e. shorter lag phase and higher cellobiose utilization, respectively), the two BGL sequences were cloned into JMY1212, thus yielding ZetaB 12. During cultivation on cellobiose, the performance of ZetaB 12 was the best among all the recombinant strains. It showed similar growth rate to that of the control grown on glucose and consumed 10 g/L of cellobiose within 40 h.

Fig. 5
figure 5

Comparison of Y. lipolytica ZetaW (control), ZetaB 1 (P TEF -BGL1) and ZetaB 2 (P TEF -BGL2) during aerobic growth on 5 g/L a glucose, b cellobiose, c cellotriose, d cellotetraose, e cellopentaose and f cellohexaose as carbon and energy source. Shown is OD600nm, optical density at 600 nm, versus time. Each data point represents the mean of at least three independent experiments and the standard deviation is less than 5%.

Fig. 6
figure 6

Comparison of Y. lipolytica a ZetaW (control), b ZetaB 1 (P TEF -BGL1), c ZetaB 2 (P TEF -BGL2) and d Zeta-B12 (P TEF -BGL1, P TEF -BGL2) during aerobic growth on 10 g/L cellobiose. Shown are OD600nm, optical density at 600 nm, and cellobiose concentration versus time. Each data point represents the mean of five independent experiments and the error bar indicates the standard deviation.

Table 4 Comparison of growth and biomass yield of Y. lipolytica JMY1212 control and recombinant strains in aerobic cellobiose cultivation on defined medium

Characterization of cellulose-based lipid production by recombinant Y. lipolytica strains

To investigate whether the BGLs described in this study could be used to construct a lipid-producing strain, a previously described strategy to increase lipid accumulation, involving the deletion of the 6 POX genes (POX1 to POX6) that encode the peroxisomal acyl-coenzyme oxidases involved in lipid β-oxydation, was adopted [29]. Accordingly, Y. lipolyticapoxB1, ∆poxB2 and ∆poxB12 were constructed and grown on cellulose in the presence of Celluclast 1.5L. Even though this cocktail is reputedly β-glucosidase-deficient, to avoid any problems (i.e. spurious results linked to the presence of β-glucosidase in Celluclast) the Celluclast loading was kept low (7.5 FPU/g cellulose), and control experiments containing the prototrophic Y. lipolyticapoxW strain grown in the presence of Celluclast 1.5L with or without β-glucosidase supplementation were performed. During the initial 6 h of cultivation an accumulation of reducing sugars was observed in all of the cultures, which was attributed to Celluclast 1.5L-mediated cellulose hydrolysis. However, further monitoring revealed that after 12-h growth, less reducing sugars were present in the Y. lipolyticapoxB12 culture (2.7 g/L) compared to the other cultures (Fig. 7a). Moreover, this observation was correlated with continued yeast growth, whereas the growth of the other cultures stagnated over the same period (Fig. 7b). After 60 h of cultivation, the growth of ∆poxB12 reached a stationary phase. At this point the amount of FAMEs obtained after the methylation of Total fatty acids produced by this yeast had reached 0.8 g/L (Fig. 7c), but further growth did not result in an increase in cellular lipid content, reflecting a limitation of the available energy source. Besides the longer lag phase, the growth of ∆poxB1 and ∆poxB2 was similar to that of the control culture supplemented with β-glucosidase. Regarding the control culture, in the absence of β-glucosidase supplementation, growth ceased after 60 h and the cell density of the culture was approximately half that of the other cultures. Moreover, continuous addition of cellulases to the control culture did not procure any obvious increase in growth. When the control was supplemented with β-glucosidase, the amount of cellulose that remained unconsumed (25 g/L) was similar to that of cultures of the ∆pox strain expressing β-glucosidases after 5 days of growth.

Fig. 7
figure 7

Growth and lipid production on cellulose medium of Y. lipolytica strains. Growth during SSF on 50 g/L cellulose supplemented with Celluclast 1.5L. a the concentration of reduced sugar versus time; b growth expressed as cell number versus time; c lipid content at 60 h. Strains are Y. lipolyticapoxB1 (P TEF -BGL1), ∆poxB2 (P TEF -BGL2), ∆poxB12 (P TEF -BGL1, P TEF -BGL2) and ∆poxW (wild type) under the same condition without (control) or with (control + Bgl) extra β-glucosidase (Novozyme 188). Each data point represents the mean of at least three independent experiments and the error bars indicate the standard deviation.


Even though fungal β-glucosidases can be produced at relatively low cost via solid-state fermentation [30], cellulases still account for almost 50% of the cost of cellulose hydrolysis processes [31]. Moreover, taking into account the fact that the T. reesei secretome is rather deficient in β-glucosidases, it is particularly relevant to engineer microorganisms that are self-sufficient with regard to this type of enzyme activity.

To date, S. cerevisiae has been the main target for engineering aimed at the creation of a cellulolytic yeast strain for consolidated bioprocessing purposes [32, 33]. However, Y. lipolytica is one of the most widely studied “non conventional” oleaginous yeast species, which is well characterized for its ability to use hydrophobic substrates (e.g. alkanes, FAs, oils), glucose and glycerol as carbon sources for the production of lipids [34]. Nevertheless, Y. lipolytica is also known for its inability to grow on cellulose or even cellobiose. Therefore, it also constitutes an attractive target for strain engineering work.

Interestingly, comparative genomics has revealed that Y. lipolytica is only distantly related to the majority of yeast species and instead shares a number of common properties with filamentous fungi [35]. Therefore, based on this observation, we investigated whether Y. lipolytica harbors genes that allow cellobiose degradation and thus the possibility to confer cellobiose-degrading ability to Y. lipolytica through the use of endogenous β-glucosidases. Gratifyingly, our data clearly revealed that at least two genes, designated BGL1 and BGL2, encode β-glycosidases, active under the control of TEF promoter, that hydrolyze cellobiose, although our results also indicate that these enzymes do not exclusively cleave β-1, 4 linkages. In this respect, it is noteworthy that many cellobiolytic yeasts, such as Debaryomyces vanrijiae [36], Candida peltata [37], Monascus purpureus [38], Kluyveromyces fragilis [39] and Metschnikowia pulcherrima [40], only produce one active Bgl and even in the case of S. fibuligera, which produces two active Bgls, only one of these (SfBgl1) actually hydrolyzes β-1,4 glycosidic bonds [41]. In this respect it is also intriguing to note that while Bgl1 (described herein) displays good activity on both cellobiose and cellodextrin, its homolog from S. fibuligera (SfBgl2, 49.5% identity) is inactive on cellobiose. Regarding the remaining four putative β-glucosidases, at this stage it is unclear why no recombinant products were revealed by western blot analysis. To pursue this further it would no doubt be useful to take a closer look at transcription, but this was not performed since the successful expression of two BGLs was more than enough to continue the present study.

The multiplicity of β-glucosidases in filamentous fungi is common and is often due to the presence of multiple genes, or differential post transcriptional modifications [42, 43]. However, our study provides the first evidence of possible β-glucosidase multiplicity in yeast. In filamentous fungi, the apparent redundancy of multiple cellulolytic enzymes can be explained by the ability of these microorganisms to adapt to different biomass resources and culture conditions [43] and is probably essential in fungal metabolism and survival. Similarly, in the case of Y. lipolytica, the presence of multiple putative β-glycosidase genes (two of which have been shown to be β-glucosidases) could be the result of adaptation of the yeast to changing environments and may help to explain the different evolutionary history of this yeast.

Considering the physicochemical characteristics of the two Y. lipolytica β-glucosidases described in this study, it is possible to tentatively correlate these with the cellular localization of the enzymes. Bgl2 was more stable in the assay conditions employed (50°C, pH 4.0), which might arguably be logical for an extracellular enzyme that needs to show a certain resilience to environmental challenges. On the other hand, the relative fragility of Bgl1 can be explained by the fact that its natural intracellular location probably protects it from major temperature and pH changes. In this respect the different cellular locations of the enzymes is no doubt beneficial to the yeast, since it provides optimal activity.

The comparison of the specific activities of Bgl1 and Bgl2 on cellobiose (108 units/mg and 25 units/mg protein, respectively) with that of the commercially available β-glucosidase from Aspergillus niger (5.2 units/mg protein), the enzyme that is generally used to complement the cellulolytic cocktail of T. reesei [44], is rather flattering for the former. Moreover, the K M values describing the cellobiolytic reactions catalyzed by Bgl1 and 2 are approximately 10 and 3.4-fold lower than those of the β-glucosidases from S. fibuligera (2.8 mM Bgl1) and A. niger (2.7 mM) [44], meaning that the minimum concentration of cellobiose required for effective catalysis to occur is much lower. Likewise, comparing the apparent performance constants, k cat/K M, of Y. lipolytica Bgls with those of other reported β-glucosidases [3641, 44] suggests that the enzymes described in this study hydrolyze cellodextrins more efficiently. In this respect, it is also interesting to consider the fact that Bgl1 apparently exhibits higher catalytic performance on low DP cellodextrins, while Bgl2 is more active on higher DP ones. This difference in substrate preference is consistent with the cellular location of the two enzymes (i.e. the extracellular enzymes deals with the longer oligosaccharides that cannot be transported into the cell) and possibly forms the basis of the superior performance of ZetaB 12 on Avicel compared to control cultures containing the β-glucosidase from Novozymes 188.

Apparently, in the lipid production experiments performed in this study, the activity of the Yarrowia β-glucosidases adequately satisfied the requirements for exo-cellulase activity, since no accumulation of reducing sugars was observed. However, lipid production was not sustained, since growth ceased before the cellulose was consumed. This failure suggests that cellulose degradation was the limiting factor, with a part of the cellulose being recalcitrant to further hydrolysis by the Celluclast cocktail, even when Bgl1 and 2 were correctly expressed and thus available.


This study has provided a clear demonstration that Y. lipolytica does not naturally use cellobiose, despite the fact that this strain contains the genetic potential to do so. This is intriguing because the protein products of these genes are active on glucose-based oligosaccharides, including cello-oligosaccharides. Moreover, our data clearly show that upon expression of BGL1 and BGL2 under the control of the constitutive well characterized TEF promoter in Y. lipolytica, growth on cellobiose becomes possible. These encouraging findings render plausible the creation of an engineered Y. lipolytica strain that could be useful in advanced generation biorefinery schemes involving the use of lignocellulosic hydrolysates as feedstock for the production of bioenergy and valuable chemicals.


Strains and media

The genotypes of the microbial strains used in the present study are summarized in Table 5. E. coli DH5 were purchased from Invitrogen (Paisley, UK) and used for plasmid construction. The Y. lipolytica strains were routinely grown in YPD (1% yeast extract, 2% bacto peptone, and 2% glucose). Solid YPD medium contains 1.5% agar. Transformants were selected on solid YNB medium (0.17% w/v YNB, 1% glucose or cellobiose w/v, 0.5% w/v ammonium chloride, with (for Ura+) or without (for Leu+) 0.2% w/v casamino acids and 50 mM sodium–potassium phosphate buffer, pH 6.8), supplemented with uracil (440 mg/L) or leucine (440 mg/L) depending on the auxotrophic requirements. The detection of β-glucosidase activity in solid YNBcasa medium was achieved by incorporating 1.0 mM p-nitrophenyl-β-d-glucoside (pNP-βGlc) [45]. For β-glucosidase characterization, enzymes were produced in YTD medium (1% w/v yeast extract, 2% w/v tryptone, 5% w/v glucose and 100 mM phosphate buffer, pH 6.8). To compare the efficiency of recombinant β-glucosidase to degrade cellobiose and cellodextrin with respect to cell growth, yeasts were aerobically cultivated in YNBcasa medium, containing 5 g/L cellobiose or cello-oligosaccharides (C3–C6), and defined medium containing vitamins, trace elements [46] and salts, including 3.5 g/L (NH4)2SO4, 3.0 g/L K2HPO4, 3.0 g/L NaH2PO4 and 1.0 g/L MgSO4∙7H2O with 10 g/L cellobiose. For lipid production using cellulose as the carbon source, Y. lipolytica strains were grown in defined media supplemented with 50 g/L Avicel PH-101.

Table 5 Microbial strains used in the present study

Plasmid constructions

The plasmids constructed in the present study are summarized in Table 6, and all primers are listed in Table 7. The plasmids used for expression of the putative β-glucosidases were constructed using the expression vectors JMP62UraTEF and JMP62LeuTEF, which are derivatives of a previously described vector [47]. Briefly, these vectors contain the Y. lipolytica TEF promoter and either the URA3ex or LEU2ex excisable selection markers, which are flanked by loxP sites and a Zeta fragment that serves as the homologous integration site [48]. Regarding β-glucosidases, six putative gene candidates (Sequences YALI0F16027g, YALI0F01672g, YALI0D18381g, YALI0B14289g, YALI0B14333g, YALI0E20185g available at Genome Resources from Yeast Chromosomes: were identified (see Additional file 1: Table S1). For the expression of wild-type and His6-tagged proteins, the genes were amplified by PCR using FA (1–6) as forward primers and RB (1–6) or RB-His (1–6) as reverse primers, respectively. The PCR fragments were digested using either BamHI/AvrII or HindIII/AvrII and inserted into the plasmid JMP62 UraTEF at the corresponding sites.

Table 6 Plasmids used or created in the present study
Table 7 The sequences of the oligonucleotide primers used in this study

After construction, all expression vectors were verified by DNA sequencing (GATC Biotech, Konstanz, Germany). For Y. lipolytica transformation, vectors were digested using NotI, thus generating a linear DNA with Zeta sequences at both extremities, and purified. Then the linear DNA fragments were introduced into the Zeta docking platform of Y. lipolytica JMY1212 Zeta or randomly into the genome of ∆pox strain using the lithium acetate method [49]. Transformants were tested for β-glucosidase activity on YNB glucose plate containing pNP-βGlc and for growth on cellobiose using solid YNB cellobiose plates. Clones displaying both activities were retained for further analysis.

Transcriptional analysis

Y. lipolytica JMY1212 wide type and recombinant strains overexpressing BGL1 and BGL2 were grown to mid-exponential phase in defined media and then transferred into fresh medium containing either glucose or cellobiose as the sole carbon source. Cells were recovered from the medium at 20 min and 1 h, respectively, and rapidly frozen in liquid nitrogen and stored at −80°C until use. Total mRNA was isolated using RNeasy Plus Mini Kit (QIAGEN) and reverse transcription was performed with iScript™ cDNA Synthesis Kit (BIO-RAD) according to the manufacturer’s instructions. Transcription of the BGLs was analyzed by PCR, using gene-specific primers and sequencing of the PCR products (Additional file 1: Table S3).

Measurement of enzyme activity

β-Glucosidase activity was measured by quantifying the release of pNP (p-nitrophenol) from pNP-βGlc as described previously [50]. One unit of pNP-βGlcase activity was defined as the amount of enzyme required to release 1 μmol pNP per min. All protein concentrations were measured using the Bradford method and bovine serum albumin as a standard [47].

Western blot analysis

Western blotting of proteins was performed as described previously [51]. Crude supernatant and cell-free extracts of Y. lipolytica JMY1212 expressing putative β-glucosidases fused with the His6 tag were concentrated 10-fold using an ultra-centrifugation filter unit (Amicon® Ultra-4 10 kDa cut-off, Merk Millipore, Bedford, MA, USA). Blots were sequentially treated with mouse non position-specific His-Tag antibody 1:2,500 (THE™ from Genscript, Piscataway, NJ, USA) and the alkaline phosphatase-conjugated goat anti-mouse IgG.

Subcellular fractionation and enzyme localization

Fractionation of yeast cells was carried out as described by Cummings and Fowler [52], with slight modifications. Briefly, yeasts were cultivated until a cellular density of 6 × 107 cells/mL was reached. Then, to quantify total β-glucosidase activity, a 50-mL sample was taken and subjected to centrifugation at 8,000×g for 5 min at 4°C thus isolating a cell pellet and supernatant. The cell pellet was disrupted in Tris–HCl buffer (50 mM, pH 7.4, 3 mM EDTA and 0.5 mM PMSF) using a MP FastPrep-24 Instrument (MP Biomedicals Inc.). β-Glucosidase activity in both the cell lysate and the supernatant was determined as described earlier to estimate total β-glucosidase activity. Using a second 50 mL yeast culture, a cell pellet containing approximately 2 × 108 cells/mL was obtained by centrifugation and then treated with zymolyase 100T at 10 mg/mL (Seikagaku corp coger) in 15 mL of sorbitol buffer (1 M sorbitol, 50 mM Tris–HCl, pH 7.4, 2 mM dithiothreitol, 10 mM MgCl2, 20 mM-sodium azide, 0.5 mM PMSF) at 30°C with gentle shaking. Protoplast formation was monitored using a microscope until ≥99% of the cells was lysed when SDS was added (1% SDS w/v). The solid protoplast fraction was then separated from the supernatant by centrifugation (1,000 rpm for 5 min at 4°C) and the latter was designated as the periplasmic fraction. The protoplasts were re-suspended in Tris–HCl buffer (50 mM Tris–HCl, pH 7.4) and disrupted by vortex in the presence of glass beads (0.4–0.45 mm). The homogenate was centrifuged (20,000×g for 2 h at 4°C) and the supernatant and solid fractions were designated as the cytoplasmic and membrane fraction, respectively. Prior to enzyme assays, the membrane fraction was suspended in citrate buffer.

Purification of β-glucosidases

Y. lipolytica JMY1212 overproducing Bgl1-His6 and Bgl2 was grown in 200 mL YTD medium at 130 rpm, 28°C for 36 h before centrifugation at 8,000×g for 5 min. For purification of His6-Bgl1, the cell pellet was washed, suspended in 50 mL phosphate buffer (50 mM, pH 7.4) and homogenized over a 3-min period using a MP FastPrep-24 Instrument. After centrifugation (8,000×g for 5 min at 4°C), the supernatant was applied to 2 mL of TALON Metal Affinity Resin (Clontech, Takara-Bio, Kyoto, Japan) and protein was eluted using imidazole buffer according to the manufacturer’s instructions.

For purification of Bgl2, the culture supernatant was concentrated fivefold using an Amicon® Ultra-4 Centrifugal Filter Unit with 30 kDa cut-off (Merk Millipore, Bedford, MA, USA). The concentrated sample was then loaded onto a Q Sepharose™ High Performance column (Hiload, 1.6 × 10 cm, Pharmacia Biotech), equilibrated with Tris-buffer (20 mM, pH8.0). The column was washed first with equilibration buffer (2 bed volumes) before applying a linear gradient of 0–1.0 M NaCl in Tris-buffer (20 mM, pH7.4) at a flow rate of 1.0 mL/min (Pharmacia Biotech ÄKTA). Eluted fractions were collected and assayed for β-glucosidase activity. All fractions displaying activity were pooled, desalted and concentrated using an Amicon ultra-filtration unit equipped with a PM-10 membrane (Millipore), before being applied to a Superdex 200 column (1.0 × 30 cm, Pharmacia Biotech) equilibrated in Tris-sodium buffer (20 mM Tris–HCl, 150 mM NaCl, pH 7.4). Protein species were separated at a flow rate of 0.5 mL/min. Fractions were collected and analyzed by SDS-PAGE to ascertain purity and estimate the approximate molecular weights of Bgl1-His6 and Bgl2. All fractions satisfying the purity criterion (>95% purity) were pooled and retained for further work.

Deglycosylation and N-terminal amino acid sequencing

Purified Bgl1-His6 and Bgl2 were treated with endoglycosidase H (New England Biolabs, Beverly, MA, USA) according to the manufacturer’s instructions. After deglycosylation, the protein species displaying M r (relative molecular mass) closest to those of the theoretical M r (predicted using Protparam, of Bgl1-His6 and Bgl2 were excised and submitted to N-terminal amino acid sequencing (PISSARO platform, Rouen, France).

Physicochemical characteristics of β-glucosidases

Optimal temperatures and pH for the activity of Bgl1 and Bgl2 were determined using pNPGlc as the substrate. Assays were either performed at pH 5.0 and various temperatures (30–70°C) or at 30°C in variable pH conditions (2.0–8.0) using either 50 mM glycine–HCl (pH 2.0), 50 mM citrate/acetate (pH 3.0–7.2) or potassium phosphate (pH 7.0–8.2) buffer. When the temperature was varied, the pH of the citrate buffer was adjusted accordingly. Stability of Bgl1 and Bgl2 depending on pH and temperature was analyzed as follows: enzymes were incubated at 30°C for up to 2 h at various pH values (2.0–8.0) or at various temperatures (30–70°C) for up to 2 h in 50 mM citrate buffer, pH 5.0. Residual glucosidase activity was then assayed at 30°C in 50 mM citrate buffer, pH 5.0.

Substrate specificity and enzyme kinetics

The substrate specificity of Bgl1-His6 and Bgl2 was investigated by assaying for activity on the aryl-glycosides pNP-β-d-glucopyranoside, pNP-α-d-glucopyranoside, pNP-β-d-galactopyranoside, pNP-β-d-xylopyranoside and pNP-β-d-cellobioside, and on the oligosaccharides cellobiose, cellotriose, cellotetraose, cellopentaose, cellohexaose, sophorose, laminaribiose, gentiobiose, methylglucoside and octylglucoside. When using aryl-substrates, the standard assay method was employed, simply replacing pNP-βGlc by another substrate as appropriate. For oligosaccharides, the release of glucose was quantified using an enzyme kit (d-Fructose/d-Glucose Assay Kit, liquid stable, Megazyme). To study the Michaelis–Menten parameters K M, V max and k cat , Bgl1 (0.120 nM) or Bgl2 (0.13 nM) were added to reaction mixtures containing different substrate concentrations: 0.25–5 mM cellobiose, 0.25–5 mM cellotriose, 0.25–5 mM cellotetraose, 0.25–5 mM cellopentaose, 0.25–5 mM cellohexaose, 0.2–4 mM sophorose, 0.1–2 mM laminaribiose, 0.1–2 mM gentiobiose, 0.5–20 mM methylglucoside and 0.2–4 mM octylglucoside. Initial rates were fitted to the Michaelis–Menten kinetic equation using a nonlinear regression (SigmaPlot 10) to extract the apparent K M and k cat [53].

Yeast growth and lipid production

Yeast growth on cellobiose and cellodextrins was performed in a 40-well microplate. A single colony from a fresh YPD plate was transferred into 5 mL of defined medium containing 10 g/L of glucose and pre-cultured until the mid-exponential phase. The cells were then harvested, washed, suspended in sterile water and used to inoculate 200 μL YNBcasa media containing 5 g/L cellobiose or cellodextrins in the microplate, achieving an initial OD600 of 0.1. This culture was grown in a microplate reader (Spectrostar Omega, BMG Labtech, Germany) at 30°C with continuous shaking (150 rpm) and automatic OD600 recording. Similarly, yeast growth on cellobiose was also performed on 30 ml defined medium containing 10 g/L cellobiose in 250 mL Erlenmeyer flasks.

For lipid production a fresh yeast culture in exponential phase was used to inoculate 50 mL defined medium containing 50 g/L Avicel in Erlenmeyer flasks, achieving an initial OD600 of 1.0. Celluclast 1.5L (60 FPU/mL, gift from Novozymes, Denmark) was added (7.5 U/g cellulose) and growth was pursued for 5 days (30°C, 150 rpm). Samples were taken at regular intervals to determine concentrations of biomass, glucose, cellobiose and citric acid. In parallel, two control experiments were conducted under the same conditions, with or without the addition of extra β-glucosidase (810 IU/mL Novozyme 188, gift from Novozyme, Denmark) at 12.0 IU/g cellulose as recommended [54].

Analysis of product formation and determination of dry cell weight

To determine the concentration of substrates and extracellular metabolites, three aliquots (1.5 mL each) of cultures were rapidly frozen in liquid nitrogen and then thawed on ice before centrifugation (8,000×g for 5 min at 4°C) to recover supernatants for analysis. Glucose, cellobiose and citric acid were measured using an Aminex HPX87-H column (Bio-Rad Laboratories, Germany), operating at 50°C using a mobile phase (5 mM H2SO4) flowing at a rate of 0.5 mL/min. Glucose and cellobiose were detected using a Shodex RI-101 refractive index detector (Showa Denko, New York, NY, USA), while citric acid was detected using an UV detector at 210 nm (Dionex, Sunnyvale, CA, USA).

To determine the dry cell weight, three aliquots (5 mL each) of cultures were filtered through pre-weighed PES filters (0.45 μm; Sartorius Biolab, Germany). The biomass retained by the filters was washed, dried in a microwave oven at 150 W for 15 min and then placed in a desiccator before weighing. The biomass yield was calculated as the ratio of the amount of biomass obtained divided by the amount of carbon source consumed.

Lipids were extracted from freeze-dried cells (~10 mg) and methylated as described previously [55]. During the lipid extraction, C17:0 (Sigma) (50 μg) was added as the internal standard and fatty acid methyl esters (FAMEs) were analyzed by gas chromatography (6890 N Network GC System, Agilent, USA). The measurements were performed in a split mode (1 μL at 250°C), with helium as the carrier gas (2 mL/min). FAMEs were separated on a HP-5 GC column (30 m × 0.32 mm I.D., 0.5-μm film thickness, Agilent, USA). The temperature program was 120°C, ramped to 180°C (10°C/min) for 6 min, 183°C (0.33°C/min) for 9 min and 250°C (15°C/min) for 5 min. Detection was performed using a flame ionization detector (FID) at 270°C (2.0 pA). FAMEs were quantified by comparing their profiles with that of standards of known concentration.


LC biomass:

lignocellulosic biomass


consolidated bioprocessing


fatty acid ethyl esters


simultaneous saccharification and fermentation


Fatty Acid Esters


dry cell weight


single cell oil


Food and Drug Administration


generally recognized as safe



M r :

relative molecular mass




  1. Pedersen M, Meyer AS (2010) Lignocellulose pretreatment severity-relating pH to biomatrix opening. N Biotechnol 27(6):739–750

    Article  CAS  Google Scholar 

  2. Wilson DB (2009) Cellulases and biofuels. Curr Opin Biotechnol 20(3):295–299

    Article  CAS  Google Scholar 

  3. Tomme P, Warren RA, Gilkes NR (1995) Cellulose hydrolysis by bacteria and fungi. Adv Microb Physiol 37:1–81

    Article  CAS  Google Scholar 

  4. Sun Y, Cheng J (2002) Hydrolysis of lignocellulosic material for ethanol production: a review. Bioresour Technol 83:1–11

    Article  CAS  Google Scholar 

  5. Lynd LR, van Zyl WH, McBride JE, Laser M (2005) Consolidated bioprocessing of cellulosic biomass: an update. Curr Opin Biotechnol 16:577–583

    Article  CAS  Google Scholar 

  6. La Grange DC, Den Haan R, Van Zyl WH (2010) Engineering cellulolytic ability into bioprocessing organisms. Appl Microbiol Biotechnol 87:1195–1208

    Article  CAS  Google Scholar 

  7. Duff SJB, Cooper DG, Fuller OM (1985) Cellulase and beta-glucosidase production by mixed culture of Trichoderma reesei Rut C30 and Aspergillus phoenicis. Biotechnol Lett 7:185–190

    Article  CAS  Google Scholar 

  8. Holtzapple M, Cognata M, Shu Y, Hendrickson C (1990) Inhibition of Trichoderma reesei cellulase by sugars and solvents. Biotechnol Bioeng 36(3):275–287

    Article  CAS  Google Scholar 

  9. Stockton DC, Mitchell DJ, Grohmann K (1991) Optimum β-glucosidase supplementation of cellulose for efficient conversion of cellulose to glucose. Biotechnol Lett 13:57–62

    Article  CAS  Google Scholar 

  10. Lee WH, Nan H, Kim HJ, Jin YS (2013) Simultaneous saccharification and fermentation by engineered Saccharomyces cerevisiae without supplementing extracellular β-glucosidase. J Biotechnol 167(3):316–322

    Article  CAS  Google Scholar 

  11. Lian J, Li Y, HamediRad M, Zhao H (2014) Directed evolution of a cellodextrin transporter for improved biofuel production under anaerobic conditions in Saccharomyces cerevisiae. Biotechnol Bioeng 111(8):1521–1531

    Article  CAS  Google Scholar 

  12. Tang H, Hou J, Shen Y, Xu L, Yang H, Fang X et al (2013) High β-glucosidase secretion in Saccharomyces cerevisiae improves the efficiency of cellulase hydrolysis and ethanol production in simultaneous saccharification and fermentation. J Microbiol Biotechnol 23(11):1577–1585

    Article  CAS  Google Scholar 

  13. Zhang F, Carothers JM, Keasling JD (2012) Design of a dynamic sensor-regulator system for production of chemicals and fuels derived from fatty acids. Nat Biotechnol 30(4):354–359

    Article  CAS  Google Scholar 

  14. Demirbaş A (2003) Biodiesel fuels from vegetable oils via catalytic and non-catalytic supercritical alcohol transesterifications and other methods: a survey. Energy Convers Manage 44:2093–2109

    Article  Google Scholar 

  15. Steen EJ, Kang Y, Bokinsky G, Hu Z, Schirmer A, McClure A et al (2010) Microbial production of fatty-acid-derived fuels and chemicals from plant biomass. Nature 463:559–562

    Article  CAS  Google Scholar 

  16. Shi S, Octavio Valle-Rodriguez J, Khoomrung S, Siewers V, Nielsen J (2012) Functional expression and characterization of five wax ester synthases in Saccharomyces cerevisiae and their utility for biodiesel production. Biotechnol Biofuels 5:7–16

    Article  CAS  Google Scholar 

  17. Runguphan W, Keasling JD (2014) Metabolic engineering of Saccharomyces cerevisiae for production of fatty acid-derived biofuels and chemicals. Metab Eng 21:103–113

    Article  CAS  Google Scholar 

  18. Nielsen J (2009) Systems biology of lipid metabolism: from yeast to human. FEBS Lett 583(24):3905–3913

    Article  CAS  Google Scholar 

  19. Valle-Rodríguez JO, Shi S, Siewers V, Nielsen J (2014) Metabolic engineering of Saccharomyces cerevisiae for production of fatty acid ethyl esters, an advanced biofuel, by eliminating non-essential fatty acid utilization pathways. Appl Energ 115:226–232

    Article  Google Scholar 

  20. Ratledge C (2005) Single cell oils for the 21th century. In: Cohen R (ed) Single cell oils. AOCS Press, Champaign, pp 1–20

    Chapter  Google Scholar 

  21. Thevenieau F, Nicaud J-M (2013) Microorganisms as sources of oils. OCL 20(6):1–8

    Google Scholar 

  22. Papanikolaou S, Aggelis G (2010) Yarrowia lipolytica: a model microorganism used for the production of tailor-made lipids. Eur J Lipid Sci Technol 112:639–654

    Article  CAS  Google Scholar 

  23. Easterling ER, French WT, Hernandez R, Licha M (2009) The effect of glycerol as a sole and secondary substrate on the growth and fatty acid composition of Rhodotorula glutinis. Bioresour Technol 100(1):356–361

    Article  CAS  Google Scholar 

  24. Blazeck J, Hill A, Liu L, Knight R, Miller J, Pan A et al (2014) Harnessing Yarrowia lipolytica lipogenesis to create a platform for lipid and biofuel production. Nat Commun 5:3131

    Article  Google Scholar 

  25. Groenewald M, Boekhout T, Neuvéglise C, Gaillardin C, van Dijck PW, Wyss M (2014) Yarrowia lipolytica: safety assessment of an oleaginous yeast with a great industrial potential. Crit Rev Microbiol 40(3):187–206

    Article  CAS  Google Scholar 

  26. Michely S, Gaillardin C, Nicaud J-M, Neuvéglise C (2013) Comparative physiology of oleaginous species from the Yarrowia clade. PLoS One 8(5):e63356

    Article  CAS  Google Scholar 

  27. Lane S, Zhang S, Wei N, Rao C, Jin YS, Lee WH et al (2014) Development and physiological characterization of cellobiose-consuming Yarrowia lipolytica. Biotechnol Bioeng (In press)

  28. Chauhan JS, Rao A, Raghava GPS (2013) In silico platform for prediction of N-, O- and C-glycosites in eukaryotic protein sequences. PLoS One 8(6):e67008

    Article  CAS  Google Scholar 

  29. Beopoulos A, Mrozova Z, Thevenieau F, Le Dall MT, Hapala I, Papanikolaou S et al (2008) Control of lipid accumulation in the yeast Yarrowia lipolytica. Appl Environ Microbiol 74(24):7779–7789

    Article  CAS  Google Scholar 

  30. Singhania RR, Patel AK, Sukumaran RK, Larroche C, Pandey A (2013) Role and significance of beta-glucosidases in the hydrolysis of cellulose for bioethanol production. Bioresour Technol 127:500–507

    Article  CAS  Google Scholar 

  31. Rodrigues AC, Leitão AF, Moreira S, Felby C, Gama M (2012) Recycling of cellulases in lignocellulosic hydrolysates using alkaline elution. Bioresour Technol 110:526–533

    Article  CAS  Google Scholar 

  32. Jeon E, Hyeon JE, Suh DJ, Suh YW, Kim SW, Song KH et al (2009) Production of cellulosic ethanol in Saccharomyces cerevisiae heterologous expressing Clostridium thermocellum endoglucanase and Saccharomycopsis fibuligera beta-glucosidase genes. Mol Cell 28:369–373

    Article  CAS  Google Scholar 

  33. Shen Y, Zhang Y, Ma T, Bao X, Du F, Zhuang G et al (2008) Simultaneous saccharification and fermentation of acid-pretreated corncobs with a recombinant Saccharomyces cerevisiae expressing beta-glucosidase. Bioresour Technol 99:5099–5103

    Article  CAS  Google Scholar 

  34. Fickers P, Benetti P-H, Waché Y, Marty A, Mauersberger S, Smit M et al (2005) Hydrophobic substrate utilisation by the yeast Yarrowia lipolytica, and its potential applications. FEMS Yeast Res 5:527–543

    Article  CAS  Google Scholar 

  35. Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I et al (2004) Genome evolution in yeasts. Nature 430:35–44

    Article  Google Scholar 

  36. Belancic A, Gunata Z, Vallier MJ, Agosin E (2003) β-Glucosidase from the grape native yeast Debaryomyces vanrijiae. Purification, characterization, and its effect on monoterpene content of a muscat grape juice. J Agric Food Chem 51:1453–1459

    Article  CAS  Google Scholar 

  37. Galas E, Romanowska I (1996) Production purification and characterization of a highly glucose-tolerant novel ß-glucosidase from Candida peltata. Appl Environ Microbiol 62:3165–3170

    Google Scholar 

  38. Daroit DJ, Simonetti A, Hertz PF, Brandelli A (2008) Purification and characterization of an extracellular beta-glucosidase from Monascus purpureus. J Microbiol Biotechnol 18:933–941

    CAS  Google Scholar 

  39. Leclerc M, Chemardin P, Arnaud A, Ratomahenina R, Galzy P, Gerbaud C et al (1987) Comparison of the properties of the purified beta-glucosidase from the transformed strain of Saccharomyces cerevisiae TYKF2 with that of the donor strain Kluyveromyces fragilis Y610. Biotechnol Appl Biochem 9:410–422

    CAS  Google Scholar 

  40. González-Pombo P, Pérez G, Carrau F, Guisán JM, Batista-Viera F, Brena BM (2008) One-step purification and characterization of an intracellular beta-glucosidase from Metschnikowia pulcherrima. Biotechnol Lett 30:1469–1475

    Article  Google Scholar 

  41. Machida M, Ohtsuki I, Fukui S, Yamashita I (1988) Nucleotide sequences of Saccharomycopsis fibuligera genes for extracellular beta-glucosidases as expressed in Saccharomyces cerevisiae. Appl Environ Microbiol 54(12):3147–3155

    CAS  Google Scholar 

  42. Collins T, Gerday C, Feller G (2005) Xylanases, xylanase families and extremophilic xylanases. FEMS Microbiol Rev 29:3–23

    Article  CAS  Google Scholar 

  43. Nazir A, Soni R, Saini HS, Kaur A, Chadha BS (2010) Profiling differential expression of cellulases and metabolite footprints in Aspergillus terreus. Appl Biochem Biotechnol 162:538–547

    Article  CAS  Google Scholar 

  44. Yan TR, Lin CL (1997) Purification and characterization of a glucose-tolerant β-glucosidase from Aspergillus niger CCRC 31494. Biosci Biotechnol Biochem 61:965–970

    Article  CAS  Google Scholar 

  45. Guo ZP, Zhang L, Ding ZY, Gu ZH, Shi GY (2011) Development of an industrial ethanol-producing yeast strain for efficient utilization of cellobiose. Enzyme Microb Technol 49(1):105–112

    Article  CAS  Google Scholar 

  46. Verduyn C, Postma E, Scheffers WA, van Dijken JP (1992) Effect of benzoic acid on metabolic fluxes in yeasts: a continuous culture study on the regulation of respiration and alcoholic fermentation. Yeast 8:501–517

    Article  CAS  Google Scholar 

  47. Bradford MM (1976) A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72:248–254

    Article  CAS  Google Scholar 

  48. Fickers P, Le Dall MT, Gaillardin C, Thonart P, Nicaud J-M (2003) New disruption cassettes for rapid gene disruption and marker rescue in the yeast Yarrowia lipolytica. J Microbiol Methods 55(3):727–737

    Article  CAS  Google Scholar 

  49. Duquesne S, Bordes F, Fudalej F, Nicaud J-M, Marty A (2012) The yeast Yarrowia lipolytica as a generic tool for molecular evolution of enzymes. Methods Mol Biol 861:301–312

    Article  CAS  Google Scholar 

  50. Bordes F, Fudalej F, Dossat V, Nicaud J-M, Marty A (2007) A new recombinant protein expression system for high-throughput screening in the yeast Yarrowia lipolytica. J Microbiol Methods 70(3):493–502

    Article  CAS  Google Scholar 

  51. Duquesne S, Bozonnet S, Bordes F, Dumon C, Nicaud J-M, Marty A (2014) Construction of a highly active xylanase displaying oleaginous yeast: comparison of anchoring systems. PLoS One 9(4):e95128

    Article  Google Scholar 

  52. Cummings C, Fowler T (1996) Secretion of Trichoderma reesei β-glucosidase by Saccharomyces cerevisiae. Curr Genet 29:227–233

    CAS  Google Scholar 

  53. Segel I (1993) Enzyme kinetics: behavior and analysis of rapid equilibrium and steady-state enzyme systems. John Wiley, New York

    Google Scholar 

  54. Lan TQ, Lou H, Zhu JY (2013) Enzymatic saccharification of lignocelluloses should be conducted at elevated pH 5.5 to 6.2. Bioenergy Res 6((2):476–485

    Article  CAS  Google Scholar 

  55. Browse J, McCourt PJ, Somerville CR (1986) Fatty acid composition of leaf lipids determined after combined digestion and fatty acid methyl ester formation from fresh tissue. Anal Biochem 152:141–145

    Article  CAS  Google Scholar 

  56. Haddouche R, Poirier Y, Delessert S, Sabirova J, Pagot Y, Neuvéglise C et al (2011) Engineering polyhydroxyalkanoate content and monomer composition in the oleaginous yeast Yarrowia lipolytica by modifying the ß-oxidation multifunctional protein. Appl Microbiol Biotechnol 91(5):1327–1340

    Article  CAS  Google Scholar 

  57. Nicaud J-M, Madzak C, van den Broek P, Gysler C, Duboc P, Niederberger P et al (2002) Protein expression and secretion in the yeast Yarrowia lipolytica. FEMS Yeast Res 2(3):371–379

    CAS  Google Scholar 

Download references

Authors’ contributions

ZPG, SB, SD, JMN, AM, MJD conceived of the study and participated in its design. ZPG designed the constructs, carried out all the experiments and drafted the manuscript. SD participated in the western blot analysis of the expressed proteins. GC participated in the protein purification. SB, SD, JMN, AM, MJD revised the manuscript. All authors read and approved the final manuscript.


The authors would like to express their gratitude to Nelly Monties for her help with chromatographic analyses. They thank the ICEO facility dedicated to enzyme screening and discovery and part of the Integrated Screening Platform of Toulouse (PICT, IBiSA) for providing access to its protein purification and analytical facilities. This work was funded by the Agence Nationale de la Recherche (Investissements d’Avenir program; reference ANR-11-BTBR-0003). Dr. Guo is a recipient of a postdoctoral fellowship from the Institut National de la Recherche Agronomique.

Compliance with ethical guidelines

Competing interests The authors declare that they have no competing interests.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Zhongpeng Guo or Alain Marty.

Additional files

Additional file 1:

Table S1. Six putative β-glucosidase coding genes identified harboring conserved glycosyl hydrolase family 3 N and/or 3C terminal domain. Table S2. The predicted N-glycosylation sites in sequence of Bgl2 by GlycoEP. Table S3. The sequences of the oligonucleotide primers used in verification of transcriptions.

Additional file 2:

Figure S1. Multiple alignments of putative conserved domains of the family 3 glycosyl hydrolases of S. fibuligera (Bgl1, GenBank Accession numbers: AAA34314.1) against Yarrowia genome. Figure S2. Screening of Y. lipolytica expressing the 6 putative β-glucosidases on (a) indication plate containing YNBcasa medium supplemented with 1 mM p-nitrophenyl-β-D-glucoside (pNP-βGlc), and (b, c) YNBC plate with cellobiose as sole carbon source. Figure S3. Transcriptional analysis of the expression of the six putative BGLs in wild type strain on glucose (A) and cellobiose (B), and recombinant strain overexpression of BGL1 and BGL2 (C). Figure S4. N-terminal amino acid sequences of Y. lipolytica Bgl1 (A) and Bgl2 (B). The first 50 N-terminal amino acid sequences are indicated with the predicted signal sequence determined with signal P (underlined), the cleavage site predicted with a*, and the N-terminal AA sequence of the purified protein determined by direct sequencing (in bold). Figure S5. Optimal pH (a) and temperature (b) of Bgl1 (square) and Bgl2 (diamond) from Y. lipolytica JMY1212. Each data point represents the mean of three independent experiments and the error bar indicates the standard deviation. Figure S6. Stability of Bgl1 (a) and Bgl2 (b) from Y. lipolytica JMY1212 at pH from 2.0–8.0 as a function of time at 30ºC, and stability of Bgl1 (c) and Bgl2 (d) at temperature from 30ºC to 60ºC as a function of time at pH 5. Each data point represents the mean of three independent experiments and the error bar indicates the standard deviation. Please note that only one curve is given to represent the stability of Bgl2 at pH 4.0, 5.0 and 6.0 (b) and at 30ºC and 40ºC (d) as 100% of enzyme activity remained for these conditions. Figure S7. The hydrolytic activity of Bgl2 on pNP-βGlc (a) and the stability of Bgl2 at 40ºC as a function of time at pH5.0 before and after deglycosylation. Each data point represents the mean of three independent experiments and the error bar indicates the standard deviation.

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

Guo, Z., Duquesne, S., Bozonnet, S. et al. Development of cellobiose-degrading ability in Yarrowia lipolytica strain by overexpression of endogenous genes. Biotechnol Biofuels 8, 109 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: