3.1. Biomass and lipid productivity
Biomass and lipid productivity was calculated at the late log phase for all the runs based on their growth curve analysis. As shown in Fig. 1, the biomass productivity was affected by different nitrogen and phosphorus concentrations. Based on the results, the highest biomass productivity (166.3 ± 5.3 mg/l.d) was obtained at the concentration of 550 and 55 mg/l nitrogen and phosphorus, respectively, which was not significantly different from the biomass productivity at the concentration of 1000 mg/l nitrogen and 100 mg/l phosphorus. Results demonstrated that the maximum biomass productivity achieved at optimum ratio of N/P (around 10) while the highest content of lipid was observed at limited concentration of nitrogen or phosphorus (N/P > 10 or N/P < 3.7) because under this condition, cellular metabolism (Biomass production) moves to fatty acid biosynthesis (Miller et al., 2010).
The lipid content of Chlorella in BG-11 medium was about 14%, which has increased in all concentrations of nitrogen and phosphorus (Fig. 2a). The maximum lipid content (40.2 ± 1.2%) was obtained in nitrogen and phosphorus concentrations of 100 and 100 mg/l, respectively. According to the Fig. 2a, in low concentration of nitrogen (100 mg/l) the lipid content raised by increasing phosphorus concentration up to 100 mg/l due to accumulate consumed phosphorus as polyphosphate in microalgae. Polyphosphate plays an important role in the metabolism of protein, RNA and ATP, these products are a source of energy for lipid accumulation in microalgae under limited nitrogen conditions (Chu et al., 2014).
The results also show that in conditions where there is a limit concentration of phosphorus (10 mg/l), the lipid content increased by increasing the concentration of nitrogen to 1000 mg/l. It has been reported that under phosphorus limitation, lipid accumulation in chlorella is higher than carbohydrate accumulation (Liang et al., 2013). Phosphate limitation in medium could lead to the breakdown of phospholipids from cell wall into neutral lipids in order to obtain required phosphate for other metabolic processes. Therefore, the lipid content increases under this condition.
According to low biomass productivity under nitrogen or phosphorus deficient, parameter of lipid productivity should be analyzed to optimize lipid production. The highest lipid productivity (41.8 ± 1 mg/l.d) was obtained at nitrogen and phosphorus concentrations of 100 and 55 mg/l, respectively. As shown in Fig. 2b, in nitrogen and phosphorus concentrations of 1000 and 10 mg/l and 100 and 100 mg / l, respectively, despite the high lipid content, the lipid production efficiency is low because in these conditions (limitation of nitrogen and phosphorus concentration), the levels of ribose 1 and 5-carboxylic bisphosphate enzymes (involved in carbon dioxide consumption) and tryptophan synthesis are reduced. These enzymes affect photosynthetic function and reduce the growth of microalgae (Msanne et al. 2012).
3.2. Nitrogen removal
As shown in Fig. 3, the highest percentage of nitrogen removal is obtained at low concentrations of nitrogen (100 mg/l). The maximum nitrogen removal efficiency was obtained at concentration of 100 and 55 mg/L nitrogen and phosphorus, respectively, due to limitation value of nitrogen in the culture medium. It has been reported when the ratio of nitrogen to phosphorus is less than 3.7, there is a limit to the concentration of nitrogen, the microalgae consume more nitrogen and this consumption is used to keep the cell alive (Alketife et al., 2017). Because in this condition, the cell produces lipids instead of producing biomass to keep itself alive that is consistent with the results of biomass production in Fig. 1. It can be concluded that biomass yield is directly proportional to N concentration in culture media while lipid accumulation of microalgal cells is inversely proportional to N concentration. In general, microalgae need both of nitrogen and carbon source to grow and do photosynthesis. In this study, the only source of carbon was carbon dioxide in aeration and organic carbon (COD), which was not enough for the consumption of microalgae.
3.3 Effect of combined nitrogen and phosphorus concentration on the expression of accD gene of C. Vulgaris
ACCase catalyzes the first rate limiting step in the fatty-acid biosynthetic pathway through the formation of malonyl-CoA from acetyl-CoA. In this study, the expression of heteromeric unit of ACCase (accD gene) has been evaluated as a function of combined nitrogen and phosphorus concentration changes.
To confirm gene expression in the samples, the PCR product was applied to the electrophoresis gel, which is shown in Fig. 4. As shown in the gel electrophoresis, the specific bands for accD gene expression were not seen in negative control and NRT control. The NRT control sample was placed for all samples and in Fig. 4–10, only the NRT gel image of the control is shown. Also, due to the high color of some bands, the results show that the rate of gene expression in samples 2, 3, 4, 5 and 7 is higher than other samples.
As shown in Fig. 5, an increase in expression level of accD gene was observed in all experimental condition as compared to control (BG11). The highest levels of gene expression were obtained at nitrogen and phosphorus concentrations of 100 and 100 mg/l, respectively, which was approximately 8.5-fold increase in expression of accD gene in the control sample (culture in BG-11 medium). This value of gene expression is completely consistent with the optimal concentration of nitrogen and phosphorus that obtained in the previous section (percentage content and efficiency of lipid production). The results show that when there is a limit to the concentration of nitrogen or phosphorus in wastewater, the expression of the accD gene is increased and the microalgae tend to produce more lipids. Nitrogen deprivation conditions could lead to reduced cell division. Reduced cell division shifts the lipid biosynthetic pathways to synthesize more neutral lipids than synthesizing membrane lipids required for the cell wall formation (Singh et al., 2016). Subsequent accumulation of NADH due to the lower photosynthetic rate inhibits enzyme citrate synthase and prevents acetyl-CoA from entering into the TCA cycle. Elevated concentrations of acetyl-CoA activate acetyl-CoA carboxylase, this leads to enhanced lipid accumulation in microalgal cells (Li et al., 2015) Similarly, Fan et al., (2015) reported that both nitrogen and phosphorus limitation resulted in increasing the expression of accD gene in C. pyrenoidosa.
3.4. Predicted model using RSM
Gene expression and nitrogen removal were taken as the response to analyze the combined effect of nitrogen and phosphorus concentration on cultivation of C.vulgaris. With the results obtained from the CCD experimental design and data analysis, reduced quadratic polynomial equations were derived. Insignificant coefficients were omitted to achieve a more suitable regression model using back-ward elimination procedure. Eqs. (3) and (4) indicates the final reduced regression model in terms of coded factors for the nitrogen removal and gene expression, respectively:
Nitrogen removal = 32.26–6.54 * A + 2.87 * B + 14.89 * A2- 5.56 * B2 (3)
Gene expression = 4.54–1.07 * A- 0.45 * B – 2.4 * AB + 0.9740 * A2 (4)
Based on previous studies, R2 should be a minimum of 0.8 to fit a proper model (Joglekar et al., 1987). In both of these equations, R2 were approximately 0.96; meaning that the models were fit. In addition, the value of the adjusted determination coefficient for both of the responses (adjusted R2 = 0.94) suggested that about 6% of the total variation could not be explained by these models.
From the Table 3 it can be concluded that both factors were significant parameters on response nitrogen removal variation, while the nutrient concentration with most notable impact on nitrogen removal was nitrogen concentration (p < 0.0001). F-Value of 52.52 and p-value less than 0.0001 stated that the reduced quadratic model was valid. In addition, there was only a 0.01% chance that a model F-value could have occurred due to the
noise. Moreover, ANOVA analysis showed that the lack of fit of dates was not significant, revealing that the indicated model terms were significant for the considered response (Muhammad et al. 2013).
Table 3
ANOVA analysis for quadratic model of nitrogen removal response
Source | Sum of Squares | df | Mean Square | F-value | p-value | |
Model | 911.71 | 4 | 227.93 | 52.52 | < 0.0001 | significant |
A-A | 249.62 | 1 | 249.62 | 57.52 | < 0.0001 | |
B-B | 49.31 | 1 | 49.31 | 11.36 | 0.0098 | |
A² | 612.74 | 1 | 612.74 | 141.20 | < 0.0001 | |
B² | 85.23 | 1 | 85.23 | 19.64 | 0.0022 | |
Residual | 34.72 | 8 | 4.34 | | | |
Lack of Fit | 28.80 | 4 | 7.20 | 4.87 | 0.0771 | not significant |
Pure Error | 5.91 | 4 | 1.48 | | | |
Cor Total | 946.43 | 12 | | | | |
The ANOVA regression for gene expression response is shown in Table 4. Based on data, the parameter influence on gene expression response could be arranged as AB < A < A2 < B. Hence, nitrogen concentration had the most considerable effect on the gene expression.
Table 4
ANOVA analysis for quadratic model of gene expression response
Source | Sum of Squares | df | Mean Square | F-value | p-value | |
Model | 34.13 | 4 | 8.53 | 50.11 | < 0.0001 | significant |
A-A | 6.81 | 1 | 6.81 | 39.97 | 0.0002 | |
B-B | 1.22 | 1 | 1.22 | 7.14 | 0.0283 | |
AB | 23.04 | 1 | 23.04 | 135.32 | < 0.0001 | |
A² | 3.07 | 1 | 3.07 | 18.00 | 0.0028 | |
Residual | 1.36 | 8 | 0.1703 | | | |
Lack of Fit | 1.12 | 4 | 0.2810 | 4.72 | 0.0810 | not significant |
Pure Error | 0.2381 | 4 | 0.0595 | | | |
Cor Total | 35.49 | 12 | | | | |
The adequacy of these models was evaluated using diagnostic plot such as a plot of predicted versus actual values. According to Fig. 6a and b, the predicted values obtained from the presented models match well with the real values. Figure 6c and d demonstrates the three dimensional surface response as a function of the nitrogen and phosphorus concentration. It also indicates the major effects and the interaction between the investigated parameters.
According to the Fig. 6c, in all nitrogen concentrations, the nitrogen removal efficiency rose with increasing the concentration of phosphorus to 55 mg/l. It can be seen that the effect of this increase in low concentrations of nitrogen is greater than high concentrations of nitrogen.
As shown in Fig. 6d, at high concentrations of nitrogen (1000 mg/L), the expression of gene decreased by increasing the concentration of phosphorus, which indicates that the limitation of phosphorus concentration in the cultivation lead the microalgae to increase the expression of accD. Also, the results show that at high concentrations of phosphorus (100 mg/l), the relative expression of the gene dropped by rising the concentration of nitrogen.
3.5. Optimization model for simultaneously nitrogen removal and accD gene expression
Numerical optimization technique was applied to investigate the optimal conditions of the process. The goal was to maximize nitrogen removal and gene expression simultaneously. Predictably, maximum value of N removal and relative accD gene expression would be reached in low nitrogen concentration at 100 mg/L and close to the maximum of phosphorus concentration at 95.7 mg/L. In this optimum values, the maximum predicted N removal and level of gene expression was 51.64% and 8.34, respectively, and the desirability was 0.909.