Regulation of amino-acid metabolism controls flux to lipid accumulation in Yarrowia lipolytica

Engineering the oleaginous yeast Yarrowia lipolytica for high-level resveratrol production

Resveratrol is a plant secondary metabolite with multiple health-beneficial properties. Microbial production of resveratrol in model microorganisms requires extensive engineering to reach commercially viable levels. Here, we explored the potential of the non-conventional yeast Yarrowia lipolytica to produce resveratrol and several other shikimate pathway-derived metabolites (p-coumaric acid, cis,cis-muconic acid, and salicylic acid). The Y. lipolytica strain expressing a heterologous pathway produced 52.1 ± 1.2 mg/L resveratrol in a small-scale cultivation. The titer increased to 409.0 ± 1.2 mg/L when the strain was further engineered with feedback-insensitive alleles of the key genes in the shikimate pathway and with five additional copies of the heterologous biosynthetic genes. In controlled fed-batch bioreactor, the strain produced 12.4 ± 0.3 g/L resveratrol, the highest reported titer to date for de novo resveratrol production, with a yield on glucose of 54.4 ± 1.6 mg/g and a productivity of 0.14 ± 0.01 g/L/h. The study showed that Y. lipolytica is an attractive host organism for the production of resveratrol and possibly other shikimate-pathway derived metabolites.

•

Oleaginous yeast Y. lipolytica was engineered for production of aromatic compounds.

•

High resveratrol production required increased activities of Aro4p and Aro7p.

•

Multiple integration of resveratrol biosynthetic genes improved production.

•

Fed-batch fermentation enabled de novo production of 12.4 g/L resveratrol.

Synthetic microbial communities of heterotrophs and phototrophs facilitate sustainable growth

Microbial communities comprised of phototrophs and heterotrophs hold great promise for sustainable biotechnology. Successful application of these communities relies on the selection of appropriate partners. Here we construct four community metabolic models to guide strain selection, pairing phototrophic, sucrose-secreting Synechococcus elongatus with heterotrophic Escherichia coli K-12, Escherichia coli W, Yarrowia lipolytica, or Bacillus subtilis. Model simulations reveae metabolic exchanges that sustain the heterotrophs in minimal media devoid of any organic carbon source, pointing to S. elongatus-E. coli K-12 as the most active community. Experimental validation of flux predictions for this pair confirms metabolic interactions and potential production capabilities. Synthetic communities bypass member-specific metabolic bottlenecks (e.g. histidine- and transport-related reactions) and compensate for lethal genetic traits, achieving up to 27% recovery from lethal knockouts. The study provides a robust modelling framework for the rational design of synthetic communities with optimized growth sustainability using phototrophic partners.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Synthetic biology, systems biology, and metabolic engineering of Yarrowia lipolytica toward a sustainable biorefinery platform

Yarrowia lipolytica is an oleaginous yeast that has been substantially engineered for production of oleochemicals and drop-in transportation fuels. The unique acetyl-CoA/malonyl-CoA supply mode along with the versatile carbon-utilization pathways makes this yeast a superior host to upgrade low-value carbons into high-value secondary metabolites and fatty acid-based chemicals. The expanded synthetic biology toolkits enabled us to explore a large portfolio of specialized metabolism beyond fatty acids and lipid-based chemicals. In this review, we will summarize the recent advances in genetic, omics, and computational tool development that enables us to streamline the genetic or genomic modification for Y. lipolytica. We will also summarize various metabolic engineering strategies to harness the endogenous acetyl-CoA/malonyl-CoA/HMG-CoA pathway for production of complex oleochemicals, polyols, terpenes, polyketides, and commodity chemicals. We envision that Y. lipolytica will be an excellent microbial chassis to expand nature's biosynthetic capacity to produce plant secondary metabolites, industrially relevant oleochemicals, agrochemicals, commodity, and specialty chemicals and empower us to build a sustainable biorefinery platform that contributes to the prosperity of a bio-based economy in the future.

Effects of a high-cultivation temperature on the physiology of three different Yarrowia lipolytica strains

Despite the increasing relevance, ranging from academic research to industrial applications, only a limited number of non-conventional, oleaginous Yarrowia lipolytica strains are characterized in detail. Therefore, we analyzed three strains in regard to their metabolic and physiological properties, especially with respect to important characteristics of a production strain. By investigating different cultivation conditions and media compositions, similarities and differences between the distinct strain backgrounds could be derived. Especially sugar alcohol production, as well as an agglomeration of cells were found to be connected with growth at high temperatures. In addition, sugar alcohol production was independent of high substrate concentrations under these conditions. To investigate the genotypic basis of particular traits, including growth characteristics and metabolite concentrations, genomic analysis were performed. We found sequence variations for one third of the annotated proteins but no obvious link to all phenotypic features.

Analysis of human metabolism by reducing the complexity of the genome-scale models using redHUMAN

Altered metabolism is associated with many human diseases. Human genome-scale metabolic models (GEMs) were reconstructed within systems biology to study the biochemistry occurring in human cells. However, the complexity of these networks hinders a consistent and concise physiological representation. We present here redHUMAN, a workflow for reconstructing reduced models that focus on parts of the metabolism relevant to a specific physiology using the recently established methods redGEM and lumpGEM. The reductions include the thermodynamic properties of compounds and reactions guaranteeing the consistency of predictions with the bioenergetics of the cell. We introduce a method (redGEMX) to incorporate the pathways used by cells to adapt to the medium. We provide the thermodynamic curation of the human GEMs Recon2 and Recon3D and we apply the redHUMAN workflow to derive leukemia-specific reduced models. The reduced models are powerful platforms for studying metabolic differences between phenotypes, such as diseased and healthy cells.

Analysis of Yarrowia lipolytica growth, catabolism, and terpenoid biosynthesis during utilization of lipid-derived feedstock

This study employs biomass growth analyses and ¹³C-isotope tracing to investigate lipid feedstock utilization by Yarrowia lipolytica. Compared to glucose, oil-feedstock in the minimal medium increases the yeast's biomass yields and cell sizes, but decreases its protein content (<20% of total biomass) and enzyme abundances for product synthesis. Labeling results indicate a segregated metabolic network (the glycolysis vs. the TCA cycle) during co-catabolism of sugars (glucose or glycerol) with fatty acid substrates, which facilitates resource allocations for biosynthesis without catabolite repressions. This study has also examined the performance of a β-carotene producing strain in different growth mediums. Canola oil-containing yeast-peptone (YP) has resulted in the best β-carotene titer (121 ± 13 mg/L), two-fold higher than the glucose based YP medium. These results highlight the potential of Y. lipolytica for the valorization of waste-derived lipid feedstock.

•

¹³C tracing was used to track Y. lipolytica metabolism of lipid-based feedstock.

•

Y. lipolytica has a segregated flux network for lipid and sugar co-utilizations.

•

