# Statistical Optimization of the Physico-Chemical Parameters for Pigment Production in Submerged Fermentation of Talaromyces albobiverticillius 30548

Mekala Venkatachalam, Alain Shum-Chéong-Sing, Laurent Dufossé, Mireille Fouillaud

► To cite this version:

Mekala Venkatachalam, Alain Shum-Chéong-Sing, Laurent Dufossé, Mireille Fouillaud. Statistical Optimization of the Physico-Chemical Parameters for Pigment Production in Submerged Fermentation of Talaromyces albobiverticillius 30548. Microorganisms, 2020, 8 (5), pp.711. 10.3390/microorganisms8050711. hal-04145603

HAL Id: hal-04145603

https://hal.science/hal-04145603

Submitted on 29 Jun 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

# Statistical Optimization of the Physico-Chemical Parameters for Pigment Production in Submerged Fermentation of *Talaromyces albobiverticillius* 30548

Mekala Venkatachalam<fnref n="1" /> ID, Alain Shum-Chéong-Sing<fnref n="1" />, Laurent Dufossé<fnref n="1" />,<fnref n="2" /> ID and
Mireille Fouillaud <fnref n="1" />,<fnref n="2" />,* ID

Received: 5 April 2020; Accepted: 8 May 2020; Published: 11 May 2020

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

check for
updates

**Abstract:** *Talaromyces albobiverticillius* 30548 is a marine-derived pigment producing filamentous fungus, isolated from the La Réunion island, in the Indian Ocean. The objective of this study was to examine and optimize the submerged fermentation (SmF) process parameters such as initial pH (4–9), temperature (21–27 °C), agitation speed (100–200 rpm), and fermentation time (0–336 h), for maximum production of pigments (orange and red) and biomass, using the Box–Behnken Experimental Design and Response Surface Modeling (BBED and RSM). This methodology allowed consideration of multifactorial interactions between a set of parameters. Experiments were carried out based on the BBED using 250 mL shake flasks, with a 100 mL working volume of potato dextrose broth (PDB). From the experimental data, mathematical models were developed to predict the pigments and biomass yields. The individual and interactive effects of the process variables on the responses were also investigated (RSM). The optimal conditions for maximum production of pigments and biomass were derived by the numerical optimization method, as follows—initial pH of 6.4, temperature of 24 °C, agitation speed of 164 rpm, and fermentation time of 149 h, respectively.

**Keywords:** *Talaromyces albobiverticillius* 30548; submerged fermentation; pigments; biomass; optimization; Box–Behnken experimental design; response surface modeling

## 1. Introduction

In recent years, development of alternative sources for the production of natural pigments has been focused on overcoming the unlimited usage of synthetic pigments, which are found to be hazardous to the human health and environment [1]. Natural colorants obtained from plants, animals, and microbial sources are alternatives to synthetic pigments [2,3]. Nowadays, the shift towards natural colorants is of great interest worldwide and is gaining its industrial importance in food, cosmetics, pharmaceuticals, paint, or dying applications. Microorganisms including bacteria, algae, yeasts, and filamentous fungi from marine or terrestrial origins are capable of producing natural dyes, which pave the way to develop the industrial production of natural compounds [4–6]. Indeed, microbial fermentations present the crucial advantages of independence from geographic or climatic constraints. Additionally, massive production of pigments using bioreactors can be totally controlled and optimized. The use of fungi for the production of commercially important products has been increasing rapidly over the past decades [4,7]. Taxonomically, it represents a large group, often prevailing in seawater, rocks, sediments,

<footnote n="1">¹ Laboratoire de Chimie et Biotechnologie des Produits Naturels-EA 2212, Université de la Réunion, 15 Avenue René Cassin, CS 92003, CEDEX 9, F-97744 Saint-Denis, Ile de la Réunion, France; mekalavenkat@gmail.com (M.V.); alain.shum@univ-reunion.fr (A.S.-C.-S.); laurent.dufosse@univ-reunion.fr (L.D.)</footnote>

<footnote n="2">² Ecole Supérieure d'Ingénieurs Réunion Océan Indien-ESIROI, 2 Rue Joseph Wetzell, F-97490 Sainte-Clotilde, Ile de la Réunion, France</footnote>

<footnote>* Correspondence: mireille.fouillaud@univ-reunion.fr; Tel.: +262-262-483-363</footnote>
sand, soils, and diverse marine macro-organisms [5]. In general, the colorants produced by fungi are associated with increased yields, and huge array of compounds due to the abundance and vast diversity of fungal strains in the environment [8]. Majority of the pigments produced by fungi are quinones, flavonoids, melanins, and azaphilones, which belong to the aromatic polyketide chemical group [9–11] and have been widely described for their medicinal uses or potential utilization as dyes [12,13]. More specifically, filamentous fungi are considered to be a promising source of natural pigments. Moreover, some of the produced molecules also possess bioactive properties like anti-cancer, immunomodulatory, anti-proliferative, antibiotic properties, and so on [14].

Fungi isolated from the marine environments reveal a plethora of known or new compounds with still relatively unexplored bioactivities [15–17]. Marine fungi of several genera, such as *Aspergillus*, *Cladosporium*, *Eurotium*, *Monascus*, *Penicillium*, *Talaromyces*, and *Trichoderma* are potent pigment producers of different classes, with wide range of hues [18–22]. For fungal pigment production, submerged fermentation (SmF) appears to be a convenient and economical fermentation method. Indeed, SmF provides the space to produce pigments, owing to their deep commitment to sophisticated control devices and easy monitoring methods during production [23–25]. Besides this, SmF processes can also be affected by various abiotic factors, such as pH of the medium, temperature, agitation speed (gas diffusion), and fermentation time. The mentioned abiotic factors create a strong pressure on the pattern of fungal growth that brings out the interesting secondary metabolites like pigments [26–29].

Over the past few years, a number of studies identified that several species of *Penicillium* and *Talaromyces* produce *Monascus*-like red pigments without the co-production of mycotoxin citrinin, which is produced in *Monascus* spp. fermentations [30]. Some species of *Talaromyces* secrete large amounts of red pigment, notably species such as *Talaromyces purpurogenus*, *T. albobiverticillius*, *T. marneffei*, and *T. minioluteus* (often referenced in the literature under earlier *Penicillium* names). Isolates identified as *T. purpurogenus* have been considered to be industrially interesting for their red pigment production. Unfortunately, the pigment production comes along with the co-production of mycotoxins such as rubratoxin a and b and luteoskyrin. Whereas, a novel strain, *Talaromyces atroroseus* CBS 133442 produces diffusible red pigments, without any mycotoxin production [19]. Certain strains of *Talaromyces albobiverticillius* were studied and proposed to produce large amount of red pigments, thus, can be used for coloring food [31]. However, considering *T. albobiverticillius* 30548, collecting large amount of pigments is a real challenge. First, at a lab scale it allows the purification, characterization, and structural elucidation of all molecules produced, including undesirable compounds, if any. In a second step, optimizing the production of compounds of interest is crucial to determining the ultimate potential of the fungus, for a successful production of pigments at an industrial scale.

The present study investigated and optimized the SmF process crucial parameters such as initial pH (4-9), temperature (21-27 °C), agitation speed (100-200 rpm), and fermentation time (24-336 h) for the maximum production of pigments (orange and red) and biomass by *Talaromyces albobiverticillius* 30548. To achieve this goal, Box-Behnken Experimental Design (BBED) and response surface modeling (RSM) were used. BBED is a methodology based on a three interlocking 2² factorial experimental design. It is also described as a spherical revolving surface surrounding a cube constructed from the responses of the center and middle points of the experimental design [32–37]. This experimental design allows consideration of the multifactorial interactions between a set of crucial parameters. Derived from this mathematical approach, it is possible to optimize the culture conditions for maximum production of pigments or biomass and its interactive effects [38]. This is the first report to investigate and optimize the SmF process parameters by BBED and RSM, for maximal pigment and biomass production, using the marine-derived *T. albobiverticillius* 30548.
## 2. Materials and Methods

### 2.1. *Microorganism and Maintenance*

The filamentous fungus used in this study was isolated from the sediments collected from the marine environment of the La Réunion island, Indian Ocean (GPS coordinates: 21°06′22.11″ S, 55°14′ 15.78″E) [39]. From taxonomic identification of the fungal biodiversity center (Westerdijk Fungal Biodiversity Institute, Utrecht, The Netherlands), the fungus was identified as *Talaromyces albobiverticillius* and named as *T. albobiverticillius* 30548 (GenBank accession number MK937814). The isolated fungus was cultured in potato dextrose agar plates (PDA, BD-Difco, Sparks, NV, USA) and maintained at 4 °C, as well as sub-cultured at regular intervals for the experiments.

### 2.2. *Preparation of Pre-Culture*

Laboratory grade potato dextrose broth (PDB) (BD-Difco)) was used to prepare the seed medium. For inoculum preparation, 100 mg of mycelial spores was scrapped from actively grown culture on potato dextrose agar petri plate (PDA), (7 days old) (composition: PDB added with 2% Bacto agar (BD-Difco) and transferred into 250 mL Erlenmeyer flask containing 80 mL of sterilized PDB medium (pre-culture). The flasks were incubated at 25 °C for 48 h in an orbital shaker, at 200 rpm (Multitron Pro, Infors HT, Bottmingen, Switzerland) before being transferred to the main fermentation medium [19].

### 2.3. *Fermentation Conditions*

After 48 h of growth, the pre-culture broth was centrifuged at 8000 rpm for 6 min at room temperature (Sigma 4K15, Sigma Laborzentrifugen GmbH, Osterode am Harz, Germany) to separate the mycelia and the culture filtrate. The harvested mycelium (200 mg) was transferred as seed inoculum into a 500 mL Erlenmeyer flask containing 200 mL of sterile fermentation medium (PDB). Twenty-nine (29) series of submerged fermentations were carried out according to the experimental design (BBED) presented in Table 1. The BBD matrix (Table 2) allowed to fix a set of selected values for initial pH, temperature, agitation speed, and fermentation time, for each fermentation. The flasks were incubated in a controlled orbital shaker with the fixed conditions (Multitron Pro, Infors HT, Bottmingen, Switzerland) during the whole culture period. All experiments were performed in triplicates for each condition. At predefined intervals of 24, 168, 336 h (Table 1), fermented broth samples were taken out from each culture flask, in order to measure the concentration of orange and red pigments, as well as the increase in biomass weight.

### 2.4. *Spectrophotometric Quantification of Extracellular Pigments*

Five milliliters of the fermented culture broth were sampled at predefined regular intervals and filtered through pre-weighed nylon mesh of 48 µm pore size (Nitex 07401, Sefar AG, Heiden, Switzerland) to separate the mycelia (Table 1). The concentration of extracellular pigments was determined using UV-visible spectrophotometer (UV-VIS spectrophotometer UV-1800, Shimadzu, Tokyo, Japan) at 470 (orange pigments) and 500 nm (red pigments), as the color in the extracellular culture filtrate of *T. albobiverticillius* 30548 was intense dark red [39]. Quinizarin (Sigma-Aldrich, St. Louis, MO, USA) and Red Yeast rice pigments (RYrp) (Wuhan Jiacheng Biotechnology Co., Ltd., Wuhan, China) were used as the reference standards, since the color and absorbance profile of the pigments produced by *T. albobiverticillius* 30548 were similar to those of these two compounds, at a wavelength of 470 and 500 nm. Thus, the OD values at 470 and 500 nm were interpolated using the standard curve equations of the two chosen standards, respectively (Figures A1 and A2 in Appendix A), to obtain the pigment concentration in terms of gram per liter equivalents (g/L equivalent). The concentrations of red and orange pigments (g/L equivalent) were consequently calculated using Equations (1) and (2), as shown below:
$$ \text{Orange pigment concentration, } C\left(\frac{g}{L}\text{equiv.}\right) = \frac{\text{Abs}_{470} - 0.1102}{0.0148} * 0.001, \qquad (1) $$

$$ \text{Red pigment concentration, } C\left(\frac{g}{L}\text{equiv.}\right) = \frac{\text{Abs}_{500} - 0.0968}{0.0069} * 0.001, \qquad (2) $$

## 2.5. Dry Biomass Concentration

The amount of wet biomass was obtained through filtration using 48 m Nitex filter cloth (SEFAR AG, Switzerland) and was weighed precisely using an analytical balance. The filters with wet biomass was dried in a hot air oven (SNB 100, Memmert GMBH, Schwabach, Germany) at 105 °C for 17 h. The dried filters were weighed after letting them cool in a desiccator for 30 min to bring it down to room temperature. The amount of dry biomass in the culture broth was inferred from the following, Equation (3):

$$ [C_{DB}]_{\frac{g}{L}} = \frac{(W_2 - W_1)}{V}, \qquad (3) $$

where W₂ is the weight of the filter with biomass obtained after drying (g), W₁ is the weight of the corresponding empty filter (g), and V is the volume of sample (L).

## 2.6. Experimental Design and Statistical Analysis (BBED and RSM)

Box-Behnken experimental design (BBED) and response surface modeling (RSM) with four factors (pH, temperature, agitation speed, and fermentation time) at three levels were used in this study to investigate and optimize the effect of process parameters on the production of pigments (orange and red pigments) and fungal biomass. The process parameters and their ranges were chosen, based on preliminary “one-factor-at-a-time” studies (Table 1) (data not shown). Each independent variable was coded at three levels as -1 (low), 0 (central point), and +1 (high), which were—pH (4–9), temperature (21–27 °C), agitation speed (100–200 rpm), and fermentation time (24–336 h). The process variables were denoted as X₁, X₂, X₃, and X₄ and the actual value of process parameters were converted into the uncoded form, using the following Equation (4) [40].

$$ x = \frac{X - ((X_{\text{max}} + X_{\text{min}})/2)}{(X_{\text{max}} - X_{\text{min}})/2}, \qquad (4) $$

<table><caption>**Table 1.** Coded and actual values of the variables for the four factor Box-Behnken experimental design.</caption><thead><tr><th rowspan="2">Variables</th><th rowspan="2">Symbol</th><th colspan="3">Coded and Actual Values</th></tr><tr><th>-1</th><th>0</th><th>+1</th></tr></thead><tbody><tr><td>pH</td><td>X<sub>1</sub></td><td>4</td><td>6.5</td><td>9</td></tr><tr><td>Temperature (°C)</td><td>X<sub>2</sub></td><td>21</td><td>24</td><td>27</td></tr><tr><td>Agitation speed (rpm)</td><td>X<sub>3</sub></td><td>100</td><td>150</td><td>200</td></tr><tr><td>Fermentation time (h)</td><td>X<sub>4</sub></td><td>24</td><td>168</td><td>336</td></tr></tbody></table>

In total, twenty-nine experiments were carried out including five center points. Design-Expert® Software (Version 9, Stat-Ease Inc., Minneapolis, MN, USA) was used in this study to construct the experimental design and statistically analyze the experimental data.

The obtained experimental results should fit well into the empirical second-order polynomial model (Equation (5)). The second-order polynomial equation correlated the relationship between independent variables and the responses. The mathematical form of the equation is given below:

$$ Y = \beta_0 + \sum_{j=1}^{k} \beta_j x_j + \sum_{j=1}^{k} \beta_{jj} x_j^2 + \sum_{i < j = 2}^{k} \beta_{ij} x_i x_j + e_i, \qquad (5) $$
where Y is the response; $x_i$ and $x_j$ are variables (i and j range from 1 to k); $\beta_0$ is the model intercept coefficient; $\beta_j$, $\beta_{jj}$ and $\beta_{ij}$ are interaction coefficients of the linear, quadratic, and second-order terms, respectively; k is the number of independent parameters (k = 4 in this study); and $e_i$ is the error [41–43].

The sequential model sum of squares and model summary statistics were carried out on the obtained experimental data to evaluate the adequacy of various models (linear, interactive, quadratic, and cubic) fitted to the experimental data. Pareto analysis of variance (ANOVA) was utilized in this present work in order to generate the ANOVA table. The experimental data were evaluated with various descriptive statistical analyses such as p value, F value, degrees of freedom (DF), coefficient of variation (CV), coefficient of determination (R²), adjusted coefficient of determination (Adj-R²), and predicted coefficient of determination (Pre-R²) to reveal the numerical consequence of the constructed quadratic mathematical model.

## 2.7. Optimization and Validation

The effect of fermentation process parameters on the responses, such as the pigment yields and the biomass were optimized by multi-response analysis, using Derringer’s desired function methodology [44]. In this method, the responses were transformed into a dimensionless individual desirability function ($d_i$) that varied from 0 to 1 (lowest to highest desirability). From the geometric means of individual desires, the overall desirability function (D) was obtained using the following equation (Equation (6)).

$$D = (d_{1}^{n_{1}} \times d_{2}^{n_{2}} \times d_{3}^{n_{3}} \times \dots \times d_{k}^{n_{k}})^{1/k} \quad (6)$$

where $d_i$ is the individual desirability ranged from 0 to 1, k is the number of considered responses, and $n_i$ is the weight of each response.

The dimensionless desirability ($d_i$) value of the response using the maximized function was calculated from the equation below (Equation (7)).

$$d_i = \frac{Y_i - Y_{min}}{Y_{max} - Y_{min}}, \quad (7)$$

where, $Y_i$ is the obtained response value, $Y_{min}$ is the obtained response minimum value, and $Y_{max}$ is the obtained response maximum value.

To determine the validity of the developed mathematical model equation, triplicate experiments were performed under the optimal condition, as predicted by the model. The average values of the experimental data were compared with the predicted values of the developed model to find out the accuracy and suitability of this model.

<table><caption>**Table 2.** Box-Behnken design matrix with the experimental responses over process parameters in *Talaromyces albobiverticillius* 30548.</caption><thead><tr><th>Exp. Run</th><th>pH</th><th>Temperature (°C)</th><th>Agitation Speed (rpm)</th><th>Fermentation Time (h)</th><th>Orange Pigment Yield (g/L<sup>1</sup>)</th><th>Red Pigment Yield (g/L<sup>2</sup>)</th><th>Dry Biomass Weight (g/L)</th></tr></thead><tbody><tr><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.30</td><td>0.83</td><td>3.35</td></tr><tr><td>2</td><td>0</td><td>0</td><td>1</td><td>-1</td><td>1.15</td><td>1.44</td><td>5.38</td></tr><tr><td>3</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0.30</td><td>0.82</td><td>1.97</td></tr><tr><td>4</td><td>-1</td><td>0</td><td>0</td><td>1</td><td>0.70</td><td>1.08</td><td>3.05</td></tr><tr><td>5</td><td>0</td><td>-1</td><td>0</td><td>-1</td><td>0.90</td><td>1.23</td><td>3.76</td></tr><tr><td>6</td><td>0</td><td>1</td><td>1</td><td>0</td><td>1.10</td><td>1.37</td><td>6.17</td></tr><tr><td>7</td><td>-1</td><td>0</td><td>0</td><td>-1</td><td>0.65</td><td>1.04</td><td>2.39</td></tr><tr><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1.25</td><td>1.47</td><td>5.47</td></tr><tr><td>9</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.65</td><td>0.95</td><td>3.34</td></tr><tr><td>10</td><td>0</td><td>1</td><td>-1</td><td>0</td><td>1.20</td><td>1.48</td><td>6.05</td></tr><tr><td>11</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0.40</td><td>0.85</td><td>4.83</td></tr><tr><td>12</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1.40</td><td>1.68</td><td>5.10</td></tr><tr><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1.00</td><td>1.31</td><td>3.33</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Exp. Run</th><th>pH</th><th>Temperature (°C)</th><th>Agitation Speed (rpm)</th><th>Fermentation Time (h)</th><th>Orange Pigment Yield (g/L<sup>1</sup>)</th><th>Red Pigment Yield (g/L<sup>2</sup>)</th><th>Dry Biomass Weight (g/L)</th></tr></thead><tbody><tr><td>14</td><td>0</td><td>-1</td><td>-1</td><td>0</td><td>-0.10</td><td>0.43</td><td>0.21</td></tr><tr><td>15</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0.65</td><td>1.03</td><td>3.75</td></tr><tr><td>16</td><td>-1</td><td>-1</td><td>0</td><td>0</td><td>1.20</td><td>1.39</td><td>5.34</td></tr><tr><td>17</td><td>0</td><td>1</td><td>0</td><td>-1</td><td>0.10</td><td>0.57</td><td>2.66</td></tr><tr><td>18</td><td>1</td><td>0</td><td>-1</td><td>0</td><td>1.20</td><td>1.48</td><td>4.28</td></tr><tr><td>19</td><td>0</td><td>-1</td><td>1</td><td>0</td><td>1.20</td><td>1.50</td><td>4.61</td></tr><tr><td>20</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1.55</td><td>1.73</td><td>6.15</td></tr><tr><td>21</td><td>1</td><td>0</td><td>0</td><td>-1</td><td>0.85</td><td>1.18</td><td>5.21</td></tr><tr><td>22</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0.70</td><td>1.12</td><td>2.10</td></tr><tr><td>23</td><td>0</td><td>0</td><td>-1</td><td>-1</td><td>0.75</td><td>1.12</td><td>3.87</td></tr><tr><td>24</td><td>1</td><td>-1</td><td>0</td><td>0</td><td>0.60</td><td>0.98</td><td>2.93</td></tr><tr><td>25</td><td>-1</td><td>1</td><td>0</td><td>0</td><td>1.70</td><td>1.92</td><td>6.67</td></tr><tr><td>26</td><td>1</td><td>0</td><td>0</td><td>1</td><td>1.71</td><td>1.92</td><td>6.7</td></tr><tr><td>27</td><td>-1</td><td>0</td><td>-1</td><td>0</td><td>1.76</td><td>1.93</td><td>6.79</td></tr><tr><td>28</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1.77</td><td>1.93</td><td>7.17</td></tr><tr><td>29</td><td>0</td><td>0</td><td>-1</td><td>1</td><td>1.75</td><td>1.93</td><td>6.65</td></tr></tbody></table>

1: g/L equiv. quinizarin 2: g/L equiv. RYrp.

# 3. Results

## 3.1. Box-Behnken Experimental Design Analysis

The selected statistical methods measured the effects of change in the operating variables and their mutual interactions on a system or a process through experimental design [45]. In this study, using the Design-Expert® Software (v9), the Box-Behnken experimental design with four factors at three levels was employed to investigate and optimize the influence of the process variables on the yields of pigments (orange and red) and fungal biomass for *T. albobiverticillius* 30548. The results are listed in Table 2.

## 3.2. Development of Second-Order Polynomial Models

Model adequacy checking was performed on the experimental data to determine whether the fitted model would give poor or misleading results. Four-degree polynomial models viz., linear, interactive (2 factors interaction, 2FI), quadratic and cubic models were fitted to the experimental data. Two different tests namely the sequential model sum of squares and model summary statistics were carried out in this study, to conclude the adequacy of models among various models to represent the responses. The adequacy of model summary output indicated that, the quadratic model was statistically highly significant due to its higher R², compared to the other models (Table 3). Hence, the quadratic model was selected in this study to investigate the effect of the process variables on the responses, such as pigments and biomass production in the *T. albobiverticillius* 30548 strain.

## 3.3. Determination of Second-Order Polynomial Equations

An empirical relationship expressed by a second-order polynomial equation with interaction terms was fitted between the experimental results obtained on the basis of the Box-Behnken experimental design and the input variables. The observed experimental results in each run were subjected to multiple regression analysis to calculate the regression coefficients of the model. Calculated regression coefficients were substituted in (Equation (4)) to obtain a model for the orange pigment yield OPY (Equation (8)), red pigment yield RPY (Equation (9)), and dry biomass weight DBW (Equation (10)). The final equations obtained in terms of the coded factors are given below:
$$ OPY = +1.74 + 0.35X_1 - 0.11X_2 + 0.27X_3 - 0.029X_4 - 0.11X_1X_2 - 0.19X_1X_3 + 0.11X_1X_4 + 0.41X_2X_3 + 2.145^{-17} X_2X_4 + 0.10X_3X_4 - 0.42X_1^2 - 0.68X_2^2 - 0.35X_3^2 - 0.38X_4^2 \quad (8) $$

$$ RPY = +1.92 + 0.28X_1 - 0.09X_2 + 0.20X_3 - 0.015X_4 - 0.087X_1X_2 - 0.17X_1X_3 + 0.075X_1X_4 + 0.31X_2X_3 - 0.019X_2X_4 + 0.071X_3X_4 - 0.33X_1^2 - 0.54X_2^2 - 0.31X_3^2 - 0.32X_4^2 \quad (9) $$

$$ DBW = +6.80 + 0.77X_1 - 0.77X_2 + 1.24X_3 - 0.17X_4 - 0.24X_1X_2 - 0.021X_1X_3 - 0.61X_1X_4 + 1.18X_2X_3 + 0.55X_2X_4 + 0.17X_3X_4 - 1.02X_1^2 - 2.31X_2^2 - 1.35X_3^2 - 0.97X_4^2 \quad (10) $$

where $X_1$, $X_2$, $X_3$, and $X_4$ are the coded values of the independent process variables, respectively, pH, temperature, agitation speed, and fermentation time.

<table><caption>**Table 3.** Model summary statistics for responses.</caption><thead><tr><th>Source</th><th>Std. Dev.</th><th>R<sup>2</sup></th><th>Adjusted R<sup>2</sup></th><th>Predicted R<sup>2</sup></th><th>Press</th><th>Remarks</th></tr></thead><tbody><tr><td colspan="7">Model summary statistics for Orange Pigment Yield</td></tr><tr><td>Linear</td><td>0.46</td><td>0.3306</td><td>0.2191</td><td>0.1387</td><td>6.54</td><td></td></tr><tr><td>2FI</td><td>0.48</td><td>0.4574</td><td>0.1560</td><td>0.0518</td><td>7.20</td><td></td></tr><tr><td>Quadratic</td><td>0.089</td><td>0.9853</td><td>0.9706</td><td>0.9172</td><td>0.63</td><td>Suggested</td></tr><tr><td>Cubic</td><td>0.072</td><td>0.9959</td><td>0.9809</td><td>0.4736</td><td>4.00</td><td>Aliased</td></tr><tr><td colspan="7">Model summary statistics for Red Pigment Yield</td></tr><tr><td>Linear</td><td>0.37</td><td>0.3146</td><td>0.2004</td><td>0.1230</td><td>4.29</td><td></td></tr><tr><td>2FI</td><td>0.39</td><td>0.4328</td><td>0.1178</td><td>0.0127</td><td>4.82</td><td></td></tr><tr><td>Quadratic</td><td>0.096</td><td>0.9734</td><td>0.9468</td><td>0.8469</td><td>0.75</td><td>Suggested</td></tr><tr><td>Cubic</td><td>0.058</td><td>0.9958</td><td>0.9806</td><td>0.402</td><td>2.92</td><td>Aliased</td></tr><tr><td colspan="7">Model summary statistics for Dry Biomass Weight</td></tr><tr><td>Linear</td><td>1.49</td><td>0.3828</td><td>0.2799</td><td>0.1902</td><td>69.91</td><td></td></tr><tr><td>2FI</td><td>1.58</td><td>0.4820</td><td>0.1943</td><td>-0.0036</td><td>86.64</td><td></td></tr><tr><td>Quadratic</td><td>0.42</td><td>0.9714</td><td>0.9429</td><td>0.8445</td><td>13.42</td><td>Suggested</td></tr><tr><td>Cubic</td><td>0.18</td><td>0.9977</td><td>0.9894</td><td>0.9790</td><td>1.81</td><td>Aliased</td></tr></tbody></table>

### 3.4. Statistical Analysis

The adequacy and fitness of the models were tested by multiple regression analysis using the least square method. Significance of the developed models could be determined through ANOVA and the results are shown in Table A1 (Appendix B). The results of ANOVA indicated that, the developed models adequately represented the actual relationship between the independent variables and responses, in the chosen range. Analysis of variance followed by Fisher's statistical test (F-test) was applied to evaluate the significance of each variable. The high F- values, 67.05 for the orange pigment yield (OPY) measured at 470 nm, 36.61 for red pigment yield (RPY) measured at 500 nm, and 34.02 for the dry biomass weight (DBW) indicated that most of the variation in the response could be explained by the developed regression equations. The associated *p*-values were used to estimate whether F is large enough to indicate statistical significance, and *p*-values lower than 0.05 indicated that the developed model and the terms were statistically significant. In our study, the *p*-values were lower than 0.0001 for all responses (OPY, RPY, DBW). It exhibited the precision and the accuracy of the developed models.

Determination coefficient (R²), adjusted R² (AdjR²), predicted R² (PreR²), and coefficient of variation (CV%) were calculated to check the adequacy (Adeq. Pre.) and accuracy of the developed models. The R² gave the proportion of total variation in the responses predicted by the models. The values of R² (0.9853 for OPY, 0.9734 for RPY, and 0.9714 for DBW) ensured a satisfactory fit of the quadratic model to the experimental data. For example, R² values indicated that the sample variation of 98.5% for OPY was attributed to the independent variables and only 1.5% of the total variations were not explained by the model. The adjusted determination coefficient (AdjR²) corrected the R² value of the sample size and the number of terms in the model. In this study, the values of AdjR² (0.9706 for OPY at 470 nm, 0.9468 for RPY at 500 nm, and 0.942 for DBW) were also high and very close to the R²
values, and indicated a better prediction of the model. However, when irrelevant variables are added to the model, the adjusted R² decreases. In general, AdjR² will always be less than or equal to R². PreR² is a measure of how good the model predicts a response value. In our case, the PreR² (0.9172 for OPY at 470 nm, 0.8469 for RPY at 500 nm, and 0.8445 for DBW) are in reasonable agreement with the AdjR².

The CV%, indicating the relative dispersion of the experimental points from the predictions of the second-order polynomial (SOP) models, were found to be 9.12 for OPY at 470 nm, 7.41 for RPY at 500 nm, and 9.41 for DBW, respectively. The low values of CV% clearly indicated a very high degree of precision and a good reliability of the experimental values. The high R² value and a small CV% value indicated that the developed model would be able to give a good estimate of response of the system over the ranges studied.

### 3.5. Effect of Process Variables on Pigment Yield

In our study, the growth of the fungus in PDB leads to the production of secondary metabolites, which are visualized through the formation of a colored broth during the fermentation process. Production of pigments started at 48 h, attained a maximum at 140 h, and showed a slight increase, thereafter, until the end of fermentation [39,46]. In pigment production, temperature and pH acts as a whole active mechanism that is probably related to the genetic and metabolic control of defense mechanisms, thus influencing the level of pigment concentration and high productivity of pigments [47].

#### 3.5.1. Combined Effect of pH, Temperature, Agitation Speed, and Time on Orange Pigment Yield (OPY)

The effect of independent variables on orange pigment yield was studied by changing the levels of any two independent variables, while keeping the other two at their constant middle level. The response surface plots (RSM) were used to locate the optimal values for this fermentation process. Therefore, three response surface plots were obtained by considering all possible combinations.

Figure 1a shows the effect of interaction between pH and temperature on the orange pigment production (maximum absorbance at 470 nm). It revealed that both factors at their lower levels had a negative impact on the orange pigment production. Increase in pH and temperature led to a gradual increase in pigment production, of up to pH 7.5 and temperature 24 °C. Increasing the value of both independent variables (pH above 7.8 and temperature above 27 °C) also showed a negative effect on orange pigment production. As shown in Figure 1b, increasing the agitation speed (up to 155 rpm) and pH (up to 7.5) led to high pigment production (1.76 g/L of orange pigments quinizarin equivalent), while agitation speed above 155 rpm and pH above 7.5 showed a negative effect on orange pigment production.

The interaction between fermentation time and pH also play an important role in pigment production (Figure 1c); pH in the range of up to 8.1 (alkaline pH) at higher levels of fermentation time, was found to be significant for orange pigment production. It was noticed that fermentation time (234 h) and pH (<8.1) up to a certain level supported orange pigment production but both factors negatively applied at higher levels and affected the pigment production. Temperature and agitation speed showed significant influence on each other and also on OPY. Both parameters showed a linear and quadratic response for pigment yield (Figure 1d). Maximum yield was observed in the above mid-values for both parameters (temperature of 25 °C and agitation speed of 165 rpm) and above that level, yields of the pigment decreased. The interactive effects of temperature and fermentation time against OPY are depicted in Figure 1e. When temperature increased from 21 to 24 °C and fermentation time was between 24 and 259 h, it resulted in a gradual increase in pigment production, up to a maximum. Further increase in temperature (above 24 °C) and fermentation time (above 259 h) led to a decrease in pigment yield.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Response surface plots showing the influence of the process variables on the orange pigment yield (OPY): (a) Interaction between temperature and pH; (b) Interaction between agitation speed and pH; (c) Interaction between time and pH; (d) Interaction between agitation speed and temperature; (e) Interaction between time and temperature; (f) Interaction between time and agitation speed.</figcaption></figure>

The three-dimensional (3D) response interactive plot of agitation speed and fermentation time illustrated that, OPY had increased with increase in agitation speed and fermentation time, up to 180 rpm and 224 h, respectively (Figure 1f). The yield of orange pigment was decreased with an increase in the agitation speed and fermentation time above 180 rpm and 247 h, respectively.

Combining the effects of all process variables investigated in this work, the initial pH and temperature of the medium played a key role in orange pigment synthesis. It was reported that a more acidic pH is favorable for most fungi for pigment synthesis in a submerged culture [4,48,49]. Similarly, acidic pH around 6.5 combined with temperature of up to 23.9 °C, agitation speed of 154.9 rpm, and fermentation time of 229 h was considered to be the optimal condition to produce maximal OPY of 1.76 g/L quinizarin equivalent.

### 3.5.2. Combined Effects of pH, Temperature, Agitation Speed, and Time on Red Pigment Yield (RPY)

Data obtained from the experiments were used to study the effect of process variables on red pigment production. From the results, it was confirmed that increasing pH and temperature from 4.0–6.5 and 22–25 °C enhanced the intensity of red pigment production in *T. albobiverticillius* 30548 and then decreased it (Figure 2a).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Response surface plots showing the influence of the process variables on the red pigment yield (RPY): (a) Interaction between temperature and pH; (b) Interaction between agitation speed and pH; (c) Interaction between time and pH; (d) Interaction between agitation speed and temperature; (e) Interaction between time and temperature; (f) Interaction between time and agitation speed.</figcaption></figure>

Additionally, increasing the agitation speed from 100–158 rpm and the pH value from 4.0–6.5 enhanced the red pigment production (Figure 2b). The maximum production of red pigments was attained at 120–168 h of fermentation and further increase in incubation time led to a decrease in the production (Figure 2c,e,f). The mycelium of the species grew rapidly at the start of fermentation (12–48 h), matured in the period of 48–96 h, and more extracellular red metabolites were released for up to 168 h, in the fermentation medium. Further increasing incubation time led to the bleaching of the extracellular red color and then a decrease in the pigment yield. To examine the influence of agitation speed on red pigment production, experiments were carried out in various agitation speeds (100–200 rpm). The results showed that, red pigment production increased from 100–158 rpm and above that, it decreased rapidly (Figure 2d,f).

Maximum pigment production was exhibited in the middle level of the process variables and strong interaction was found between the variables, since the shape of the response surface and contour plot was elliptical in nature. By applying the numerical optimization method, the optimal condition was attained as follows—initial pH of 6.5, temperature of 24 °C, agitation speed of 158.4 rpm, and fermentation time of 198.6 h, with a maximal red pigment yield of 1.92 g/L, in terms of RYrp equivalent, respectively.

### 3.5.3. Combined Effects of pH, Temperature, Agitation Speed, and Time on Dry Biomass Weight (DBW)

The influence of the process variables over the dry biomass weight was similarly examined and the outcomes were depicted in Figure 3. The graph (Figure 3a) demonstrated the biomass yield was
linearly increased with increasing levels of pH up to 7.8 and temperature up to 25 °C, and then the
yields decreased gradually.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Response surface plots showing the influence of the process variables on dry biomass yield: (a) Interaction between temperature and pH; (b) Interaction between agitation speed and pH; (c) Interaction between time and pH; (d) Interaction between agitation speed and temperature; (e) Interaction between time and temperature; (f) Interaction between time and agitation speed.</figcaption></figure>

Figure 3b shows the effect of interaction between pH and agitation speed on biomass production, which revealed that both factors at their lower level had no significant effect on the biomass production, but increase in pH (up to pH 7.7) and agitation speed (up to 180 rpm) led to a gradual increase in the biomass production, while increase in the value of both independent variables beyond pH 7.7 and agitation speed of 180 rpm showed a negative effect on biomass production. From Figure 3c, it can be seen that the biomass yield increased with increasing pH and fermentation time, up to its middle level. However, further increasing pH and fermentation time would decrease the biomass yield. The 3D response surface plots based on independent variables such as temperature and agitation speed were developed, while other variables were kept at middle levels (Figure 3d) and exhibited that, the biomass yield was enhanced with increasing temperature (24.5 °C) and agitation speed (185 rpm), and the yield
decreased beyond these levels. Biomass production was also increased with increase in both variables (temperature (up to 25 °C) and fermentation time (up to 245 h)) (Figure 3e). However, there was a sharp convergence of the curve near the boundary, explaining that in the presence of temperature and fermentation time above certain limit (above temperature of 25 °C and fermentation time of 245 h), it would not contribute to increasing biomass production.

The interaction between agitation speed and fermentation time also played an important role in biomass production (Figure 3f). Agitation speed and fermentation time at above middle level was found to be significant for biomass production. It was noticed that agitation speed up to a certain level (190 rpm), as well as fermentation time (245 h) supported the biomass production but both affected the biomass production negatively, at higher levels. The analysis of response surface was performed in order to determine the optimal condition to produce maximal biomass yield.

The optimal condition was an initial pH of 6.6, temperature of 23.9 °C, agitation speed of 168.9 rpm, and fermentation time of 145.1 h, respectively, and the maximum predicted yield of biomass (7.16 g/L) was attained in the optimal condition.

## 3.6. Multi Response Optimization

The second-order polynomial models developed in this study were utilized for each response, in order to obtain specified optimum conditions. An optimal condition for the maximum pigment production (orange and red) and biomass yield in *T. albobiverticillius* 30548 under submerged fermentation was derived by Derringer's desired function methodology. This function searches for a combination of factor levels that simultaneously satisfies the requirements for each response in the design. Individual desirability ($d_i$) evaluates how the process parameters optimize a single response. This numerical optimization evaluates a point that maximizes the desirability function. One main objective of optimization could be to maximize the final pigment yield and biomass production and recalculating all responsible independent factors by using desirability functions, therefore, the goal for pH, temperature, agitation speed, and fermentation time was assigned, as in range.

The goal for the orange pigment yield (OPY), red pigment yield (RPY), and biomass production (DBW) was assigned to get the maximum. A weight factor of 1 was chosen for all individual desirability in this work. The "importance" of a goal could be changed in relation to the other goals. It could range from 1 (least importance) to 5 (most important). The default is for all goals to be equally important in a setting of 3. The optimization procedure was conducted under these settings and boundaries. Under the optimal conditions, the predicted orange pigment yield was 1.76 g/L quinizarin equivalent, red pigment yield was 1.92 g/L RYrp equivalent, and dry biomass yield was 7.16 g/L, with a desirability value of 0.983. The maximized overall desirability ($D = 0.983$) was calculated from the geometric means of the individual desirability functions ($d_i$) of each response.

## 3.7. Validation of the Optimized Condition

The validation was carried out in shake flasks in triplicates under the optimized conditions of the media, predicted by the polynomial model. The mean of experimental values of OPY, RPY, and DBW were compared with the predicted values. The existence of a good correlation between the experimental and predicted values indicated the reliability and validity of the proposed model (Figure 4a-c). Due to the applicability of optimal extraction condition in a practical manner, the attained optimal condition was altered as follows—initial pH of 6.4, temperature of 24 °C, agitation speed of 164 rpm, and fermentation time of 149 h, respectively. The experimental efficiency of the production of pigments and dry biomass under the optimum condition was found to be 1.76 ± 0.58 g/L quinizarin equivalent for OPY, 1.92 ± 0.76 g/L RYrp equivalent for RPY, and 7.18 ± 0.41 g/L for DBW, respectively.
<figure><img src="image_7.png" /></figure>

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

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Validation of the polynomial model for (a) OPY, (b) red pigment yield (RPY), and (c) dry biomass weight (DBW) productions from *Talaromyces albobiverticillius* 30548.</figcaption></figure>

**4. Discussion**

The results obtained in this work for the cultivation of *T. albobiverticillius* 30548 were consistent with Babitha et al. (2007), who reported that maximum production of pigments was noticed in *M. purpureus* at pH 4.5 to 7.5 [50]. Studies by many researchers have revealed that among the numerous environmental factors, medium pH and the source of nitrogen determine pigment production in a submerged culture. The optimal pH for *T. albobiverticillius* 30548 was found to be 6.4 in this study, which is close to the optimum pH for pigment production in *Talaromyces verruculosus* mentioned by Chadni et al. (2017) [51]. Indeed, acidic pH might enhance the hydrolysis of the substrates and subsequently favors the generation of the metabolite production [52]. From Yongsmith et al. (1993) it was demonstrated that the compounds excreted in the culture broth could react with ammonium
ion or free amino groups to transform to red amine derivatives towards neutral pH and beyond [53]. This might happen in *T. albobiverticillius* 30548 as there is an evolution of the color hue in the culture broth (from yellow to red), along the entire fermentation period [39].

Temperature, another key parameter, controls the growth and metabolic activity of fungi. To find the optimal temperature for mycelial growth and pigment production, *T. albobiverticillius* 30548 was cultivated under various temperatures (21–27 °C). Fermentation at 24 °C is regarded as favorable for fungal growth in combination with an acidic pH and a medium agitation speed, which aid good production of pigments. Mendez et al. (2011) found a similar result in obtaining the highest level of pigment production at 24 °C in *P. purpurogenum* GH2, combined with an acidic pH [28]. However, in other studies with *Monascus spp.*, *Penicillium spp.*, temperature at 30 °C is considered to be the best temperature for red pigment production. However, the microbial growth was hindered at 30 °C but a maximum growth was attained by increasing the temperature to 34 °C [26,54,55]. For the marine-derived *T. albobiverticillius* 30548, it is interesting to note that the optimal temperature for both mycelial growth and pigment production was found to be 24 °C.

Fermentation of any substrate is influenced by the time of incubation, so the pigment synthesis was monitored at predefined intervals. Pigment production, as measured by absorbance, was slightly delayed when compared to biomass growth, and was higher at 149 h (day 6) of fermentation time, until the stationary phase [39,46]. After 149 h, a decrease in pigment production was observed, which might be due to the decomposition of pigments (degradation of the chromophore pigment group) or changes in the pigment structure. A similar outcome was observed in the pigment production by *Monascus purpureus* [56].

Basically, in liquid fermentation, agitation speed increases the amount of dissolved oxygen and makes the oxygen more accessible to cells, leading to greater growth [57]. Similarly, in this present work, increasing agitation speed from 100–158 rpm enhanced the supply of oxygen in the growth phase and thus enhanced the biomass and pigment yield. However, above 158 rpm, agitation speed might damage the mycelial growth by high shear force and lead to a decrease in the pigment yields. Mohamed et al. (2012) reported similar results, as the highest pigment production was obtained at agitation speed of 150 rpm using *Monascus purpureus* [58]. Additionally, the highest level of yellow pigments (1.38 g/L) was produced by *Penicillium aculeatum* ATCC 10409 at an agitation speed of 100 and 150 rpm [47].

**5. Conclusions**

Box-Behnken experimental design (BBED) coupled with response surface modeling (RSM) was successfully applied in the present work. The methodology allowed to locate the optimal process parameters in the chosen range (initial pH (4-9), temperature (21-27 °C), agitation speed (100-200 rpm), and fermentation time (0-336 h)), to maximize pigment and biomass production in submerged fermentation of *T. albobiverticillius* 30548. The outcomes of this study indicated that, all process variables had a significant influence on the responses. The actual values obtained through the experimental studies were used to construct second-order quadratic model for the responses to predict the observed data. The optimal condition was attained through Derringer's desired function methodology (initial pH of 6.4, temperature of 24 °C, agitation speed of 164 rpm, and fermentation time of 149 h). It was materially validated by the experiment (1.76 ± 0.58 g/L quinizarin equivalent for OPY, 1.92 ± 0.76 g/L RYrp equivalent for RPY, and 7.18 ± 0.41 g/L for DBW).

The results of this study revealed that *T. albobiverticillius* 30548 could be an avenue for industrial application of pigment production. Considering the potential toxicity of compounds produced by *T. albobiverticillius* 30548 in SmF, complementary studies need to be conducted, based on the characterization of the 12 major compounds synthesized [59], and on the elucidation of the biosynthetic pathway, before using them for applications. In addition, the mathematical models constructed in this study and the methodology applied could be useful for the selection of appropriate parameters for the production of pigments at a large scale.
This study allowed for the consideration of other factors of pigment production, besides the
physical process factors. Generally, carbon and nitrogen sources in the culture medium play a major
role in the growth and production of metabolites. Further work on optimizing various nutrient types
and their concentration in the culture media is in progress and will be developed in a further study.

**Author Contributions:** Conceptualization, methodology, writing—original draft preparation, M.V. and M.F.; investigation and visualization, M.V.; formal analysis, M.V. and A.S.-C.-S.; writing—review and editing, M.V., M.F., L.D., and A.S.-C.-S.; supervision, L.D. and M.F.; funding acquisition and project administration, M.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** We would like to thank the Conseil Régional de La Réunion, Reunion Island, France, for the financial support of research activities dedicated to microbial pigments (DIRED 20140704 - COLORMAR Program and allocation régionale de recherche n° 2016020577/181733).

**Acknowledgments:** Authors are thankful to J. Prakash Maran for his valuable advices about the design of experiments and Cathie Milhau for technical support.

**Conflicts of Interest:** The authors declare no conflict 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.

Appendix A

<figure><img src="image_10.png" /><figcaption>Figure A1. Standard curve for quinizarin at 470 nm.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure A2. Standard curve for red yeast rice pigment (RYrp) at 500 nm.</figcaption></figure>
## Appendix B

<table><caption>**Table A1.** Regression coefficient (RC) of the models and their statistical parameters on the responses (Xi: coded values of independent variables (i = 1: pH; i = 2: temperature, i = 3: agitation speed; i = 4: fermentation time) and their interactions.</caption><thead><tr><th rowspan="2">Source</th><th rowspan="2">DF</th><th colspan="2">Orange Pigment Yield</th><th colspan="2">Red Pigment Yield</th><th colspan="2">Dry Biomass Weight</th></tr><tr><th>RC</th><th><i>p</i> Value *</th><th>RC</th><th><i>p</i> Value *</th><th>RC</th><th><i>p</i> Value *</th></tr></thead><tbody><tr><td>Model</td><td>14</td><td>1.74</td><td>&lt;0.0001</td><td>1.92</td><td>&lt;0.0001</td><td>6.80</td><td>&lt;0.0001</td></tr><tr><td>X<sub>1</sub></td><td>1</td><td>0.35</td><td>&lt;0.0001</td><td>0.28</td><td>&lt;0.0001</td><td>0.77</td><td>&lt;0.0001</td></tr><tr><td>X<sub>2</sub></td><td>1</td><td>-0.11</td><td>0.0009</td><td>-0.090</td><td>0.0058</td><td>-0.77</td><td>&lt;0.0001</td></tr><tr><td>X<sub>3</sub></td><td>1</td><td>0.27</td><td>&lt;0.0001</td><td>0.20</td><td>&lt;0.0001</td><td>1.24</td><td>&lt;0.0001</td></tr><tr><td>X<sub>4</sub></td><td>1</td><td>-0.029</td><td>0.2767</td><td>-0.015</td><td>0.5901</td><td>-0.17</td><td>0.1837</td></tr><tr><td>X<sub>12</sub></td><td>1</td><td>-0.11</td><td>0.0245</td><td>-0.087</td><td>0.0927</td><td>-0.24</td><td>0.2788</td></tr><tr><td>X<sub>13</sub></td><td>1</td><td>-0.19</td><td>0.0009</td><td>-0.17</td><td>0.0031</td><td>-0.02</td><td>0.9224</td></tr><tr><td>X<sub>14</sub></td><td>1</td><td>0.11</td><td>0.0245</td><td>0.075</td><td>0.1394</td><td>-0.61</td><td>0.0117</td></tr><tr><td>X<sub>23</sub></td><td>1</td><td>0.41</td><td>&lt;0.0001</td><td>0.31</td><td>&lt;0.0001</td><td>1.18</td><td>&lt;0.0001</td></tr><tr><td>X<sub>24</sub></td><td>1</td><td>0.000</td><td>1.000</td><td>-0.019</td><td>0.7011</td><td>0.55</td><td>0.0210</td></tr><tr><td>X<sub>34</sub></td><td>1</td><td>0.10</td><td>0.0418</td><td>0.071</td><td>0.1615</td><td>0.17</td><td>0.4359</td></tr><tr><td>X<sub>1</sub><sup>2</sup></td><td>1</td><td>-0.42</td><td>&lt;0.0001</td><td>-0.33</td><td>&lt;0.0001</td><td>-1.02</td><td>&lt;0.0001</td></tr><tr><td>X<sub>2</sub><sup>2</sup></td><td>1</td><td>-0.68</td><td>&lt;0.0001</td><td>-0.54</td><td>&lt;0.0001</td><td>-2.31</td><td>&lt;0.0001</td></tr><tr><td>X<sub>3</sub><sup>2</sup></td><td>1</td><td>-0.35</td><td>&lt;0.0001</td><td>-0.31</td><td>&lt;0.0001</td><td>-1.35</td><td>&lt;0.0001</td></tr><tr><td>X<sub>4</sub><sup>2</sup></td><td>1</td><td>-0.38</td><td>&lt;0.0001</td><td>-0.32</td><td>&lt;0.0001</td><td>-0.97</td><td>&lt;0.0001</td></tr><tr><td>R<sup>2</sup></td><td></td><td>0.9853</td><td></td><td>0.9734</td><td></td><td>0.9714</td><td></td></tr><tr><td>AdjR<sup>2</sup></td><td></td><td>0.9706</td><td></td><td>0.9468</td><td></td><td>0.9429</td><td></td></tr><tr><td>PreR<sup>2</sup></td><td></td><td>0.9172</td><td></td><td>0.8469</td><td></td><td>0.8445</td><td></td></tr><tr><td>CV%</td><td></td><td>9.12</td><td></td><td>7.41</td><td></td><td>9.41</td><td></td></tr><tr><td>Adeq. Pre.</td><td></td><td>28.30</td><td></td><td>21.04</td><td></td><td>22.71</td><td></td></tr></tbody></table>

R<sup>2</sup>: Determination coefficient; AdjR<sup>2</sup>: adjusted R<sup>2</sup>; PreR<sup>2</sup>: predicted R<sup>2</sup>; CV%: coefficient of variation; and Adeq. Pre.: predicted adequacy. *p* < 0.01 highly significant; 0.01 < *p* < 0.05 significant; and *p* > 0.05 not significant.

## References

1. Weiss, B. Synthetic Food Colors and Neurobehavioral Hazards: The View from Environmental Health Research. *Environ. Health Perspect.* **2012**, *120*, 1–5. [CrossRef]

2. Dufossé, L. Microbial production of food grade pigments. *Food Technol. Biotechnol.* **2006**, *44*, 313–323.

3. Dufosse, L. Anthraquinones, the Dr Jekyll and Mr Hyde of the food pigment family. *Food Res. Int.* **2014**, *65*, 132–136. [CrossRef]

4. Kumar, A.; Vishwakarma, H.S.; Singh, J.; Dwivedi, S.; Kumar, M. Microbial pigments: Production and their applications in various industries. *IJPCBS* **2015**, *5*, 203–212.

5. Kamarudin, K.R. Microbial production of food grade pigments-screening and metabolic pathway analysis. *Food Technol. Biotechnol.* **2006**, *44*, 13–23.

6. Caro, Y.; Anamale, L.; Fouillaud, M.; Laurent, P.; Petit, T.; Dufosse, L. Natural hydroxyanthraquinoid pigments as potent food grade colorants: An overview. *Nat. Prod. Bioprospect.* **2012**, *2*, 174–193. [CrossRef]

7. Duran, N.; Teixeira, M.F.; De Conti, R.; Esposito, E. Ecological-friendly pigments from fungi. *Crit. Rev. Food Sci. Nutr.* **2002**, *42*, 53–66. [CrossRef]

8. Ventura, S.P.; Santos-Ebinuma, V.C.; Pereira, J.F.; Teixeira, M.F.; Pessoa, A.; Coutinho, J.A. Isolation of natural red colorants from fermented broth using ionic liquid-based aqueous two-phase systems. *J. Ind. Microbiol. Biotechnol.* **2013**, *40*, 507–516. [CrossRef]

9. Pastre, R.; Marinho, A.M.; Rodrigues-Filho, E.; Souza, A.Q.; Pereira, J.O. Diversity of polyketides produced by Penicillium species isolated from Melia azedarach and *Murraya paniculata*. *Química Nova* **2007**, *30*, 1867–1871. [CrossRef]

10. Bugni, T.S.; Ireland, C.M. Marine-derived fungi: A chemically and biologically diverse group of microorganisms. *Nat. Prod. Rep.* **2004**, *21*, 143–163. [CrossRef] [PubMed]
11. Osmanova, N.; Schultze, W.; Ayoub, N. Azaphilones: A class of fungal metabolites with diverse biological activities. *Phytochem. Rev.* **2010**, *9*, 315–342. [CrossRef]

12. Kongruang, S. Growth kinetics of biopigment production by Thai isolated *Monascus purpureus* in a stirred tank bioreactor. *J. Ind. Microbiol. Biotechnol.* **2011**, *38*, 93–99. [CrossRef] [PubMed]

13. Teixeira, M.F.; Martins, M.S.; Da Silva, J.; Kirsch, L.S.; Fernandes, O.C.; Carneiro, A.L.; De Conti, R.; Durán, N. Amazonian biodiversity: Pigments from *Aspergillus* and *Penicillium*-characterizations, antibacterial activities and their toxicities. *Curr. Trends Biotechnol. Pharm.* **2012**, *6*, 300–311.

14. Torres, F.A.E.; Zaccarim, B.R.; De Lencastre Novaes, L.C.; Jozala, A.F.; Santos, C.A.d.; Teixeira, M.F.S.; Santos-Ebinuma, V.C. Natural colorants from filamentous fungi. *Appl. Microbiol. Biotechnol.* **2016**, *100*, 2511–2521. [CrossRef]

15. Bhatnagar, I.; Kim, S.-K. Immense Essence of Excellence: Marine Microbial Bioactive Compounds. *Mar. Drugs* **2010**, *8*, 2673–2701. [CrossRef]

16. Kjer, J.; Debbab, A.; Aly, A.H.; Proksch, P. Methods for isolation of marine-derived endophytic fungi and their bioactive secondary products. *Nat. Protoc.* **2010**, *5*, 479–490. [CrossRef]

17. Imhoff, J.F. Natural products from marine fungi—Still an underrepresented resource. *Mar. drugs* **2016**, *14*, 19. [CrossRef]

18. Yilmaz, N.; Houbraken, J.; Hoekstra, E.; Frisvad, J.C.; Visagie, C.; Samson, R. Delimitation and characterisation of *Talaromyces purpurogenus* and related species. *Pers. Mol. Phylogeny Evol. Fungi* **2012**, *29*, 39–54. [CrossRef]

19. Frisvad, J.C.; Yilmaz, N.; Thrane, U.; Rasmussen, K.B.; Houbraken, J.; Samson, R.A. Talaromyces atroroseus, a new species efficiently producing industrially relevant red pigments. *PLoS ONE* **2013**, *8*, e84102. [CrossRef]

20. Mapari, S.A.; Meyer, A.S.; Thrane, U.; Frisvad, J.C. Identification of potentially safe promising fungal cell factories for the production of polyketide natural food colorants using chemotaxonomic rationale. *Microbial Cell Factories* **2009**, *8*, 1. [CrossRef]

21. Nicoletti, R.; Trincone, A. Bioactive compounds produced by strains of *Penicillium* and *Talaromyces* of marine origin. *Mar. Drugs* **2016**, *14*, 37. [CrossRef]

22. Caro, Y.; Venkatachalam, M.; Lebeau, J.; Fouillaud, M.; Dufossé, L. Pigments and colorants from filamentous fungi. In *Fungal Metabolites*, 1st ed.; Merillon, J.-M., Ramawat, K.G., Eds.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 1–70.

23. Musoni, M.; Destain, J.; Thonart, P.; Bahama, J.-B.; Delvigne, F. Bioreactor design and implementation strategies for the cultivation of filamentous fungi and the production of fungal metabolites: From traditional methods to engineered systems/Conception de bioréacteurs et mise en oeuvre de stratégies pour la culture de champignons filamenteux et la production de métabolites d'origine fongique: Des méthodes traditionnelles aux technologies actuelles. *Biotechnol. Agron. Société Environ.* **2015**, *19*, 430.

24. General, T.; Kim, H.-J.; Prasad, B.; Ngo, H.T.A.; Vadakedath, N.; Cho, M.-G. Fungal utilization of a known and safe macroalga for pigment production using solid-state fermentation. *J. Appl. Phycol.* **2014**, *26*, 1547–1555. [CrossRef]

25. Hajjaj, H.; Blanc, P.; Goma, G.; Francois, J. Sampling techniques and comparative extraction procedures for quantitative determination of intra-and extracellular metabolites in filamentous fungi. *FEMS Microbiol. Lett.* **1998**, *164*, 195–200. [CrossRef]

26. Santos-Ebinuma, V.C.; Roberto, I.C.; Teixeira, M.F.S.; Pessoa, A., Jr. Improvement of submerged culture conditions to produce colorants by *Penicillium purpurogenum*. *Braz. J. Microbiol.* **2014**, *45*, 731–742. [CrossRef]

27. Lee, B.-K.; Park, N.-H.; Piao, H.Y.; Chung, W.-J. Production of red pigments by *Monascus purpureus* in submerged culture. *Biotechnol. Bioprocess Eng.* **2001**, *6*, 341–346. [CrossRef]

28. Méndez, A.; Pérez, C.; Montañéz, J.C.; Martínez, G.; Aguilar, C.N. Red pigment production by *Penicillium purpurogenum* GH2 is influenced by pH and temperature. *J. Zhejiang Univ. Sci. B* **2011**, *12*, 961–968. [CrossRef]

29. Orozco, S.F.B.; Kilikian, B.V. Effect of pH on citrinin and red pigments production by *Monascus purpureus* CCT3802. *World J. Microbiol. Biotechnol.* **2008**, *24*, 263–268. [CrossRef]

30. Liu, L.; Zhao, J.; Huang, Y.; Xin, Q.; Wang, Z. Diversifying of chemical structure of native *Monascus* pigments. *Front. Microbiol.* **2018**, *9*, 3143. [CrossRef]

31. Yilmaz, N.; Visagie, C.M.; Houbraken, J.; Frisvad, J.C.; Samson, R.A. Polyphasic taxonomy of the genus *Talaromyces*. *Studies Mycol.* **2014**, *78*, 175–341. [CrossRef]
32. Zhang, G.; He, L.; Hu, M. Optimized ultrasonic-assisted extraction of flavonoids from *Prunella vulgaris* L. and evaluation of antioxidant activities in vitro. *Innov. Food Sci. Emerg. Technol.* **2011**, *12*, 18–25. [CrossRef]

33. Yang, Z.; Zhai, W. Optimization of microwave-assisted extraction of anthocyanins from purple corn (Zea mays L.) cob and identification with HPLC-MS. *Innov. Food Sci. Emerg. Technol.* **2010**, *11*, 470–476. [CrossRef]

34. Parajó, J.; Santos, V.; Domínguez, H.; Vázquez, M. NH4OH-Based pretreatment for improving the nutritional quality of single-cell protein (SCP). *Appl. Biochem. Biotechnol.* **1995**, *55*, 133–149. [CrossRef]

35. Huang, W.; Xue, A.; Niu, H.; Jia, Z.; Wang, J. Optimised ultrasonic-assisted extraction of flavonoids from *Folium eucommiae* and evaluation of antioxidant activity in multi-test systems in vitro. *Food Chem.* **2009**, *114*, 1147–1154. [CrossRef]

36. Hayouni, E.A.; Abedrabba, M.; Bouix, M.; Hamdi, M. The effects of solvents and extraction method on the phenolic contents and biological activities in vitro of Tunisian *Quercus coccifera* L. and *Juniperus phoenicea* L. fruit extracts. *Food Chem.* **2007**, *105*, 1126–1134. [CrossRef]

37. Ge, Y.; Ni, Y.; Yan, H.; Chen, Y.; Cai, T. Optimization of the supercritical fluid extraction of natural vitamin E from wheat germ using response surface methodology. *J. Food Sci.* **2002**, *67*, 239–243. [CrossRef]

38. Sehrawat, R.; Panesar, P.S.; Swer, T.L.; Kumar, A. Response surface methodology (RSM) mediated interaction of media concentration and process parameters for the pigment production by *Monascus purpureus* MTCC 369 under solid state fermentation. *Pigment. Resin Technol.* **2017**, *46*, 14–20. [CrossRef]

39. Venkatachalam, M.; Magalon, H.; Dufossé, L.; Fouillaud, M. Production of pigments from the tropical marine-derived fungi *Talaromyces albobiverticillius*: New resources for natural red-colored metabolites. *J. Food Compos. Anal.* **2018**, *70*, 35–48. [CrossRef]

40. Bas, D.; Boyaci, I.H. Modeling and optimization I: Usability of response surface methodology. *J. Food Eng.* **2007**, *78*, 836–845. [CrossRef]

41. Bezerra, M.A.; Santelli, R.E.; Oliveira, E.P.; Villar, L.S.; Escaleira, L.A. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. *Talanta* **2008**, *76*, 965–977. [CrossRef]

42. Ferreira, S.C.; Bruns, R.; Ferreira, H.; Matos, G.; David, J.; Brandao, G.; Da Silva, E.P.; Portugal, L.; Dos Reis, P.; Souza, A. Box-Behnken design: An alternative for the optimization of analytical methods. *Anal. Chim. Acta* **2007**, *597*, 179–186. [CrossRef]

43. Karichappan, T.; Venkatachalam, S.; Jeganathan, P.M. Analysis of efficiency of *Bacillus subtilis* to treat bagasse based paper and pulp industry wastewater-A novel approach. *J. Korean Chem. Soc.* **2014**, *58*, 198–204. [CrossRef]

44. Derringer, G.; Suich, R. Simultaneous Optimization of Several Response Variables. *J. Qual. Technol.* **1980**, *12*, 214–219. [CrossRef]

45. Ekren, O.; Ekren, B.Y. Size optimization of a PV/wind hybrid energy conversion system with battery storage using response surface methodology. *Appl. Energy* **2008**, *85*, 1086–1101. [CrossRef]

46. Venkatachalam, M.; Gérard, L.; Milhau, C.; Vinale, F.; Dufossé, L.; Fouillaud, M. Salinity and Temperature Influence Growth and Pigment Production in the Marine-Derived Fungal Strain *Talaromyces albobiverticillius* 30548. *Microorganisms* **2019**, *7*, 10. [CrossRef] [PubMed]

47. Afshari, M.; Shahidi, F.; Mortazavi, S.A.; Tabatabai, F.; Es' haghi, Z. Investigating the influence of pH, temperature and agitation speed on yellow pigment production by *Penicillium aculeatum* ATCC 10409. *Nat. Prod. Res.* **2015**, *29*, 1300–1306. [CrossRef] [PubMed]

48. Cho, Y.; Park, J.; Hwang, H.; Kim, S.; Choi, J.; Yun, J. Production of red pigment by submerged culture of *Paecilomyces sinclairii*. *Lett. Appl. Microbiol.* **2002**, *35*, 195–202. [CrossRef]

49. Bae, J.-T.; Sinha, J.; Park, J.-P.; Song, C.-H.; Yun, J.-W. Optimization of submerged culture conditions for exo-biopolymer production by *Paecilomyces japonica*. *J. Microbiol. Biotechnol.* **2000**, *10*, 482–487.

50. Babitha, S. Microbial Pigments. In *Biotechnology for Agro-Industrial Residues Utilisation: Utilisation of Agro-Residues*; Singh nee' Nigam, P., Pandey, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 147–162.

51. Chadni, Z.; Rahaman, M.H.; Jerin, I.; Hoque, K.; Reza, M.A. Extraction and optimisation of red pigment production as secondary metabolites from *Talaromyces verruculosus* and its potential use in textile industries. *Mycology* **2017**, *8*, 48–57. [CrossRef]

52. Wong, H.; Lin, Y.; Koehler, P. Regulation of growth and pigmentation of *Monascus purpureus* by carbon and nitrogen concentrations. *Mycologia* **1981**, *649–654*. [CrossRef]
53. Yongsmith, B.; Tabloka, W.; Yongmanitchai, W.; Bavavoda, R. Culture conditions for yellow pigment formation by *Monascus* sp. KB 10 grown on cassava medium. *World J. Microbiol. Biotechnol.* **1993**, *9*, 85–90.
[CrossRef]

54. Babitha, S.; Soccol, C.R.; Pandey, A. Solid-state fermentation for the production of *Monascus* pigments from jackfruit seed. *Bioresour. Technol.* **2007**, *98*, 1554–1560. [CrossRef]

55. Gunasekaran, S.; Poorniammal, R. Optimization of fermentation conditions for red pigment production from *Penicillium* sp. under submerged cultivation. *Afr. J. Biotechnol.* **2008**, *7*, 1894–1898. [CrossRef]

56. Chen, M.-H.; Johns, M.R. Effect of pH and nitrogen source on pigment production by *Monascus purpureus*. *Appl. Microbiol. Biotechnol.* **1993**, *40*, 132–138. [CrossRef]

57. Jafari, A.; Sarrafzadeh, M.; Alemzadeh, I.; Vosoughi, M. Effect of stirrer speed and aeration rate on the production of glucose oxidase by *Aspergillus niger*. *J. Biol. Sci.* **2007**, *7*, 270–275.

58. Mohamed, M.S.; Mohamad, R.; Manan, M.A.; Ariff, A.B. Enhancement of red pigment production by *Monascus purpureus* FTC 5391 through retrofitting of helical ribbon impeller in stirred-tank fermenter. *Food Bioprocess. Technol.* **2012**, *5*, 80–91. [CrossRef]

59. Venkatachalam, M.; Zelena, M.; Cacciola, F.; Ceslova, L.; Girard-Valenciennes, E.; Clerc, P.; Dugo, P.; Mondello, L.; Fouillaud, M.; Rotondo, A. Partial characterization of the pigments produced by the marine-derived fungus *Talaromyces albobiverticillius* 30548. Towards a new fungal red colorant for the food industry. *J. Food Compos. Anal.* **2018**, *67*, 38–47. [CrossRef]

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

© 2020 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 (http://creativecommons.org/licenses/by/4.0/).