Comparative metabolic profiling reveals the key role of amino acids metabolism in the rapamycin overproduction by Streptomyces hygroscopicus

Integrated Metabolomic and Network Analysis to Explore the Potential Mechanism of Three Chemical Elicitors in Rapamycin Overproduction

Rapamycin is a polyketide macrocyclic antibiotic with exceptional pharmacological potential. To explore the potential mechanism of rapamycin overproduction, the intracellular metabolic differences of three chemical elicitor treatments were first investigated by combining them with dynamic metabolomics and network analysis. The metabolic response characteristics of each chemical elicitor treatment were identified by a weighted gene co-expression network analysis (WGCNA) model. According to the analysis of the identified metabolic modules, the changes in the cell membrane permeability might play a key role in rapamycin overproduction for dimethyl sulfoxide (DMSO) treatment. The enhancement of the starter unit of 4,5-dihydroxycyclohex-1-ene carboxylic acid (DHCHC) and the nicotinamide adenine dinucleotide phosphate (NADPH) availability were the main functions in the LaCl3 treatment. However, for sodium butyrate (SB), the improvement of the methylmalonyl-CoA and NADPH availability was a potential reason for the rapamycin overproduction. Further, the responsive metabolic pathways after chemical elicitor treatments were selected to predict the potential key limiting steps in rapamycin accumulation using a genome-scale metabolic network model (GSMM). Based on the prediction results, the targets within the reinforcement of the DHCHC and NADPH supply were selected to verify their effects on rapamycin production. The highest rapamycin yield improved 1.62 fold in the HT-aroA/zwf2 strain compared to the control.

Evaluation and optimization of analytical procedure and sample preparation for polar Streptomyces albus J1074 metabolome profiling

Metabolomics is an essential discipline in omics technology that promotes research on the biology of microbial systems. Streptomyces albus J1074 is a model organism used in fundamental research and industrial microbiology. Nevertheless, a comprehensive and standardized method for analyzing the metabolome of S. albus J1074 is yet to be developed. Thus, we comprehensively evaluated and optimized the analytical procedure and sample preparation for profiling polar metabolites using hydrophilic interaction liquid chromatography (HILIC) coupled with high-resolution mass spectrometry (HRMS). We systematically examined the HILIC columns, quenching solutions, sample-to-quenching ratios, and extraction methods. Then, the optimal protocol was used to investigate the dynamic intracellular polar metabolite profile of the engineered S. albus J1074 strains during spinosad (spinosyn A and spinosyn D) fermentation. A total of 3648 compounds were detected, and 83 metabolites were matched to the standards. The intracellular metabolomic profiles of engineered S. albus J1074 strains (ADE-AP and OE3) were detected; furthermore, their metabolomes in different stages were analyzed to reveal the reasons for their differences in their spinosad production, as well as the current metabolic limitation of heterologous spinosad production in S. albus J1074. The HILIC-HRMS method is a valuable tool for investigating polar metabolomes, and provides a reference methodology to study other Streptomyces metabolomes.

•

A HILIC-HRMS method was developed for polar metabolome profiling.

•

Sample preparation protocol for Streptomyces albus J1074 intracellular metabolites was studied for the first time.

•

This study revealed the possible reasons for different production of spinosad of engineered S. albus J1074 strains.

Comparative genomic analysis of Streptomyces rapamycinicus NRRL 5491 and its mutant overproducing rapamycin

Streptomyces rapamycinicus NRRL 5491 is a well-known producer of rapamycin, a secondary metabolite with useful bioactivities, including antifungal, antitumor, and immunosuppressive functions. For the enhanced rapamycin production, a rapamycin-overproducing strain SRMK07 was previously obtained as a result of random mutagenesis. To identify genomic changes that allowed the SRMK07 strain’s enhanced rapamycin production, genomes of the NRRL 5491 and SRMK07 strains were newly sequenced in this study. The resulting genome sequences of the wild-type and SRMK07 strains showed the size of 12.47 Mbp and 9.56 Mbp, respectively. Large deletions were observed at both end regions of the SRMK07 strain’s genome, which cover 17 biosynthetic gene clusters (BGCs) encoding secondary metabolites. Also, genes in a genomic region containing the rapamycin BGC were shown to be duplicated. Finally, comparative metabolic network analysis using these two strains’ genome-scale metabolic models revealed biochemical reactions with different metabolic fluxes, which were all associated with NADPH generation. Taken together, the genomic and computational approaches undertaken in this study suggest biological clues for the enhanced rapamycin production of the SRMK07 strain. These clues can also serve as a basis for systematic engineering of a production host for further enhanced rapamycin production.

