# Scale-Up to Pilot of a Non-Axenic Culture of Thraustochytrids Using Digestate from Methanization as Nitrogen Source

Denis de La Broise, Mariana Ventura, Luc Chauchat, Maurean Guerreiro, Teo Michez, Thibaud Vinet, Nicolas Gautron, Fabienne Le Grand, Nelly Le Goïc, Antoine Bideau, et al.

► To cite this version:

Denis de La Broise, Mariana Ventura, Luc Chauchat, Maurean Guerreiro, Teo Michez, et al. Scale-Up to Pilot of a Non-Axenic Culture of Thraustochytrids Using Digestate from Methanization as Nitrogen Source. Marine drugs, 2022, 20 (8), pp.499. 10.3390/md20080499. hal-04237349

HAL Id: hal-04237349

https://hal.univ-brest.fr/hal-04237349

Submitted on 12 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Scale-Up to Pilot of a Non-Axenic Culture of Thraustochytrids Using Digestate from Methanization as Nitrogen Source

Denis de la Broise <fnref n="1" />,*, Mariana Ventura <fnref n="1" />, Luc Chauchat <fnref n="1" />, Maurean Guerreiro <fnref n="1" />,†, Teo Michez <fnref n="1" />, Thibaud Vinet <fnref n="2" />, Nicolas Gautron <fnref n="1" />,†, Fabienne Le Grand <fnref n="1" />,†, Antoine Bideau <fnref n="1" />, Nelly Le Goïc <fnref n="1" />, Adeline Bidault <fnref n="1" />, Christophe Lambert <fnref n="1" />,† and Philippe Soudant <fnref n="1" />,*,†

**Abstract:** The production of non-fish based docosahexaenoic acid (DHA) for feed and food has become a critical need in our global context of over-fishing. The industrial-scale production of DHA-rich Thraustochytrids could be an alternative, if costs turned out to be competitive. In order to reduce production costs, this study addresses the feasibility of the non-axenic (non-sterile) cultivation of *Aurantiochytrium mangrovei* on industrial substrates (as nitrogen and mineral sources and glucose syrup as carbon and energy sources), and its scale-up from laboratory (250 mL) to 500 L cultures. Pilot-scale reactors were airlift cylinders. Batch and fed-batch cultures were tested. Cultures over 38 to 62 h achieved a dry cell weight productivity of 3.3 to 5.5 g.L⁻<fnref n="1" />.day⁻<fnref n="1" />, and a substrate to biomass yield of up to 0.3. DHA productivity ranged from 10 to 0.18 mg.L⁻<fnref n="1" />.day⁻<fnref n="1" />. Biomass productivity appears linearly related to oxygen transfer rate. Bacterial contamination of cultures was low enough to avoid impacts on fatty acid composition of the biomass. A specific work on microbial risks assessment (in supplementary files) showed that the biomass can be securely used as feed. However, to date, there is a law void in EU legislation regarding the recycling of nitrogen from digestate from animal waste for microalgae biomass and its usage in animal feed. Overall, the proposed process appears similar to the industrial yeast production process (non-axenic heterotrophic process, dissolved oxygen supply limiting growth, similar cell size). Such similarity could help in further industrial developments.

**Keywords:** Thraustochytrids; *Aurantiochytrium*; PUFA; DHA; digestate; oxygen; scale-up; pilot

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** de la Broise, D.; Ventura, M.; Chauchat, L.; Guerreiro, M.; Michez, T.; Vinet, T.; Gautron, N.; Le Grand, F.; Bideau, A.; Goïc, N.L.; et al. Scale-Up to Pilot of a Non-Axenic Culture of Thraustochytrids Using Digestate from Methanization as Nitrogen Source. *Mar. Drugs* **2022**, *20*, 499. https://doi.org/10.3390/ md20080499

Academic Editor: Ipek Kurtboke

Received: 24 June 2022

Accepted: 28 July 2022

Published: 2 August 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/ 4.0/).

## 1. Introduction

Oceans are today exposed to dramatic pressure on fishing resources. The unsustain-ability of marine fish catches rose from 10–15% in the 1970s to 35% in 2017 [1]. From the approximately 85 million tons of marine fishes caught in 2018, 18 million were used for fish meal and fish oil production, mainly in aquaculture and breeding.

Fish oil is part of a large panel of the aquaculture diet, as it supplies long-chain n-3 polyunsaturated fatty acids (LC n-3 PUFA), a prerequisite for fish growth and for the nutritional quality of the final products [2].

Docosahexaenoic acid (DHA) is the most abundant LC n-3 PUFA in the mammalian central and peripheral nervous system and plays a crucial role in brain and retina development, cognitive function and visual acuity [3]. DHA is involved in pre- and postnatal development and is particularly important in the later stages of pregnancy and early infancy. Recent studies have also demonstrated that food intake of DHA has a beneficial effect on a number of psychological, behavioral, and neural disorders, such as Alzheimer's

<footnote n="1">¹ CNRS, IRD, Ifremer, Université de Bretagne Occidentale, UMR6539 LEMAR, F-29280 Plouzané, France; mariana.ventura@univ-brest.fr (M.V.); luc.chauchat@univ-brest.fr (L.C.); guerreiro.maurean@outlook.com (M.G.); michez.teo@gmail.com (T.M.); nicolas.gautron.37@gmail.com (N.G.); fabienne.legrand@univ-brest.fr (F.L.G.); bideau.antoine56@gmail.com (A.B.); nelly.legoic@univ-brest.fr (N.L.G.); adeline.bidault@univ-brest.fr (A.B.); christophe.lambert@univ-brest.fr (C.L.)</footnote>

<footnote n="2">² DENITRAL, F-22400 Lamballe, France; thibaud.vinet@orange.fr</footnote>

<footnote>* Correspondence: denis.de-labroise@univ-brest.fr (D.d.l.B.); philippe.soudant@univ-brest.fr (P.S.)</footnote>
disease [4]. DHA might have played a critical role in human evolution, as DHA is a determinant of neuronal migration, neurogenesis and the expression of several genes involved in brain growth and function [5]. Although human beings can synthesize DHA from n-3 PUFA precursors, it is more efficient to have a dietary supply of DHA to meet our health requirements [2]. This high demand results in a clear risk that DHA supply, mostly from small pelagic fishes, may become insufficient for a growing human population by 2040 [6].

As the supply of PUFA from fish catches is becoming a hurdle [7], alternative sources are being investigated, and among them, research has focused on Thraustochytrids [8], planktonic microorganisms that exhibit high DHA content [9,10] and has led to recent industrial application [11]. Obviously, as annual fish oil production in the world accounts for as high as 1 million tons per year [1], the Thraustochytrids production could be a significant and realistic alternative only if developed up to a large industrial scale and at a competitive market price. Then, production cost at industrial scale for DHA-containing biomass would be equal or lower than the current market price for an equivalent DHA supply in products obtained from small pelagic fisheries.

For this to happen, the Thraustochytrids industrial production process should fulfill the requirements for both high productivity and low production cost, as observed for example, in the yeast industry.

One of the usual means to reduce microorganism production cost is, when possible, to run the microbial cultivation under non-axenic (or "non-sterile") conditions, as it leads to significant reductions in both investment and running costs. Such conditions were pointed out as key to competitive process development in "next generation biotechnology" [12]. When this approach is usual for large-scale phototrophic cultivations, its application to heterotrophic cultures can be more difficult, as massive bacterial contamination can rapidly appear in organic nutrient-rich culture media. However, this can be solved by establishing a competitive advantage for the cultivated strain. A competitive advantage can be obtained by means of high inoculation rate and short-term culture as applied for yeast production [13], by harsh environmental conditions (i.e., pH, temperature, salinity) that the cultivated strain can competitively withstand [14], or even by strain engineering for its growth on a specific substrate that competitors cannot use [15]. High inoculation rates and short-term culture could be applied to Thraustochytrids.

Another approach to reduce production cost can be the use of fed-batch or continuous cultures [16,17], as they allow one to increase the productivity and thus reduce investment and process costs.

Substrates for culture medium often strongly affect the biomass production cost. As Thraustochytrids exhibit a heterotrophic metabolism, the organic substrates, sources of carbon and energy, are a key point in the final production cost. Therefore, they must be cheap and allow a high substrate-to-product conversion yield. Low-cost by-products from agriculture and from the food industry are available in large amounts and were investigated for Thraustochytrids production [18,19]. These include by-products from the food industry (whey and derivatives, plant process by-products) and from the anaerobic digestion (AD) process (digestates). They could be used in non-sterile microbial processes [12,20]. Using such by-products also contributes to the development of a circular economy as part of sustainability. The digestate is an emerging by-product, produced in large amounts in northwestern Europe. Its disposal as agricultural fertilizer leads, when in excess, to nitrogen and phosphorus pollution, resulting in eutrophication of natural waters.

