Rapid metabolic analysis of Rhodococcus opacus PD630 via parallel 13 C‐metabolite fingerprinting

Unleashing the capacity of Rhodococcus for converting lignin into lipids

Bioconversion of bioresources/wastes (e.g., lignin, chemical pulping byproducts) represents a promising approach for developing a bioeconomy to help address growing energy and materials demands. Rhodococcus, a promising microbial strain, utilizes numerous carbon sources to produce lipids, which are precursors for synthesizing biodiesel and aviation fuels. However, compared to chemical conversion, bioconversion involves living cells, which is a more complex system that needs further understanding and upgrading. Various wastes amenable to bioconversion are reviewed herein to highlight the potential of Rhodococci for producing lipid-derived bioproducts. In light of the abundant availability of these substrates, Rhodococcus' metabolic pathways converting them to lipids are analyzed from a "beginning-to-end" view. Based on an in-depth understanding of microbial metabolic routes, genetic modifications of Rhodococcus by employing emerging tools (e.g., multiplex genome editing, biosensors, and genome-scale metabolic models) are presented for promoting the bioconversion. Co-solvent enhanced lignocellulose fractionation (CELF) strategy facilitates the generation of a lignin-derived aromatic stream suitable for the Rhodococcus' utilization. Novel alkali sterilization (AS) and elimination of thermal sterilization (ETS) approaches can significantly enhance the bioaccessibility of lignin and its derived aromatics in aqueous fermentation media, which promotes lipid titer significantly. In order to achieve value-added utilization of lignin, biodiesel and aviation fuel synthesis from lignin and lipids are further discussed. The possible directions for unleashing the capacity of Rhodococcus through synergistically modifying microbial strains, substrates, and fermentation processes are proposed toward a sustainable biological lignin valorization.

A High-Quality Genome-Scale Model for Rhodococcus opacus Metabolism

Rhodococcus opacus is a bacterium that has a high tolerance to aromatic compounds and can produce significant amounts of triacylglycerol (TAG). Here, we present iGR1773, the first genome-scale model (GSM) of R. opacus PD630 metabolism based on its genomic sequence and associated data. The model includes 1773 genes, 3025 reactions, and 1956 metabolites, was developed in a reproducible manner using CarveMe, and was evaluated through Metabolic Model tests (MEMOTE). We combine the model with two Constraint-Based Reconstruction and Analysis (COBRA) methods that use transcriptomics data to predict growth rates and fluxes: E-Flux2 and SPOT (Simplified Pearson Correlation with Transcriptomic data). Growth rates are best predicted by E-Flux2. Flux profiles are more accurately predicted by E-Flux2 than flux balance analysis (FBA) and parsimonious FBA (pFBA), when compared to 44 central carbon fluxes measured by 13C-Metabolic Flux Analysis (13C-MFA). Under glucose-fed conditions, E-Flux2 presents an R2 value of 0.54, while predictions based on pFBA had an inferior R2 of 0.28. We attribute this improved performance to the extra activity information provided by the transcriptomics data. For phenol-fed metabolism, in which the substrate first enters the TCA cycle, E-Flux2's flux predictions display a high R2 of 0.96 while pFBA showed an R2 of 0.93. We also show that glucose metabolism and phenol metabolism function with similar relative ATP maintenance costs. These findings demonstrate that iGR1773 can help the metabolic engineering community predict aromatic substrate utilization patterns and perform computational strain design.

In situ diversity of metabolism and carbon use efficiency among soil bacteria

The central carbon (C) metabolic network harvests energy to power the cell and feed biosynthesis for growth. In pure cultures, bacteria use some but not all of the network's major pathways, such as glycolysis and pentose phosphate and Entner-Doudoroff pathways. However, how these pathways are used in microorganisms in intact soil communities is unknown. Here, we analyzed the incorporation of ¹³C from glucose isotopomers into phospholipid fatty acids. We showed that groups of Gram-positive and Gram-negative bacteria in an intact agricultural soil used different pathways to metabolize glucose. They also differed in C use efficiency (CUE), the efficiency with which a substrate is used for biosynthesis. Our results provide experimental evidence for diversity among microbes in the organization of their central carbon metabolic network and CUE under in situ conditions. These results have important implications for our understanding of how community composition affects soil C cycling and organic matter formation.

Current challenges of hydrothermal treated wastewater (HTWW) for environmental applications and their perspectives: A review

Hydrothermal treatment (HT) is an emerged thermochemical approach for the utilization of biomass. In the last decade, intense research has been conducted on bio-oil and hydrochar, during which extensive amount of hydrothermal treated wastewater (HTWW) is produced, containing large amount of organic compounds along with several toxic chemicals. The composition of HTWW is highly dependent on the process conditions and organic composition of biomass, which determines its further utilization. The current study provides a comprehensive overview of recent advancements in HTWW utilization and its properties which can be changed by varying different parameters like temperature, residence time, solid concentration, mass ratio and catalyst including types of biomasses. HTWW characterization, parameters, reaction mechanism and its application were also summarized. By considering the challenges of HTWW, some suggestions and proposed methodology to overcome the bottleneck are provided.

Rhodococcus as Biofactories for Microbial Oil Production

