Next Article in Journal
Classification of Environmental Strains from Order to Genus Levels Using Lipid and Protein MALDI-ToF Fingerprintings and Chemotaxonomic Network Analysis
Next Article in Special Issue
Biofertilizer: The Future of Food Security and Food Safety
Previous Article in Journal
Influence of Different Vegetation Types on Soil Physicochemical Parameters and Fungal Communities
Previous Article in Special Issue
Effects of Rhizobia Isolated from Coffee Fields in the High Jungle Peruvian Region, Tested on Phaseolus vulgaris L. var. Red Kidney
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of Bacillus amyloliquefaciens BLB369 Culture Medium by Response Surface Methodology for Low Cost Production of Antifungal Activity

1
Laboratory of Biopesticides, Centre of Biotechnology of Sfax, University of Sfax, P.O. Box 1177, Sfax 3018, Tunisia
2
Department of Health and Medical Sciences, Khawarizmi International College, Abu Dhabi P.O. Box 25669, United Arab Emirates
3
Laboratory of Probability and Statistics, Faculty of Sciences of Sfax, University of Sfax, P.O. Box 1171, Sfax 3000, Tunisia
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(4), 830; https://doi.org/10.3390/microorganisms10040830
Submission received: 5 March 2022 / Revised: 31 March 2022 / Accepted: 6 April 2022 / Published: 16 April 2022
(This article belongs to the Special Issue Biofertilizer Microorganism in Agriculture)

Abstract

:
Bacillus amyloliquefaciens BLB369 is an important plant growth-promoting bacterium, which produces antifungal compounds. A statistics-based experimental design was used to optimize a liquid culture medium using inexpensive substrates for increasing its antifungal activity. A Plackett–Burman design was first applied to elucidate medium components having significant effects on antifungal production. Then the steepest ascent method was employed to approach the experimental design space, followed by an application of central composite design. Three factors were retained (candy waste, peptone, and sodium chloride), and polynomial and original trigonometric models fitted the antifungal activity. The trigonometric model ensured a better fit. The contour and surface plots showed concentric increasing levels pointing out an optimized activity. Hence, the polynomial and trigonometric models showed a maximal antifungal activity of 251.9 (AU/mL) and 255.5 (AU/mL) for (19.17, 19.88, 3.75) (g/L) and (19.61, 20, 3.7) (g/L) of candy waste, peptone, and NaCl, respectively. This study provides a potential strategy for improving the fermentation of B. amyloliquefaciens BLB369 in low-cost media for large-scale industrial production.

1. Introduction

Phytopathogenic fungi cause several plant diseases responsible for various crop losses and agricultural products deterioration. Fusarium graminearum produces various pathogenicity and virulence factors allowing it to enter into the plant and advance within the interior of the infected tissue. It causes Fusarium head blight (FHB), one of the most economically disastrous diseases of wheat, barley, rice, and other grain crops worldwide [1,2]. It produces deoxynivalenol (DON) and zearalenone, two harmful toxins to humans and animals. Consequently, chemical fungicides were used. Triazoles (e.g., tebuconazole, metconazole, and prothioconazole [3,4]) and benzimidazole obstruct sterol biosynthesis and are the most effective to suppress FHB symptoms and decrease mycotoxin concentration [5,6,7]. Prochloraz, an imidazole derivate acting similarly to triazole fungicides [8], is widely applied to control fungal growth in cereals in several European countries [9]. The widespread use of chemical pesticides (more than 97% of control measures) overexposed nature to their ecotoxicity and led to a loss of efficiency due to adaptation of the targeted phytopathogens and disadvantages on non-target populations sharing the ecosystem [10]. High doses of triazole fungicides strongly affect the structure of the microbial communities in soil and usually decrease the soil microbial population and the activities of enzymes found in soil [11]. The use of non-pathogenic microorganisms as biopesticides is an emerging technology, is ecologically compatible, and is considered a promising alternative to synthetic pesticides [12,13,14]. Numerous Bacillus species offer many advantages for agricultural biotechnology such as Bacillus amyloliquefaciens, B. subtilis, B. licheniformis, and B. pumilus. They synthesize several secondary metabolites, essentially, the cyclic lipopeptides surfactin, iturin, and fengycin with antifungal activities [15,16,17,18]. They also produce primary metabolites with antifungal activities like the hydrolytic enzymes (chitinase, glucanase, and protease enzymes) acting on fungal cell walls [19,20,21]. The growth of cells and metabolite concentrations is influenced strongly by medium composition such as the carbon source, nitrogen source, and inorganic salts. Efforts must therefore be redirected to improve production efficiency and recovery bioprocesses to optimize yields [22,23].
The statistical methods using screening and response surface methodology or artificial neural networks offer several advantages over conventional methods to optimize numerous multi-factorial processes or formulations [24,25,26,27,28]. They shortlist significant nutrients for culture media, help understand the interactions among the nutrients at various concentrations, and reduce the total number of experiments leading to saving time and resources [29,30,31]. They allow the production increase of antagonist compounds, spores, and enzymes by Bacillus spp. strains [32,33,34,35]. Many researchers have used central composite design (CCD) to identify optimal reaction conditions [36,37]. Therefore, the present study was undertaken to optimize a medium for economical production of B. amyloliquefaciens BLB369 antifungal activity. We applied a Plackett–Burman design (PBD) to screen the significant factors, the steepest ascent method (SAM) to approach the experimental design space, and the (CCD) to optimize the concentrations of selected variables using two different regression models.

