Next Article in Journal
Bioactive Capnosanes and Cembranes from the Soft Coral Klyxum flaccidum
Next Article in Special Issue
Effects of Temperature, Growth Media, and Photoperiod on Growth and Toxin Production of Azadinium spinosum
Previous Article in Journal
Microalgal Enzymes with Biotechnological Applications
Previous Article in Special Issue
Biorefining and the Functional Properties of Proteins from Lipid and Pigment Extract Residue of Chlorella pyrenoidosa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Light Intensity and Nitrogen Concentration Impact on the Biomass and Phycoerythrin Production by Porphyridium purpureum

by
Juan Eduardo Sosa-Hernández
1,
Laura Isabel Rodas-Zuluaga
1,
Carlos Castillo-Zacarías
1,
Magdalena Rostro-Alanís
1,
Reynaldo de la Cruz
1,
Danay Carrillo-Nieves
2,
Carmen Salinas-Salazar
1,
Claudio Fuentes Grunewald
3,
Carole A. Llewellyn
3,
Eugenia J. Olguín
4,
Robert W. Lovitt
5,6,
Hafiz M. N. Iqbal
1,* and
Roberto Parra-Saldívar
1,*
1
Tecnologico de Monterrey, Escuela de Ingenieria y Ciencias, Campus Monterrey, Ave. Eugenio Garza Sada 2501, Monterrey, Nuevo Leon 64849, Mexico
2
Tecnologico de Monterrey, Escuela de Ingenieria y Ciencias, Av. General Ramón Corona 2514, Nuevo México, Zapopan C.P. 45138, Jalisco, Mexico
3
Department of Biosciences, Singleton Park, Swansea University, Swansea, Wales SA2 8PP, UK
4
Environmental Biotechnology Group, Institute of Ecology(INECOL), Carretera Antigua a Coatepec No. 351, Xalapa, Veracruz 91070, Mexico
5
College of Engineering, Swansea University, Swansea SA2 8PP, UK
6
Membranology Ltd., Unit D5 Rainbow Business Centre, Swansea SA7 9FP, UK
*
Authors to whom correspondence should be addressed.
Mar. Drugs 2019, 17(8), 460; https://doi.org/10.3390/md17080460
Submission received: 21 June 2019 / Revised: 31 July 2019 / Accepted: 1 August 2019 / Published: 7 August 2019
(This article belongs to the Special Issue Bioactive Compounds Derived from Marine Microalgae)

Abstract

:
Several factors have the potential to influence microalgae growth. In the present study, nitrogen concentration and light intensity were evaluated in order to obtain high biomass production and high phycoerythrin accumulation from Porphyridium purpureum. The range of nitrogen concentrations evaluated in the culture medium was 0.075–0.450 g L−1 and light intensities ranged between 30 and 100 μmol m−2 s−1. Surprisingly, low nitrogen concentration and high light intensity resulted in high biomass yield and phycoerythrin accumulation. Thus, the best biomass productivity (0.386 g L−1 d−1) and biomass yield (5.403 g L−1) were achieved with NaNO3 at 0.075 g L−1 and 100 μmol m−2 s−1. In addition, phycoerythrin production was improved to obtain a concentration of 14.66 mg L−1 (2.71 mg g−1 of phycoerythrin over dry weight). The results of the present study indicate that it is possible to significantly improve biomass and pigment production in Porphyridium purpureum by limiting nitrogen concentration and light intensity.

Graphical Abstract

1. Introduction

Phycoerythrin (PE) is a bioactive pigment mainly found in Rhodophyta, but it can also be recovered from Cyanobacteria [1]. This water-soluble pigment is a bioactive molecule with interesting properties such as anti-inflammatory, immunosuppressive, antitumor, antioxidant, antidiabetic, and antihypertensive effects [2,3,4,5]. In addition, PE moiety can be exploited as a fluorescent probe useful in cell tracking and imaging assays such as flow cytometry, fluorescent immunoassays, and immunophenotyping. [6,7]. As well as its application in the biotechnological and medicine areas, PE can be used in the food and cosmetic industries [8,9]. Commercial PE from microalgae has been produced by Spirulina, Rhodella, and Porphyridium [10,11,12]. The market price of this purified pigment is around $200 USD mg−1 (Sigma-Aldrich, product number: 52412-1MG-F), this being why PE is regarded as a molecule of high value [13]. Porphyridium purpureum has the ability to grow in different culture systems and locations [14], coupled with the capacity to produce and accumulate a large amount of PE, which production can be enhanced when the microalgae are cultured under stress conditions [11]. The effect of different physical and nutritional factors on the culture conditions of different microalgae groups has been investigated for several purposes, including the production of biomass, pigments, oil, and other metabolites [11,15,16,17,18]. Light intensity is one of the most critical factors that influence microalgal pigment composition. This is extremely important in those microalgae who have a complex harvest light system such as Porphiridium because, for instance, PE is the major accessory light-harvesting compound in Rhodophyta [19,20]. On the other hand, nitrogen is a fundamental constituent for protein synthesis, cell division and it also is required for the synthesis of photosynthetic pigments [21,22]. The effect of these two factors on Porphyridium strains has been reported in some studies [23,24,25]. However, in most of the cases, the variables of light intensity and nitrogen concentrations have been evaluated separately. Furthermore, most of the investigations regarding light influence on the production of PE focused only in the evaluation of different wavelengths [26,27] or illumination periods [28], but scarce information is related to light intensity.
Complementing biotechnology work, mathematical models have been widely used to explore microalgae growth. The application of computational power allows us to understand, predict and optimize cellular processes present in microalgae, without the need for experimental tests [29]. Mathematical models have been used in the cultivation of P. purpureum, however, because of the complexity of these models, it is difficult to apply them [30,31]. The Gompertz model allows us to predict parameters such as maximum specific growth rate (µmax), latency phase (λ), duplication time of microalgae cells (G), productivity and because of its simplicity, it can be easily applied; however, this model has not been applied in the growth of P. purpureum [32,33,34].
Due to the high value and remarkable diverse bioactivities of PE, in-depth exploration about the abiotic parameters management of microalgae cultures is needed in order to improve its production. In the present study, the effect of different light intensities and nitrogen concentrations were evaluated with the aim to obtain high biomass production and to accumulate a high PE concentration in P. purpureum. Besides the identification of the best levels of the abiotic factors evaluated, the correlation between the response variables was aimed as well.

2. Results

2.1. Experimental Results and Mathematical Model Fitness

The Gompertz model was used with the same conditions from the experimental results. In this manner, it was possible to get parameters that are not directly quantifiable. To fit this model, all the data were used in hours instead of days to have a better resolution in the model, then the time parameter was returned to days. The model was used to compare all experiments performed. The mathematical modeling to determine the biological parameters of P. purpureum culture, as well as the productivity and determination coefficient of the linear regression model obtained between the experimental and simulated results, are shown in Table 1.
The model fitness parameter used was the determination coefficient R2, which are over 90% for almost all curves. The two experiments under 90% were in the light intensity of 65 µmol m−2 s−1. Both experiments had 0.075 and 0.450 g L−1 concentration of NaNO3 and were in the range of 85% which is reasonably expected to be the biological variability of growth [35].
According to the productivity of the systems that were subjected to one way ANOVA analysis and paired two-sample mean with t-Test analysis, the values of productivity were higher (p = 0.105) for 65 µmol m−2 s−1 and significantly (p = 0.047) higher 100 µmol m−2 s−1 compared to the value registered at 30 µmol m−2 s−1, which were 350 and 386 mg L−1 d−1, respectively.

2.2. Influence of Light Intensity and NaNO3 Concentration on the Production of P. purpureum Biomass

The kinetics of the biomass production from the culture of P. purpureum under different light intensities and NaNO3 concentration is presented in Figure 1A–C. The lowest biomass production was observed at 30 µmol m−2 s−1 regardless of the variation of nitrogen source (Figure 1A and Table 1). Under these culture conditions (30 µmol m−2 s−1), in the tests at 0.075 and 0.225 g L−1 of NaNO3, a delay in the latency phase of five days was observed, while at 0.45 g L−1 of NaNO3, it was reduced to three days (also, according to the fit parameters shown in Table 1). This could be associated with the fact that the medium had a higher concentration of this nutrient. In addition, the maximum specific growth rate and productivity values were higher for the kinetics with nitrogen source of 0.075 and 0.225 g L−1 of NaNO3 and the duplication time of microalgae cells was shorter, compared with the condition of 0.45 g L−1 of NaNO3 (see Table 1).
Under a higher light intensity of 65 µmol m−2 s−1, a lag phase of six days and exponential phase of eight days of exponential were observed for the three concentrations of NaNO3 evaluated. Overall, the culture reached a maximum biomass concentration at day 14, with 4.896, 2.914, and 2.248 g L−1 for 0.075, 0225, and 0.450 g L−1 of NaNO3, respectively (see Figure 1).
Furthermore, high values of biomass were obtained under an increased light intensity at 100 µmol m−2 s−1 (Figure 1C). The maximum biomass production of 5.403 g L−1 was obtained at 14 days when the concentration of NaNO3 was fixed at 0.075 g L−1. The maximal production for the treatments with NaNO3 concentrations of 0.225 and 0.450 g L−1 were 2.579 and 2.598 g L−1, respectively at 10 days of culture. Culture conditions fixed at 100 µmol m−2 s−1 of light intensity and 0.075 g L−1 of NaNO3 concentration resulted in the highest productivity and maximum specific growth rate (μmax) with values of 0.386 g L−1 d−1 and 0.015 d−1, respectively. Similar productivity of 0.349 g L−1 d−1 was obtained under the concentration of NaNO3 0.075 g L−1 by applying a light intensity of 65 µmol m−2 s−1. The productivity under a light intensity of 30 µmol m−2 s−1 was lower (and independent) than that of the NaNO3 concentration. The low levels of NaNO3 resulted in a biomass of 0.950 g L−1, while the biomass obtained under the high level of NaNO3 was 0.340 g L−1. In other words, 2.79 times more biomass was obtained maintaining of NaNO3 at a low concentration.

2.3. Influence of Light Intensity and NaNO3 Concentration on the Production of PE by P. purpureum

The kinetics of PE production by P. purpureum under different light intensities and NaNO3 concentrations is presented in Figure 2. The experiments under a light intensity of 30 µmol m−2 s−1 are shown in Figure 2A. Under these culture conditions, the three experiments showed similar lag phases (6 days). A maximum PE concentration of 11.24 mg L−1 (5.70 mg g−1 PE/DW) was observed at 0.225 g L−1 of NaNO3 at the 13th day of culture. However, high PE content was detected since the ninth day (10.37 mg L−1). Similarly, 10.58 mg L−1 PE was determined at 0.075 g L−1 of NaNO3, but obtained at the 10th day. The PE content was lower (maximum value of 4.29 mg L−1) when the concentration of NaNO3 was 0.450 g L−1. The results of PE concentration under a light intensity of 65 µmol m−2 s−1 are shown in Figure 2B. Six days of a lag phase were observed when the concentration of NaNO3 was set at 0.075 g L−1. This experiment showed a long exponential phase of seven days, reaching a maximum PE concentration of 15.18 mg L−1 (3.10 mg g−1 PE/DW). Lower concentrations of PE (4.74 and 4.71 mg L−1) were observed when the concentration of NaNO3 was set at 0.225 and 0.450 g L−1, respectively, at day 10. The experiments under a light intensity of 100 µmol m−2 s−1 are shown in Figure 2C. Maximum PE production of 14.65 mg L−1 (2.71 mg g−1 PE/DW) was obtained at day 14, with NaNO3 at 0.075 g L−1. When the concentration of NaNO3 was 0.450 g L−1, it resulted in a maximum PE production of 10.73 mg L−1 (4.13 mg g−1 PE/DW) after 10 days. A significantly lower value of maximum PE production (4.35 mg L−1) was observed with NaNO3 concentration of 0.225 g L−1. The best value of PE productivity (1.08 mg L−1 d−1) was observed under a light intensity of 65 µmol m−2 s−1 and 0.075 g L−1 of NaNO3. These same conditions were used to obtain the highest PE accumulation.
In addition, the nitrogen and phosphorus ratio was changed between the three experiments. The N:P ratio calculated for each experiment is represented by 11:1, 33:1 and 66:1 respectively with 0.075, 0.225 and 0.450 g L−1 concentrations of NaNO3. The best PE production is at 11:1 N:P ratio in combination with high light intensity. The higher N:P ratios did not promote PE production, except the 66:1 N:P ratio at the high light intensity with a PE production of 10.73 mg L−1 (4.13 mg g−1 PE/DW) and the 33:1 ratio at the low light intensity with a PE production of 11.25 mg L−1 (5.70 mg g−1 PE/DW).
Similar to the biomass behavior, the highest values of PE content were observed under light intensities of 65 and 100 µmol m−2 s−1 (15.18 and 14.66 mg L−1, respectively) that correspond to 3.10 and 2.71 mg g−1 PE/DW when the experiments were maintained at 0.075 g L−1 of NaNO3. PE production under 65 µmol m−2 s−1 was reduced by 69% and 71% when the concentration of NaNO3 was increased to 0.225 and 0.450 mg L−1. Under the experiment at 100 µmol m−2 s−1, a similar effect occurred reducing PE in 70% and 27%, respectively. Experiments under a light intensity of 30 µmol m−2 s−1 did not affect PE when the nitrate concentration increased from 0.075 to 0.225 g L−1.
The best values of PE production were observed in the experiments maintained at 0.075 g L−1. An increment of the light intensities of 65 and 100 µmol m−2 s−1 increased PE production by 43% and 30%, respectively. However, a contrary effect was observed when the concentration of NaNO3 was fixed at 0.225 g L−1, causing a PE reduction of 57% and 61%, respectively, due to the light intensities of 65 and 100 µmol m−2 s−1. Experiments with NaNO3 at 0.450 g L−1 did not show any effect when light intensity increased from 30 to 65 µmol m−2 s−1. However, when cultures were exposed to 100 µmol m−2 s−1, PE production increased by 149%.

3. Discussion

3.1. Mathematical Model Parameter Analysis

The model fitness to the experimental data was high for almost all curves (over 90%). Even the two that were about 80% was good enough considering the biological variability of the system. The Gompertz model helped to accurately estimate the productivity and latency parameters of the P. purpureum culture. The mentioned parameters were used to distinguish the kinematics of culture growth and differentiate the actual potential to select the best conditions for biomass production.
Regarding the selection of the best set of culture conditions to maximize biomass, higher productions were found in 65 and 100 µmol m−2 s−1, compared to 30 µmol m−2 s−1. For instance, the best NaNO3 concentration was 0.075 g L−1, with no significant difference (p = 0.084) between light intensities of 65 and 100 µmol m−2 s−1. Although, it is possible that the higher light intensity can be better for larger volume scale with a longer light path.

3.2. Effect of Light Intensity on Biomass and PE Production in P. purpureum

The high impact of light intensity on biomass production by P. purpureum was evident in this study. Light intensities of 65 and 100 µmol m−2 s−1 resulted in increases of 253% and 390% of biomass production, respectively, compared to that obtained with 30 µmol m−2 s−1. However, surprisingly, this increased biomass production was observed when NaNO3 was lower than 0.075 g L−1. Experiments carried out at 0.225 and 0.450 g L−1 of NaNO3 also showed lower increments on biomass production (up to 187%) at high light intensities. In the literature, light intensities ranging from 18.75 to 50 µmol m−2 s−1 resulted in biomass production in P. purpureum of up to 0.950 g L−1 [36,37,38]. On the other hand, higher levels of light intensity (100 to 140 µmol m−2 s−1) yielded biomass production above 2 g L−1 [26,39]. A study with Porphyridium marinum was able to produce 2 g L−1 with a light intensity of 150 µmol m−2 s−1 to later change the condition to 70 µmol m−2 s−1 [40]. The trends in the present study and in the literature mentioned above indicate that high biomass production occurs at a high light intensity. A similar trend has been observed in P. purpureum SCS-2 [41] and several other microalgae strains, such as Chlorella vulgaris [42], Parachlorella sp. [19], and Spirulina platensis [43]. It is important to compare these results with other similar biomass yields found in other studies with P. purpureum (shown in Table 2). The yield of biomass production by weight is within the range of the presented reports. Our hypothesis for this effect is that low nutrient availability in the culture medium may induce polysaccharide, starch and/or carbohydrate accumulation, [27,44,45,46], and coupled with low light intensity, which promotes the accumulation of lipids 19.4% dry weight [37] and/or dark light cycle [47]. By using the mechanisms of photo-adaption and photo-protection with different light wavelengths sustained 17% dry weight [27]. Guihéneuf and Stengel in 2015 performed a study with several culture conditions, where the nitrogen was maintained replete, limited and depleted. The results of their work show that in nitrogen depleted cultures, P. purpureum accumulates carbohydrates (57% dry weight), and interestingly found that the biomass productivity reached 0.15 g L−1 d−1.
Experiments under low light intensity (30 µmol m−2 s−1) when the nitrate concentration increased from 0.075 to 0.225 g L−1 did not affect PE production. However, 60% of PE decrease was observed as a result of an increment of NaNO3 concentration at 0.450 mg L−1. A similar effect was observed by Sánchez–Saavedra since they reported a PE production of 44.45 µg L−1 under a light intensity of 200 µmol m−2 s−1 [25]. In contrast, the evaluation at a light intensity of 50 µmol m−2 s−1 caused a reduction of 60% in the PE concentration. For instance, we can hypothesize that at high light intensity more PE molecules are produced in order to harvest the maximum amount of photons reaching a limit.
It is reported in the literature that an extremely high light intensity causes low biomass yield, due to a phenomenon of photoinhibition in microalgae especially at the beginning of the curve when the culture is non-dense enough [48]. The research done by Guihéneuf and Stengel in 2015, where they explored different conditions to growth P. purpureum the best light intensity was at 100 µmol m−2 s−1 compared to results with 30 and 200 µmol m−2 s−1, where the biomass and other compounds production were decreased substantially. In accordance with Guihéneuf and Stengel, in the present study, the light intensity of 100 µmol m−2 s−1 was the highest level evaluated and it resulted in the best biomass yield. Therefore, these results suggest that higher light levels, on the range of evaluation above of 100 µmol m−2 s−1 could improve biomass yield as a result of the direct relationship between the increase of biomass produced and the light intensity, but at an intensity below 200 µmol m−2 s−1 similar to the report of Gargouch et al. (2018). This effect may be due to the evaluated strain of P. purpureum which showed high photosynthetic efficiency to transform this energy into high biomass production [49]. The photosynthetic efficiency of P. purpureum should be related to the accumulation of PE since this molecule is responsible for the light absorption [50]. Furthermore, studies on the macroalgae Gracilaria domingensis have shown this trend where high irradiance stimulated the accumulation of proteins [51].
A possibility to select a lower light intensity could arise due to the reduction of cost by using less energy for light in an upscale process. This option could be tested by the cost-benefit analysis of a larger-volume culture set-up against the possible decrease in the penetration of light through it. Moreover, other similar parameters should be tested, such as mixing, aeration, and temperature as is recommended in a previous study [27].

3.3. Effect of Nitrogen Concentration on Biomass and Phycoerythrin from P. purpureum

At low light intensities, the biomass production was independent of nitrate concentrations used in this study. At a high light intensity, the systems became sensitive to higher concentrations of nitrate, reducing the biomass yield. In contrast, biomass production was reduced by 50% when the concentration of NaNO3 was increased in 0.225 and 0.450 g L−1 under a light intensity of 65 and 100 µmol m−2 s−1. Therefore, the higher biomass production resulted when the lower value of NaNO3 was evaluated (0.075 g L−1) only when high light intensities were applied. A similar effect is reported in a study done by Kavitha et al. (2016), where they evaluated 1 and 2 g L−1 of NaNO3 for production of P. purpureum biomass. High levels of biomass (3.4–9.4 g L−1) from P. purpureum were reported by Kathiresan et al. (2007), who suggested that low concentrations of chloride, nitrate, and phosphate did not produce any significant effect on the biomass production. However, an increment of NaNO3 above 1 g L−1 provokes a negative impact on biomass yield as well as in PE production. It is important to notice, that the main effect caused by a disproportionate increase in the concentration of nitrogen source is the possible limitation of other nutrients, such as phosphorus, trace metals and vitamins [46].
In the present study, increments of the nitrogen concentration were evaluated to obtain higher biomass production and higher PE accumulation. However, under the culture conditions evaluated, the best results were obtained at a low nitrogen concentration (0.075 g L−1). Phosphorous is also an important element for microalgae growth, which in conjunction with nitrogen allows biomass production; still, microalgae cannot absorb nitrogen (N) without the presence of phosphorus (P) in the culture, or vice versa, because a proper N:P ratio is essential for their growth [52]. The phosphorous concentration was maintained at 0.005 g L−1 (original concentration of the culture medium). The increase of nitrogen concentration changed the N:P ratio. Therefore, results suggest that nitrogen concentration could be increased, but if the N:P ratio of the original culture medium (15:1) is maintained. As studied previously by Razaghi in 2014, the results presented here show similar behavior for N:P ratio change impact the P. purpureum growth and PE production. In contrast, the report of Gargouch et al. (2018), where the strain P. marinum is used to produce PE (a maximum of 40 mg/g of PE/DW). The culture medium was optimized to produce high biomass and PE and the formulation contains 3.4 g L−1 of NaNO3 and 0 g L−1 K2HPO4. This reflects a complete different behavior to P. purpureum regarding the N:P ratio in terms of this study and others [23,24,37,41,46,52].

4. Materials and Methods

4.1. Strain, Medium and Culture Conditions

The marine strain Porphyridium purpureum UTEX LB 2757 was obtained from the University of Texas, United States of America. The strain was maintained in a mixture (1:1) of Bold 1NV media (UTEX, Weimar, TX, USA) and Erdshreiber media (UTEX, Weimar, TX, USA) in 500 mL flasks containing 300 mL of culture medium. The inoculum (1:10 v/v) was cultured at 20 °C in a flask containing 300 mL of culture medium mixed by bubbling sterile air at a flow rate of 25 mL min−1. The culture was continuously illuminated with white fluorescent tubes with an average intensity of 65 µmol m−2 s−1 for 9 days. Then, aliquots of the culture were taken to seed 500 mL Erlenmeyer flasks filled with 270 mL of medium and 30 mL of inoculum coupled with sealed-in tubes to pass a stream of air. Two key factors, the light intensity and nitrogen limitation influencing the production of PE pigment were evaluated. The light intensity was designated as high, medium or low by adjusting different levels of the intensity (100, 65, and 30 µmol m−2 s−1) with a combined light with wavelengths of green and blue (520–525 and 465–468 nm). The effects of nitrogen concentration were evaluated at three levels of NaNO3 concentration (0.450, 0.225 and 0.075 g L−1), keeping the F/2 Phosphorous content fixed. Culture nutrient composition was based on the medium F/2 (UTEX, TX, USA), temperature and air flow rate were constant at 20 °C and 25 mL min−1. Then, samples were collected each day within 15 days including the days 0 to 3, 6 to 10 and 13 to 14 for the analysis of total biomass and PE produced by microalgae. The microalgae biomass was washed three times with MilliQ water. The samples were centrifuged for 15 min at 4500 rpm. The supernatant was discarded and the pellet was resuspended in MilliQ water. Then, the sample was filtered through a Whatman regenerated cellulose membrane filter, pore size: 0.45 µm, diameter 4.7 cm followed by drying in an oven at 105 °C for 24 h.

4.2. Mathematical Model and Kinetic Parameters Calculation to Determine Growth Parameters

In this study, the Gompertz model (Equation (1)) was used to determine maximum specific growth rate (μmax, Equation (2)), duration of the latency phase (Equation (3)), duplication time of microalgae cells (Equation (4)) and biomass productivity (Equation (5)) in experiments at different conditions of light intensity and NaNO3 concentration.
Y = a e ( e ( b c t ) )
where:
  • Y = log ( N / N 0 ) , N 0 . Initial biomass ( t = 0 ) in the sample (mg L−1), biomass counts for each time.
  • t : Time (h).
  • a : Constant determined from the experimental data. Represents the maximum population of biomass counts when time increases indefinitely.
  • b : Constant determined from the experimental data.
  • c : Constant determined from the experimental data.
Maximum specific growth rate:
μmax = a × c × (h − 1)
Duration of the latency phase:
λ = (b − 1)/c × (h)
Generation or duplication time:
G = ln (2)/μmax (h)
The biological parameters determined by the Gompertz model were estimated using the catalog of non-linear regression functions offered by Minitab Inc. 17.1.0 software (LEAD Technologies, Inc., Charlotte, North Carolina, USA) with the Levenberg—Marquardt algorithm.
The determination coefficient (R2) of linear regression model obtained between the experimental and simulated results allow us to estimate the precision of the experimental data with the model obtained, as well as the quality of the model to replicate the results. Values above 95% conclude that the model manages to replicate 95% of the experiments were carried out with very good quality.
Biomass productivity ( P x ) was calculated (Equation (5)) as follows:
P x = X m a x X i t m a x
where X m a x and X i are maximum and initial biomass concentrations, respectively, and t m a x is the time corresponding to X m a x .

4.3. PE Extraction Methodology

A sample of 2 mL from each P. purpureum culture was centrifuged at 4500 rpm for 15 min and the supernatant was discarded. The sediment was mixed with 2 mL of distilled water and was vortexed until homogenization. The suspension was frozen at −20 °C for one hour, and then it was thawed and sonicated in an ice bath for 10 min. The sample was vortexed until homogenization and repeat steps until completing five cycles. After that, the samples were centrifuged at 4500 rpm for 15 min. The pink supernatant was recovered and measured in a spectrophotometer at 564, 592 and 455 nm. The obtained data were substituted in the Beer and Eshel (1985) equation (Equation (6)) to calculate the concentration of PE:
P E   ( m g   m L 1 ) = [ ( O D 564 n m O D 592 n m ) ( O D 455 n m O D 592 n m )   0.2 ] × 0.12

5. Conclusions

The comparison between experimental and modeled results shows that the Gompertz model can be used correctly for modeling the Phorphyridium purpureum culture within the presented ranges of light intensities and NaNO3. This model has the potential to be a simple and powerful tool for the prediction of biomass and PE production. Although, corroboration is needed to include wider ranges of culture conditions.
Light intensity and nitrogen content resulted in strong impact factors to the biomass production of P. purpureum and its PE accumulation. The values obtained in the present study are promising since the production is within those reported in the literature, despite a low nitrogen concentration. In addition, this study is the basis for forthcoming optimization and scaling-up studies. The strategy to follow in future experiments to obtain high PE production by P. purpureum is to first promote biomass growth by using low nitrogen and high light intensity. After this nitrogen can be increased and a high light intensity can be sustained.

Author Contributions

J.E.S.-H., L.I.R.-Z., C.C.-Z., H.M.N.I. and R.P.-S. conceptualized the study; L.I.R.-Z. and C.C.-Z. perform the experiments and collected the data; J.E.S.-H., L.I.R.-Z., C.C.-Z., M.R.-A., R.d.l.C. and D.C.-N. analyzed the results; J.E.S.-H. and D.C.-N. constructed the mathematical model and evaluated parameters; J.E.S.-H., L.I.R.-Z., C.C.-Z., R.d.l.C., D.C.-N. drafted the manuscript; C.S.-S., C.F.G., C.A.L., E.J.O., R.W.L., H.M.N.I. and R.P.-S. revised the manuscript; C.F.G., R.W.L., E.J.O., H.M.N.I. and R.P.-S. coordinated the study. All authors read and approved the final manuscript.

Funding

The authors would like to thank the Sustainable Applied Biotechnology Group from Tecnologico de Monterrey, Mexico, the Institute of Ecology, Mexico and Membranology, Co, Swansea, UK. This work was totally funded by CONACYT (Mexican Council of Science and Technology) and INNOVATE UK project Phycopigments (grant number 268792).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hirose, Y.; Misawa, N.; Yonekawa, C.; Nagao, N.; Watanabe, M.; Ikeuchi, M.; Eki, T. Characterization of the genuine type 2 chromatic acclimation in the two Geminocystis cyanobacteria. DNA Res. 2017, 24, 387–396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sosa-Hernández, J.E.; Romero-Castillo, K.D.; Parra-Arroyo, L.; Aguilar-Aguila-Isaías, M.A.; García-Reyes, I.E.; Ahmed, I.; Parra-Saldivar, R.; Bilal, M.; Iqbal, H. Mexican Microalgae Biodiversity and State-Of-The-Art Extraction Strategies to Meet Sustainable Circular Economy Challenges: High-Value Compounds and Their Applied Perspectives. Mar. Drugs 2019, 17, 174. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, D.; Nishizawa, M.; Shimizu, Y.; Saeki, H. Anti-inflammatory effects of dulse (Palmaria palmata) resulting from the simultaneous water-extraction of phycobiliproteins and chlorophyll a. Food Res. Int. 2017, 100, 514–521. [Google Scholar] [CrossRef] [PubMed]
  4. Li, P.; Ying, J.; Chang, Q.; Zhu, W.; Yang, G.; Xu, T.; Yi, H.; Pan, R.; Zhang, E.; Zeng, X.; et al. Effects of phycoerythrin from Gracilaria lemaneiformis in proliferation and apoptosis of SW480 cells. Oncol. Rep. 2016, 36, 3536–3544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Patel, S.N.; Sonani, R.R.; Jakharia, K.; Bhastana, B.; Patel, H.M.; Chaubey, M.G.; Singh, N.K.; Madamwar, D. Antioxidant activity and associated structural attributes of Halomicronema phycoerythrin. Int. J. Biol. Macromol. 2018, 111, 359–369. [Google Scholar] [CrossRef]
  6. Al-Essa, M.K.; Melzer, S.; Tarnok, A.; Hadidi, K.A.; El-Khateeb, M. Fast RBC loading by fluorescent antibodies and nuclei staining dye and their potential bioanalytical applications. Z. Für Nat. C 2018, 73, 95–105. [Google Scholar] [CrossRef]
  7. Wu, J.; Chen, H.; Zhao, J.; Jiang, P. Fusion proteins of streptavidin and allophycocyanin alpha subunit for immunofluorescence assay. Biochem. Eng. J. 2017, 125, 97–103. [Google Scholar] [CrossRef]
  8. Chen, X.; Wu, M.; Yang, Q.; Wang, S. Preparation, characterization of food grade phycobiliproteins from Porphyra haitanensis and the application in liposome-meat system. LWT 2017, 77, 464–474. [Google Scholar] [CrossRef]
  9. Intrarapuk, A.; Awakairt, S.; Thammapalerd, N.; Mahannop, P.; Pattanawong, U.; Suppasiri, T. Comparison between R-phycocyanin-labeled and R-phycoerythrin-labeled monoclonal antibody (Mab) probes for the detection of Entamoeba histolytica trophozoites. Southeast Asian J. Trop. Med. Public Health 2001, 32, 165–171. [Google Scholar]
  10. Chandra, R.; Parra, R.; MN Iqbal, H. Phycobiliproteins: A novel green tool from marine origin blue-green algae and red algae. Protein Pept. Lett. 2017, 24, 118–125. [Google Scholar] [CrossRef]
  11. Koch, M.; Kniggendorf, A.K.; Meinhardt-Wollweber, M.; Roth, B. In vivo determination of carotenoid resonance excitation profiles of Chlorella vulgaris, Haematococcus pluvialis, and Porphyridium purpureum. J. Raman Spectrosc. 2018, 49, 404–411. [Google Scholar] [CrossRef]
  12. Rodrigues, R.D.P.; de Castro, F.C.; de Santiago-Aguiar, R.S.; Rocha, M.V.P. Ultrasound-assisted extraction of phycobiliproteins from Spirulina (Arthrospira) platensis using protic ionic liquids as solvent. Algal Res. 2018, 31, 454–462. [Google Scholar] [CrossRef]
  13. Kersen, P.; Paalme, T.; Pajusalu, L.; Martin, G. Biotechnological applications of the red alga Furcellaria lumbricalis and its cultivation potential in the Baltic Sea. Bot. Mar. 2017, 60, 207–218. [Google Scholar]
  14. Fuentes-Grünewald, C.; Bayliss, C.; Zanain, M.; Pooley, C.; Scolamacchia, M.; Silkina, A. Evaluation of batch and semi-continuous culture of Porphyridium purpureum in a photobioreactor in high latitudes using Fourier Transform Infrared spectroscopy for monitoring biomass composition and metabolites production. Bioresour. Technol. 2015, 189, 357–363. [Google Scholar] [CrossRef]
  15. Chang, J.; Le, K.; Song, X.; Jiao, K.; Zeng, X.; Ling, X.; Shi, T.; Tang, X.; Sun, Y.; Lin, L. Scale-up cultivation enhanced arachidonic acid accumulation by red microalgae Porphyridium purpureum. Bioprocess Biosyst. Eng. 2017, 40, 1763–1773. [Google Scholar]
  16. Pancha, I.; Chokshi, K.; George, B.; Ghosh, T.; Paliwal, C.; Maurya, R.; Mishra, S. Nitrogen stress triggered biochemical and morphological changes in the microalgae Scenedesmus sp. CCNM 1077. Bioresour. Technol. 2014, 156, 146–154. [Google Scholar] [CrossRef]
  17. Contreras-Angulo, J.R.; Mata, T.M.; Cuellar-Bermudez, S.P.; Caetano, N.S.; Chandra, R.; Garcia-Perez, J.S.; Muylaert, K.; Parra-Saldivar, R. Symbiotic Co-Culture of Scenedesmus sp. and Azospirillum brasilense on N-Deficient Media with Biomass Production for Biofuels. Sustainability 2019, 11, 707. [Google Scholar] [CrossRef]
  18. Fuentes-Grünewald, C.; Garcés, E.; Alacid, E.; Sampedro, N.; Rossi, S.; Camp, J. Improvement of lipid production in the marine strains Alexandrium minutum and Heterosigma akashiwo by utilizing abiotic parameters. J. Ind. Microbiol. Biotechnol. 2012, 39, 207–216. [Google Scholar]
  19. Heo, J.; Shin, D.S.; Cho, K.; Cho, D.H.; Lee, Y.J.; Kim, H.S. Indigenous microalga Parachlorella sp. JD-076 as a potential source for lutein production: Optimization of lutein productivity via regulation of light intensity and carbon source. Algal Res. 2018, 33, 1–7. [Google Scholar] [CrossRef]
  20. Sciandra, A.; Lazzara, L.; Claustre, H.; Babin, M. Responses of growth rate, pigment composition and optical properties of Cryptomonas sp. to light and nitrogen stresses. Mar. Ecol. Prog. Ser. 2000, 201, 107–120. [Google Scholar] [CrossRef] [Green Version]
  21. Gigova, L.G.; Ivanova, N.J. Microalgae respond differently to nitrogen availability during culturing. J. Biosci. 2015, 40, 365–374. [Google Scholar] [CrossRef]
  22. Singh, P.; Guldhe, A.; Kumari, S.; Rawat, I.; Bux, F. Investigation of combined effect of nitrogen, phosphorus and iron on lipid productivity of microalgae Ankistrodesmus falcatus KJ671624 using response surface methodology. Biochem. Eng. J. 2015, 94, 22–29. [Google Scholar] [CrossRef]
  23. Guihéneuf, F.; Stengel, D.B. Towards the biorefinery concept: Interaction of light, temperature and nitrogen for optimizing the co-production of high-value compounds in Porphyridium purpureum. Algal Res. 2015, 10, 152–163. [Google Scholar] [CrossRef]
  24. Kathiresan, S.; Sarada, R.; Bhattacharya, S.; Ravishankar, G.A. Culture media optimization for growth and phycoerythrin production from Porphyridium purpureum. Biotechnol. Bioeng. 2007, 96, 456–463. [Google Scholar] [CrossRef]
  25. Del Pilar Sánchez-Saavedra, M.; Castro-Ochoa, F.Y.; Nava-Ruiz, V.M.; Ruiz-Güereca, D.A.; Villagómez-Aranda, A.L.; Siqueiros-Vargas, F.; Molina-Cárdenas, C.A. Effects of nitrogen source and irradiance on Porphyridium cruentum. J. Appl. Phycol. 2018, 30, 783–792. [Google Scholar] [CrossRef]
  26. Baer, S.; Heining, M.; Schwerna, P.; Buchholz, R.; Hübner, H. Optimization of spectral light quality for growth and product formation in different microalgae using a continuous photobioreactor. Algal Res. 2016, 14, 109–115. [Google Scholar] [CrossRef]
  27. Coward, T.; Fuentes-Grünewald, C.; Silkina, A.; Oatley-Radcliffe, D.L.; Llewellyn, G.; Lovitt, R.W. Utilising light-emitting diodes of specific narrow wavelengths for the optimization and co-production of multiple high-value compounds in Porphyridium purpureum. Bioresour. Technol. 2016, 221, 607–615. [Google Scholar] [CrossRef] [Green Version]
  28. Munier, M.; Jubeau, S.; Wijaya, A.; Morancais, M.; Dumay, J.; Marchal, L.; Jaouen, P.; Fleurence, J. Physicochemical factors affecting the stability of two pigments: R-phycoerythrin of Grateloupia turuturu and B-phycoerythrin of Porphyridium cruentum. Food Chem. 2014, 150, 400–407. [Google Scholar] [CrossRef]
  29. Darvehei, P.; Bahri, P.A.; Moheimani, N.R. Model development for the growth of microalgae: A review. Renew. Sustain. Energy Rev. 2018, 97, 233–258. [Google Scholar] [CrossRef]
  30. Fleck, P.; Posten, C. Modelling of growth and product formation of Porphyridium purpureum under defined conditions. IFAC Proc. Vol. 2004, 37, 281–286. [Google Scholar] [CrossRef]
  31. Fleck-Schneider, P.; Lehr, F.; Posten, C. Modelling of growth and product formation of Porphyridium purpureum. J. Biotechnol. 2007, 132, 134–141. [Google Scholar] [CrossRef] [PubMed]
  32. Nayak, B.K.; Das, D. Improvement of carbon dioxide biofixation in a photobioreactor using Anabaena sp. PCC 7120. Process Biochem. 2013, 48, 1126–1132. [Google Scholar] [CrossRef]
  33. Zhen, G.; Lu, X.; Kobayashi, T.; Kumar, G.; Xu, K. Anaerobic co-digestion on improving methane production from mixed microalgae (Scenedesmus sp., Chlorella sp.) and food waste: Kinetic modeling and synergistic impact evaluation. Chem. Eng. J. 2016, 299, 332–341. [Google Scholar] [CrossRef]
  34. Dammak, M.; Hadrich, B.; Barkallah, M.; Hentati, F.; Hlima, H.B.; Pichon, C.; Denis, M.; Fendri, I.; Michaud, P.; Abdelkafi, S. Modelling Tetraselmis sp. growth-kinetics and optimizing bioactive-compound production through environmental conditions. Bioresour. Technol. 2018, 249, 510–518. [Google Scholar] [CrossRef] [PubMed]
  35. Buchanan, R.L.; Whiting, R.C.; Damert, W.C. When is simple good enough: A comparison of the Gompertz, Baranyi, and three-phase linear models for fitting bacterial growth curves. Food Microbiol. 1997, 14, 313–326. [Google Scholar] [CrossRef]
  36. Assunção, M.F.; Varejão, J.M.; Santos, L.M. Nutritional characterization of the microalga Ruttnera lamellosa compared to Porphyridium purpureum. Algal Res. 2017, 26, 8–14. [Google Scholar] [CrossRef]
  37. Kavitha, M.D.; Kathiresan, S.; Bhattacharya, S.; Sarada, R. Culture media optimization of Porphyridium purpureum: Production potential of biomass, total lipids, arachidonic and eicosapentaenoic acid. J. Food Sci. Technol. 2016, 53, 2270–2278. [Google Scholar] [CrossRef]
  38. Schulze, C.; Wetzel, M.; Reinhardt, J.; Schmidt, M.; Felten, L.; Mundt, S. Screening of microalgae for primary metabolites including β-glucans and the influence of nitrate starvation and irradiance on β-glucan production. J. Appl. Phycol. 2016, 28, 2719–2725. [Google Scholar] [CrossRef]
  39. Sandefur, H.N.; Asgharpour, M.; Mariott, J.; Gottberg, E.; Vaden, J.; Matlock, M.; Hestekin, J. Recovery of nutrients from swine wastewater using ultrafiltration: Applications for microalgae cultivation in photobioreactors. Ecol. Eng. 2016, 94, 75–81. [Google Scholar] [CrossRef]
  40. Gargouch, N.; Karkouch, I.; Elleuch, J.; Elkahoui, S.; Michaud, P.; Abdelkafi, S.; Laroche, C.; Fendri, I. Enhanced B-phycoerythrin production by the red microalga Porphyridium marinum: A powerful agent in industrial applications. Int. J. Biol. Macromol. 2018, 120, 2106–2114. [Google Scholar] [CrossRef]
  41. Li, T.; Xu, J.; Wu, H.; Jiang, P.; Chen, Z.; Xiang, W. Growth and Biochemical Composition of Porphyridium purpureum SCS-02 under Different Nitrogen Concentrations. Mar. Drugs 2019, 17, 124. [Google Scholar] [CrossRef] [PubMed]
  42. Bazdar, E.; Roshandel, R.; Yaghmaei, S.; Mardanpour, M.M. The effect of different light intensities and light/dark regimes on the performance of photosynthetic microalgae microbial fuel cell. Bioresour. Technol. 2018, 261, 350–360. [Google Scholar] [CrossRef] [PubMed]
  43. Qiang, H.; Richmond, A. Productivity and photosynthetic efficiency of Spirulina platensis as affected by light intensity, algal density and rate of mixing in a flat plate photobioreactor. J. Appl. Phycol. 1996, 8, 139–145. [Google Scholar] [CrossRef]
  44. Köst, H.P.; Senser, M.; Wanner, G. Effect of nitrate and sulphate starvation on Porphyridium cruentum cells. Z. Pflanzenphysiol. 1984, 113, 231–249. [Google Scholar] [CrossRef]
  45. Adda, M.; Merchuk, J.C.; Arad, S.M. Effect of nitrate on growth and production of cell-wall polysaccharide by the unicellular red alga Porphyridium. Biomass 1986, 10, 131–140. [Google Scholar] [CrossRef]
  46. Razaghi, A.; Godhe, A.; Albers, E. Effects of nitrogen on growth and carbohydrate formation in Porphyridium cruentum. Open Life Sci. 2014, 9, 156–162. [Google Scholar] [CrossRef]
  47. Oh, S.H.; Han, J.G.; Kim, Y.; Ha, J.H.; Kim, S.S.; Jeong, M.H.; Jeong, H.S.; Kim, N.Y.; Cho, J.S.; Yoon, W.B.; et al. Lipid production in Porphyridium cruentum grown under different culture conditions. J. Biosci. Bioeng. 2009, 108, 429–434. [Google Scholar] [CrossRef]
  48. George, B.; Pancha, I.; Desai, C.; Chokshi, K.; Paliwal, C.; Ghosh, T.; Mishra, S. Effects of different media composition, light intensity and photoperiod on morphology and physiology of freshwater microalgae Ankistrodesmus falcatus–A potential strain for bio-fuel production. Bioresour. Technol. 2014, 171, 367–374. [Google Scholar] [CrossRef]
  49. Ji, F.; Hao, R.; Liu, Y.; Li, G.; Zhou, Y.; Dong, R. Isolation of a novel microalgae strain Desmodesmus sp. and optimization of environmental factors for its biomass production. Bioresour. Technol. 2013, 148, 249–254. [Google Scholar] [CrossRef]
  50. Wiltbank, L.B.; Kehoe, D.M. Two cyanobacterial photoreceptors regulate photosynthetic light harvesting by sensing teal, green, yellow, and red light. MBio 2016, 7, e02130-15. [Google Scholar] [CrossRef]
  51. Ramlov, F.; Souza, J.; Faria, A.V.; Maraschin, M.; Horta, P.A.; Yokoya, N.S. Growth and accumulation of carotenoids and nitrogen compounds in Gracilaria domingensis (Kütz.) Sonder ex Dickie (Gracilariales, Rhodophyta) cultured under different irradiance and nutrient levels. Rev. Bras. Farmacogn. 2011, 21, 255–261. [Google Scholar] [CrossRef]
  52. Beuckels, A.; Smolders, E.; Muylaert, K. Nitrogen availability influences phosphorus removal in microalgae-based wastewater treatment. Water Res. 2015, 77, 98–106. [Google Scholar] [CrossRef]
Figure 1. Biomass concentration of P. purpureum (g L−1) during cultivation under three different light intensities (A) 30, (B) 65 and (C) 100 (µmol m−2 s−1); and three different NaNO3 concentrations 0.075 (●), 0.225 (■) and 0.450 (♦) (g L−1). The error is less than 5% or in general, so it is in the range of the symbol size.
Figure 1. Biomass concentration of P. purpureum (g L−1) during cultivation under three different light intensities (A) 30, (B) 65 and (C) 100 (µmol m−2 s−1); and three different NaNO3 concentrations 0.075 (●), 0.225 (■) and 0.450 (♦) (g L−1). The error is less than 5% or in general, so it is in the range of the symbol size.
Marinedrugs 17 00460 g001
Figure 2. Phycoerythrin concentration (mg L−1) during cultivation under three different light intensities (A) 30, (B) 65 and (C) 100 (µmol m−2 s−1); and three different NaNO3 concentrations 0.075 (●), 0.225 (■) and 0.450 (♦) (g L−1).
Figure 2. Phycoerythrin concentration (mg L−1) during cultivation under three different light intensities (A) 30, (B) 65 and (C) 100 (µmol m−2 s−1); and three different NaNO3 concentrations 0.075 (●), 0.225 (■) and 0.450 (♦) (g L−1).
Marinedrugs 17 00460 g002
Table 1. Kinetic parameters of P. purpureum cultivation under different light intensities and NaNO3 concentrations (maximum specific growth rate (µmax), latency phase (λ), duplication time of microalgae cells (G), determination coefficient (R2) and productivity).
Table 1. Kinetic parameters of P. purpureum cultivation under different light intensities and NaNO3 concentrations (maximum specific growth rate (µmax), latency phase (λ), duplication time of microalgae cells (G), determination coefficient (R2) and productivity).
Light Intensity (µmol m−2 s−1)
3065100
NaNO3 (g L−1)NaNO3 (g L−1)NaNO3 (g L−1)
0.0750.2250.450.0750.2250.450.0750.2250.45
Gompertz model constantsa0.740.882.051.331.741.041.591.032.27
b9.314.261.511.761.163.111.421.281.27
c0.070.090.0070.0130.0060.0170.0090.0050.008
µmax (d−1) 0.050.080.0150.0180.0110.0180.0150.0050.018
λ (h) 122.3139.170.6456.125.72124.745.135634.31
G (h) 13.748.2246.9638.2362.2739.3746.6445.7839.2
R2 (%) 93.8394.7195.4288.691.4986.294.8291.4795.59
Productivity (mg L−1 d−1) 173.2151.8138.6349.8208.2160.6386258259.9
Table 2. Comparison with recent reports of P. purpureum biomass growth and yields.
Table 2. Comparison with recent reports of P. purpureum biomass growth and yields.
DW (g L−1)N (g L−1)Nitrogen SourceY (N/DW)Time (d)Light Intensity (μmol m−2 s−1)Reference
5.541.78KNO30.04516350[41]
3.41NaNO30.04910100[23]
9.120.97NaNO30.0181418.85[24]
3.0–5.4030.075–0.45NaNO30.0023–0.02514100This report

Share and Cite

MDPI and ACS Style

Sosa-Hernández, J.E.; Rodas-Zuluaga, L.I.; Castillo-Zacarías, C.; Rostro-Alanís, M.; de la Cruz, R.; Carrillo-Nieves, D.; Salinas-Salazar, C.; Fuentes Grunewald, C.; Llewellyn, C.A.; Olguín, E.J.; et al. Light Intensity and Nitrogen Concentration Impact on the Biomass and Phycoerythrin Production by Porphyridium purpureum. Mar. Drugs 2019, 17, 460. https://doi.org/10.3390/md17080460

AMA Style

Sosa-Hernández JE, Rodas-Zuluaga LI, Castillo-Zacarías C, Rostro-Alanís M, de la Cruz R, Carrillo-Nieves D, Salinas-Salazar C, Fuentes Grunewald C, Llewellyn CA, Olguín EJ, et al. Light Intensity and Nitrogen Concentration Impact on the Biomass and Phycoerythrin Production by Porphyridium purpureum. Marine Drugs. 2019; 17(8):460. https://doi.org/10.3390/md17080460

Chicago/Turabian Style

Sosa-Hernández, Juan Eduardo, Laura Isabel Rodas-Zuluaga, Carlos Castillo-Zacarías, Magdalena Rostro-Alanís, Reynaldo de la Cruz, Danay Carrillo-Nieves, Carmen Salinas-Salazar, Claudio Fuentes Grunewald, Carole A. Llewellyn, Eugenia J. Olguín, and et al. 2019. "Light Intensity and Nitrogen Concentration Impact on the Biomass and Phycoerythrin Production by Porphyridium purpureum" Marine Drugs 17, no. 8: 460. https://doi.org/10.3390/md17080460

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