Bacteria belonging to the Rhodococcus genus are frequent components of microbial communities in diverse natural environments. Some rhodococcal species exhibit the outstanding ability to produce significant amounts of triacylglycerols (TAG) (>20% of cellular dry weight) in the presence of an excess of the carbon source and limitation of the nitrogen source. For this reason, they can be considered as oleaginous microorganisms. As occurs as well in eukaryotic single-cell oil (SCO) producers, these bacteria possess specific physiological properties and molecular mechanisms that differentiate them from other microorganisms unable to synthesize TAG. In this review, we summarized several of the well-characterized molecular mechanisms that enable oleaginous rhodococci to produce significant amounts of SCO. Furthermore, we highlighted the ability of these microorganisms to degrade a wide range of carbon sources coupled to lipogenesis. The qualitative and quantitative oil production by rhodococci from diverse industrial wastes has also been included. Finally, we summarized the genetic and metabolic approaches applied to oleaginous rhodococci to improve SCO production. This review provides a comprehensive and integrating vision on the potential of oleaginous rhodococci to be considered as microbial biofactories for microbial oil production.

An Improved CRISPR Interference Tool to Engineer Rhodococcus opacus

Rhodococcus opacus is a nonmodel bacterium that is well suited for valorizing lignin. Despite recent advances in our systems-level understanding of its versatile metabolism, studies of its gene functions at a single gene level are still lagging. Elucidating gene functions in nonmodel organisms is challenging due to limited genetic engineering tools that are convenient to use. To address this issue, we developed a simple gene repression system based on CRISPR interference (CRISPRi). This gene repression system uses a T7 RNA polymerase system to express a small guide RNA, demonstrating improved repression compared to the previously demonstrated CRISPRi system (i.e., the maximum repression efficiency improved from 58% to 85%). Additionally, our cloning strategy allows for building multiple CRISPRi plasmids in parallel without any PCR step, facilitating the engineering of this GC-rich organism. Using the improved CRISPRi system, we confirmed the annotated roles of four metabolic pathway genes, which had been identified by our previous transcriptomic analysis to be related to the consumption of benzoate, vanillate, catechol, and acetate. Furthermore, we showed our tool's utility by demonstrating the inducible accumulation of muconate that is a precursor of adipic acid, an important monomer for nylon production. While the maximum muconate yield obtained using our tool was 30% of the yield obtained using gene knockout, our tool showed its inducibility and partial repressibility. Our CRISPRi tool will be useful to facilitate functional studies of this nonmodel organism and engineer this promising microbial chassis for lignin valorization.

Investigating the potential use of an oleaginous bacterium, Rhodococcus opacus PD630, for nano-TiO2 remediation

The occurrence of titanium dioxide nanoparticles (nTiO₂), in the effluents released from wastewater treatment plants, has raised concerns. The fate of nTiO₂ and their potential impact on organisms from different ecosystems are widely investigated. For the first time, in this work, we report the responses of an oleaginous bacteria Rhodococcus opacus PD630, belonging to an ecologically important genus Rhodococcus to environmentally relevant concentrations of nTiO_2, under dark and UV light conditions. We observed a dose-dependent increase in nTiO₂ uptake by the bacteria that reached a maximum of 1.4 mg nTiO₂ (g cell)^-1 under mid-log UV exposure, corresponding to 97% uptake. The nTiO₂ induced oxidative stress in bacteria that increased from 25.1 to a maximum of 100.3, 44.1, and 51.7 μmol ^.OH (g cell)^-1 under dark, continuous, and mid-log UV, respectively. However, nTiO₂ did not affect bacterial viability. Further, due to oxidative stress, the triacylglycerol (biodiesel) content from bacteria increased from 30% to a maximum of 54% CDW. Based on our findings, we propose an application of R. opacus PD 630 in nTiO₂ remediation due to their high nTiO₂ uptake and resistance.

High-resolution 13C metabolic flux analysis

Precise quantification of metabolic pathway fluxes in biological systems is of major importance in guiding efforts in metabolic engineering, biotechnology, microbiology, human health, and cell culture. ¹³C metabolic flux analysis (¹³C-MFA) is the predominant technique used for determining intracellular fluxes. Here, we present a protocol for ¹³C-MFA that incorporates recent advances in parallel labeling experiments, isotopic labeling measurements, and statistical analysis, as well as best practices developed through decades of experience. Experimental design to ensure that fluxes are estimated with the highest precision is an integral part of the protocol. The protocol is based on growing microbes in two (or more) parallel cultures with ¹³C-labeled glucose tracers, followed by gas chromatography-mass spectrometry (GC-MS) measurements of isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose. Fluxes are then estimated using software for ¹³C-MFA, such as Metran, followed by comprehensive statistical analysis to determine the goodness of fit and calculate confidence intervals of fluxes. The presented protocol can be completed in 4 d and quantifies metabolic fluxes with a standard deviation of ≤2%, a substantial improvement over previous implementations. The presented protocol is exemplified using an Escherichia coli ΔtpiA case study with full supporting data, providing a hands-on opportunity to step through a complex troubleshooting scenario. Although applications to prokaryotic microbial systems are emphasized, this protocol can be easily adjusted for application to eukaryotic organisms.

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.

Construction of Genetic Logic Gates Based on the T7 RNA Polymerase Expression System in Rhodococcus opacus PD630

Rhodococcus opacus PD630 (R. opacus) is a nonmodel, Gram-positive bacterium that holds promise as a biological catalyst for the conversion of lignocellulosic biomass to value-added products. In particular, it demonstrates both a high tolerance for and an ability to consume inhibitory lignin-derived aromatics, generates large quantities of lipids, exhibits a relatively rapid growth rate, and has a growing genetic toolbox for engineering. However, the availability of genetic parts for tunable, high-activity gene expression is still limited in R. opacus. Furthermore, genetic logic circuits for sophisticated gene regulation have never been demonstrated in Rhodococcus spp. To address these shortcomings, two inducible T7 RNA polymerase-based expression systems were implemented for the first time in R. opacus and applied to the construction of AND and NAND genetic logic gates. Additionally, three isopropyl β-d-1-thiogalactopyranoside (IPTG)-inducible promoters were created by inserting LacI binding sites into newly characterized constitutive promoters. Furthermore, four novel aromatic sensors for 4-hydroxybenzoic acid, vanillic acid, sodium benzoate, and guaiacol were developed, expanding the gene expression toolbox. Finally, the T7 RNA polymerase platform was combined with a synthetic IPTG-inducible promoter to create an IMPLY logic gate. Overall, this work represents the first demonstration of a heterologous RNA polymerase system and synthetic genetic logic in R. opacus, enabling complex and tunable gene regulation in this promising nonmodel host for bioproduction.

Mapping Salmonella typhimurium pathways using 13C metabolic flux analysis

In the last years, Salmonella has been extensively studied not only due to its importance as a pathogen, but also as a host to produce pharmaceutical compounds. However, the full exploitation of Salmonella as a platform for bioproduct delivery has been hampered by the lack of information about its metabolism. Genome-scale metabolic models can be valuable tools to delineate metabolic engineering strategies as long as they closely represent the actual metabolism of the target organism. In the present study, a ¹³C-MFA approach was applied to map the fluxes at the central carbon pathways of S. typhimurium LT2 growing at glucose-limited chemostat cultures. The experiments were carried out in a 2L bioreactor, using defined medium enriched with 20% ¹³C-labeled glucose. Metabolic flux distributions in central carbon pathways of S. typhimurium LT2 were estimated using OpenFLUX2 based on the labeling pattern of biomass protein hydrolysates together with biomass composition. The results suggested that pentose phosphate is used to catabolize glucose, with minor fluxes through glycolysis. In silico simulations, using Optflux and pFBA as simulation method, allowed to study the performance of the genome-scale metabolic model. In general, the accuracy of in silico simulations was improved by the superimposition of estimated intracellular fluxes to the existing genome-scale metabolic model, showing a better fitting to the experimental extracellular fluxes, whereas the intracellular fluxes of pentose phosphate and anaplerotic reactions were poorly described.

Development of Rhodococcus opacus as a chassis for lignin valorization and bioproduction of high-value compounds

The current extraction and use of fossil fuels has been linked to extensive negative health and environmental outcomes. Lignocellulosic biomass-derived biofuels and bioproducts are being actively considered as renewable alternatives to the fuels, chemicals, and materials produced from fossil fuels. A major challenge limiting large-scale, economic deployment of second-generation biorefineries is the insufficient product yield, diversity, and value that current conversion technologies can extract from lignocellulose, in particular from the underutilized lignin fraction. Rhodococcus opacus PD630 is an oleaginous gram-positive bacterium with innate catabolic pathways and tolerance mechanisms for the inhibitory aromatic compounds found in depolymerized lignin, as well as native or engineered pathways for hexose and pentose sugars found in the carbohydrate fractions of biomass. As a result, R. opacus holds potential as a biological chassis for the conversion of lignocellulosic biomass into biodiesel precursors and other value-added products. This review begins by examining the important role that lignin utilization will play in the future of biorefineries and by providing a concise survey of the current lignin conversion technologies. The genetic machinery and capabilities of R. opacus that allow the bacterium to tolerate and metabolize aromatic compounds and depolymerized lignin are also discussed, along with a synopsis of the genetic toolbox and synthetic biology methods now available for engineering this organism. Finally, we summarize the different feedstocks that R. opacus has been demonstrated to consume, and the high-value products that it has been shown to produce. Engineered R. opacus will enable lignin valorization over the coming years, leading to cost-effective conversion of lignocellulose into fuels, chemicals, and materials.

13C-Fingerprinting and Metabolic Flux Analysis of Bacterial Metabolisms

¹³C-assisted metabolism analysis provides rigorous calculations of the intracellular reaction rates (i.e., fluxes) within the central metabolism of microbial hosts. This mapping of the intracellular network within microbes has proven to be essential for understanding the cell physiology. The approach is also a key to identifying central metabolic nodes, probing the rigidity of a metabolic network, revealing cofactor balances, and delineating hidden pathways. Here we present the methodology of using stable isotopic carbon substrates for both qualitative (¹³C-fingerprinting of functional pathways) and quantitative (Metabolic Flux Analysis) metabolism studies on bacterial species. In this methodology, we include step-by-step instructions to use the open source WUflux software for the steady-state flux calculations based on labeling information of amino acids or free metabolites.

Discovery of potential pathways for biological conversion of poplar wood into lipids by co-fermentation of Rhodococci strains

BACKGROUND

Biological routes for utilizing both carbohydrates and lignin are important to reach the ultimate goal of bioconversion of full carbon in biomass into biofuels and biochemicals. Recent biotechnology advances have shown promises toward facilitating biological transformation of lignin into lipids. Natural and engineered Rhodococcus strains (e.g., R. opacus PD630, R. jostii RHA1, and R. jostii RHA1 VanA^-) have been demonstrated to utilize lignin for lipid production, and co-culture of them can promote lipid production from lignin.

RESULTS

In this study, a co-fermentation module of natural and engineered Rhodococcus strains with significant improved lignin degradation and/or lipid biosynthesis capacities was established, which enabled simultaneous conversion of glucose, lignin, and its derivatives into lipids. Although Rhodococci sp. showed preference to glucose over lignin, nearly half of the lignin was quickly depolymerized to monomers by these strains for cell growth and lipid synthesis after glucose was nearly consumed up. Profiles of metabolites produced by Rhodococcus strains growing on different carbon sources (e.g., glucose, alkali lignin, and dilute acid flowthrough-pretreated poplar wood slurry) confirmed lignin conversion during co-fermentation, and indicated novel metabolic capacities and unexplored metabolic pathways in these organisms. Proteome profiles suggested that lignin depolymerization by Rhodococci sp. involved multiple peroxidases with accessory oxidases. Besides the β-ketoadipate pathway, the phenylacetic acid (PAA) pathway was another potential route for the in vivo ring cleavage activity. In addition, deficiency of reducing power and cellular oxidative stress probably led to lower lipid production using lignin as the sole carbon source than that using glucose.

CONCLUSIONS

This work demonstrated a potential strategy for efficient bioconversion of both lignin and glucose into lipids by co-culture of multiple natural and engineered Rhodococcus strains. In addition, the involvement of PAA pathway in lignin degradation can help to further improve lignin utilization, and the combinatory proteomics and bioinformatics strategies used in this study can also be applied into other systems to reveal the metabolic and regulatory pathways for balanced cellular metabolism and to select genetic targets for efficient conversion of both lignin and carbohydrates into biofuels.

Lipid metabolism of phenol-tolerant Rhodococcus opacus strains for lignin bioconversion

BACKGROUND

Lignin is a recalcitrant aromatic polymer that is a potential feedstock for renewable fuel and chemical production. Rhodococcus opacus PD630 is a promising strain for the biological upgrading of lignin due to its ability to tolerate and utilize lignin-derived aromatic compounds. To enhance its aromatic tolerance, we recently applied adaptive evolution using phenol as a sole carbon source and characterized a phenol-adapted R. opacus strain (evol40) and the wild-type (WT) strain by whole genome and RNA sequencing. While this effort increased our understanding of the aromatic tolerance, the tolerance mechanisms were not completely elucidated.

RESULTS

We hypothesize that the composition of lipids plays an important role in phenol tolerance. To test this hypothesis, we applied high-resolution mass spectrometry analysis to lipid samples obtained from the WT and evol40 strains grown in 1 g/L glucose (glucose), 0.75 g/L phenol (low phenol), or 1.5 g/L phenol (high phenol, evol40 only) as a sole carbon source. This analysis identified > 100 lipid species of mycolic acids, phosphatidylethanolamines (PEs), phosphatidylinositols (PIs), and triacylglycerols. In both strains, mycolic acids had fewer double bond numbers in phenol conditions than the glucose condition, and evol40 had significantly shorter mycolic acid chain lengths than the WT strain in phenol conditions. These results indicate that phenol adaptation affected mycolic acid membrane composition. In addition, the percentage of unsaturated phospholipids decreased for both strains in phenol conditions compared to the glucose condition. Moreover, the PI content increased for both strains in the low phenol condition compared to the glucose condition, and the PI content increased further for evol40 in the high phenol condition relative to the low phenol condition.

CONCLUSIONS

This work represents the first comprehensive lipidomic study on the membrane of R. opacus grown using phenol as a sole carbon source. Our results suggest that the alteration of the mycolic acid and phospholipid membrane composition may be a strategy of R. opacus for phenol tolerance.

Selection of stable reference genes for RT-qPCR in Rhodococcus opacus PD630

Rhodococcus opacus PD630 is a gram-positive bacterium with promising attributes for the conversion of lignin into valuable fuels and chemicals. To develop an organism as a cellular factory, it is necessary to have a deep understanding of its metabolism and any heterologous pathways being expressed. For the purpose of quantifying gene transcription, reverse transcription quantitative PCR (RT-qPCR) is the gold standard due to its sensitivity and reproducibility. However, RT-qPCR requires the use of reference genes whose expression is stable across distinct growth or treatment conditions to normalize the results. Unfortunately, no in-depth analysis of stable reference genes has been conducted in Rhodococcus, inhibiting the utilization of RT-qPCR in R. opacus. In this work, ten candidate reference genes, chosen based on previously collected RNA sequencing data or literature, were examined under four distinct growth conditions using three mathematical programs (BestKeeper, Normfinder, and geNorm). Based on this analysis, the minimum number of reference genes required was found to be two, and two separate pairs of references genes were identified as optimal normalization factors for when ribosomal RNA is either present or depleted. This work represents the first validation of reference genes for Rhodococcus, providing a valuable starting point for future research.

Molecular Toolkit for Gene Expression Control and Genome Modification in Rhodococcus opacus PD630

Rhodococcus opacus PD630 is a non-model Gram-positive bacterium that possesses desirable traits for lignocellulosic biomass conversion. In particular, it has a relatively rapid growth rate, exhibits genetic tractability, produces high quantities of lipids, and can tolerate and consume toxic lignin-derived aromatic compounds. Despite these unique, industrially relevant characteristics, R. opacus has been underutilized because of a lack of reliable genetic parts and engineering tools. In this work, we developed a molecular toolbox for reliable gene expression control and genome modification in R. opacus. To facilitate predictable gene expression, a constitutive promoter library spanning ∼45-fold in output was constructed. To improve the characterization of available plasmids, the copy numbers of four heterologous and nine endogenous plasmids were determined using quantitative PCR. The molecular toolbox was further expanded by screening a previously unreported antibiotic resistance marker (HygR) and constructing a curable plasmid backbone for temporary gene expression (pB264). Furthermore, a system for genome modification was devised, and three neutral integration sites were identified using a novel combination of transcriptomic data, genomic architecture, and growth rate analysis. Finally, the first reported system for targeted, tunable gene repression in Rhodococcus was developed by utilizing CRISPR interference (CRISPRi). Overall, this work greatly expands the ability to manipulate and engineer R. opacus, making it a viable new chassis for bioproduction from renewable feedstocks.