2. Materials and Methods

2.1. Microorganisms and Cultivation

The B. amyloliquefaciens BLB369 strain was previously characterized by Zalila-Kolsi et al. [38]. In fact, it was isolated in our laboratory from Tunisian rhizosphere soil sample and then identified by using API50CH and API20E strips and partially sequencing 16S rDNA and gyrA genes. It could produce the extracellular cyclic lipopeptides iturin and surfactin harboring antifungal activities. Its antagonist effect against F. graminearum for protection of durum wheat was demonstrated in vivo [38]. The phytopathogenic fungus Fusarium graminearum was kindly provided by the Agricultural Culture Collection of Biopesticides laboratory, Centre of Biotechnology of Sfax. It was maintained on potato dextrose agar (PDA) and used as the target pathogen for testing antifungal activity. For spore suspension preparation, the fungus was incubated on a PDA plate for 5 days at 28 °C; then, 0.9% NaCl solution was added, and the top of the mold was scraped with a sterile loop to release spores. The collected fungus suspension was filtered with sterile cotton to remove mycelial fragments. The obtained spore suspension was enumerated on Malassez cell, adjusted to a concentration of approximately 105 spores/mL, and stored at 4 °C.
The B. amyloliquefaciens BLB369 was grown in a 250 mL flask containing 50 mL of MOLB medium [39] at 30 °C for 14 h, on a rotary shaker set 200 rpm in order to prepare the culture inoculums. For antifungal production, the inoculums were served to inoculate 50 mL of culture medium in a 250 mL flask with an initial OD600 of 0.15, and the culture was incubated at 30 °C for 48 h, shaking at 200 rpm.

2.2. Determination of the Antifungal Activity

The antagonistic activity of B. amyloliquefaciens BLB369 strain against F. graminearum was assessed by the well diffusion method. The BLB369 culture broth of 48 h was centrifuged in 2 mL Eppendorf tube at 10,000 rpm for 15 min, and the culture medium was recovered. Several dilutions of the bacterial culture medium were prepared in the MOLB medium with various dilution factors (DF) (values: 1, 2, 3, 4, 5, …, 30) and used for antifungal activity determination. The dilution factor may be expressed as the ratio of the volume of the final diluted bacterial culture medium to the initial volume removed from the original bacterial culture medium. Then, aliquots (v = 100 µL) were filled in wells of 5 mm diameter made in PDA + chloramphenicol (30 µg/mL) pH 7.0 previously plated with 100 µL of collected F. graminearum suspension (105 spores/mL). After incubation at 28 °C for 3–5 days, the inhibition zones were observed. Antifungal activity (AU/mL) = (DFh * 1000/v), with AU: Arbitrary Unit, DFh: higher dilution factor of the bacterial culture medium able to inhibit fungal growth. v: volume (µL) of the diluted bacterial culture medium used for well test. Each experience was repeated three times.

2.3. Identification of Significant Factors Using Plackett–Burman Design

The PBD [40] was used to select the most significant factors of medium components for biofungicides production where the interactions between the factors were considered negligible. An orthogonal matrix was generated using seven factors, and each factor was represented by a high level (+1) and a low level (−1). Carbon sources, nitrogen sources, and salts could influence the production of antifungal activity. The candy waste was an agro-industrial by-product corresponding to an aqueous solution rich in carbohydrates. Its reducing sugar content was estimated to be about 56.5 g/L using the colorimetric method with 3.5-dinitrosalicylic acid reagent (DNS) [41]. The fish extract was obtained from tuna canning industry waste after being treated with NaOH solution at pH 10-12 for 2 h. Its nitrogen content (non-protein nitrogen and nitrogen in proteins) was estimated to be about 40 g/L using the Kjeldahl method. Yeast extract is essentially composed of amino acids, peptides, carbohydrates, and soluble vitamins. Peptone from casein enzymatic digest is a rich source of peptides and amino acids. Hence, the candy waste expressed as reducing sugar (g/L) (10, 20), fish extract expressed as nitrogen content (g/L) (0, 16), peptone from casein enzymatic digest (Fluka) (g/L) (0, 10), yeast extract (Sigma-Aldrich) (g/L) (0, 5), NaCl (g/L) (0, 4), MgSO4 (g/L) (0, 0.5), and MnSO4 (g/L) (0, 0.006) corresponding to (−1, +1) levels were screened for biofungicides production. Each of the 8 experiences reported in the PBD is repeated twice (Table 1). A linear approach presented by the following equation (Equation (1)) was considered to be sufficient for screening:
Y 1 = α 0 + α 1   A + α 2   B + α 3   C + α 4   D + α 5   E + α 6   F + α 7   G
where Y1 is the predicted target response (the antifungal activity); αi are the regression coefficients; and A, B, C, D, E, F, and G are dimensionless coded values of the independent variables candy waste, fish extract, peptone, yeast extract, NaCl, MgSO4, and MnSO4, respectively. The experimental data were fitted using MATLAB software. The most significant factors were then investigated more thoroughly in subsequent experiments (Table 2).

2.4. Optimization by Steepest Ascent Method

The path of the SAM was achieved to set up basal concentrations of media components selected from the PBD to be used in a CCD. It permitted rapid movement towards the most favorable of variable concentrations. Increments are direct ratios of regression coefficients αi (Equation (1)). Experiments were performed along with the SAM until the response did not increase anymore, and the starting point was the center of PBD, i.e., the medium level of factors reported in Table 1 and Table 3.

2.5. Central Composite Design and Response Surface

CCD Matrix and Antifungal Activity

Once the critical factors were identified via screening and the experimental design space was approached by SAM, the CCD was used to define the level of the significant parameters and the interactions between them, which significantly influence the antifungal activity. Each parameter was analyzed at five levels coded as (−2, −1, 0, +1, +2) (Table 4).

2.6. Regression Models and Statistical Analysis

The experimental data were fitted using MATLAB and Eureqa software. To determine the antifungal activity relation to input variables according to Table 4, a polynomial regression and then a trigonometric model were used. The second-order polynomial regression was first applied:
Y 1 = β 0 + i = 1 3 β i X i + i < j = 2 3 β i j X i X j + i = 1 3 β i i X i 2
where Y1 is the antifungal activity; Xi are input variables (3 variables retained); β 0 , β i , β i j , and β i i are, respectively, the regression coefficients for the intercept, linear, interaction, and quadratic effects. The model was then adjusted using the stepwise technique [42]. Given that the three retained variables took limited levels denoted by   X ¯ 1 ,   X ¯ 2 , and X ¯ 3 , an alternative trigonometric model with oscillating behavior was investigated. For a limited number of regression coefficients, we considered the following trigonometric model optimized by Eureqa software:
Y 1 = a + b   cos X ¯ 2   cos X ¯ 3 +   c   cos d e X ¯ 1
We were motivated by finding the model that guarantees high values of the coefficient of determination (R2), the adjusted R2 (AR2), and the corrected Akaike information criterion (AICc). The R2 increases every time we add predictors, even those insignificant. Therefore, the AR2 was introduced to compare the explanatory power of regression models and was correlated to R2 [43]. The AR2 increases only if the additional term is a good predictor, and it decreases with poor quality predictors. The Akaike information criterion rewards goodness of fit (using the likelihood); however, it penalizes increasing the number of estimated coefficients, which may differ from the number of predictors [44]. The AICc was introduced by Hurvich and Tsai [45] as a correction to the Akaike information criterion for small samples.

3. Results and Discussion

3.1. Screening of the Significant Medium Components Using PBD

The PBD was used to screen the seven factors candy waste, fish extract, peptone, yeast extract, NaCl, MgSO4, and MnSO4 with main effects on antifungal activity production by B. amyloliquefaciens BLB369 expressed in arbitrary unit per milliliter (AU/mL) (Table 1). The magnitudes and signs of each effect of the variables were shown. Candy waste (A), peptone (C), and NaCl (E) were the main variables that positively affected antifungal activity production as their p-values were lower than 1‰ (Table 2). The obtained linear model is given by
Y1 = 101.563 + 17.188 A + 48.438 C + 14.063 E
The peptone was the most significant factor due to its pronounced coefficient effect. It was manufactured by controlled enzymatic hydrolysis of casein, which is an excellent organic nitrogen source necessary for the bacteria to synthesize proteins and nucleic acids. Additionally, the industrial by-product candy waste constitutes a good carbon source for the growth of BLB369 strain and could then be valued for antifungal production. NaCl impacts the antifungal activity. It constitutes a nutriment for bacteria and has an osmotic effect in the culture medium. It could affect enzymatic activities and bacterial growth as it influences water and salts transport across bacterial membrane and could affect secretion and stability of the antifungal compounds. The fish extract did not affect the antifungal activity, probably due to fish extract instability. Therefore, the fish extract, yeast extract, MgSO4, and MnSO4 not retained after PBD screening were not added to the media in the coming runs as their concentrations at lower levels (−1) correspond to zero. The steepest ascent method further investigated the significant factors to optimize the experimental design space.

3.2. Optimization by Steepest Ascent Method

The direction of the steepest ascent method was determined by Equation (4). Concentrations of candy waste, peptone, and NaCl were increased because they presented positive effects on antifungal production. Hence, increments were direct ratios to regression coefficients αi, corresponding to 0.354 and 0.289 units of coded variables   X ¯ 1 and X ¯ 3 for each unit of X ¯ 2 , respectively. When the concentrations of candy waste, peptone, and NaCl were 20.25, 20, and 3.74 (g/L), respectively, the production of antifungal activity reached its maximal value of 250 (AU/mL) at run number 4 (Table 3). This point was chosen as a clue to set up basal concentrations for further optimization by CCD (the center point for optimization by CCD).
Table 3. Experimental design and response of the SAM experiments.
Table 3. Experimental design and response of the SAM experiments.
RunCandy WastePeptoneNaClAntifungal Activity (AU/mL)
X ¯ 1
0.354 #
X1 (g/L)
1.77 #
X ¯ 2 1 #X2 (g/L)
5 #
X ¯ 3
0.289 #
X3 (g/L)
1.58 #
1015.000502.00100
20.35416.771100.2892.58100
30.71018.552150.5783.16175
41.06520.253200.8673.74250
51.42022.004251.1564.30250
61.77523.755301.4454.90225
X ¯ i : variables in coded levels; Xi: variables in real values; #: increment.

3.3. Optimization of the Selected Medium Components Using the CCD

The CCD was used to define the optimum levels of the significant factors and study their interactions. The candy waste, peptone, and NaCl (denoted in what follows by X1, X2, X3, respectively) were studied at five levels (−2, −1, 0, +1, +2). The experimental design and the experimental responses of antifungal activity (Y1) were reported (Table 4).
Table 4. Response surface of CCD and results for antifungal activity.
Table 4. Response surface of CCD and results for antifungal activity.
RunCandy WastePeptoneNaClAntifungal Activity
X ¯ 1 X1
(g/L)
X ¯ 2 X2
(g/L)
X ¯ 3 X3
(g/L)
Y1
(AU/mL)
1+122−115+14.3175
2−118.4+125+14.3150
3−118.4−115−13.1175
4020.202003.7250
5020.202003.7250
6+122+125−13.1150
7−118.4−115+14.3175
8+122+125+14.3150
9−118.4+125−13.1175
10020.202003.7250
11+122−115−13.1125
12020.202003.7250
13020.2020−22.5125
14020.202003.7250
15+223.802003.7250
16020.2+23003.7125
17−216.602003.7200
18020.202003.7250
19020.2020+24.9125
20020.2−21003.7125

3.3.1. Regression Models for Antifungal Activity and Their Comparison

The regression coefficients were calculated using the least square technique and reported in Table 5. The standard errors and p-values were also reported. The standard error of the coefficient indicates the precision of the coefficient estimates. The p-value indicates the significance of each coefficient (low p-values indicate statistically significant terms). The significant variables are highlighted by asterisks: those with p-values less than 0.05 (Table 5). Moreover, to compare the resulting models, R2, AR2, AICc, and F-value for regression model and lack of fit are reported (Table 5 and Table 6).

Polynomial Regression

The polynomial fit shows only five significant terms (Table 5). However, before postulating such a statement, the model should be adjusted using the stepwise technique [42]: adding or removing variables one by one based on their p-value (the term with the highest p-value must be removed first). As additional terms, we tested including cubic (e.g., X13) and second-order interaction terms (e.g., X12X2). The retained polynomial model is given by the following equation (Table 5):
Y 1 = 12312 + 1624 X 1 + 53 X 2 + 351 X 3 78 X 1 2 1.33 X 2 2 + 1.25 X 1 3 8.3 X 3 3
Comparing results before and after applying the stepwise technique (Table 5), although a lower number of predictors (8 versus 10), the coefficient of determination is slightly improved after applying the stepwise technique from 0.92 to 0.928. The AR2 increased from 0.848 to 0.887, indicating more suitable predictors. The AICc criterion gets also improved (a lower value is observed of 190 versus 208). The F-value of the retained polynomial model was also computed (Table 6) and evaluated to 22, with a low p-value exhibiting the significance of the model. However, the hypothesis of lack of fit could not be rejected because of the high F-value associated with the lack of fit. This result is coherent with the relatively low AR2, and an alternative regression model should be examined.

Trigonometric Model

We started by estimating the coefficients through iterative optimization process. Then, we carried linear regression to check the coefficients of all terms and predict their significance (Table 5). The resulting equation is
Y 1 = 145 + 87   cos X ¯ 2   cos X ¯ 3 + 23.8   cos 5.5 2.3 X ¯ 1
The trigonometric model guarantees satisfactory values of both R2 and AR2 around 0.96 (96% of the variability in the response could be explained by this model), which are better than those obtained with the polynomial fit. The AICc criterion also shows the superiority of the trigonometric model with the lowest value of 160, which means a good compromise between the model fit and slight complexity. Moreover, the trigonometric model guarantees the lowest mean absolute error and root mean squared error showing a good agreement between the experimental and predicted values (Table 5). Finally, the ANOVA results demonstrate that the model is significant with an F-value of 212 and an associated p-value less than 10−12. The p-value associated with the lack of fit is equal to 0.14, indicating that the lack of fit is not significantly associated with the pure error (Table 6). All these statistical parameters illustrate the adequacy of the trigonometric model.

Predicted Versus Actual Plot and Residuals Versus Fits Plot

Figure 1 displays the predicted values by the retained polynomial and the trigonometric models, calculated from Equations (5) and (6), as a function of the actual values. Most markers (triangles for the polynomial model and squares for the trigonometric model) are scattered near the first bisector. This indicates that the predicted values are in good agreement with the actual ones. Figure 2 presents the internally studentized residuals versus the predicted values. All values are in theinterval (−3; +3) and almost equally scatter above and below the x-axis. Hence, both models produce random residuals, which means they are unbiased models.

Response Surface and Contour Plots

For a graphical illustration of the regression Equations (Equations (5) and (6)), contour plots and 3D response surfaces are represented (Figure 3 and Figure 4). In these plots, one of the three factors is set to its mean level (its value at the central point in the design space) while the two others are varied. Contour lines are concentric curves with increasing levels. Consequently, the antifungal activity should reach its maximal value for a specific combination of factors inside the considered domain.

Optimization of Production Conditions of the Antifungal Activity

For optimization of the antifungal activity, we carry a triple loop search sweeping the possible values of (X1, X2, X3) in the domain (18.4; 22) × (15; 25) × (3.1; 4.3) (g/L) with a step of 0.01 (g/L). The polynomial model (Equation (5)) yields a maximum of antifungal activity of 251.9 (AU/mL) for X1 = 19.17, X2 = 19.88, and X3 = 3.75 (g/L). Using the trigonometric model (Equation (6)), the maximal value of antifungal activity of 255.5 (AU/mL) is achieved for X1 = 19.61, X2 = 20, and X3 = 3.7 (g/L). The optimized values with both models are very close to the center point of CCD, and only a slight alteration of X1 should be applied.

4. Conclusions

PBD, SAM, and CCD were applied to maximize the antifungal production by B. amyloliquefaciens BLB369. Candy waste, peptone, and NaCl, the most significant among seven factors, were used for regression in CCD. Regression analysis showed the supremacy of a new trigonometric model over the usually used polynomial model, which may encourage using the trigonometric model in CDD. The residuals plot showed the adequacy of both models (unbiased). Moreover, contour plots and numerical optimization revealed a maximal activity for a point close to the CCD center. This low-cost medium could be investigated to further optimize BLB369 and similar strains for industrial and agricultural applications to control fungal diseases.

Author Contributions

Conceptualization, K.J.; methodology, K.J. and I.Z.-K.; software, S.K.; validation, I.Z.-K. and K.J.; formal analysis, S.K. and K.J.; investigation, I.Z.-K.; resources, S.T.; writing—original draft preparation, K.J. and S.K.; writing—review and editing, K.J. and S.K.; supervision, K.J.; project administration, K.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This study was supported by grants from the Tunisian “Ministry of Higher Education and Scientific Research”.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Kazan, K.; Gardiner, D.M. Transcriptomics of cereal-Fusarium graminearum interactions: What we have learned so far. Mol. Plant Pathol. 2018, 19, 764–778. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Rauwane, M.E.; Ogugua, U.V.; Kalu, C.M.; Ledwaba, L.K.; Woldesemayat, A.A.; Ntushelo, K. Pathogenicity and Virulence Factors of Fusarium graminearum Including Factors Discovered Using Next Generation Sequencing Technologies and Proteomics. Microorganisms 2020, 8, 305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pirgozliev, S.R.; Edwards, S.G.; Hare, M.C.; Jenkinson, P. Strategies for the control of Fusarium head blight in cereals. Eur. J. Plant Pathol. 2003, 109, 731–742. [Google Scholar] [CrossRef]
  4. Paul, P.A.; Lipps, P.E.; Hershman, D.E.; McMullen, M.P.; Draper, M.A.; Madden, L.V. Efficacy of triazole-based fungicides for fusarium head blight and deoxynivalenol control in wheat: A multivariate meta-analysis. Phytopathology 2008, 98, 999–1011. [Google Scholar] [CrossRef] [Green Version]
  5. Edwards, S.G.; Pirgozliev, S.R.; Hare, M.C.; Jenkinson, P. Quantification of trichothecene-producing Fusarium species in harvested grain by competitive PCR to determine efficacies of fungicides against fusarium head blight of winter wheat. Appl. Environ. Microbiol. 2001, 67, 1575–1580. [Google Scholar] [CrossRef] [Green Version]
  6. Paul, P.A.; McMullen, M.P.; Hershman, D.E.; Madden, L.V. Meta-analysis of the effects of triazole-based fungicides on wheat yield and test weight as influenced by fusarium head blight intensity. Phytopathology 2010, 100, 160–171. [Google Scholar] [CrossRef] [Green Version]
  7. Machado, F.J.; Santana, F.M.; Lau, D.; Del Ponte, E.M. Quantitative review of the effects of triazole and benzimidazole fungicides on fusarium head blight and wheat yield in Brazil. Plant Dis. 2017, 101, 1633–1641. [Google Scholar] [CrossRef] [Green Version]
  8. Kotowicz, N.K.; Frąc, M.; Lipiec, J. The importance of Fusarium fungi in wheat cultivation—Pathogenicity and mycotoxins production: A review. J. Anim. Plant Sci. 2014, 21, 3326–3343. [Google Scholar]
  9. Mateo, E.M.; Valle-Algarra, F.M.; Mateo, R.; Jiménez, M.; Magan, N. Effect of fenpropimorph, prochloraz and tebuconazole on growth and production of T-2 and HT-2 toxins by Fusarium langsethiae in oat-based medium. Int. J. Food Microbiol. 2011, 151, 289–298. [Google Scholar] [CrossRef] [Green Version]
  10. Kafle., S.; Vaidya, A.; Pradhan., B.; Jørs, E.; Onta, S. Factors Associated with Practice of Chemical Pesticide Use and Acute Poisoning Experienced by Farmers in Chitwan District Nepal. Int. J. Environ. Res. Public Health 2021, 18, 4194. [Google Scholar] [CrossRef]
  11. Roman, D.L.; Voiculescu, D.I.; Filip, M.; Ostafe, V.; Isvoran., A. Effects of Triazole Fungicides on Soil Microbiota and on the Activities of Enzymes Found in Soil: A Review. Agriculture 2021, 11, 893. [Google Scholar] [CrossRef]
  12. Kumar, J.; Ramlal, A.; Mallick, D.; Mishra, V. An Overview of Some Biopesticides and Their Importance in Plant Protection for Commercial Acceptance. Plants 2021, 10, 1185. [Google Scholar] [CrossRef] [PubMed]
  13. Asari, S.; Ongena, M.; Debois, D.; Pauw, E.; Chen, K.; Bejai, S.; Meijer, J. Insights into the molecular basis of biocontrol of Brassica pathogens by Bacillus amyloliquefaciens UCMB5113 lipopeptides. Ann. Bot. 2017, 20, 551–562. [Google Scholar] [CrossRef] [PubMed]
  14. Kunova, A.; Bonaldi, M.; Saracchi, M.; Pizzatti, C.; Chen, X.; Cortesi, P. Selection of Streptomyces against soil borne fungal pathogens by a standardized dual culture assay and evaluation of their effects on seed germination and plant growth. BMC Microbiol. 2016, 16, 272. [Google Scholar] [CrossRef] [Green Version]
  15. Gu, Q.; Yang, Y.; Yuan, Q.; Shi, G.; Wu, L.; Lou, Z.; Huo, R.; Wu, H.; Borriss, R.; Gao, X. Bacillomycin D Produced by Bacillus amyloliquefaciens Is Involved in the Antagonistic Interaction with the Plant-Pathogenic Fungus Fusarium graminearum. Appl. Environ. Microbiol. 2017, 83, e01075-17. [Google Scholar] [CrossRef] [Green Version]
  16. Penha, R.O.; Vandenberghe, L.P.S.; Faulds, C.; Soccol, V.; Soccol, C. Bacillus lipopeptides as powerful pest control agents for a more sustainable and healthy agriculture: Recent studies and innovations. Planta 2020, 251, 70. [Google Scholar] [CrossRef] [Green Version]
  17. Mora, I.; Cabrefiga, J.; Montesinos, E. Cyclic lipopeptide biosynthetic genes and products, and inhibitory activity of plant-associated Bacillus against phytopathogenic bacteria. PLoS ONE 2015, 10, e0127738. [Google Scholar] [CrossRef] [Green Version]
  18. Abdallah, R.A.B.; Stedel, C.; Garagounis, C.; Nefzi, A.; Jabnoun-Khiareddine, H.; Papadopoulou, K.K.; Daami-Remadi, M. Involvement of lipopeptide antibiotics and chitinase genes and induction of host defense in suppression of Fusarium wilt by endophytic Bacillus spp. in tomato. Crop Prot. 2017, 99, 45–58. [Google Scholar] [CrossRef]
  19. Han, J.H.; Shim, H.; Shin, J.H.; Kim, K.S. Antagonistic activities of Bacillus spp. strains isolated from tidal flat sediment towards anthracnose pathogens Colletotrichum acutatum and C. gloeosporioides in South Korea. Plant Pathol. J. 2015, 31, 165–175. [Google Scholar] [CrossRef] [Green Version]
  20. Hamid, R.; Khan, M.A.; Ahmad, M.; Ahmad, M.M.; Abdin, M.Z.; Musarrat, J.; Javed, S. Chitinases: An update. J. Pharm. Bioallied Sci. 2013, 5, 21–29. [Google Scholar]
  21. Zalila-Kolsi, I.; Sallemi, S.; Tounsi, S.; Jamoussi, K. Heterologous expression and periplasmic secretion of an antifungal Bacillus amyloliquefaciens BLB369 endo-β-1,3-1,4-glucanase in Escherichia coli. J. Phytopathol. 2017, 166, 28–33. [Google Scholar] [CrossRef]
  22. Mulatu, A.; Alemu, T.; Megersa, N.; Vetukuri, R.R. Optimization of Culture Conditions and Production of Bio-Fungicides from Trichoderma Species under Solid-State Fermentation Using Mathematical Modeling. Microorganisms 2021, 9, 1675. [Google Scholar] [CrossRef] [PubMed]
  23. Zhan, Y.; Zhu, P.; Liang, J.; Xu, Z.; Feng, X.; Liu, Y.; Xu, H.; Li, S. Economical production of isomaltulose from agricultural residues in a system with sucrose isomerase displayed on Bacillus subtilis spores. Bioprocess Biosyst. Eng. 2020, 43, 75–84. [Google Scholar] [CrossRef] [PubMed]
  24. Aguirre, A.M.; Bassi, A. Investigation of biomass concentration, lipid production, and cellulose content in Chlorella vulgaris cultures using response surface methodology. Biotechnol. Bioeng. 2013, 110, 2114–2122. [Google Scholar] [CrossRef]
  25. Bhatia, S.K.; Mehta, P.K.; Bhatia, R.K.; Bhalla, T.C. Optimization of arylacetonitrilase production from Alcaligenes sp. MTCC 10675 and its application in mandelic acid synthesis. Appl. Microbiol. Biotechnol. 2014, 98, 83–94. [Google Scholar] [CrossRef]
  26. Hsieh, S.C.; Liu, J.M.; Pua, X.H.; Ting, Y.; Hsu, R.J.; Cheng, K.C. Optimization of Lactobacillus acidophilus cultivation using taro waste and evaluation of its biological activity. Appl. Microbiol. Biotechnol. 2016, 100, 2629–2639. [Google Scholar] [CrossRef]
  27. Samaram, S.; Mirhosseini, H.; Tan, C.P.; Ghazali, H.M.; Bordbar, S.; Serjouie, A. Optimisation of ultrasound-assisted extraction of oil from papaya seed by response surface methodology: Oil recovery, radical scavenging antioxidant activity, and oxidation stability. Food Chem. 2015, 172, 7–17. [Google Scholar] [CrossRef]
  28. Samavati, V. Central composite rotatable design for investigation of microwave-assisted extraction of okra pod hydrocolloid. Int. J. Biol. Macromol. 2013, 61, 142–149. [Google Scholar] [CrossRef]
  29. Caldeira, A.T.; Arteiro, J.M.; Roseiro, J.C.; Neves, J.; Vicente, H. An artificial intelligence approach to Bacillus amyloliquefaciens CCMI 1051 cultures: Application to the production of anti-fungal compounds. Bioresour. Technol. 2011, 102, 1496–1502. [Google Scholar] [CrossRef]
  30. Tang, X.J.; He, G.Q.; Chen, Q.H.; Zhang, X.Y.; Ali, M.A. Medium optimization for the production of thermal stable β-glucanase by Bacillus subtilis ZJF-1A5 using response surface methodology. Bioresour. Technol. 2004, 93, 175–181. [Google Scholar] [CrossRef]
  31. Wang, L.; Zhang, M.; Li, Y.; Cui, Y.; Zhang, Y.; Wang, Z.; Wang, M.; Huang, Y. Application of response surface methodology to optimize the production of antimicrobial metabolites by Micromonospora Y15. Biotechnol. Biotechnol. Equip. 2017, 31, 1016–1025. [Google Scholar] [CrossRef] [Green Version]
  32. Jung, J.; Yu, K.O.; Ramzi, A.B.; Choe, S.H.; Kim, S.W.; Han, S.O. Improvement of surfactin production in Bacillus subtilis using synthetic wastewater by overexpression of specific extracellular signaling peptides, comX and phrC. Biotechnol. Bioeng. 2012, 109, 2349–2356. [Google Scholar] [CrossRef] [PubMed]
  33. Raza, W.; Hongsheng, W.; Qirong. Use of response surface methodology to evaluate the effect of metal ions (Ca2+, Ni2+, Mn2+, Cu2+) on production of antifungal compounds by Paenibacillus polymyxa. Bioresour. Technol. 2010, 101, 1904–1912. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, X.; Han, Y.; Tan, X.Q.; Wang, J.; Zhou, Z.J. Optimization of antifungal lipopeptide production from Bacillus sp. BH072 by response surface methodology. J. Microbiol. 2014, 52, 324–332. [Google Scholar] [CrossRef] [PubMed]
  35. Dai, J.; Wang, Z.; Xiu, Z. High production of optically pure (3R)-acetoin by a newly isolated marine strain of Bacillus subtilis CGMCC 13141. Bioprocess Biosyst. Eng. 2019, 42, 475–483. [Google Scholar] [CrossRef]
  36. De Oliveira, C.C.; de Souza, A.K.S.; de Castro, R.J.S. Bioconversion of Chicken Feather Meal by Aspergillus niger: Simultaneous Enzymes Production Using a Cost-Effective Feedstock Under Solid State Fermentation. Indian J. Microbiol. 2019, 59, 209–216. [Google Scholar] [CrossRef]
  37. Sreedharan, S.M.; Singh, S.P.; Singh, R. Flower Shaped Gold Nanoparticles: Biogenic Synthesis Strategies and Characterization. Indian J. Microbiol. 2019, 59, 321–327. [Google Scholar] [CrossRef]
  38. Zalila-Kolsi, I.; Ben Mahmoud, A.; Ali, H.; Sellami, S.; Nasfi, Z.; Tounsi, S.; Jamoussi, K. Antagonist effects of Bacillus spp. strains against Fusarium graminearum for protection of durum wheat (Triticum turgidum L. subsp. durum). Microbiol. Res. 2016, 192, 148–158. [Google Scholar] [CrossRef]
  39. Ahimou, F.; Jacques, P.; Daleu, M. Surfactin and iturin A effects on Bacillus subtilis surface hydrophobicity. Enzyme Microb. Technol. 2000, 27, 749–754. [Google Scholar] [CrossRef]
  40. Plackett, R.L.; Burman, J.P. The design of optimum multifactorial experiments. Biometrika 1946, 33, 305–325. [Google Scholar] [CrossRef]
  41. Miller, G.L. Use of dinitrosalicylic acid reagent for determination of reducing sugars. Anal. Chem. 1959, 31, 426–428. [Google Scholar] [CrossRef]
  42. Hocking, R.R. The analysis and selection of variables in linear regression. Biometrics 1976, 32, 1–49. [Google Scholar] [CrossRef]
  43. Theil, H. Economic Forecasts and Policy; North-Holland: Amsterdam, The Netherlands, 1961. [Google Scholar]
  44. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  45. Hurvich, C.M.; Tsai, C.L. Model selection for extended quasi-likelihood models in small samples. Biometrics 1995, 5, 1077–1084. [Google Scholar] [CrossRef]
Figure 1. Plot of the actual versus predicted values.
Figure 1. Plot of the actual versus predicted values.
Microorganisms 10 00830 g001
Figure 2. Plot of the internally studentized residuals versus the predicted values.
Figure 2. Plot of the internally studentized residuals versus the predicted values.
Microorganisms 10 00830 g002
Figure 3. Contour plots (a,c,e) and response surface curves (b,d,f) predicted by the retained polynomial model.
Figure 3. Contour plots (a,c,e) and response surface curves (b,d,f) predicted by the retained polynomial model.
Microorganisms 10 00830 g003
Figure 4. Contour plots (a,c,e) and response surface curves (b,d,f) predicted by the trigonometric model.
Figure 4. Contour plots (a,c,e) and response surface curves (b,d,f) predicted by the trigonometric model.
Microorganisms 10 00830 g004
Table 1. The Plackett–Burman experiments design matrix with factors given in coded levels and biofungicide production values.
Table 1. The Plackett–Burman experiments design matrix with factors given in coded levels and biofungicide production values.
RunA:
Candy Waste
B:
Fish Extract
C:
Peptone
D:
Yeast Extract
E:
NaCl
F:
MgSO4
G:
MnSO4
Antifungal
Activity
(AU/mL)
1+1−1−1+1−1+1+175
2+1+1−1−1+1−1+175
3+1+1+1−1−1+1−1150
4−1+1+1+1−1−1+1125
5+1−1+1+1+1−1−1175
6−1+1−1+1+1+1−162.5
7−1−1+1−1+1+1+1150
8−1−1−1−1−1−1−10
Variables in real values (g/L): A (10, 20), B (0, 16), C (0, 10), D (0, 5), E (0, 4), F (0, 0.5), and G (0, 0.006).
Table 2. Statistical analysis of factors using Plackett–Burman design.
Table 2. Statistical analysis of factors using Plackett–Burman design.
CoefficientValuep-ValueSignificance
α0101.5633.49 × 10−12***
α117.1884.15 × 10−6***
α21.5630.347
α348.4381.27 × 10−9***
α47.8131.05 × 10−3**
α514.0631.85 × 10−5***
α67.8131.05 × 10−3**
α74.6880.017*
*** Significance level 99.9%; ** significance level 99%; * significance level 95%. R2 = 0.9935; AR2 = 0.9878.
Table 5. Regression models: regression coefficients, their significance, and some statistical parameters.
Table 5. Regression models: regression coefficients, their significance, and some statistical parameters.
ModelTermCoefficientp-ValueSignificance
Polynomial model #Intercept−2050.50.03693*
X164.5340.2974
X257.4180.0151*
X35730.0062**
X12−2.54280.0703
X22−1.32959.75 × 10−6***
X32−92.339.75 × 10−6***
X1 X20.34720.6739
X1 X38.68060.2227
X2 X3−3.1250.2227
Retained polynomial model after applying the stepwise technique §Intercept−123120.0108*
X11623.80.0211*
X252.8497.74 × 10−7***
X3351.17.21 × 10−7***
X12−78.3060.0245*
X22−1.3296.38 × 10−7***
X131.25030.0285*
X33−8.30346.34 × 10−7***
Trigonometric model £Intercept144.523.1 × 10−19***
cos X ¯ 2   cos X ¯ 3 87.2026.6 × 10−13***
cos 5.5358 2.2737 X ¯ 1 23.7752.9 × 10−6***
*** Significance level 99.9%; ** significance level 99%; * significance level 95%. # R2 = 0.920; AR2 = 0.848; AICc = 208. § R2 = 0.928; AR2 = 0.887; AICc = 190. £ R2 = 0.962; AR2 = 0.957; AICc = 161.
Table 6. ANOVA for significance of regression and lack of fit of the retained polynomial model and the trigonometric model.
Table 6. ANOVA for significance of regression and lack of fit of the retained polynomial model and the trigonometric model.
Model Mean of SquareF-Valuep-Value
Retained polynomial modelTotal2728.6
Model6876.222.245.717 × 10−6
Residual309.17
Lack of fit530.01Inf0
Pure error0
Trigonometric modelTotal2728.6
Model24926212.719.32 × 10−13
Residual117.18
Lack of fit175.772.06240.1413
Pure error85.227
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zalila-Kolsi, I.; Kessentini, S.; Tounsi, S.; Jamoussi, K. Optimization of Bacillus amyloliquefaciens BLB369 Culture Medium by Response Surface Methodology for Low Cost Production of Antifungal Activity. Microorganisms 2022, 10, 830. https://doi.org/10.3390/microorganisms10040830

AMA Style

Zalila-Kolsi I, Kessentini S, Tounsi S, Jamoussi K. Optimization of Bacillus amyloliquefaciens BLB369 Culture Medium by Response Surface Methodology for Low Cost Production of Antifungal Activity. Microorganisms. 2022; 10(4):830. https://doi.org/10.3390/microorganisms10040830

Chicago/Turabian Style

Zalila-Kolsi, Imen, Sameh Kessentini, Slim Tounsi, and Kaïs Jamoussi. 2022. "Optimization of Bacillus amyloliquefaciens BLB369 Culture Medium by Response Surface Methodology for Low Cost Production of Antifungal Activity" Microorganisms 10, no. 4: 830. https://doi.org/10.3390/microorganisms10040830

APA Style

Zalila-Kolsi, I., Kessentini, S., Tounsi, S., & Jamoussi, K. (2022). Optimization of Bacillus amyloliquefaciens BLB369 Culture Medium by Response Surface Methodology for Low Cost Production of Antifungal Activity. Microorganisms, 10(4), 830. https://doi.org/10.3390/microorganisms10040830

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop