Regulatory properties of malic enzyme in the oleaginous yeast, Yarrowia lipolytica, and its non-involvement in lipid accumulation

Malic enzyme (EC 1.1.1.40) converts l-malate to pyruvate and CO₂ providing NADPH for metabolism especially for lipid biosynthesis in oleaginous microorganisms. However, its role in the oleaginous yeast, Yarrowia lipolytica, is unclear. We have cloned the malic enzyme gene (YALI0E18634g) from Y. lipolytica into pET28a, expressed it in Escherichia coli and purified the recombinant protein (YlME). YlME used NAD⁺ as the primary cofactor. K_m values for NAD⁺ and NADP⁺ were 0.63 and 3.9 mM, respectively. Citrate, isocitrate and α-ketoglutaric acid (>5 mM) were inhibitory while succinate (5–15 mM) increased NADP⁺- but not NAD⁺-dependent activity. To determine if fatty acid biosynthesis could be increased in Y. lipolytica by providing additional NADPH from an NADP⁺-dependent malic enzyme, the malic enzyme gene (mce2) from an oleaginous fungus, Mortierella alpina, was expressed in Y. lipolytica. No significant changes occurred in lipid content or fatty acid profiles suggesting that malic enzyme is not the main source of NADPH for lipid accumulation in Y. lipolytica.     

Comprehensive Metabolomic,Lipidomic and Microscopic Profiling of Yarrowia lipolytica during Lipid Accumulation Identifies Targets for Increased Lipogenesis

Yarrowia lipolytica is an oleaginous ascomycete yeast that accumulates large amounts of lipids and has potential as a biofuel producing organism. Despite a growing scientific literature focused on lipid production by Y. lipolytica, there remain significant knowledge gaps regarding the key biological processes involved. We applied a combination of metabolomic and lipidomic profiling approaches as well as microscopic techniques to identify and characterize the key pathways involved in de novo lipid accumulation from glucose in batch cultured, wild-type Y. lipolytica. We found that lipids accumulated rapidly and peaked at 48 hours during the five day experiment, concurrent with a shift in amino acid metabolism. We also report that exhaustion of extracellular sugars coincided with thickening of the cell wall, suggesting that genes involved in cell wall biogenesis may be a useful target for improving the efficiency of lipid producing yeast strains.     

Characterization of a lipid droplet protein from Yarrowia lipolytica that is required for its oleaginous phenotype

Oleaginous microorganisms are characterized by their ability to store high amounts of triacylglycerol (TAG) in intracellular lipid droplets (LDs). In this work, we characterized a protein of the oleaginous yeast Yarrowia lipolytica that is associated with LD and plays a role in the regulation of TAG storage. This protein is required for the oleaginous phenotype of Y. lipolytica because deletion of the coding gene results in a strongly reduced TAG content of the mutant. Therefore, we named it Oleaginicity Inducing LD protein, Oil1. Furthermore, a mutant overexpressing OIL1 accumulates more TAG than the wild type and is delayed in TAG lipolysis when this process is stimulated. We found that Oil1p plays a role in protecting the TAG content of the LD from degradation through lipases under conditions where the cell aims at building up its TAG reserves. Heterologous expression studies showed that Oil1p rescued the phenotype of a Saccharomyces cerevisiae mutant deleted for the perilipin-like protein Pln1p and that its expression in COS-7 cells resulted in increased TAG accumulation, similar to the phenotype of a perilipin 1 expressing control strain. Despite this phenotypical parallels to mammalian perilipins, Oil1p is not a member of this protein family and its activity does not depend on phosphorylation. Rather, our results suggest that ubiquitination might contribute to the function of Oil1p in Y. lipolytica and that a different mechanism evolved in this species to regulate TAG homeostasis.     

Control of Lipid Accumulation in the Yeast Yarrowia lipolytica

A genomic comparison of Yarrowia lipolytica and Saccharomyces cerevisiae indicates that the metabolism of Y. lipolytica is oriented toward the glycerol pathway. To redirect carbon flux toward lipid synthesis, the GUT2 gene, which codes for the glycerol-3-phosphate dehydrogenase isomer, was deleted in Y. lipolytica in this study. This Δgut2 mutant strain demonstrated a threefold increase in lipid accumulation compared to the wild-type strain. However, mobilization of lipid reserves occurred after the exit from the exponential phase due to β-oxidation. Y. lipolytica contains six acyl-coenzyme A oxidases (Aox), encoded by the POX1 to POX6 genes, that catalyze the limiting step of peroxisomal β-oxidation. Additional deletion of the POX1 to POX6 genes in the Δgut2 strain led to a fourfold increase in lipid content. The lipid composition of all of the strains tested demonstrated high proportions of FFA. The size and number of the lipid bodies in these strains were shown to be dependent on the lipid composition and accumulation ratio.     

Yarrowia lipolytica: A model and a tool to understand the mechanisms implicated in lipid accumulation

The oleaginous yeast Yarrowia lipolytica is known to inhabit various lipid-containing environments. One of the most striking features in this yeast is the presence of several multigene families involved in the metabolic pathways of hydrophobic substrate utilization. The complexity and the multiplicity of these genes give Y. lipolytica a wide capability range towards hydrophobic substrate (HS) utilization and storage. The combination of the increasing knowledge of this yeast's metabolism and the development of more efficient genetic tools is offering new perspectives in using Y. lipolytica as a model organism to study the mechanisms involved in lipid metabolism associated to fat uptake, storage, deposition, mobilization and regulation. Nutrient status and culture conditions seem to play a major role in obesity.     

Identification and characterization of DGA2, an acyltransferase of the DGAT1 acyl-CoA:diacylglycerol acyltransferase family in the oleaginous yeast Yarrowia lipolytica. New insights into the storage lipid metabolism of oleaginous yeasts

Triacylglycerols (TAG) and steryl esters (SE) are the principal storage lipids in all eukaryotic cells. In yeasts, these storage lipids accumulate within special organelles known as lipid bodies (LB). In the lipid accumulation-oriented metabolism of the oleaginous yeast Yarrowia lipolytica, storage lipids are mostly found in the form of TAG, and only small amounts of SE accumulate. We report here the identification of a new DAG acyltransferase gene, DGA2, homologous to the ARE genes of Saccharomyces cerevisiae. This gene encodes a member of the type 1 acyl-CoA:diacylglycerol acyltransferase family (DGAT1), which has not previously been identified in yeasts, but is commonly found in mammals and plants. Unlike the Are proteins in S. cerevisiae, Dga2p makes a major contribution to TAG synthesis via an acyl-CoA-dependent mechanism and is not involved in SE synthesis. This enzyme appears to affect the size and morphology of LB, suggesting a direct role of storage lipid proteins in LB formation. We report that the Are1p of Y. lipolytica was essential for sterol esterification, as deletion of the encoding gene (ARE1) completely abolished SE synthesis. Unlike its homologs in yeasts, YlARE1 has no DAG acyltransferase activity. We also reconsider the role and function of all four acyltransferase enzymes involved in the final step of neutral lipid synthesis in this oleaginous yeast.     

Hexokinase—A limiting factor in lipid production from fructose in Yarrowia lipolytica

Microbial biolipid production has become an important part of making biofuel production economically feasible. Genetic engineering has been used to improve the ability of Yarrowia lipolytica, an oleaginous yeast, to produce lipids using glucose-based media. However, few studies have examined lipid accumulation by Y. lipolytica׳s ability to utilize other hexose sugars, and as of yet, the rate-limiting steps in this process are unidentified. In this study, we investigated the de novo accumulation of lipids by Y. lipolytica when grown in glucose, fructose, and sucrose. Three Y. lipolytica wild-type (WT) strains of varied origin differed significantly in their lipid production, growth, and fructose utilization. Hexokinase (ylHXK1p) activity partially explained these differences. Overexpression of the ylHXK1 gene led to increased hexokinase activity (6.5–12 times higher) in the mutants versus the WT strains; a pronounced reduction in cell filamentation in mutants grown in fructose-based media; and improved biomass production, particularly in the mutant whose parent had shown the lowest growth capacity in fructose (French strain W29). All mutants showed improved lipid yield and production when grown on fructose, although the effect was strain dependent (23–55% improvement). Finally, we overexpressed ylHXK1 in a highly modified strain of Y. lipolytica W29 engineered to optimize oil production. This modification was combined with Saccharomyces cerevisiae invertase gene expression to evaluate the resulting mutant׳s ability to produce lipids using cheap industrial substrates, namely sucrose (a major component of molasses). Sucrose turned out to be a better substrate than either of its building blocks, glucose or fructose. Over its 96 h of growth in the bioreactors, this highly modified strain produced 9.15 g L⁻¹ of lipids, yielding 0.262 g g⁻¹ of biomass.     

Engineering polyhydroxyalkanoate content and monomer composition in the oleaginous yeast <Emphasis Type="Italic">Yarrowia lipolytica</Emphasis> by modifying the ß-oxidation multifunctional protein

Recombinant strains of the oleaginous yeast Yarrowia lipolytica expressing the PHA synthase gene (PhaC) from Pseudomonas aeruginosa in the peroxisome were found able to produce polyhydroxyalkanoates (PHA). PHA production yield, but not the monomer composition, was dependent on POX genotype (POX genes encoding acyl-CoA oxidases) (Haddouche et al. FEMS Yeast Res 10:917–927, 2010). In this study of variants of the Y. lipolytica β-oxidation multifunctional enzyme, with deletions or inactivations of the R-3-hydroxyacyl-CoA dehydrogenase domain, we were able to produce hetero-polymers (functional MFE enzyme) or homo-polymers (with no 3-hydroxyacyl-CoA dehydrogenase activity) of PHA consisting principally of 3-hydroxyacid monomers (>80%) of the same length as the external fatty acid used for growth. The redirection of fatty acid flux towards β-oxidation, by deletion of the neutral lipid synthesis pathway (mutant strain Q4 devoid of the acyltransferases encoded by the LRO1, DGA1, DGA2 and ARE1 genes), in combination with variant expressing only the enoyl-CoA hydratase 2 domain, led to a significant increase in PHA levels, to 7.3% of cell dry weight. Finally, the presence of shorter monomers (up to 20% of the monomers) in a mutant strain lacking the peroxisomal 3-hydroxyacyl-CoA dehydrogenase domain provided evidence for the occurrence of partial mitochondrial β-oxidation in Y. lipolytica.     

Reconstruction and In Silico Analysis of Metabolic Network for an Oleaginous Yeast,Yarrowia lipolytica

With the emergence of energy scarcity, the use of renewable energy sources such as biodiesel is becoming increasingly necessary. Recently, many researchers have focused their minds on Yarrowia lipolytica, a model oleaginous yeast, which can be employed to accumulate large amounts of lipids that could be further converted to biodiesel. In order to understand the metabolic characteristics of Y. lipolytica at a systems level and to examine the potential for enhanced lipid production, a genome-scale compartmentalized metabolic network was reconstructed based on a combination of genome annotation and the detailed biochemical knowledge from multiple databases such as KEGG, ENZYME and BIGG. The information about protein and reaction associations of all the organisms in KEGG and Expasy-ENZYME database was arranged into an EXCEL file that can then be regarded as a new useful database to generate other reconstructions. The generated model iYL619_PCP accounts for 619 genes, 843 metabolites and 1,142 reactions including 236 transport reactions, 125 exchange reactions and 13 spontaneous reactions. The in silico model successfully predicted the minimal media and the growing abilities on different substrates. With flux balance analysis, single gene knockouts were also simulated to predict the essential genes and partially essential genes. In addition, flux variability analysis was applied to design new mutant strains that will redirect fluxes through the network and may enhance the production of lipid. This genome-scale metabolic model of Y. lipolytica can facilitate system-level metabolic analysis as well as strain development for improving the production of biodiesels and other valuable products by Y. lipolytica and other closely related oleaginous yeasts.     

Improving CRISPR/Cas9-mediated genome editing efficiency in Yarrowia lipolytica using direct tRNA-sgRNA fusions

Yarrowia lipolytica is an important oleaginous yeast currently used in the production of specialty chemicals and has a great potential for further applications in lipid biotechnology. Harnessing the full potential of Y. lipolytica is, however, limited by its inherent recalcitrance to genetic manipulation. In contrast to Saccharomyces cerevisiae, Y. lipolytica is poor in homology-mediated DNA repair and thus in homologous recombination, which limits site-specific gene editing in this yeast. Recently developed CRISPR/Cas9-based methods using tRNA-sgRNA fusions succeeded in editing some genomic loci in Y. lipolytica. Nonetheless, the majority of other tested loci either failed editing or editing was achieved but at very low efficiency using these methods. Using tools of secondary RNA structure prediction, we were able to improve the design of the tRNA-sgRNA fusions used for the expression of single guide RNA (sgRNA) in such methods. This resulted in high efficiency CRISPR/cas9 gene editing at chromosomal loci that failed gene editing or were edited at very low efficiencies with previous methods. In addition, we characterized the gene editing performance of our newly designed tRNA-sgRNA fusions for both chromosomal gene integration and deletion. As such, this study presents an efficient CRISPR/Cas9-mediated gene-editing tool for efficient genetic engineering of Yarrowia lipolytica.     

Yarrowia lipolytica AAL genes are involved in peroxisomal fatty acid activation

In yeast, β-oxidation of fatty acids (FAs) essentially takes place in peroxisomes, and FA activation must precede FA oxidation. In Saccharomyces cerevisiae, a single fatty-acyl–CoA-synthetase, ScFaa2p, mediates peroxisomal FA activation. We have previously shown that this reaction also exists in the oleaginous yeast Yarrowia lipolytica; however, the protein involved in this process remains unknown. Here, we found that proteins, named Aal proteins (Acyl/Aryl-CoA-ligases), resembling the 4-coumarate–CoA-ligase-like enzymes found in plants are involved in peroxisomal FA activation in Y. lipolytica; Y. lipolytica has 10 AAL genes, eight of which are upregulated by oleate. All the Aal proteins contain a PTS1-type peroxisomal targeting sequence (A/SKL), suggesting a peroxisomal localization. The function of the Aal proteins was analyzed using the faa1Δant1Δ mutant strain, which demonstrates neither cytoplasmic FA activation (direct result of FAA1 deletion) nor peroxisomal FA activation (indirect result of ANT1 deletion, a gene coding an ATP transporter). This strain is thus highly sensitive to external FA levels and unable to store external FAs in lipid bodies (LBs). Whereas the overexpression of (cytoplasmic) AAL1ΔPTS1 was able to partially complement the growth defect observed in the faa1Δant1Δ mutant on short-, medium- and long-chain FA media, the presence of Aal2p to Aal10p only allowed growth on the short-chain FA medium. Additionally, partial LB formation was observed in the oleate medium for strains overexpressing Aal1ΔPTS1p, Aal4ΔPTS1p, Aal7ΔPTS1p, and Aal8ΔPTS1p. Finally, an analysis of the FA content of cells grown in the oleate medium suggested that Aal4p and Aal6p present substrate specificity for C16:1 and/or C18:0.     

Combination of cell disruption technologies for lipid recovery from dry and wet biomass of Yarrowia lipolytica and using green solvents

This work aims to investigate the efficiency of bead milling combined or not with high-pressure homogenization on the cell disruption and lipid recovery from Yarrowia lipolytica oleaginous yeast. First, a simulation study involving the use of the Hansen solubility parameters’ approach was performed in order to identify, among 41 conventional and “green” solvents, the most promising ones that are able to replace n-hexane for lipid recovery from Y. lipolytica biomass. The results obtained showed that the pre-treatment involving both high-pressure homogenization and bead milling applied sequentially was more performant than that involving bead milling alone. In addition, bead-milling parameters were optimized showing an optimal bead size of 4.9 mm and a processing time of 30 s. Among the tested solvents, isoamyl acetate was selected as the most appropriate “green” solvent, maximizing the lipid extraction, compared to n-hexane. Despite the better performance of the dry route compared to the wet one, promising results were obtained towards 1) minimizing the energy consumed and 2) replacing n-hexane by “green” solvents for lipid recovery from Y. lipolytica yeast.     

Genome-scale model-driven strain design for dicarboxylic acid production in <Emphasis Type="Italic">Yarrowia lipolytica</Emphasis>

Background

Recently, there have been several attempts to produce long-chain dicarboxylic acids (DCAs) in various microbial hosts. Of these, Yarrowia lipolytica has great potential due to its oleaginous characteristics and unique ability to utilize hydrophobic substrates. However, Y. lipolytica should be further engineered to make it more competitive: the current approaches are mostly intuitive and cumbersome, thus limiting its industrial application.

Results

In this study, we proposed model-guided metabolic engineering strategies for enhanced production of DCAs in Y. lipolytica. At the outset, we reconstructed genome-scale metabolic model (GSMM) of Y. lipolytica (iYLI647) by substantially expanding the previous models. Subsequently, the model was validated using three sets of published culture experiment data. It was finally exploited to identify genetic engineering targets for overexpression, knockout, and cofactor modification by applying several in silico strain design methods, which potentially give rise to high yield production of the industrially relevant long-chain DCAs, e.g., dodecanedioic acid (DDDA). The resultant targets include (1) malate dehydrogenase and malic enzyme genes and (2) glutamate dehydrogenase gene, in silico overexpression of which generated additional NADPH required for fatty acid synthesis, leading to the increased DDDA fluxes by 48% and 22% higher, respectively, compared to wild-type. We further investigated the effect of supplying branched-chain amino acids on the acetyl-CoA turn-over rate which is key metabolite for fatty acid synthesis, suggesting their significance for production of DDDA in Y. lipolytica.

Conclusion

In silico model-based strain design strategies allowed us to identify several metabolic engineering targets for overproducing DCAs in lipid accumulating yeast, Y. lipolytica. Thus, the current study can provide a methodological framework that is applicable to other oleaginous yeasts for value-added biochemical production.

     

Engineering 4-coumaroyl-CoA derived polyketide production in Yarrowia lipolytica through a β-oxidation mediated strategy

Polyketides are a diverse class of molecules sought after for their valuable properties, including as potential pharmaceuticals. Previously, we demonstrated that the oleaginous yeast Yarrowia lipolytica is an optimal host for production of the simple polyketide, triacetic acid lactone (TAL). We here expand the capacities of this host by overcoming previous media challenges and enabling production of more complex polyketides. Specifically, we employ a β-oxidation related strategy to improve polyketide production directly from defined media. Beyond TAL production, we establish biosynthesis of the 4-coumaroyl-CoA derived polyketides: naringenin, resveratrol, and bisdemethoxycurcumin, as well as the diketide intermediate, (E)-5-(4-hydroxyphenyl)-3-oxopent-4-enoic acid. In this background, we enable high-level de novo production of naringenin through import of both a heterologous pathway and a mutant Y. lipolytica allele. In doing so, we generated an averaged maximum titer of 898 mg/L naringenin, the highest titer reported to date in any host. These results demonstrate that Y. lipolytica is an ideal polyketide production host for more complex 4-coumaroyl-CoA derived products.