Mixomics analysis of Bacillus subtilis: effect of oxygen availability on riboflavin production
- Junlang Hu†1,
- Pan Lei†1,
- Ali Mohsin1,
- Xiaoyun Liu1,
- Mingzhi Huang1Email author,
- Liang Li1,
- Jianhua Hu2,
- Haifeng Hang1,
- Yingping Zhuang1 and
- Meijin Guo1Email author
© The Author(s) 2017
Received: 24 April 2017
Accepted: 7 September 2017
Published: 12 September 2017
Riboflavin, an intermediate of primary metabolism, is one kind of important food additive with high economic value. The microbial cell factory Bacillus subtilis has already been proven to possess significant importance for the food industry and have become one of the most widely used riboflavin-producing strains. In the practical fermentation processes, a sharp decrease in riboflavin production is encountered along with a decrease in the dissolved oxygen (DO) tension. Influence of this oxygen availability on riboflavin biosynthesis through carbon central metabolic pathways in B. subtilis is unknown so far. Therefore the unveiled effective metabolic pathways were still an unaccomplished task till present research work.
In this paper, the microscopic regulation mechanisms of B. subtilis grown under different dissolved oxygen tensions were studied by integrating 13C metabolic flux analysis, metabolomics and transcriptomics. It was revealed that the glucose metabolic flux through pentose phosphate (PP) pathway was lower as being confirmed by smaller pool sizes of metabolites in PP pathway and lower expression amount of ykgB at transcriptional level. The latter encodes 6-phosphogluconolactonase (6-PGL) under low DO tension. In response to low DO tension in broth, the glucose metabolic flux through Embden–Meyerhof–Parnas (EMP) pathway was higher and the gene, alsS, encoding for acetolactate synthase was significantly activated that may result due to lower ATP concentration and higher NADH/NAD+ ratio. Moreover, ResE, a membrane-anchored protein that is capable of oxygen regulated phosphorylase activity, and ResD, a regulatory protein that can be phosphorylated and dephosphorylated by ResE, were considered as DO tension sensor and transcriptional regulator respectively.
This study shows that integration of transcriptomics, 13C metabolic flux analysis and metabolomics analysis provides a comprehensive understanding of biosynthesized riboflavin’s regulatory mechanisms in B. subtilis grown under different dissolved oxygen tension conditions. The two-component system, ResD–ResE, was considered as the signal receiver of DO tension and gene regulator that led to differences between biomass and riboflavin production after triggering the shifts in gene expression, metabolic flux distributions and metabolite pool sizes.
KeywordsRiboflavin Bacillus subtilis Oxygen availability 13C metabolic flux analysis Metabolomics Transcriptomics
Riboflavin (vitamin B2) is a water-soluble vitamin that is produced by microorganisms and plants. It is a precursor of flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN), which are cofactors involved in several important metabolic reactions. Nutritional deficiency of riboflavin can cause a variety of fundamental and clinical abnormalities from mitochondrial dysfunction and hemolytic anemia to growth retardation and neurological disorders . Low blood counts, dizziness, hair loss and poor digestion can also be attributed as symptoms to the riboflavin deficiency. It can also be used in the treatment of neonatal jaundice, corneal disorder  and migraine . Beside these nutritive and medicinal values, riboflavin as a supplement is essentially required in animal raising .
As a model organism, Bacillus subtilis has played an important role in the medicinal industry  and possess the potential to produce riboflavin. In order to have commercial scale production of riboflavin by B. subtilis, many researches have focused on the establishment of productive strains [5, 6] and optimal fermentation strategies [7, 8]. In industrial scale fermentation processes, we found that the yield of riboflavin decreased sharply from 1.08 to 0.59 g riboflavin/g glucose while the dissolved oxygen (DO) tension was below 30% of saturated concentration in water (37 °C, 0.15 MPa). However, to maintain DO tension above 30% during the later phase of fermentation is difficult and highly energy-consumptive especially for production level of industrial scale. Hence, enhanced production at lower DO tension can economize the overall process. As we already know, oxygen is an essential substrate in most of aerobic organisms’ metabolic processes because of its significance as electron receptor. Oxygen deficiency may lead to low efficiency of NADH and NADPH oxidation that are the essential reactions to produce ATP and low flux through tricarboxylic acid (TCA) cycle . It is reasonable to hypothesize that a low energy conversion rate can affect the cell multiplication that leads to lower riboflavin production. In past decade, the development of OMICS and high-throughput technology enabled the further study on cellular molecular regulatory mechanism during different biological processes. For examples, the whole genome sequence of B. subtilis 168 was completed in 1997 , since then gene re-sequencing has been used to find out the distinction of different B. subtilis strains based on the genomic DNA data. Global gene expression and metabolomics analysis were used to identify the response to different oxygen availability in the riboflavin fermentation. Furthermore, 13C metabolic flux analysis (MFA) was also employed to reveal the metabolic flux shift between these two conditions. By integrating these omics data and analysis of 13C MFA, it is possible to get a comprehensive understanding of regulatory mechanisms of riboflavin-producing strain from gene expression and metabolic levels.
In the previous studies, many researchers developed the fermentation process of riboflavin by B. subtilis, but the majority of them focused on the macroscopic process behavior  or gene manipulation for strain development using metabolic engineering technology . However, it was the first time that Sauer et al. on quantifying the metabolic fluxes of a wild-type strain and an industrial riboflavin-producing B. subtilis, found that the metabolic flux through PP pathway was growth rate-dependent and higher in the riboflavin-producing strain . Then Shi et al. integrated a comparative transcriptome profiling and metabolic pool analysis to explore the differences between a riboflavin-producing and the wild-type strain, and revealed that the lower expression of PurR-regulated genes was the bottleneck for increasing riboflavin production in their riboflavin-producing strain . These Pur-regulated genes are one set that are involved in purine metabolism and glycine biosynthesis. The expression of these genes can be inhibited by PurR and may lead to lower riboflavin production. To the best of our knowledge, no one ever has revealed the mechanism of riboflavin’s biosynthesis by B. subtilis under the influence of oxygen availability. Fast-growing aerobic microorganisms like B. subtilis require large amount of oxygen for oxidation of NAD(P)H or FADH2 to generate enough ATP for the growth and targeted bio-product formation . Despite its occurrence during fermentation, 12 mol ATP per mol of synthesized product is required that makes this process as highly oxygen consuming. Therefore, this work aims to bring insights into the possible regulatory roles of oxygen in riboflavin biosynthesis by an industrial strain B. subtilis cultivated under different DO tension conditions using integrated 13C MFA, metabolomics and transcriptomics analysis.
Strain and cultivations
Bacillus subtilis was stored at Aceibright Company (Shanghai, China), which was mutated by chemicals like NTG and physical method like UV and 60Co radiation, and genomic DNA was shuffled for many rounds .
Fermentation was carried out as follows: a 100 μL of seed vial, stored in 50% glycerol at −80 °C, it was spread on a Petri dish containing solid medium, and cultured at 37 °C for 24 h. Then, a single colony grown on Petri dish was picked into a 500 mL shaking flask containing 50 mL seed medium, and incubated in a rotary shaker at 37 °C and 220 rpm for 24 h.
Based on the fermentation process, sudden decrease of riboflavin production was observed as DO concentration decreased below 20%. So we considered 20% of saturated DO concentration under 37 °C, 0.15 MPa as critical oxygen concentration and set high DO level above 30% and low DO level in range of 0–10%. Batch cultures in a 5 L stirred bioreactor (FUS-5L, Shanghai, China) were conducted with 2.0% of inoculum for initial 3.0 L working volume at 37 °C. Aeration is kept at 2.0 VVM and process pH was controlled at 7.2 with 2.0 M sodium hydroxide. During the fermentation, agitation speed (maximum 800 rpm) was adopted to maintain dissolved oxygen (DO) tension above 30% of air saturation in the first 24 h, then shift to the high DO level of 35.0 ± 5.0% (Batch H) or low DO level of below 10% (Batch L) till the end of fermentation, respectively. The cultivations for each procedure (high DO level or low DO level) were performed in three independent cultivation replicates.
The solid medium contained (g/L) peptone 15, yeast extract 7.5, NaCl 6.5, maltose 25, agar 20, erythromycin 0.01 and chloramphenicol 0.01. The seed medium contained (g/L) glucose·H2O 10, corn steep liquor 15, molasses 15, yeast extra 7.5, MgSO4·2H2O 0.5, (NH4)2SO4 7.5. The fermentation medium contained (g/L) glucose·H2O 20, K2HPO4·3H2O 15, KH2PO4 6, (NH4)2SO4 8, MgSO4·2H2O 0.2, erythromycin 0.01 and chloramphenicol 0.01. Besides, 10 mL/L trace element buffer that contained (g/L) CaCl2·2H2O 0.55, FeCl3 1.0, MnCl2·4H2O 0.1, ZnCl2 0.17, CuCl2·2H2O 0.043, CoCl2·6H2O 0.06, Na2MoO4 0.06 was added into the fermentation medium.
Quantification of biomass, riboflavin and extracellular glucose
For biomass weight, 5 mL culture was taken and filtered through pre-dried and pre-weighed filter paper. To remove solutes, the cells were rinsed three times with deionized water. Then the wet filters with the biomass were put in the 80 °C oven and dried for 24 h. The dried filter was re-weighed immediately.
Regarding riboflavin determination, 0.05% sodium hydroxide solution was used to dissolve the riboflavin in the broth. After appropriate dilution, absorbance of sample was detected by spectrophotometer (A444).
A glucose analyzing Kit (Shanghai Kexin Biotechnology Institute, China) was used for quantitative analysis for glucose concentration. After centrifugation (4000 rpm, 10 min) and appropriate dilution of samples, they were mixed with reagents provided by the Kit and maintained 37 °C for 20 min before determination by spectrophotometer (A510).
Analysis of extracellular organic acids
The concentrations of extracellular organic acids in the culture were analyzed using a high-performance liquid chromatography (HPLC) system. An ion-exclusion column (Hi-Plex H, Agilent) was eluted at 50 °C with 10 mM H2SO4 at a flow rate of 0.5 mL/min and connected with an absorbance detector spectrophotometer at 210 nm. 1 mL broth was centrifuged at 12,000 rpm, 4 °C for 5 min and the supernatant can be analyzed after appropriate dilution and filtration by filter membrane of diameter 0.22 μm.
Sampling and quantification of intracellular metabolites concentration
Because of the massive leaking of metabolites into the quenching reagent in our previous experiments and other studies , we adapted a rapid filtration protocol for sampling instead of quenching. 1 mL broth was filtered through a vacuum filtration device, then the filter cake was washed by precooled 2.3% (m/m) sodium chloride solution. This step was performed as soon as possible. Subsequently, the filter cake with filter paper was transferred rapidly to 10 mL 75% (v/v) ethanol solution and the extraction continued for 4 min at 95 °C. To accurately determine metabolite’s concentration, 30 μL of 13C-labeled cell extract was added to every sample as internal standard for isotope dilution mass spectrometry (IDMS) before extraction. The supernatant was evaporated through a rotary evaporator at −101 kPa (relative pressure) and then stored at −80 °C before determination by UPLC-MS/MS (Thermo Fisher Scientific Corporation, Waltham, USA).
13C-glucose labeling experiments
13C-labeling experiments were conducted in a 250 mL bioreactor (NCBio, Shanghai, China) with the working volume of 150 mL using glucose as sole carbon source. The labeled glucose contained 30% natural glucose, 20% [1-13C]-glucose and 50% [U-13C]-glucose. The labeling experiments were performed in two independent cultivation bioreactors.
Sampling and extraction of intracellular free amino acids
5 mL broth was taken during the 13C-labeling experiments for analysis of mass isotopomer distributions (MID) of intracellular free amino acids and rapidly mixed with 15 mL pre-cooled quenching reagent (60% methanol, −40 °C). The cells were washed twice with deionized water after 5 min centrifugation at 4 °C, 7500 rpm. Then the intracellular free amino acids were extracted with 75% (v/v) ethanol solution at 95 °C for 4 min. The supernatant was evaporated through a rotary evaporator at −101 kPa (relative pressure) before determination by gas chromatography–mass spectrometry (GC–MS).
13C-labeled metabolic flux analysis
A stoichiometric model of B. subtilis with compartments and reversible reactions was constructed based on the models as previously described . To determine the goodness-of-fit, the 13C-metabolic flux analysis (MFA) fitting results were subjected to a χ2-statistical test. The estimated fluxes were considered acceptable only when the obtained minimal weighted residual was below the χ2 at a 95% confidence level. At convergence, accurate 68 and 95% confidence intervals were computed for all estimated fluxes by evaluating the sensitivity of the minimized SSR to flux variations.
To study the gene expressing distinctions between the high DO level mode and the low DO level mode, during whole fermentation process, three sampling points were set in each culture mode. The supernatant was removed after centrifugation at 4 °C, 4000 rpm for 5 min. Then the cell pellets were stored at −80 °C till analysis. The extraction of mRNA, reverse transcription and sequencing were done by Shanghai Majorbio Bio-Pharm Technology Co., Ltd (Shanghai, China).
Primer sequences used for real-time PCR
Fermentation performance of B. subtilis cultivated under different DO tension conditions
Growth kinetic parameters of B. subtilis in two modes
Rib (36 h) (g/L)
13C metabolic flux analysis
Carbon recovery during 13C-labeling experiments under different dissolved oxygen concentrations
Carbon recovery (%)
Obvious differences were found in the central carbon fluxes between two batches. The relative flux (fluxes normalized to glucose uptake rate) through PP pathway was higher in Batch H (64.85%) than that in Batch L (50.79%). PP pathway offered precursors for nucleotide and amino acid biosynthesis which are essential for the cell growth . Furthermore, Ru-5P is a precursor of riboflavin, also generated by PP pathway. This remarkable flux reduction through PP pathway may be responsible for the lower specific growth rate and riboflavin production rate (Table 2) in Batch L. Compared to Batch H, the fluxes through EMP pathway and TCA cycle were much higher in Batch L. The increased flux through EMP pathway (48.67% vs 33.85%) was predictable because cells would need ATP generated by substrate phosphorylation through EMP pathway under low oxygen availability environment. However, the higher flux through TCA cycle (12.08% vs 7.50%) was not expected. To explain this phenomenon, other analyses were also needed.
It can be speculated based on MFA that in Batch L, lower flux through PP pathway led to lower specific growth rate and riboflavin production rate. Higher flux through EMP pathway and TCA cycle may cause higher NADH generating rate but lower ATP producing rate because of deficiency of electron receptors like oxygen. These results can be explained by following metabolomics and transcriptomics analysis.
Metabolic pools in central carbon metabolism
When it came to TCA cycle, the metabolite pool sizes were larger in the Batch H, except for succinate and fumarate (Fig. 6). This phenomenon was agreed with the trends of organic acids those were secreted into culture medium (Fig. 3). The concentration of 6-PG, the metabolite connecting PP and EMP pathway, was found significantly larger in Batch H and that may be the consequence of larger flux through PP pathway in Batch H.
Gene expression profiling
Differential expression analysis of genes
To study the regulatory mechanism of B. subtilis grown under low dissolved oxygen tension stress, the differential expression analysis of genes between the same phases under two conditions was performed. It was found that there were 224 and 169 differential expression between two conditions in the second and third phase respectively; and 53 genes were overlapped. Interactive pathways explorer 2 (http://pathways.embl.de/iPath2.cgi#) was used to show the results of the analysis (Figures shown in Additional file 2). It is obvious that the reaction from pyruvate to acetoin (alsS) was up-regulated under low oxygen availability. This can be explained that this reaction was proved to be propitious to stabilize the ratio of NAD+/NADH and maintain the redox balance . In addition, the genes that related to the reactions that generate cytochrome bd, from maltose to oxaloacetate and from polysaccharide like starch and sucrose to glucose were also up-regulated. Besides, ResD–ResE, one two-component regulatory factor which has been proved to be responsible for the gene regulation to adapt oxygen tension shifting in B. subtilis [19, 20], was up-regulated significantly in Batch L based on transcriptomic analysis. It was reasonable to suppose that this two-component regulatory factor received the signal of low DO tension and regulated the gene expression globally in cells tested here.
From transcriptomics analysis, global gene expression profiling was performed and the genes involved in purine metabolism were picked out as potential manipulating targets to optimize riboflavin fermentation. Activated ResD–ResE system was supposed to be the sensor of DO tension and regulator of relative genes for the adaptation of DO shift. Gene alsS and ykgB were thought to be the key genes that redistribute the fluxes through central carbon metabolism.
Central carbon metabolism
Based on the result of 13C MFA, the relative flux through PP pathway was significantly larger in Batch H than Batch L, while the situation was vice versa in EMP pathway and TCA cycle (Fig. 5). Various studies proved that, gene zwf is the key gene which can regulate the flux through PP pathway [17, 21]. The manners that regulate the flux shift between PP and EMP pathway are different among organisms and aren’t controlled by the demand of NADPH or pentose [17, 22]. Various signal transduction pathway, in response to different factors in organisms, may affect different genes. It is speculated that when DO concentration decreases, a number of regulatory factors are activated by ResE–ResD leading to a repressed expression of ykgB. Beside the latter phenomenon, cellular physiological status can also affect the gene expression, i.e. higher NADPH/NADP+ ratio repressed expression of zwf in Batch H. In our study, the gene ykgB which encodes 6-phosphogluconolactonase catalyzing the formation of 6-phosphogluconate from 6-phosphogluconolactone was expressed higher in Batch H than Batch L and this may be responsible for the increased flux through PP pathway in Batch H. Obvious difference of expression in zwf between two conditions was not observed in this case. The transcription of zwf in Batch H might be suppressed by higher NADPH/NADP+ ratio. Because of the larger flux through PP pathway, it is not surprising to find out the pool sizes of main metabolites in PP pathway were larger in Batch H. PP pathway is the main source of ribulose 5-phosphate which is one of the precursor of riboflavin, and links with purine metabolism responsible for producing GTP, one another precursor of riboflavin, by PRPP. Besides, PP pathway is a major source of reducing power and metabolic intermediates for cell growth . To summarize, the rise in flux through PP pathway would be attributed as the main reason for higher biomass and riboflavin production.
Higher flux through EMP pathway was observed in Batch L compared with Batch H. It was hypothesized that due to lack of oxygen, NADH/NAD+ ratio (Fig. 7) was increased which may suppress the conversion from pyruvate to acetyl-CoA . Hence, pyruvate was difficult to be consumed and subsequently led to accumulation of EMP downstream metabolites (Fig. 5). The genes that encode acetolactate synthase were significantly up-regulated both at the second and third time point in Batch L compared with Batch H. This reaction converts pyruvate to acetoin which was proved to play an important role in maintaining the NAD+/NADH balance . The vast amount of pyruvate (both extracellular and intracellular) and absence of electron acceptor in Batch L may be the main reason for this stronger expression.
It was surprised to find out that the flux through TCA cycle was higher in Batch L than Batch H, because TCA cycle is a redox process that relies on oxygen. The only reasonable interpretation was the large accumulation of pyruvate that links the EMP pathway and TCA cycle. The high driving force enlarged the flux through TCA cycle, but the NADH produced by this process was accumulated instead of oxidizing to NAD+ by oxygen. Even so, the accumulated NADH can be converted into NAD+ by substrate phosphorylation and that eventually supported the higher flux through EMP pathway in Batch L.
Purine metabolism is important for riboflavin production, because it provides GTP for formation of riboflavin. Comparing the expression of genes that were involved in purine metabolism at all three time points in between Batch H and Batch L; no significant difference was found between the same time points of two batches. However, obvious expression changing trend of genes in purine metabolism (Fig. 10) was found while clustering the whole genes. In two batches, one similar trend was showed that most of the pur genes expression decreased along with the bioprocess both in Batch H and Batch L, although with different degree. This trend indicated that the purine metabolism may be affected by glucose concentration and oxygen availability as those genes were down-regulated more significantly under low DO condition (Fig. 10). And this may be the potential target manipulation for explaining the productivity loss of riboflavin during later phase.
How to compensate lack of electron acceptor
Current study focuses on momentous effect of DO concentration upon riboflavin production. Carbon flux distributions, intracellular metabolite pool size and gene expression profile aided in obtaining major conclusions for this effect. Furthermore, covering of different supporting omics for the conducted study validates the consideration of above stated parameters. However, there lie certain future prospects that can further unveil extent of such promising DO concentration effect for enhancing biosynthesis. Specifically, if pur-related genes, ykgB overexpression or alsS knocking can improve the riboflavin production under low DO condition? If expressions of ykgB or alsS could be regulated by FNR regulator. If providing more nitrate can reduce the NADH/NAD+ ratio, improve ATP generation and riboflavin production? If co-feeding strategy on using different carbon source as substrates could be applied to alleviate the formation of by-products? . All these hypotheses on exploration can open new avenues for effective resource utilization in biosynthesis.
In this paper, integration of 13C metabolic flux analysis, metabolomics and transcriptomics unveils the responses of B. subtilis towards low oxygen tension. ResDE system was verified to be the signal transduction between DO concentration and gene expression by transcriptomics. In addition, the shift in flux from PP to EMP pathway is considered as the main reason of declined riboflavin production under low oxygen availability due to less precursors and reducing power supply as per metabolomics and transcriptomics. The genes involved in purine metabolism would be regarded as genetically possible manipulated targets to maintain high riboflavin production in the later phase of fermentation.
JH performed transcriptomics and metabolomics analysis. PL carried out measurement of extracellular metabolites concentration and 13C metabolic flux analysis. MG and MH participated in design of whole study and analysis of results. AM helped in drafting the manuscript. JH, HH and LL participated in fermentation control. XL and YZ participated in the data discussion. All authors read and approved the final manuscript.
JH, PL, AM, XL, MH, LL, HH, YZ and MG declare no financial interests. JH is employee of Shanghai Acebright Pharmaceuticals Group Co., Ltd.
Availability of data and materials
All data generated or analyzed during this study are included in this article.
Consent for publication
The present manuscript has been seen and proved by all authors. They agree with its content and submission to Microbial Cell Factories.
Ethics approval and consent to participate
The authors declare that this is not a study involving human participants and reporting health related outcomes.
This study was funded by National High Technology Research and Development Program of China (863 Program, 2014AA021705).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pinto J, Rivlin R. Riboflavin (vitamin B2). In: Zempleni J, Suttie JW, Gregory JF III, Stover PJ, editors. Handbook of vitamins. Boca Raton: CRC Press; 2013. p. 191–266.Google Scholar
- Makdoumi K, Mortensen J, Sorkhabi O, Malmvall BE, Crafoord S. UVA-riboflavin photochemical therapy of bacterial keratitis: a pilot study. Graefes Arch Clin Exp Ophthalmol. 2012;250:95–102.View ArticleGoogle Scholar
- Tang J, Wen ZG, Guo ZB, Huang W, Guo YM, Xie M, et al. Dietary riboflavin supplementation improve the growth performance and antioxidant status of starter white Pekin ducks fed a corn–soybean meal diets. Livest Sci. 2014;170:131–6.View ArticleGoogle Scholar
- Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, Azevedo V, et al. The complete genome sequence of the Gram-positive bacterium Bacillus subtilis. Nature. 1997;390:249–56.View ArticleGoogle Scholar
- Sauer U, Hatzimanikatis V, Hohmann H, Manneberg M, Loon A, Bailey J. Physiology and metabolic fluxes of wild-type and riboflavin-producing Bacillus subtilis. Appl Environ Microbiol. 1996;62:3687–96.Google Scholar
- Zhu Y, Chen X, Chen T, Zhao X. Enhancement of riboflavin production by overexpression of acetolactate synthase in a pta mutant of Bacillus subtilis. FEMS Microbiol Lett. 2007;266:224–30.View ArticleGoogle Scholar
- Wu QL, Chen T, Gan Y, Chen X, Zhao XM. Optimization of riboflavin production by recombinant Bacillus subtilis RH44 using statistical designs. Appl Microbiol Biotechnol. 2007;76:783–94.View ArticleGoogle Scholar
- Abd-Alla MH, Bagy MMK, Nafady NA, Morsy FM, Mahmoud AE. Activation of riboflavin production by Bacillus subtilis (KU559874) and Bacillus tequilensis (KU559876). EC Bacteriol Virol Res. 2016;2:131–50.Google Scholar
- Zhang Q, Xiu Z. Metabolic pathway analysis of glycerol metabolism in Klebsiella pneumoniae incorporating oxygen regulatory system. Biotechnol Prog. 2009;25:103–15.View ArticleGoogle Scholar
- Horiuchi J, Hiraga K. Industrial application of fuzzy control to large-scale recombinant vitamin B2 production. J Biosci Bioeng. 1999;87:365–71.View ArticleGoogle Scholar
- Wang Z, Chen T, Ma X, Shen Z, Zhao X. Enhancement of riboflavin production with Bacillus subtilis by expression and site-directed mutagenesis of zwf and gnd gene from Corynebacterium glutamicum. Bioresour Technol. 2011;102:3934–40.View ArticleGoogle Scholar
- Shi S, Chen T, Zhang Z, Chen X, Zhao X. Transcriptome analysis guided metabolic engineering of Bacillus subtilis for riboflavin production. Metab Eng. 2009;11:243–52.View ArticleGoogle Scholar
- Xu H, Dou W, Xu H, Zhang X, Rao Z, Shi Z, et al. A two-stage oxygen supply strategy for enhanced l-arginine production by Corynebacterium crenatum based on metabolic fluxes analysis. Biochem Eng J. 2009;43:41–51.View ArticleGoogle Scholar
- Chen T, Wang J, Zhou S, Chen X, Ban R, Zhao X. Trait improvement of riboflavin-producing Bacillus subtilis by genome shuffling and metabolic flux analysis. J Chem Ind Eng. 2004;55:1842–8.Google Scholar
- Wittmann C, Kromer JO, Kiefer P, Binz T, Heinzle E. Impact of the cold shock phenomenon on quantification of intracellular metabolites in bacteria. Anal Biochem. 2004;327:135–9.View ArticleGoogle Scholar
- Dauner M, Sauer U. Stoichiometric growth model for riboflavin-producing Bacillus subtilis. Biotechnol Bioeng. 2001;76:132–43.View ArticleGoogle Scholar
- Stincone A, Prigione A, Cramer T, Wamelink MM, Campbell K, Cheung E, et al. The return of metabolism: biochemistry and physiology of the pentose phosphate pathway. Biol Rev Camb Philos Soc. 2015;90:927–63.View ArticleGoogle Scholar
- Johansen L, Bryn K, Stormer FC. Physiological and biochemical role of the butanediol pathway in Aerobacter (Enterobacter) aerogenes. J Bacteriol. 1975;123:1124–30.Google Scholar
- Vetter SM, Schlievert PM. The two-component system Bacillus respiratory response A and B (BrrA-BrrB) is a virulence factor regulator in Bacillus anthracis. Biochemistry. 2007;46:7343–52.View ArticleGoogle Scholar
- Nakano MM, Zhu Y. Involvement of ResE phosphatase activity in down-regulation of ResD-controlled genes in Bacillus subtilis during aerobic growth. J Bacteriol. 2001;183:1938–44.View ArticleGoogle Scholar
- Stanton RC. Glucose-6-phosphate dehydrogenase, NADPH, and cell survival. IUBMB Life. 2012;64:362–9.View ArticleGoogle Scholar
- Dauner M, Sonderegger M, Hochuli M, Szyperski T, Wuthrich K, Hohmann HP, et al. Intracellular carbon fluxes in riboflavin-producing Bacillus subtilis during growth on two-carbon substrate mixtures. Appl Environ Microbiol. 2002;68:1760–71.View ArticleGoogle Scholar
- Kruger NJ, von Schaewen A. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6:236–46.View ArticleGoogle Scholar
- Ojima Y, Suryadarma P, Tsuchida K, Taya M. Accumulation of pyruvate by changing the redox status in Escherichia coli. Biotechnol Lett. 2012;34:889–93.View ArticleGoogle Scholar
- And MMN, Zuber P. Anaerobic growth of a “strict aerobe” (Bacillus subtilis). Annu Rev Microbiol. 1998;52:165–90.View ArticleGoogle Scholar
- Yarwood JM, McCormick JK, Schlievert PM. Identification of a novel two-component regulatory system that acts in global regulation of virulence factors of Staphylococcus aureus. J Bacteriol. 2001;183:1113–23.View ArticleGoogle Scholar
- Zhang X, Hulett FM. ResD signal transduction regulator of aerobic respiration in Bacillus subtilis: ctaA promoter regulation. Mol Microbiol. 2000;37:1208–19.View ArticleGoogle Scholar
- Shalel-Levanon S, San KY, Bennett GN. Effect of oxygen, and ArcA and FNR regulators on the expression of genes related to the electron transfer chain and the TCA cycle in Escherichia coli. Metab Eng. 2005;7:364–74.View ArticleGoogle Scholar
- Shalel-Levanon S, San KY, Bennett GN. Effect of ArcA and FNR on the expression of genes related to the oxygen regulation and the glycolysis pathway in Escherichia coli under microaerobic growth conditions. Biotechnol Bioeng. 2005;92:147–59.View ArticleGoogle Scholar