Next Article in Journal
Identification of More Feasible MicroRNA–mRNA Interactions within Multiple Cancers Using Principal Component Analysis Based Unsupervised Feature Extraction
Previous Article in Journal
Exchange Transfusion in the Treatment of Neonatal Septic Shock: A Ten-Year Experience in a Neonatal Intensive Care Unit
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 5
10.3390/ijms17050694
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Metabolic Regulation of the Snow Alga Chlamydomonas nivalis in Response to Nitrate or Phosphate Deprivation by a Metabolome Profile Analysis

by
Na Lu
1,†,
Jun-Hui Chen
1,†,
Dong Wei
1,*,
Feng Chen
1,2 and
Gu Chen
1
1
School of Food Science and Engineering, South China University of Technology, Guangzhou 510641, China
2
Institute for Food and Bioresource Engineering, College of Engineering, Peking University, Beijing 100871, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the work and should be regarded as co-first authors.
Int. J. Mol. Sci. 2016, 17(5), 694; https://doi.org/10.3390/ijms17050694
Submission received: 1 April 2016 / Revised: 29 April 2016 / Accepted: 29 April 2016 / Published: 10 May 2016
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
In the present work, Chlamydomonas nivalis, a model species of snow algae, was used to illustrate the metabolic regulation mechanism of microalgae under nutrient deprivation stress. The seed culture was inoculated into the medium without nitrate or phosphate to reveal the cell responses by a metabolome profile analysis using gas chromatography time-of-flight mass spectrometry (GC/TOF-MS). One hundred and seventy-one of the identified metabolites clustered into five groups by the orthogonal partial least squares discriminant analysis (OPLS-DA) model. Among them, thirty of the metabolites in the nitrate-deprived group and thirty-nine of the metabolites in the phosphate-deprived group were selected and identified as “responding biomarkers” by this metabolomic approach. A significant change in the abundance of biomarkers indicated that the enhanced biosynthesis of carbohydrates and fatty acids coupled with the decreased biosynthesis of amino acids, N-compounds and organic acids in all the stress groups. The up- or down-regulation of these biomarkers in the metabolic network provides new insights into the global metabolic regulation and internal relationships within amino acid and fatty acid synthesis, glycolysis, the tricarboxylic acid cycle (TCA) and the Calvin cycle in the snow alga under nitrate or phosphate deprivation stress.

Graphical Abstract

1. Introduction

Chlamydomonas nivalis, a typical unicellular green microalga, is distributed worldwide in the snowfield of polar regions and similar extreme environments, and it has attracted attention recently due to its application as a model species of snow algae for illustrating the cellular response mechanism to stress conditions [1]. Nutrient deprivation is the most important stress factor that significantly impacts algal survival in harsh habitats, and it induces various cell responses through different mechanisms in the microalgae to adapt to the stress [2]. Of all the nutrients, nitrogen and phosphorus are two macronutrients with broad effects on algal growth and survival through complex biochemical processes, involving metabolic energy generation and transformation [3,4], membrane stability and signaling [5] and photosynthetic metabolism [6]. When a specific nutrient is deprived, the corresponding responses of the algal cells tend to regulate the metabolic pathways to compensate for the nutrient limitation [7]. For example, the rapid removal of nitrogen from the medium activates the assimilation of inorganic nitrogen, making nitrogen available through the production and re-use of NH4+ [8]. Phosphorus limitation leads to the up-regulation of protein digestion associated with lipid accumulation, while the down-regulation of protein biosynthesis is associated with amino acid metabolism in the marine diatom Phaeodactylum tricornutum [9]. The responses to sulfur starvation are quite similar to those of phosphorus starvation. The uptake and assimilations of cellular sulfate-containing metabolites are induced immediately, which directly affect the metabolism of vitamins and polyamines in Arabidopsis [10]. Most previous reports indicate that the responses of microalgae, corresponding to nutrient limitation as stress conditions, are achieved through various molecular regulations in different metabolic pathways.
The metabolome profile reflects the cellular metabolic stage at a specific physiological phase in response to environmental changes [11]. A metabolome profile analysis is a valid approach towards a quantitative description of the regulated and controllable process under stress conditions [12], and therefore, it provides unique perspectives of the fundamental regulations in algal cells. Gas chromatography time-of-flight mass spectrometry (GC/TOF-MS) is an effective method for metabolome profile analyses obtaining a global vision of the small molecules of metabolites (molecular weight <1000 dalton) due to its convenience and high accuracy in identifying and quantifying more than 100 metabolites in biological samples by a notably fast scan from a complicated matrix [13]. These small molecules, which include carbohydrates, organic acids, alcohols, aldehydes and amines, provide a comprehensive overview in the primary metabolic pathways [14]. The detailed procedures for metabolite profiling based on GC/TOF-MS have been established using Chlamydomonas reinhardtii as the model species [15,16,17,18], which are crucial for investigating the metabolism of microalgae. Based on these approaches and methods, the metabolic pathways of Chlamydomonas reinhardtii, changing accordingly under the deficiencies of nutrients (nitrogen, phosphorus, sulfur and iron), were investigated recently using GC/TOF-MS [19,20]. However, due to complications in establishing stable technical procedures of the sampling and extraction, only a few algal studies using GC-MS have been reported, which are worthy of further investigation to explore the algal metabolisms by GC/TOF-MS [15,21,22].
The Kyoto encyclopedia of genes and genomes (KEGG) [23] is an online tool to reconstruct the metabolic network from a similar organism and to manually draw pathway maps based on the representing knowledge about the molecular reactions, the regulation of metabolic pathways and the genomic information in various organisms. Pandhal et al. [24] investigated the metabolism of the cyanobacterium Synechocystis sp. PCC6803 (Synechocystis 6803) and its adaption to the high salinity stress by connecting the metabolomic and proteomic data to the KEGG pathway. The nitrogen-starved metabolism of the diatom was investigated by linking the change in the protein biosynthesis to the carbon metabolisms in the KEGG pathway [25].
In our previous studies, the responses of C. nivalis to salt stress and nutrient starvation (nitrate or phosphate deprivation stress) were investigated on a lipidomic level using UPLC/QTOF-MS (ultra-performance liquid chromatography/quadrupole time-of-flight mass spectrometry) coupled with OPLS-DA (the orthogonal partial least squares discriminant analysis), and more than twenty kinds of lipid/fatty acid biomarkers were selected and identified [5,26,27,28]. However, the biological functions of these biomarkers, involving membrane stability, signaling and photosynthesis, are under speculation according to the literature. It is difficult to link the changes of lipid biomarkers in the membrane system to the intercellular metabolic network. Therefore, the aim of this work is to elucidate the metabolic changes in certain stress conditions and to reveal the possible regulation mechanisms under these stress conditions by coupling GC/TOF-MS with a multivariate statistical analysis.
In the present study, GC/TOF-MS was first applied to detect the whole metabolome in the different stress groups of C. nivalis cells, and a multivariate statistical analysis (OPLS-DA) was employed to recognize patterns of the metabolome data from the groups. Then, the potential “responding biomarkers” were selected and identified. Significant changes in the abundance of the responding biomarkers were discussed subsequently, and finally they were visualized in the map of the KEGG pathway to reveal the global regulation of the metabolic network under these stress conditions.

2. Results and Discussion

2.1. Identification and Overview of the Metabolites

By rapid sampling, extraction and standard analysis using GC/TOF-MS, 670 metabolite peaks were detected and imported to the BinBase database, but only 171 metabolites were identified with a chemical structure and name. Overall, they belonged to several classes, including carbohydrates, organic acids, polyols, amino acids, N-compounds (nitrogen-containing compounds except amino acids), phosphates, fatty acids and other compounds.
Figure 1 shows the overview of the identified metabolites in classes and their abundance in the stress groups and the control. The composition in the classes looks similar in all the groups, but the abundance of each class was significantly different compared to the control. Similar changes in the trends of their abundance in all of the classes were found in both of the stress groups; except that the class of phosphates increased significantly (p < 0.01) in the nitrate-deprived group but decreased significantly (p < 0.01) in the phosphate-deprived group.
Furthermore, comparing the abundance with that of the control, the carbohydrates significantly increased by 11.78% and 7.62% in the nitrate- and phosphate-deprived groups (p < 0.01), respectively. This could be due to the decreased catabolism and reutilization of the stored carbon compounds for the biosynthesis under both stress conditions [29]. Amino acids decreased by 12.47% and 2.09% (p < 0.01) in both of the stress conditions, and N-compounds also decreased significantly by 14.06% and 10.35%, respectively (p < 0.01). Protein biosynthesis is slowed down by low nitrate and phosphate concentrations [29]. A nitrogen limitation reduces the protein synthesis, mainly due to the limited translation of mRNA without the efficient supply of amino acids, while a phosphorus limitation causes a short supply of adenosine triphosphate (ATP), which leads to the accumulation of carbohydrates in preference to proteins due to the fact that the enzymes for the protein biosynthesis have a low affinity for ATP [30]. Additionally, organic acids, including the central intermediates from glycolysis, the citric acid cycle and pyruvate metabolism, decreased by 9.2% and 4.8%, respectively (p < 0.01). However, fatty acids increased 20.9% and 7.65%, respectively (p < 0.01). Excess carbon resources are preferred to biosynthesize fatty acids during nutrient limitation, which can explain the increased triglyceride in response to the nitrogen or phosphate limitation [25]. Therefore, following evidence of a global metabolic rebalance in the stress groups, the investigation in this work was extended toward a metabolome profile analysis.

2.2. Discrimination of Metabolomic Data by OPLS-DA (the Orthogonal Partial Least Squares Discriminant Analysis)

To assess the differentiation between the stress groups and the control on a metabolomic level, the OPLS-DA model was employed to discriminate the origins of the metabolite profiles [17]. The score scatter plot in Figure 2 was generated from the peak area of the identified metabolites in the different groups. Five experimental clusters with time points were clearly separated by the OPLS-DA model, where each spot represents a sample with the parameters of R2X, R2Y, and Q2 (0.944, 0.981, and 0.971, respectively). The samples from the different experimental groups were analyzed repeatedly to validate the original model, which showed that this procedure demonstrated great repeatability from the stress treatments to the data analysis. In the “y-scrambling” validation that was applied for this OPLS-DA model, random permutation tests were performed, and the parameter of Q2 (−0.364 < 0) indicated that the valid model had high predictability without overfitting [31].
In Figure 2, the samples that were stressed for 6 and 12 h clustered together in the nutrient-deprived groups and separated with another cluster (stressed for 24, 48, and 72 h), indicating that 6 and 24 h might be the turning point of this alga for the essential metabolic regulation responding to the nutrient deprivation stress. All the clusters in the phosphate-deprived group located closer to the control, and thus, their metabolite profiles are more similar to the control than the nitrate-deprived group. This could be due to the fact that the algae were not yet fully P-deprived at the beginning of the P-stressed treatment. The organisms often store substantial amounts of excess P as polyphosphates to maintain the phosphorus use efficiency during P deprivation [32]. These results also indicated that the metabolism in this alga was affected more effectively by nitrate than phosphate, and the OPLS-DA model was perfect for obtaining metabolic readouts in investigating the influence of nutrients in the culture of C. nivalis.

2.3. Selection of the Responding Biomarkers

To illustrate the pattern of metabolic regulation that responds to nutrient-deprived stress in C. nivalis, key contributors among the metabolites that clustered needed to be identified as responding biomarkers. For this purpose, the related plots of the OPLS-DA models shown in Figure 3 and Figure 4 were generated from both of the stress groups vs. the control, respectively. The performance of R2X, R2Y, and Q2 in the OPLS-DA model generated from the nitrated-deprived group vs. the control (Figure 3a) was 0.953, 0.999, 0.998, respectively. The intercept value of Q2 in the permutation test was −0.34. All the results showed that the established model is robust and not overfitting. The procedure for the selection of the responding biomarkers was the same as described in our previous work [28]. Briefly, the potential biomarkers had to meet all of the three criteria (i.e., VIP (variable influence on the projection) > 1, ｜p (corr(correlation coefficient))｜> 0.6 and the span of CIJKJK(jack-knifed confidence interval) excluding 0) for further investigation. Overall, thirty metabolites were selected as “responding biomarkers” among which the largest VIP value was 5.48 from the nitrate-deprived group vs. the control. These included one carbohydrate, five N-compounds, fourteen amino acids, six organic acids and four fatty acids (shown in Table 1).
The performance of R2X, R2Y, and Q2 in the OPLS-DA model generated from the phosphate-deprived group vs. the control (Figure 4a) was 0.877, 0.987, 0.971, respectively, and the intercept value of Q2 in the permutation test was −0.348, indicating that the established model is also robust and not overfitting. According to the same procedure above, thirty-nine of the identified metabolites were selected as “responding biomarkers” among which the largest VIP value is 3.42 from the phosphate-deprived group vs. the control. These included five carbohydrates, four N-compounds, twelve amino acids, ten organic acids and eight fatty acids (shown in Table 2).
Based on the same criteria, the results showed that less responding biomarkers were selected from the nitrate-deprived group, but the largest VIP value was higher than that from the phosphate-deprived group, indicating that nitrate has a greater effect than the phosphate on the metabolomic profiling. Especially, the unique biomarkers discovered only from the nitrate-deprived group vs. the control were tryptophan, lysine, isoleucine, glutamine and citric acid. Most of these are amino acids (in Table 1), suggesting that nitrate deprivation has a greater influence on amino acid metabolism. While, the unique biomarkers discovered only from the phosphate-deprived group vs. the control were trehalose, sucrose, cellobiose, threitol, tyrosine, threonine, threonic acid, gluconic acid, erythronic acid lactone, benzoic acid, 3-phosphoglycerate, pyrophosphate, inositol-4-monophosphate, fructose-6-phosphate and adenosine-5-phosphate. Most of these are carbohydrates and organic acids related to central carbon metabolism (Table 2).

2.4. Abundance of the Responding Biomarkers in the Different Groups