Metabolism of aromatics by Trichosporon oleaginosus while remaining oleaginous

BACKGROUND

The oleaginous yeast, Trichosporon oleaginosus, has been extensively studied for its ability to metabolize non-conventional feedstocks. These include phenol-containing waste streams, such as distillery wastewater, or streams consisting of non-conventional sugars, such as hydrolyzed biomass and various bagasse. An initial BLAST search suggests this yeast has putative aromatic metabolizing genes. Given the desirability to valorize underutilized feedstocks such as lignin, we investigated the ability of T. oleaginosus to tolerate and metabolize lignin-derived aromatic compounds.

RESULTS

Trichosporon oleaginosus can tolerate and metabolize model lignin monoaromatics and associated intermediates within funneling pathways. Growth rates and biomass yield were similar to glucose when grown in 4-hydroxybenzoic acid (pHBA) and resorcinol, but had an increased lag phase when grown in phenol. Oleaginous behavior was observed using resorcinol as a sole carbon source. Fed-batch feeding resulted in lipid accumulation of 69.5% on a dry weight basis.

CONCLUSIONS

Though the exact pathway of aromatic metabolism remains to be determined for T. oleaginosus, the results presented in this work motivate use of this organism for lignin valorization and phenolic wastewater bioremediation. Trichosporon oleaginosus is the first yeast shown to be oleaginous while growing on aromatic substrates, and shows great promise as a model industrial microbe for biochemical and biofuel production from depolymerized lignin.

Development of Chemical and Metabolite Sensors for Rhodococcus opacus PD630

Rhodococcus opacus PD630 is a nonmodel, Gram-positive bacterium that possesses desirable traits for biomass conversion, including consumption capabilities for lignocellulose-based sugars and toxic lignin-derived aromatic compounds, significant triacylglycerol accumulation, relatively rapid growth rate, and genetic tractability. However, few genetic elements have been directly characterized in R. opacus, limiting its application for lignocellulose bioconversion. Here, we report the characterization and development of genetic tools for tunable gene expression in R. opacus, including: (1) six fluorescent reporters for quantifying promoter output, (2) three chemically inducible promoters for variable gene expression, and (3) two classes of metabolite sensors derived from native R. opacus promoters that detect nitrogen levels or aromatic compounds. Using these tools, we also provide insights into native aromatic consumption pathways in R. opacus. Overall, this work expands the ability to control and characterize gene expression in R. opacus for future lignocellulose-based fuel and chemical production.

Synthetic biology for manufacturing chemicals: constraints drive the use of non-conventional microbial platforms

Genetically modified microbes have had much industrial success producing protein-based products (such as antibodies and enzymes). However, engineering microbial workhorses for biomanufacturing of commodity compounds remains challenging. First, microbes cannot afford burdens with both overexpression of multiple enzymes and metabolite drainage for product synthesis. Second, synthetic circuits and introduced heterologous pathways are not yet as "robust and reliable" as native pathways due to hosts' innate regulations, especially under suboptimal fermentation conditions. Third, engineered enzymes may lack channeling capabilities for cascade-like transport of metabolites to overcome diffusion barriers or to avoid intermediate toxicity in the cytoplasmic environment. Fourth, moving engineered hosts from laboratory to industry is unreliable because genetic mutations and non-genetic cell-to-cell variations impair the large-scale fermentation outcomes. Therefore, synthetic biology strains often have unsatisfactory industrial performance (titer/yield/productivity). To overcome these problems, many different species are being explored for their metabolic strengths that can be leveraged to synthesize specific compounds. Here, we provide examples of non-conventional and genetically amenable species for industrial manufacturing, including the following: Corynebacterium glutamicum for its TCA cycle-derived biosynthesis, Yarrowia lipolytica for its biosynthesis of fatty acids and carotenoids, cyanobacteria for photosynthetic production from its sugar phosphate pathways, and Rhodococcus for its ability to biotransform recalcitrant feedstock. Finally, we discuss emerging technologies (e.g., genome-to-phenome mapping, single cell methods, and knowledge engineering) that may facilitate the development of novel cell factories.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.

Biological conversion of the aqueous wastes from hydrothermal liquefaction of algae and pine wood by Rhodococci

In this study, R. opacus PD630, R. jostii RHA1, R. jostii RHA1 VanA^-, and their co-culture were employed to convert hydrothermal liquefaction aqueous waste (HTLAW) into lipids. After 11days, the COD reduction of algal-HTLAW reached 93.4% and 92.7% by R. jostii RHA1 and its mutant VanA^-, respectively. Woody-HTLAW promoted lipid accumulation of 0.43glipid/gcell dry weight in R. opacus PD630 cells. Additionally, the total number of chemicals in HTLAW decreased by over 1/3 after 7days of coculture, and 0.10g/L and 0.46g/L lipids were incrementally accumulated in the cellular mass during the fermentation of wood- and algal-HTLAW, respectively. The GC-MS data supported that different metabolism pathways were followed when these Rhodococci strains degraded algae- and woody-HTLAW. These results indicated promising potential of bioconversion of under-utilized carbon and toxic compounds in HTLAW into useful products by selected Rhodococci.

Recent applications of metabolomics to advance microbial biofuel production