Lipid feedstock and nitrogen sources affect cell morphology and optical density.

•

Lipid feedstock benefits both Y. lipolytica growth and carotenoid biosynthesis.

Dynamic resource allocation drives growth under nitrogen starvation in eukaryotes

Cells can sense changes in their extracellular environment and subsequently adapt their biomass composition. Nutrient abundance defines the capability of the cell to produce biomass components. Under nutrient-limited conditions, resource allocation dramatically shifts to carbon-rich molecules. Here, we used dynamic biomass composition data to predict changes in growth and reaction flux distributions using the available genome-scale metabolic models of five eukaryotic organisms (three heterotrophs and two phototrophs). We identified temporal profiles of metabolic fluxes that indicate long-term trends in pathway and organelle function in response to nitrogen depletion. Surprisingly, our calculations of model sensitivity and biosynthetic cost showed that free energy of biomass metabolites is the main driver of biosynthetic cost and not molecular weight, thus explaining the high costs of arginine and histidine. We demonstrated how metabolic models can accurately predict the complexity of interwoven mechanisms in response to stress over the course of growth.

C/N ratio and carbon source-dependent lipid production profiling in Rhodotorula toruloides

Microbial oils are lipids produced by oleaginous microorganisms, which can be used as a potential feedstock for oleochemical production. The oleaginous yeast Rhodotorula toruloides can co-produce microbial oils and high-value compounds from low-cost substrates, such as xylose and acetic acid (from hemicellulosic hydrolysates) and raw glycerol (a byproduct of biodiesel production). One step towards economic viability is identifying the best conditions for lipid production, primarily the most suitable carbon-to-nitrogen ratio (C/N). Here, we aimed to identify the best conditions and cultivation mode for lipid production by R. toruloides using various low-cost substrates and a range of C/N ratios (60, 80, 100, and 120). Turbidostat mode was used to achieve a steady state at the maximal specific growth rate and to avoid continuously changing environmental conditions (i.e., C/N ratio) that inherently occur in batch mode. Regardless of the carbon source, higher C/N ratios increased lipid yields (up to 60% on xylose at a C/N of 120) but decreased the specific growth rate. Growth on glycerol resulted in the highest specific growth and lipid production (0.085 g lipids/gDW*h) rates at C/Ns between 60 and 100. We went on to study lipid production using glycerol in both batch and fed-batch modes, which resulted in lower specific lipid production rates compared with turbisdostat, however, fed batch is superior in terms of biomass production and lipid titers. By combining the data we obtained in these experiments with a genome-scale metabolic model of R. toruloides, we identified targets for improvements in lipid production that could be carried out either by metabolic engineering or process optimization.

Application of Stable Isotope Tracing to Elucidate Metabolic Dynamics During Yarrowia lipolytica α-Ionone Fermentation

Targeted metabolite analysis in combination with ¹³C-tracing is a convenient strategy to determine pathway activity in biological systems; however, metabolite analysis is limited by challenges in separating and detecting pathway intermediates with current chromatographic methods. Here, a hydrophilic interaction chromatography tandem mass spectrometry approach was developed for improved metabolite separation, isotopologue analysis, and quantification. The physiological responses of a Yarrowia lipolytica strain engineered to produce ∼400 mg/L α-ionone and temporal changes in metabolism were quantified (e.g., mevalonate secretion, then uptake) indicating bottleneck shifts in the engineered pathway over the course of fermentation. Dynamic labeling results indicated limited tricarboxylic acid cycle label incorporation and, combined with a measurable ATP shortage during the high ionone production phase, suggested that electron transport and oxidative phosphorylation may limit energy supply and strain performance. The results provide insights into terpenoid pathway metabolic dynamics of non-model yeasts and offer guidelines for sensor development and modular engineering.

•

A HILIC method is demonstrated for efficient separation of 57 cellular metabolites

•

Production of α-ionone was ∼400 mg/L in bench-top bioreactors

•

Engineered Y. lipolytica secreted then consumed mevalonate during fermentation

•

Oxidative phosphorylation may limit performance in high-cell-density fermentations

Induction of triacylglycerol synthesis in yeast by cell cycle arrest

In this study, we found that cell cycle arrest induced by alpha-factor mating pheromone (G1), hydroxyurea (S) or nocodazole (G2/M) was associated to increased lipid droplet (LD) content. To identify novel cell cycle genes involved in LD homeostasis, we screened a deletion library for strains with altered LD levels. Among the mutants related to mitotic cell cycle, we found 24 hits that displayed a significantly higher LD content. Ontology mapping showed that neither a biological process nor a specific cell cycle phase was enriched among the hits. We decided to further study the role of SWI4 on LD homeostasis as it is involved in G1/S transition, a stage where lipolysis is active. The high LD content of swi4Δ mutant was not due to inhibition of lipolysis, but due to an increase in triacylglycerol (TAG) synthesis. In addition, deletion of the AMP kinase gene SNF1 or inhibition of TORC1 activity, both known regulators of LD homeostasis, further increased the LD content of a swi4Δ mutant. These findings highlight a role of the cell cycle regulator SWI4 in the coordination of lipid metabolism which is independent of the TORC1 and SNF1/AMPK pathways.

A comprehensive genome-scale model for Rhodosporidium toruloides IFO0880 accounting for functional genomics and phenotypic data

Rhodosporidium toruloides is a red, basidiomycetes yeast that can accumulate a large amount of lipids and produce carotenoids. To better assess this non-model yeast's metabolic capabilities, we reconstructed a genome-scale model of R. toruloides IFO0880's metabolic network (iRhto1108) accounting for 2204 reactions, 1985 metabolites and 1108 genes. In this work, we integrated and supplemented the current knowledge with in-house generated biomass composition and experimental measurements pertaining to the organism's metabolic capabilities. Predictions of genotype-phenotype relations were improved through manual curation of gene-protein-reaction rules for 543 reactions leading to correct recapitulations of 84.5% of gene essentiality data (sensitivity of 94.3% and specificity of 53.8%). Organism-specific macromolecular composition and ATP maintenance requirements were experimentally measured for two separate growth conditions: (i) carbon and (ii) nitrogen limitations. Overall, iRhto1108 reproduced R. toruloides's utilization capabilities for 18 alternate substrates, matched measured wild-type growth yield, and recapitulated the viability of 772 out of 819 deletion mutants. As a demonstration to the model's fidelity in guiding engineering interventions, the OptForce procedure was applied on iRhto1108 for triacylglycerol overproduction. Suggested interventions recapitulated many of the previous successful implementations of genetic modifications and put forth a few new ones.

Remodeling of Chlamydomonas Metabolism Using Synthetic Inducers Results in Lipid Storage during Growth

