Proteomics analysis of high lipid-producing strain Mucor circinelloides WJ11: an explanation for the mechanism of lipid accumulation at the proteomic level
© Tang et al. 2016
Received: 16 October 2015
Accepted: 19 January 2016
Published: 11 February 2016
The oleaginous fungus, Mucor circinelloides, is attracting considerable interest as it produces oil rich in γ-linolenic acid. Nitrogen (N) deficiency is a common strategy to trigger the lipid accumulation in oleaginous microorganisms. Although a simple pathway from N depletion in the medium to lipid accumulation has been elucidated at the enzymatic level, global changes at protein levels upon N depletion have not been investigated. In this study, we have systematically analyzed the changes at the levels of protein expression in M. circinelloides WJ11, a high lipid-producing strain (36 %, lipid/cell dry weight), during lipid accumulation.
Proteomic analysis demonstrated that N depletion increased the expression of glutamine synthetase, involved in ammonia assimilation, for the supply of cellular nitrogen but decreased the metabolism of amino acids. Upon N deficiency, many proteins (e.g., fructose-bisphosphate aldolase, glyceraldehyde-3-phosphate dehydrogenase, enolase, pyruvate kinase) involved in glycolytic pathway were up-regulated while proteins involved in the tricarboxylic acid cycle (e.g., isocitrate dehydrogenase, succinyl-CoA ligase, succinate dehydrogenase, fumarate hydratase) were down-regulated, indicating this activity was retarded thereby leading to a greater flux of carbon into fatty acid biosynthesis. Moreover, glucose-6-phosphate dehydrogenase, transaldolase and transketolase, which participate in the pentose phosphate pathway, were up-regulated, leading to the increased production of NADPH, the reducing power for fatty acid biosynthesis. Furthermore, protein and nucleic acid metabolism were down-regulated and some proteins involved in energy metabolism, signal transduction, molecular chaperone and redox homeostasis were up-regulated upon N depletion, which may be the cellular response to the stress produced by the onset of N deficiency.
N limitation increased those expressions of the proteins involved in ammonia assimilation but decreased that involved in the biosynthesis of amino acids. Upon N deprivation, the glycolytic pathway was up-regulated, while the activity of the tricarboxylic acid cycle was retarded, thus, leading more carbon flux to fatty acid biosynthesis. Moreover, the pentose phosphate pathway was up-regulated, then this would increase the production of NADPH. Together, coordinated regulation of central carbon metabolism upon N limitation, provides more carbon flux to acetyl-CoA and NADPH for fatty acid biosynthesis.
Oleaginous microorganisms include fungi, yeasts, microalgae and bacteria that can accumulate oil to more than 20 % of their cell dry weight (CDW) . The metabolism of lipid accumulation in oleaginous microorganisms has been extensively studied as microbial oils can be used as commercial sources of several nutritionally-important polyunsaturated fatty acids (PUFAs) and as potential sources of biofuels [2, 3]. Among oleaginous filamentous fungi, Mucor circinelloides has been considered as an important model organism for lipid accumulation studies due to its ability to produce an oil rich in γ-linolenic acid (GLA, 18:3; n-6), that may have beneficial effects for the treatment of premenstrual tension, atopic dermatitis and some other diseases  and also due to the availability of genome data and genetic tools.
Lipid accumulation in oleaginous microorganisms is triggered by a nutrient imbalance in the culture medium. When cells run out of a key nutrient, usually nitrogen (N), excess carbon substrate continues to be assimilated by the cells and converted into storage lipids . The biochemistry of lipid accumulation in oleaginous microorganisms triggered by N deficiency has been widely investigated. Previous data indicated that the activity of isocitrate dehydrogenase decreases rapidly and even ceases completely under N-deficient conditions as the diminishing concentration of its allosteric regulator AMP, and then results in a shift in carbon flux through the citric acid cycle and into lipid biosynthesis [5, 6]. In addition, ATP: citrate lyase (ACL), which generates acetyl-CoA as the precursor of fatty acids via the cleavage of citric acid, is an essential enzyme for fatty acid biosynthesis and possibly catalyzes the rate-limiting reaction for fatty acid biosynthesis in some oleaginous organism [7–9].
Beside acetyl-CoA, the provision of reducing power in the form of NADPH is another critical process for fatty acid biosynthesis. The key roles of malic enzyme and the pentose phosphate pathway (glucose-6-phosphate dehydrogenase coupled with 6-phosphogluconate dehydrogenase) to supply NADPH for fatty acid synthesis during lipid accumulation have been proposed [2, 10–17]. However, lipid accumulation is a complicated process involving many metabolic pathways and thus it is impossible to achieve maximal lipid production simply by regulating just one or two genes. An investigation of lipid metabolism at systematic level is therefore required to gain insights into the molecular mechanism of lipid accumulation.
Cellular responses to N deficiency are the subtle behaviors of living organisms. Although there are alternative approaches to understand the molecular mechanisms of cellular response under N deficiency, microbial proteomics has become a powerful tool to investigate the complex cellular processes. In addition, it can also determine new functions of gene products as it represents not only the gene product, but also translational rate and post-translational modifications. Identification of proteins that are up/down-regulated under N deficiency is important in studying the mechanisms of lipid accumulation. Accordingly, comparative proteomics have been explored to gain insights into the lipid metabolism under N deficiency in microalgae and yeast [18–21].
M. circinelloides was the first microorganism to be used commercially to produce an oil for human consumption—an oil rich in GLA . However, the process began in 1985 but lasted only 6 years as the commercial strain only produced 25 % lipid of its cell biomass and at that time high GLA-producing plant species came on to the market [3, 4]. In this study, we performed a comparative proteomic study on the oleaginous fungus M. circinelloides WJ11, which is the highest lipid-producing strain (up to 36 % lipid, w/w) of this species as far and its lipid content is much higher than that commercial strain . M. circinelloides WJ11 could be the potential strain to produce the commercial oil rich in GLA and its mechanism of lipid accumulation at proteomic level will provide a foundation of restarting the commercial production of GLA by microorganism. We compared proteomes from three growth stages (the balanced growth stage, the fast lipid accumulation stage and the slow lipid accumulation stage) to provide new insights into the mechanism of lipid metabolism in this fungus. To the best of our knowledge, this is the first proteomic study on the lipid overproduction process of the oleaginous filamentous fungus, M. circinelloides WJ11, under N deficiency.
Results and discussion
Cell growth and lipid accumulation in M. circinelloides WJ11
Average lipid biosynthesis rate in M. circinellodies WJ11 during the bioprocess
Average lipid biosynthesis ratea (mg g−1 h−1)
From 6 to 9 h
From 9 to 12 h
From 12 to 24 h
From 24 to 36 h
From 36 to 48 h
From 48 to 60 h
From 60 to 72 h
From 72 to 84 h
From 84 to 96 h
Proteome analysis of M. circinelloides during lipid accumulation upon N deficiency
Identification of differentially expressed proteins under N deficiency in M. circinellodies WJ11 by MALDI-TOF/TOF analysis
Theoretical MW (Da)/pI
Fold change 24 h/6 h
Fold change 60 h/6 h
Amino acid metabolism
Branched-chain amino acid aminotransferase
Serine hydroxymethyltransferase 2
Serine hydroxymethyltransferase, mitochondrial
Glycine cleavage system T protein
Glyceraldehyde-3-phosphate dehydrogenase 1
Glyceraldehyde-3-phosphate dehydrogenase 1
Glyceraldehyde-3-phosphate dehydrogenase 1
Glyceraldehyde-3-phosphate dehydrogenase 2
Pyruvate dehydrogenase complex
Pyruvate dehydrogenase complex
Pyruvate dehydrogenase E1 component subunit alpha
Pyruvate dehydrogenase E2 component
Pyruvate dehydrogenase E2 component
2,3-Bisphosphoglycerate-independent phosphoglycerate mutase
2,3-Bisphosphoglycerate-independent phosphoglycerate mutase
NAD+:isocitrate dehydrogenase, mitochondrial
NADP+:isocitrate dehydrogenase, mitochondrial
Succinyl-CoA ligase, mitochondrial
Fumarate hydratase, mitochondrial
Thiazole biosynthetic enzyme
40S ribosomal protein S21
40S ribosomal protein S7
Eukaryotic translation initiation factor 5A
Elongation factor Tu
Nascent polypeptide-associated complex subunit beta
26S protease regulatory subunit 6A
26S protease regulatory subunit 6B
20S proteasome subunit Alpha 6
Nucleic acid metabolism
Dihydroorotase, homodimeric type
Adenylate kinase 1
Adenylate kinase 1
V-type H+-transporting ATPase subunit I
V-type proton ATPase catalytic subunit A
14-3-3 family protein epsilon
Heat shock 70 kDa protein 1/8
Peptidyl-prolyl cis–trans isomerase cyp5
Ferritin heavy chain
Hypothetical protein HMPREF1544_00199
Hypothetical protein HMPREF1544_01177
Hypothetical protein HMPREF1544_01386
Hypothetical protein HMPREF1544_01888
Hypothetical protein HMPREF1544_01888
Hypothetical protein HMPREF1544_02101
Hypothetical protein HMPREF1544_02181
Hypothetical protein HMPREF1544_02260
Hypothetical protein HMPREF1544_03035
Hypothetical protein HMPREF1544_03136
Hypothetical protein HMPREF1544_03262
Hypothetical protein HMPREF1544_03549
Hypothetical protein HMPREF1544_03639
Hypothetical protein HMPREF1544_03790
Hypothetical protein HMPREF1544_03919
Hypothetical protein HMPREF1544_03956
Hypothetical protein HMPREF1544_04707
Hypothetical protein HMPREF1544_05024
Hypothetical protein HMPREF1544_05024
Hypothetical protein HMPREF1544_05863
Hypothetical protein HMPREF1544_06871
Hypothetical protein HMPREF1544_06986
Hypothetical protein HMPREF1544_07306
Hypothetical protein HMPREF1544_07438
Hypothetical protein HMPREF1544_07752
Hypothetical protein HMPREF1544_07893
Hypothetical protein HMPREF1544_08616
Hypothetical protein HMPREF1544_08758
Hypothetical protein HMPREF1544_08758
Hypothetical protein HMPREF1544_09175
Hypothetical protein HMPREF1544_09235
Hypothetical protein HMPREF1544_09235
Hypothetical protein HMPREF1544_10038
Hypothetical protein HMPREF1544_10287
Hypothetical protein HMPREF1544_10622
Hypothetical protein HMPREF1544_10686
Hypothetical protein HMPREF1544_10852
Hypothetical protein HMPREF1544_10852
Hypothetical protein HMPREF1544_10852
Hypothetical protein HMPREF1544_11327
Hypothetical protein HMPREF1544_11393
Nitrogen and amino acid metabolism
Glutamine synthase (GS) is the key enzyme involved in ammonia assimilation in both plants and microorganisms [25, 26]. Previous proteomic analysis has indicated that its expression is induced immediately following N deprivation in both Phaeodactylum tricornutum and Rhodosporidium toruloides [18, 21]. Indeed, our result showed its expression (spots 1 and 2) in M. circinelloides was up-regulated upon N depletion. In oleaginous yeasts, upon N deficiency and at the beginning of lipid accumulation, AMP is deaminated to release ammonium for further cell use. Thus, the substrate NH4 + for GS might be from the degradation of AMP or other N compounds. Saccharopepsin (spots 3 and 4), which links to vacuolar protein degradation, was also induced upon N depletion in this study. The elevated expression of this protein would enhance degradation of proteins suggesting that, immediately after N exhaustion, the cells initiate the turnover and recycling of intracellular components-especially proteins that are no longer needed for anabolic reactions .
Most of the proteins associated with amino acid metabolism showed decreased expression level after N depletion. Acetolactate synthase (spot 5), dihydroxy-acid dehydratase (spot 6) and ketol-acid reductoisomerase (spot 7) catalyze the synthesis of the branched-chain amino acids (BCAA, e.g., valine, leucine, and isoleucine) [27–29], and branched-chain amino acid aminotransferase (spot 8) participates in the degradation BCAA . These proteins related to BCAA metabolism were all down-regulated upon N depletion. Acetylornithine aminotransferase (spot 9) and asparagine synthetase (spot 10) are involved in the asparagine biosynthesis . Saccharopine dehydrogenase (spot 11) and homoisocitrate dehydrogenase (spot 12) participate in lysine biosynthesis [32, 33]. Serine hydroxymethyltransferase (spots 13 and 15) take part in serine degradation and glycine cleavage system T protein (spot 15) is involved in glycine degradation . 3-deoxy-7-phosphoheptulonate synthase (spot 16) is responsible for the biosynthesis of phenylalanine, tyrosine, and tryptophan. These above proteins related to metabolism of amino acids were also down-regulated upon N depletion. Decreased expression of the proteins associated with amino acid metabolism indicated that amino acid biosynthesis was at least partially inhibited due to the absence of nitrogen. In these proteins related to amino acid metabolism, one exception to the overall decreased expression was S-adenosylmethionine synthase (SAS, spots 17 and 18) which participates in the biosynthesis of S-adenosylmethionine, which is a precursor of glutathione and other methylated cell components .
N deprivation can result in a metabolic imbalance of reactive oxygen species (ROS) and an excess of oxygen radicals , while glutathione can quench free radicals and then improve the resistance of the cells to this stress . Therefore, the increased expression of SAS could play a role in providing glutathione to improve stress resistance arising during N deficiency.
The glycolytic pathway provides pyruvate, a key precursor for acetyl-CoA, which is the substrate for fatty acid biosynthesis. Fructose-bisphosphate aldolase (FBA, spots 19 and 20), an essential enzyme involved in glycolysis that catalyzes a reversible cleavage reaction of fructose 1,6-bisphosphate into glyceraldehyde 3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) , was up-regulated under N deficiency. The increased FBA activity could stimulate glycolysis and triacylglycerol biosynthesis. Glyceraldehyde-3-phosphate dehydrogenase (spots 21, 22, 23 and 24), enolase (spots 25) and pyruvate kinase (spots 26 and 27), which are key enzymes in the glycolytic pathway, were also up-regulated under N deficiency. Up-regulation of these proteins may also lead to the increased production of pyruvate. In addition, the expression of pyruvate dehydrogenase (PDH, spots 28, 29, 30, 31 and 32) was increased and thus lead to increased conversion of pyruvate into acetyl-CoA in the mitochondrion. Taken together, under N deficiency, the expression of many proteins involved in glycolytic pathway was increased, which would then provide more acetyl-CoA in the mitochondrion.
Phosphoglycerate mutase (PGM, spots 33 and 34) is also involved in glycolysis, catalyzing the interconversion of 3-phosphoglycerate (3PGA) to 2-phosphoglycerate (2PGA). In this study, PGM was down-regulated when the fungus was grown under N deficient condition, as FBA was up-regulated, the expression changes of these proteins could lead to accumulation of some intermediate between fructose-1,6-bisphosphate and 3-phosphoglyceric acid and this could be glycerol.
The pentose phosphate pathway (PPP) generates NADPH, which is an important source besides malic enzyme (ME) for providing reducing power for fatty acid biosynthesis . In this study, glucose-6-phosphate dehydrogenase (G6PDH, spot 35), which is a key enzyme in the PPP that generates NADPH, was up-regulated under N deficiency. This is in accordance with our previous report showing that lipid accumulation in M. circinelloides is accompanied by increased G6PDH activity after N depletion from the medium . Furthermore, transaldolase (spot 36) and transketolase (spot 37), as part of the PPP, were both up-regulated. These results indicated that PPP was more active during the lipid accumulation phase after N depletion, and this would provide more NADPH for fatty acid biosynthesis. For ME, the critical enzyme plays a key role in supplying NADPH for fatty acid synthesis and desaturation in M. circinelloides [10, 38, 40, 41], we did not find any change in its expression at proteomic level upon N deficiency, which indicates the ME gene is being expressed all the time, irrespective of the status of the cells, and regulation of ME maybe complicated.
Oleaginous yeasts deaminate AMP to release ammonium and IMP upon nitrogen limitation . NAD+:isocitrate dehydrogenase (NAD+:ICDH) requires AMP for activity. With ICDH activity being curtailed by the lack of AMP, isocitrate accumulates and equilibrates back to citrate which is then transported out of the mitochondrion into the cytosol, and thus provides the increased carbon flux to acetyl-CoA for fatty acid synthesis . In this study, the NAD+:ICDH (spot 38) was down-regulated upon N deficiency which would increase citrate accumulation and then lead to greater carbon flux to acetyl-CoA for fatty acid synthesis. Mitochondrial NADP+:ICDH (spot 39) was also down-regulated, and this regulation was similar to NAD+:ICDH. In addition, the expression of some other proteins (e.g., succinyl-CoA ligase, spot 40; succinate dehydrogenase, spot 41; fumarate hydratase, spot 42) involved in TCA cycle were all decreased, and this further suggests that the TCA cycle is retarded upon N deficiency which will therefore lead to a greater carbon flux to lipids biosynthesis. These results are well in accordance with our previous studies of lipid accumulation in M. circinelloides and the multi-omic analysis of lipid accumulation in Rhodosporidium toruloides by Zhu et al. [6, 12].
Aldehyde dehydrogenase (ALDH, spots 43 and 44), which converts acetaldehyde into acetate, was up-regulated upon N deficiency. Although in oleaginous fungi and higher eukaryotes, the major route of acetyl-CoA production is by the cleavage of citrate by ATP:citrate lyase (ACL) , it can also be produced by cytoplasmic acetyl-CoA synthase (ACS) from acetate when acetate is being produced by the degradation of proteins and amino acids that are surplus to the requirements of the cell. The level of acetyl-CoA C-acetyltransferase (ACAT, spot 187), which is the branch point enzyme for acetyl-CoA to synthesize isoprenoids, was decreased under N deficiency which indicates that the acetyl-CoA flux is being preferentially switched into the synthesis of fatty acids but not into isoprenoids.
Galactokinase (spot 49) catabolizes β-D-galactose to glucose 1-phosphate. Phosphoglucomutase (spot 46) facilitates the interconversion of glucose 1-phosphate and glucose 6-phosphate. UTP-glucose-1-phosphate uridylyltransferase (spot 47) synthesizes UDP-glucose from glucose-1-phosphate. UDP-glucose 4-epimerase (spot 48) catalyzes the reversible conversion of UDP-galactose to UDP-glucose. These proteins are involved in galactose and glucose conversion into polysaccharides and were all up-regulated upon N deficiency, suggesting that glucose utilization of the fungus might be sequentially proceeded during lipid accumulation triggered by N depletion.
S-Formylglutathione hydrolase (spot 50) can produce glutathione from hydrolysis of S-formylglutathione and its expression was up-regulated upon N deficiency. N deprivation can result in an excess of oxygen radicals and glutathione can quench free radicals . Therefore the up-regulated protein could be a stress response to N deficiency. Thiazole biosynthetic enzyme (TBE, spot 51) associated with thiamine metabolism was also up-regulated upon N depletion. It plays additional roles in adaptation to various stress conditions and in DNA damage tolerance [43, 44]. Thus the increased expression of TBE in the fungus is likely the cellular adaptation to N depletion.
As de novo protein biosynthesis is no longer occurring after N deficiency, then the cell must conserve its key proteins for as long as possible. Thus, the expression of proteins (spots 52–59) related to protein metabolism were diminished in the fungus under N deficient conditions. In addition, the expression of proteins involved in nucleic acid metabolism (spots 60, 61 and 62) were also decreased and this further indicates that cell reproduction and growth are decreased upon N depletion .
The activity of adenylate kinase has been found to be stimulated after N-exhaustion in M. circinelloides  and we found that, in accordance with this, the expression of adenylate kinase (spots 63 and 64) was increased upon N deficiency. Furthermore, ATPase (spots 65 and 66) was also up-regulated. The metabolic pattern of M. circinelloides is therefore readjusted when N becomes exhausted in the medium. In M. circinelloides, the concentrations of ATP, ADP and AMP decrease at the time of N-exhaustion so that energy, in the form of ATP, is now produced at a lower level [20, 45]. Therefore the increased expression of adenylate kinase and ATPase might play a role in helping to maintain energy production.
14-3-3 Family protein epsilon (spot 67) can restrain cell apoptosis and promote cell survival under stress condition  and the increased expression in N limitation could be the cellular response to N deficiency. Heat shock proteins (HSPs) are a group of functionally related proteins responsible for protein folding and unfolding. High-level expression of HSPs can be triggered by exposure to different environmental stress conditions, including exposure of the cell to nitrogen deficiency . We identified similar trends in N-deprived cells, including two HSPs (spots 68 and 69), which were significantly up-regulated. Peptidyl-prolyl cis–trans isomerase plays roles in protein folding and transport, RNA splicing and the regulation of multi-protein complexes in cells . In this study, expression of this protein (spot 70) was up-regulated in N deprivation, and this might also be part of the cellular adaption to stress conditions engendered by N deficiency.
Excessive generation of ROS or oxidative stress is an integral part of many stress situations, including N limitation . Peroxiredoxin (Prx) is a ubiquitous family of antioxidant enzymes and glutathione peroxidase (GPx) has the biological role in protecting organism from oxidative damage. Indeed, the Prx (spots 71, 72 and 73) and GPx (spot 74) were up-regulated upon N deficiency. The increased expression of an oxidoreductase (spots 75 and 76) that catalyzes the transfer of electrons from reductant to oxidant is also a response to oxidative stress upon N limitation. Ferritin heavy chain (spot 77), a ubiquitous and highly conserved protein, which plays a major role in iron homeostasis, was also up-regulated upon N deficiency. Many haem-proteins would be degraded as being surplus to the requirements of the cell during N deprivation. The iron being released from these proteins and other non-haem iron proteins will therefore be scavenged by the cell and stored intracellularly in ferritin. This may lead to the increased expression of ferritin heavy chain.
Analysis of the transcription of selected genes by quantitative RT-PCR
The mRNA of the genes encoding PDH were increased upon N deficiency. This is coincident with expression level of the protein, suggesting that conversion of pyruvate to acetyl CoA for TCA cycle was increased upon N depletion. Transcription levels of the genes encoding NAD+:ICDH were decreased upon N depletion indicating that TCA cycle is retarded after N exhaustion. Thus, the combined coordinated regulation of PDH and NAD+:ICDH upon N depletion results in an increased cytosolic acetyl-CoA production for fatty acid biosynthesis. Furthermore, the mRNA of the genes encoding ACAT were also decreased upon N limitation, which is in according with its protein level. This will direct the flux of acetyl-CoA into fatty acid biosynthesis rather than into isoprenoid biosynthesis. Taken together, the quantitative RT-PCR analysis demonstrated the transcription level of some proteins related to lipid accumulation under N deficiency is consistant with their protein level from proteomic data.
M. circinelloides is a GLA-producing filamentous fungus and fatty acid delta-6 desaturase converts linoleic acid into GLA, which is an key enzyme for GLA biosynthesis . Thus we also investigated transcription levels of the genes (d61 and d62) encoding fatty acid delta-6 desaturase (Fig. 4). The result showed the transcription level of d61 was significantly lower at 24 h and 60 h than that at 6 h, while the level of d62 had no obvious changes in whole bioprocess. The fatty acid composition showed the GLA content in total fatty acid at 24 h and 60 h were both less than half of that at 6 h (Additional file 1: Table S1), which was in according with the transcription level of d61.
This study represents a proteomic analysis of lipid accumulation in M. circinelloides WJ11, a higher lipid-producing strain (36 % lipid, w/w), grown under N limited condition. N limitation increased the expression of the proteins involved in ammonia assimilation for the supply of cellular nitrogen but decreased the expression of proteins involved in the biosynthesis of amino acids. Carbon metabolism is crucial for lipid accumulation, which was significantly affected upon N deficiency. Upon N deprivation, the glycolysis pathway, which provides pyruvate, a key precursor for acetyl-CoA, was up-regulated, while the activity of the TCA cycle was retarded, thus, leading more carbon flux to fatty acid biosynthesis. Moreover, there is some evidence that the PPP was up-regulated and then this would increase the production of NADPH needed for fatty acid biosynthesis. Together, coordinated regulation of central carbon metabolism upon N limitation, provide more carbon flux to acetyl-CoA and NADPH for fatty acid biosynthesis. In addition, protein and nucleic acid metabolism was down-regulated upon N limitation, this shifts the cellular metabolism to lipid biosynthesis. Also, some protein involved in energy metabolism, signal transduction, molecular chaperone and redox homeostasis were up-regulated in response to stress condition of N deficiency.
Microorganism and cultivation
M. circinelloides WJ11 isolated in our laboratory from soil at Jiangnan University was used in this study [6, 49]. 100 µl spore suspension (approx. 107 spores/mL) of M. circinelloides WJ11 was cultivated in 150 mL K & R medium  held in 1 L flasks equipped with baffles for 24 h with shaking at 150 rpm and 30 °C, and then used at 10 % (v/v) to inoculate 2 L fermenters containing 1.5 L modified K & R medium (80 g glucose and 2 g diammonium tartrate per liter). Fermenters were controlled at 30 °C with stirring at 700 rpm and aeration at 0.5 v/v min−1. The pH was maintained at 6.0 by auto-addition of 4 M KOH or 2 M H2SO4.
Analysis of cell dry weight (CDW), culture supernatant and lipid accumulation
Biomass was periodically harvested by filtration using a dried and pre-weighed filter paper and a Buchner funnel under vacuum and washed three times with distilled water, frozen overnight at −80 °C and then freeze-dried. The weight of dry cells was determined gravimetrically. Glucose concentration in the culture medium was measured using a glucose oxidase kit according to the manufacturer’s instructions. Ammonium concentration in the culture filtrate was determined using the indophenol test . Total lipid was extracted and analyzed by the procedure reported in our previous work .
Protein extraction for two-dimensional electrophoresis (2-DE)
Cells taken from the balanced phase of growth (at 6 h), the fast lipid accumulation stage (at 24 h) and the stable stage of lipid accumulation (at 60 h) were filtered as above and washed with cold distilled water at 4 °C. The collected mycelia were flash-frozen in liquid N2 and stored at −80 °C. For each condition, protein extraction was performed according to a modified version of the method of Liu et al.  and Chen et al. . The frozen mycelia were ground in liquid N2, resuspended in cold 10 % (w/v) trichloracetic acid/acetone and allowed to precipitate at −20 °C for 1 h. The samples were centrifuged at 15,000g for 15 min at 4 °C and the supernatant was discarded. The pellets were resuspended in cold acetone, kept at −20 °C for 1 h and then centrifuged at 15,000g for 15 min at 4 °C. This procedure was repeated twice. The pellets were freeze-dried, suspended in approx. 10 mL extraction buffer [0.7 M sucrose, 0.1 M NaCl, 0.5 M Tris/HCl (pH 7.5), 50 mM EDTA and 0.2 % DTT], mixed with an equal volume of Tris/saturated phenol (1 g/100 mL) (pH 7.5) and then homogenized for 30 min at 4 °C. The homogenate was centrifuged at 5000 g for 10 min at 4 °C to collect the phenol phase and the phenol extraction was repeated three time. The combined phenol phases were mixed with five volumes of a precipitation buffer (0.1 M ammonium acetate in methanol). Precipitation was carried out at −20 °C for 1 h and the pellets were washed three times with cold methanol followed by rinsing three time with cold acetone. The protein pellets were dissolved in a lysis buffer [9 M urea, 4 % CHAPS, 1 % immobilized pH gradient buffer (pH 3–10) and 1 % DTT] and then centrifuged at 15,000g for 15 min. The supernatants were collected and protein concentrations were determined according to the Bradford method with BSA as a standard.
Before 2-DE, 1200 µg protein solution was mixed with rehydration buffer [9 M urea, 4 % CHAPS, 1 % immobilized pH gradient buffer (pH 3–10), 1 % DTT and 0.002 % Bromphenol Blue] and then loaded onto the IPG strip (pH 3–10, nonlinear, 24 cm). Strips were focused on a IPGhor Isoelectric Focusing System (GE Healthcare) at 20 °C with the following program: 50 V for 12 h (for rehydration), 100 V for 1 h, 200 V for 1 h, 500 V for 1 h, 1000 V for 1 h, gradient from 1000 to 10,000 V within 1 h and 10,000 V for 11 h. After that, strips were equilibrated at room temperature in two steps: 15 min in equilibration buffer [50 mM Tris/HCl (pH 8.8), 6 M urea, 30 % w/v glycerol and 2 % w/v SDS] with the addition of 1 % DTT followed by 15 min in equilibration buffer with the addition of 2.5 % (w/v) iodoacetamide. The equilibrated strips were transferred onto 12 % SDS-PAGE for the second dimension electrophoresis by using DALT-SIX SDS-PAGE Vertical System (GE Healthcare) at 15 °C with two steps: 100 V for 45 min and 200 V until the Bromophenol Blue reached the bottom of the gel. The gels were fixed with 10 % (w/v) trichloroacetic acid and stained with Coomassie Brilliant Blue G-250.
Nine 2-DE gels (three independent analytical replicate gels for each growth stage) were scanned at 300 dpi using Image Scanner LabScan (GE Healthcare). Spot detection, gel matching and group analysis of the gels were performed using PDQuest 8.0 software (Bio-Rad). Quantitative analyses were carried out after normalizing the quantities of spots in all gels in order to compensate for non-expression related variations, and quantity of each spot was normalized by total valid spot intensity. For each spot, the mean quantity was computed at every stage, and the spots showing a mean value that changed more than 1.5-fold or less than 0.67-fold (P < 0.05) in different stages were considered differentially expressed proteins.
Protein identification and database search
Protein spots with different expression levels were manually excised from the gels, washed with Millipore pure water for three times, destained three times with 100 mM NH4HCO3 in 30 % (v/v) acetonitrile, and then vacuum dried. Every protein spot was digested overnight with 50 ng trypsin (Promega) in 30 μl 25 mM NH4HCO3 containing 10 % (v/v) acetonitrile at 37 °C. The supernatants were transferred into another tube followed by vacuum dried. The dried peptides were dissolved in 0.1 % trifluoroacetic acid and mixed with an equal volume of 0.7 mg α-cyano-4-hydroxy-trans-cinnamic acid/mL in acetonitrile/trifluoroacetic acid (85:0.1 v/v), and then spotted on the sample target plate for analysis using UltrafleXtrem MALDI-TOF/TOF mass spectrometer (Bruker-Daltonics). Tryptic peptides were analyzed in the positive ion reflector mode, and spectra were calibrated using Bruker peptide calibration standard II (Bruker-Daltonics). At least 10 peptide fragments were selected to be analyzed in lift mode. Then spectra were processed by FlexAnalysis software and analyzed by BioTools software (Bruker-Daltonics). An in-house Mascot server (http://www.matrixscience.com) was used for database search and the following parameters were used in the search: NCBInr fungi database; trypsin as the digestion enzyme; monoisotopic peptide values; a maximum of one missed cleavage per peptide; fragment mass tolerance of 0.99 Da and peptide mass tolerance of 300 ppm, together with the acceptance of cysteine carbamidomethylation (fixed modifications) and methionine oxidation (variable modifications). For a positive identification, a score calculated by the Mowse scoring algorithm in MASCOT was considered as significant (P < 0.05).
Quantitative RT-PCR analysis was performed to quantify the transcriptional levels of genes. Total RNA was extracted from samples taken at the three growth stages with TRIzol and then converted to cDNA using a PrimeScrip 1st strand cDNA Synthesis Kit (Takara) according to the manufacturer’s instructions. The quantitative RT-PCR reaction was performed with CFX Connect Real-Time System (Bio-Rad) and iTaq Universal SYBR Green PCR Supermix (Bio-Rad) was used to identify mRNA level. The primer pairs are listed in Additional file 2: Table S2 using the 18S rRNA of M. circinelloides WJ11 as an internal control in PCR amplification.
The mean values and the standard error of the mean were calculated from the data obtained from three biological replicates. A statistical analysis of the obtained data was carried out using SPSS 16.0 for Windows (SPSS Inc., Chicago, IL). One way analysis of variance (ANOVA) with Tukey’s test was conducted on the data, and P < 0.05 was considered significantly different.
cell dry weight
polyunsaturated fatty acids
- ACL, ATP:
total fatty acids
branched-chain amino acids
reactive oxygen species
pentose phosphate pathway
thiazole biosynthetic enzyme
heat shock proteins
XT carried out the experiments and drafted the manuscript. XZ and LZ participated in protein identification and database search. HC, YQC, WC, YS, and CR participated in the experimental design and reviewed the manuscript. YS and CR conceived the study and reviewed the final manuscript. All authors read and approved the final manuscript.
This work was supported by the National Natural Science Foundation of China (31271812, 31530056, 21276108), the National High Technology Research and Development Program of China (2012AA022105C), the Program for Changjiang Scholars and Innovative Research Team in University (IRT1249), the Program for New Century Excellent Talents (NCET-13-0831), and the Strategic Merieux Research Grant.
The authors declare that they have no competing interests.
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.
- Thorpe R, Ratledge C. Fatty acid distribution in triglycerides of yeasts grown on glucose or n-alkanes. J Gen Microbiol. 1972;72:151–63.View ArticleGoogle Scholar
- Ratledge C. The role of malic enzyme as the provider of NADPH in oleaginous microorganisms: a reappraisal and unsolved problems. Biotech Lett. 2014;36:1557–68.View ArticleGoogle Scholar
- Ratledge C, Wynn JP. The biochemistry and molecular biology of lipid accumulation in oleaginous microorganisms. Adv Appl Microbiol. 2002;51:1–51.View ArticleGoogle Scholar
- Ratledge C. Microbial production of gamma-linolenic acid. In: Akoh CC, editor. Handbook of functional lipids. Boca Raton: CRC Press; 2005. p. 19–45.View ArticleGoogle Scholar
- Botham PA, Ratledge C. A biochemical explanation for lipid accumulation in Candida 107 and other oleaginous micro-organisms. J Gen Microbiol. 1979;114:361–75.View ArticleGoogle Scholar
- Tang X, Chen H, Chen YQ, Chen W, Garre V, Song Y, Ratledge C. Comparison of biochemical activities between high and low lipid-producing strains of Mucor circinelloides: an explanation for the high oleaginicity of strain WJ11. PLoS ONE. 2015;10:e0128396.View ArticleGoogle Scholar
- Boulton CA, Ratledge C. Correlation of lipid accumulation in yeasts with possession of ATP: citrate lyase. J Gen Microbiol. 1981;127:169–76.Google Scholar
- Tamano K, Bruno KS, Karagiosis SA, Culley DE, Deng S, Collett JR, Umemura M, Koike H, Baker SE, Machida M. Increased production of fatty acids and triglycerides in Aspergillus oryzae by enhancing expressions of fatty acid synthesis-related genes. Appl Microbiol Biotechnol. 2013;97:269–81.View ArticleGoogle Scholar
- Zhang H, Zhang L, Chen H, Chen YQ, Chen W, Song Y, Ratledge C. Enhanced lipid accumulation in the yeast Yarrowia lipolytica by over-expression of ATP: citrate lyase from Mus musculus. J Biotechnol. 2014;192:78–84.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Wynn JP, Ratledge C. Malic enzyme is a major source of NADPH for lipid accumulation by Aspergillus nidulans. Microbiology. 1997;143:253–7.View ArticleGoogle Scholar
- Zhu Z, Zhang S, Liu H, Shen H, Lin X, Yang F, Zhou YJ, Jin G, Ye M, Zou H. A multi-omic map of the lipid-producing yeast Rhodosporidium toruloides. Nat Commun. 2012;3:1112.View ArticleGoogle Scholar
- Liu Z, Gao Y, Chen J, Imanaka T, Bao J, Hua Q. Analysis of metabolic fluxes for better understanding of mechanisms related to lipid accumulation in oleaginous yeast Trichosporon cutaneum. Bioresour Technol. 2013;130:144–51.View ArticleGoogle Scholar
- Paula Alonso A, Dale VL, Shachar Hill Y. Understanding fatty acid synthesis in developing maize embryos using metabolic flux analysis. Metab Eng. 2010;12:488–97.View ArticleGoogle Scholar
- Xiong W, Liu L, Wu C, Yang C, Wu Q. 13C-tracer and gas chromatography-mass spectrometry analyses reveal metabolic flux distribution in the oleaginous microalga Chlorella protothecoides. Plant Physiol. 2010;154:1001–11.View ArticleGoogle Scholar
- Wasylenko TM, Ahn WS, Stephanopoulos G. The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica. Metab Eng. 2015;30:27–39.View ArticleGoogle Scholar
- Zhao L, Zhang H, Wang L, Chen H, Chen YQ, Chen W, Song Y. 13 C-metabolic flux analysis of lipid accumulation in the oleaginous fungus Mucor circinelloides. Bioresour Technol. 2015;197:23–9.View ArticleGoogle Scholar
- Yang ZK, Ma YH, Zheng JW, Yang WD, Liu JS, Li HY. Proteomics to reveal metabolic network shifts towards lipid accumulation following nitrogen deprivation in the diatom Phaeodactylum tricornutum. J Appl Phycol. 2014;26:73–82.View ArticleGoogle Scholar
- Garnier M, Carrier G, Rogniaux H, Nicolau E, Bougaran G, Saint-Jean B, Cadoret J-P. Comparative proteomics reveals proteins impacted by nitrogen deprivation in wild-type and high lipid-accumulating mutant strains of Tisochrysis lutea. J Proteomics. 2014;105:107–20.View ArticleGoogle Scholar
- Song P, Li L, Liu J. Proteomic analysis in nitrogen-deprived Isochrysis galbana during lipid accumulation. PLoS ONE. 2013;8:e82188.View ArticleGoogle Scholar
- Liu H, Zhao X, Wang F, Li Y, Jiang X, Ye M, Zhao ZK, Zou H. Comparative proteomic analysis of Rhodosporidium toruloides during lipid accumulation. Yeast. 2009;26:553–66.View ArticleGoogle Scholar
- Chen H, Hao G, Wang L, Wang H, Gu Z, Liu L, Zhang H, Chen W, Chen YQ. Identification of a critical determinant that enables efficient fatty acid synthesis in oleaginous fungi. Sci Rep. 2015;5:11247.View ArticleGoogle Scholar
- Xia C, Zhang J, Zhang W, Hu B. A new cultivation method for microbial oil production: cell pelletization and lipid accumulation by Mucor circinelloides. Biotechnol Biofuels. 2011;4:1–10.View ArticleGoogle Scholar
- Tang X, Zhang H, Chen H, Chen YQ, Chen W, Song Y. Effects of 20 standard amino acids on the growth, total fatty acids production, and γ-linolenic acid yield in Mucor circinelloides. Curr Microbiol. 2014;69:899–908.View ArticleGoogle Scholar
- Stacey G, Van Baalen C, Tabita FR. Nitrogen and ammonia assimilation in the cyanobacteria: regulation of glutamine synthetase. Arch Biochem Biophys. 1979;194:457–67.View ArticleGoogle Scholar
- Magasanik B, Kaiser CA. Nitrogen regulation in Saccharomyces cerevisiae. Gene. 2002;290:1–18.View ArticleGoogle Scholar
- Chipman D, Ze Barak, Schloss JV. Biosynthesis of 2-aceto-2-hydroxy acids: acetolactate synthases and acetohydroxyacid synthases. BBA Protein Struct Mol Enzymol. 1998;1385:401–19.View ArticleGoogle Scholar
- Myers JW. Dihydroxy acid dehydrase: an enzyme involved in the biosynthesis of isoleucine and valine. J Biol Chem. 1961;236:1414–8.Google Scholar
- Chunduru SK, Mrachko GT, Calvo K. Mechanism of ketol acid reductoisomerase. Steady-state analysis and metal ion requirement. Biochemistry. 1989;28:486–93.View ArticleGoogle Scholar
- Ichihara A, Koyama E. Transaminase of branched chain amino acids. J Biochem. 1966;59:160–9.Google Scholar
- Huber TA, Streeter JG. Asparagine biosynthesis in soybean nodules. Plant Physiol. 1984;74:605–10.View ArticleGoogle Scholar
- Fujioka M, Nakatani Y. Saccharopine dehydrogenase. Eur J Biochem. 1972;25:301–7.View ArticleGoogle Scholar
- Miyazaki J, Kobashi N, Nishiyama M, Yamane H. Characterization of homoisocitrate dehydrogenase involved in lysine biosynthesis of an extremely thermophilic bacterium, Thermus thermophilus HB27, and evolutionary implication of β-decarboxylating dehydrogenase. J Biol Chem. 2003;278:1864–71.View ArticleGoogle Scholar
- Kikuchi G. The glycine cleavage system: composition, reaction mechanism, and physiological significance. Mol Cell Biochem. 1973;1:169–87.View ArticleGoogle Scholar
- Lu SC. S-adenosylmethionine. Int J Biochem cell B. 2000;32:391–5.View ArticleGoogle Scholar
- Shin R, Berg RH, Schachtman DP. Reactive oxygen species and root hairs in Arabidopsis root response to nitrogen, phosphorus and potassium deficiency. Plant Cell Physiol. 2005;46:1350–7.View ArticleGoogle Scholar
- Guy C, Carter J. Effect of low temperature on the glutathione status of plant cells. Plant Cold Hardiness Freez Stress. 1982;2:169–79.View ArticleGoogle Scholar
- Kendrick A, Ratledge C. Desaturation of polyunsaturated fatty acids in Mucor circinelloides and the involvement of a novel membrane-bound malic enzyme. Eur J Biochem. 1992;209:667–73.View ArticleGoogle Scholar
- Marsh JJ, Lebherz HG. Fructose-bisphosphate aldolases: an evolutionary history. Trends Biochem Sci. 1992;17:110–3.View ArticleGoogle Scholar
- Song Y, Wynn JP, Li Y, Grantham D, Ratledge C. A pre-genetic study of the isoforms of malic enzyme associated with lipid accumulation in Mucor circinelloides. Microbiology. 2001;147:1507–15.View ArticleGoogle Scholar
- Zhang Y, Adams IP, Ratledge C. Malic enzyme: the controlling activity for lipid production? Overexpression of malic enzyme in Mucor circinelloides leads to a 2.5-fold increase in lipid accumulation. Microbiology. 2007;153:2013–25.View ArticleGoogle Scholar
- Noctor G, Mhamdi A, Chaouch S, Han Y, Neukermans J, Marquez-Garcia B, Queval G, Foyer CH. Glutathione in plants: an integrated overview. Plant, Cell Environ. 2012;35:454–84.View ArticleGoogle Scholar
- Machado C, Oliveira RCL, Boiteux S, Praekelt UM, Meacock PA, Menck CFM. Thi1, a thiamine biosynthetic gene in Arabidopsis thaliana, complements bacterial defects in DNA repair. Plant Mol Biol. 1996;31:585–93.View ArticleGoogle Scholar
- Machado CR, Praekelt UM, Oliveira RCL, Barbosa ACC, Byrne KL, Meacock PA, Menck CF. Dual role for the yeast THI4 gene in thiamine biosynthesis and DNA damage tolerance. J Mol Biol. 1997;273:114–21.View ArticleGoogle Scholar
- Wynn JP, Hamid AA, Li Y, Ratledge C. Biochemical events leading to the diversion of carbon into storage lipids in the oleaginous fungi Mucor circinelloides and Mortierella alpina. Microbiology. 2001;147:2857–64.View ArticleGoogle Scholar
- Hermeking H, Benzinger A. 14-3-3 proteins in cell cycle regulation. Semin Cancer Biol. 2006;16:183–92.View ArticleGoogle Scholar
- Choi YE, Kwon KW, Lee JC, Woo SY. Expression of the rice cytoplasmic cysteine synthase gene in tobacco reduces ozone-induced damage. Plant Biotechnol Rep. 2007;1:93–100.View ArticleGoogle Scholar
- Bell A, Monaghan P, Page AP. Peptidyl-prolyl cis–trans isomerases (immunophilins) and their roles in parasite biochemistry, host-parasite interaction and antiparasitic drug action. Int J Parasitol. 2006;36:261–76.View ArticleGoogle Scholar
- Tang X, Zhao L, Chen H, Chen YQ, Chen W, Song Y, Ratledge C. Complete genome sequence of a high lipid-producing strain of Mucor circinelloides WJ11 and comparative genome analysis with a low lipid-producing strain CBS 277.49. PLoS ONE. 2015;10:e0137543.View ArticleGoogle Scholar
- Chaney AL, Marbach EP. Modified reagents for determination of urea and ammonia. Clin Chem. 1962;8:130–2.Google Scholar
- Liu X, Wu H, Ji C, Wei L, Zhao J, Yu J. An integrated proteomic and metabolomic study on the chronic effects of mercury in Suaeda salsa under an environmentally relevant salinity. PLoS ONE. 2013;8:e64041.View ArticleGoogle Scholar
- Chen Y, Pang Q, Dai S, Wang Y, Chen S, Yan X. Proteomic identification of differentially expressed proteins in Arabidopsis in response to methyl jasmonate. J Plant Physiol. 2011;168:995–1008.View ArticleGoogle Scholar