Enhanced ascomycin production in Streptomyces hygroscopicus var. ascomyceticus by employing polyhydroxybutyrate as an intracellular carbon reservoir and optimizing carbon addition

BACKGROUND

Ascomycin is a multifunctional antibiotic produced by Streptomyces hygroscopicus var. ascomyceticus. As a secondary metabolite, the production of ascomycin is often limited by the shortage of precursors during the late fermentation phase. Polyhydroxybutyrate is an intracellular polymer accumulated by prokaryotic microorganisms. Developing polyhydroxybutyrate as an intracellular carbon reservoir for precursor synthesis is of great significance to improve the yield of ascomycin.

RESULTS

The fermentation characteristics of the parent strain S. hygroscopicus var. ascomyceticus FS35 showed that the accumulation and decomposition of polyhydroxybutyrate was respectively correlated with cell growth and ascomycin production. The co-overexpression of the exogenous polyhydroxybutyrate synthesis gene phaC and native polyhydroxybutyrate decomposition gene fkbU increased both the biomass and ascomycin yield. Comparative transcriptional analysis showed that the storage of polyhydroxybutyrate during the exponential phase accelerated biosynthesis processes by stimulating the utilization of carbon sources, while the decomposition of polyhydroxybutyrate during the stationary phase increased the biosynthesis of ascomycin precursors by enhancing the metabolic flux through primary pathways. The comparative analysis of cofactor concentrations confirmed that the biosynthesis of polyhydroxybutyrate depended on the supply of NADH. At low sugar concentrations found in the late exponential phase, the optimization of carbon source addition further strengthened the polyhydroxybutyrate metabolism by increasing the total concentration of cofactors. Finally, in the fermentation medium with 22 g/L starch and 52 g/L dextrin, the ascomycin yield of the co-overexpression strain was increased to 626.30 mg/L, which was 2.11-fold higher than that of the parent strain in the initial medium (296.29 mg/L).

CONCLUSIONS

Here we report for the first time that polyhydroxybutyrate metabolism is beneficial for cell growth and ascomycin production by acting as an intracellular carbon reservoir, stored as polymers when carbon sources are abundant and depolymerized into monomers for the biosynthesis of precursors when carbon sources are insufficient. The successful application of polyhydroxybutyrate in increasing the output of ascomycin provides a new strategy for improving the yields of other secondary metabolites.

Metabolomics analysis of Actinosynnema pretiosum with improved AP-3 production by enhancing UDP-glucose biosynthesis

Ansamitocin P-3 (AP-3) shows strong anticancer effects and has used as a payload for antibody-drug conjugates. Our previous study have shown that although genetically engineered Actinosynnema pretiosum strains with enhanced UDP-glucose (UDPG) biosynthesis displayed improved AP-3 production compared to the wild-type strain, the increase in yield was far from meeting the industrial demand. In this study, comparative metabolomics analysis complemented with quantitative real-time PCR analysis was performed for the wild-type strain and two mutants (OpgmOugp, ΔzwfΔgnd) to identify possible metabolic bottlenecks and non-intuitive targets for further enhancement of AP-3 production. We observed that enhancing intracellular UDPG availability facilitated the accumulation of intracellular N-demethyl-AP-3 and AP-3, where the transporting of them outside the cell still needs to be developed. We also found that the UDPG biosynthesis was closely associated with the availability of fructose in the medium and a suitable fructose feeding strategy could promote the further improvement of AP-3 titer. In addition, pathway abundance analysis revealed that undesired fatty acid accumulation and down-regulation of amino acid metabolism may be unfavorable for ansamitocin biosynthesis in later stage of production. These results indicate that genetic modification of the UDPG biosynthetic pathways may have pleiotropic effects on AP-3 production. Efforts must be made to eliminate these newly identified metabolic bottlenecks to boost AP-3 production in A. pretiosum.