Biofuel production from plant biomass is a promising source of renewable energy [1]. However, efficient biofuel production involves the complex task of engineering high-performance microorganisms, which requires detailed knowledge of metabolic function and regulation. This review highlights the potential of mass-spectrometry-based metabolomic analysis to guide rational engineering of biofuel-producing microbes. We discuss recent studies that apply knowledge gained from metabolomic analyses to increase the productivity of engineered pathways, characterize the metabolism of emerging biofuel producers, generate novel bioproducts, enable utilization of lignocellulosic feedstock, and improve the stress tolerance of biofuel producers.

Recent advances in high-throughput 13C-fluxomics

The rise of high throughput (HT) strain engineering tools accompanying the area of synthetic biology is supporting the generation of a large number of microbial cell factories. A current bottleneck in process development is our limited capacity to rapidly analyze the metabolic state of the engineered strains, and in particular their intracellular fluxes. HT ¹³C-fluxomics workflows have not yet become commonplace, despite the existence of several HT tools at each of the required stages. This includes cultivation and sampling systems, analytics for isotopic analysis, and software for data processing and flux calculation. Here, we review recent advances in the field and highlight bottlenecks that must be overcome to allow the emergence of true HT ¹³C-fluxomics workflows.

WUFlux: an open-source platform for 13C metabolic flux analysis of bacterial metabolism

Background

Flux analyses, including flux balance analysis (FBA) and ¹³C-metabolic flux analysis (¹³C-MFA), offer direct insights into cell metabolism, and have been widely used to characterize model and non-model microbial species. Nonetheless, constructing the ¹³C-MFA model and performing flux calculation are demanding for new learners, because they require knowledge of metabolic networks, carbon transitions, and computer programming. To facilitate and standardize the ¹³C-MFA modeling work, we set out to publish a user-friendly and programming-free platform (WUFlux) for flux calculations in MATLAB^®.

Results

We constructed an open-source platform for steady-state ¹³C-MFA. Using GUIDE (graphical user interface design environment) in MATLAB, we built a user interface that allows users to modify models based on their own experimental conditions. WUFlux is capable of directly correcting mass spectrum data of TBDMS (N-tert-butyldimethylsilyl-N-methyltrifluoroacetamide)-derivatized proteinogenic amino acids by removing background noise. To simplify ¹³C-MFA of different prokaryotic species, the software provides several metabolic network templates, including those for chemoheterotrophic bacteria and mixotrophic cyanobacteria. Users can modify the network and constraints, and then analyze the microbial carbon and energy metabolisms of various carbon substrates (e.g., glucose, pyruvate/lactate, acetate, xylose, and glycerol). WUFlux also offers several ways of visualizing the flux results with respect to the constructed network. To validate our model’s applicability, we have compared and discussed the flux results obtained from WUFlux and other MFA software. We have also illustrated how model constraints of cofactor and ATP balances influence fluxome results.

Conclusion

Open-source software for ¹³C-MFA, WUFlux, with a user-friendly interface and easy-to-modify templates, is now available at http://www.13cmfa.org/or (http://tang.eece.wustl.edu/ToolDevelopment.htm). We will continue documenting curated models of non-model microbial species and improving WUFlux performance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1314-0) contains supplementary material, which is available to authorized users.

Engineering levoglucosan metabolic pathway in Rhodococcus jostii RHA1 for lipid production

Oleaginous strains of Rhodococcus including R. jostii RHA1 have attracted considerable attention due to their ability to accumulate triacylglycerols (TAGs), robust growth properties and genetic tractability. In this study, a novel metabolic pathway was introduced into R. jostii by heterogenous expression of the well-characterized gene, lgk encoding levoglucosan kinase from Lipomyces starkeyi YZ-215. This enables the recombinant R. jostii RHA1 to produce TAGs from the anhydrous sugar, levoglucosan, which can be generated efficiently as the major molecule from the pyrolysis of cellulose. The recombinant R. jostii RHA1 could grow on levoglucosan as the sole carbon source, and the consumption rate of levoglucosan was determined. Furthermore, expression of one more copy of lgk increased the enzymatic activity of LGK in the recombinant. However, the growth performance of the recombinant bearing two copies of lgk on levoglucosan was not improved. Although expression of lgk in the recombinants was not repressed by the glucose present in the media, glucose in the sugar mixture still affected consumption of levoglucosan. Under nitrogen limiting conditions, lipid produced from levoglucosan by the recombinant bearing lgk was up to 43.54 % of the cell dry weight, which was comparable to the content of lipid accumulated from glucose. This work demonstrated the technical feasibility of producing lipid from levoglucosan, an anhydrosugar derived from the pyrolysis of lignocellulosic materials, by the genetically modified rhodococci strains.

Current state and challenges for dynamic metabolic modeling

While the stoichiometry of metabolism is probably the best studied cellular level, the dynamics in metabolism can still not be well described, predicted and, thus, engineered. Unknowns in the metabolic flux behavior arise from kinetic interactions, especially allosteric control mechanisms. While the stoichiometry of enzymes is preserved in vitro, their activity and kinetic behavior differs from the in vivo situation. Next to this challenge, it is infeasible to test the interaction of each enzyme with each intracellular metabolite in vitro exhaustively. As a consequence, the whole interacting metabolome has to be studied in vivo to identify the relevant enzymes properties. In this review we discuss current approaches for in vivo perturbation experiments, that is, stimulus response experiments using different setups and quantitative analytical approaches, including dynamic carbon tracing. Next to reliable and informative data, advanced modeling approaches and computational tools are required to identify kinetic mechanisms and their parameters.

Optimal tracers for parallel labeling experiments and 13C metabolic flux analysis: A new precision and synergy scoring system