Microalgae accumulate lipids during stress such as that of nutrient deprivation, concomitant with cessation of growth and depletion of chloroplasts. By contrast, certain small chemical compounds selected by high-throughput screening in Chlamydomonas reinhardtii can induce lipid accumulation during growth, maintaining biomass. Comprehensive pathway analyses using proteomics, transcriptomics, and metabolomics data were acquired from Chlamydomonas cells grown in the presence of one of two structurally distinct lipid activators. WD10784 stimulates both starch and lipid accumulation, whereas WD30030-treated cells accumulate only lipids. The differences in starch accumulation are largely due to differential effects of the two compounds on substrate levels that feed into starch synthesis and on genes encoding starch metabolic enzymes. The compounds had differential effects on photosynthesis, respiration, and oxidative stress pathways. Cells treated with WD10784 showed slowed growth over time and reduced abundance of photosynthetic proteins, decreased respiration, and increased oxidative stress proteins, glutathione, and reactive oxygen species specific to this compound. Both compounds maintained central carbon and nitrogen metabolism, including the tricarboxylic acid cycle, glycolysis, respiration, and the Calvin-Benson-Bassham cycle. There were few changes in proteins and transcripts related to fatty acid biosynthesis, whereas proteins and transcripts for triglyceride production were elevated, suggesting that lipid synthesis is largely driven by substrate availability. This study reports that the compound WD30030 and, to a lesser extent WD10784, increases lipid and lipid droplet synthesis and storage without restricting growth or biomass accumulation by mechanisms that are substantially different from nutrient deprivation.

Multi-Omics Analysis of Fatty Alcohol Production in Engineered Yeasts Saccharomyces cerevisiae and Yarrowia lipolytica

Fatty alcohols are widely used in various applications within a diverse set of industries, such as the soap and detergent industry, the personal care, and cosmetics industry, as well as the food industry. The total world production of fatty alcohols is over 2 million tons with approximately equal parts derived from fossil oil and from plant oils or animal fats. Due to the environmental impact of these production methods, there is an interest in alternative methods for fatty alcohol production via microbial fermentation using cheap renewable feedstocks. In this study, we aimed to obtain a better understanding of how fatty alcohol biosynthesis impacts the host organism, baker’s yeast Saccharomyces cerevisiae or oleaginous yeast Yarrowia lipolytica. Producing and non-producing strains were compared in growth and nitrogen-depletion cultivation phases. The multi-omics analysis included physiological characterization, transcriptome analysis by RNAseq, ¹³Cmetabolic flux analysis, and intracellular metabolomics. Both species accumulated fatty alcohols under nitrogen-depletion conditions but not during growth. The fatty alcohol–producing Y. lipolytica strain had a higher fatty alcohol production rate than an analogous S. cerevisiae strain. Nitrogen-depletion phase was associated with lower glucose uptake rates and a decrease in the intracellular concentration of acetyl–CoA in both yeast species, as well as increased organic acid secretion rates in Y. lipolytica. Expression of the fatty alcohol–producing enzyme fatty acyl–CoA reductase alleviated the growth defect caused by deletion of hexadecenal dehydrogenase encoding genes (HFD1 and HFD4) in Y. lipolytica. RNAseq analysis showed that fatty alcohol production triggered a cell wall stress response in S. cerevisiae. RNAseq analysis also showed that both nitrogen-depletion and fatty alcohol production have substantial effects on the expression of transporter encoding genes in Y. lipolytica. In conclusion, through this multi-omics study, we uncovered some effects of fatty alcohol production on the host metabolism. This knowledge can be used as guidance for further strain improvement towards the production of fatty alcohols.

Dimorphism of Trichosporon cutaneum and impact on its lipid production

BACKGROUND

Compared to the oleaginous yeast Yarrowia lipolytica, Trichosporon cutaneum can metabolize pentose sugars more efficiently, and in the meantime is more tolerant to inhibitors, which is suitable for lipid production from lignocellulosic biomass. However, this species experiences dimorphic transition between yeast-form cells and hyphae during submerged fermentation, which consequently affects the rheology and mass transfer performance of the fermentation broth and its lipid production.

RESULTS

The strain T. cutaneum B3 was cultured with medium composed of yeast extract, glucose and basic minerals. The experimental results indicated that yeast-form morphology was developed when yeast extract was supplemented at 1 g/L, but hyphae were observed when yeast extract supplementation was increased to 3 g/L and 5 g/L, respectively. We speculated that difference in nitrogen supply to the medium might be a major reason for the dimorphic transition, which was confirmed by the culture with media supplemented with yeast extract at 1 g/L and urea at 0.5 g/L and 1.0 g/L to maintain total nitrogen at same levels as that detected in the media with yeast extract supplemented at 3 g/L and 5 g/L. The morphological change of T. cutaneum B3 affected not only the content of intracellular lipids but also their composition, due to its impact on the rheology and oxygen mass transfer performance of the fermentation broth, and more lipids with less polyunsaturated fatty acids such as linoleic acid (C18:2) were produced by the yeast-form cells. When T. cutaneum B3 was cultured at an aeration rate of 1.5 vvm for 72 h with the medium composed of 60 g/L glucose, 3 g/L yeast extract and basic minerals, 27.1 g (dry cell weight)/L biomass was accumulated with the lipid content of 46.2%, and lipid productivity and yield were calculated to be 0.174 g/L/h and 0.21 g/g, respectively. Comparative transcriptomics analysis identified differently expressed genes for sugar metabolism and lipid synthesis as well as signal transduction for the dimorphic transition of T. cutaneum B3.

CONCLUSIONS

Assimilable nitrogen was validated as one of the major reasons for the dimorphic transition between yeast-form morphology and hyphae with T. cutaneum, and the yeast-form morphology was more suitable for lipid production at high content with less polyunsaturated fatty acids as feedstock for biodiesel production.

Insights into the Metabolism of Oleaginous Rhodococcus spp

Some species belonging to the Rhodococcus genus, such as Rhodococcus opacus, R. jostii, and R. wratislaviensis, are known to be oleaginous microorganisms, since they are able to accumulate triacylglycerols (TAG) at more than 20% of their weight (dry weight). Oleaginous rhodococci are promising microbial cell factories for the production of lipids to be used as fuels and chemicals. Cells could be engineered to create strains capable of producing high quantities of oils from industrial wastes and a variety of high-value lipids. The comprehensive understanding of carbon metabolism and its regulation will contribute to the design of a reliable process for bacterial oil production. Bacterial oleagenicity requires an integral configuration of metabolism and regulatory processes rather than the sole existence of an efficient lipid biosynthesis pathway. In recent years, several studies have been focused on basic aspects of TAG biosynthesis and accumulation using R. opacus PD630 and R. jostii RHA1 strains as models of oleaginous bacteria. The combination of results obtained in these studies allows us to propose a metabolic landscape for oleaginous rhodococci. In this context, this article provides a comprehensive and integrative view of different metabolic and regulatory attributes and innovations that explain the extraordinary ability of these bacteria to synthesize and accumulate TAG. We hope that the accessibility to such information in an integrated way will help researchers to rationally select new targets for further studies in the field.

Role of Adenosine Monophosphate Deaminase during Fatty Acid Accumulation in Oleaginous Fungus Mortierella alpina

In oleaginous micro-organisms, nitrogen limitation activates adenosine monophosphate deaminase (AMPD) and promotes lipogenesis via the inhibition of isocitrate dehydrogenase. We found that the overexpression of homologous AMPD in Mortierella alpina favored lipid synthesis over cell growth. Total fatty acid content in the recombinant strain was 15.0-34.3% higher than that in the control, even though their biomass was similar. During the early fermentation stage, the intracellular AMP level reduced by 40-60%, together with a 1.9-2.7-fold increase in citrate content compared with the control, therefore provided more precursors for fatty acid synthesis. Moreover, the decreased AMP level resulted in metabolic reprogramming, reflected by the blocked TCA cycle and reduction of amino acids, distributing more carbon to lipid synthesis pathways. By coupling the energy balance with lipogenesis, this study provides new insights into cell metabolism under nitrogen-limited conditions and targets the regulation of fatty acid accumulation in oleaginous micro-organisms.

Common aspects in the engineering of yeasts for fatty acid- and isoprene-based products

The biosynthetic pathways for most lipophilic metabolites share several common principles. These substances are built almost exclusively from acetyl-CoA as the donor for the carbon scaffold and NADPH is required for the reductive steps during biosynthesis. Due to their hydrophobicity, the end products are sequestered into the same cellular compartment, the lipid droplet. In this review, we will summarize the efforts in the metabolic engineering of yeasts for the production of two major hydrophobic substance classes, fatty acid-based lipids and isoprenoids, with regard to these common aspects. We will compare and discuss the results of genetic engineering strategies to construct strains with enhanced synthesis of the precursor acetyl-CoA and with modified redox metabolism for improved NADPH supply. We will also discuss the role of the lipid droplet in the storage of the hydrophobic product and review the strategies to either optimize this organelle for higher capacity or to achieve excretion of the product into the medium.

Glutamate dehydrogenases in the oleaginous yeast Yarrowia lipolytica

Glutamate dehydrogenases (GDHs) are fundamental to cellular nitrogen and energy balance. Yet little is known about these enzymes in the oleaginous yeast Yarrowia lipolytica. The YALI0F17820g and YALI0E09603g genes, encoding potential GDH enzymes in this organism, were examined. Heterologous expression in gdh-null Saccharomyces cerevisiae and examination of Y. lipolytica strains carrying gene deletions demonstrate that YALI0F17820g (ylGDH1) encodes a NADP-dependent GDH whereas YALI0E09603g (ylGDH2) encodes a NAD-dependent GDH enzyme. The activity encoded by these two genes accounts for all measurable GDH activity in Y. lipolytica. Levels of the two enzyme activities are comparable during logarithmic growth on rich medium, but the NADP-ylGDH1p enzyme activity is most highly expressed in stationary and nitrogen starved cells by threefold to 12-fold. Replacement of ammonia with glutamate causes a decrease in NADP-ylGdh1p activity, whereas NAD-ylGdh2p activity is increased. When glutamate is both carbon and nitrogen sources, the activity of NAD-ylGDH2p becomes dominant up to 18-fold compared with that of NADP-ylGDH1p. Gene deletion followed by growth on different carbon and nitrogen sources shows that NADP-ylGdh1p is required for efficient nitrogen assimilation whereas NAD-ylGdh2p plays a role in nitrogen and carbon utilization from glutamate. Overexpression experiments demonstrate that ylGDH1 and ylGDH2 are not interchangeable. These studies provide a vital basis for future consideration of how these enzymes function to facilitate energy and nitrogen homeostasis in Y. lipolytica.

Monitoring of the enzymatic activity of intracellular lipases of

Under nitrogen starvation, Ustilago maydis forms lipid droplets (LDs). Although the dynamics of these organelles are known in the literature, the identity of the lipases implicated in their degradation is unknown. We determined lipase activity and identified the intracellular lipases expressed during growth under nitrogen starvation and YPD media by zymograms. The results showed that cytosolic extracts exhibited higher lipase activity when cells were grown in YPD. Under nitrogen starvation, lipase activity was not detected after 24 h of culture, resulting in lipid accumulation in LDs. This suggests that these lipases could be implicated in LD degradation. In the zymogram, two bands, one of 25 and the other of 37 kDa, presented lipase activity. The YPD extracts showed lipase activity in olive and almond oils, which contain triacylglycerols with mono and polyunsaturated fatty acids. This is the first report about U. maydis cytosolic lipases involved in LD degradation.

Transcriptomic analysis of the oleaginous yeast Lipomyces starkeyi during lipid accumulation on enzymatically treated corn stover hydrolysate

BACKGROUND

Efficient and economically viable production of biofuels from lignocellulosic biomass is dependent on mechanical and chemical pretreatment and enzymatic hydrolysis of plant material. These processing steps yield simple sugars as well as plant-derived and process-added organic acids, sugar-derived dehydration products, aldehydes, phenolics and other compounds that inhibit the growth of many microorganisms. Lipomyces starkeyi is an oleaginous yeast capable of robust growth on a variety of sugars and lipid accumulation on pretreated lignocellulosic substrates making it attractive as an industrial producer of biofuels. Here, we examined gene expression during batch growth and lipid accumulation in a 20-L bioreactor with either a blend of pure glucose and xylose or pretreated corn stover (PCS) that had been enzymatically hydrolyzed as the carbon sources.

RESULTS

We monitored sugar and ammonium utilization as well as biomass accumulation and found that growth of L. starkeyi is inhibited with PCS hydrolysate as the carbon source. Both acetic acid and furfural are present at concentrations toxic to L. starkeyi in PCS hydrolysate. We quantified gene expression at seven time-points for each carbon source during batch growth and found that gene expression is similar at physiologically equivalent points. Analysis of promoter regions revealed that gene expression during the transition to lipid accumulation is regulated by carbon and nitrogen catabolite repression, regardless of carbon source and is associated with decreased expression of the translation machinery and suppression of the cell cycle. We identified 73 differentially expressed genes during growth phase in the bioreactor that may be involved in detoxification of corn stover hydrolysate.

CONCLUSIONS

Growth of L. starkeyi is inhibited by compounds present in PCS hydrolysate. Here, we monitored key metabolites to establish physiologically equivalent comparisons during a batch bioreactor run comparing PCS hydrolysate and purified sugars. L. starkeyi's response to PCS hydrolysate is primarily at the beginning of the run during growth phase when inhibitory compounds are presumably at their highest concentration and inducing the general detoxification response by L. starkeyi. Differentially expressed genes identified herein during growth phase will aid in the improvement of industrial strains capable of robust growth on substrates containing various growth inhibitory compounds.

Metabolic traits specific for lipid-overproducing strain of Mucor circinelloides WJ11 identified by genome-scale modeling approach

The genome-scale metabolic model of a lipid-overproducing strain of Mucor circinelloides WJ11 was developed. The model (iNI1159) contained 1,159 genes, 648 EC numbers, 1,537 metabolites, and 1,355 metabolic reactions, which were localized in different compartments of the cell. Using flux balance analysis (FBA), the iNI1159 model was validated by predicting the specific growth rate. The metabolic traits investigated by phenotypic phase plane analysis (PhPP) showed a relationship between the nutrient uptake rate, cell growth, and the triacylglycerol production rate, demonstrating the strength of the model. A putative set of metabolic reactions affecting the lipid-accumulation process was identified when the metabolic flux distributions under nitrogen-limited conditions were altered by performing fast flux variability analysis (fastFVA) and relative flux change. Comparative analysis of the metabolic models of the lipid-overproducing strain WJ11 (iNI1159) and the reference strain CBS277.49 (iWV1213) using both fastFVA and coordinate hit-and-run with rounding (CHRR) showed that the flux distributions between these two models were significantly different. Notably, a higher flux distribution through lipid metabolisms such as lanosterol, zymosterol, glycerolipid and fatty acids biosynthesis in iNI1159 was observed, leading to an increased lipid production when compared to iWV1213. In contrast, iWV1213 exhibited a higher flux distribution across carbohydrate and amino acid metabolisms and thus generated a high flux for biomass production. This study demonstrated that iNI1159 is an effective predictive tool for the pathway engineering of oleaginous strains for the production of diversified oleochemicals with industrial relevance.

Proteome analysis of xylose metabolism in Rhodotorula toruloides during lipid production

BACKGROUND

Rhodotorula toruloides is a promising platform organism for production of lipids from lignocellulosic substrates. Little is known about the metabolic aspects of lipid production from the lignocellolosic sugar xylose by oleaginous yeasts in general and R. toruloides in particular. This study presents the first proteome analysis of the metabolism of R. toruloides during conversion of xylose to lipids.

RESULTS

Rhodotorula toruloides cultivated on either glucose or xylose was subjected to comparative analysis of its growth dynamics, lipid composition, fatty acid profiles and proteome. The maximum growth and sugar uptake rate of glucose-grown R. toruloides cells were almost twice that of xylose-grown cells. Cultivation on xylose medium resulted in a lower final biomass yield although final cellular lipid content was similar between glucose- and xylose-grown cells. Analysis of lipid classes revealed the presence of monoacylglycerol in the early exponential growth phase as well as a high proportion of free fatty acids. Carbon source-specific changes in lipid profiles were only observed at early exponential growth phase, where C18 fatty acids were more saturated in xylose-grown cells. Proteins involved in sugar transport, initial steps of xylose assimilation and NADPH regeneration were among the proteins whose levels increased the most in xylose-grown cells across all time points. The levels of enzymes involved in the mevalonate pathway, phospholipid biosynthesis and amino acids biosynthesis differed in response to carbon source. In addition, xylose-grown cells contained higher levels of enzymes involved in peroxisomal beta-oxidation and oxidative stress response compared to cells cultivated on glucose.

CONCLUSIONS

The results obtained in the present study suggest that sugar import is the limiting step during xylose conversion by R. toruloides into lipids. NADPH appeared to be regenerated primarily through pentose phosphate pathway although it may also involve malic enzyme as well as alcohol and aldehyde dehydrogenases. Increases in enzyme levels of both fatty acid biosynthesis and beta-oxidation in xylose-grown cells was predicted to result in a futile cycle. The results presented here are valuable for the development of lipid production processes employing R. toruloides on xylose-containing substrates.

Cellular and metabolic engineering of oleaginous yeast Yarrowia lipolytica for bioconversion of hydrophobic substrates into high-value products

The non-conventional oleaginous yeast Yarrowia lipolytica is able to utilize both hydrophilic and hydrophobic carbon sources as substrates and convert them into value-added bioproducts such as organic acids, extracellular proteins, wax esters, long-chain diacids, fatty acid ethyl esters, carotenoids and omega-3 fatty acids. Metabolic pathway analysis and previous research results show that hydrophobic substrates are potentially more preferred by Y. lipolytica than hydrophilic substrates to make high-value products at higher productivity, titer, rate, and yield. Hence, Y. lipolytica is becoming an efficient and promising biomanufacturing platform due to its capabilities in biosynthesis of extracellular lipases and directly converting the extracellular triacylglycerol oils and fats into high-value products. It is believed that the cell size and morphology of the Y. lipolytica is related to the cell growth, nutrient uptake, and product formation. Dimorphic Y. lipolytica demonstrates the yeast-to-hypha transition in response to the extracellular environments and genetic background. Yeast-to-hyphal transition regulating genes, such as YlBEM1, YlMHY1 and YlZNC1 and so forth, have been identified to involve as major transcriptional factors that control morphology transition in Y. lipolytica. The connection of the cell polarization including cell cycle and the dimorphic transition with the cell size and morphology in Y. lipolytica adapting to new growth are reviewed and discussed. This review also summarizes the general and advanced genetic tools that are used to build a Y. lipolytica biomanufacturing platform.

Eliciting the impacts of cellular noise on metabolic trade-offs by quantitative mass imaging

Optimal metabolic trade-offs between growth and productivity are key constraints in strain optimization by metabolic engineering; however, how cellular noise impacts these trade-offs and drives the emergence of subpopulations with distinct resource allocation strategies, remains largely unknown. Here, we introduce a single-cell strategy for quantifying the trade-offs between triacylglycerol production and growth in the oleaginous microorganism Yarrowia lipolytica. The strategy relies on high-throughput quantitative-phase imaging and, enabled by nanoscale secondary ion mass spectrometry analyses and dedicated image processing, allows us to image how resources are partitioned between growth and productivity. Enhanced precision over population-averaging biotechnologies and conventional microscopy demonstrates how cellular noise impacts growth and productivity differently. As such, subpopulations with distinct metabolic trade-offs emerge, with notable impacts on strain performance and robustness. By quantifying the self-degradation of cytosolic macromolecules under nutrient-limiting conditions, we discover the cell-to-cell heterogeneity in protein and fatty-acid recycling, unmasking a potential bet-hedging strategy under starvation.

Understanding lipogenesis by dynamically profiling transcriptional activity of lipogenic promoters in Yarrowia lipolytica

