Genome-scale metabolic model of oleaginous yeast Papiliotrema laurentii

Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used industrially for production of biofuels and bioproducts. Developing strains with better performance in biomanufacturing contexts requires a systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated for several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. Here, we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test-learn cycles, ultimately advancing the field of fungal biomanufacturing.

Effect of carbon and nitrogen concentrations on lipid accumulation and regulation of acetyl-CoA carboxylase in Papiliotrema laurentii

Biodiesel is an interesting alternative to petroleum diesel as it is renewable, biodegradable, and has a low pollutant content. Yeast oils can be used for biodiesel production instead of edible oils, mitigating the use of arable land and water for biodiesel production. Maximum lipid accumulation is reached at 48 h of cultivation by the oleaginous yeast Papiliotrema laurentii UFV-1. Nevertheless, the effects of carbon and nitrogen concentrations on lipid accumulation, as well as the regulation of lipid metabolism in this yeast are still not well-characterised. Therefore, this work evaluated the effects of carbon and nitrogen concentrations on the lipid accumulation in P. laurentti, the expression of the ACC gene, and the activity of the enzyme acetyl-CoA carboxylase (ACCase) in different carbon:nitrogen ratios (C:N) and glucose concentrations. The variation of ammonium sulfate concentration did not affect the growth and lipid accumulation in P. laurentii UFV-1. On the other hand, glucose concentration remarkably influenced biomass and lipid production by this yeast. Therefore, the carbon concentration is more important than the nitrogen concentration for lipid production by P. laurentii UFV-1. Importantly, the levels of both ACC gene expression and ACCase activity were maximum during the late-exponential growth phase and decreased after reaching the highest lipid contents, which was easier evidenced during the accumulation and maximum lipid levels. As such, the reduction of ACCase enzyme activity seems to be related to the decrease in the expression level of the ACC gene.

New Papiliotrema laurentii UFV-1 strains with improved acetic acid tolerance selected by adaptive laboratory evolution

The production of yeast oil from lignocellulosic biomasses is impaired by inhibitors formed during the pretreatment step, mainly acetic acid. Herein, we applied Adaptive Laboratory Evolution (ALE) to select three Acetic acid Tolerant Strains (ATS) of P. laurentii UFV-1. Different phenotypes emerged alongside evolution. The ATS II presented trade-offs in the absence of acetic acid, suggesting that it displays a specialized phenotype of tolerance to growth on organic acids. On the other hand, ATS I and ATS III presented phenotypes associated with the behavior of generalists. ATS I was considered the most promising evolved strain as it displayed the oleaginous phenotype in all conditions tested. Thus, we applied whole-genome sequencing to detect the mutations that emerged in this strain during the ALE. We found alterations in genes encoding proteins involved in different cellular functions, including multidrug resistance (MDR) transporters, energy metabolism, detoxification, coenzyme recycling, and cell envelope remodeling. To evaluate acetic acid stress responses, both parental and ATS I strains were cultivated in chemostat mode in the absence and presence of acetic acid. In contrast to ATS I, the parental strain presented alterations in the cell envelope and cell size under acetic acid stress conditions. Furthermore, the parental strain and the ATS I presented differences regarding acetic acid assimilation. Contrary to the parental strain, the ATS I displayed an increase in unsaturated fatty acid content irrespective of acetic acid stress, which might be related to improved tolerance to acetic acid. Altogether, these results provided insights into the mechanisms involved with the acetic acid tolerance displayed by ATS I and the responses of P. laurentii to this stressful condition.

Recent Development of Advanced Biotechnology in the Oleaginous Fungi for Arachidonic Acid Production

Arachidonic acid is an essential ω-6 polyunsaturated fatty acid, which plays a significant role in cardiovascular health and neurological development, leading to its wide use in the food and pharmaceutical industries. Traditionally, ARA is obtained from deep-sea fish oil. However, this source is limited by season and is depleting the already threatened global fish stocks. With the rapid development of synthetic biology in recent years, oleaginous fungi have gradually attracted increasing attention as promising microbial sources for large-scale ARA production. Numerous advanced technologies including metabolic engineering, dynamic regulation of fermentation conditions, and multiomics analysis were successfully adapted to increase ARA synthesis. This review summarizes recent advances in the bioengineering of oleaginous fungi for ARA production. Finally, perspectives for future engineering approaches are proposed to further improve the titer yield and productivity of ARA.