Transcriptome analysis guided metabolic engineering of Bacillus subtilis for riboflavin production

Bacillus subtilis and Escherichia coli essential genes and minimal cell factories after one decade of genome engineering

Investigation of essential genes, besides contributing to understanding the fundamental principles of life, has numerous practical applications. Essential genes can be exploited as building blocks of a tightly controlled cell ‘chassis’. Bacillus subtilis and Escherichia coli K-12 are both well-characterized model bacteria used as hosts for a plethora of biotechnological applications. Determination of the essential genes that constitute the B. subtilis and E. coli minimal genomes is therefore of the highest importance. Recent advances have led to the modification of the original B. subtilis and E. coli essential gene sets identified 10 years ago. Furthermore, significant progress has been made in the area of genome minimization of both model bacteria. This review provides an update, with particular emphasis on the current essential gene sets and their comparison with the original gene sets identified 10 years ago. Special attention is focused on the genome reduction analyses in B. subtilis and E. coli and the construction of minimal cell factories for industrial applications.

Metabolic engineering of Escherichia coli for the production of riboflavin

BACKGROUND

Riboflavin (vitamin B2), the precursor of the flavin cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD), is used commercially as an animal feed supplement and food colorant. E. coli is a robust host for various genetic manipulations and has been employed for efficient production of biofuels, polymers, amino acids, and bulk chemicals. Thus, the aim of this study was to understand the metabolic capacity of E. coli for the riboflavin production by modification of central metabolism, riboflavin biosynthesis pathway and optimization of the fermentation conditions.

RESULTS

The basic producer RF01S, in which the riboflavin biosynthesis genes ribABDEC from E. coli were overexpressed under the control of the inducible trc promoter, could accumulate 229.1 mg/L of riboflavin. Further engineering was performed by examining the impact of expression of zwf (encodes glucose 6-phosphate dehydrogenase) and gnd (encodes 6-phosphogluconate dehydrogenase) from Corynebacterium glutamicum and pgl (encodes 6-phosphogluconolactonase) from E. coli on riboflavin production. Deleting pgi (encodes glucose-6-phosphate isomerase) and genes of Entner-Doudoroff (ED) pathway successfully redirected the carbon flux into the oxidative pentose phosphate pathway, and overexpressing the acs (encodes acetyl-CoA synthetase) reduced the acetate accumulation. These modifications increased riboflavin production to 585.2 mg/L. By further modulating the expression of ribF (encodes riboflavin kinase) for reducing the conversion of riboflavin to FMN in RF05S, the final engineering strain RF05S-M40 could produce 1036.1 mg/L riboflavin in LB medium at 37°C. After optimizing the fermentation conditions, strain RF05S-M40 produced 2702.8 mg/L riboflavin in the optimized semi-defined medium, which was a value nearly 12-fold higher than that of RF01S, with a yield of 137.5 mg riboflavin/g glucose.

CONCLUSIONS

The engineered strain RF05S-M40 has the highest yield among all reported riboflavin production strains in shake flask culture. This work collectively demonstrates that E. coli has a potential to be a microbial cell factory for riboflavin bioproduction.

Deregulation of purine pathway in Bacillus subtilis and its use in riboflavin biosynthesis

Background

Purine nucleotides are essential metabolites for living organisms because they are involved in many important processes, such as nucleic acid synthesis, energy supply, and biosynthesis of several amino acids and riboflavin. Owing to the pivotal roles of purines in cell physiology, the pool of intracellular purine nucleotides must be maintained under strict control, and hence the de novo purine biosynthetic pathway is tightly regulated by transcription repression and inhibition mechanism. Deregulation of purine pathway is essential for this pathway engineering in Bacillus subtilis.

Results

Deregulation of purine pathway was attempted to improve purine nucleotides supply, based on a riboflavin producer B. subtilis strain with modification of its rib operon. To eliminate transcription repression, the pur operon repressor PurR and the 5’-UTR of pur operon containing a guanine-sensing riboswitch were disrupted. Quantitative RT-PCR analysis revealed that the relative transcription levels of purine genes were up-regulated about 380 times. Furthermore, site-directed mutagenesis was successfully introduced into PRPP amidotransferase (encoded by purF) to remove feedback inhibition by homologous alignment and analysis. Overexpression of the novel mutant PurF (D293V, K316Q and S400W) significantly increased PRPP amidotransferase activity and triggered a strong refractory effect on purine nucleotides mediated inhibition. Intracellular metabolite target analysis indicated that the purine nucleotides supply in engineered strains was facilitated by a stepwise gene-targeted deregulation. With these genetic manipulations, we managed to enhance the metabolic flow through purine pathway and consequently increased riboflavin production 3-fold (826.52 mg/L) in the purF-VQW mutant strain.

Conclusions

A sequential optimization strategy was applied to deregulate the rib operon and purine pathway of B. subtilis to create genetic diversities and to improve riboflavin production. Based on the deregulation of purine pathway at transcription and metabolic levels, an extended application is recommended for the yield of products, like inosine, guanosine, adenosine and folate which are directly stemming from purine pathway in B. subtilis.

Improved production of a heterologous amylase in Saccharomyces cerevisiae by inverse metabolic engineering

The increasing demand for industrial enzymes and biopharmaceutical proteins relies on robust production hosts with high protein yield and productivity. Being one of the best-studied model organisms and capable of performing posttranslational modifications, the yeast Saccharomyces cerevisiae is widely used as a cell factory for recombinant protein production. However, many recombinant proteins are produced at only 1% (or less) of the theoretical capacity due to the complexity of the secretory pathway, which has not been fully exploited. In this study, we applied the concept of inverse metabolic engineering to identify novel targets for improving protein secretion. Screening that combined UV-random mutagenesis and selection for growth on starch was performed to find mutant strains producing heterologous amylase 5-fold above the level produced by the reference strain. Genomic mutations that could be associated with higher amylase secretion were identified through whole-genome sequencing. Several single-point mutations, including an S196I point mutation in the VTA1 gene coding for a protein involved in vacuolar sorting, were evaluated by introducing these to the starting strain. By applying this modification alone, the amylase secretion could be improved by 35%. As a complement to the identification of genomic variants, transcriptome analysis was also performed in order to understand on a global level the transcriptional changes associated with the improved amylase production caused by UV mutagenesis.

Overexpression of a non-native deoxyxylulose-dependent vitamin B6 pathway in Bacillus subtilis for the production of pyridoxine

Vitamin B6 is a designation for the vitamers pyridoxine, pyridoxal, pyridoxamine, and their respective 5'-phosphates. Pyridoxal 5'-phosphate, the biologically most-important vitamer, serves as a cofactor for many enzymes, mainly active in amino acid metabolism. While microorganisms and plants are capable of synthesizing vitamin B6, other organisms have to ingest it. The vitamer pyridoxine, which is used as a dietary supplement for animals and humans is commercially produced by chemical processes. The development of potentially more cost-effective and more sustainable fermentation processes for pyridoxine production is of interest for the biotech industry. We describe the generation and characterization of a Bacillus subtilis pyridoxine production strain overexpressing five genes of a non-native deoxyxylulose 5'-phosphate-dependent vitamin B6 pathway. The genes, derived from Escherichia coli and Sinorhizobium meliloti, were assembled to two expression cassettes and introduced into the B. subtilis chromosome. in vivo complementation assays revealed that the enzymes of this pathway were functionally expressed and active. The resulting strain produced 14mg/l pyridoxine in a small-scale production assay. By optimizing the growth conditions and co-feeding of 4-hydroxy-threonine and deoxyxylulose the productivity was increased to 54mg/l. Although relative protein quantification revealed bottlenecks in the heterologous pathway that remain to be eliminated, the final strain provides a promising basis to further enhance the production of pyridoxine using B. subtilis.

Kinetic modeling of riboflavin biosynthesis in Bacillus subtilis under production conditions

To study the network dynamics of the riboflavin biosynthesis pathway and to identify potential bottlenecks in the system, an ordinary differential equation-based model was constructed using available literature data for production strains. The results confirmed that the RibA protein is rate limiting in the pathway. Under the conditions investigated, we determined a potential limiting order of the remaining enzymes under increased RibA concentration (>0.102 mM) and therefore higher riboflavin production (>0.045 mmol g(CDW)(-1) h(-1) and 0.0035 mM s(-1), respectively). The reductase activity of RibG and lumazine synthase (RibH) might be the next most limiting steps. The computational minimization of the enzyme concentrations of the pathway suggested the need for a greater RibH concentration (0.251 mM) compared with the other enzymes (RibG: 0.188 mM, RibB: 0.023 mM).

Construction and fed-batch cultivation of Candida famata with enhanced riboflavin production

Riboflavin (vitamin B2) is an essential nutrition component serving as a precursor of coenzymes FMN and FAD that are involved mostly in reactions of oxidative metabolism. Riboflavin is produced in commercial scale and is used in feed and food industries, and in medicine. The yeast Candida famata (Candida flareri) belongs to the group of so called "flavinogenic yeasts" which overproduce riboflavin under iron limitation. Three genes SEF1, RIB1 and RIB7 coding for a putative transcription factor, GTP cyclohydrolase II and riboflavin synthase, respectively were simultaneously overexpressed in the background of a non-reverting riboflavin producing mutant AF-4, obtained earlier in our laboratory using methods of classical selection (Dmytruk et al. (2011), Metabolic Engineering 13, 82-88). Cultivation conditions of the constructed strain were optimized for shake-flasks and bioreactor cultivations. The constructed strain accumulated up to 16.4g/L of riboflavin in optimized medium in a 7L laboratory bioreactor during fed-batch fermentation.

Multiplex iterative plasmid engineering for combinatorial optimization of metabolic pathways and diversification of protein coding sequences

Engineering complex biological systems typically requires combinatorial optimization to achieve the desired functionality. Here, we present Multiplex Iterative Plasmid Engineering (MIPE), which is a highly efficient and customized method for combinatorial diversification of plasmid sequences. MIPE exploits ssDNA mediated λ Red recombineering for the introduction of mutations, allowing it to target several sites simultaneously and generate libraries of up to 10(7) sequences in one reaction. We also describe "restriction digestion mediated co-selection (RD CoS)", which enables MIPE to produce enhanced recombineering efficiencies with greatly simplified coselection procedures. To demonstrate this approach, we applied MIPE to fine-tune gene expression level in the 5-gene riboflavin biosynthetic pathway and successfully isolated a clone with 2.67-fold improved production in less than a week. We further demonstrated the ability of MIPE for highly multiplexed diversification of protein coding sequence by simultaneously targeting 23 codons scattered along the 750 bp sequence. We anticipate this method to benefit the optimization of diverse biological systems in synthetic biology and metabolic engineering.

Comparative transcriptome analysis for metabolic engineering

Transcriptome profiling allows massively parallel analysis of the dynamic expression of all genes and captures the cell physiology and regulatory mechanism in a holistic manner. Compared to other "omic" techniques, transcriptome is more tractable and sensitive. Transcriptomics has profoundly promoted development and applications of metabolic engineering by analysis of cell metabolism at a system level. Our recent effort was performed on a comparative transcriptome profiling between a riboflavin-producing Bacillus subtilis strain RH33 and the wild-type strain B. subtilis 168 to rationally identify new targets for improving riboflavin production. This transcriptome analysis-guided method improved the riboflavin titer by 32 ± 3 %. Herein, we describe the detailed experimental protocols for predicting new engineering targets using comparative transcriptome analysis.

In silico aided metabolic engineering of Streptomyces roseosporus for daptomycin yield improvement

In silico metabolic network models are valuable tools for strain improvement with desired properties. In this work, based on the comparisons of each pathway flux under two different objective functions for the reconstructed metabolic network of Streptomyces roseosporus, three potential targets of zwf2 (code for glucose-6-phosphate hydrogenase), dptI (code for α-ketoglutarate methyltransferase), and dptJ (code for tryptophan oxygenase) were identified and selected for the genetic modifications. Overexpression of zwf2, dptI, and dptJ genes increased the daptomycin concentration up to 473.2, 452.5, and 489.1 mg/L, respectively. Furthermore, co-overexpression of three genes in series resulted in a 34.4% higher daptomycin concentration compared with the parental strain, which ascribed to the synergistic effect of the enzymes responsible for daptomycin biosynthesis. Finally, the engineered strain enhanced the yield of daptomycin up to 581.5 mg/L in the fed-batch culture, which was approximately 43.2% higher than that of the parental strain. These results demonstrated that the metabolic network based on in silico prediction would be accurate, reasonable, and practical for target gene identification and strain improvement.

One-step of tryptophan attenuator inactivation and promoter swapping to improve the production of L-tryptophan in Escherichia coli

BACKGROUND

L-tryptophan is an aromatic amino acid widely used in the food, chemical and pharmaceutical industries. In Escherichia coli, L-tryptophan is synthesized from phosphoenolpyruvate and erythrose 4-phosphate by enzymes in the shikimate pathway and L-tryptophan branch pathway, while L-serine and phosphoribosylpyrophosphate are also involved in L-tryptophan synthesis. In order to construct a microbial strain for efficient L-tryptophan production from glucose, we developed a one step tryptophan attenuator inactivation and promoter swapping strategy for metabolic flux optimization after a base strain was obtained by overexpressing the tktA, mutated trpE and aroG genes and inactivating a series of competitive steps.

RESULTS

The engineered E. coli GPT1002 with tryptophan attenuator inactivation and tryptophan operon promoter substitution exhibited 1.67 ~ 9.29 times higher transcription of tryptophan operon genes than the control GPT1001. In addition, this strain accumulated 1.70 g l(-1) L-tryptophan after 36 h batch cultivation in 300-mL shake flask. Bioreactor fermentation experiments showed that GPT1002 could produce 10.15 g l(-1) L-tryptophan in 48 h.

CONCLUSIONS

The one step inactivating and promoter swapping is an efficient method for metabolic engineering. This method can also be applied in other bacteria.

Wild-type and feedback-resistant phosphoribosyl pyrophosphate synthetases from Bacillus amyloliquefaciens: purification, characterization, and application to increase purine nucleoside production

Bacillus strains are used for the industrial production of the purine nucleosides inosine and guanosine, which are raw materials for the synthesis of the flavor enhancers disodium inosinate and disodium guanylate. An important precursor of purine nucleosides is 5-phospho-α-D: -ribosyl-1-pyrophosphate, which is synthesized by phosphoribosyl pyrophosphate synthetase (PRS, EC 2.7.6.1). Class I PRSs are widespread in bacteria and mammals, are highly conserved among different organisms, and are negatively regulated by two end products of purine biosynthesis, adenosine 5'-diphosphate (ADP) and guanosine 5'-diphosphate (GDP). The D52H, N114S, and L129I mutations in the human PRS isozyme I (PRS1) have been reported to cause uric acid overproduction and gout due to allosteric deregulation and enzyme superactivity. In this study, to find feedback-resistant Bacillus amyloliquefaciens PRS, the influence of the D58H, N120S, and L135I mutations (corresponding to the D52H, N114S, and L129I mutations in PRS1, respectively) on PRS enzymatic properties has been studied. Recombinant histidine-tagged wild-type PRS and three mutant PRSs were expressed in Escherichia coli, purified, and characterized. The N120S and L135I mutations were found to release the enzyme from ADP and GDP inhibition and significantly increase its sensitivity to inorganic phosphate (P(i)) activation. In contrast, PRS with the D58H mutation exhibited nearly identical sensitivity to ADP and GDP as the wild-type protein and had a notably greater P(i) requirement for activation. The N120S and L135I mutations improved B. amyloliquefaciens and Bacillus subtilis purine nucleoside-producing strains.

Synthetic biology: an emerging research field in China

Synthetic biology is considered as an emerging research field that will bring new opportunities to biotechnology. There is an expectation that synthetic biology will not only enhance knowledge in basic science, but will also have great potential for practical applications. Synthetic biology is still in an early developmental stage in China. We provide here a review of current Chinese research activities in synthetic biology and its different subfields, such as research on genetic circuits, minimal genomes, chemical synthetic biology, protocells and DNA synthesis, using literature reviews and personal communications with Chinese researchers. To meet the increasing demand for a sustainable development, research on genetic circuits to harness biomass is the most pursed research within Chinese researchers. The environmental concerns are driven force of research on the genetic circuits for bioremediation. The research on minimal genomes is carried on identifying the smallest number of genomes needed for engineering minimal cell factories and research on chemical synthetic biology is focused on artificial proteins and expanded genetic code. The research on protocells is more in combination with the research on molecular-scale motors. The research on DNA synthesis and its commercialisation are also reviewed. As for the perspective on potential future Chinese R&D activities, it will be discussed based on the research capacity and governmental policy.

Genetic control of biosynthesis and transport of riboflavin and flavin nucleotides and construction of robust biotechnological producers

Riboflavin [7,8-dimethyl-10-(1'-d-ribityl)isoalloxazine, vitamin B₂] is an obligatory component of human and animal diets, as it serves as the precursor of flavin coenzymes, flavin mononucleotide, and flavin adenine dinucleotide, which are involved in oxidative metabolism and other processes. Commercially produced riboflavin is used in agriculture, medicine, and the food industry. Riboflavin synthesis starts from GTP and ribulose-5-phosphate and proceeds through pyrimidine and pteridine intermediates. Flavin nucleotides are synthesized in two consecutive reactions from riboflavin. Some microorganisms and all animal cells are capable of riboflavin uptake, whereas many microorganisms have distinct systems for riboflavin excretion to the medium. Regulation of riboflavin synthesis in bacteria occurs by repression at the transcriptional level by flavin mononucleotide, which binds to nascent noncoding mRNA and blocks further transcription (named the riboswitch). In flavinogenic molds, riboflavin overproduction starts at the stationary phase and is accompanied by derepression of enzymes involved in riboflavin synthesis, sporulation, and mycelial lysis. In flavinogenic yeasts, transcriptional repression of riboflavin synthesis is exerted by iron ions and not by flavins. The putative transcription factor encoded by SEF1 is somehow involved in this regulation. Most commercial riboflavin is currently produced or was produced earlier by microbial synthesis using special selected strains of Bacillus subtilis, Ashbya gossypii, and Candida famata. Whereas earlier RF overproducers were isolated by classical selection, current producers of riboflavin and flavin nucleotides have been developed using modern approaches of metabolic engineering that involve overexpression of structural and regulatory genes of the RF biosynthetic pathway as well as genes involved in the overproduction of the purine precursor of riboflavin, GTP.

One-step production of lactate from cellulose as the sole carbon source without any other organic nutrient by recombinant cellulolytic Bacillus subtilis

Although intensive efforts have been made to create recombinant cellulolytic microorganisms, real recombinant cellulose-utilizing microorganisms that can produce sufficient secretory active cellulase, hydrolyze cellulose, and utilize released soluble sugars for supporting both cell growth and cellulase synthesis without any other organic nutrient (e.g., yeast extract, peptone, amino acids), are not available. Here we demonstrated that over-expression of Bacillus subtilis endoglucanase BsCel5 enabled B. subtilis to grow on solid cellulosic materials as the sole carbon source for the first time. Furthermore, two-round directed evolution was conducted to increase specific activity of BsCel5 on regenerated amorphous cellulose (RAC) and enhance its expression/secretion level in B. subtilis. To increase lactate yield, the alpha-acetolactate synthase gene (alsS) in the 2,3-butanediol pathway was knocked out. In the chemically defined minimal M9/RAC medium, B. subtilis XZ7(pBscel5-MT2C) strain (ΔalsS), which expressed a BsCel5 mutant MT2C, was able to hydrolyze RAC with cellulose digestibility of 74% and produced about 3.1g/L lactate with a yield of 60% of the theoretical maximum. When 0.1% (w/v) yeast extract was added in the M9/RAC medium, cellulose digestibility and lactate yield were enhanced to 92% and 63% of the theoretical maximum, respectively. The recombinant industrially safe cellulolytic B. subtilis would be a promising consolidated bioprocessing platform for low-cost production of biocommodities from cellulosic materials.

Enhancement of riboflavin production with Bacillus subtilis by expression and site-directed mutagenesis of zwf and gnd gene from Corynebacterium glutamicum

Zwf (code for glucose-6-phosphate dehydrogenase) and gnd (code for 6-phosphogluconate dehydrogenase) genes from Corynebacterium glutamicum were firstly cloned, and then site-directed mutagenesis was successfully introduced to remove allosteric inhibition by intracellular metabolites. Expression of the mutant zwf and gnd in Bacillus subtilis RH33 resulted in significant enhancement of riboflavin productivity, while the specific growth rate decreased slightly and the specific glucose uptake rate was unchanged. Introduction of the mutant zwf and gnd led to approximately 18% and 22% increased riboflavin production, respectively. An improvement by 31% and 39% of the riboflavin production was obtained by co-expression of the mutated dehydrogenases in shaker flask and fed-batch cultivation. Intracellular metabolites analysis indicated that metabolites detected in pentose phosphate pathway or riboflavin synthesis pathway of engineered strains showed higher concentration, while TCA cycle and glycolysis metabolites detected were lower abundance than that of parent strain.

Control of α-amylase production by Bacillus subtilis

This study proposes two adaptive control algorithms for the fed-batch production of α-amylase. The first one uses online information from hardware measuring glucose. Online information of both biomass and glucose concentrations measured with different frequency is used in the second algorithm. Hardware measuring variables are inputs for software sensors of glucose concentration and (specific) glucose consumption rate. Either of the algorithms do not require any kinetic coefficients. This is a benefit, because the kinetic coefficients can vary during cultivation and between cultivations, leading to low process reproducibility and the non-stationary state of the bioprocess. The results of simulation investigations show good performance of the proposed control schemes.

Metabolic engineering and classic selection of the yeast Candida famata (Candida flareri) for construction of strains with enhanced riboflavin production

Currently, the mutant of the flavinogenic yeast Candida famata dep8 isolated by classic mutagenesis and selection is used for industrial riboflavin production. Here we report on construction of a riboflavin overproducing strain of C. famata using a combination of random mutagenesis based on the selection of mutants resistant to different antimetabolites as well as rational approaches of metabolic engineering. The conventional mutagenesis involved consecutive selection for resistance to riboflavin structural analog 7-methyl-8-trifluoromethyl-10-(1'-d-ribityl)isoalloxazine), 8-azaguanine, 6-azauracil, 2-diazo-5-oxo-L-norleucine and guanosine as well as screening for yellow colonies at high pH. The metabolic engineering approaches involved introduction of additional copies of transcription factor SEF1 and IMH3 (coding for IMP dehydrogenase) orthologs from Debaryomyces hansenii, and the homologous genes RIB1 and RIB7, encoding GTP cyclohydrolase II and riboflavin synthetase, the first and the last enzymes of riboflavin biosynthesis pathway, respectively. Overexpression of the aforementioned genes in riboflavin overproducer AF-4 obtained by classical selection resulted in a 4.1-fold increase in riboflavin production in shake-flask experiments. D. hansenii IMH3 and modified ARO4 genes conferring resistance to mycophenolic acid and fluorophenylalanine, respectively, were successfully used as new dominant selection markers for C. famata.