Comparative Metabolomic and Network Analysis of the Hyperaccumulation of Astaxanthin in Haematococcus pluvialis under the Different Stress Conditions

1. Introduction

Astaxanthin (C₄₀H₅₂O₄) is a ketocarotenoid belonging to the terpene compounds, which is recognized as one of the most powerful antioxidants in nature [1]. The ability to neutralize free radicals and scavenge singlet oxygen of astaxanthin is much stronger than a range of common antioxidants such as zeaxanthin, lycopene, β-carotene or vitamin E [2]. Moreover, excellent biosafety properties of natural astaxanthin has been confirmed, which made it apply to health products, dietory supplements and aquaculture. According to the recent research, astaxanthin has been showed obviously beneficial effects on human health and used to the adjuvant treatment for a range of diseases, such as cataract, atherosclerosis and cardiovascular disease [3]. In the field of aquaculture, astaxanthin has been widely used as the colorant and immunopotentiator, which maintained body color of ornamental fish, and provided immune enhancement ability for aquaculture economic animals. Given the excellent biological activity and a wide range of applications, the astaxanthin market will be estimated up to $2.57 billion by 2025 [4].

Haematococcus pluvialis is a type of unicellular green algae, which is distributed in variable freshwater, saline and semi-saline habitats worldwide [5,6,7,8,9]. Furthermore, H. pluvialis is widely recognized as one of the best natural sources of astaxanthin. Under the stressed conditions, the planktonic cells of H. pluvialis gradually lose their flagella and then formed immotile spores, in which astaxanthin would be generated and accumulated up to 4.0% of the dry weight [10,11,12]. Moreover, variable strategies to accelerate the astaxanthin accumulation in H. pluvialis have been conducted, contained regulation of the environmental stress factors, artificial mutation breeding, external carbon source or nutrient enrichment, and induction of exogenous stimuli [13,14,15,16].

Genomics, transcriptomics and metabolomics techniques have been well proved to be the effective approach to analyse the biosynthesis mechanism of the secondary metabolites in H. pluvialis [17,18,19,20]. However, metabolomics is considered to be the most reflective analysis method for the physiology of organisms among all the “omics” techniques, because metabolites are the terminal products after a series of complex physiological metabolic processes. Currently, nuclear magnetic resonance (NMR), mass spectrometry (MS), chromatography spectroscopy (GC), GC-MS and LC-MS techniques have been used to metabolomics study [21,22,23]. Thus, a series of research results related to a global insight into the regulation networks involved in the secondary metabolites such as β-carotene, astaxanthin and lipids in H. pluvialis has been showed [24,25,26].

In this study, the optimum NaCl and NaNO₃ concentrations for the astaxanthin accumulation in H. pluvialis were determined, and then under which the content of intracellular astaxanthin were measured. Subsequently, the metabolomics analysis using LC-MS/MS technique and enrichment analysis based on KEGG database were conducted. Furthermore, the differential metabolic patterns induced by high salinity and nitrogen deficiency were interpreted. The results provided the comparative metabolomic view of astaxanthin accumulation in H. pluvialis under the different stress conditions, moreover showed a novel insights into the astaxanthin production induced by the environmental stress factors.

2. Materials and Methods

2.1. Experimental Material and Treatment

H. pluvialis was provided by Tianjin Key Lab for Aquaculture Ecology and Cultivation. The microalgal was cultivated in BBM medium at 22±1℃ in a light incubator with light intensity of 55μmol/(m²·s) and light-dark cycle of 12hL: 12hD. During the culture period, the flasks were shaken six times each day to avoid microalgal cells to attach to the walls or sink. H. pluvialis in the exponential growth stage were used for the experiments.

75mL cultures of H. pluvialis in the exponential growth stage (ca. 5×10⁵ cells·mL-1) was subcultured into 250-mL conical flasks containing 75mL medium. According to the Table 1., the control group and stress treatment groups were arranged in the experiment. Parallel triplicates were set for each experimental group. As the accumulation of astaxanthin in H. pluvialis was observed obviously from the 12th day, thus the measurement of astaxanthin was carried out every 2 days from the 12th day to the 26th day throughout the experiment. The high salinity and low nitrogen treatment groups in which the content of intracellular astaxanthin reached the highest level, compared with any other treatment group, were named as group HS and group LN respectively. In addition, the control group was abbreviated as group CK. Subsequently, the microalgae suspension on the day when the intracellular astaxanthin reached the highest level, from the group CK, group HS and group LN, were sampled to execute metabolomics analysis.

2.2. Determination of Astaxanthin in H. pluvialis

Astaxanthin in H. pluvialis was extracted and analyzed in the light of the experimental procedure of Boussiba et al.[27]. Firstly, microalgal cells were collected and treated with 5% KOH solution in 30% methanol in order to eliminate the interference of chlorophyll. Subsequently, astaxanthin was extracted by dimethylsulfoxide (DMSO) repeatedly until the precipitation turned into colourless. Finally, the absorbance of extracts was measured at 490 nm using a spectrophotometer. The blank contained DMSO only.

2.3. Metabolite Extraction and LC-MS/MS-based Metabolomic Analysis

Sixty milligrams of microalgal samples were thawed at 4°C and 100μL aliquots were mixed with 400μL pre-cooled methanol/acetonitrile (1:1, v/v) to eliminate the protein. And then the mixture was collected and centrifuged for 20 min (14000r·min^-1, 4°C). Furthermore, the samples were re-dissolved in 100μL acetonitrile/water (1:1, v/v) and fully shaked, and then centrifuged for 15min (14000r·min^-1, 4°C). The supernatants were collected to conduct the LC-MS/MS analysis.

For HILIC separation, all the samples were detected using a ACQUIY UPLC BEH column. In both ESI positive and negative modes, the mobile phase contained A and B. A was a mixed solution of ammonium acetate and ammonium hydroxide with the concentration of 25 mmol·L^-1 respectively, and B was a acetonitrile solution. For RPLC separation, a ACQUIY UPLC HSS T3 column was used. In ESI positive mode, the mobile phase contained A and B. A was a solution with 0.1% formic acid dissolved in water, and B was a solution of acetonitrile mixed with 0.1% formic acid. While in ESI negative mode, the mobile phase included A and B. A was a solution of 0.5 mmol·L^-1 ammonium fluoride dissolved in water, and B was a solution of acetonitrile. The ESI source parameters were made the following settings: Ion Source Gas1 and Gas2 as 60 respectively, curtain gas (CUR) as 30, source temperature as 600℃, IonSpray Voltage Floating (ISVF) as ± 5500 V.

2.4. Data Process and Analysis

Extraction of raw MRM data for all the metabolites was conducted, and then the raw data was conversed to netCDF format. Subsequently, the data with netCDF format was imported to XCMS based on R 3.5.2 software, and the data extraction, peak matching, baseline correction and filling gaps were conducted, therefore a file containing information such as retention time, charge-to-mass ratio and peak area that can be used for subsequent data analysis was obtained. In addition, one-way ANOVA was used to analyze the differences among the metabolites. Moreover the metabolites with P<0.05 and fold change≥1.5 were defined as the differential metabolites with statistical significance among the CK, HS and LN group.

Principal component analysis (PCA), a multivariate statistical analysis method, was used to judge the differences among the samples belonging to the variable treatments through R 3.5.2 software. And then the hierarchical clustering was conducted based on the euclidean distance algorithm. Furthermore, the heatmap of hierarchical clustering was obtained with the aim to visualize the differential metabolites. KEGG (Kyoto Encyclopedia of Genes and Genomes) is the major public database related to the pathways in the organism, which contains a series of genes and metabolites information. The metabolites detected in this study were mapped to KEGG metabolic pathways for pathway analysis and enrichment analysis, then the enriched metabolic pathways and signal transduction pathways were identified significantly in differential metabolites by pathway enrichment analysis, compared with the whole background [28].

The metabolic networks were constructed using the Advanced Network Graph Online Tools on the GENE DENOVO website (http://www.omicshare.com). In the network gragh, the metabolic pathways and the differential metabolites were represented using the green filled circles and the pink filled circles respectively, and the quatity of metabolites enriched in the specific pathway was reflected by the size of the green circle. Otherwise, the associated pathways and metabolites were linked up by the black solid lines, which exhibited the regulatory relationship among the metabolic pathways and the metabolites.

The contents of astaxanthin in the different treatments on the each sampling time were analyzed using one-way ANOVA (One-way Analysis of Variance) and Duncan’s multiple comparison based on the software SPSS 22.0. The significant level was set at P=0.05.

3. Results

3.1. Effect of NaCl Concentration on Astaxanthin Accumulation in H. pluvialis

The level of astaxanthin in H. pluvialis raised gradually, and then tended to keep steady under the high salinity stress during the experiment. The contents of astaxanthin in 5, 7.5 and 10g·L^-1 treatment groups were higher than the CK group significantly (P<0.05), while which in 12.5g·L^-1 treatment group was insignificant compared with the CK group (P>0.05). According to the Figure 1., the level of astaxanthin accumulated in H. pluvialis, in 7.5g·L^-1 NaCl group, was higher than any other groups roughly thoughout the experiment (P<0.05), furthermore in which the highest level of astaxanthin reached up to 7.51mg·L^-1 on the 24th day. However, the excessive NaCl concentration (12.5g·L^-1) was not beneficial to the accumlation of astaxanthin in H. pluvialis. This result indicated that there was a certain NaCl concentration threshold, above or below which the level of astaxanthin accumulated in H. pluvialis would drop.

3.2. Effect of NaNO₃ Concentration on Astaxanthin Accumulation in H. pluvialis

The concentration of astaxanthin in H. pluvialis increasd gradually, and then reduced slowly under the nitrogen limitation stress during the experiment. According to the Figure 2., the content of astaxanthin was the highest among the treatment groups under the nitrogen-deficient stress (0g·L^-1) thoughout the experiment (P<0.05), while those of 0.032、0.065 and 0.13g·L^-1 treatment groups were insignificant compared with the CK group (P>0.05). The astaxanthin accumulation in H. pluvialis was relatively slower before the 14th day, however the accumulation rate speeded up from the 14th day to the 22th day. The astaxanthin content in H. pluvialis reached up to 5.60mg·L^-1 at the highest level in the nitrogen-deficient group, which was 28% higher than the CK group.

3.3. Comparison of the Metabolic Characteristics in H. pluvialis under the Different Stress Conditions

PCA plot of the metabolomic profiles under the different stress conditions was generated, in which the distribution and concentration of each sample was showed clearly. According to the Figure 3., the metabolomic results of the environmental stress treatments and control group among all samples were distinct in the experiment, which suggested that the effects of high salinity and nitrogen-deficient on the global metabolomic profiles related to the astaxanthin accumulation in H. pluvialis were differential.

As shown in Table 2., a total of 132 metabolites in H. pluvialis, under the different stress conditions, were identified by LC-MS/MS technigue. In addition, one-way ANOVA was used to analyze the differences among the metabolites. Metabolites with P<0.05 and fold change≥1.5 were defined as the differential ones with statistical significance among the CK, HS and LN group. And then the heatmap of hierarchical clustering was obtained with the aim to visualize the differential metabolites (Figure 4.).

According to the judging criteria, a total of 30 differential metabolites containing 3 fatty acids and intermediates, 6 amino acids and derivatives, 1 carbohydrates, 12 nucleoside derivatives and 8 other metabolites were identified and clustered among the CK, HS and LN group. The results showed distinct downregulation of 18 metabolites in the environmental stress treatments, compared with the CK group. Especially, the levels of 3 fatty acids and intermediates, including Glycerol3P, Glyceric acid and Glycerophosphocholine, were much lower than the CK Group. Moreover the downregulation of the intermediate products of nucleotides metabolism such as AMP, UMP and 5'-Deoxyadenosine in the stress treatments was observed. However, 8 metabolites including L-Aspartic Acid, Thymine, Glycine, Guanosine, Sucrose, L-Valine, L-2-Hydroxygluterate and L-Citrulline were upregulated in the HS group. Additionally, the significant upregulation of GDP, IMP, N4-Acetylcytidine and 4-Pyridoxic acid was detected in the LN group.

3.4. KEGG Pathway Analysis

To further explore the crucial metabolism pathways related to the accumulation of astaxanthin in H. pluvialis, stressed under 7.5g·L^-1 NaCl, KEGG pathways enrichment analysis was performed for the CK and HS group. 51 differential metabolic pathways were enriched and annotated, among which 18 pathways were statistically significant (P<0.05) (Figure 5.). The results showed that the metabolites involved in key metabolic pathways were regulated and modified under the nitrogen-deficient stress. Furthermore, the significant pathways mainly focused on carbohydrate metabolism, amino acid metabolism and transport, glycerol and derivatives metabolism, nucleotide and derivative metabolism and energy metabolism.

To further reveal the dominant role in the key metabolic pathways related to the astaxanthin accumulation in H. pluvialis, under nitrogen-deficient stress, KEGG pathways enrichment analysis was performed for the CK and LN group. 51 differential metabolic pathways were enriched and annotated, among which 10 pathways were statistically significant (P<0.05) (Figure 6.). Moreover, the significant pathways mainly focused on amino acid metabolism and transport, nucleotide and derivative metabolism and energy metabolism.

3.5. Metabolimic Network Construction and Analysis

To exhibit the metabolic characteristics related to the astaxanthin accumulation in H. pluvialis under the high salinity and nitrogen-deficient conditions respectively, the metabolic networks were constructed based on the crucial pathways and metabolites (Figure 7 and Figure 8.).

3.6. Major Carbon Metabolism and the Changes in Identified Metabolites Related to Astaxanthin Accumulation in H. pluvialis

In aims to reveal and compare to the characteristics of the changes in identified metabolites involved in the major carbon metabolism related to the astaxanthin accumulation, under the different environmental stress conditions, the global profiling was exhibited (Figure 9.). The major carbon metabolism contained the carbohydrate and derivatives metabolism, glycerol and derivatives metabolism, amino acid metabolism and energy metabolism. The series of metabolic processes provided not only the precursors and intermediates related to the astaxanthin biosynthesis, but also the energy required in the metabolism scheme. As shown in Figure 9., the contents of the metabolites in LN group were almost lower than those both in the CK and HS groups, while which in HS group were substantially at the higher or close levels contrast to the CK group. However, the distinct patterns of several metabolites were exhibited according to the gragh. Specially, the content of cis-Aconitate was higher in the LN group than which both in the CK and HS groups. Otherwise, the levels of L-Glutamine in the three groups were almost identical. Moreover, the concentrations of Glycerol3P and DHAP were close to each other in the LN and HS group respectively, while both lower than the CK group. Additionally, the levels of L-Aspartate and GABA increased in order of CK, LN and HS group.

4. Discussion

4.1. Metabolic Characteristics of Astaxanthin Accumulation in H. pluvialis under the High Salinity Stress

When the fluctuation of the environmental conditions occurs, the physiological and metabolic processes of microalgae would change obviously. In response to the external environmental stress better, the energy flow and substances synthesis pathways in the cells would be adjusted accordingly. Many organic molecules such as carbohydrates, lipids and pigment synthesized in microalgae have been reported to cope with the variable environmental stress factors. Induced by the environmental stress, the biosynthesis and accumulation of astaxanthin which was the most important substance against stress in immotile cysts of H. pluvialis occurred gradually [29,30,31].

In the present study, the contents of astaxanthin in H. pluvialis stressed by the high salinity reached up to 7.51mg·L^-1, higher than the CK group significantly(P<0.05). It was highlighted that most of the astaxanthin in H. pluvialis did not exist in the free state, but the form of mono- and diesters conmbined with fatty acid (FA) after dehydration. Therefore, the synthesis and accumulation of astaxanthin and FA showed a obvious synchronization [32,33]. This view was proved by the result that upregulation of “glycerophospholipid metabolism” and “glycerolipid metabolism” pathways along with increasing of astaxanthin in the HS group, in the present study. Furthermore, sn-Glycerol-3-phosphate, one of the crucial intermediates of glycerophospholipid metabolism, is a key compound for lipid synthesis, providing the carbon skeleton to formate astaxanthin and lipids in microalgae cells [34]. In this study, the relative contents of sn-glycerol 3-phosphate in H. pluvialis were only 36.54% in the HS group of the control group, which suggested the biosynthesis process of astaxanthin was quite active to withstand the high salinity stress condition. In this research, activation of amino acid metabolism is beneficial to astaxanthin accumulation, because amino acid metabolism process can provide the necessary precursors and intermediates for astaxanthin biosynthesis [35]. It was reported that the contents of amino acid, protein and astaxanthin in H. pluvialis increased simultaneously under sodium chloride stress condition [36]. However, the metabolic pathways associated with different amino acids showed an obvious complexity in the present study. Specifically, “glycine, serine and threonine metabolism”, “valine, leucine and isoleucine biosynthesis” and “tyrosine metabolism” were activated in H. pluvialis treated by high salinity stress, leading to significant increase of glycine, L-serine, L-valine and L-tyrosine. This was consistent with the result of previous research [37]. In contrast, “arginine and proline metabolism” was downregulated under the same condition, accompanied with a drop of L-proline. The results implied that the different amino acids maybe play the viarable roles in astaxanthin biosynthesis. Moreover, simultaneous accumulation of GABA and astaxanthin in H. pluvialis stressed by the high salinity was observed in the present study, which was demonstrated in several previous studies [38,39]. In microalgae, GABA metabolism plays an important role in the organic acid and amino acid biosynthesis, and GABA is considered as a temporary nitrogen storage under stressful conditions [40]. Therefore, rising GABA would provide more intermediates for astaxanthin accumulation and substance reserves to cope with high salinity stress. As one precursor involved into GABA biosynthesis, glutamate would be transformed into GABA through GAD pathway [36]. In this research, the significantly opposite trends between GABA and glutamate were observed, indicated that the biotransformation of glutamate to GABA was quite active.

Carbohydrates are not only the important energy providers for microalgae, but also the raw materials for many physiological and metabolic process. It was demonstrated that the carbohydrates depletion in microalgae might be promoted because of the secondary metabolites accumulation, containing lipids and astaxanthin, stressed by the abiotic factors [37]. Yet, the level of sucrose in H. pluviali, the only carbohydrate detected in the present study, rised in tandem with astaxanthin under high salinity stress. Nevertheless, although “fructose and mannose metabolism” and “starch and sucrose metabolism” were enriched significantly between the CK and HS group, their involvement in astaxanthin accumulation were limited according to Figure 7. These inconsistent results suggested that the carbohydrate metabolism would be closely related to the strain of microalgal and the specific state of stress, and the raw materials required for astaxanthin synthesis were probably variable in different situations. Maltose would be turned into the carbon skeleton for the biosynthesis of astaxanthin and FA [18]. Although maltose and sucrose are isomers of each other, the interconversion of them in microalgae has not been found. In this research, increasing sucrose content indicated that it has not been consumed appreciably, implying that the biotransformation between sucrose and maltose would not occur in H. pluvialis. Additionally, inositol phosphate is an important signal molecule involved in several physiological metabolic processes, while inositol is the product of dephosphorylation of inositol phosphate. The relative content of inositol in the HS group was 1.33-fold of the control, hinting that the inositol phosphate signal transduction system inhibited under high salinity stress.

4.2. Metabolic Characteristics of Astaxanthin Accumulation in H. pluvialis under the Nitrogen-deficient Stress

According to the previous studies, nitrogen-deficient culture was an effective approach to facilitate the astaxanthin accumulation in H. pluvialis. However, variation of substance metabolism and flow occurred under the nitrogen-deficient stress condition, compared with the normal status [35,37,39]. Furthermore, to cope with the nitrogen constraint stress, resource allocation in cells would be adjusted accordingly. In particular, the processes of carbon metabolism and nitrogen metabolism were reintegrated, not only ensured the nitrogen balance in the cells, but the biosynthesis of anti-stress substances.

Under the nitrogen-deficient condition, proteins, one of the intracellular nitrogen storage, serve as a crucial role in the modulation of nitrogen metabolism, the levels of intracellular amino acids and protein reduced significantly, and nitrogen was perhaps recycled into amino acids synthesis in response to the resource constraints [40]. The proteins related to the catabolism would be re-distributed after the fluctuation of carbon skeletons during the biosynthesis of organic compounds [41]. In this study, “glycine, serine and threonine metabolism” and “arginine biosynthesis” were regulated significantly under the nitrogen-deficient condition. Subsequently, the contents of most amino acids including L-proline, glycine, L-glutamate, L-methionine, L-serine, L-leucine, L-phenylalanine, L-Tyrosine, arginine and L-valine in the LN group decreased significantly, which implied that quite a lot of amino acids in H. pluvialis were converted to endogenous nitrogen sources in response to nitrogen limitation. In particular, L-glutamate which is the crucial precursor in GAD pathway and serves as the hub compound in the carbon and nitrogen metabolism, specifically in nitrogen assimilation in microalgae [34]. In this research, the low level of L-glutamate in H. pluvialis demonstrated it has been heavily consumed during the readjustment between carbon metabolism and nitrogen metabolism, coping with the nitrogen-deficient stress. However, significant upregulation of L-aspartate was observed in the LN group, 2.42 fold of the control. The results indicated that the amino acids related to “alanine, aspartate and glutamate metabolism” had the different effects on astaxanthin accumulation and stress response in H. pluvialis. Although “lysine biosynthesis” pathway was enriched significantly between the CK and LN group, the content of lysine in the LN group was insignificant compared with the CK group. The results implied that the effects of lysine and “lysine biosynthesis” on astaxanthin accumulation in H. pluvialis were not obvious.

As a crucial procedure in the respiratory metabolism, the TCA cycle is considered as not only the hub of carbon-nitrogen metabolism but also the primary energy provider to life activities [42]. The intermediates involved in TCA cycle would provide the carbon skeletons components for the biosynthesis of many vital compounds, including astaxanthin and lipids [43]. In the present research, the contents of the pivotal compounds involved in TCA cycle, including succinate, fumarate and malate, showed a sharp decline in the LN group. The results suggested that TCA cycle was fully activated under the nitrogen constraint stress condition. Furthermore, malate from TCA cycle could enter into the FA biosynthesis pathway directly through NADP-dependent malic enzyme [44]. And then most of the astaxanthin in H. pluvialis was esterified with FA to generate the astaxanthin monoesters and diesters [33]. Therefore, the decline of malate and increase of astaxanthin were observed simultaneously in this study, which demonstrated that the synchronous accumulation of FA and astaxanthin were promoted by the consumption of malate. Furthermore, in order to provide the intermediates and precursors for the recycle of certain amino acids, the upregulation of TCA cycle could be triggered following the massive degradation of protein and amino acid [45]. The view were confirmed by the phenomenon, that the synchronization between sharp decline of most amino acids and significant promotion of TCA cycle in H. pluvialis under the nitrogen constraint stress, observed in the present study.

Yet, a series of small molecule organics such as glycolic acid, glyceric acid and formic acid, were supplied by “glyoxylate and dicarboxylate metabolism”, with participation in stress resistance and other metabolic pathways. Thus the obvious upregulation of “glyoxylate and dicarboxylate metabolism” pathway was observed in H. pluvialis under the nitrogen-deficient stress condition. It was reported that the endogenous nitrogen in H. pluvialis would be utilized for crucial amino acids and NAD(P)H production, which supplied the astaxanthin biosynthesis depended on NAD(P)H [26]. In the present study, the level of NADP in the LN group was only 54.21% of the control, which indicated that the biosynthesis of astaxanthin in H. pluvialis was activated effectively, leading to a sharp decline of NADP consumed heavily in the process.

4.3. Comparative Analysis of the major carbon metabolism in H. pluvialis under the different stress conditions

According to Figure 9., the patterns of the metabolic features involved in the major carbon metabolism related to the astaxanthin accumulation in H. pluvialis, under the high salinity and nitrogen-deficient conditions, varied wildly. Specially, the contents of the metabolites in LN group were almost lower than those both in the CK and HS groups, which indicated that the overall metabolic level stressed by nitrogen-deficient was inhibited, and the astaxanthin was accumulated to response to the stress. This was similar to the previous research [46]. In contrast, the metabolites in HS group were substantially at the higher or close levels contrast to the CK group, which demonstrated the promotion effects of the high salinity on the carbon metabolism related to the astaxanthin accumulation. However, the former study has been shown the reactive oxygen species (ROS) level would increased, accompanied by the astaxanthin accumulation [47].

In fact, the content of astaxanthin accumulated in H. pluvialis, meanwhile the level of oxidative stress rised induced by the high salinity unilaterally. The resent research indicated that the transcript levels of the dxs, lcy, bkt and chy genes related to the astaxanthin biosynthesis in H. pluvialis increased, but the ROS content decreased synchronously under the comprehensive induced conditions containing the high salinity, high-light and endogenous GABA [48]. The results implied that we could use the composite measures to neutralize the side effects of high salinity on the microalgae life actitivies in the astaxanthin production process.

5. Conclusions

In the present study, the differential metabolites and metabolic pathways under the high salinity and nitrogen-deficient stress conditions were highlighted respectively. Furthermore, the metabolic networks related to the astaxanthin accumulation were constucted based on the crucial pathways and metabolites. Moreover, the major carbon metabolism in H. pluvialis under the high salinity and nitrogen deficient stress were presented, which exhibited the characteristics in the identified metabolites changes under the different stress conditions.