¹³C-Metabolic flux analysis (¹³C-MFA) is a widely used approach in metabolic engineering for quantifying intracellular metabolic fluxes. The precision of fluxes determined by ¹³C-MFA depends largely on the choice of isotopic tracers and the specific set of labeling measurements. A recent advance in the field is the use of parallel labeling experiments for improved flux precision and accuracy. However, as of today, no systemic methods exist for identifying optimal tracers for parallel labeling experiments. In this contribution, we have addressed this problem by introducing a new scoring system and evaluating thousands of different isotopic tracer schemes. Based on this extensive analysis we have identified optimal tracers for ¹³C-MFA. The best single tracers were doubly ¹³C-labeled glucose tracers, including [1,6-¹³C]glucose, [5,6-¹³C]glucose and [1,2-¹³C]glucose, which consistently produced the highest flux precision independent of the metabolic flux map (here, 100 random flux maps were evaluated). Moreover, we demonstrate that pure glucose tracers perform better overall than mixtures of glucose tracers. For parallel labeling experiments the optimal isotopic tracers were [1,6-¹³C]glucose and [1,2-¹³C]glucose. Combined analysis of [1,6-¹³C]glucose and [1,2-¹³C]glucose labeling data improved the flux precision score by nearly 20-fold compared to widely use tracer mixture 80% [1-¹³C]glucose +20% [U-¹³C]glucose.

On the Kinetic and Allosteric Regulatory Properties of the ADP-Glucose Pyrophosphorylase from Rhodococcus jostii: An Approach to Evaluate Glycogen Metabolism in Oleaginous Bacteria

Rhodococcus spp. are oleaginous bacteria that accumulate glycogen during exponential growth. Despite the importance of these microorganisms in biotechnology, little is known about the regulation of carbon and energy storage, mainly the relationship between glycogen and triacylglycerols metabolisms. Herein, we report the molecular cloning and heterologous expression of the gene coding for ADP-glucose pyrophosphorylase (EC 2.7.7.27) of Rhodococcus jostii, strain RHA1. The recombinant enzyme was purified to electrophoretic homogeneity to accurately characterize its oligomeric, kinetic, and regulatory properties. The R. jostii ADP-glucose pyrophosphorylase is a homotetramer of 190 kDa exhibiting low basal activity to catalyze synthesis of ADP-glucose, which is markedly influenced by different allosteric effectors. Glucose-6P, mannose-6P, fructose-6P, ribose-5P, and phosphoenolpyruvate were major activators; whereas, NADPH and 6P-gluconate behaved as main inhibitors of the enzyme. The combination of glucose-6P and other effectors (activators or inhibitors) showed a cross-talk effect suggesting that the different metabolites could orchestrate a fine regulation of ADP-glucose pyrophosphorylase in R. jostii. The enzyme exhibited some degree of affinity toward ATP, GTP, CTP, and other sugar-1P substrates. Remarkably, the use of glucosamine-1P was sensitive to allosteric activation. The relevance of the fine regulation of R. jostii ADP-glucose pyrophosphorylase is further analyzed in the framework of proteomic studies already determined for the bacterium. Results support a critical role for glycogen as a temporal reserve that provides a pool of carbon able of be re-routed to produce long-term storage of lipids under certain conditions.

Engineering of an L-arabinose metabolic pathway in Rhodococcus jostii RHA1 for biofuel production

The oleaginous bacterium, Rhodococcus jostii RHA1 has attracted considerable attention due to its capability to accumulate significant levels of triacylglycerol as renewable hydrocarbon. To enable the strain to utilize arabinose derived from lignocellulosic biomass, the metabolic pathway of L-arabinose utilization was introduced into R. jostii RHA1 by heterogenous expression of the operon, araBAD from Escherichia coli. The results showed that recombinant bearing araBAD could grow on L-arabinose as the sole carbon source, and additional expression of araFGH encoding the arabinose transporter from E. coli could improve the cell biomass yield from high contents of arabinose. We further increased the content of lipid produced from arabinose in the recombinants from 47.9 to 56.8 % of the cell dry weight (CDW) by overexpression of a gene, atf1 encoding a diglyceride acyltransferase from R. opacus PD630. This work demonstrated the feasibility of producing lipid from arabinose by genetic modification of the rhodococci strain.

Comparison of gas chromatography-mass spectrometry and gas chromatography-combustion-isotope ratio mass spectrometry analysis for in vivo estimates of metabolic fluxes

Gas chromatography-mass spectrometry (GC-MS) was compared with gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS) for measurements of cholesterol (13)C enrichment after infusion of labeled precursor ([(13)C1,2]acetate). Paired results were significantly correlated, although GC-MS was less accurate than GC-C-IRMS for higher enrichments. Nevertheless, only GC-MS was able to provide information on isotopologue distribution, bringing new insights to lipid metabolism. Therefore, we assessed the isotopologue distribution of cholesterol in humans and dogs known to present contrasted cholesterol metabolic pathways. The labeled tracer incorporation was different in both species, highlighting the subsidiarity of GC-MS and GC-C-IRMS to analyze in vivo stable isotope studies.

Comparative transcriptomics elucidates adaptive phenol tolerance and utilization in lipid-accumulating Rhodococcus opacus PD630

