Search

Sunscreen has been used for thousands of years to protect skin from ultraviolet radiation. However, the use of modern commercial sunscreen containing oxybenzone, ZnO, and TiO2 has raised concerns due to their negative effects on human health and the environment. In this study, we aim to establish an efficient microbial platform for production of shinorine, a UV light absorbing compound with anti-aging properties. First, we methodically selected an appropriate host for shinorine production by analyzing central carbon flux distribution data from prior studies alongside predictions from genome-scale metabolic models (GEMs). We enhanced shinorine productivity through CRISPRi-mediated downregulation and utilized shotgun proteomics to pinpoint potential competing pathways. Simultaneously, we improved the shinorine biosynthetic pathway by refining its design, optimizing promoter usage, and altering the strength of ribosome binding sites. Finally, we conducted amino acid feeding experiments under various conditions to identify the key limiting factors in shinorine production. The study combines meta-analysis of 13C-metabolic flux analysis, GEMs, synthetic biology, CRISPRi-mediated gene downregulation, and omics analysis to improve shinorine production, demonstrating the potential of Pseudomonas putida KT2440 as platform for shinorine production.

Sustainable aviation fuel (SAF) will significantly impact global warming in the aviation sector, and important SAF targets are emerging. Isoprenol is a precursor for a promising SAF compound DMCO (1,4-dimethylcyclooctane) and has been produced in several engineered microorganisms. Recently, Pseudomonas putida has gained interest as a future host for isoprenol bioproduction as it can utilize carbon sources from inexpensive plant biomass. Here, we engineer metabolically versatile host P. putida for isoprenol production. We employ two computational modeling approaches (Bilevel optimization and Constrained Minimal Cut Sets) to predict gene knockout targets and optimize the "IPP-bypass" pathway in P. putida to maximize isoprenol production. Altogether, the highest isoprenol production titer from P. putida was achieved at 3.5 g/L under fed-batch conditions. This combination of computational modeling and strain engineering on P. putida for an advanced biofuels production has vital significance in enabling a bioproduction process that can use renewable carbon streams.

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.

Multidimensional measurements using state-of-the-art separations and mass spectrometry provide advantages in untargeted metabolomics analyses for studying biological and environmental bio-chemical processes. However, the lack of rapid analytical methods and robust algorithms for these heterogeneous data has limited its application. Here, we develop and evaluate a sensitive and high-throughput analytical and computational workflow to enable accurate metabolite profiling. Our workflow combines liquid chromatography, ion mobility spectrometry and data-independent acquisition mass spectrometry with PeakDecoder, a machine learning-based algorithm that learns to distinguish true co-elution and co-mobility from raw data and calculates metabolite identification error rates. We apply PeakDecoder for metabolite profiling of various engineered strains of Aspergillus pseudoterreus, Aspergillus niger, Pseudomonas putida and Rhodosporidium toruloides. Results, validated manually and against selected reaction monitoring and gas-chromatography platforms, show that 2683 features could be confidently annotated and quantified across 116 microbial sample runs using a library built from 64 standards.

The global regulator LaeA controls secondary metabolism in diverse Aspergillus species. Here we explored its role in regulation of itaconic acid production in Aspergillus pseudoterreus. To understand its role in regulating metabolism, we deleted and overexpressed laeA, and assessed the transcriptome, proteome, and secreted metabolome prior to and during initiation of phosphate limitation induced itaconic acid production. We found that secondary metabolite clusters, including the itaconic acid biosynthetic gene cluster, are regulated by laeA and that laeA is required for high yield production of itaconic acid. Overexpression of LaeA improves itaconic acid yield at the expense of biomass by increasing the expression of key biosynthetic pathway enzymes and attenuating the expression of genes involved in phosphate acquisition and scavenging. Increased yield was observed in optimized conditions as well as conditions containing excess nutrients that may be present in inexpensive sugar containing feedstocks such as excess phosphate or complex nutrient sources. This suggests that global regulators of metabolism may be useful targets for engineering metabolic flux that is robust to environmental heterogeneity.