Crystal structures and biochemical analyses of the bacterial arginine dihydrolase ArgZ suggests a "bond rotation" catalytic mechanism

A recently discovered ornithine-ammonia cycle (OAC) serves as a conduit in the nitrogen storage and remobilization machinery in cyanobacteria. The OAC involves an arginine catabolic reaction catalyzed by the arginine dihydrolase ArgZ whose catalytic mechanism is unknown. Here we determined the crystal structures at 1.2-3.0 Å of unliganded ArgZ from the cyanobacterium Synechocystis sp. PCC6803 and of ArgZ complexed with its substrate arginine, a covalently linked reaction intermediate, or the reaction product ornithine. The structures reveal that a key residue, Asn⁷¹, in the ArgZ active center functions as the determinant distinguishing ArgZ from other members of the guanidino group-modifying enzyme superfamily. The structures, along with biochemical evidence from enzymatic assays coupled with electrospray ionization MS techniques, further suggest that ArgZ-catalyzed conversion of arginine to ornithine, ammonia, and carbon dioxide consists of two successive cycles of amine hydrolysis. Finally, we show that arginine dihydrolases are broadly distributed among bacteria and metazoans, suggesting that the OAC may be frequently used for redistribution of nitrogen from arginine catabolism or nitrogen fixation.

Coupling metabolomics analysis and DOE optimization strategy towards enhanced IBDV production by chicken embryo fibroblast DF-1 cells

Infectious bursal disease (IBD) caused by IBD virus (IBDV) is highly contagious viral and vaccination in chicken embryo has been an effective mean to prevent acute infection. However, the current production of IBDV vaccine faces serious batch instability and external contamination. The chicken embryonic fibroblast cell line DF-1 is widely used for the proliferation of avian viruses and vaccine production. Thus, optimizing the production of IBDV by DF-1 cells has an important application value. Combining metabolomics analysis and a Design of Experiments (DOE) statistical strategy, this study successfully optimized the process of IBDV production by DF-1 cells. Differential analysis and time series analysis of metabolite data in both IBDV-infected and uninfected DF-1 cells were performed by multivariate statistical analysis. The results showed that the intracellular metabolite intensities of glycolysis, the pentose phosphate pathway, the nucleoside synthesis pathway, lipid metabolism, and glutathione metabolism were upregulated, and the TCA cycle underwent a slight downregulation after IBDV infection of DF-1 cells. Based on the metabolome results and DOE statistical optimization method, the additive components suitable for IBDV proliferation were determined. The IBDV titer increased by 20.7 times upon exogenous addition of cysteine, methionine, lysine and nucleosides in the control medium, which is consistent with the predicted result (20.0 times) by a multivariate quadratic equation. This study provides a strategy for the efficient production of IBDV vaccines and could potentially be utilized to improve the production of other viral vaccines and biologics.

Optimized sampling protocol for mass spectrometry-based metabolomics in Streptomyces

In quantitative metabolomics studies, the most crucial step was arresting snapshots of all interesting metabolites. However, the procedure customized for Streptomyces was so rare that most studies consulted the procedure from other bacteria even yeast, leading to inaccurate and unreliable metabolomics analysis. In this study, a base solution (acetone: ethanol = 1:1, mol/mol) was added to a quenching solution to keep the integrity of the cell membrane. Based on the molar transition energy (ET) of the organic solvents, five solutions were used to carry out the quenching procedures. These were acetone, isoamylol, propanol, methanol, and 60% (v/v) methanol. To the best of our knowledge, this is the first report which has utilized a quenching solution with ET values. Three procedures were also adopted for extraction. These were boiling, freezing–thawing, and grinding ethanol. Following the analysis of the mass balance, amino acids, organic acids, phosphate sugars, and sugar alcohols were measured using gas chromatography with an isotope dilution mass spectrometry. It was found that using isoamylol with a base solution (5:1, v/v) as a quenching solution and that freezing–thawing in liquid nitrogen within 50% (v/v) methanol as an extracting procedure were the best pairing for the quantitative metabolomics of Streptomyces ZYJ-6, and resulted in average recoveries of close to 100%. The concentration of intracellular metabolites obtained from this new quenching solution was between two and ten times higher than that from 60% (v/v) methanol, which until now has been the most commonly used solution. Our findings are the first systematic quantitative metabolomics tools for Streptomyces ZYJ-6 and, therefore, will be important references for research in fields such as 13C based metabolic flux analysis, multi-omic research and genome-scale metabolic model establishment, as well as for other Streptomyces.