As partner in the INTERREG North West Europe ALG-AD project, we investigated the scale-up of Thraustochytrids production using AD digestate as a source of nitrogen and minerals. As recycling animal waste is a sensitive issue in Europe, a regulatory review on the use of digestate to cultivate algal biomass for animal feed was realized within the framework of the ALG-AD project (see report in supplementary files). With the aim of industrial scalability, we focused on economically sound process development by complying with the above-described goals and constraints. The objective was thus to assess the feasibility of the biomass production process at pilot scale in non-axenic conditions, and
to achieve high productivity and cell density on AD digestate and corn syrup hydrolysate as main substrates, the ultimate objective being the industrial production of biomass for DHA and protein supplementation in feed. We describe in this paper the successive steps, up to a non-axenic process at 500 L scale. We present the identified key points of the process and the productivity and conversion yield observed.

## 2. Results

### 2.1. Cultures Submitted to Air Bubbling and Shaking Exhibit Linear Growth

The first scale-up step was to run 8 L cultures in a modified yeast extract and peptone (YEP) medium with continuous air supply (Cf. Section 4.4.1). The specific growth rate (Figure 1) exhibited an unusual linear shape, as compared to common microbial exponential growth rates. In usual exponential growth conditions, the biomass production per time unit increase as the cell density increases, usually expressed as:

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Biomass production kinetics of *A. mangrovei* cultivation in a modified YEP medium, with air bubbling 0.4 VVM. Volume of culture was 8 L before day 3, and 10 L from day 3 to day 10, due to addition of 2 L of a 500 g.L⁻¹ glucose solution on day 3. Three replicates, error bars represent S.D.</figcaption></figure>

$$X = X^{\circ} \times \text{Exp} \mu \times (t - t^{\circ}),$$

where $X^{\circ}$ is the biomass at time $t^{\circ}$, X is the biomass at time t and $\mu$ is the growth rate.

In this experiment, we observed a different, stable biomass production rate over time, which can be expressed as:

$$X = X^{\circ} + A \times (t - t^{\circ}).$$

where productivity is $A = 0.05 \text{ g.L}^{-1}\text{.h}^{-1}$.

In exponential growth, as described by Monod [21] the growth rate is determined by the limiting nutrient concentration. In this experiment, glucose was in large excess (16 g/L total added), and initial proteins from yeast extract (Y.E.) and peptones (approximately 3.5 g/L) represented 50% of the 7 g/L biomass dry cell weight (DCW) measured after 6 h of cultivation. This implies that the limiting factor for the growth rate was likely not related to nutrients. The constant biomass production per time unit suggests that the limitation is related to the constant and limiting supply of an external component. Preliminary experiments showed that growth remained poor when bubbling was lacking and dissolved oxygen values were below 10% after 2 days in bubbling cultures. This critical role of oxygen supply in Thraustochytrids cultures was already reported [11]. The volumetric oxygen mass transfer coefficient can be described as:

$$d\text{Cl}/dt = K\text{L.a} \times (C^{*} - \text{Cl})$$

where $d\text{Cl}/dt$ is the oxygen transfer rate, $K\text{L.a}$ is the volumetric oxygen mass transfer coefficient, $C^*$ is the dissolved oxygen value when medium is saturated, and $\text{Cl}$ is dissolved oxygen value at time $t$ [22]. Thus, in this experiment, the volumetric oxygen mass transfer coefficient $K\text{L.a}$ that drove oxygen transfer at constant rate, as dissolved oxygen was low and stable, was likely the limiting factor that drove the linear biomass production rate.

Then, the productivity can be expressed as a linear function of the oxygen transfer $K\text{L.a}$, and the biomass production rate can be expressed as:
$$X = X^{\circ} + \alpha \times Kl \cdot a \times (t - t^{\circ})$$

where α is a strain related constant value when the substrate is not limiting.

This model points out the oxygen transfer rate as the factor that drives the productivity in such cultures, when nutrients are not limiting. Therefore, during the further scale-up steps, special attention was dedicated to oxygen transfer efficiency. Further experiments should be run to confirm this proposed model.

## 2.2. Comparison of Performance of Cultures in Digestates from Three Different Sources

The performance of digestates, made from three different anaerobic digestion biomass sources, was compared over 136 h in 800 mL of aerated and agitated cultures in flasks. Based on previous NH₄ content analysis, the digestate volumes added in cultures were adjusted in order to start the culture with 100 mg.L⁻¹ of NH₄ in each culture. The added volume of NH₄Cl solution was calculated to obtain 100 mg.L⁻¹ in control culture (Cf. Section 4.1.2). Biomass production (Figure 2A) and glucose consumption (Figure 2B) showed the Cooperl (pig manure) digestate gave the highest biomass (3.2 g.L⁻¹ DCW), while Langage digestate (agricultural and food by-products) and plant-based digestate (U. Ghent) allowed 2 g/L and 1.6 g/L, respectively. However, all digestates allowed a biomass production higher than the NH₄Cl control (1.5 g.L⁻¹ DCW) (as digestates were filtered at 0.2 µm, their contribution to biomass was negligible). This suggests that digestates did not only supply N-NH₄, but also other growth-promoting nutrients. Lipid content, as illustrated by fluorescence level of BODIPY-stained cells (Figure 2C), dropped at the end of culture. Together with the increase in N-NH₄ values (Figure 2D), these drops could be due to cell stress as other nutrient availability dropped (Figure 2B). The following experiments were all performed with Cooperl digestate.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Cultures in bubbling and shaking 800 mL supplemented with digestates from three different sources and with 0.297 g.L⁻¹ NH₄Cl (final concentration) as control source of nitrogen (see Section 4.1.2 for detailed description of culture medium). Mean values (n = 3, error bars = S.D.). (A): biomass production kinetics; (B): glucose concentration in cultures; (C): relative lipid content in cells as measured by BODIPY-stained cell fluorescence; (D): NH₄ concentrations in cultures.</figcaption></figure>
## 2.3. The Digestate Concentration Influences the Growth Kinetics in Flask Cultures

The influence of digestate supply quantity was investigated in flask cultures (Figure 3). The addition of 5% digestate allowed, as compared to YEP modified control condition, a slightly higher cell concentration between 7 and 10 days, and the same final cell concentration after 13 days.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Biomass growth kinetics of *A. mangrovei* in 150 mL cultures in 500 mL shaken flasks. Digestate as % v/v. Mean values (n = 3, error bars = S.D.). Cell concentrations were estimated using O.D. measurement at 620 nm. Relationship between cell concentration and O.D. at 620 nm was previously established by the following equation: cell concentration = $1 \times 10^8 \times (\text{O.D. at 620 nm}) - 1 \times 10^7 (R^2 = 0.98)$.</figcaption></figure>

Higher digestate concentrations induced a reduction in cell concentration (10% digestate), or worse, inhibited growth (20% digestate). We concluded that digestate exhibited a positive effect only for the lowest percentage tested (5%), while a toxic effect seemed to appear for higher percentages. This is possibly due to the toxicity of ammonium in solution [23], and/or possibly to other unidentified deterrent compounds. We then decided to run the following experiments with 2.5% digestate in batch cultures and up to 5% in fed-batch cultures.

## 2.4. Batch Cultures Can Be Run in Non-Axenic Airlift Cylinders

In order to study the scalability of the culture on digestate and in non-axenic conditions, cultures were run in 10 L PMMA cylinders. Two options were tested for air supply: large bubbles (4.5 mm diameter bottom tube) similar to carboy cultivation, and smaller bubbles (five holes, 1.5 mm diameter, in a bottom PVC plate). Three flow rates (0.8, 0.4 and 0.2 VVM) were tested with the small bubbles. As an axenic control, a cultivation was run in an agitated carboy with the large bubbles for aeration. These cultures were preliminary tests prior to pilot scale-up and were not replicated.

Results (Figure 4a) showed that the highest biomass concentration value ($3.7 \text{ g.L}^{-1}$) in non-axenic cylinder was obtained after 66 h of culture, with small bubbles at the maximum tested flow rate (0.8 VVM). Both cylinders with small bubbles at lower flow rates (0.4 and 0.2 VVM) allowed only poor cell and biomass production (Figure 4a,b). As the dissolved oxygen values were also the lowest in these two cylinders (Figure 4c), oxygen transfer limitation likely explained this lower growth performance. At 0.4 VVM, large bubbles, as compared to small bubbles, allowed higher dissolved oxygen values, with higher cell growth and biomass values. We hypothesized that the airlift effect was less efficient when small bubbles were used, then reducing oxygen transfer and mixing.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Effect of air bubbling systems on cultures (6 L) in non-axenic airlift cylinders and in an axenic carboy with air bubbling and 100 rpm shaking as control. (a): biomass production kinetics; (b): Cell concentration kinetics; (c): Dissolved oxygen concentration; (d): relative lipid content in cells as measured by BODIPY stained cell fluorescence.</figcaption></figure>

Intracellular lipid content was measured by flow cytometry (green fluorescence level in arbitrary unit after staining with BODIPY [20]). The BODIPY values were the highest (Figure 4d) at 66 h in the three cultures that exhibited the highest cell density, and were similar in axenic and non-axenic cultures. This suggests that cell size and/or cell lipid content were increased in these highly aerated cultures and that the potential bacterial presence did not dramatically modify the kinetics as compared to axenic cultures. Bacterial presence, as observed by microscopy (with a 40× objective and 10× eyepiece), was detected but low (few bacteria in the field) in all cylinders. Because we were using a large culture
volume (6 L), the experiment was performed without replicates. Therefore, it should only
be considered as indicative. Nevertheless, it showed that non-axenic cultivation appeared
feasible and permitted running further tests at pilot scale.

2.5. *Non-Axenic Cultivation in Pilot-Scale (500 L) Culture in Batch and Fed-Batch Modes*
*Appears Feasible*

Cultures were performed in 800 L cylinders, in 500 L culture medium. Three culturing modes were tested: (1) simple batch; (2) batch with feeding after 16 h of cultivation with 30% of the initial media nutrients (glucose, peptone, yeast extract and digestate); or (3) batch with feeding after 16 h with 100% of initial media nutrients. The cultures were harvested after 38 h for (1) and (2) and 62 h for (3). Before 16 h, growth rates appeared similar (Figure 5). Dissolved oxygen measurement showed, in fed cultures after 24 h, values ranging from 2 to 9% of saturation (range values obtained from 11 simple batches; kinetics of dissolved oxygen values during 24 h of 11 simple batches are presented in Figure S1, supplementary material). Stability of dissolved oxygen level after 9 h of culture revealed that the oxygen transfer rate was equal to the oxygen uptake rate. On the other hand, as aeration rate does not change, and as the physical properties of the liquid medium, and its volume can be considered as stable, then the oxygen transfer rate does not change. These stable oxygen transfer rate and uptake rate could induce a linear growth, as already discussed above (Cf. Section 2.1).

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Effect of medium supplementation on biomass production in 500 L reactor cultures (2 replicates). 1X-1, 1X-2: simple batch (no feeding); 1.3X-1, 1.3X-2: feeding with 30% of initial medium organic concentrate after 16 h; 2X-1, 2X-2: feeding with 100% of initial medium organic concentrate after 16 h.</figcaption></figure>

However, as the number of samples collected over time was too scarce, data only
allowed suggesting that a similar linear growth would occur in this pilot experiment after
15–20 h. Productivity and yields are discussed below (Cf. Section 2.9).

2.6. *Bacterial Contamination Assessment in Non-Axenic Cultures*

Bacterial contaminations could be critical for the quality and the security of use of
biomass from non-axenic cultures. A series of 14 cultures at pilot scale were performed from
July 2020 to January 2021. Among these cultures, temperatures ranged from 26 to 30 °C
at 0 h and from 27 to 34 °C at 24 h (Figure S2, supplementary material), pH ranged from
7.1 to 7.8 at 0 h and 3.8 to 6.3 at 24 h (Figure S3, supplementary material), culture volumes
were 450L or 500L, mean inoculum concentration was $2.9 \times 10^6 \pm 5.5 \times 10^5$ cells.mL$^{-1}$ (min = 2.1 × 10$^6$ and max = 3.9 × 10$^6$), and culture harvest was performed after 24 to 62 h of cultivation. The contaminations were first assessed during the process by microscopic observation of fresh cultures. The only detected contaminations were bacteria, which were always present. It must be noted that in some cultures, these contaminants were less present by the end of the culture. Some lab tests were run to study a putative antibacterial effect of culture supernatant, but were unsuccessful. We thus hypothesized a bacterivory impact of *A. mangrovei*, as such a function was previously reported in Thraustochytrids [24]. In order to allow the use of the biomass for feed trials, a microbial risk assessment was performed by ANSES (the food and environmental safety agency of France) on a September production run (and on a previous production run), within the framework of the ALG-AD program. The report established that products were safe for use as feed (Cf. Report Interreg NWE ALG-AD project on Safety analysis digestate and algal biomass produced by the three ALG-AD pilots in Supplementary Files).

In order to further estimate bacterial contamination, a detailed analysis of fatty acids distribution was performed for each culture endpoint sample. A focus was placed on iso 15:0, known as a specific proxy of marine bacterial contamination [25] in favorable growth conditions [26]. The fatty acid composition of an axenic culture (Cf. Section 2.1) was also measured and used as a reference. Results (Figure 6) showed that iso 15:0 was detected at a level below 0.04% in the axenic control. Two cultures exhibited values of 0.15–0.20%, and two cultures showed contamination in the range of 2–4% (December A and B). These cultures concomitantly exhibited an increase in 16:1n-7, suggesting that this fatty acid could be also a proxy for bacterial contamination in our culture conditions. However, even in the most contaminated cultures, the percentage of the main *Aurantiochytrium mangrovei* fatty acids remained within the range of the other cultures and axenic control. Only 15:0 and 17:0, generally present at low concentrations, were reduced when bacterial fatty acids were at their highest.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Relative distribution of the main fatty acids from 14 cultures of 450–500 L (black bars are an average of 3 axenic cultures in carboys on modified YEP medium. Error bars = S.D). January A and January B are the fed batch cultures (+30% of the initial media; see Section 2.5 for details).</figcaption></figure>
Such results illustrate that this non-axenic process can be performed at pilot scale without
significant impact of bacterial contamination on the lipid characteristics of the biomass.

2.7. *Glucose Conversion and Nitrogen Trapping Efficiency*

The most productive pilot-scale culture (fed-batch 2X, 5% digestate, 40 g.L⁻¹ glucose) was assessed for glucose conversion and for N-NH₄ depuration efficiency. Figure 7 shows that from 20 g.L⁻¹ of glucose added at 0 h followed by a second addition of 20 g.L⁻¹ at 20 h, only 3.3 g.L⁻¹ remained at the end of culture, and that it was reduced progressively as biomass concentration increased. From the 0.30 g.L⁻¹ N-NH₄ added in medium, only 0.025 g.L⁻¹ remained at the end of the culture (more than 90% depuration ratio). Ammonia stripping due to aeration could occur in cultures when pH values are above 7 [27]. However, as pH was maintained below pH 6.5 after 14 h (data not shown), the rapid drop in N-NH₄ observed after feeding (16 h) could be assigned to cell trapping of N-NH₄, illustrating the high efficiency of *A. mangrovei* in depuration of N-NH₄ and the use of NH₄ from digestate as part of the nitrogen source in medium.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Biomass (DCW g.L⁻¹), glucose (g.L⁻¹) and nitrogen (g.L⁻¹) concentrations in 500L airlift fed-batch culture; 100% of initial organic substrate was added at 16 h. Mean (n = 2, error bars = S.D.).</figcaption></figure>

2.8. *Lipids and Proteins in Pilot-Scale Production*

In order to assess the characteristics of the biomass produced, analyses were run on a
set of 280, 450 and 500 L pilot cultures, with samples collected over culture time. Figure 8a
shows a correlation between the cell protein concentration and the DCW concentration
in samples. Then, in spite of an important dispersion of values over the curve, it seems
that the proportion of protein in cells remained stable in cultures. In contrast, the lipid
concentration in samples appeared highly variable and could not confidently be related to
DCW concentration (Figure 8b, R²= 0.51). This variability could be explained by a strong
qualitative and quantitative change in lipids over the culture, as illustrated from the 2X
fed-batch culture (Figure 9): both polar and neutral lipids increased (as mg.g⁻¹ DCW)
during the first phase of the culture. After 23 h for polar lipids and 38 h for neutral lipids,
these values started to reduce when the neutral lipids/polar lipids ratio increased. This
could be explained by a shift in the lipid synthesis process from membrane lipids (polar
lipids) produced during the early phase to mostly reserve lipids (neutral lipids) synthesized
during the late phase. The decrease in neutral lipids by the end of the culture could be due
to the catabolism of this energy reserve when glucose was almost fully depleted.
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Cellular protein (a) and total fatty acid (b) concentrations versus biomass concentration, from a data set of 14 kinetics in 450 and 500 L cultures and 3 kinetics in 260 L cultures, over a 6 month period.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Mean concentration (n = 2, error bars = S.D.) over a 2X-fed-batch kinetic in a 500 L culture, of neutral and polar fatty acids in biomass, and percentage (w/w) of neutral lipids over total lipids.</figcaption></figure>

## 2.9. Productivity and Conversion Yields of Substrate to Biomass

In Table 1, we compared the productivities from various experiments. The lowest productivities were observed when no air bubbling was applied. Highest productivities were observed in 500 L cultures. They were 3 times higher or more, as compared to productivities in 10 L cylinders. This resulted from higher cell densities obtained in a shorter time. These large cylinders allow a high water column (1 m), and their large diameter and center-focused bubbling allow an efficient airlift effect and, therefore, a high mixing and oxygen transfer efficiency. Productivities were further increased when fed batch mode was applied, both in axenic and non-axenic cultures, and with full substrate or
glucose-only addition. This could be explained both by a longer part of the culture spent at maximum linear cell production rate (the less productive initial phase of the culture represented a shorter part of the whole culture time), and by a possible positive influence of the high concentration of organic substrates (Y.E and peptone and/or glucose) and/or digestate on this maximum cell production rate. Such positive effect of feeding on cultures was already reported [11,16,17].

<table><caption>**Table 1.** Yields and productivities observed in cultures. Yields were calculated as kg of biomass obtained per kg of initial organic substrate (see Section 4.6.6 for calculation description). When S.D. is indicated, biomass is mean of 3 replicates for flasks and carboys, and 2 replicates for pilot cylinders. DCW = dry cell weight.</caption><thead><tr><th></th><th>Culture Condition</th><th>Initial-Culture Medium</th><th>Feed</th><th>Culture Volume (L)</th><th>Culture Time (Days)</th><th>Biomass Concentration (g.L<sup>-1</sup> ± S.D.)</th><th>Productivity (gDCW.L<sup>-1</sup>. day<sup>-1</sup>)</th><th>Yield</th></tr></thead><tbody><tr><td>Flask-500 mL axenic</td><td>100 rpm</td><td>YEP modified</td><td>/</td><td>0.25</td><td>14</td><td>4.3</td><td>0.31</td><td>0.18</td></tr><tr><td>Flask-1 L axenic</td><td>100 rpm + 0.4 VVM</td><td>Glucose + Salts</td><td>/</td><td>0.8</td><td>5.67</td><td>3.21 ± 0.16</td><td>0.57</td><td>0.32</td></tr><tr><td>Carboy-20 L axenic</td><td>100 rpm + 0.4 VVM</td><td>YEP modified + 2.5% digestate</td><td>/</td><td>6</td><td>2.75</td><td>2.49</td><td>0.91</td><td>0.10</td></tr><tr><td>Cylinder-10 L non-axenic</td><td>Air-lift + 0.4 VVM-large bubbles</td><td>YEP modified + 2.5% digestate</td><td>/</td><td>6</td><td>2.75</td><td>2.82</td><td>1.03</td><td>0.12</td></tr><tr><td>Cylinder-10 L non-axenic</td><td>Air-lift + 0.8 VVM-small bubbles</td><td>YEP modified + 2.5% digestate</td><td>/</td><td>6</td><td>2.75</td><td>3.71</td><td>1.35</td><td>0.15</td></tr><tr><td>Cylinder-800 L non-axenic</td><td>Air-lift + 0.4 VVM</td><td>YEP modified + 2.5% digestate</td><td>/</td><td>500</td><td>1.58</td><td>5.21 ± 0.66</td><td>3.29</td><td>0.22</td></tr><tr><td>Cylinder-800 L non-axenic</td><td>Air-lift + 0.4 VVM</td><td>YEP modified + 2.5% digestate</td><td>+0.3X medium</td><td>500</td><td>1.58</td><td>7.15 ± 0.35</td><td>4.52</td><td>0.23</td></tr><tr><td>Cylinder-800 L non-axenic</td><td>Air-lift + 0.4 VVM</td><td>YEP modified + 2.5% digestate</td><td>+1X medium</td><td>500</td><td>2.57</td><td>14.34 ± 0.72</td><td>5.57</td><td>0.30</td></tr></tbody></table>

Yields were calculated as DCW produced per initial organic substrate added in medium. Then, when substrates were remaining in media, the substrate to biomass conversion efficiencies were underestimated. Despite this, a yield (Table 1) of 0.30 was observed in 2X fed-batch 500 L culture, 0.32 was observed in glucose and salts medium, and 0.28 in YEP-modified medium. The two last yields were obtained when productivities were low. This suggests that the energy budget for long-term cell maintenance, needed when productivity is low, does not influence the yield. With such yields, 3.3 kg of organic substrate is needed to produce 1 kg of DCW Thraustochytrids.

### 3. Discussion

#### 3.1. *Non-Axenic Cultivation*

In our goal to develop an economically realistic industrial process for Thraustochytrids production, a key point was to succeed in running cultivation under non-axenic conditions. To our knowledge, such an approach has not been investigated to date. Microbial sanitary risks evaluation for feed use of our product did not point out any sanitary risks and confirmed the possibility of using this biomass as feed. Furthermore, lipid quality assessment over a 6 month survey of productions (Figure 6) showed that even when bacterial contamination was detected, the fatty acid profile was not significantly impacted. More details on seasonal variations in production can be found in the supplementary file “Interreg NWE ALG-AD project Seasonal variation of algal biomass cultivated using nutrient rich digestate”. Together with positive feed trials on this product [28], these results allow considering further investigations of performance and scale-up of non-axenic industrial production of Thraustochytrids.
## 3.2. Substrates

The use of by-products in order to reduce the costs of process and the environmental impacts of the agro-industry is a major concern. We have shown here that digestate can be used in non-axenic conditions, at pilot scale, in fed-batch mode. The fed-batch mode was already successfully investigated in Thraustochytrids studies [16]. We have shown here that, to avoid a negative impact on growth, digestate from anaerobic digestion (AD) should be added only in small amounts (Figure 3). However, the addition of higher amounts under progressive fed-batch mode should be investigated, as digestate could be progressively metabolized and toxicity avoided. Ultimately, it will allow reducing the usage of yeast extract and peptone to culture medium.

## 3.3. Limiting Factor for Growth and Productivity

We showed (Figures 1 and 4) that the controlling factor for growth was the oxygen availability and that biomass productivity increased when culture volume increased (Table 1). This could be explained by an increase in the oxygen transfer coefficient (Kl.a) when increasing the bioreactor size. We can hypothesize that a further increase in reactor size (and height) could similarly allow a further increase in oxygen transfer, and then in productivity.

The maximum productivity (5.6 g.L⁻¹.day⁻¹), and maximum biomass concentration (14.3 g.L⁻¹) observed in our 500 L fed-batch culture can be compared to reported values in the literature. We calculated, from a recent review from Du et al. [11], Thraustochytrids productivities ranging from 2.5 to 38 g.L⁻¹.day⁻¹ with biomass concentration from 4.8 to 152 g.L⁻¹. It is difficult to comment on such data, as each value is related to a specific process. Furthermore, such values are only a part of the economic assessment. Productivities in our process were drastically increased from lab to pilot scale. Further scale-up should be investigated and could lead to further productivity improvements.

This process shows similarities with the yeast production processes [13,29]. They include non-axenic last step(s) after axenic cultivation, fed-batch mode, oxygen as major limiting factor for growth and high substrate to biomass conversion yield.

The yeast production industry uses low-cost substrates, whereas our process still used small amounts of costly peptone and yeast extract. However, we showed at lab scale (Figure 2 and Table 1) that such compounds are avoidable when maintaining yield.

## 3.4. Lipids

This work was focused on biomass production and process feasibility, and no specific attempt was made to increase lipids or DHA content. DHA productivity ranged from 0.10 to 0.18 mg.L⁻¹.day⁻¹. Fatty acids content in our cultures (i.e., 280 mg.g⁻¹ DCW at maximum and 110 mg.g⁻¹ DCW at harvest, in 2X fed-batch, Figure 9) was far below reported maximum values obtained when, for example, the C/N ratio or dissolved oxygen profiles over the culture were specifically designed [11]. However, harvested biomass from pilot productions was tested in feed trials [28] and proved suitable for fish feed.

In order to further improve productivity and yield together with DHA content, further investigations should include tests on the C/N ratio over the culture phases, eventually coupled with air supply modification during late phase [30–32] and further increases in feeding rate.

## 3.5. AD Digestate Regulation

The use of AD digestate as nitrogen and mineral substrate in industrial cultivation would be an opportunity to reduce adverse impacts such as eutrophication related to excessive agricultural spreading. Unfortunately, recent investigations about the applicable European regulation on AD digestate showed that animal-based digestates could not be used for feed production (see the report of Interreg NWE ALG-AD project "A Regulatory Review on the Use of Digestate to Cultivate Algal Biomass for Animal Feed" included in supplementary file). This limits the potential industrial scale-up of this process to only crop-based digestates. However, other by-products from the food industry, such as
milk derivatives and vegetable processing by-products, could be investigated as potential alternatives to digestate.

This paper showed that the use of a process of the same kind as the widely used low-cost yeast technology, together with the use of low-cost substrates and by-products, could lead to the development of a low-cost process for industrial-scale production of biomass from Thraustochytrids. Low-cost substrates and by-products include dairy waste water, spent brewery yeast, and spent osmotic solution from the candied fruit industry. This process still needs an optimization step for DHA production. For instance, in order to be competitive in the fish feed industry, the production cost of our Thraustochytrid biomass (in dry matter) should range between EUR 3 and 6 per kg according to the targeted fish species and stages. Then, such economically sound development could turn to reality the hope for a reduction in industrial fishing for feed and food.

# 4. Materials and Methods

## 4.1. Culture Media

### 4.1.1. Modified YEP Medium for Lab Experiments and Pilot Precultures

The YEP-modified culture medium was prepared as follows: Yeast extract (VWR™ 84601.5000): 2 g; peptone from casein (VWR™, Rosny-sous-bois, France) peptone from casein 84610.0500): 2 g; sea salt (Sigma™, Saint Louis, MO, USA; Sigma S9883):15 g; dextrose (TITOL chimica SpA™, Pontecchio Polesine, Italy, Glucose anhydre pur CA:50-99-7): 20 g; water: 1 L, autoclaved for 20 min at 120 °C. When added, the digestate obtained from anaerobic digestion of pig manure was stored for settling of large particles, adjusted to pH 7 by H₂SO₄ addition, and autoclaved for 20 min at 120 °C.

### 4.1.2. Modified Guillard and Ryther F/2 Medium for Digestate Comparison Experiment

The medium [33] was designed to allow the supply of nitrogen either from digestate or from NH₄Cl as control. A digestate produced by anaerobic digestion (AD) from agriculture and food by-products was obtained from the Langage AD site (Plymouth, Great Britain), an animal waste-free AD digestate was obtained from University of Ghent (Belgium), and a pig manure digestate was obtained from the Cooperl site (Lamballe, France).

* Digestates and NH₄Cl (297 g.L⁻¹) solution were filter-sterilized through 0.2 µm filter.

* The vitamin solution was prepared as follows: cyanocobalamin: 5 mg; biotin: 5 mg; thiamine HCl: 1 g; distilled water: 500 mL, autoclaved 20 min at 120 °C.

* The trace element solution was prepared as follows: stock solution: CuSO₄·5H₂O: 100 mg, ZnSO₄·7H₂O: 220 mg, CoCl₂·6H₂O: 100 mg, MnCl₂·4H₂O: 1800 mg and Na₂MoO₄·2H₂O: 60 mg. Distilled water: 200 mL. The final trace element solution contained stock solution: 20 mL; Na₂EDTA: 4.36 g; FeCl₃·6H₂O: 3.15 g; distilled water: 980 mL. Solution was filter-sterilized through 0.2 µm.

* Each 800 mL culture was prepared independently as follows:

Twelve grams sea salt (Le Saulnier de Camargue™, Clichy, France) was added in a volume of distilled water calculated to obtain 800 mL of final volume, and autoclaved for 20 min at 120 °C. Solutions were then added in a sterile manner: Na₂HPO₄·7H₂O (56.8 g.L⁻¹, autoclaved for 20 min at 120 °C): 20 mL; MgSO₄·7H₂O (42 g.L⁻¹, filtered through 0.2 µm): 8 mL; glucose (500 g.L⁻¹ autoclaved for 20 min at 120 °C): 16 mL; trace elements solution: 0.8 mL; vitamin solution: 0.8 mL. The nitrogen source was then added: either NH₄Cl (0.297 g.L⁻¹, final concentration): 0.8 mL, Cooperl digestate: 25.8 mL, Langage digestate: 16.7 mL, or U. Ghent digestate: 15.1 mL; pH was finally adjusted to 6.5 in each flask by successive sterile sampling and acetic acid addition.

### 4.1.3. Culture Media for Digestate Concentration Tests

Yeast and peptone concentrate (4 g.L⁻¹ each) in sea water and glucose concentrate (200 g.L⁻¹) in sea water were prepared and autoclaved for 20 min at 120 °C. Digestate (unfiltered) was autoclaved for 20 min at 120 °C. Chemical characteristics of digestate are
provided in Table S1. Solutions were mixed in a 500 mL sterile flask and adjusted with sterile sea water in order to obtain 150 mL of a solution with: 2 g.L⁻¹ peptone; 2 g.L⁻¹ Y.E.; 60 g.L⁻¹ glucose and either 0 (control), 5, 10 or 20% digestate (v/v).

### 4.1.4. Medium for Non-Axenic Culture in 10 L Cylinder Experiments at Lab

Media (6 liters each) for non-axenic cultures in 10 L cylinder and for axenic culture in 20 L carboy were prepared as follows:

Solution 1: Sea salts from Le Saulnier de CamargueTM untreated: 90 g; digestate from Cooperl AD pilot plant, centrifuged in order to eliminate large particles and then filtered through 300 KDa tangential flow filtration membrane: 150 mL (2.5% v/v); tap water: 5.12 L. For carboy culture, the medium was split into five 2 L bottles and autoclaved for 20 min at 120 °C.

Solution 2: Glucose: 150 g; water: 60 mL. Autoclaved for 20 min at 120 °C.

Solution 3: Peptone: 12 g; yeast extract: 12 g; water: 600 mL. Autoclaved for 20 min at 120 °C.

Solutions 1, 2, and 3 were mixed in each cylinder or mixed sterilely in a previously autoclaved (20 min at 120 °C) carboy. Autoclaved (20 min at 120 °C) silicone-based anti-foam was added (final concentration: 1 mL.L⁻¹): Antimousse 426R, LabogrosTM (Coueron, France).

### 4.1.5. Medium for Non-Axenic 400 or 500 L Culture in 800 L Cylinders at Pilot Plant

The medium was prepared by addition of the following solutions:

Solution 1: Filtered digestate (2.5% v/v), containing 48.3 mg nitrogen/L, 2.6 mg phosphorus/L with an N/P ratio of 18.6.

Solution 2: Glucose syrup Isosweet 470 (Tereos™, Moussy-le-vieux, France): 24 g.L⁻¹ of culture medium (as 20 g.L⁻¹ dextrose replacement).

Solution 3: Yeast extract and peptone at 90 g.L⁻¹ each, autoclaved for 20 min at 120 °C (for initial concentration in culture medium: yeast extract 2 g.L⁻¹ and peptone 2 g/L).

Solution 4: Sea salt from Le SaulnierTM (15 g.L⁻¹ of culture medium); Mg SO₄·7H₂O (0.46 g.L⁻¹ of culture medium), tap water as needed.

Silicone-based anti-foam was added as soon as the yeast extract and peptone solutions were poured into the cylinder (final concentration: 1 mL.L⁻¹, Antimousse 426R, LabogrosTM, Coueron, France). For fed-batch feeding, only digestate and organic substrate (Isosweet, peptone and yeast extract) concentrates were added as 30% of initial volumes for 1.3X, and 100% of initial volumes for 2X.

## 4.2. *Thraustochytrid Strain*

The thraustochytrid *A. mangrovei* was obtained from the Roscoff culture collection (Strain RCC 893). Cryovials for pre-culture were prepared as follows:

A 150 mL modified YEP culture was run in a 500 mL round flask on a shaking table at 100 rpm, 25 °C. One milliliter of culture at approximately 2 × 10⁷ cells/mL was centrifuged. The pellet was resuspended in 600 µL of fresh medium and transferred to a cryovial, and 600 µL of autoclaved (20 min at 120 °C) glycerol 40% in water was added. Vials were vortexed and, after 30 min at room temperature, were stored at -80°C.

## 4.3. Pre-Cultures

### 4.3.1. Pre-Cultures for Lab Experiments

Two 500 mL round flasks with 250 mL of modified YEP medium (Cf. Section 4.1.1) were inoculated with the content of a strain cryovial freshly thawed at room temperature. Cultures were maintained for 2 days at 23–25 °C on a shaking table (100 rpm). Two milliliters of this first culture were then used to inoculate four new flasks in the same conditions for 64 h. Flasks were used as inoculum for culture experiments.

### 4.3.2. Pre-Cultures for Pilot-Scale Experiments

Eight liters of modified YEP medium (Cf. Section 4.1.1) were prepared, dispatched as approximately 1.5 L in 2 L bottles, and autoclaved for 20 min at 120 °C. The medium was
then transferred in a previously autoclaved (120 °C, 20 min) 20 L carboy equipped with inlet and outlet 0.2 µm gas filters (Millipore™) and silicone tubing connected to brass end tubing (4.5 mm internal diameter). This brass tubing was in contact with the inside bottom of the carboy in order to allow bubbling in the medium. A sampling bottle (250 mL) was also connected to the silicone tubing in order to allow sterile sampling.

A 250 mL amount of pre-culture (Cf. Section 4.3.1) was inoculated in the 8 L of modified YEP medium in a 20 L carboy. Cultures were run on a shaking table at 100 rpm, with air bubbling at 0.4 volume of air per volume of culture per min (VVM). After 20 h of cultivation, they were transported to the pilot site, maintained under bubbling. Inoculation of pilot cultures occurred within 6 h.

## 4.4. Cultures at Laboratory Scale

### 4.4.1. Scale-Up of Axenic 8 L Cultures in 20 L Carboys

Cultures were prepared following the same protocol as described above (Cf. Section 4.3.2).

Three carboy replicates were prepared and maintained at 25 °C on a shaking table (100 rpm), and air was supplied at 2.4 L.min⁻¹ (0.4 VVM). Samples for analysis were collected under a sterile laminar flow hood. After 3 days, 2 L of a 500 g.L⁻¹ autoclaved (20 min at 120 °C) dextrose solution was added in each carboy (total glucose in final medium: 216 g.L⁻¹).

### 4.4.2. Comparison of Three Sources of Digestate

Cultures in triplicate were run in 800 mL culture medium (Cf. Section 4.1.5), in 1 L flasks with a stainless steel, 2 mm ID tube dispensing 1VVM air through a 0.2 µm inlet filter and maintained on a shaking table at 90 rpm, 25 ± 1 °C. Three different sources of digestate and NH₄Cl as control nitrogen source were tested.

A cryovial (−80 °C) was inoculated in 250 mL medium in a 500 mL flask and run for 72 h, 90 rpm, 25 ± 1 °C. This culture (10 mL) was then inoculated in a new YEP-modified medium flask for 24 h. The resulting culture was used as inoculum for culture (40 mL in each 800 mL culture). Cultures were run in triplicate.

### 4.4.3. Effect of Digestate Concentration on Growth Kinetics in Flask Cultures

Cultures were prepared in 500 mL flasks with 150 mL culture medium (Cf. Section 4.1.2).

Cultures were run in triplicate at 16 ± 1 °C and 90 rpm on a shaking table. In this experiment, cell concentrations were estimated using O.D. measurement at 620 nm. The relationship between cell concentration and O.D. at 620 nm was previously established by the following equation: cell concentration = $1 \times 10^8 \times (\text{O.D. at 620 nm}) - 1 \times 10^7 (R^2 = 0.98)$.

### 4.4.4. Scale-up to Non-Axenic 6 L Cultures in 10 L Cylinders

Cultures (6 liters) were run in 10 L cylinders. Cylinders were home-designed, made of PMMA (100 mm ID) tubing closed with a flat PVC bottom. They were built by AFE consulting, France. The bottom was equipped with inlet tubing for sampling and for air bubbling. Two bubbling systems were tested: a 50 mm PVC plate with 5 holes (1.5 mm) and a 4.5 mm I.D. brass tube. Before each culture, the cylinders were filled up to the top with water, and two drops of 3.3 g of dichloroisocyanuric disodium salt at 81% (N° CAS: 5180-86-0; disinfectant bleach) were added. Cylinders were kept under bubbling for 24 h and then thoroughly rinsed with tap water. The silicone-based anti-foam was added (1 mL.L⁻¹ of culture) immediately after inoculation in cylinders and in one of the carboys.

## 4.5. Cultures at Pilot Scale

The two cylinders were home designed, made of PMMA (800 mm ID, 1700 mm height) tubing closed with a flat 20 mm thick PVC bottom. They were built by AFE consulting, France. Tap water (for process and for cleaning) was supplied by a pump and delivered at the top of the cylinders through a rotating nozzle. Agitation and oxygen supply in each cylinder were provided by airflow bubbling from the bottom of the cylinder, at a rate of
0.4 volume of air per volume of culture per minute (Air-lift system). The air was supplied through a 400 mm diameter ring of 16 mm diameter PVC tubing pierced with 500 holes (1.5 mm diameter). All other necessary inputs were introduced manually at the top of the two cylinders.

Each cylinder was filled up with 500 L of the pilot-scale medium, and each inoculated with 8L culture as prepared above (Cf. Section 4.1.1). The temperature of the culture was maintained between 28 and 30 °C using automated control from an immersed stainless steel electric heater (2 kW for 500 L) and a temperature probe, while pH was maintained above 4.5 by manual pH measurement and addition of 10 N NaOH solution as needed. Three cultures were also run in 330 L cylinders (260 L cultures) in similar conditions. Protein and total lipids analysis were obtained from these three cultures.

## 4.6. Analysis

### 4.6.1. Dry Cell Weight (DCW)

Culture was filtered through Whatman GF/F glass fiber filter, rinsed with the same volume of an ammonium formate solution (0.26 M) and dried at 105 °C for 24 h.

### 4.6.2. Monitoring of Culture Concentration and Cellular Parameters

Concentration and cellular parameters (size, complexity, lipid content) of *A. mangrovei* were measured using an Easy-Cyte Plus 6HT flow cytometer (Guava Merck Millipore®, Darmstadt, Germany) equipped with a 488 nm blue laser, detectors of forward (FSC) and side (SSC) light scatter, and three fluorescence detectors: green (525/30 nm), yellow (583/26 nm) and red (680/30 nm). Cell morphological variables, i.e., forward scatter (FSC) and side scatter (SSC), were used to identify and select the *A. mangrovei* cell population. FSC and SSC give, respectively, information on the relative size and complexity of cells [34]. The flow cytometry measurements were performed on fresh (living) samples. The BODIPY probe (BODIPY 505/515 FL; Molecular Probes, Invitrogen, Eugene OR, USA, final concentration of 10 mM), which stains lipid droplets/bodies within Thraustochytrid cells, was used as a proxy of lipid reserves [34]. The green fluorescence emitted was proportional to the quantity of lipid reserve present in the cells. Measurements were performed at a flow rate of 59 µL.min⁻¹. Concentration of Thraustochytrid was given as cells per mL, and cellular parameters (FSC, SCC, and lipid reserve) were expressed as fluorescence intensity in arbitrary units (A.U).

### 4.6.3. Nutrient Content Analysis

Glucose concentration was measured using the DNS colorimetric method. Dissolved ammonia in medium was measured following Hansen and Koroleff [35] and Eppley et al. [36]. Dissolved phosphate in medium was measured following Murphy and Riley [37].

### 4.6.4. Lipid Content Analysis

*   **Lipid Extraction:**

Fifty milligrams of freeze-dried Thraustochytrid biomass were extracted in 6 mL of solvent mixture (CHCl₃:MeOH, 2:1, v:v) directly added into glass vials [38]. Extracts were flushed with nitrogen gas, vortexed, and sonicated for 10 min to ensure complete lipid extraction and stored at -20 °C overnight for further analysis.

*   **Separation of polar lipids (PL) and neutral lipids (NL):**

To fractionate NL and PL, 1 mL of total lipid (TL) extract was evaporated to dryness under nitrogen, recovered with 3 washes of 0.5 mL of CHCl₃:MeOH (98:2, v:v; final volume 1.5 mL) and spotted at the top of a silica gel column (40 mm × 4 mm, silica gel 60A 63–200 µm rehydrated with 6% H₂O, 70–230 mesh, Sigma-Aldrich, Darmstadt, Germany). The NL was then eluted using CHCl₃:MeOH (98:2, v:v; 10 mL) and PL fraction with methanol (20 mL). Both were collected in glass vials containing an internal standard (C23:0, 2.3 µg). Lipid fractions were then stored at -20 °C under nitrogen atmosphere until further analysis.
* Fatty acid analysis in total lipids (TL), polar lipids (PL) and neutral lipids (NL):

Polar and neutral fractions or total lipid extracts were dried under vacuum with an evaporator (Genevac). Dried lipid fractions were saponified with 1 mL of KOH:MeOH (0.5 M) heated for 30 min at 80 °C, then they were transesterified with 800 µL of MeOH:H₂SO₄ (3.4%; *v*/*v*) heated at 100 °C for 10 min. The fatty acid methyl esters (FAME) formed were recovered in hexane. FAME of the NL of Thraustochytrid biomass were separated from other unwanted compounds (e.g., sterols and alcohols, potentially present in the NL fraction) using high performance liquid chromatography (HPLC) equipped with two columns (LiChrospher Si 60 and LiChrospher 100 DIOL, both 5 µm). A Dionex HPLC system (P680 pump AS-100 auto sampler, UVD170U UV detector with deuterium lamp, Foxy fraction collector), was used. Details on the analytical methods (i.e., solvent proportions, flow rate) can be found in [39]. The purified FAME were recovered in new vials for gas chromatography analysis.

FAME composition was analyzed by gas chromatography coupled to a flame ionization
detector (GC-FID; Varian CP8400 gas chromatograph, Agilent). Samples (2 µL) were
injected at 250 °C in splitless mode at an oven temperature of 60 °C, with hydrogen as the
carrier gas. Detector temperature was set at 280 °C. The GC was equipped with a ZB-WAX
column (30 m in length, 0.25 mm internal diameter, 0.25 µm film thickness, Phenomenex).
The oven temperature was raised to 150 °C at 50 °C.min⁻¹, to 170 °C at 3.5 °C.min⁻¹,
to 185 °C at 1.5 °C.min⁻¹, to 225 °C at 2.4 °C.min⁻¹ and then to 250 °C at 5.5 °C.min⁻¹.
FAME were identified by comparing their retention times to those of external commercial
standard mixtures (S37 FAME Mix, PUFA n°1, and PUFA n°3, Supelco) using the software
Galaxie 1.9.3.2 (Agilent). FAME peak area was converted into µg of FA based on the peak
area of the internal standard C23:0. Concentrations of total neutral and polar lipids were
expressed in mg.g⁻¹ of biomass dry weight, and FA compositions of pilot-scale cultures
were expressed in percentage of total fatty acids.

4.6.5. Protein Analysis

Two milliliters of culture (cells+ supernatant) and of supernatant previously filtered through 0.2 µm were thawed and hydrolyzed with 1 N NaOH for 1 h at room temperature. Two hundred microliters of resulting protein extracts were analyzed following the Lowry method, using a BIO-RAD™ kit (Ref: 5000116 DC protein assays reagent package). Standard curve was prepared using bovine serum albumin between 0.1 and 1.5 mg.mL⁻¹. Whenever necessary, protein extracts were diluted to fit within the standard curve.

4.6.6. Conversion Yield

Yields were calculated as follows: [biomass concentration dry cell weight (g.L⁻¹)]/[sum of initial organic substrate concentrations (g.L⁻¹)].

4.6.7. On Line Measurements

For pilot scale, the temperature, pH and polarographic dissolved oxygen probes,
connected on a multi-parameters system KM 3000, were obtained from Sensortechnik
Meinsberg GmbH, Germany.

For 10 L cylinders in the lab, temperature, pH and polarographic oxygen probes
from IngoldTM were connected to an ApplikonTM controller ADI 1030 (Getinge AB, Göte-
borg, Sweden).

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/md20080499/s1. It includes three reports named: Interreg NWE ALG-AD project Seasonal variation of algal biomass cultivated using nutrient rich digestate; Interreg NWE ALG-AD project Safety analysis digestate and algal biomass produced by the three ALG-AD pilots; Interreg NWE ALG-AD project A Regulatory Review on the Use of Digestate to Cultivate Algal Biomass for Animal Feed. One table (S1) and three figures (S1, S2, and S3) are also in a supplementary file with the following captions: Table S1: Chemical characterization of digestate from COOPERL organic biological waste (pig manure) in Lamballe (Brittany, France); Figure S1: Percentages of
dissolved oxygen of *A. mangrovei* cultivated in 800 L cylinder at 0, 3, 6, 9, and 24 h (n=11 batches, error bars= S.D.); Figure S2: Mean values of *A. mangrovei* culture temperature in 800 L cylinder at 0, 3, 6, 9 and 24 h (n=11 batches, error bars= S.D.); Figure S3: Mean values of *A. mangrovei* culture pH in 800 L cylinder at 0, 3, 6, 9 and 24 h (n=11 batches, error bars= S.D.).

**Author Contributions:** Conceptualization and methodology, D.d.l.B., L.C., M.V., P.S.; investigation, D.d.l.B., L.C., M.V., M.G., T.M., T.V., N.G., F.L.G., A.B. (Antoine Bideau), N.L.G., A.B. (Adeline Bidault), C.L., P.S.; writing—original draft preparation, D.d.l.B.; writing—review and editing, C.L., L.C., M.V., P.S.; project administration and funding acquisition, D.d.l.B., P.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the European Union, "Interreg North-West Europe" under the acronym ALG-AD.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not Applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** Authors thank Barbara Clément-Larosière and all his team at DENITRAL for its technical and material support at pilot plant.

**Conflicts of Interest:** The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Food and Agriculture Organization of the United Nations. *The State of World Fisheries and Aquaculture 2020: Sustainability in Action*; Food and Agriculture Organization of the United Nations FAO: Rome, Italy, 2020; ISBN 978-92-5-132692-3.

2. Tocher, D.R. Omega-3 Long-Chain Polyunsaturated Fatty Acids and Aquaculture in Perspective. *Aquaculture* **2015**, *449*, 94–107. [CrossRef]

3. Arts, M.T.; Ackman, R.G.; Holub, B.J. "Essential Fatty Acids" in Aquatic Ecosystems: A Crucial Link between Diet and Human Health and Evolution. *Can. J. Fish. Aquat. Sci.* **2001**, *58*, 122–137. [CrossRef]

4. Amminger, G.P.; Schäfer, M.R.; Schlögelhofer, M.; Klier, C.M.; McGorry, P.D. Longer-Term Outcome in the Prevention of Psychotic Disorders by the Vienna Omega-3 Study. *Nat. Commun.* **2015**, *6*, 7934. [CrossRef] [PubMed]

5. Crawford, M.A.; Broadhurst, C.L. The Role of Docosahexaenoic and the Marine Food Web as Determinants of Evolution and Hominid Brain Development: The Challenge for Human Sustainability. *Nutr. Health* **2012**, *21*, 17–39. [CrossRef] [PubMed]

6. Food and Agriculture Organization of the United Nations. *The State of World Fisheries and Aquaculture*; Food and Agriculture Organization of the United Nations FAO: Rome, Italy, 2014.

7. Lotze, H.K.; Tittensor, D.P.; Bryndum-Buchholz, A.; Eddy, T.D.; Cheung, W.W.L.; Galbraith, E.D.; Barange, M.; Barrier, N.; Bianchi, D.; Blanchard, J.L.; et al. Global Ensemble Projections Reveal Trophic Amplification of Ocean Biomass Declines with Climate Change. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 12907–12912. [CrossRef]

8. Fossier Marchan, L.; Lee Chang, K.J.; Nichols, P.D.; Mitchell, W.J.; Polglase, J.L.; Gutierrez, T. Taxonomy, Ecology and Biotechnological Applications of Thraustochytrids: A Review. *Biotechnol. Adv.* **2018**, *36*, 26–46. [CrossRef] [PubMed]

9. Patel, A.; Karageorgou, D.; Katapodis, P.; Sharma, A.; Rova, U.; Christakopoulos, P.; Matsakas, L. Bioprospecting of Thraustochytrids for Omega-3 Fatty Acids: A Sustainable Approach to Reduce Dependency on Animal Sources. *Trends Food Sci. Technol.* **2021**, *115*, 433–444. [CrossRef]

10. Orozco Colonia, B.S.; Vinícius de Melo Pereira, G.; Soccol, C.R. Omega-3 Microbial Oils from Marine Thraustochytrids as a Sustainable and Technological Solution: A Review and Patent Landscape. *Trends Food Sci. Technol.* **2020**, *99*, 244–256. [CrossRef]

11. Du, F.; Wang, Y.-Z.; Xu, Y.-S.; Shi, T.-Q.; Liu, W.-Z.; Sun, X.-M.; Huang, H. Biotechnological Production of Lipid and Terpenoid from Thraustochytrids. *Biotechnol. Adv.* **2021**, *48*, 107725. [CrossRef] [PubMed]

12. Liu, W.; Mao, W.; Zhang, C.; Liu, L.; Zhang, Z.; Guo, C.; Lin, J. Effective and Economic Microbial Lipid Biosynthesis for Biodiesel Production by Two-Phase Whole-Cell Biocatalytic Process. *J. Clean. Prod.* **2021**, *298*, 126798. [CrossRef]

13. Loiez, A. *Production de La Levure de Panification Par Biotechnologie*; Techniques de L'ingénieur Bioprocédés et Bioproductions 2003, Base Documentaire: TIP140WEB; Techniques de L'ingénieur: Paris, France, 2003.

14. Gahlawat, G.; Kumari, P.; Bhagat, N.R. Technological Advances in the Production of Polyhydroxyalkanoate Biopolymers. *Curr. Sustain. Renew. Energy Rep.* **2020**, *7*, 73–83. [CrossRef]

15. Lei, M.; Peng, X.; Sun, W.; Zhang, D.; Wang, Z.; Yang, Z.; Zhang, C.; Yu, B.; Niu, H.; Ying, H.; et al. Nonsterile L-Lysine Fermentation Using Engineered Phosphite-Grown Corynebacterium Glutamicum. *ACS Omega* **2021**, *6*, 10160–10167. [CrossRef] [PubMed]

16. Qu, L.; Ren, L.-J.; Sun, G.-N.; Ji, X.-J.; Nie, Z.-K.; Huang, H. Batch, Fed-Batch and Repeated Fed-Batch Fermentation Processes of the Marine *Thraustochytrid Schizochytrium* sp. For Producing Docosahexaenoic Acid. *Bioprocess Biosyst. Eng.* **2013**, *36*, 1905–1912. [CrossRef]
17. Ye, H.; He, Y.; Xie, Y.; Sen, B.; Wang, G. Fed-Batch Fermentation of Mixed Carbon Source Significantly Enhances the Production of Docosahexaenoic Acid in *Thraustochytriidae* sp. PKU#Mn16 by Differentially Regulating Fatty Acids Biosynthetic Pathways. *Bioresour. Technol.* **2020**, *297*, 122402. [CrossRef] [PubMed]

18. Silva, D.; Villarroel, M.P.; Roa, A.L.; Quilodrán, B.H. Use of Waste from Agroindustrial Sources as Substrate for Polyunsaturated Fatty Acids Production by Thraustochytrium Kinney VAL-B1. *Int. J. Eng. Res. Afr.* **2017**, *33*, 50–55. [CrossRef]

19. Yamasaki, T.; Aki, T.; Shinozaki, M.; Taguchi, M.; Kawamoto, S.; Ono, K. Utilization of Shochu Distillery Wastewater for Production of Polyunsaturated Fatty Acids and Xanthophylls Using Thraustochytrid. *J. Biosci. Bioeng.* **2006**, *102*, 323–327. [CrossRef] [PubMed]

20. Smerilli, M.; Neureiter, M.; Wurz, S.; Haas, C.; Frühauf, S.; Fuchs, W. Direct Fermentation of Potato Starch and Potato Residues to Lactic Acid by Geobacillus Stearothermophilus under Non-Sterile Conditions. *J. Chem. Technol. Biotechnol.* **2015**, *90*, 648–657. [CrossRef] [PubMed]

21. Monod, J. The Growth of Bacterial Cultures. *Annu. Rev. Microbiol.* **1949**, *3*, 371–394. [CrossRef]

22. Garcia-Ochoa, F.; Gomez, E. Bioreactor Scale-up and Oxygen Transfer Rate in Microbial Processes: An Overview. *Biotechnol. Adv.* **2009**, *27*, 153–176. [CrossRef]

23. Collos, Y.; Harrison, P.J. Acclimation and Toxicity of High Ammonium Concentrations to Unicellular Algae. *Mar. Pollut. Bull.* **2014**, *80*, 8–23. [CrossRef]

24. Raghukumar, S. Bacterivory: A Novel Dual Role for Thraustochytrids in the Sea. *Mar. Biol.* **1992**, *113*, 165–169. [CrossRef]

25. Dalsgaard, J.; St. John, M.; Kattner, G.; Müller-Navarra, D.; Hagen, W. Fatty Acid Trophic Markers in the Pelagic Marine Environment. In *Advances in Marine Biology*; Academic Press: Cambridge, MA, USA, 2003; Volume 46, pp. 225–340.

26. Najdek, M.; Debobbis, D.; Mioković, D.; Ivančić, I. Fatty Acid and Phytoplankton Compositions of Different Types of Mucilaginous Aggregates in the Northern Adriatic. *J. Plankton Res.* **2002**, *24*, 429–441. [CrossRef]

27. Bonmatí, A.; Flotats, X. Air Stripping of Ammonia from Pig Slurry: Characterisation and Feasibility as a Pre- or Post-Treatment to Mesophilic Anaerobic Digestion. *Waste Manag.* **2003**, *23*, 261–272. [CrossRef]

28. Soudant, P.; Ventura, M.; Chauchat, L.; Guerreiro, M.; Le Grand, F.; Simon, V.; Collet, S.; Zambonino-Infante, J.-L.; Le Goïc, N.; Lambert, C.; et al. Evaluation of Aurantiochytrium mangrovei Biomass Grown on Digestate as a Feed Ingredient of Sea Bass, Dicentrarchus labrax, Juveniles and Larvae. 2022; submitted.

29. Lisičar, J.; Scheper, T.; Barbe, S. Turning Industrial Baker's Yeast Manufacture into a Powerful Zero Discharge Multipurpose Bioprocess. *Ind. Biotechnol.* **2017**, *13*, 184–191. [CrossRef]

30. Jakobsen, A.N.; Aasen, I.M.; Josefsen, K.D.; Strøm, A.R. Accumulation of Docosahexaenoic Acid-Rich Lipid in *Thraustochytrid Aurantiochytrium* sp. Strain T66: Effects of N and P Starvation and O₂ Limitation. *Appl. Microbiol. Biotechnol.* **2008**, *80*, 297–306. [CrossRef]

31. Chi, Z.; Liu, Y.; Frear, C.; Chen, S. Study of a Two-Stage Growth of DHA-Producing Marine Algae Schizochytrium Limacinum SR21 with Shifting Dissolved Oxygen Level. *Appl. Microbiol. Biotechnol.* **2009**, *81*, 1141–1148. [CrossRef] [PubMed]

32. Ling, X.; Guo, J.; Liu, X.; Zhang, X.; Wang, N.; Lu, Y.; Ng, I.-S. Impact of Carbon and Nitrogen Feeding Strategy on High Production of Biomass and Docosahexaenoic Acid (DHA) by *Schizochytrium* sp. LU310. *Bioresour. Technol.* **2015**, *184*, 139–147. [CrossRef] [PubMed]

33. Guillard, R.R.; Ryther, J.H. Studies of Marine Planktonic Diatoms. I. Cyclotella Nana Hustedt, and *Detonula confervacea* (Cleve) Gran. *Can. J. Microbiol.* **1962**, *8*, 229–239. [CrossRef] [PubMed]

34. Lelong, A.; Jolley, D.F.; Soudant, P.; Hégaret, H. Impact of Copper Exposure on *Pseudo-Nitzschia* spp. Physiology and Domoic Acid Production. *Aquat. Toxicol.* **2012**, *118–119*, 37–47. [CrossRef]

35. Hansen, H.P.; Koroleff, F. Determination of Nutrients. In *Methods of Seawater Analysis: Third, Completely Revised and Extended Edition*; John and Wiley and Sons: Hoboken, NJ, USA, 2007; pp. 159–228. ISBN 978-3-527-61398-4.

36. Eppley, R.W.; Coatsworth, J.L.; Solórzano, L. Studies of Nitrate Reductase in Marine Phytoplankton. *Linnol. Oceanogr.* **1969**, *14*, 194–205. [CrossRef]

37. Murphy, J.; Riley, J.P. A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. *Anal. Chim. Acta* **1962**, *27*, 31–36. [CrossRef]

38. Sardenne, F.; Bodin, N.; Metral, L.; Crottier, A.; Le Grand, F.; Bideau, A.; Brisset, B.; Bourjea, J.; Saraux, C.; Bonhommeau, S.; et al. Effects of Extraction Method and Storage of Dry Tissue on Marine Lipids and Fatty Acids. *Anal. Chim. Acta* **2019**, *1051*, 82–93. [CrossRef] [PubMed]

39. Marty, Y.; Soudant, P.; Perrotte, S.; Moal, J.; Dussauze, J.; Samain, J.F. Identification and Occurrence of a Novel Cis-4,7,10,Trans-13-Docosatetraenoic Fatty Acid in the Scallop *Pecten maximus* (L.). *J. Chromatogr. A* **1999**, *839*, 119–127. [CrossRef]