Lipogenesis is a complicated process involving global transcriptional reprogramming of lipogenic pathways. It is commonly believed that nitrogen starvation triggers a metabolic shift that reroutes carbon flux from Krebs cycles to lipogenesis. In this study, we systematically surveyed and dynamically profiled the transcriptional activity of 22 lipogenic promoters aiming to delineate a picture how nitrogen starvation regulates lipogenesis in Y. lipolytica. These lipogenic promoters drive the expression of critical pathways that are responsible for the generation of reducing equivalents (NADPH), carbon backbones (acetyl-CoA, malonyl-CoA, DHAP, etc.), synthesis and degradation of fatty acids. Specifically, our investigated promoters span across an array of metabolic pathways, including glycolysis, Krebs cycle, pentose phosphate pathway, mannitol cycle, glutamine-GABA cycle, fatty acid and lipid synthesis, glyoxylate, β-oxidation, and POM (pyruvate-oxaloacetate-malate) cycle. Our work provides evidences that mannitol cycle, glutamine-GABA cycle and amino acid degradation, pyruvate oxidation, and acetate assimilation pathways are lipogenesis-related steps involved in generating cytosolic NADPH and acetyl-CoA precursors. This systematic investigation and dynamic profiling of lipogenic promoters may help us better understand lipogenesis, facilitate the formulation of structure-based kinetic models, as well as develop efficient cell factories for fuels and chemical production in oleaginous species.

SLIMEr: probing flexibility of lipid metabolism in yeast with an improved constraint-based modeling framework

BACKGROUND

A recurrent problem in genome-scale metabolic models (GEMs) is to correctly represent lipids as biomass requirements, due to the numerous of possible combinations of individual lipid species and the corresponding lack of fully detailed data. In this study we present SLIMEr, a formalism for correctly representing lipid requirements in GEMs using commonly available experimental data.

RESULTS

SLIMEr enhances a GEM with mathematical constructs where we Split Lipids Into Measurable Entities (SLIME reactions), in addition to constraints on both the lipid classes and the acyl chain distribution. By implementing SLIMEr on the consensus GEM of Saccharomyces cerevisiae, we can represent accurate amounts of lipid species, analyze the flexibility of the resulting distribution, and compute the energy costs of moving from one metabolic state to another.

CONCLUSIONS

The approach shows potential for better understanding lipid metabolism in yeast under different conditions. SLIMEr is freely available at https://github.com/SysBioChalmers/SLIMEr .

In silico identification of metabolic engineering strategies for improved lipid production in Yarrowia lipolytica by genome-scale metabolic modeling

BACKGROUND

Yarrowia lipolytica, an oleaginous yeast, is a promising platform strain for production of biofuels and oleochemicals as it can accumulate a high level of lipids in response to nitrogen limitation. Accordingly, many metabolic engineering efforts have been made to develop engineered strains of Y. lipolytica with higher lipid yields. Genome-scale model of metabolism (GEM) is a powerful tool for identifying novel genetic designs for metabolic engineering. Several GEMs for Y. lipolytica have recently been developed; however, not many applications of the GEMs have been reported for actual metabolic engineering of Y. lipolytica. The major obstacle impeding the application of Y. lipolytica GEMs is the lack of proper methods for predicting phenotypes of the cells in the nitrogen-limited condition, or more specifically in the stationary phase of a batch culture.

RESULTS

In this study, we showed that environmental version of minimization of metabolic adjustment (eMOMA) can be used for predicting metabolic flux distribution of Y. lipolytica under the nitrogen-limited condition and identifying metabolic engineering strategies to improve lipid production in Y. lipolytica. Several well-characterized overexpression targets, such as diglyceride acyltransferase, acetyl-CoA carboxylase, and stearoyl-CoA desaturase, were successfully rediscovered by our eMOMA-based design method, showing the relevance of prediction results. Interestingly, the eMOMA-based design method also suggested non-intuitive knockout targets, and we experimentally validated the prediction with a mutant lacking YALI0F30745g, one of the predicted targets involved in one-carbon/methionine metabolism. The mutant accumulated 45% more lipids compared to the wild-type.

CONCLUSION

This study demonstrated that eMOMA is a powerful computational method for understanding and engineering the metabolism of Y. lipolytica and potentially other oleaginous microorganisms.

Yarrowia lipolytica: a beneficious yeast in biotechnology as a rare opportunistic fungal pathogen: a minireview

Yarrowia lipolytica is one of the most studied "non-conventional" yeast species capable of synthesizing a wide group of valuable metabolites, in particular lipases and other hydrolytic enzymes, microbial oil, citric acid, erythritol and γ-decalactone. Processes based on the yeast have GRAS status ("generally recognized as safe") given by Food and Drug Administration. The majority of research communications regarding to Y. lipolytica claim that the yeast species is non-pathogenic. In spite of that, Y. lipolytica, like other fungal species, can cause infections in immunocompromised and critically ill patients. The yeast possess features that facilitate invasion of the host cell (particularly production of hydrolytic enzymes), as well as the protection of the own cells, such as biofilm formation. The aim of this study was to present well-known yeast species Y. lipolytica as a rare opportunistic fungal pathogen. Possible pathogenicity and epidemiology of this yeast species were discussed. Antifungal drugs susceptibility and increasing resistance to azoles in Y. lipolytica yeasts were also presented.

Synthetic biology tools for engineering Yarrowia lipolytica

The non-conventional oleaginous yeast Yarrowia lipolytica shows great industrial promise. It naturally produces certain compounds of interest but can also artificially generate non-native metabolites, thanks to an engineering process made possible by the significant expansion of a dedicated genetic toolbox. In this review, we present recently developed synthetic biology tools that facilitate the manipulation of Y. lipolytica, including 1) DNA assembly techniques, 2) DNA parts for constructing expression cassettes, 3) genome-editing techniques, and 4) computational tools.

Oleaginicity of the yeast strain Saccharomyces cerevisiae D5A

BACKGROUND

The model yeast, Saccharomyces cerevisiae, is not known to be oleaginous. However, an industrial wild-type strain, D5A, was shown to accumulate over 20% storage lipids from glucose when growth is nitrogen-limited compared to no more than 7% lipid accumulation without nitrogen stress.

METHODS AND RESULTS

To elucidate the mechanisms of S. cerevisiae D5A oleaginicity, we compared physiological and metabolic changes; as well as the transcriptional profiles of the oleaginous industrial strain, D5A, and a non-oleaginous laboratory strain, BY4741, under normal and nitrogen-limited conditions using analytic techniques and next-generation sequencing-based RNA-Seq transcriptomics. Transcriptional levels for genes associated with fatty acid biosynthesis, nitrogen metabolism, amino acid catabolism, as well as the pentose phosphate pathway and ethanol oxidation in central carbon (C) metabolism, were up-regulated in D5A during nitrogen deprivation. Despite increased carbon flux to lipids, most gene-encoding enzymes involved in triacylglycerol (TAG) assembly were expressed at similar levels regardless of the varying nitrogen concentrations in the growth media and strain backgrounds. Phospholipid turnover also contributed to TAG accumulation through increased precursor production with the down-regulation of subsequent phospholipid synthesis steps. Our results also demonstrated that nitrogen assimilation via the glutamate-glutamine pathway and amino acid metabolism, as well as the fluxes of carbon and reductants from central C metabolism, are integral to the general oleaginicity of D5A, which resulted in the enhanced lipid storage during nitrogen deprivation.

CONCLUSION

This work demonstrated the disequilibrium and rebalance of carbon and nitrogen contribution to the accumulation of lipids in the oleaginous yeast S. cerevisiae D5A. Rather than TAG assembly from acyl groups, the major switches for the enhanced lipid accumulation of D5A (i.e., fatty acid biosynthesis) are the increases of cytosolic pools of acetyl-CoA and NADPH, as well as alternative nitrogen assimilation.

Mitotic defects in fission yeast lipid metabolism 'cut' mutants are suppressed by ammonium chloride

Fission yeast 'cut' mutants show defects in temporal coordination of nuclear division with cytokinesis, resulting in aberrant mitosis and lethality. Among other causes, the 'cut' phenotype can be triggered by genetic or chemical perturbation of lipid metabolism, supposedly resulting in shortage of membrane phospholipids and insufficient nuclear envelope expansion during anaphase. Interestingly, penetrance of the 'cut' phenotype in mutants of the transcription factor cbf11 and acetyl-coenzyme A carboxylase cut6, both related to lipid metabolism, is highly dependent on growth media, although the specific nutrient(s) affecting 'cut' occurrence is not known. In this study, we set out to identify the growth media component(s) responsible for 'cut' phenotype suppression in Δcbf11 and cut6-621 cells. We show that mitotic defects occur rapidly in Δcbf11 cells upon shift from the minimal EMM medium ('cut' suppressing) to the complex YES medium ('cut' promoting). By growing cells in YES medium supplemented with individual EMM components, we identified ammonium chloride, an efficiently utilized nitrogen source, as a specific and potent suppressor of the 'cut' phenotype in both Δcbf11 and cut6-621. Furthermore, we found that ammonium chloride boosts lipid droplet formation in wild-type cells. Our findings suggest a possible involvement of nutrient-responsive signaling in 'cut' suppression.

Morphological response of Yarrowia lipolytica under stress of heavy metals

The marine dimorphic yeast Yarrowia lipolytica has been proposed as a suitable model for the dimorphism study. In this study, the morphological behaviour of two marine strains of Y. lipolytica (NCIM 3589 and NCIM 3590) was studied under stress of different heavy metals. Scanning electron microscopy was used to investigate the morphological features of yeast cells. This study revealed that the normal ellipsoidal shape of yeast cells was changed into oval, rounded, or elongated in response to different heavy-metal stress. Light microscopy was also used to investigate individual properties of yeast cells. The average cell length and radius of both marine strains was increased with increasing concentrations of heavy-metal ions. In addition, the elongation factor was calculated and was increased in the presence of heavy metals like Pb(II), Co(II), Cr(III), Cr(VI), and Zn(II) under the static conditions.

Metabolic engineering in the host Yarrowia lipolytica

The nonconventional, oleaginous yeast, Yarrowia lipolytica is rapidly emerging as a valuable host for the production of a variety of both lipid and nonlipid chemical products. While the unique genetics of this organism pose some challenges, many new metabolic engineering tools have emerged to facilitate improved genetic manipulation in this host. This review establishes a case for Y. lipolytica as a premier metabolic engineering host based on innate metabolic capacity, emerging synthetic tools, and engineering examples. The metabolism underlying the lipid accumulation phenotype of this yeast as well as high flux through acyl-CoA precursors and the TCA cycle provide a favorable metabolic environment for expression of relevant heterologous pathways. These properties allow Y. lipolytica to be successfully engineered for the production of both native and nonnative lipid, organic acid, sugar and acetyl-CoA derived products. Finally, this host has unique metabolic pathways enabling growth on a wide range of carbon sources, including waste products. The expansion of carbon sources, together with the improvement of tools as highlighted here, have allowed this nonconventional organism to act as a cellular factory for valuable chemicals and fuels.

Holistic Approaches in Lipid Production by Yarrowia lipolytica

Concerns about climate change have driven research on the production of lipid-derived biofuels as an alternative and renewable liquid fuel source. Using oleaginous yeasts for lipid synthesis creates the potential for cost-effective industrial-scale operations due to their ability to reach high lipid titer, yield, and productivity resulting from their unique metabolism. Yarrowia lipolytica is the model oleaginous yeast, with the best-studied lipid metabolism, the greatest number of genetic tools, and a fully sequenced genome. In this review we highlight multiomics studies that elucidate the mechanisms allowing this yeast to achieve lipid overaccumulation and then present several major metabolic engineering efforts that enhanced the production metrics in Y. lipolytica. Recent achievements that applied novel engineering strategies are emphasized.

Enhanced squalene biosynthesis in Yarrowia lipolytica based on metabolically engineered acetyl-CoA metabolism

As a bioactive triterpenoid, squalene is widely used in the food industry, cosmetics, and pharmacology. Squalene's major commercial sources are the liver oil of deep-sea sharks and plant oils. In this study, we focused on the enhancement of squalene biosynthesis in Yarrowia lipolytica, with particular attention to the engineering of acetyl-CoA metabolism based on genome-scale metabolic reaction network analysis. Although the overexpression of the rate-limiting endogenous ylHMG1 (3-hydroxy-3-methylglutaryl-CoA reductase gene) could improve squalene synthesis by 3.2-fold over that by the control strain, the availability of the key intracellular precursor, acetyl-CoA, was found to play a more significant role in elevating squalene production. Analysis of metabolic networks with the newly constructed genome-scale metabolic model of Y. lipolytica iYL_2.0 showed that the acetyl-CoA pool size could be increased by redirecting carbon flux of pyruvate dehydrogenation towards the ligation of acetate and CoA or the cleavage of citrate to form oxaloacetate and acetyl-CoA. The overexpression of either acetyl-CoA synthetase gene from Salmonella enterica (acs*) or the endogenous ATP citrate lyase gene (ylACL1) resulted in a more than 50% increase in the cytosolic acetyl-CoA level. Moreover, iterative chromosomal integration of the ylHMG1, asc*, and ylACL1 genes resulted in a significant improvement in squalene production (16.4-fold increase in squalene content over that in the control strain). We also found that supplementation with 10 mM citrate in a flask culture further enhanced squalene production to 10 mg/g DCW. The information obtained in this study demonstrates that rationally engineering acetyl-CoA metabolism to ensure the supply of this key metabolic precursor is an efficient strategy for the enhancement of squalene biosynthesis.

Prediction of Metabolite Concentrations, Rate Constants and Post-Translational Regulation Using Maximum Entropy-Based Simulations with Application to Central Metabolism of Neurospora crassa

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin’s 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Genome-scale model-driven strain design for dicarboxylic acid production in Yarrowia lipolytica

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.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0542-5) contains supplementary material, which is available to authorized users.

Functional genomics of lipid metabolism in the oleaginous yeast Rhodosporidium toruloides

The basidiomycete yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) accumulates high concentrations of lipids and carotenoids from diverse carbon sources. It has great potential as a model for the cellular biology of lipid droplets and for sustainable chemical production. We developed a method for high-throughput genetics (RB-TDNAseq), using sequence-barcoded Agrobacterium tumefaciens T-DNA insertions. We identified 1,337 putative essential genes with low T-DNA insertion rates. We functionally profiled genes required for fatty acid catabolism and lipid accumulation, validating results with 35 targeted deletion strains. We identified a high-confidence set of 150 genes affecting lipid accumulation, including genes with predicted function in signaling cascades, gene expression, protein modification and vesicular trafficking, autophagy, amino acid synthesis and tRNA modification, and genes of unknown function. These results greatly advance our understanding of lipid metabolism in this oleaginous species and demonstrate a general approach for barcoded mutagenesis that should enable functional genomics in diverse fungi.

The fungus Rhodosporidium toruloides can grow on substances extracted from plant matter that is inedible to humans such as corn stalks, wood pulp, and grasses. Under some growth conditions, the fungus can accumulate massive stores of hydrocarbon-rich fats and pigments. A community of scientists and engineers has begun genetically modifying R. toruloides to convert these naturally produced fats and pigments into fuels, chemicals and medicines. These could form sustainable replacements for products made from petroleum or harvested from threatened animal and plant species.

Fungi, plants, animals and other eukaryotes store fat in specialized compartments called lipid droplets. The genes that control the metabolism – the production, use and storage – of fat in lipid bodies have been studied in certain eukaryotes, including species of yeast. However, R. toruloides is only distantly related to the most well-studied of these species. This means that we cannot be certain that a gene will play the same role in R. toruloides as in those species.

To assemble the most comprehensive list possible of the genes in R. toruloides that affect the production, use, or storage of fat in lipid bodies, Coradetti, Pinel et al. constructed a population of hundreds of thousands of mutant fungal strains, each with its own unique DNA ‘barcode’. The effects that mutations in over 6,000 genes had on growth and fat accumulation in these fungi were measured simultaneously in several experiments. This general approach is not new, but technical limitations had, until now, restricted its use in fungi to a few species.

Coradetti, Pinel et al. identified hundreds of genes that affected the ability of R. toruloides to metabolise fat. Many of these genes were related to genes with known roles in fat metabolism in other eukaryotes. Other genes are involved in different cell processes, such as the recycling of waste products in the cell. Their identification adds weight to the view that the links between these cellular processes and fat metabolism are deep and widespread amongst eukaryotes. Finally, some of the genes identified by Coradetti, Pinel et al. are not closely related to any well-studied genes. Further study of these genes could help us to understand why R. toruloides can accumulate much larger amounts of fat than most other fungi.

The methods developed by Coradetti, Pinel et al. should be possible to implement in many species of fungi. As a result these techniques may eventually contribute to the development of new treatments for human fungal diseases, the protection of important food crops, and a deeper understanding of the roles various fungi play in the broader ecosystem.

Metabolomic elucidation of the effects of media and carbon sources on fatty acid production by Yarrowia lipolytica

Lipid production by oleaginous Yarrowia lipolytica depends highly on culture environments, such as carbon sources, carbon/nitrogen (C/N) ratios, types of media, and cellular growth phases. In this study, the effects of media and carbon sources on lipid and metabolite production were investigated by profiling fatty acids and intracellular metabolites of Y. lipolytica grown in various media. The highest total fatty acid yield 114.04 ± 6.23 mg/g dry cell weight was achieved by Y. lipolytica grown in minimal medium with glycerol (SCG) in the exponential phase. The high lipid production by Y. lipolytica in SCG was presumed to be due to the higher C/N ratio in SCG than in the complex media. Moreover, glycerol promoted lipid production better than glucose in both complex and minimal media because glycerol can easily incorporate into the core of triglycerides. Metabolite profiling revealed that levels of long-chain fatty acids, such as stearic acid, palmitic acid, and arachidic acid, increased in SCG medium. Meanwhile, in complex media supplemented with either glucose or glycerol, levels of amino acids, such as cysteine, methionine, and glycine, highly increased. This metabolomic approach could be applied to modulate the global metabolic network of Y. lipolytica for producing lipids and other valuable products.

Metabolic Engineering of Oleaginous Yeasts for Production of Fuels and Chemicals

Oleaginous yeasts have been increasingly explored for production of chemicals and fuels via metabolic engineering. Particularly, there is a growing interest in using oleaginous yeasts for the synthesis of lipid-related products due to their high lipogenesis capability, robustness, and ability to utilize a variety of substrates. Most of the metabolic engineering studies in oleaginous yeasts focused on Yarrowia that already has plenty of genetic engineering tools. However, recent advances in systems biology and synthetic biology have provided new strategies and tools to engineer those oleaginous yeasts that have naturally high lipid accumulation but lack genetic tools, such as Rhodosporidium, Trichosporon, and Lipomyces. This review highlights recent accomplishments in metabolic engineering of oleaginous yeasts and recent advances in the development of genetic engineering tools in oleaginous yeasts within the last 3 years.

A synthetic biology approach to transform Yarrowia lipolytica into a competitive biotechnological producer of β-carotene

The increasing market demands of β-carotene as colorant, antioxidant and vitamin precursor, requires novel biotechnological production platforms. Yarrowia lipolytica, is an industrial organism unable to naturally synthesize carotenoids but with the ability to produce high amounts of the precursor Acetyl-CoA. We first found that a lipid overproducer strain was capable of producing more β-carotene than a wild type after expressing the heterologous pathway. Thereafter, we developed a combinatorial synthetic biology approach base on Golden Gate DNA assembly to screen the optimum promoter-gene pairs for each transcriptional unit expressed. The best strain reached a production titer of 1.5 g/L and a maximum yield of 0.048 g/g of glucose in flask. β-carotene production was further increased in controlled conditions using a fed-batch fermentation. A total production of β-carotene of 6.5 g/L and 90 mg/g DCW with a concomitant production of 42.6 g/L of lipids was achieved. Such high titers suggest that engineered Y. lipolytica is a competitive producer organism of β-carotene.

Elucidation of complexity and prediction of interactions in microbial communities

Microorganisms engage in complex interactions with other members of the microbial community, higher organisms as well as their environment. However, determining the exact nature of these interactions can be challenging due to the large number of members in these communities and the manifold of interactions they can engage in. Various omic data, such as 16S rRNA gene sequencing, shotgun metagenomics, metatranscriptomics, metaproteomics and metabolomics, have been deployed to unravel the community structure, interactions and resulting community dynamics in situ. Interpretation of these multi-omic data often requires advanced computational methods. Modelling approaches are powerful tools to integrate, contextualize and interpret experimental data, thus shedding light on the underlying processes shaping the microbiome. Here, we review current methods and approaches, both experimental and computational, to elucidate interactions in microbial communities and to predict their responses to perturbations.