BackgroundMitigation of climate change requires that new routes for the production of fuels and chemicals be as oil-independent as possible. The microbial conversion of lignocellulosic feedstocks into terpene-based biofuels and bioproducts represents one such route. This work builds upon previous demonstrations that the single-celled carotenogenic basidiomycete, Rhodosporidium toruloides, is a promising host for the production of terpenes from lignocellulosic hydrolysates.ResultsThis study focuses on the optimization of production of the monoterpene 1,8-cineole and the sesquiterpene α-bisabolene in R. toruloides. The α-bisabolene titer attained in R. toruloides was found to be proportional to the copy number of the bisabolene synthase (BIS) expression cassette, which in turn influenced the expression level of several native mevalonate pathway genes. The addition of more copies of BIS under a stronger promoter resulted in production of α-bisabolene at 2.2 g/L from lignocellulosic hydrolysate in a 2-L fermenter. Production of 1,8-cineole was found to be limited by availability of the precursor geranylgeranyl pyrophosphate (GPP) and expression of an appropriate GPP synthase increased the monoterpene titer fourfold to 143 mg/L at bench scale. Targeted mevalonate pathway metabolite analysis suggested that 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGR), mevalonate kinase (MK) and phosphomevalonate kinase (PMK) may be pathway bottlenecks are were therefore selected as targets for overexpression. Expression of HMGR, MK, and PMK orthologs and growth in an optimized lignocellulosic hydrolysate medium increased the 1,8-cineole titer an additional tenfold to 1.4 g/L. Expression of the same mevalonate pathway genes did not have as large an impact on α-bisabolene production, although the final titer was higher at 2.6 g/L. Furthermore, mevalonate pathway intermediates accumulated in the mevalonate-engineered strains, suggesting room for further improvement.ConclusionsThis work brings R. toruloides closer to being able to make industrially relevant quantities of terpene from lignocellulosic biomass.

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides. High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

Biological engineering of microorganisms to produce value-added chemicals is a promising route to sustainable manufacturing. However, overproduction of metabolic intermediates at high titer, rate, and yield from inexpensive substrates is challenging in non-model systems where limited information is available regarding metabolic flux and its control in production conditions. Integrated multi-omic analyses of engineered strains offers an in-depth look at metabolites and proteins directly involved in growth and production of target and non-target bioproducts. Here we applied multi-omic analyses to overproduction of the polymer precursor 3-hydroxypropionic acid (3HP) in the filamentous fungus Aspergillus pseudoterreus. A synthetic pathway consisting of aspartate decarboxylase, beta-alanine pyruvate transaminase, and 3HP dehydrogenase was designed and built for A. pseudoterreus. Strains with single- and multi-copy integration events were isolated and multi-omics analysis consisting of intracellular and extracellular metabolomics and targeted and global proteomics was used to interrogate the strains in shake-flask and bioreactor conditions. Production of a variety of co-products (organic acids and glycerol) and oxidative degradation of 3HP were identified as metabolic pathways competing with 3HP production. Intracellular accumulation of nitrogen as 2,4-diaminobutanoate was identified as an off-target nitrogen sink that may also limit flux through the engineered 3HP pathway. Elimination of the high-expression oxidative 3HP degradation pathway by deletion of a putative malonate semialdehyde dehydrogenase improved the yield of 3HP by 3.4 × after 10 days in shake-flask culture. This is the first report of 3HP production in a filamentous fungus amenable to industrial scale biomanufacturing of organic acids at high titer and low pH.

Biology has changed radically in the past two decades, growing from a purely descriptive science into also a design science. The availability of tools that enable the precise modification of cells, as well as the ability to collect large amounts of multimodal data, open the possibility of sophisticated bioengineering to produce fuels, specialty and commodity chemicals, materials, and other renewable bioproducts. However, despite new tools and exponentially increasing data volumes, synthetic biology cannot yet fulfill its true potential due to our inability to predict the behavior of biological systems. Here, we showcase a set of computational tools that, combined, provide the ability to store, visualize, and leverage multiomics data to predict the outcome of bioengineering efforts. We show how to upload, visualize, and output multiomics data, as well as strain information, into online repositories for several isoprenol-producing strain designs. We then use these data to train machine learning algorithms that recommend new strain designs that are correctly predicted to improve isoprenol production by 23%. This demonstration is done by using synthetic data, as provided by a novel library, that can produce credible multiomics data for testing algorithms and computational tools. In short, this paper provides a step-by-step tutorial to leverage these computational tools to improve production in bioengineered strains.

Machine learning provides researchers a unique opportunity to make metabolic engineering more predictable. In this review, we offer an introduction to this discipline in terms that are relatable to metabolic engineers, as well as providing in-depth illustrative examples leveraging omics data and improving production. We also include practical advice for the practitioner in terms of data management, algorithm libraries, computational resources, and important non-technical issues. A variety of applications ranging from pathway construction and optimization, to genetic editing optimization, cell factory testing, and production scale-up are discussed. Moreover, the promising relationship between machine learning and mechanistic models is thoroughly reviewed. Finally, the future perspectives and most promising directions for this combination of disciplines are examined.

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides. High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

BackgroundRhodosporidium toruloides is a promising host for the production of bioproducts from lignocellulosic biomass. A key prerequisite for efficient pathway engineering is the availability of robust genetic tools and resources. However, there is a lack of characterized promoters to drive expression of heterologous genes for strain engineering in R. toruloides.ResultsThis data describes a set of native R. toruloides promoters, characterized over time in four different media commonly used for cultivation of this yeast. The promoter sequences were selected using transcriptional analysis and several of them were found to drive expression bidirectionally. Promoter expression strength was determined by measurement of EGFP and mRuby2 reporters by flow cytometry. A total of 20 constitutive promoters (12 monodirectional and 8 bidirectional) were found, and are expected to be of potential value for genetic engineering of R. toruloides.ConclusionsA set of robust and constitutive promoters to facilitate genetic engineering of R. toruloides is presented here, ranging from a promoter previously used for this purpose (P7, glyceraldehyde 3-phosphate dehydrogenase, GAPDH) to stronger monodirectional (e.g., P15, mitochondrial adenine nucleotide translocator, ANT) and bidirectional (e.g., P9 and P9R, histones H3 and H4, respectively) promoters. We also identified promoters that may be useful for specific applications such as late-stage expression (e.g., P3, voltage-dependent anion channel protein 2, VDAC2). This set of characterized promoters significantly expands the range of engineering tools available for this yeast and can be applied in future metabolic engineering studies.

Isopentenyl pyrophosphate (IPP) toxicity presents a challenge in engineered microbial systems since its formation is unavoidable in terpene biosynthesis. In this work, we develop an experimental platform to study IPP toxicity in isoprenol-producing Escherichia coli. We first characterize the physiological response to IPP accumulation, demonstrating that elevated IPP levels are linked to growth inhibition, reduced cell viability, and plasmid instability. We show that IPP toxicity selects for pathway "breakage", using proteomics to identify a reduction in phosphomevalonate kinase (PMK) as a probable recovery mechanism. Next, using multi-omics data, we demonstrate that endogenous E. coli metabolism is globally impacted by IPP accumulation, which slows nutrient uptake, decreases ATP levels, and perturbs nucleotide metabolism. We also observe the extracellular accumulation of IPP and present preliminary evidence that IPP can be transported by E. coli, findings that might be broadly relevant for the study of isoprenoid biosynthesis. Finally, we discover that IPP accumulation leads to the formation of ApppI, a nucleotide analog of IPP that may contribute to observed toxicity phenotypes. This comprehensive assessment of IPP stress suggests potential strategies for the alleviation of prenyl diphosphate toxicity and highlights possible engineering targets for improved IPP flux and high titer isoprenoid production.

Plant biomass, once reduced to its composite sugars, can be converted to fuel substitutes. One means of overcoming the recalcitrance of lignocellulose is pretreatment followed by enzymatic hydrolysis. However, currently available commercial enzyme cocktails are inhibited in the presence of residual pretreatment chemicals. Recent studies have identified a number of cellulolytic enzymes from bacteria that are tolerant to pretreatment chemicals such as ionic liquids. The challenge now is generation of these enzymes in copious amounts, an arena where fungal organisms such as Aspergillus niger have proven efficient. Fungal host strains still need to be engineered to increase production titers of heterologous protein over native enzymes, which has been a difficult task. Here, we developed a forward genetics screen coupled with whole-genome resequencing to identify specific lesions responsible for a protein hyper-production phenotype in A. niger. This strategy successfully identified novel targets, including a low-affinity glucose transporter, MstC, whose deletion significantly improved secretion of recombinant proteins driven by a glucoamylase promoter.

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 ., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Czajka, Han, Kim, Mondo, Hofstad, Robles, Haridas, Riley, LaButti, Pangilinan, Andreopoulos, Lipzen, Yan, Wang, Ng, Grigoriev, Spatafora, Magnuson, Baker and Pomraning.)