Omics-based analyses revealed metabolic responses of Clostridium acetobutylicum to lignocellulose-derived inhibitors furfural, formic acid and phenol stress for butanol fermentation

BACKGROUND

Clostridium acetobutylicum is a model fermentative anaerobe for consolidated bioprocessing of lignocellulose hydrolysates into acetone-butanol-ethanol (ABE). However, the main inhibitors (acids, furans and phenols) ubiquitous in lignocellulose hydrolysates strictly limit the conversion efficiency. Thus, it is essential to understand the underlying mechanisms of lignocellulose hydrolysate inhibitors to identify key industrial bottlenecks that undermine efficient biofuel production. The recently developed omics strategy for intracellular metabolites and protein quantification now allow for the in-depth mapping of strain metabolism and thereby enable the resolution of the underlying mechanisms.

RESULTS

The toxicity of the main inhibitors in lignocellulose hydrolysates against C. acetobutylicum and ABE production was systematically investigated, and the changes in intracellular metabolism were analyzed by metabolomics and proteomics. The toxicity of the main lignocellulose hydrolysate inhibitors at the same dose was ranked as follows: formic acid > phenol > furfural. Metabolomic analysis based on weighted gene coexpression network analysis (WGCNA) revealed that the three inhibitors triggered the stringent response of C. acetobutylicum. Proteomic analysis based on peptide mass spectrometry (MS) supported the above results and provided more comprehensive conclusions. Under the stress of three inhibitors, the metabolites and key enzymes/proteins involved in glycolysis, reductive tricarboxylic acid (TCA) cycle, acetone-butanol synthesis and redox metabolism were lower than those in the control group. Moreover, proteins involved in gluconeogenesis, the oxidative TCA cycle, thiol peroxidase (TPX) for oxidative stress were significantly upregulated, indicating that inhibitor stress induced the stress response and metabolic regulation. In addition, the three inhibitors also showed stress specificity related to fatty acid synthesis, ATP synthesis, nucleic acid metabolism, nicotinic acid metabolism, cell wall synthesis, spore synthesis and flagellum synthesis and so on.

CONCLUSIONS

Integrated omics platforms provide insight into the cellular responses of C. acetobutylicum to cytotoxic inhibitors released during the deconstruction of lignocellulose. This insight allows us to fully improve the strain to adapt to a challenging culture environment, which will prove critical to the industrial efficacy of C. acetobutylicum.

Medium optimization based on comparative metabolomic analysis of chicken embryo fibroblast DF-1 cells

Chicken embryo fibroblast DF-1 cells are increasingly being used in the production of avian virus vaccines. However, the relatively low proliferative capacity does not meet the requirements of industrial production. In this study, we attempted to improve the proliferative capacity of DF-1 cells. The results of intracellular metabolomics showed that 28 types of metabolites could play roles in DF-1 cell growth based on the variance and timing analysis of intracellular metabolites from DF-1 cells grown in two media with distinct growth difference, DMEM/F12 (1 : 1) and DMEM. By examining the differences in the components in the two media, DOE was used to screen and optimize the growth medium for DF-1 cells. The maximum cell density was 40.72% higher, and the infectious bursal disease virus (IBDV) titer was 2.68 times higher, in the optimized medium than in the control. This study proposes a complete solution from metabolomics to media optimization.

This study proposes a medium optimal method based on DOE methodology and metabolomics of chicken embryo fibroblasts DF-1 cells.

Metabolic network model guided engineering ethylmalonyl-CoA pathway to improve ascomycin production in Streptomyces hygroscopicus var. ascomyceticus

Background

Ascomycin is a 23-membered polyketide macrolide with high immunosuppressant and antifungal activity. As the lower production in bio-fermentation, global metabolic analysis is required to further explore its biosynthetic network and determine the key limiting steps for rationally engineering. To achieve this goal, an engineering approach guided by a metabolic network model was implemented to better understand ascomycin biosynthesis and improve its production.

Results

The metabolic conservation of Streptomyces species was first investigated by comparing the metabolic enzymes of Streptomyces coelicolor A3(2) with those of 31 Streptomyces strains, the results showed that more than 72% of the examined proteins had high sequence similarity with counterparts in every surveyed strain. And it was found that metabolic reactions are more highly conserved than the enzymes themselves because of its lower diversity of metabolic functions than that of genes. The main source of the observed metabolic differences was from the diversity of secondary metabolism. According to the high conservation of primary metabolic reactions in Streptomyces species, the metabolic network model of Streptomyces hygroscopicus var. ascomyceticus was constructed based on the latest reported metabolic model of S. coelicolor A3(2) and validated experimentally. By coupling with flux balance analysis and using minimization of metabolic adjustment algorithm, potential targets for ascomycin overproduction were predicted. Since several of the preferred targets were highly associated with ethylmalonyl-CoA biosynthesis, two target genes hcd (encoding 3-hydroxybutyryl-CoA dehydrogenase) and ccr (encoding crotonyl-CoA carboxylase/reductase) were selected for overexpression in S. hygroscopicus var. ascomyceticus FS35. Both the mutants HA-Hcd and HA-Ccr showed higher ascomycin titer, which was consistent with the model predictions. Furthermore, the combined effects of the two genes were evaluated and the strain HA-Hcd-Ccr with hcd and ccr overexpression exhibited the highest ascomycin production (up to 438.95 mg/L), 1.43-folds improvement than that of the parent strain FS35 (305.56 mg/L).

Conclusions

The successful constructing and experimental validation of the metabolic model of S. hygroscopicus var. ascomyceticus showed that the general metabolic network model of Streptomyces species could be used to analyze the intracellular metabolism and predict the potential key limiting steps for target metabolites overproduction. The corresponding overexpression strains of the two identified genes (hcd and ccr) using the constructed model all displayed higher ascomycin titer. The strategy for yield improvement developed here could also be extended to the improvement of other secondary metabolites in Streptomyces species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0787-5) contains supplementary material, which is available to authorized users.

Integrating multi-omics analyses of Nonomuraea dietziae to reveal the role of soybean oil in [(4'-OH)MeLeu]4-CsA overproduction

Background

Nonomuraea dietziae is a promising microorganism to mediate the region-specific monooxygenation reaction of cyclosporine A (CsA). The main product [(4′-OH)MeLeu]⁴-CsA possesses high anti-HIV/HCV and hair growth-stimulating activities while avoiding the immunosuppressive effect of CsA. However, the low conversion efficiency restricts the clinical application. In this study, the production of [(4′-OH)MeLeu]⁴-CsA was greatly improved by 55.6% from 182.8 to 284.4 mg/L when supplementing soybean oil into the production medium, which represented the highest production of [(4′-OH)MeLeu]⁴-CsA so far.

Results

To investigate the effect of soybean oil on CsA conversion, some other plant oils (corn oil and peanut oil) and the major hydrolysates of soybean oil were fed into the production medium, respectively. The results demonstrated that the plant oils, rather than the hydrolysates, could significantly improve the [(4′-OH)MeLeu]⁴-CsA production, suggesting that soybean oil might not play its role in the lipid metabolic pathway. To further unveil the mechanism of [(4′-OH)MeLeu]⁴-CsA overproduction under the soybean oil condition, a proteomic analysis based on the two-dimensional gel electrophoresis coupled with MALDI TOF/TOF mass spectrometry was implemented. The results showed that central carbon metabolism, genetic information processing and energy metabolism were significantly up-regulated under the soybean oil condition. Moreover, the gas chromatography-mass spectrometry-based metabolomic analysis indicated that soybean oil had a great effect on amino acid metabolism and tricarboxylic acid cycle. In addition, the transcription levels of cytochrome P450 hydroxylase (CYP) genes for CsA conversion were determined by RT-qPCR and the results showed that most of the CYP genes were up-regulated under the soybean oil condition.

Conclusions

These findings indicate that soybean oil could strengthen the primary metabolism and the CYP system to enhance the mycelium growth and the monooxygenation reaction, respectively, and it will be a guidance for the further metabolic engineering of this strain.

Electronic supplementary material

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

Comparative metabolomics reveals the mechanism of avermectin production enhancement by S-adenosylmethionine

It was found that S-adenosylmethionine (SAM) could effectively improve avermectin titer with 30–60 μg/mL addition to FH medium. To clearly elucidate the mechanism of SAM on intracellular metabolites of Streptomyces avermitilis, a GC–MS-based comparative metabolomics approach was carried out. First, 230 intracellular metabolites were identified and 14 of them remarkably influenced avermectin biosynthesis were discriminative biomarkers between non-SAM groups and SAM-treated groups by principal components analysis (PCA) and partial least squares (PLS). Based on further key metabolic pathway analyses, these biomarkers, such as glucose, oxaloacetic acid, fatty acids (in soybean oil), threonine, valine, and leucine, were identified as potentially beneficial precursors and added in medium. Compared with single-precursor feeding, the combined feeding of the precursors and SAM markedly increased the avermectin titer. The co-feeding approach not only directly verified our hypothesis on the mechanism of SAM by comparative metabolomics, but also provided a novel strategy to increase avermectin production.

Uncertainty budgeting in fold change determination and implications for non-targeted metabolomics studies in model systems

Uncertainty budgeting provides error intervals for fold change values and complements significance testing in non-targeted metabolomics.

Omics-based approaches reveal phospholipids remodeling of Rhizopus oryzae responding to furfural stress for fumaric acid-production from xylose

In order to relieve the toxicity of furfural on Rhizopus oryzae fermentation, the molecular mechanism of R. oryzae responding to furfural stress for fumaric acid-production was investigated by omics-based approaches. In metabolomics analysis, 29 metabolites including amino acid, sugars, polyols and fatty acids showed significant changes for maintaining the basic cell metabolism at the cost of lowering fumaric acid production. To further uncover the survival mechanism, lipidomics was carried out, revealing that phosphatidylcholine, phosphatidylglycerol, phosphatidylinositol and polyunsaturated acyl chains might be closely correlated with R. oryzae's adapting to furfural stress. Based on the above omics analysis, lecithin, inositol and soybean oil were exogenously supplemented separately with an optimized concentration in the presence of furfural, which increased fumaric acid titer from 5.78g/L to 10.03g/L, 10.05g/L and 12.13g/L (increased by 73.5%, 73.8% and 110%, respectively). These findings provide a methodological guidance for hemicellulose-fumaric acid development.

An overview of rapamycin: from discovery to future perspectives

Rapamycin is an immunosuppressive metabolite produced from several actinomycete species. Besides its immunosuppressive activity, rapamycin and its analogs have additional therapeutic potentials, including antifungal, antitumor, neuroprotective/neuroregenerative, and lifespan extension activities. The core structure of rapamycin is derived from (4R,5R)-4,5-dihydrocyclohex-1-ene-carboxylic acid that is extended by polyketide synthase. The resulting linear polyketide chain is cyclized by incorporating pipecolate and further decorated by post-PKS modification enzymes. Herein, we review the discovery and biological activities of rapamycin as well as its mechanism of action, mechanistic target, biosynthesis, and regulation. In addition, we introduce the many efforts directed at enhancing the production of rapamycin and generating diverse analogs and also explore future perspectives in rapamycin research. This review will also emphasize the remarkable pilot studies on the biosynthesis and production improvement of rapamycin by Dr. Demain, one of the world's distinguished scientists in industrial microbiology and biotechnology.