Lignin-derived (e.g. phenolic) compounds can compromise the bioconversion of lignocellulosic biomass to fuels and chemicals due to their toxicity and recalcitrance. The lipid-accumulating bacterium Rhodococcus opacus PD630 has recently emerged as a promising microbial host for lignocellulose conversion to value-added products due to its natural ability to tolerate and utilize phenolics. To gain a better understanding of its phenolic tolerance and utilization mechanisms, we adaptively evolved R. opacus over 40 passages using phenol as its sole carbon source (up to 373% growth improvement over wild-type), and extensively characterized two strains from passages 33 and 40. The two adapted strains showed higher phenol consumption rates (∼20 mg/l/h) and ∼2-fold higher lipid production from phenol than the wild-type strain. Whole-genome sequencing and comparative transcriptomics identified highly-upregulated degradation pathways and putative transporters for phenol in both adapted strains, highlighting the important linkage between mechanisms of regulated phenol uptake, utilization, and evolved tolerance. Our study shows that the R. opacus mutants are likely to use their transporters to import phenol rather than export them, suggesting a new aromatic tolerance mechanism. The identified tolerance genes and pathways are promising candidates for future metabolic engineering in R. opacus for improved lignin conversion to lipid-based products.

Elucidation of the co-metabolism of glycerol and glucose in Escherichia coli by genetic engineering, transcription profiling, and (13)C metabolic flux analysis

BACKGROUND

Glycerol, a byproduct of biodiesel, has become a readily available and inexpensive carbon source for the production of high-value products. However, the main drawback of glycerol utilization is the low consumption rate and shortage of NADPH formation, which may limit the production of NADPH-requiring products. To overcome these problems, we constructed a carbon catabolite repression-negative ΔptsGglpK* mutant by both blocking a key glucose PTS transporter and enhancing the glycerol conversion. The mutant can recover normal growth by co-utilization of glycerol and glucose after loss of glucose PTS transporter. To reveal the metabolic potential of the ΔptsGglpK* mutant, this study examined the flux distributions and regulation of the co-metabolism of glycerol and glucose in the mutant.

RESULTS

By labeling experiments using [1,3-(13)C]glycerol and [1-(13)C]glucose, (13)C metabolic flux analysis was employed to decipher the metabolisms of both the wild-type strain and the ΔptsGglpK* mutant in chemostat cultures. When cells were maintained at a low dilution rate (0.1 h(-1)), the two strains showed similar fluxome profiles. When the dilution rate was increased, both strains upgraded their pentose phosphate pathway, glycolysis and anaplerotic reactions, while the ΔptsGglpK* mutant was able to catabolize much more glycerol than glucose (more than tenfold higher). Compared with the wild-type strain, the mutant repressed its flux through the TCA cycle, resulting in higher acetate overflow. The regulation of fluxomes was consistent with transcriptional profiling of several key genes relevant to the TCA cycle and transhydrogenase, namely gltA, icdA, sdhA and pntA. In addition, cofactor fluxes and their pool sizes were determined. The ΔptsGglpK* mutant affected the redox NADPH/NADH state and reduced the ATP level. Redox signaling activated the ArcA regulatory system, which was responsible for TCA cycle repression.

CONCLUSIONS

This work employs both (13)C-MFA and transcription/metabolite analysis for quantitative investigation of the co-metabolism of glycerol and glucose in the ΔptsGglpK* mutant. The ArcA regulatory system dominates the control of flux redistribution. The ΔptsGglpK* mutant can be used as a platform for microbial cell factories for the production of biofuels and biochemicals, since most of fuel molecule (e.g., alcohols) synthesis requires excess reducing equivalents.

Examining Escherichia coli glycolytic pathways, catabolite repression, and metabolite channeling using Δpfk mutants

BACKGROUND

Glycolysis breakdowns glucose into essential building blocks and ATP/NAD(P)H for the cell, occupying a central role in its growth and bio-production. Among glycolytic pathways, the Entner Doudoroff pathway (EDP) is a more thermodynamically favorable pathway with fewer enzymatic steps than either the Embden-Meyerhof-Parnas pathway (EMPP) or the oxidative pentose phosphate pathway (OPPP). However, Escherichia coli do not use their native EDP for glucose metabolism.

RESULTS

Overexpression of edd and eda in E. coli to enhance EDP activity resulted in only a small shift in the flux directed through the EDP (~20 % of glycolysis flux). Disrupting the EMPP by phosphofructokinase I (pfkA) knockout increased flux through OPPP (~60 % of glycolysis flux) and the native EDP (~14 % of glycolysis flux), while overexpressing edd and eda in this ΔpfkA mutant directed ~70 % of glycolytic flux through the EDP. The downregulation of EMPP via the pfkA deletion significantly decreased the growth rate, while EDP overexpression in the ΔpfkA mutant failed to improve its growth rates due to metabolic burden. However, the reorganization of E. coli glycolytic strategies did reduce glucose catabolite repression. The ΔpfkA mutant in glucose medium was able to cometabolize acetate via the citric acid cycle and gluconeogenesis, while EDP overexpression in the ΔpfkA mutant repressed acetate flux toward gluconeogenesis. Moreover, ¹³C-pulse experiments in the ΔpfkA mutants showed unsequential labeling dynamics in glycolysis intermediates, possibly suggesting metabolite channeling (metabolites in glycolysis are pass from enzyme to enzyme without fully equilibrating within the cytosol medium).

CONCLUSIONS

We engineered E. coli to redistribute its native glycolytic flux. The replacement of EMPP by EDP did not improve E. coli glucose utilization or biomass growth, but alleviated catabolite repression. More importantly, our results supported the hypothesis of channeling in the glycolytic pathways, a potentially overlooked mechanism for regulating glucose catabolism and coutilization of other substrates. The presence of channeling in native pathways, if proven true, would affect synthetic biology applications and metabolic modeling.

13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production

Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many ¹³C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms