- Research article
- Open Access
Reconstruction and analysis of a genome-scale metabolic model of the oleaginous fungus Mortierella alpina
BMC Systems Biologyvolume 9, Article number: 1 (2015)
Mortierella alpina is an oleaginous fungus used in the industrial scale production of arachidonic acid (ARA). In order to investigate the metabolic characteristics at a systems level and to explore potential strategies for enhanced lipid production, a genome-scale metabolic model of M. alpina was reconstructed.
This model included 1106 genes, 1854 reactions and 1732 metabolites. On minimal growth medium, 86 genes were identified as essential, whereas 49 essential genes were identified on yeast extract medium. A series of sequential desaturase and elongase catalysed steps are involved in the synthesis of polyunsaturated fatty acids (PUFAs) from acetyl-CoA precursors, with concomitant NADPH consumption, and these steps were investigated in this study. Oxygen is known to affect the degree of unsaturation of PUFAs, and robustness analysis determined that an oxygen uptake rate of 2.0 mmol gDW−1 h−1 was optimal for ARA accumulation. The flux of 53 reactions involving NADPH was significantly altered at different ARA levels. Of these, malic enzyme (ME) was confirmed as a key component in ARA production and NADPH generation. When using minimization of metabolic adjustment, a knock-out of ME led to a 38.28% decrease in ARA production.
The simulation results confirmed the model as a useful tool for future research on the metabolism of PUFAs.
Mortierella alpina is an oleaginous zygomycete, that can accumulate lipids to 50% of its dry weight in the form of triacylglycerols . The important ω-6 polyunsaturated fatty acid (PUFA) arachidonic acid (ARA) can account for over 50% of the lipid content . M. alpina is nonpathogenic and nonallergenic, including the spores produced during the industrial production of ARA  which is widely used in food ingredients . ARA has been produced at levels up to 19.8 g/L in 5 L cultures grown over 7 days . Various methods have been attempted in order to improve ARA production including screening potentially higher yielding mutant strains following treatment with N-methyl-N’-nitro-N-nitrosoguanidine (MNNG) . This work led to the generation of M. alpina strain Y11 which possessed lowered ω-3 desaturation activity, and ARA production was 2.1 fold (2.21 g/L) higher than strain 1S-4. Optimizing the culture medium and fixing the ratio of defatted soybean meal to potassium nitrate at 2:1 gave a fourfold increase (6.0 g/L) in ARA production . In another study the fermentation process was optimized, and a two-stage temperature-shift strategy increased ARA production by 26.1% (9.2 g/L) . Additionally, overexpressing GLELO gene using a genetic manipulation approach increased ARA production by 101.2% (5.05 g/L) .
Genetic engineering of M. alpina for enhanced ARA production remains an attractive proposition and for further research. Early work identified desaturases as key enzymes in PUFA synthesis. Specifically, Δ5 desaturase , which catalyzes dehydrogenation of dihomo-γ-linolenic acid (DGLA) to form ARA, was isolated and functionally characterized. The rate-limiting step for ARA biosynthesis is catalyzed by elongase which converts γ-linolenic acid (GLA) to DGLA [11,12]. NADH-cytochrome b5 reductase (Cb5R), an electron carrier and a major component of the cytochrome b5-dependent electron transport system, is also crucial. Cb5R catalyzes several different reduction reactions, including the desaturation and elongation of acyl chains built from acetyl-CoA during PUFA synthesis . Despite various studies that have identified the importance of these enzymes in PUFA synthesis and metabolism [2,14], their exact roles are not completely understood, and neither are the pathways through which glucose relates to PUFA biosynthesis and metabolism.
A genome-scale metabolic model (GSMM) is an indispensable tool for the study of metabolism that adopts a systems biology approach to integrate data from genomics, transcriptomics, proteomics and metabolomics. It has been widely used in the analysis of the network properties of metabolism , prediction and analysis of organism growth phenotypes , model-based interpretation of experimental data , and metabolic engineering . Oleaginous organisms such as M. alpina can accumulate large quantities of lipids, but maximizing lipid production is complicated by the complexity of the regulatory mechanisms associated with lipid metabolism. It is generally difficult to identify key metabolic modules contributing to lipid physiology. Using reconstruction GSMM, we can systematically analyze the function of each gene and metabolic reaction and model the effects using flux balance analysis (FBA). Specific pathways can be understood based on the model of the whole metabolic network, and strain design strategies can also be used to guide metabolic engineering experiments. Two GSMM studies on Yarrowia lipolytica (iNL895  and iYL619_PCP ) have been published along with recent modeled networks of Mucor circinelloides and M. alpina . GSMM studies therefore provide a new approach to investigating the complex lipid metabolism in M. alpina. Vongsangnak et al. (2013)  previously published a M. alpina network model, however, this was a refined network that could only be used to investigate genome annotation and metabolic routes, and not flux distribution or phenotypic behaviors . To systematically study flux distribution and the mechanism of lipid accumulation, we reconstructed a new M. alpina GSMM and used the COBRA Toolbox  for subsequent research.
In this study, a genome-scale metabolic model (iCY1106) of M. alpina ATCC 32222 was reconstructed based on sequencing results . Using this model, we first identified essential genes and reactions in fermentation medium containing sodium nitrate as a nitrogen source . The de novo synthesis pathways of PUFAs such as ARA and eicosapentaenoic acid (EPA) were subsequently resolved. The roles of acetyl-CoA and NADPH in the regulation of PUFA biosynthesis were probed, and important genes and reactions were systematically analyzed using FBA, flux variability analysis (FVA), and minimization of metabolic adjustment (MOMA).
Results and discussions
Characterization of the M. alpina GSMM iCY1106
The genome-scale metabolic model iCY1106 reconstructed in this study included 1106 genes representing 9.51% of the 11631 protein-coding genes in the genome of M. alpina ATCC 32222. The model consisted of 1854 reactions and 1732 metabolites (Table 1). Compared with the previous model constructed by Vongsangnak , iCY1106 is more detailed and contains a higher degree of functional information that can be used in in silico simulations. During the construction of iCY1106, some missing gene were added from the Vongsangnak’s model such as MA-090-247, which encodes an asparaginase (EC: 18.104.22.168), and catalyzes the transformation of L-asparagine into L-aspartate. When compared with the existing model, many improvements in iCY1106 are apparent. (1) A total of 139 transport proteins were annotated based on TCDB and TransportDB. (2) A total of 216 exchange reactions were added into model iCY1106. (3) Biomass composition of M. alpina was calculated based on literature mining. (4) Additional chemical reactions were also included. (5) The reversibility of some reactions was redefined according to MetaCyc and BioPath databases. For example, in the published model, the reaction catalyzed by aldehyde dehydrogenase (EC: 22.214.171.124) was reversible, whereas this reaction is irreversible according to BioPath. The flux distribution of iCY1106 was normal following this change.
In model iCY1106, four compartments (extracellular, cytoplasmic, peroxisome and mitochondrial compartments) were linked by 141 trans-plasma membrane transport reactions, 98 cytoplasmic-mitochondrial transport reactions and eight cytoplasmic-peroxisome transport reactions (Additional file 1). The extracellular compartment contained 267 reactions, including 216 exchange reactions and 51 metabolic reactions. These metabolic reactions catalyzed by extracellular enzymes may enable M. alpina to adapt to different environmental conditions. For example, M. alpina grows on different carbohydrates such as glycerol , sucrose , d-mannose , and raffinose  using various extracellular galactosidases. The peroxisome compartment contained 34 reactions mostly related to glyoxylate metabolism.
Metabolic reactions in the model iCY1106 were divided into ten metabolic subsystems (Figure 1). The largest subsystem was lipid metabolism (fatty acid biosynthesis, fatty acid elongation, steroid, glycerolipid, glycerophospholipid and sphingolipid metabolism) that accounted for 23.53% of the total reactions. Reactions involving the synthesis of PUFAs such as DGLA, ARA, and EPA are listed in model iCY1106. For most subsystems, most categories were more abundant in iCY1106 than in the Vongsangnak’s model, with the exception of Amino Acid Metabolism and Nucleotide Metabolism (Figure 1a). Since there were six compartments (cytosol, mitochondria, peroxisome, extracellular, plasma membrane and unclassified compartment) in Vongsangnak’s model, some reactions occur in multiple compartments at the same time. When ignoring compartments, not including transport reactions and exchange reactions, these two models shared 566 reactions (Additional file 2). Vongsangnak’s model and iCY1106 included 179 and 621 unique reactions, respectively (Figure 1b). Compared with model iYL619_PCP of the oleaginous yeast Y. lipolytica , lipid metabolism accounted for a higher proportion (16.7%, 191 of total 1142 reactions) of reactions in iCY1106. Amino acid and carbohydrate metabolism subsystems were 13.17% and 11.66%, respectively, in iCY1106.
Identification of essential genes and reactions in iCY1106
When cultured on minimal growth medium (MG) containing glucose, KH2PO4, MgSO4 and NaNO3 , 86 M. alpina genes (7.78% of the total) were identified as essential for growth using the FBA method (Additional file 3). In contrast, on yeast extract (YE) medium (based on MG but including all 20 regular amino acids) only 49 genes were identified as essential. On MG medium, over half of the essential genes were involved in amino acid (36.05%) and nucleotide metabolism (23.26%; Figure 2). However on YE medium, most of the essential genes were associated with nucleotide metabolism (38.78%).
The identified genes were compared with the database of essential genes (DEG), and 40 genes were found to share high homology with genes identified by DEG (identity ≥ 40%, E-value ≤ 1.0E-30; Table 2). One gene, UTP-glucose-1-phosphate uridylyltransferase (MA-090-485), is involved in the synthesis of the cell wall component 1,3-beta-D-glucan (Table 2). MA-090-485 is highly homologous with the Saccharomyces cerevisiae enzyme (sequence identity = 62.58%, e-value = 0) identified in the DEG search. Similarly, nucleoside diphosphate kinase (MA-120-261), which is involved in the synthesis of CTP and GTP, also shared high homology with one of the essential genes identified by DEG (identity = 66.89%, e-value = 2.0E-68). Besides, in MG, there were 46 genes that did not match with DEG (Additional file 3). Of these, when nitrate was used as nitrogen source, nitrate reductase (MA-291-2) and nitrite reductase (MA-291-3), which together convert nitrate into NH3 for cell growth, were essential under this condition. Another two essential genes, delta 5 desaturase (MA-326-160) and delta 12 desaturase (MA-334-414) are vital for the synthesis of ARA in M. alpina. In addition, most of these essential genes were distributed in the Amino Acid Metabolism category (28.3%). This was understandable since essential gene identified by DEG were in a rich medium in which genes involving amino acid synthesis may not be essential.
Simulation and verification of model iCY1106
FBA with constraints was used to investigate the metabolic properties of model iCY1106. YE medium was used to simulate the growth rate and ARA production, the maximum uptake rate of all 20 amino acids were set at 0.01 mmol gDW−1 h-1 , and the glucose consumption rate was set at 0.8 mmol gDW−1 h-1 , and 2.0 mmol gDW−1 h-1  respectively, and the in silico growth rate were closed to experimental values (differences were only 0.8% and 9.48%; Table 3). To investigate ARA production in silico, the growth rate was fixed at 0.03 h-1 , and the ARA exchange reaction was set as the objective function. Under these conditions, the ARA rate was 0.128 mmol gDW−1 h−1 by FBA. This in silico value was lower than the experimentally determined value of 0.149 mmol gDW−1 h−1 (14.09%) , and this may due to the use of two different M. alpina strains (ATCC 32222 and ME-1) for modeling and experimental components, with ME-1 being the higher ARA producing strain . Furthermore, the difficulties associated with accurately modeling the yeast extract may have resulted in a nutritional deficiency in the model.
The capability of utilizing 14 different carbon sources for cell growth was qualitatively predicted by the model iCY1106 using FBA (Additional file 4). Each substance was the sole carbon source and the rate of uptake was 1.0 mmol gDW−1 h−1. Initially, there were four types of carbon source (ethanol, xylose, maltose, and rhamnose), and the model could not achieve growth with any individual carbon source. For ethanol, the reaction catalyzed by alcohol dehydrogenase (EC: 126.96.36.199) was irreversible, but should be reversible according to the MetaCyc database. For rhamnose, the lack of L-rhamnulose-1-phosphate lactaldehyde-lyase (EC: 188.8.131.52) meant that rhamnose could not generate glycerone phosphate, which could otherwise be used by the glycolysis pathway. After elimination of initial discrepancies using the continuous gap filling process, a 100% match was acquired. This indicated that model iCY1106 could predict catabolic pathways of various carbon sources including common sugars and alcohols. A similar FBA simulation of nitrogen sources (ammonium, nitrate, urea, and amino acids) also generated results that were consistent with experimental fermentation data. Both glutamate and glycine supported M. alpina growth, confirming that the model performed well for predicting the utilization of different carbon and nitrogen sources.
PUFA biosynthetic pathways based on GSMM
To further investigate the pathways associated with PUFA biosynthesis and metabolism, genome annotation, literature mining and KEGG database analysis were performed. ARA is a ω-6 PUFA synthesized by M. alpina (Figure 3). ARA (c204(6)) is synthesized from acetyl-CoA in 38 enzymatic steps, including the de novo synthesis of palmitic acid (c160) and the synthesis of very long chain fatty acids. During PUFA biosynthesis, all fatty acids should be in the form of acyl-CoA . Desaturase and elongase were necessary for the conversion of octadecanoyl-CoA (c180coa) to ARA. Δ9 desaturase  (MA-055-51, MA-101-533, MA-184-235), Δ12 desaturase  (MA-334-414) and Δ6 desaturase [30,31] (MA-268-54, MA-101-36) catalyzed the conversion of c180coa into γ-linolenoyl-CoA (c183(6)coa) which was subsequently lengthened to 8,11,14-eicosatrienoyl-CoA (c203(6)coa) by an elongase such as ELO/GLELO  (MA-189-257), MAELO  (MA-184-206) or MALCE1  (MA-073-327, MA-320-221). Δ5 desaturase  (MA-326-160) catalyzed c203(6)coa to arachidonyl-CoA (c204(6)coa), which was finally hydrolyzed to c204(6) by choloyl-CoA hydrolase (MA-049-4, MA-139-47).
M. alpina is capable of producing ARA and EPA, and EPA can be produced via the ω-3 or the ω-6 PUFA biosynthetic pathways through ARA desaturation . During EPA biosynthesis, Δ15 (ω3) desaturase (MA-072-10) is the key enzyme and was the first fungal desaturase known that uses both 18-carbon and 20-carbon ω-6 PUFAs as substrates . In the ω-3 pathway, (11Z,14Z,17Z)-icosatrienoyl-CoA (c204(3)coa) was used as substrate, while the ω-6 pathway utilized c204(6)coa.
Both the maximum and minimum flux of oxygen absorption of DGLA, ARA and EPA production in silico were found to be 1.552, 1.674 and 1.891 mmol gDW−1 h−1 respectively, using FVA. This indicated that the degree of PUFA unsaturation was dependent on oxygen content. Robustness analysis was used to elucidate the effects of oxygen uptake rate on PUFA biosynthesis. The oxygen uptake rate was perturbed by constraints-based flux analysis between 0 and 20 mmol gDW−1 h−1. The minimum growth rate was set at 0.03 h−1, with ARA exchange reaction as the objective function. When oxygen absorption was lower than 2.0 mmol gDW−1 h−1, ARA production increased with oxygen absorption rate (Additional file 5a). However, beyond this range, ARA production gradually declined with increasing oxygen uptake rate. As reported by Higashiyama et al., when the DO concentration range was 10-15 ppm, the ARA yield was enhanced 1.6-fold over a DO of 7 ppm, and a DO between 20 and 50 ppm drastically decreased ARA production . This indicated that excessive oxygen exposure during fermentation could impair ARA overproduction. During ARA biosynthesis, high oxygen concentrations may affect metabolism and cell growth in the filamentous mycelia, and beta-oxidation of fatty acids may be needed to obtain the extra energy required for adaptation to the high oxygen conditions.
Regulation of PUFA biosynthesis and metabolism based on GSMM
Acetyl-CoA is the essential two-carbon donor molecule for fatty acid synthesis. In model iCY1106, there were ten pathways that could produce acetyl-CoA, including amino acid degradation (from aspartate, isoleucine, leucine, and lysine), fatty acid beta-oxidation, and other acetyl-CoA generation processes (phosphoenolpyruvate [PEP], malate, L-lactate, acetate, and citrate) (Figure 4). However, only the pathway from PEP to pyruvate could generate sufficient flux for synthesis of PUFAs on the fermentation medium. During the growth stage, the flux of acetyl-CoA generated by pyruvate was 1.22 mmol gDW−1 h−1, while this value was 1.56 mmol gDW−1 h−1 during the product synthesis stage. This corresponds to a 27.87% increase, which would ensure acetyl-CoA was available for ARA production due to the enhanced activity of the pyruvate dehydrogenase complex (EC: 184.108.40.206, 220.127.116.11, 18.104.22.168) . Although aspartate, isoleucine, leucine, and lysine could all be used to generate acetyl-CoA for ARA production. Asparate has four carbon, while isoleucine, leucine, and lysine all contains six carbon. When the maximum uptake rate of asparate and the other amino acids were set at 0.15 mmol gDW−1 h−1 and 0.1 mmol gDW−1 h−1, respectively. FBA showed only that leucine and lysine could increase ARA production by 14.06% and 13.28%.
In model iCY1106, acetyl-CoA can be consumed by amino acid synthesis, amino sugar metabolism and fatty acid synthesis as well as ARA production (Additional file 5b). During the synthesis of six amino acids (threonine, methionine, lysine, cysteine, ornithine, and leucine), acetyl-CoA was required. FBA results showed that during the growth stage, 96.5% of the acetyl-CoA flux was used to synthesize fatty acids and 2.58% was used for amino acid synthesis. In contrast, during the product synthesis stage, the flux of acetyl-CoA used for fatty acid synthesis accounted for 99.52%, which indicates that amino acid synthesis was inhibited during ARA production. Limiting the nitrogen source can be an effective strategy to control amino acids biosynthesis , in which the lipid yield may be increased despite decreases in mycelia concentration.
During the growth stages, the flux of acetyl-CoA used to synthesize malonyl-CoA, catalyzed by acetyl-CoA carboxylase (ACC, EC 22.214.171.124), was 0.79 mmol gDW−1 h−1. In contrast during the ARA production stage, the corresponding flux was 1.43 mmol gDW−1 h−1 (81.0%). ACC catalyzes the first committal step of the fatty acid biosynthetic pathway, and should be overexpressed to maximize ARA production . Glutamate has been shown to activate ACC, and adding glutamate to the culture medium led to an increase in total lipid (21.81%) and ARA yield (66.07%) . As the M. alpina growth rate in batch culture ranged from 0.06 h−1 to 0.19 h-1 [1,7,8,41,42], the glutamate uptake rates used for the in silico analysis should be set between 0 and 1.8 mmol gDW−1 h−1. When the glutamate uptake rate was increased, ARA production increased from 0.128 mmol gDW−1 h−1 to 0.355 mmol gDW−1 h−1 (Additional file 5c).
During PUFA biosynthesis, NADPH is a necessary cofactor for de novo fatty acid synthesis involving desaturation and elongation. In model iCY1106, there were 172 reactions involving NADPH (Additional file 6). With ARA production fixed at 0.059 mmol gDW−1 h−1 (μ = 0.06 h−1) and 0.128 mmol gDW−1 h−1 (μ = 0.03 h−1), the flux of NADPH reactions were investigated. In total, 53 NADPH reactions exhibited flux changes of > 10−6 mmol gDW−1 h−1 (Figure 5).
Further research showed that only five reactions (R2, R3, R4, R5, and R6) involving accumulation of NADPH were significantly altered. When ARA production was increased, three reactions (R4, R5, and R6) underwent flux increases, which corresponded to enhanced enzyme activities due to up-regulation of the genes encoding the enzymes. Of these, malic enzyme (ME, EC 126.96.36.199), which catalyzes the conversion of malate to pyruvate (R4), is considered a key enzyme in lipogenesis in M. alpina. ME expression and enzyme activity are enhanced during ARA production , and the minimization of metabolic adjustment (MOMA) method was used to investigate the role of this enzyme. With the biomass and ARA flux maximized by ME deletion using MOMA, 410 (22.11%) reactions showed flux changes (Additional file 6) compared with the wild type model. The ARA exchange reaction flux was lowered from 0.128 mmol gDW−1 h−1 to 0.079 mmol gDW−1 h−1, a decrease of 38.28%. Additionally, all reactions involving NADPH consumption reactions other than R13, R16, R51, and R53 were associated with lipid metabolism. These results indicated that the increase in ARA production was directly correlated with the NADPH consumption rate.
In this study, a genome-scale metabolic model (iCY1106) for Mortierella alpina ATCC 32222 was successfully reconstructed using knowledge from the scientific literature and publicly accessible databases. The model comprised 1106 genes, 1854 reactions and 1732 metabolites, and included 247 transport reactions and 216 exchange reactions. Following a series of simulations and verification by measurement of growth rates and substrate usage, the model was found to agree with published literature. The model was used to investigate the effects of the important precursors acetyl-CoA and NADPH on the biosynthesis of PUFAs such as ARA and EPA. FBA results showed that enhancing the pyruvate dehydrogenase complex increased acetyl-CoA availability for increasing ARA production. Limiting the nitrogen source was an effective method for reducing acetyl-CoA consumption, and malic enzyme was found to be a key node in the regulation of NADPH in ARA biosynthesis. Model iCY1106 could serve as a useful predictive tool for future systems biology studies to guide the genetic engineering of M. alpina to improve the production of industrially important metabolites.
Construction of the M. alpina GSMM
The draft model was constructed by amassing reactions from the genome-scale metabolic models of genetically related organisms Aspergillus terreus , Pichia pastoris  and Yarrowia lipolytica . Reactions were chosen based on orthologs shared between M. alpina and the three reference organisms identified by protein sequence similarity searches using BLAST. Protein sequences from A. terreus NIH 2624, P. pastoris GS 115 and Y. lipolytica CLIB 122 were downloaded from UniProt . Open reading frame (ORF) information for M. alpina ATCC 32222 was provided by Yong Q. Chen based on sequencing results. The iJL1454 GSMM was used as a reference. Additionally, iLC915  and iYL619_PCP  were also used for comparison since both models were concerned with lipid metabolism and overproduction. To ensure accuracy, only sufficiently similar orthologs with e-values ≤10−30 and sequence identities ≥40% were included . To expand and update the draft model, the genome of M. alpina was re-annotated by submitting ORFs to the KAAS online annotation server . Metabolic reactions absent from the draft model were added from the KEGG database  based on KAAS annotation results. MetaCyc  and BioPath  databases were used to judge reaction reversibility. Compartmentalization information assigned to reactions was determined by subcellular localization prediction tools CELLO  and WoLF PSORT.  BaCelLo  was also used for proteins that were difficult to determine with the other tools. Transport information was obtained by cross-referencing BLATSp searches and the Transporter Classification Database TCDB .
To refine the draft model, the gapFind  program in the COBRA Toolbox  was used to identify metabolic gaps in draft model and literature data were used to fill these gaps. The metabolites in each reaction were characterized based on their chemical formulae and neutral charges, which were obtained using CHEBI  and PubChem .
The biomass equation of M. alpina was assumed to have six components: proteins, DNA, RNA, lipids, the cell wall and the small molecule pool [59,60]. Since no detailed information on M. alpina DNA and RNA was available, the ratio was assumed to be the same as in the related Aspergillus niger . The nucleotide and amino acid composition were calculated based on the M. alpina ATCC 32222 genome , as no specific experimental data were available. Similarly, the cell wall composition was calculated based on the typical fungal cell wall structure . The lipid composition was calculated based on the current literature . For calculation of energetic parameters, the growth and non-growth associated ATP maintenance values (GAM and NGAM, respectively) were assumed to be the same as those in the central carbon metabolic model of A. niger . Detailed biomass composition information can be found in Additional file 7.
Simulation and analysis
The reconstructed metabolic network was converted into stoichiometric matrix (S = M * N) using the Matlab program, where M represents metabolites and N represents reactions . The basic tools used for model analysis were flux balance analysis (FBA) and flux variability analysis (FVA). GLPK was used for linear programming , and Gurobi was used for quadratic programming . In silico analysis included growth simulation, gene and reaction essentiality analysis, robustness analysis, and minimization of metabolic adjustment. Analyses were performed according to the instructions for the COBRA Toolbox . Constraints used in each analysis are mentioned in the results section. To analyze model parameters relevant to cell growth, the biomass equation was selected as the objective function. For analysis of ARA production, the exchange reaction of ARA was the objective function.
Polyunsaturated fatty acids
Genome-scale metabolic model
NADH-cytochrome b5 reductase
Flux balance analysis
Flux variability analysis
Minimization of metabolic adjustment
Growth associated ATP maintenance
Non-growth associated ATP maintenance
Ho S-Y, Jiang Y, Chen F. Polyunsaturated fatty acids (PUFAs) content of the fungus Mortierella alpina isolated from soil. J Agric Food Chem. 2007;55:3960–6.
Wang L, Chen W, Feng Y, Ren Y, Gu Z, Chen H, et al. Genome characterization of the oleaginous fungus Mortierella alpina. PLoS One. 2011;6:e28319.
Shimizu S, Ogawa J, Kataoka M, Kobayashi M. Screening of novel microbial enzymes for the production of biologically and chemically useful compounds. Adv Biochem Eng Biotechnol. 1997;58:45–87.
Streekstra H. On the safety of Mortierella alpina for the production of food ingredients, such as arachidonic acid. J Biotechnol. 1997;56:153–65.
Jin MJ, Huang H, Xiao AH, Zhang K, Liu X, Li S, et al. A novel two-step fermentation process for improved arachidonic acid production by Mortierella alpina. Biotechnol Lett. 2008;30:1087–91.
Sakuradani E, Hirano Y, Kamada N, Nojiri M, Ogawa J, Shimizu S. Improvement of arachidonic acid production by mutants with lower n-3 desaturation activity derived from Mortierella alpina 1S-4. Appl Microbiol Biotechnol. 2004;66:243–8.
Zhang J, Zhao JX, Chen HQ, Chen W, Zhang H. The preliminary study of fermentation conditions for polyunsaturated fatty acids production by Mortierella aplina. Sci Technol Food Int. 2011;11:216–9.
Peng C, Huang H, Ji X, Liu X, You J, Lu J, et al. A temperature-shift strategy for efficient arachidonic acid fermentation by Mortierella alpina in batch culture. Biochem Eng J. 2010;53:92–6.
Takeno S, Sakuradani E, Tomi A, Inohara-Ochiai M, Kawashima H, Shimizu S. Transformation of oil-producing fungus, Mortierella alpina 1S-4, using zeocin, and application to arachidonic acid production. J Biosci Bioeng. 2005;100:617–22.
Abe T, Sakuradani E, Ueda T, Shimizu S. Identification of mutation sites on delta5 desaturase genes from Mortierella alpina 1S-4 mutants. J Biosci Bioeng. 2005;99:296–9.
Wynn JP, Ratledge C. Evidence that the rate-limiting step for the biosynthesis of arachidonic acid in Mortierella alpina is at the level of the 18:3 to 20:3 elongase. Microbiology. 2000;146:2325–31.
Takeno S, Sakuradani E, Murata S, Inohara-Ochiai M, Kawashima H, Ashikari T, et al. Molecular evidence that the rate-limiting step for the biosynthesis of arachidonic acid in Mortierella alpina is at the level of an elongase. Lipids. 2005;40:25–30.
Certik M, Sakuradani E, Kobayashi M, Shimizu S. Characterization of the second form of NADH-cytochrome b(5) reductase gene from arachidonic acid-producing fungus Mortierella alpina 1S-4. J Biosci Bioeng. 1999;88:667–71.
Sakuradani E, Ando A, Ogawa J, Shimizu S. Improved production of various polyunsaturated fatty acids through filamentous fungus Mortierella alpina breeding. Appl Microbiol Biotechnol. 2009;84:1–10.
Mahadevan R, Schilling CH. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003;5:264–76.
Varma A, Palsson BO. Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol. 1994;60:3724–31.
Barrett CL, Herring CD, Reed JL, Palsson BO. The global transcriptional regulatory network for metabolism in Escherichia coli exhibits few dominant functional states. Proc Natl Acad Sci U S A. 2005;102:19103–8.
Pharkya P, Burgard AP, Maranas CD. OptStrain: A computational framework for redesign of microbial production systems. Genome Res. 2004;14:2367–76.
Loira N, Dulermo T, Nicaud J-M, Sherman DJ. A genome-scale metabolic model of the lipid-accumulating yeast Yarrowia lipolytica. BMC Syst Biol. 2012;6:35–43.
Pan P, Hua Q. Reconstruction and in silico analysis of metabolic network for an oleaginous yeast, Yarrowia lipolytica. PLoS One. 2012;7:e51535.
Vongsangnak W, Ruenwai R, Xin T, Xinjie H, Hao Z, Shen B, et al. Genome-scale analysis of the metabolic networks of oleaginous Zygomycete fungi. Gene. 2013;521:180–90.
Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.
Lu J, Peng C, Ji XJ, You J, Cong L, Ouyang P, et al. Fermentation characteristics of Mortierella alpina in response to different nitrogen sources. Appl Biochem Biotechnol. 2011;164:979–90.
Yu A-Q, Zhu J-C, Zhang B, Xing L-J, Li M. Effects of different carbon sources on the growth, fatty acids production, and expression of three desaturase genes of mortierella alpina ATCC 16266. Curr Microbiol. 2011;62:1617–22.
Nisha A, Venkateswaran G. Effect of culture variables on mycelial arachidonic acid production by Mortierella alpina. Food Bioprocess Tech. 2011;4:232–40.
Boghigian BA, Armando J, Salas D, Pfeifer BA. Computational identification of gene over-expression targets for metabolic engineering of taxadiene production. Appl Microbiol Biotechnol. 2012;93:2063–73.
Cong L, Peng C, Ji X, Li Z, You J, Lu J, et al. Progress in production of arachidonic acid by Mortierella alpina and genetic modification. Chin J Biotechnol. 2010;26:1232–8.
Abe T, Sakuradani E, Asano T, Kanamaru H, Shimizu S. Functional characterization of Δ9 and ω9 desaturase genes in Mortierella alpina 1S-4 and its derivative mutants. Appl Microbiol Biotechnol. 2006;70:711–9.
Sakuradani E, Abe T, Matsumura K, Tomi A, Shimizu S. Identification of mutation sites on delta 12 desaturase genes from Mortierella alpina 1S-4 mutants. J Biosci Bioeng. 2009;107:99–101.
Chen R, Matsui K, Ogawa M, Oe M, Ochiai M, Kawashima H, et al. Expression of delta 6, delta 5 desaturase and GLELO elongase genes from Mortierella alpina for production of arachidonic acid in soybean glycine [max (L.) merrill] seeds. Plant Sci. 2006;170:399–406.
Sakuradani E, Shimizu S. Gene cloning and functional analysis of a second delta 6-fatty acid desaturase from an arachidonic acid-producing Mortierella fungus. Biosci Biotechnol Biochem. 2003;67:704–11.
Sakuradani E, Murata S, Kanamaru H, Shimizu S. Functional analysis of a fatty acid elongase from arachidonic acid-producing Mortierella alpina 1S-4. Appl Microbiol Biotechnol. 2008;81:497–503.
Sakuradani E, Nojiri M, Suzuki H, Shimizu S. Identification of a novel fatty acid elongase with a wide substrate specificity from arachidonic acid-producing fungus Mortierella alpina 1S-4. Appl Microbiol Biotechnol. 2009;84:709–16.
Bajpai P, Bajpai PK, Ward OP. Eicosapentaenoic acid (EPA) Production by Mortierella alpina ATCC 32222. Appl Biochem Biotechnol. 1991;31:267–72.
Sakuradani E, Abe T, Iguchi K, Shimizu S. A novel fungal omega 3-desaturase with wide substrate specificity from arachidonic acid-producing Mortierella alpina 1S-4. Appl Microbiol Biotechnol. 2005;66:648–54.
Higashiyama K, Murakami K, Tsujimura H, Matsumoto N, Fujikawa S. Effects of dissolved oxygen on the morphology of an arachidonic acid production by Mortierella alpina 1S-4. Biotechnol Bioeng. 1999;63:442–8.
Pettit FH, Pelley JW, Reed LJ. Regulation of pyruvate dehydrogenase kinase and phosphatase by acetyl-CoA/CoA and NADH/NAD ratios. Biochem Biophys Res Commun. 1975;65:575–82.
Koike Y, Cai HJ, Higashiyama K, Fujikawa S, Park EY. Effect of consumed carbon to nitrogen ratio of mycelial morphology and arachidonic acid production in cultures of Mortierella alpina. J Biosci Bioeng. 2001;91:382–9.
Davis MS, Solbiati J, Cronan Jr JE. Overproduction of acetyl-CoA carboxylase activity increases the rate of fatty acid biosynthesis in Escherichia coli. J Biol Chem. 2000;275:28593–8.
Yu L, Qin W, Lan W, Zhou P, Zhu M. Improved arachidonic acids production from the fungus Mortierella alpina by glutamate supplementation. Bioresour Technol. 2003;88:265–8.
Eroshin V, Satroutdinov A, Dedyukhina E, Chistyakova T. Arachidonic acid production by Mortierella alpina with growth-coupled lipid synthesis. Process Biochem. 2000;35:1171–5.
Zeng Y, Ji X, Chang S, Nie Z, Huang H. Improving arachidonic acid accumulation in Mortierella alpina through B-group vitamin addition. Bioprocess Biosystems Eng. 2012;35:683–8.
Wynn JP, Bin Abdul Hamid A, Ratledge C. The role of malic enzyme in the regulation of lipid accumulation in filamentous fungi. Microbiology. 1999;145:1911–7.
Liu J, Gao Q, Xu N, Liu L. Genome-scale reconstruction and in silico analysis of Aspergillus terreus metabolism. Mol Biosyst. 2013;9:1939–48.
Caspeta L, Shoaie S, Agren R, Nookaew I, Nielsen J. Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials. BMC Syst Biol. 2012;6:24–37.
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, et al. UniProt: the Universal Protein knowledgebase. Nucleic Acids Res. 2004;32:D115–9.
Tian W, Skolnick J. How well is enzyme function conserved as a function of pairwise sequence identity? J Mol Biol. 2003;333:863–82.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–5.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, et al. The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 2008;36:D623–31.
Reitz M, Sacher O, Tarkhov A, Trumbach D, Gasteiger J. Enabling the exploration of biochemical pathways. Org Biomol Chem. 2004;2:3226–37.
Yu CS, Lin CJ, Hwang JK. Predicting subcellular localization of proteins for gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004;13:1402–6.
Horton P, Park K-J, Obayashi T, Fujita N, Harada H, Adams-Collier C, et al. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35:W585–7.
Pierleoni A, Martelli PL, Fariselli P, Casadio R. BaCelLo: a balanced subcellular localization predictor. Bioinformatics. 2006;22:e408–16.
Saier Jr MH, Tran CV, Barabote RD. TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006;34:D181–6.
Kumar VS, Dasika MS, Maranas CD. Optimization based automated curation of metabolic reconstructions. BMC Bioinformatics. 2007;8:212–27.
Degtyarenko K, de Matos P, Ennis M, Hastings J, Zbinden M, McNaught A, et al. ChEBI: a database and ontology for chemical entities of biological interest. Nucleic Acids Res. 2008;36:D344–50.
Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH. PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009;37:W623–33.
Nisha A, Muthukumar SP, Venkateswaran G. Safety evaluation of arachidonic acid rich Mortierella alpina biomass in albino rats-A subchronic study. Regul Toxicol Pharmacol. 2009;53:186–94.
Shinmen Y, Kawashima H, Shimizu S, Yamada H. Concentration of eicosapentaenoic acid and docosahexaenoic acid in an arachidonic acid-producing fungus, Mortierella alpina IS-4, grown with fish oil. Appl Microbiol Biotechnol. 1992;38:301–4.
Andersen MR, Nielsen ML, Nielsen J. Metabolic model integration of the bibliome, genome, metabolome and reactome of Aspergillus niger. Mol Syst Biol. 2008;4:178.
Bowman SM, Free SJ. The structure and synthesis of the fungal cell wall. Bioessays. 2006;28:799–808.
Jareonkitmongkol S, Sakuradani E, Shimizu S. Isolation and characterization of a Delta 9-desaturation-defective mutant of an arachidonic acid-producing fungus, Mortierella alpina 1S-4. J Am Oil Chem Soc. 2002;79:1021–6.
David H, Akesson M, Nielsen J. Reconstruction of the central carbon metabolism of Aspergillus niger. Eur J Biochem. 2003;270:4243–53.
Thiele I, Palsson BO. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5:93–121.
Yin W: Gurobi Mex: a MATLAB interface for Gurobi. Online at http://convexoptimization.com/wikimization/index.php/gurobi_mex 2009.
The author would like to thank Dr Jens Nielsen, Professor at Chalmers University of Technology, for technological help on the model reconstruction. This work was supported by the Program for Innovative Research Team in University (IRT1249), the Provincial Outstanding Youth Foundation of Jiangsu Province (BK2012002), the National Natural Science Foundation of China (21422602), the Program for Young Talents in China, the National High-Tech R & D Program of China (No. 2014AA021701). The Program for New Century Excellent Talents (NCET-13-0831), the National Natural Science Foundation of China (No. 31125021, No. 21276108), and the Fundamental Research Funds for the Central Universities (No. JUSRP51320B).
The authors declare that they have no competing interests.
CY carried out the construction of iCY1106, performed the statistical analysis, and drafted the manuscript. NX carried out literature mining for M. alpina, provided help for the construct of model. HC, YQC and WC provided the refined genome sequence of M. alpina. LL participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.
Detailed information about model i CY1106.
Shared reactions between Vongsangnak’s model and i CY1106 (not including transport and exchange reactions, no compartments information).
Result of essential genes identified by FBA using minimum growth medium and yeast extract medium.
Simulation result for the capability of usage different carbon and nitrogen sources for i CY1106.
a Figure for the robustness analysis about oxygen and ARA. b Figure about the acetyl-CoA consumption pathways in iCY1106. c Figure for the robustness analysis about glutamate and ARA.
FBA simulation result for the reactions involving NADPH at different ARA levels.
Biomass composition information for M. alpina.