A framework for accelerated phototrophic bioprocess development: integration of parallelized microscale cultivation, laboratory automation and Kriging-assisted experimental design
Biotechnology for Biofuels volume 10, Article number: 26 (2017)
Even though microalgae-derived biodiesel has regained interest within the last decade, industrial production is still challenging for economic reasons. Besides reactor design, as well as value chain and strain engineering, laborious and slow early-stage parameter optimization represents a major drawback.
The present study introduces a framework for the accelerated development of phototrophic bioprocesses. A state-of-the-art micro-photobioreactor supported by a liquid-handling robot for automated medium preparation and product quantification was used. To take full advantage of the technology’s experimental capacity, Kriging-assisted experimental design was integrated to enable highly efficient execution of screening applications. The resulting platform was used for medium optimization of a lipid production process using Chlorella vulgaris toward maximum volumetric productivity. Within only four experimental rounds, lipid production was increased approximately threefold to 212 ± 11 mg L−1 d−1. Besides nitrogen availability as a key parameter, magnesium, calcium and various trace elements were shown to be of crucial importance. Here, synergistic multi-parameter interactions as revealed by the experimental design introduced significant further optimization potential.
The integration of parallelized microscale cultivation, laboratory automation and Kriging-assisted experimental design proved to be a fruitful tool for the accelerated development of phototrophic bioprocesses. By means of the proposed technology, the targeted optimization task was conducted in a very timely and material-efficient manner.
By virtue of significant advantages offered over agricultural crops [1–6], microalgae are generally accepted as promising feedstock for bio-economy applications [7–9]. However, until now, their industrial exploitation remains mostly uneconomic, especially when lower-priced products like biofuels are targeted . Currently, the integrated utilization of biomass is intensively investigated as a promising concept to improve the overall efficiency in terms of cost and energy [9, 11, 12]. In this context, intracellular lipids represent a compound class of special interest as they can be either transesterificated to biodiesel  or boost the nutritional quality of algae for functional food applications .
Regarding phototrophic bioprocess development, early-stage strain and parameter screening are of crucial importance to the successful set-up of economic processes . Today, these aspects are typically studied by means of only marginally parallelized reactor systems like shake flasks, test tubes or even single-vessel reactors [15–18]. Consequently, experimental throughput is fairly limited rendering screening tasks rather laborious and highly time consuming. Only recently a strong demand for high throughput micro-photobioreactors has been identified, based on which some prototype systems have been developed [19–27]. To take full advantage of phototrophic microscale cultivation, supporting methodologies and technologies, such as simplified strain maintenance , high throughput analytics [29, 30], and automated processing [20, 31], are needed.
In the medium term, these initial developments and especially further progress in high throughput technology and laboratory automation will clearly boost the efficiency of phototrophic process development. Phototrophic processes are characterized by their intrinsic complexity induced by a high number of potentially interacting input variables. Hence, experimental capacities, i.e., mainly cultivation, will be always a crucial factor due to the trade-off between throughput and the necessary laboratory resources. Current micro-photobioreactors mainly rely on standardized microtiter plates. Thus, a further rise of cultivation capacity by intensifying parallelization would need to be based on scale out and inevitably be accompanied by increasing cost. Hence, an alternative strategy to focus cultivation activities on only the most informative experiments is the ultimate solution to tackle the omnipresent challenge of restricted experimental throughput.
One approach to achieve this efficiently, which is already well-established for microbial bioprocess development , is the use of Design of Experiments (DoE) to focus on experiments providing the highest information content in a targeted parameter space. Despite having been established during the early twentieth century , there is still ongoing research into this methodology . This approach is regarded to be particularly suitable to deal with the combinatorial explosion typically occurring when investigating multi-parameter relations . Moreover, DoE overcomes a critical limitation of “conventional” one-factor-at-a-time experiments, as such approaches often fail in locating global optima by not taking potentially synergistic or antagonistic interactions of input variables into account . Regarding bioprocess development, the most prominent application of DoE is the culture media optimization . During such tasks, the omnipresent interactions between single compounds render locating a global optimum by “conventional” experimental planning to chance.
In the above context, the current study aims at the combination of emerging technologies for parallelized microscale cultivation and analytics to phototrophic microorganisms with elaborate experimental design as has previously been fruitfully applied for heterotrophic systems by . Thereby, an integrated framework for the accelerated development of phototrophic bioprocesses is to be set up. Optimizing medium composition toward maximized lipid productivity of the unicellular microalga Chlorella vulgaris was chosen as a model process for the above purpose.
All chemicals were purchased either from Sigma-Aldrich (Steinheim/Germany) or Roth (Karlsruhe/Germany) and were of analytical grade. The unicellular microalga C. vulgaris 211-11b , purchased from the Culture Collection of Algae at the University of Göttingen (Germany), was used throughout all cultivation experiments.
Cultivations were carried out in variations of an enriched Bold’s Basal Medium  prepared from stock solutions. The previously established reference medium [27, 28] was composed of chemicals as follows: 9.76 g L−1 2-(N-morpholino)ethanesulfonic acid (MES), 0.6 g L−1 K2HPO4, 1.4 g L−1 KH2PO4, 1.5 g L−1 NaNO3, 187.5 mg L−1 MgSO4·7 H2O, 6.25 mg L−1 NaCl, 125 mg L−1 CaCl2·2 H2O, 17.64 mg L−1 ZnSO4·7 H2O, 2.88 mg L−1 MnCl2·4 H2O, 2.4 mg L−1 Na2MoO4·2 H2O, 3.14 mg L−1 CuSO4·5 H2O, 0.94 mg L−1 CoSO4·7 H2O, 22.8 mg L−1 H3BO3, 9.96 mg L−1 FeSO4·7 H2O, 3.68 mg L−1 H2SO4, 100 mg L−1 Na2EDTA·2 H2O, 62 mg L−1 KOH and 100 mg L−1 penicillin-G sodium salt. The pH value was set to 6.5 with 5 M NaOH. Ultrapure water (type 1) was used for the preparation of all cultivation media.
During optimization experiments, the medium composition was varied according to the respective experimental plan by adjusting the applied volumes of the individual stock solutions. These media variants were prepared by a liquid-handling platform as previously described in literature. Medium preparation was carried out in a fully automated manner, while a surrounding laminar flow hood ensured sterile conditions [34, 40, 41]. Media were prepared at 2.5 mL scale in an MTP-R-48-B “Round Well Plate” (m2p-labs, Baesweiler/Germany) under continuous shaking at 500 rpm on an integrated Teleshake 95 (Inheco, Martinsried/Germany), while the minimum volume to be pipetted was set to 10 µL. Thereby, achieving sufficiently high accuracy (±0.3%) and precision (±0.3%) could be ensured . Subsequently, 950 µL of each medium was transferred to a well of an MTP-48-B “FlowerPlate®” (m2p-labs, Baesweiler/Germany) in which the cultivation took place (see “Main cultivation” section).
Strain maintenance and pre-cultivation
Chlorella vulgaris was maintained as glucose-adapted cryocultures. Preserved cells were re-adapted to light during phototrophic pre-cultivation in illuminated shake flasks. A detailed description of strain maintenance and pre-cultivation strategy can be obtained from . After 60 h of incubation, the cells were harvested by 5 min centrifugation at 3939×g and 4 °C in a Labofuge 400R (Heraeus Instruments, Hanau/Germany). The supernatant was discarded, and the pellet re-suspended in 0.9% (w v−1) NaCl to a biovolume of 2 µL mL−1 to generate the stock solution required for the inoculation of subsequent main cultivations.
Main cultivations were conducted in pre-sterilized, disposable 48-well MTP-48-B “FlowerPlates®” (m2p-labs, Baesweiler/Germany). Each well was filled with 950 µL of medium and inoculated to a biovolume of 0.1 µL mL−1 with 50 µL of the inoculation stock solution generated as described in “Strain maintenance and pre-cultivation” section. The plates were sealed using an F-R48-10 “perforated sealing foil for evaporation reduction” (m2p-labs, Baesweiler/Germany), pasted over with an F-GP-AB10 “gas-permeable seal” (m2p-labs, Baesweiler/Germany).
The microtiter plates were incubated using a micro-photobioreactor prototype. The system relies on bottom-side illumination with a set of blue and white LEDs and indirect temperature control via placement of the plates in a tempered incubation chamber. A detailed description and the schematic representation of the system are given in . The following cultivation conditions were applied: 25 °C, continuous shaking at 1200 rpm, 3 mm shaking diameter, 2.5% (v v−1) CO2, 200 µmol m−2 s−1 photon flux density (constant), and ≥85% relative humidity.
Optical density (OD) was acquired using 10-mm polystyrene semi-micro cuvettes (ratiolab, Dreieich/Germany) and an UV-1800 photometer (Shimadzu, Duisburg/Germany) at 750 nm, while desalted water served as a blank. If needed, samples were diluted to OD750 ≤ 0.3 using 0.9% (w v−1) NaCl solution to fit the linear range of the photometer.
The biovolume was measured taking advantage of a particle counter (MultiSizer 3, Beckman Coulter, Krefeld/Germany) using the “Coulter principle” . The device was equipped with a 30 µm capillary which had been calibrated using a suspension of 3 µm latex beads (Beckman Coulter, Krefeld/Germany) according to the manufacturer’s specification and was operated in volumetric control mode. Prior to measurement, cell suspensions were diluted to OD750 ≤ 0.025 in CASYton buffer (Omni Life Science, Bremen/Germany), and only particles in the range of 1.8–14 µm were analyzed.
The cell dry weight was determined by means of gravimetry. Culture liquid from two replicate wells of a microtiter plate was pooled to obtain sufficient sample amounts for the analysis. Cells were spun down in pre-dried and weighed 2-mL reaction tubes for 5 min at 16,060×g (Biofuge Pico, Heraeus Instruments, Hanau/Germany). The supernatants were discarded and the pellets freeze-dried in an LT-105 freeze dryer (Christ Gefriertrocknungsanlagen, Osterode am Harz/Germany) until attaining a constant weight. After acclimatization to room temperature in a desiccator, weighing was repeated, and the cell dry weight was derived from the resulting mass difference.
The intracellular accumulation of neutral lipids was quantitatively monitored by means of an automated high throughput Nile red staining assay as previously described in .
Cells were removed by filtration using 0.2 µm cellulose acetate syringe filters (DIA-Nielsen, Düren/Germany), and the cell-free supernatant was stored at −20 °C prior to analysis, if needed. Nitrate was quantified using the Spectroquant 1.09713.0002 nitrate test (Merck, Darmstadt/Germany) according to the manufacturer’s specifications, scaled down to one quarter of the recommended volume. Supernatants were pre-diluted with desalted water to fit the linear range of the assay, if needed. The measurements were conducted in UV semi-micro cuvettes (Brand, Wertheim/Germany) using an UV-1800 photometer (Shimadzu, Duisburg/Germany).
Acquisition of fatty acid fingerprints
Lyophilized biomass from cell dry weight determination (see Sect. 2.5) was in-situ transesterificated using acidic methanol [10% (w w−1) H2SO4], and the resulting fatty acid methyl esters were subsequently extracted with heptane. Semi-quantitative fingerprints were accessed by gas chromatography time-of-flight mass spectrometry of the extracts. A detailed description of the methodology can be obtained from .
Media composition was optimized with respect to lipid productivity using a Design of Experiments methodology. The applied optimization strategy was adopted from . Initially, fractional and full factorial experimental designs were applied for estimating single component effects and combinatorial interactions. Myers et al.  provide a good overview of these classical DoE methods.
Based on the initially collected data, the statistically more advanced concept of Kriging was applied for data analysis, visualization, and for designing further experiments with potentially improved lipid productivity. Kriging is an interpolation method that provides unbiased approximations of the underlying nonlinear functional relationships between media composition and lipid productivity with minimal prediction error. This method originates in geostatistics and has recently been adapted for optimizing biotechnology processes . Further mathematical details of the Kriging method can be found in the monograph of Cressie . The statistical analysis tools applied in this study are part of the open source Kriging toolkit “KriKit”, which can be freely downloaded at https://github.com/modsim/KriKit.
Given a Kriging model of the current dataset, further experiments were designed to maximize the Expected Improvement (EI). This experimental design strategy seeks a compromise between maximizing lipid productivity and reducing prediction uncertainty of the Kriging approximation in relevant regions of the parameter space . In a comparative study, EI has been found to outperform other sampling strategies in Kriging-based optimization .
In sequential optimization, new experiments are typically planned at maximal EI. Parallel experiments, as in the present study, are most efficiently planned by sampling from the EI distribution. In a non-deterministic sampling process, using the Markov Chain Monte Carlo (MCMC) method, new experiments are selected with probability proportional to their EI. Naturally, experiments with high EI are preferred over experiments with lower EI, which nonetheless have a reduced chance of being selected, while experiments with zero EI are strictly excluded. Freier et al. have demonstrated that MCMC sampling can significantly reduce the number of required experiments in process optimization . In the present study, the Delay Rejection Adaptive Metropolis algorithm  was applied with a chain length of 10,000 elements, of which the first 1000 are discarded (burn in phase of the MCMC method).
Results and discussion
Choice of relevant media components
The medium targeted for optimization incorporates 17 different components (see “Medium”) with phosphate salts counted as one compound due to their pH-dependent equilibrium. This number is too high to efficiently perform the experimental study with a manageable number of experiments, since a full factorial design with two concentration levels would result in 217 ≈ 130,000 experiments. In order to keep the number of components of interest, preselection was completed based on the literature information. Table 1 summarizes the known biological effects of the individual components. Penicillin-G concentration was kept constant under all conditions, and all trace elements were clustered to one single input variable as a similar effect on cultivation was expected. Sulfuric acid and potassium hydroxide had to be varied together with FeSO4 and Na2EDTA, respectively, as they were needed to keep the latter two components dissolved in their stock solutions. Thereby, the number of input variables was reduced by almost 50% from 17 to 9.
Starting with the nine remaining media components of interest, a full factorial design would require 29 ≈ 500 experiments. Making full use of 48-fold parallelized microtiter plate cultivation (see “Main cultivation” section), this leads to a total of 11 experimental runs, equivalent to 4 months of cultivation time. In 12-fold parallelized shake flasks, the experiments would even take 14 months. Yet, such time scales are clearly far from feasible, underlining the necessity to effectively reduce the experimental effort.
Fractional factorial designs allow the reduction of the number of experiments by estimating only single component effects and a subset of combinatorial effects . The chosen design (see Additional file 1 for both, design and corresponding measurement data) comprises 37 experiments, five of which represent the reference point using the enBBMref medium (see Additional file 2 for medium composition). Taking reference points into account allows for the investigation of measurement noise and normalization. The other experiments allowed for a statistical analysis of the effect of single components, as well as the interaction with magnesium ions. The interaction with this divalent metal ion was analyzed, as it is reported to be an effector of the acetyl-CoA carboxylase, an enzyme essential for lipid biosynthesis responsible for the initial step of carbon dioxide fixation to malonyl-CoA (see Table 1). An overview of the functionality of this enzyme complex and its regulation is given by Ohlrogge and Browse . Thus, any interactions with this input variable are of special interest with respect to product accumulation in the cells.
Figure 1a shows the resulting statistical analysis of the fractional factorial experiments. The green bars indicate the expected effect of varying the medium concentrations between their minimal and maximal values (see Additional file 1). The error bars indicate the uncertainty of the estimations. In the following, the main and combinatorial effects of the components are checked for significance using a t test with a significance level of p = 0.1. Using a lower significance level would increase the risk of false negatives, i.e., excluding relevant components from the remaining study. The diagram shows that an increase in the concentration of NaNO3 has a significant negative effect on lipid productivity. On the other hand, an increase in the trace element’s concentration results in a significant (p < 0.1, t test) productivity improvement. Furthermore, the analysis indicates a positive tendency with the increasing CaCl2 concentration and the lowering EDTA concentration. However, because of measurement noise and the comparably low number of experiments, the uncertainty of the estimation is relatively high and leads to no reliable statement about the effects of CaCl2 and EDTA, respectively. Similarly, this holds true for MgSO4, but here, the pairwise interaction with another component was additionally investigated. As shown in Additional file 3, a significant negative combinatorial effect was identified with the sodium salts, NaNO3 and NaCl.
For visual inspection of the negative combinatorial effect, a Kriging model was constructed based on the given data. The predicted functional relationship between MgSO4, NaNO3, and the lipid productivity is displayed in Fig. 1b. In case of low NaNO3 concentration, the interpolation reveals a positive correlation between an increase in MgSO4 and that of the performance indicator. With the increasing NaNO3 concentration, this positive effect is weakened.
In conclusion, significant effects of NaNO3 and the trace elements were identified, as well as positive tendencies of MgSO4 and CaCl2. Furthermore, the effect of MgSO4 appears to depend on the sodium salts, NaNO3 and NaCl. The remaining components have only low potential to affect the lipid productivity and were thus excluded from further analysis.
In order to verify the observed tendencies and to investigate potential pairwise or higher combinatorial effects, a full factorial design was constructed for the remaining five input variables: NaNO3, MgSO4, CaCl2, NaCl, and the clustered trace elements. This design again comprises five reference points and 32 experiments with minimal/maximal concentration (see Additional file 4 for the individual designs and the corresponding measurement data).
Figure 2a shows the updated statistical results after performing the full factorial design. The previously observed effects of NaNO3 and the trace elements were confirmed. The positive tendency of CaCl2 turned out to be significant, while the effect of NaCl remained insignificant. However, the interaction of MgSO4 with the sodium salts could be investigated in more detail. Figure 2b shows the opposing effect of MgSO4 dependent on NaNO3. This interaction leads to a non-distinguishable single component effect of MgSO4, as indicated in Fig. 2a. The analysis also revealed a negative interaction between CaCl2 and the trace elements, as indicated in the screening plot shown in Additional file 5.
Locating optimal medium composition
In “Full factorial” section, single and combinatorial effects of the media components were investigated on the basis of a full factorial design providing a rough estimate about optimal medium. The goal of the next step was to examine limitations of the particular effects and to identify potential optimal media compositions. To achieve this, the minimum and maximum concentrations were adjusted, and a more complex experimental design scheme was applied, comprising several nested factorial designs (see Additional file 6 for the full experimental design including the corresponding measurement data).
The maximum concentration of NaNO3 was lowered from 1.7 × Ref to 1 × Ref. The upper bound of the concentration of the clustered trace elements was increased by 50% to 3.75 × Ref. The concentration of CaCl2 could not be increased, as various types of precipitation effects were observed that distorted lipid analysis (see Additional file 7).
However, MgSO4 was varied over three levels, as illustrated in Fig. 3a. For each level, the concentrations of NaNO3, CaCl2, and trace elements were distributed using a full factorial design. For the intermediate concentration of MgSO4, the remaining components were varied only over half of their total ranges. A center point was located in each of these full factorial cubes. An additional nine points were space filling distributed over the edges of the cubes. In total, 39 experiments were performed and analyzed, including four reference replicates.
Figure 3b shows the Kriging interpolation based on all the data available after the third round of experiments. The figure shows three contour plots where the third component was fixed to the front, bottom–left corner of the inner cube as shown in Fig. 3a. The contour plots clearly show an interaction of the trace elements with NaNO3 and CaCl2, whereas MgSO4 influences the lipid productivity only slightly positively. Moreover, an optimal region for the medium composition can be identified around MgSO4 = 3.25 × Ref, CaCl2 = 1.5 × Ref, Trace = 2 × Ref, and NaNO3 = 0.3 × Ref (A three-dimensional plot of the Kriging model together with the measured data can be obtained from Additional file 8).
Refining the optimum
In the fourth and last round of the experiments, twelve experiments were placed around the optimum predicted by the Kriging interpolation. These experiments were planned by sampling the EI distribution, as described in “Expected Improvement” section, for maximizing the lipid productivity and minimizing the prediction uncertainty of the Kriging model. In addition, 23 experiments were uniformly distributed over the parameter space in a random manner, in order to improve prediction accuracy also in non-optimal regions. In total, 39 experiments were performed, including the four reference experiments (The full experimental design including the respective measurement data can be obtained from Additional file 9).
Figure 4 shows predictions of the updated Kriging model in the same fashion as described in “Locating optimal medium composition” section. The location of the optimum shifted toward MgSO4 = 3.25 × Ref, CaCl2 = 1.25 × Ref, trace elements = 2.5 × Ref, and NaNO3 = 0.45 × Ref (A three-dimensional plot of the Kriging model together with the measured data can be obtained from Additional file 10). For the optimal medium composition, the Kriging model predicts an increase by a factor of 3.03 ± 0.81 in lipid productivity compared with the reference medium.
Validating the optimal medium composition
In order to validate the determined optimal medium composition (see “Kriging-assisted optimization” section) and to highlight potential changes of process kinetics, cultivations using enBBMref and enBBMopt were carried out (see Additional file 2 for medium composition). Both processes were monitored in-depth by sequential harvest of replicate wells from microtiter plate cultivations (see Fig. 5). To maximize comparability with the literature reports, biomass concentration at harvest was acquired as cell dry weight rather than biovolume in this context.
Medium optimization resulted in a series of significant changes in process performance as summarized in Table 2. While the exponential growth rates in both media did not differ significantly (p < 0.05, t test), times to nitrogen depletion were 84 h and 52 h for enBBMref and enBBMopt, respectively. This was due to the reduction of nitrate concentration during medium optimization down to 0.45 × Ref. In the reference process, exponential growth shifted to linear kinetics reaching an optical density of 4.94 ± 0.06 typically indicating the onset of light limitation and in clear accordance with prior experiments . This effect was not observed for the optimized medium before nitrogen depletion. Neutral lipid accumulations started within 36 h (enBBMref) and 20 h (enBBMopt) after nitrogen limitation which corresponds to a reduction of approx. 45%. Moreover, the biomass-specific lipid accumulation rate (estimated by linear fit) increased by approx. 32% from 4.87 ± 0.53% (w w−1) d−1 to 6.43 ± 0.17% (w w−1) d−1 due to medium optimization. Most probably, both effects are attributable to the increased availability of magnesium and calcium ions, as well as trace elements in the medium. This might result in a boost of the enzymatic turnover of lipid synthesis, especially regarding acetyl-CoA carboxylase (see Table 1).
Alternatively, a kinetic limitation of ion import into the cells at the low concentrations in the reference medium could be an explanation. Regarding downstream processing, the increased magnesium concentration offers another positive aspect, as it was previously reported to assist flocculation of the cells at high pH . This mechanism is currently being investigated as an alternative to the comparably costlier biomass separation by centrifugation.
Most strikingly, cell dry weight at harvest did not differ significantly (p < 0.05, t test) for both media, despite the nitrate concentration being reduced to 45% in enBBMopt. This indicates the nitrate-specific biomass yield as being a function of the initial nitrate availability, a phenomenon that has recently been recognized and discussed for a fairly comparable Chlorella process . Together with an increase in the neutral lipid content from 10.55 ± 0.35% (w w−1) to 23.9 ± 1.2% (w w−1), this translated into a 2.3-fold increase of volumetric productivity up to 169 ± 7 mg L−1 d−1.
Besides the evaluation of productivity-related issues, the relative composition of the fatty acids from the neutral lipid product fraction was compared by gas chromatography time-of-flight mass spectrometry (see Fig. 6).
The obtained fingerprints were in clear agreement with the previous literature reports  as palmitic, oleic, linoleic, and α-linolenic acids made up the major product fractions of 85% (enBBMref) and 89% (enBBMopt). There are indications that the lipid fingerprint largely depends on cultivation conditions such as temperature , illumination , etc. However, our results demonstrate that changes in the medium composition can also lead to differences in the fatty acid fingerprint. The fractions of palmitoleic (16:1 Δ9), hexadecadienoic (16:2 Δ7,10), hexadecatrienoic (16:3 Δ7,10,13), stearic (18:0), and linoleic (18:2 Δ9,12) remained nearly unchanged. On the contrary, the proportions of palmitic (16:0) and α-linolenic (18:3 Δ9,12,15) acids shrank by 22 and 42%, respectively, while linoleic (18:1 Δ9) acid increased by 92% to a total share of 48 ± 1.8% using enBBMopt. With respect to biodiesel synthesis, this reduction in the polyunsaturated fatty acids’ fraction is clearly advantageous, increasing the fuel’s oxidative stability .
Final medium simplification
In “Kriging-assisted optimization” section, several input variables were identified to be ‘non-relevant’ and thus kept at the respective reference values throughout the whole study. However, for MES and especially for EDTA, a negative, but still non-significant (p < 0.05, t test) trend was observed. Besides economic aspects, culture media should only contain the necessary ingredients in appropriate concentrations to ensure high nutrient usage efficiency. Thus, an additional variant, in the following denoted as enBBMopt,min, was investigated. Here, the concentrations of all ‘non-relevant’ components were set to the respective minimum values during the screening analysis. In particular, this included the complete omissions of MES buffer and the chelator EDTA, as well as NaCl (see Additional file 2 for medium composition), while phosphate availability was reduced to 0.125 × Ref.
In comparison with the results using enBBMref, as well as enBBMopt, these adaptations did not change the overall obtained cell dry weight significantly (p < 0.05, t test) but led to an increase of the neutral lipid to 30.1 ± 1.6% (w w−1), while the respective lipid fingerprint remained unchanged in comparison with enBBMopt (see Additional file 11). The resulting volumetric productivity of 212 ± 11 mg L−1 d−1 represents a total 2.9-fold improvement compared with the reference. Leaving out EDTA and especially the MES buffer drastically reduces the medium costs, so that the price per liter is lowered by 96%. Most probably, MES is not required as the phosphate salts offer sufficient pH stabilization capacity. Although EDTA is commonly used as a metal chelator to improve long-term stability of algae cultivation media, the results clearly indicate that its usage is not beneficial for this specific application. Moreover, the reduction of phosphate concentration to 12.5% is advantageous for large-scale application where the recovery of excess nutrients to prevent overfertilization by wastewater is an important economic aspect. Yet, these results clearly confirm the validity of the initial screening analysis.
Assessment of achieved volumetric productivity
In the last decade, numerous studies addressed the lipid production of diverse C. vulgaris strains in different laboratory-scale batch processes [6, 17, 58, 70, 71, 73–81]. Among these, the average volumetric productivity was approximately 51 ± 36 mg L−1 d−1 and thus was comparable to the achieved value of 74 ± 1 mg L−1 d−1 using the enBBM reference medium. However, the reported values exhibit a wide spread, and it has to be assumed that these differences do not only originate from the different strains used, but from process conditions and reactor design as well. Some studies report productivities in the range of 130 mg L−1 d−1 when cultivating C. vulgaris in laboratory-scale batch processes with optimized nitrogen availability [17, 80]. Unfortunately, it is not generally clarified if productivities refer to the neutral lipid or the total lipid content. In this study, the volumetric productivity of neutral lipids of up to 212 ± 11 mg L−1 d−1 clearly exceeds previous reports and thereby underlines the importance of medium optimization not only for nitrate as commonly done, but especially for the concentrations of further salts and trace ions.
In this study, a blueprint strategy for the accelerated development of phototrophic bioprocesses is presented. This strategy is very efficient in terms of time and material, by incorporating state-of-the-art phototrophic cultivation and analytics with higher throughput that is closely linked to sophisticated experimental design strategies.
Taking neutral lipid production by the unicellular microalga C. vulgaris as a model process, the cultivation medium was optimized toward volumetric productivity. Fractional and full factorial designs in combination with Kriging-based approaches for data analysis, visualization, and experimental design allowed for an efficient and effective optimization in terms of time and cost. The optimized process has an approximately threefold increased lipid productivity of 212 ± 11 mg L−1 d−1, which was achieved with only four experimental rounds with one microtiter plate each.
Besides the commonly addressed concentration of the nitrogen source (here nitrate), especially magnesium, calcium, and various trace elements were shown to be of crucial importance. Analysis tools furthermore revealed multi-parameter interactions that could have been overlooked otherwise. Over and above this, the concentration of non-relevant medium components was successfully minimized, contributing to reducing medium cost. Taking all the above results together, a smart combination of microscale phototrophic cultivation with sophisticated design of experiments led to a tremendous improvement of neutral lipid production with C. vulgaris, at the same time reducing cost for media components by 96%, while all other process performance indicators were kept constant.
Design of Experiments
enriched Bold’s Basal Medium
Markov Chain Monte Carlo
- nLC [% (w w−1)]:
neutral lipid content
- OD750 [−]:
optical density at 750 nm
- Pvol [mg L−1 d−1]:
- µ [d−1]:
exponential growth rate
Mata TM, Martins AA, Caetano NS. Microalgae for biodiesel production and other applications: a review. Renew Sust Energ Rev. 2010;14:217–32.
Weyer KM, Bush DR, Darzins A, Willson BD. Theoretical maximum algal oil production. Bioenerg Res. 2010;3:204–13.
Gouveia L, Oliveira A. Microalgae as a raw material for biofuels production. J Ind Microbiol Biotechnol. 2009;36:269–74.
Gui MM, Lee KT, Bhatia S. Feasibility of edible oil vs. non-edible oil vs. waste edible oil as biodiesel feedstock. Energy. 2008;33:1646–53.
Yang J, Xu M, Zhang X, Hu Q, Sommerfeld M, Chen Y. Life-cycle analysis on biodiesel production from microalgae: water footprint and nutrients balance. Bioresour Technol. 2011;102:159–65.
Rodolfi L, Chini Zittelli G, Bassi N, Padovani G, Biondi N, Bonini G, et al. Microalgae for oil: strain selection, induction of lipid synthesis and outdoor mass cultivation in a low-cost photobioreactor. Biotechnol Bioeng. 2009;102:100–12.
Chisti Y. Biodiesel from microalgae. Biotechnol Adv. 2007;25:294–306.
Wijffels RH, Barbosa MJ. An outlook on microalgal biofuels. Science. 2010;329:796–9.
Wijffels RH, Barbosa MJ, Eppink MHM. Microalgae for the production of bulk chemicals and biofuels. Biofuel Bioprod Biorefin. 2010;4:287–95.
Norsker N-H, Barbosa MJ, Vermuë MH, Wijffels RH. Microalgal production—a close look at the economics. Biotechnol Adv. 2011;29:24–7.
Hariskos I, Posten C. Biorefinery of microalgae—opportunities and constraints for different production scenarios. Biotechnol J. 2014;9:739–52.
Brennan L, Owende P. Biofuels from microalgae—a review of technologies for production, processing, and extractions of biofuels and co-products. Renew Sust Energ Rev. 2010;14:557–77.
Adarme-Vega TC, Lim DKY, Timmins M, Vernen F, Li Y, Schenk PM. Microalgal biofactories: a promising approach towards sustainable omega-3 fatty acid production. Microb Cell Fact. 2012;11:96.
Van Wagenen J, Holdt SL, De Francisci D, Valverde-Pérez B, Plósz BG, Angelidaki I. Microplate-based method for high-throughput screening of microalgae growth potential. Bioresour Technol. 2014;169:566–72.
Radmann EM, Camerini FV, Santos TD, Costa JAV. Isolation and application of SOX and NOX resistant microalgae in biofixation of CO2 from thermoelectricity plants. Energy Convers Manag. 2011;52:3132–6.
de Morais MG, Costa JAV. Isolation and selection of microalgae from coal fired thermoelectric power plant for biofixation of carbon dioxide. Energy Convers Manag. 2007;48:2169–73.
Breuer G, Lamers PP, Martens DE, Draaisma RB, Wijffels RH. The impact of nitrogen starvation on the dynamics of triacylglycerol accumulation in nine microalgae strains. Bioresour Technol. 2012;124:217–26.
Debska D, Potvin G, Lan C, Zhang Z. Effects of medium composition on the growth of Chlorella vulgaris during photobioreactor batch cultivations. J Biobased Mater Bio. 2010;4:68–72.
Radzun KA, Wolf J, Jakob G, Zhang E, Stephens E, Ross I, et al. Automated nutrient screening system enables high-throughput optimisation of microalgae production conditions. Biotechnol Biofuels. 2015;8:1–17.
Tillich UM, Wiolter N, Schulze K, Kramer D, Brödel O, Frohme M. High-throughput cultivation and screening platform for unicellular phototrophs. BMC Microbiol. 2014;14:239.
Chen M, Mertiri T, Holland T, Basu AS. Optical microplates for high-throughput screening of photosynthesis in lipid-producing algae. Lab Chip. 2012;12:3870–4.
Han W, Li C, Miao X, Yu G. A novel miniature culture system to screen CO2-sequestering microalgae. Energies. 2012;5:4372–89.
Heo J, Cho D-H, Ramanan R, Oh H-M, Kim H-S. PhotoBiobox: a tablet sized, low-cost, high throughput photobioreactor for microalgal screening and culture optimization for growth, lipid content and CO2 sequestration. Biochem Eng J. 2015;103:193–7.
Ojo EO, Auta H, Baganz F, Lye GJ. Design and parallelisation of a miniature photobioreactor platform for microalgal culture evaluation and optimisation. Biochem Eng J. 2015;103:93–102.
Kim HS, Weiss TL, Thapa HR, Devarenne TP, Han A. A microfluidic photobioreactor array demonstrating high-throughput screening for microalgal oil production. Lab Chip. 2014;14:1415–25.
Graham PJ, Riordon J, Sinton D. Microalgae on display: a microfluidic pixel-based irradiance assay for photosynthetic growth. Lab Chip. 2015;15:3116–24.
Morschett H, Schiprowski D, Müller C, Mertens K, Felden P, Meyer J, et al. Design and validation of a parallelized micro-photobioreactor enabling phototrophic bioprocess development at elevated throughput. Biotechnol Bioeng. 2017;114:122–31.
Morschett H, Reich S, Wiechert W, Oldiges M. Simplified cryopreservation of the microalga Chlorella vulgaris integrating a novel concept for cell viability estimation. Eng Life Sci. 2016;16:36–44.
Morschett H, Wiechert W, Oldiges M. Automation of a Nile red staining assay enables high throughput quantification of microalgal lipid production. Microb Cell Fact. 2016;15:34.
Pereira H, Barreira L, Mozes A, Florindo C, Polo C, Duarte C, et al. Microplate-based high throughput screening procedure for the isolation of lipid-rich marine microalgae. Biotechnol Biofuels. 2011;4:1–12.
Tillich UM, Wolter N, Franke P, Dühring U, Frohme M. Screening and genetic characterization of thermo-tolerant Synechocystis sp. PCC6803 strains created by adaptive evolution. BMC Biotechnol. 2014;14:1–15.
Mandenius C-F, Brundin A. Bioprocess optimization using design-of-experiments methodology. Biotechnol Progr. 2008;24:1191–203.
Fisher RA. The design of experiments. Edinborough: Oliver and Boyd; 1935.
Freier L, Hemmerich J, Schöler K, Wiechert W, Oldiges M, von Lieres E. Framework for Kriging-based iterative experimental analysis and design: optimization of secretory protein production in Corynebacterium glutamicum. Eng Life Sci. 2016. doi:10.1002/elsc.201500171.
Montgomery DC. Design and analysis of experiments. New York: Wiley; 2012.
Lee K-M, Gilmore DF. Statistical experimental design for bioprocess modeling and optimization analysis. Appl Biochem Biotechnol. 2006;135:101–15.
Islam RS, Tisi D, Levy MS, Lye GJ. Framework for the rapid optimization of soluble protein expression in Escherichia coli combining microscale experiments and statistical experimental design. Biotechnol Progr. 2007;23:785–93.
Beijerinck MW. Culturversuche mit Zoochlorellen und anderen niederen Algen. Btg Ztg. 1890;45:725–40.
Bold HC. The morphology of Chlamydomonas chlamydogama, sp. nov. Bull Torrey Bot Club. 1949;76:101–8.
Rohe P, Venkanna D, Kleine B, Freudl R, Oldiges M. An automated workflow for enhancing microbial bioprocess optimization on a novel microbioreactor platform. Microb Cell Fact. 2012;11:144.
Unthan S, Radek A, Wiechert W, Oldiges M, Noack S. Bioprocess automation on a Mini Pilot Plant enables fast quantitative microbial phenotyping. Microb Cell Fact. 2015;14:32.
Graham MD. The Coulter principle: foundation of an industry. JALA. 2003;8:72–81.
Paczia N, Nilgen A, Lehmann T, Gätgens J, Wiechert W, Noack S. Extensive exometabolome analysis reveals extended overflow metabolism in various microorganisms. Microb Cell Fact. 2012;11:122.
Myers RH, Montgomery DC, Anderson-Cook CM. Response surface methodology: process and product optimization using designed experiments. Hoboken: Wiley; 2016.
Cressie N. Statistics for spatial data. Hoboken: Wiley-Intersciences; 2015.
Jones DR, Schonlau M, Welch WJ. Efficient global optimization of expensive black-box functions. J Glob Optim. 1998;13:455–92.
Jones DR. A taxonomy of global optimization methods based on response surfaces. J Glob Optim. 2001;21:345–83.
Freier L, von Lieres E. Multi-objective global optimization (MOGO) Algorithm and case study in gradient elution chromatography. Biotechnol J. doi:10.1002/biot.201600613.
Haario H, Laine M, Mira A, Saksman EDRAM. Efficient adaptive MCMC. Stat Comput. 2006;16:339–54.
Chen H, Zhang Y, He C, Wang Q. Ca2+ signal transduction related to neutral lipid synthesis in an oil-producing green alga Chlorella sp. C2. Plant Cell Physiol. 2014;55:634–44.
Ren H-Y, Liu B-F, Kong F, Zhao L, Xie G-J, Ren N-Q. Enhanced lipid accumulation of green microalga Scenedesmus sp. by metal ions and EDTA addition. Bioresour Technol. 2014;169:763–7.
Jones RGW, Lunt OR. The function of calcium in plants. Bot Rev. 1967;33:407–26.
Liu ZY, Wang GC, Zhou BC. Effect of iron on growth and lipid accumulation in Chlorella vulgaris. Bioresour Technol. 2008;99:4717–22.
Ruangsomboon S, Ganmanee M, Choochote S. Effects of different nitrogen, phosphorus, and iron concentrations and salinity on lipid production in newly isolated strain of the tropical green microalga Scenedesmus dimorphus KMITL. J Appl Phycol. 2013;25:867–74.
Concas A, Steriti A, Pisu M, Cao G. Comprehensive modeling and investigation of the effect of iron on the growth rate and lipid accumulation of Chlorella vulgaris cultured in batch photobioreactors. Bioresour Technol. 2014;153:340–50.
Guckert JB, Cooksey KE. Triglyceride accumulation and fatty acid profile changes in Chlorella (Chlorophyta) during high pH-induced cell cycle inhibition. J Phycol. 1990;26:72–9.
Deng X, Fei X, Li Y. The effects of nutritional restriction on neutral lipid accumulation in Chlamydomonas and Chlorella. Afr J Microbiol Res. 2011;5:260–70.
Lv JM, Cheng LH, Xu XH, Zhang L, Chen HL. Enhanced lipid production of Chlorella vulgaris by adjustment of cultivation conditions. Bioresour Technol. 2010;101:6797–804.
Alyabyev AJ, Loseva NL, Gordon LK, Andreyeva IN, Rachimova GG, Tribunskih VI, et al. The effect of changes in salinity on the energy yielding processes of Chlorella vulgaris and Dunaliella maritima cells. Thermochim Acta. 2007;458:65–70.
Duan X, Ren GY, Liu LL, Zhu WX. Salt-induced osmotic stress for lipid overproduction in batch culture of Chlorella vulgaris. Afr J Biotechnol. 2012;11:7072–8.
Mallick N. Biotechnological potential of immobilized algae for wastewater N, P and metal removal: a review. Biometals. 2002;15:377–90.
Pittman JK, Dean AP, Osundeko O. The potential of sustainable algal biofuel production using wastewater resources. Bioresour Technol. 2011;102:17–25.
Muñoz R, Guieysse B. Algal–bacterial processes for the treatment of hazardous contaminants: a review. Wat Res. 2006;40:2799–815.
Chen M, Tang H, Ma H, Holland TC, Ng KYS, Salley SO. Effect of nutrients on growth and lipid accumulation in the green algae Dunaliella tertiolecta. Bioresour Technol. 2011;102:1649–55.
Yang J, Cao J, Xing G, Yuan H. Lipid production combined with biosorption and bioaccumulation of cadmium, copper, manganese and zinc by oleaginous microalgae Chlorella minutissima UTEX2341. Bioresour Technol. 2015;175:537–44.
Clijsters H, Van Assche F. Inhibition of photosynthesis by heavy metals. Photosynth Res. 1985;7:31–40.
Ohlrogge J, Browse J. Lipid biosynthesis. Plant Cell. 1995;7:957–70.
Vandamme D, Foubert I, Fraeye I, Meesschaert B, Muylaert K. Flocculation of Chlorella vulgaris induced by high pH: role of magnesium and calcium and practical implications. Bioresour Technol. 2012;105:114–9.
Petkov G, Garcia G. Which are fatty acids of the green alga Chlorella? Biochem Syst Ecol. 2007;35:281–5.
Converti A, Casazza AA, Ortiz EY, Perego P, Del Borghi M. Effect of temperature and nitrogen concentration on the growth and lipid content of Nannochloropsis oculata and Chlorella vulgaris for biodiesel production. Chem Eng Process. 2009;48:1146–51.
Atta M, Idris A, Bukhari A, Wahidin S. Intensity of blue LED light: a potential stimulus for biomass and lipid content in fresh water microalgae Chlorella vulgaris. Bioresour Technol. 2013;148:373–8.
Wadumesthrige K, Salley SO, Ng KYS. Effects of partial hydrogenation, epoxidation, and hydroxylation on the fuel properties of fatty acid methyl esters. Fuel Process Technol. 2009;90:1292–9.
Widjaja A, Chien C-C, Ju Y-H. Study of increasing lipid production from fresh water microalgae Chlorella vulgaris. J Taiwan Inst Chem E. 2009;40:13–20.
Yeh K-L, Chang J-S. Effects of cultivation conditions and media composition on cell growth and lipid productivity of indigenous microalga Chlorella vulgaris ESP-31. Bioresour Technol. 2012;105:120–7.
Mallick N, Mandal S, Singh AK, Bishai M, Dash A. Green microalga Chlorella vulgaris as a potential feedstock for biodiesel. J Chem Technol Biotechnol. 2012;87:137–45.
Liang Y, Sarkany N, Cui Y. Biomass and lipid productivities of Chlorella vulgaris under autotrophic, heterotrophic and mixotrophic growth conditions. Biotechnol Lett. 2009;31:1043–9.
Hsieh C-H, Wu W-T. Cultivation of microalgae for oil production with a cultivation strategy of urea limitation. Bioresour Technol. 2009;100:3921–6.
Griffiths M, Hille R, Harrison SL. The effect of nitrogen limitation on lipid productivity and cell composition in Chlorella vulgaris. Appl Microbiol Biotechnol. 2014;98:2345–56.
Gorain PC, Bagchi SK, Mallick N. Effects of calcium, magnesium and sodium chloride in enhancing lipid accumulation in two green microalgae. Environ Technol. 2013;34:1887–94.
Yeh K-L, Chang J-S. Nitrogen starvation strategies and photobioreactor design for enhancing lipid content and lipid production of a newly isolated microalga Chlorella vulgaris ESP-31: implications for biofuels. Biotechnol J. 2011;6:1358–66.
Yoo C, Jun SY, Lee JY, Ahn CY, Oh HM. Selection of microalgae for lipid production under high levels carbon dioxide. Bioresour Technol. 2010;101:71–4.
HM and LF designed the study, evaluated the results, and prepared the manuscript. HM carried out the experimental work with assistance from JR, while LF performed the computational analysis of the acquired data. WW helped to finalize the manuscript and EvL and MO helped to finalize the manuscript and supervised the computational and experimental work, respectively. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of supporting data
The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.
The authors thank the Federal Ministry for Economic Affairs and Energy on the basis of a grant by the German Bundestag for the support and funding (grant no. KF2519304CS3) for Holger Morschett, as well as the Ministry of Innovation, Science and Research of North Rhine-Westphalia and the Heinrich Heine University Düsseldorf for the scholarship awarded to Lars Freier within the purview of CLIB-Graduate Cluster Industrial Biotechnology.
Holger Morschett and Lars Freier contributed equally to this work
Additional file 1. Fractional factorial design for the initial screening analysis and corresponding measurement data of the individual cultivations.
Additional file 3. Estimated effect of two factor interactions with MgSO4. Estimations are based on the experiments in section “Fractional factorial” using the fractional factorial design and the corresponding measurement data given in Additional material 1.
Additional file 4. Experimental design and corresponding measurement data evaluated in section “Full factorial”.
Additional file 5. Screening plot around reference point (enBBMref medium). Estimation of the functional relationship between media components and lipid productivity was done by Kriging. The Kriging model is based on the experiments in section “Kriging-assisted optimization” using the open source software KriKit.
Additional file 6. Experimental design and corresponding measurement data evaluated in section “Locating optimal medium composition”.
Additional file 7. Analysis of calcium precipitation by means of optical density measurements. Shaded area represents the parameter space covered during cultivation experiments. Error bars deviated from technical replicates (n = 3).
Additional file 8. Measured data and Kriging interpolation after the third round of experiments in three-dimensional representation. The non-varied components were fixed at MgSO4=1.7xRef, CaCl2=0.82xRef, Trace=1xRef, or NaNO3=0.37xRef. Red dots indicate measurement data.
Additional file 9. Experimental design and corresponding measurement data evaluated in section “Refining the optimum”.
Additional file 10. Measurement data and updated Kriging interpolation after the fourth round of experiments in three-dimensional representation. The non-varied components were fixed at MgSO4=1.7xRef, CaCl2=0.82xRef, Trace=1xRef, or NaNO3=0.37xRef. Red dots indicate measurement data.
Additional file 11. Relative composition of the fatty acids from the neutral lipid product fraction for enBBMopt and enBBMopt,min. Error bars represent min/max from biological replicates (n = 2).
About this article
Cite this article
Morschett, H., Freier, L., Rohde, J. et al. A framework for accelerated phototrophic bioprocess development: integration of parallelized microscale cultivation, laboratory automation and Kriging-assisted experimental design. Biotechnol Biofuels 10, 26 (2017). https://doi.org/10.1186/s13068-017-0711-6