The abundance of the responding biomarkers in the classes in the different groups is shown in Figure 5. Although the abundance changes in all the classes were significantly different compared to the control, the changes in the biomarker classes have the same trends as that of the whole identified metabolite classes as shown in Figure 1. Therefore, the cellular response and adaptation to both nutrient deprivations could be represented by the changes in the abundance in the responding biomarkers during the certain stress.
As shown in Figure 5, the abundance of the carbohydrate biomarkers significantly increased by 11.46% (p < 0.01, Figure 5a), which was mainly trisaccharide (Table 1), and by 31.34% (p < 0.05, Figure 5b), which was mainly trisaccharide, trehalose, sucrose and cellobiose (Table 2), in the N and P stressed groups. However, a significant decrease in the abundance of the N-compound (p < 0.01) and amino acid (p < 0.05) biomarkers was found in the two stress groups. Moreover, the decrease in the abundance of the amino acid biomarkers in the N deprivation group (p < 0.01) was much higher than in the P deprivation group (p < 0.05), indicating that N deprivation affects the metabolism of amino acids more effectively. The abundance of the fatty acid biomarkers significantly increased by 46.12% (p < 0.01, Figure 5a) and 11.45% (p < 0.01, Figure 5b) in the two stress groups, which is in accordance with that of the carbohydrate biomarkers. Linolenic acid (18:3), stearic acid (C18:0), palmitic acid (16:0), and pentadecanoic acid (15:0) were the fatty acid biomarkers with the largest variation in both of the stress groups and also increased in all the nitrate-deprived groups (Table 1), but the linolenic acid decrease was observed only in the phosphate-deprived group (Table 2).
It was proposed that changes in carbon and nitrogen metabolism to regulate their balance for storing excess carbon in molecular pools are crucial under conditions of limited or no nitrogen. When higher amounts of polysaccharides are released, this was caused by a switch in metabolism from protein to carbohydrate [29]. However, the metabolic pathways and the enzymes activated by the nutrient deprivations were quite different in various microalgae. For example, the overexpression of acetyl-CoA carboxylase, which increases oil production in green algae, does not increase oil accumulation in diatom [33]. Although nitrogen affects the biosynthesis of fatty acids and lipids, the direct explanation for the effects in algae is not clearly understood at this time. The assimilation of carbon under nitrogen starvation is tuned to accumulate carbon skeletons and reduce equivalents for the fatty acid biosynthesis instead of nitrogenous compounds, such as amino acids [25]. In the freshwater green alga Chlamydomonas reinhardtii, excess carbon is directed to the biosynthesis of fatty acids and lipids under nitrogen-deprived conditions [34] and the up-regulation of the expression of the related genes [35].

2.4.1. Responding Biomarkers in the Nitrate-Deprived Group

In the nitrate-deprived group, five N-compounds and fourteen amino acids were selected as responding biomarkers, and the change rates in abundance are shown in Table 1. Among these, eleven amino acid biomarkers (valine, serine, proline, phenylalanine, lysine, leucine, glutamine, glutamic acids, aspartic acid, asparagine and alanine) decreased from 1.03% to 58.9%, contributing to the most declines of the total amino acids. Asn, Asp, Pro and Gln are well known to be involved in the biosynthesis of nitrogen compounds and the metabolism of nitrogen transportation. Reducing the supply of these amino acids suddenly led to the reduction of the de novo biosynthesis of several amino acids. Conversely, three amino acid biomarkers (tryptophan, isoleucine and histidine) increased by 85.84%, 75.56% and 86.94%, respectively. This might be due to the fact that the enzymes that enhance the nitrogen uptake were induced in response to N deprivation, while the other protein biosynthesis was reduced due to the decreased availability of the nitrogen sources [35]. These results indicate that C. nivalis is incapable of maintaining the homeostasis of amino acid metabolism in vivo under this stress due to changes in the biosynthesis of several amino acids and amino group mediators.
The ratio of glutamic acid to glutamine (Glu/Gln) decreased 19.69% (from 0.52 to 0.42), while the proportion of glutamine and α-ketoglutaric acid increased 68.94% (from 0.96 to 1.63) in the nitrate-deprived group (Table 1). The proportion of glutamine/α-ketoglutaric acid is crucial to affect the bioactivities of glutamine synthetase and glutamate synthase under nitrogen stress, and it remains stable by regulating the capability of utilizing alternative nitrogen sources in several metabolisms [36]. Specially, the increase of this ratio (glutamine/α-ketoglutaric acid) leads to a decreased capability of nitrogen transport and metabolism and protein turnover (N-recycling) and an increased mobilization of C-reserves. In addition, citric acid, α-ketoglutaric acid and the downstream metabolites (fumaric acid and malic acid) in the tricarboxylic acid cycle (TCA) decreased 33.14%, 55.11%, 34.74% and 11.30%, respectively in the nitrate-deprived group (Table 1). These decrease indicated that an inhibition effect occurred in the TCA cycle in the N-stress condition [37]. However, a small increase of 8.36% for succinic acid was found in our experiments as shown in Table 1. As discussed before, reduced nitrate assimilation leads to the accumulation of reducing equivalents (NADH), and thus, less malic acid is formed from succinic acid.

2.4.2. Responding Biomarkers in the Phosphate-Deprived Group

The declines of 43.47% for intracellular adenosine and 4.98% for adenosine-5-phosphate (AMP) were the distinct metabolic changes observed in response to phosphate deprivation (Table 2). The same trend was found in Selenastrum minutum [38]. Because the activity of certain glycolytic enzymes was affected by the supplementation of ATP, ADP (adenosine diphosphate) or Pi (inorganic phosphate), the reduced biosynthesis of these metabolites had a great impact on the function of glycolysis in the algal cells. The reduced respiration rates and the consequential reduction in the Calvin cycle and light utilization were also observed in Selenastrum sp. due to the cellular metabolite changes under a phosphorus limitation [39].
Inositol-4-monophosphate decreased 37.26% due to the Pi limitation (Table 2). Inositol-4-monophosphate is a phosphorylated inositol and plays key roles in cellular metabolic processes, especially those involved in the signaling pathway of phosphophatidylinositol. Thus, it could partly explain the decreased cell growth of C. nivalis in the phosphate-deprived stress condition [5]. Interestingly, pyrophosphate increased 41.9% in this stress group. In general, pyrophosphate is presumably used to regulate the concentration of Pi in the cytoplasm, and it may play an important role in the storage of internal calcium and calcium mediated signal transduction [40]. Polyphosphate is mainly biosynthesized from the assimilation of the inorganic phosphate in the cells of Chlamydomonas under phosphate stress [41]. Park et al. [8] reported that high-phosphorus treatments induced a higher cellular consumption of phosphorus, which led to a significant increase of pyrophosphate in algae, but the treatments did not significantly affect the biosynthesis of polyphosphates, usually as the subsequent intracellular storage of phosphorus. The large accumulation of polyphosphates and pyrophosphate is also toxic to the microalgal cells because it depresses the key metabolite reactions related to proteins, nucleic acids and lipids [42]. Additionally, aspartic acid increased 42.85% in the phosphate-deprived group. Aspartic acid is an inhibitor of pyruvate carboxylase, which catalyzes the carboxylation of pyruvate to oxaloacetate in order to supply the intermediates of the citrate cycle [43]. Therefore, we speculated that the aspartic acid increase led to the inhibition of the activities of pyruvate carboxylase, which further demonstrated an inhibitory effect on the TCA cycle. These data are consistent with a previous study that aspartate inhibits pyruvate carboxylase competitively with respect to acetyl-CoA and decreases the levels of the TCA cycle intermediates [44]. Based on the above results, phosphorus deprivation increases the intermediates of the TCA cycle and reduces the substrates in the Calvin cycle and respiration.

2.5. Metabolic Pathway and Regulation Analysis under Stresses

To visualize the global regulation of the metabolic network in which the responding biomarkers are involved, the IDs (Identity) of the biomarkers in the Binbase database were matched with the IDs in KEGG and were then mapped on the metabolic pathways based on the Chlamydomonas reinhardtii pathway network in the KEGG database [45]. Almost all the responding biomarkers (threonic acid, threitol, erythronic acid-1,4-lactone and pentadecanoic acid were not found in the KEGG pathways) in the major pathways and their regulation referred as the change rates in abundance under stresses are shown in Figure 6 and Figure 7. Seven central pathways in microalgae are included, such as photorespiration, the Calvin cycle, glycolysis, amino acid synthesis, the TCA cycle, fatty acid metabolism, succinate production and glutamate synthesis.
Two major amino acids (serine and glutamic acid) and α-ketoglutaric acid were identified as biomarkers participating in photorespiration (Figure 6 and Figure 7). In the nitrate-deprived group, a decline of their abundance indicated a decrease in photorespiration (Table 1); while, in the phosphate-deprived group, glutamic acid increased 10.29% and 3-phosphoglycerate increased 297.76% (Table 2), which could be synthesized through the enhanced photorespiration, the Calvin cycle and glycolysis for carbon cycling.
Three fatty acid biomarkers (linolenic acid, stearic acid and palmitic acid) participate in fatty acid metabolism. Their abundance changes showed the same trends in the two stress groups (Figure 6 and Figure 7) with increased contents of palmitic acid and stearic acid and a decreased content of linolenic acid, leading to the reduced abundance of unsaturated lipids [17] to enhance the stress tolerance of nutrient deprivation by changing the membrane fluidity and permeability [27]. The results are in accordance with that fact that the fatty acid biomarkers increased significantly in the two stress groups as described in Figure 5.
The metabolisms of several amino acid biomarkers are relatively complicated under different stress conditions. In the nitrate-deprived group (Figure 6), all the amino acid biomarkers mapped on the network were decreased except for isoleucine [34]. In the phosphate-deprived group (Figure 7), some of the amino acid biomarkers (aspartic acid, threonine, glutamic acid and tyrosine) increased, while others decreased, and thus, the regulatory mechanisms were not obvious. Five organic acids (succinic acid, citric acid, fumaric acid, malic acid, and α-ketoglutaric acid) were detected as responding biomarkers in the TCA cycle. In the nitrate-deprived group (Figure 6), succinic acid increased slightly at 8.36%, but this small increase cannot counteract the decline of the other organic acid biomarkers (Table 1). While in the phosphate-deprived group (Figure 7) all the organic acid biomarkers decreased, leading to a down-regulation of TCA activity.
Compared with the widely studied model strain, Chlamydomonas reinhardtii, the metabolic responses of C. nivalis under nutrient stresses were not well known. C. nivalis was originated from polar region and similar extreme environments, which were different from that of C. reinhardtii. It was reported that the snow alga Chlamydomonas nivalis is cryotolerant in comparison with the mesophilic C. reinhardtii, and it has the special mechanisms for adapting photosynthesis to low temperatures and changing thylakoid lipid composition to enhance the membrane fluidity [46]. The metabolic responses of C. nivalis to nitrogen and phosphate stresses in our study were investigated and compared with that of C. reinhardtii in previous studies. In our case, the metabolic pathways of nitrogen uptake and lipid biosynthesis were enhanced significantly in C. nivalis under nitrogen deprivation. The increased mobilization of C-reservers and the accumulation of reducing equivalents for the fatty acid biosynthesis were also observed in our study. These were similar with that of C. reinhardtii [6,8]. However, it’s worth to note that, for example, the arginine synthesis pathway in C. reinhardtii was up-regulated significantly after nitrogen deprivation in order to take exogenous available nitrogen rapidly [8]. This phenomenon was not found in C. nivalis in our study, which indicated the different regulation of the nitrogen uptake between the two alga species.
In contrast to nitrogen deprivation, phosphate deprivation has more effects on central carbon and energy metabolism in C. nivalis. Previous studies have shown that under phosphate starvation conditions, C. reinhardtii has evolved the adaptive mechanisms to recycle phosphate from the intracellular stocks, such as polyphosphate, pyrophosphate and other metabolites containing phosphate [47]. The increasing content of pyrophosphate was also found in C. nivalis under phosphate deprivation; however, the large accumulation of these metabolites depressed the cell growth to some extent. Phosphate deprivation resulted in the decrease of the organic acid biomarkers which affected TCA cycle in C. nivalis. Interestingly, it was found in plants that the release of organic acid is helpful for mobilizing phosphate and increasing its concentrations in soil solution [48]. As to the metabolism of amino acids, the distinct difference between the two algae was that cysteine and its pool were significantly increased in C. reinhardtii under phosphate deprivation [49] but not in C. nivalis in this study. Cysteine and its derivatives were the major metabolites of sulfur assimilation in C. reinhardtii [47]. The continual supply of sulfur is necessary to cellular functions and metabolic responses to nutrient-deprived stress. These results indicated sulfur metabolisms in the two algae might be different, which were worthy of further investigating. Therefore, it could be concluded that our work provides new insight for a more fundamental understanding of the metabolome and the regulation of metabolism in Chlamydomonas nivalis.
Overall, the change rates in the abundance of the responding biomarkers only represent the instantaneous state fixed in multiple metabolic pathways at the same time point. The regulatory strategy in each pathway becomes extremely complex by connecting it to a global metabolic network. By determining the limited responding biomarkers under stresses, it was very feasible to realize the visualization of the metabolic network. However, there are still some limitations on the interpretation due to the complications of the metabolism. In general, the metabolites are not only involved in multiple pathways, but are also presented in different subcellular pools. Another technical factor is that the current system does not offer enough precision to detect all the metabolites and/or to identify all the detected metabolites. Despite that, the metabolic pathway analysis and reconstruction, based on the responding biomarkers in the metabolomic analysis, is still a very useful tool for providing substantial information to elucidate the key regulatory points in the network and to reveal the global regulation model in microalgae under environmental stress.

3. Experimental Section

3.1. Materials

The snow alga Chlamydomonas nivalis LB-2824 were purchased from the Culture Collection of Algae at the University of Texas at Austin (UTEX), Austin, TX, USA. HPLC-grade methanol and chloroform were purchased from Merck Co. (Darmstadt, Germany), and the other analytical-grade chemicals and reagents were purchased from Jinke (Shanghai, China) Co., Ltd. Deionized water was obtained by using a Milli-Q water system (Millipore, Bedford, MA, USA). Methoxyamine hydrochloride, N-Methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) and pyridine were purchased from Sigma-Aldrich (St. Louis, MO, USA).

3.2. Algal Culture, Stress Treatment and Biomass Harvest

The seed culture procedure, the nitrate- or phosphate-deprived stress treatment and the biomass harvest are reported in our previous work [5]. All the cultures were performed in triplicate for the experiments.

3.3. Extraction and Derivatization of the Cellular Compounds

Four milligrams of lyophilized algal powder and 1 mL of solvent (methanol:chloroform:water, 10:3:1) were mixed to extract the cellular compounds [50]. The extracts were derivatized by adding 10 µL of a solution of methoxyamine hydrochloride in pyridine (40 mg/mL) followed by shaking for 90 min at low temperature (30 °C) to prevent damage to the chemical groups (aldehyde and ketone). Sequentially, MSTFA (100 µL) was added into each reaction, and it was shaken at 37 °C for 30 min. As a negative control, one blank group without algal powder was performed in the whole procedure.

3.4. Metabolome Profile Analyses by Gas Chromatography Time-of-Flight Mass Spectrometry (GC/TOF-MS)

3.4.1. Gas Chromatography Procedure

Agilent 6890 gas chromatography (Agilent, Santa Clara, CA, USA) equipped with a Gerstel CIS4-dual MPS Injector (Gerstel, Muehlheim, Germany) was employed for the metabolome profile analysis. Aliquot of the samples (0.5-µL) were splitless injected at 10 µL/s into Rtx-5 Sil MS column (30 m × 250 µm, 95% dimethyl and 5% diphenyl polysiloxane, 0.25 µm; Restek, Bellefonte, PA, USA). The injector temperature was raised from 50 to 275 °C by 12 °C/s and was held for 3 min. The temperature of the oven was held constant at 50 °C for 1 min, was raised to 330 °C by 20 °C/min and was then held for 5 min. The carrier gas was 99.9999% pure helium at a constant flow rate of 1 mL/min.

3.4.2. Mass Spectrometry

A Leco Pegasus III time-of-flight mass spectrometer (LECO Corporation, St. Joseph, MI, USA) was used to analyze the metabolome profile. The shift temperature was set at 280 °C, and the ion source temperature was set to 250 °C in order to achieve the electron impact ionization at 70 eV. The acquisition rate was 17 spectra/s, and the scan mass ranged from 85 to 500 Da.

3.4.3. Data Processing and Normalization

All the peaks of the detected metabolites were analyzed by Leco ChromaTOF software V2.32 (LECO Corporation, St. Joseph, MI, USA) followed by processing with the programmed database BinBase according to a former report [51]. The metabolite data were normalized according to our previous study [5]. Briefly, the ion intensities of each peak detected were normalized to the sum of the peak intensities and were then normalized to the cell number in the sample. The normalized data were multiplied by 10,000 and were transformed to log10 for establishing a normal distribution in each sample for further multivariate statistical analyses.

3.5. Multivariate Statistical Analyses

SIMCA-P+ V12.0 software (Umetrics AB, Umea, Sweden) was used to perform the OPLS-DA model. Pareto scaling was first applied to pretreat the data sets to make them more comparable in magnitude for the analysis. OPLS-DA is a useful statistical procedure to discriminate the different classes by comparing and selecting the informative variables, which adopt to class membership [5]. Permutation tests with 200 iterations were performed to evaluate the possible over-fitting of data. R2X, R2Y and Q2, varied between 0 and 1, are the primary parameters of the model, which indicates the fitting degree of the data. An intercept value of Q2 > 0.05 means that overfitting in the original model exists [31]. Potential responding biomarkers were identified based on the S-plot, the column loading plot and the VIP value [26,27,28]. The change rate in the percentage of certain biomarkers was calculated by comparing the abundance of the corresponding metabolites in the stress groups with the control. The significance test for the variation of the biomarker’s abundance was implemented by a student’s t-test.

Acknowledgments

This work was supported by National Sciences Foundation of China (Grant No. 31370383, 31270085), National Hi-tech Research and Development Program (863 Project) (Grant No. 2013AA065802) and the Major State Basic Research Development Program of China (973 Project) (Grant Nos. 2011CB200904). We appreciate Oliver Fiehn in Genome Center, University of California at Davis, USA for providing GC/TOF-MS, BinBase database and compound identification.

Author Contributions

Na Lu and Dong Wei conceived, designed and performed the experiments, participated in the statistical work, and wrote the paper; Jun-Hui Chen participated in creating figures, performed the data analysis and discussion, and modified the paper; Feng Chen and Gu Chen provided the helpful discussion; All authors contributed to review of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wiencke, C.; Clayton, M.N. Biology of polar benthic algae. Bot. Mar. 2009, 52, 479–481. [Google Scholar] [CrossRef]
  2. Schwarz, R.; Forchhammer, K. Acclimation of unicellular cyanobacteria to macronutrient deficiency: Emergence of a complex network of cellular responses. Microbiol. SGM 2005, 151, 2503–2514. [Google Scholar] [CrossRef] [PubMed]
  3. Gaude, N.; Brehelin, C.; Tischendorf, G.; Kessler, F.; Doermann, P. Nitrogen deficiency in arabidopsis affects galactolipid composition and gene expression and results in accumulation of fatty acid phytyl esters. Plant J. 2007, 49, 729–739. [Google Scholar] [CrossRef] [PubMed]
  4. El-Sheek, M.M.; Rady, A.A. Effect of phosphorus starvation on growth, photosynthesis and some metabolic processes in the unicellular green alga chlorella kessleri. Phyton 1995, 35, 139–151. [Google Scholar]
  5. Lu, N.; Wei, D.; Chen, F.; Yang, S.-T. Lipidomic profiling reveals lipid regulation in the snow alga Chlamydomonas nivalis in response to nitrate or phosphate deprivation. Process Biochem. 2013, 48, 605–613. [Google Scholar] [CrossRef]
  6. Schmollinger, S.; Muehlhaus, T.; Boyle, N.R.; Blaby, I.K.; Casero, D.; Mettler, T.; Moseley, J.L.; Kropat, J.; Sommer, F.; Strenkert, D.; et al. Nitrogen-sparing mechanisms in Chlamydomonas affect the transcriptome, the proteome, and photosynthetic metabolism. Plant Cell 2014, 26, 1410–1435. [Google Scholar] [CrossRef] [PubMed]
  7. Soga, T.; Ohashi, Y.; Ueno, Y.; Naraoka, H.; Tomita, M.; Nishioka, T. Quantitative metabolome analysis using capillary electrophoresis mass spectrometry. J. Proteome Res. 2003, 2, 488–494. [Google Scholar] [CrossRef] [PubMed]
  8. Park, J.-J.; Wang, H.; Gargouri, M.; Deshpande, R.R.; Skepper, J.N.; Holguin, F.O.; Juergens, M.T.; Shachar-Hill, Y.; Hicks, L.M.; Gang, D.R. The response of Chlamydomonas reinhardtii to nitrogen deprivation: A systems biology analysis. Plant J. 2015, 81, 611–624. [Google Scholar] [CrossRef] [PubMed]
  9. Feng, T.-Y.; Yang, Z.-K.; Zheng, J.-W.; Xie, Y.; Li, D.-W.; Murugan, S.B.; Yang, W.-D.; Liu, J.-S.; Li, H.-Y. Examination of metabolic responses to phosphorus limitation via proteomic analyses in the marine diatom phaeodactylum tricornutum. Sci. Rep. 2015, 5. [Google Scholar] [CrossRef] [PubMed]
  10. Bielecka, M.; Watanabe, M.; Morcuende, R.; Scheible, W.-R.; Hawkesford, M.J.; Hesse, H.; Hoefgen, R. Transcriptome and metabolome analysis of plant sulfate starvation and resupply provides novel information on transcriptional regulation of metabolism associated with sulfur, nitrogen and phosphorus nutritional responses in arabidopsis. Front. Plant Sci. 2015, 5, 1–18. [Google Scholar] [CrossRef] [PubMed]
  11. Fiehn, O. Metabolomics—The link between genotypes and phenotypes. Plant Mol. Biol. 2002, 48, 155–171. [Google Scholar] [CrossRef] [PubMed]
  12. Beisken, S.; Eiden, M.; Salek, R.M. Getting the right answers: Understanding metabolomics challenges. Expert Rev. Mol. Diagn. 2015, 15, 97–109. [Google Scholar] [CrossRef] [PubMed]
  13. Peterson, A.C.; Balloon, A.J.; Westphall, M.S.; Coon, J.J. Development of a GC/quadrupole-orbitrap mass spectrometer, part II: New approaches for discovery metabolomics. Anal. Chem. 2014, 86, 10044–10051. [Google Scholar] [CrossRef] [PubMed]
  14. Obata, T.; Fernie, A.R. The use of metabolomics to dissect plant responses to abiotic stresses. Cell. Mol. Life Sci. 2012, 69, 3225–3243. [Google Scholar] [CrossRef] [PubMed]
  15. Veyel, D.; Erban, A.; Fehrle, I.; Kopka, J.; Schroda, M. Rationales and approaches for studying metabolism in eukaryotic microalgae. Metabolites 2014, 4, 184–217. [Google Scholar] [CrossRef] [PubMed]
  16. Lee, D.Y.; Fiehn, O. High quality metabolomic data for Chlamydomonas reinhardtii. Plant Methods 2008, 4. [Google Scholar] [CrossRef] [PubMed]
  17. Wiklund, S.; Johansson, E.; Sjostrom, L.; Mellerowicz, E.J.; Edlund, U.; Shockcor, J.P.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using opls class models. Anal. Chem. 2008, 80, 115–122. [Google Scholar] [CrossRef] [PubMed]
  18. Cajka, T.; Fiehn, O. Toward merging untargeted and targeted methods in mass spectrometry-based metabolomics and lipidomics. Anal. Chem. 2016, 88, 524–545. [Google Scholar] [CrossRef] [PubMed]
  19. Wienkoop, S.; Weiss, J.; May, P.; Kempa, S.; Irgang, S.; Recuenco-Munoz, L.; Pietzke, M.; Schwemmer, T.; Rupprecht, J.; Egelhofer, V.; et al. Targeted proteomics for Chlamydomonas reinhardtii combined with rapid subcellular protein fractionation, metabolomics and metabolic flux analyses. Mol. Biosyst. 2010, 6, 1018–1031. [Google Scholar] [CrossRef] [PubMed]
  20. Diaz, G.; Melis, M.; Batetta, B.; Angius, F.; Falchi, A.M. Hydrophobic characterization of intracellular lipids in situ by nile red red/yellow emission ratio. Micron 2008, 39, 819–824. [Google Scholar] [CrossRef] [PubMed]
  21. Fernie, A.R.; Obata, T.; Allen, A.E.; Araujo, W.L.; Bowler, C. Leveraging metabolomics for functional investigations in sequenced marine diatoms. Trends Plant Sci. 2012, 17, 395–403. [Google Scholar] [CrossRef] [PubMed]
  22. Vidoudez, C.; Pohnert, G. Comparative metabolomics of the diatom skeletonema marinoi in different growth phases. Metabolomics 2012, 8, 654–669. [Google Scholar] [CrossRef]
  23. Kanehisa, M. Toward pathway engineering: A new database of genetic and molecular pathways. Sci. Technol. Jpn. 1996, 59, 34–38. [Google Scholar]
  24. Pandhal, J.; Noirel, J.; Wright, P.C.; Biggs, C.A. A systems biology approach to investigate the response of Synechocystis sp. PCC6803 to a high salt environment. Saline Syst. 2009, 5. [Google Scholar] [CrossRef] [PubMed]
  25. Hockin, N.L.; Mock, T.; Mulholland, F.; Kopriva, S.; Malin, G. The response of diatom central carbon metabolism to nitrogen starvation is different from that of green algae and higher plants. Plant Physiol. 2012, 158, 299–312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lu, N.; Wei, D.; Jiang, X.-L.; Chen, F.; Yang, S.-T. Regulation of lipid metabolism in the snow alga Chlamydomonas nivalis in response to NACL stress: An integrated analysis by cytomic and lipidomic approaches. Process Biochem. 2012, 47, 1163–1170. [Google Scholar] [CrossRef]
  27. Lu, N.; Wei, D.; Jiang, X.-L.; Chen, F.; Yang, S.-T. Fatty acids profiling and biomarker identification in snow alga Chlamydomonas nivalis by NACL stress using GC/MS and multivariate statistical analysis. Anal. Lett. 2012, 45, 1172–1183. [Google Scholar] [CrossRef]
  28. Lu, N.; Wei, D.; Chen, F.; Yang, S.-T. Lipidomic profiling and discovery of lipid biomarkers in snow alga Chlamydomonas nivalis under salt stress. Eur. J. Lipid Sci. Technol. 2012, 114, 253–265. [Google Scholar] [CrossRef]
  29. Guerrini, F.; Cangini, M.; Boni, L.; Trost, P.; Pistocchi, R. Metabolic responses of the diatom Achnanthes brevipes (bacillariophyceae) to nutrient limitation. J. Phycol. 2000, 36, 882–890. [Google Scholar] [CrossRef]
  30. Lai, J.; Yu, Z.; Song, X.; Cao, X.; Han, X. Responses of the growth and biochemical composition of prorocentrum donghaiense to different nitrogen and phosphorus concentrations. J. Exp. Mar. Biol. Ecol. 2011, 405, 6–17. [Google Scholar] [CrossRef]
  31. Kang, J.; Choi, M.-Y.; Kang, S.; Kwon, H.N.; Wen, H.; Lee, C.H.; Park, M.; Wiklund, S.; Kim, H.J.; Kwon, S.W.; et al. Application of a 1H nuclear magnetic resonance (NMR) metabolomics approach combined with orthogonal projections to latent structure-discriminant analysis as an efficient tool for discriminating between korean and chinese herbal medicines. J. Agric. Food Chem. 2008, 56, 11589–11595. [Google Scholar] [CrossRef] [PubMed]
  32. Van de Wiel, C.C.M.; van der Linden, C.G.; Scholten, O.E. Improving phosphorus use efficiency in agriculture: Opportunities for breeding. Euphytica 2016, 207, 1–22. [Google Scholar] [CrossRef]
  33. Blatti, J.L.; Michaud, J.; Burkart, M.D. Engineering fatty acid biosynthesis in microalgae for sustainable biodiesel. Curr. Opin. Chem. Biol. 2013, 17, 496–505. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, Z.T.; Ullrich, N.; Joo, S.; Waffenschmidt, S.; Goodenough, U. Algal lipid bodies: Stress induction, purification, and biochemical characterization in wild-type and starchless Chlamydomonas reinhardtii. Eukaryot. Cell 2009, 8, 1856–1868. [Google Scholar] [CrossRef] [PubMed]
  35. Miller, R.; Wu, G.; Deshpande, R.R.; Vieler, A.; Gaertner, K.; Li, X.; Moellering, E.R.; Zaeuner, S.; Cornish, A.J.; Liu, B.; et al. Changes in transcript abundance in Chlamydomonas reinhardtii following nitrogen deprivation predict diversion of metabolism. Plant Physiol. 2010, 154, 1737–1752. [Google Scholar] [CrossRef] [PubMed]
  36. Flynn, K. Algal carbon-nitrogen metabolism: A biochemical basis for modelling the interactions between nitrate and ammonium uptake. J. Plankton Res. 1991, 13, 373–387. [Google Scholar] [CrossRef]
  37. Ito, T.; Tanaka, M.; Shinkawa, H.; Nakada, T.; Ano, Y.; Kurano, N.; Soga, T.; Tomita, M. Metabolic and morphological changes of an oil accumulating trebouxiophycean alga in nitrogen-deficient conditions. Metabolomics 2013, 9, 178–187. [Google Scholar] [CrossRef] [PubMed]
  38. Theodorou, M.E.; Elrifi, I.R.; Turpin, D.H.; Plaxton, W.C. Effects of phosphorus limitation on respiratory metabolism in the green alga selenastrum minutum. Plant Physiol. 1991, 95, 1089–1095. [Google Scholar] [CrossRef] [PubMed]
  39. Juneja, A.; Ceballos, R.M.; Murthy, G.S. Effects of environmental factors and nutrient availability on the biochemical composition of algae for biofuels production: A review. Energies 2013, 6, 4607–4638. [Google Scholar] [CrossRef]
  40. Siderius, M.; Musgrave, A.; Ende, H.; Koerten, H.; Cambier, P.; Meer, P. Chlamydomonas eugametos (Chlorophyta) stores phosphate in polyphosphate bodies together with calcium. J. Phycol. 1996, 32, 402–409. [Google Scholar] [CrossRef]
  41. Hebeler, M.; Hentrich, S.; Mayer, A.; Leibfritz, D.; Grimme, L. Phosphate regulation and compartmentation in Chlamydomonas reinhardtii studied by in vivo p-31-NMR. Photosynth. Res. 1992, 34, 199. [Google Scholar]
  42. Loera-Quezada, M.M.; Leyva-Gonzalez, M.A.; Lopez-Arredondo, D.; Herrera-Estrella, L. Phosphite cannot be used as a phosphorus source but is non-toxic for microalgae. Plant Sci. 2015, 231, 124–130. [Google Scholar] [CrossRef] [PubMed]
  43. Muratsubaki, H. Regulation of reductive production of succinate under anaerobic conditions in baker's yeast. J. Biochem. 1987, 102, 705–714. [Google Scholar] [PubMed]
  44. Sirithanakorn, C.; Adina-Zada, A.; Wallace, J.C.; Jitrapakdee, S.; Attwood, P.V. Mechanisms of inhibition of Rhizobium etli pyruvate carboxylase by l-aspartate. Biochemistry 2014, 53, 7100–7106. [Google Scholar] [CrossRef] [PubMed]
  45. Kanehisa, M.; Goto, S.; Sato, Y.; Furumichi, M.; Tanabe, M. Kegg for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40, 109–114. [Google Scholar] [CrossRef] [PubMed]
  46. Lukes, M.; Prochazkova, L.; Shmidt, V.; Nedbalova, L.; Kaftan, D. Temperature dependence of photosynthesis and thylakoid lipid composition in the red snow alga Chlamydomonas cf. nivalis (Chlorophyceae). FEMS Microbiol. Ecol. 2014, 89, 303–315. [Google Scholar] [CrossRef] [PubMed]
  47. Irihimovitch, V.; Yehudai-Resheff, S. Phosphate and sulfur limitation responses in the chloroplast of Chlamydomonas reinhardtii. FEMS Microbiol. Lett. 2008, 283, 1–8. [Google Scholar] [CrossRef] [PubMed]
  48. Plaxton, W.C.; Tran, H.T. Metabolic adaptations of phosphate-starved plants. Plant Physiol. 2011, 156, 1006–1015. [Google Scholar] [CrossRef] [PubMed]
  49. Bolling, C.; Fiehn, O. Metabolite profiling of Chlamydomonas reinhardtii under nutrient deprivation. Plant Physiol. 2005, 139, 1995–2005. [Google Scholar] [CrossRef] [PubMed]
  50. Fiehn, O.; Kind, T. Metabolite profiling in blood plasma. Methods Mol. Biol. 2007, 358, 3–17. [Google Scholar] [PubMed]
  51. Kind, T.; Wohlgemuth, G.; Lee, D.Y.; Lu, Y.; Palazoglu, M.; Shahbaz, S.; Fiehn, O. Fiehnlib: Mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-fight gas chromatography/mass spectrometry. Anal. Chem. 2009, 81, 10038–10048. [Google Scholar] [CrossRef] [PubMed]
Ijms 17 00694 g001 1024
Figure 1. Overview of the identified metabolites in the classes and their abundance in the stress groups and the control group. *, ** mean significant difference (p < 0.05 or p < 0.01).
Figure 1. Overview of the identified metabolites in the classes and their abundance in the stress groups and the control group. *, ** mean significant difference (p < 0.05 or p < 0.01).
Ijms 17 00694 g001
Ijms 17 00694 g002 1024
Figure 2. Score scatter plot generated by the orthogonal partial least squares discriminant analysis (OPLS-DA) model from the peak area of the identified metabolites in the different groups. (▲) the control; (■) nitrate deprivation for 6 and 12 h; (□) nitrate deprivation for 24, 48 and 72 h; (◆) phosphate deprivation for 6 and 12 h; (◇) phosphate deprivation for 24, 48 and 72 h. t[1], the score of the first predictive component to explain the largest variation; t[2], the score of the first orthogonal component to explain the largest orthogonal variation.
Figure 2. Score scatter plot generated by the orthogonal partial least squares discriminant analysis (OPLS-DA) model from the peak area of the identified metabolites in the different groups. (▲) the control; (■) nitrate deprivation for 6 and 12 h; (□) nitrate deprivation for 24, 48 and 72 h; (◆) phosphate deprivation for 6 and 12 h; (◇) phosphate deprivation for 24, 48 and 72 h. t[1], the score of the first predictive component to explain the largest variation; t[2], the score of the first orthogonal component to explain the largest orthogonal variation.
Ijms 17 00694 g002
Ijms 17 00694 g003 1024
Figure 3. Related plots of the OPLS-DA model between the nitrate-deprived group and the control. (a) Score scatter plot for metabolite discrimination. (▲) the control; (■) nitrate-deprived group for 6 and 12 h; (□) nitrate-deprived group for 24, 48 and 72 h; (b) VIP plot (with VIP > 1 metabolites); (c) s-plot; square marked metabolites are those with |p(corr)| > 0.6; (d) Loading column plot of the potential metabolites biomarkers.
Figure 3. Related plots of the OPLS-DA model between the nitrate-deprived group and the control. (a) Score scatter plot for metabolite discrimination. (▲) the control; (■) nitrate-deprived group for 6 and 12 h; (□) nitrate-deprived group for 24, 48 and 72 h; (b) VIP plot (with VIP > 1 metabolites); (c) s-plot; square marked metabolites are those with |p(corr)| > 0.6; (d) Loading column plot of the potential metabolites biomarkers.
Ijms 17 00694 g003
Ijms 17 00694 g004 1024
Figure 4. Related plots of the OPLS-DA model between the phosphate-deprived group and the control. (a) Score scatter plot for metabolite discrimination. (▲) the control; (◆) phosphate-deprived group for 6 and 12 h; (◇) phosphate-deprived group for 24, 48 and 72 h; (b) VIP plot (with VIP > 1 metabolites); (c) s-plot; square marked metabolites are those with |p(corr)| > 0.6; (d) Loading column plot of the potential metabolites biomarkers.
Figure 4. Related plots of the OPLS-DA model between the phosphate-deprived group and the control. (a) Score scatter plot for metabolite discrimination. (▲) the control; (◆) phosphate-deprived group for 6 and 12 h; (◇) phosphate-deprived group for 24, 48 and 72 h; (b) VIP plot (with VIP > 1 metabolites); (c) s-plot; square marked metabolites are those with |p(corr)| > 0.6; (d) Loading column plot of the potential metabolites biomarkers.
Ijms 17 00694 g004
Ijms 17 00694 g005 1024
Figure 5. The abundance in the biomarker classes responding to the stresses compared with the control. *, ** mean significant difference (p < 0.05 or p < 0.01).
Figure 5. The abundance in the biomarker classes responding to the stresses compared with the control. *, ** mean significant difference (p < 0.05 or p < 0.01).
Ijms 17 00694 g005
Ijms 17 00694 g006 1024
Figure 6. Metabolic pathway and regulation visualized by Kyoto encyclopedia of genes and genomes (KEGG) under nitrate-deprived stress. Direct connections between the metabolites are indicated with solid line arrows, and putative connections between the metabolites are indicated with dash line arrows. The red words in red solid line boxes represent an increase in the abundance of the responding biomarkers, the blue words in blue dash line boxes represent a decrease in the abundance, and the black words represent no significant change. The percentages represent the change rates of the biomarkers by up- or down-regulation.
Figure 6. Metabolic pathway and regulation visualized by Kyoto encyclopedia of genes and genomes (KEGG) under nitrate-deprived stress. Direct connections between the metabolites are indicated with solid line arrows, and putative connections between the metabolites are indicated with dash line arrows. The red words in red solid line boxes represent an increase in the abundance of the responding biomarkers, the blue words in blue dash line boxes represent a decrease in the abundance, and the black words represent no significant change. The percentages represent the change rates of the biomarkers by up- or down-regulation.
Ijms 17 00694 g006
Ijms 17 00694 g007 1024
Figure 7. Metabolic pathway and regulation visualized by KEGG under phosphate-deprived stress. Direct connections between the metabolites are indicated with solid line arrows, and putative connections between the metabolites are indicated with dash line arrows. The red words in red solid line boxes represent an increase in the abundance of the responding biomarkers, the blue words in blue dash line boxes represent a decrease in the abundance, and the black words represent no significant change. The percentages represent the change rates of the responding biomarkers by up- or down-regulation.
Figure 7. Metabolic pathway and regulation visualized by KEGG under phosphate-deprived stress. Direct connections between the metabolites are indicated with solid line arrows, and putative connections between the metabolites are indicated with dash line arrows. The red words in red solid line boxes represent an increase in the abundance of the responding biomarkers, the blue words in blue dash line boxes represent a decrease in the abundance, and the black words represent no significant change. The percentages represent the change rates of the responding biomarkers by up- or down-regulation.
Ijms 17 00694 g007
Table
Table 1. Thirty of the selected metabolites with a VIP > 1 and a |p (corr)| > 0.6 as responding biomarkers from the nitrate-deprived group vs. the control group.
Table 1. Thirty of the selected metabolites with a VIP > 1 and a |p (corr)| > 0.6 as responding biomarkers from the nitrate-deprived group vs. the control group.
NumberIdentified MetabolitesVIPp[1]p(corr)[1]Change Rate in Abundance, % (N Deprivation vs. the Control)
1Linolenic acid5.480.220.655.13
2α-Ketoglutaric acid3.240.280.94−55.11
3Stearic acid2.60−0.20−0.8867.15
4Putrescine2.480.240.97−29.31
5Glutamine2.420.220.95−21.44
6Asparagine2.330.200.93−58.90
7Palmitic acid2.01−0.20−0.9767.00
8Uracil1.920.200.99−11.27
9Adenosine1.870.160.93−28.98
10Glutamic acid1.860.180.97−35.24
11Nicotinamide1.750.190.99−8.17
12Adenine1.740.170.96−30.60
13Succinic acid1.710.170.948.36
14Valine1.650.170.99−6.68
15Isoleucine1.61−0.10−0.8075.65
16Trisaccharide1.590.171.00−2.16
17Citric acid1.520.130.91−33.14
18Pentadecanoic acid1.45−0.13−0.9575.54
19Alanine1.430.150.98−56.06
20Malic acid1.410.150.99−11.30
21Fumaric acid1.360.130.94−34.74
22Aspartic acid1.350.110.92−5.22
23Proline1.340.130.97−1.03
24Serine1.300.110.93−9.06
25Phenylalanine1.190.120.97−29.42
26Histidine1.13−0.12−0.9986.94
27Tryptophan1.11−0.12−1.0085.84
28Shikimic acid1.110.121.00−35.76
29Leucine1.090.110.97−27.64
30Lysine1.010.110.99−4.38
VIP: variable influence on the projection; p[1]: the loading of the first predictive component to explain the largest variation; p(corr)[1]: the correlation coefficient between the metabolite and the score of the first predictive component to explain the largest variation; N: nitrogen.
Table
Table 2. Thirty-nine of the selected metabolites with a VIP > 1 and a |p(corr)| > 0.6 as responding biomarkers from the phosphate-deprived group vs. the control group.
Table 2. Thirty-nine of the selected metabolites with a VIP > 1 and a |p(corr)| > 0.6 as responding biomarkers from the phosphate-deprived group vs. the control group.
NumberIdentified MetabolitesVIPp[1]p(corr)[1]Change Rate in Abundance, % (P deprivation vs. the Control)
1Linolenic acid3.420.300.98−22.77
2Stearic acid2.89−0.25−0.9647.74
3Trisaccharide2.75−0.12−0.8824.77
4Adenosine2.530.210.98−43.47
5Glutamic acid2.490.210.9710.29
6α-Ketoglutaric acid2.340.250.99−22.15
7Putrescine2.160.241.00−19.51
8Fumaric acid1.840.170.99−55.15
9Trehalose1.810.100.93−20.32
10Threonic acid1.790.110.921.05
11Shikimic acid1.750.140.97−42.50
12Nicotinamide1.630.150.94−3.00
13Valine1.600.150.98−4.15
14Succinic acid1.510.160.94−1.38
15Proline1.490.110.95−1.51
16Aspartic acid1.49−0.13−0.9742.85
17Benzoic acid1.480.050.646.81
18Histidine1.460.130.97−23.75
19Uracil1.420.150.99−1.08
20Sucrose1.410.040.739.17
21Malic acid1.390.100.89−0.39
22Serine1.370.100.93−5.55
23Threitol1.360.050.71−1.11
24Phenylalanine1.350.150.99−39.56
25Erythronic acid lactone1.300.080.801.74
26Leucine1.290.130.99−38.43
27Adenosine-5-phosphate1.270.120.97−4.98
28Palmitic acid1.23−0.13−0.7925.12
293-Phosphoglycerate1.22−0.12−0.99297.76
30Fructose-6-phosphate1.210.130.99−21.05
31Inositol-4-monophosphate1.210.120.99−37.26
32Alanine1.190.130.99−27.64
33Pyrophosphate1.13−0.12−0.9941.90
34Gluconic acid1.090.080.93−35.19
35Asparagine1.090.120.99−3.23
36Pentadecanoic acid1.07−0.10−0.9632.98
37Tyrosine1.020.100.971.73
38Threonine1.00−0.08−0.9423.50
39Cellobiose1.000.100.948.64
VIP: variable influence on the projection; p[1]: the loading of the first predictive component to explain the largest variation; p(corr)[1]: the correlation coefficient between the metabolite and the score of the first predictive component to explain the largest variation; P: phosphate.

Share and Cite

MDPI and ACS Style

Lu, N.; Chen, J.-H.; Wei, D.; Chen, F.; Chen, G. Global Metabolic Regulation of the Snow Alga Chlamydomonas nivalis in Response to Nitrate or Phosphate Deprivation by a Metabolome Profile Analysis. Int. J. Mol. Sci. 2016, 17, 694. https://doi.org/10.3390/ijms17050694

AMA Style

Lu N, Chen J-H, Wei D, Chen F, Chen G. Global Metabolic Regulation of the Snow Alga Chlamydomonas nivalis in Response to Nitrate or Phosphate Deprivation by a Metabolome Profile Analysis. International Journal of Molecular Sciences. 2016; 17(5):694. https://doi.org/10.3390/ijms17050694

Chicago/Turabian Style

Lu, Na, Jun-Hui Chen, Dong Wei, Feng Chen, and Gu Chen. 2016. "Global Metabolic Regulation of the Snow Alga Chlamydomonas nivalis in Response to Nitrate or Phosphate Deprivation by a Metabolome Profile Analysis" International Journal of Molecular Sciences 17, no. 5: 694. https://doi.org/10.3390/ijms17050694

APA Style

Lu, N., Chen, J. -H., Wei, D., Chen, F., & Chen, G. (2016). Global Metabolic Regulation of the Snow Alga Chlamydomonas nivalis in Response to Nitrate or Phosphate Deprivation by a Metabolome Profile Analysis. International Journal of Molecular Sciences, 17(5), 694. https://doi.org/10.3390/ijms17050694

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop