Articles Service
Research article
Optimization of Medium Composition for Biomass Production of Lactobacillus plantarum 200655 Using Response Surface Methodology
Department of Food Science and Biotechnology of Animal Resources, Konkuk University, Seoul 05029, Republic of Korea
Correspondence to:J. Microbiol. Biotechnol. 2021; 31(5): 717-725
Published May 28, 2021 https://doi.org/10.4014/jmb.2103.03018
Copyright © The Korean Society for Microbiology and Biotechnology.
Abstract
Keywords
Graphical Abstract
Introduction
The consumption of probiotics as functional foods or supplements has increased with the growing interest of consumers. Due to various health benefits such as balancing microbiota composition, improving immunity, protecting against intestinal pathogens, and controlling bowel diseases, the popularity of probiotics has grown remarkably in recent years [1]. As the importance of probiotics is consistently emphasized, research related to large-scale industrial production of probiotics has become necessary to satisfy the increasing demand.
Lactic acid bacteria (LAB), including the genera
Among the strains of lactobacilli, which consist of more than 50 species,
The commercial medium used for
To find effective ways to increase the yield of desired products, we optimized the medium using various strategies. One of the classical methods is one-factor-at-a-time (OFAT) in which only one factor is changed while the other factors are fixed constant. Although it is simple and convenient, OFAT involves a large number of experiments with an increasing number of factors. Moreover, the method disregards the interactions between factors [12, 13]. Thus, statistical methods such as Plackett-Burman design (PBD), factorial design, and response surface methodology (RSM) are preferred owing to their efficiency. Furthermore, RSM includes factorial designs and multiple regression analysis and considers the effects of factors, relationships between variables, and optimal conditions [14].
Although many studies have reported medium optimization of
Materials and Methods
Bacterial Strain and Medium
The basal medium was derived from MRS medium and consisted of 5 g/l sodium acetate (Sigma, USA), 2 g/l K2HPO4 (Samchun, Korea), 1 g/l Tween 80 (Yakuri Pure Chemicals Co. Ltd., Japan), 0.1 g/l MgSO4·7H2O (Shinyo Pure Chemicals Co., Japan), and 0.05 g/l MnSO4·H2O (Samchun). Basal medium was added to all media used in this study. The initial pH of all media was adjusted to 6.5 ± 0.05 with 1 M NaOH and 1 M HCl before sterilization. The MRS medium was used as a control (unoptimized medium) for comparison with the optimized medium.
Each fermentation condition for optimized medium was as follows: A colony was inoculated in MRS broth and incubated to an optical density (OD) of 0.5 ± 0.05 at 600 nm. One hundred milliliters of the medium were inoculated with 1% (v/v) of seed culture in a 250-ml Erlenmeyer flask and incubated at 37°C for 24 h without shaking.
Measurement of Biomass
Culture broth (50 ml) was collected and separated by centrifugation at 4,000 ×
Determination of Medium Compositions Using OFAT
To investigate the effects of carbon and nitrogen sources on bacterial cell mass, the OFAT method was used. This method can estimate the effects of each variable, but it does not consider the interactions between variables. Before proceeding with the OFAT method, the API 50 CHL Medium Kit (BioMerieux, France) was used to assess the availability of metabolizing carbohydrates. Glucose, sucrose, maltose, fructose, lactose, and galactose were selected, and six carbon sources were added individually to 100 ml of basal medium containing 10 g/l yeast extract at a concentration of 20 g/l. Likewise, six nitrogen sources (peptone, soytone, tryptone, yeast extract, beef extract, and malt extract) were added to 100 ml of basal medium containing 20 g/l glucose at a concentration of 10 g/l.
Plackett-Burman Design (PBD)
PBD was used to screen for significant factors affecting the biomass production by
-
Table 2 . Experimental design and response values of Plackett-Burman design.
Run Variablesa Biomass (g/l) X1 (g/l) X2 (g/l) X3 (g/l) X4 (g/l) X5 (g/l) X6 (g/l) 1 1 (30) -1 (10) 1 (30) -1 (5) -1 (5) -1 (5) 2.036 2 1 (30) 1 (30) -1 (10) 1 (10) -1 (5) -1 (5) 2.407 3 -1 (10) 1 (30) 1 (30) -1 (5) 1 (10) -1 (5) 2.109 4 1 (30) -1 (10) 1 (30) 1 (10) -1 (5) 1 (10) 2.407 5 1 (30) 1 (30) -1 (10) 1 (10) 1 (10) -1 (5) 2.561 6 1 (30) 1 (30) 1 (30) -1 (5) 1 (10) 1 (10) 2.360 7 -1 (10) 1 (30) 1 (30) 1 (10) -1 (5) 1 (10) 2.232 8 -1 (10) -1 (10) 1 (30) 1 (10) 1 (10) -1 (5) 2.181 9 -1 (10) -1 (10) -1 (10) 1 (10) 1 (10) 1 (10) 2.432 10 1 (30) -1 (10) -1 (10) -1 (5) 1 (10) 1 (10) 2.514 11 -1 (10) 1 (30) -1 (10) -1 (5) -1 (5) 1 (10) 2.122 12 -1 (10) -1 (10) -1 (10) -1 (5) -1 (5) -1 (5) 2.079 Actual values are presented in parentheses.
aX1, maltose; X2, sucrose; X3, lactose; X4, yeast extract; X5, soytone; X6, tryptone.
Central Composite Design and Response Surface Methodology
The RSM with CCD was conducted to optimize the concentration of medium components and estimate the effects of each variable and the interactions between variables. Maltose, yeast extract, and soytone were used as the independent variables. Variables were set at five different levels (-α, -1, 0, 1, α), and the range of actual values is presented in Table 4. The experimental design of the CCD was made up of a full factorial design with six center points (Table 5). The media were prepared according to a combination of variables in the experimental runs, and all experiments were conducted under static conditions. The obtained biomass was used to establish the regression model, and the quadratic polynomial equation for the variables was as follows:
-
Table 4 . Coded and real values of independent variables used in the central composite design.
Independent variables (g/l) Actual levels of coded values -α -1 0 1 α Maltose (X1) 8.18 15 25 35 41.8 Yeast extract (X2) 8.18 15 25 35 41.8 Soytone (X3) 4.77 15 30 45 55.2
-
Table 5 . Central composite design and response values.
Run Independent variablesa Biomass (g/l) X1 X2 X3 1 -1 -1 -1 2.484 2 1 -1 -1 3.216 3 -1 1 -1 2.620 4 1 1 -1 3.486 5 -1 -1 1 2.896 6 1 -1 1 3.714 7 -1 1 1 3.196 8 1 1 1 3.892 9 -α 0 0 2.348 10 α 0 0 3.570 11 0 - α 0 3.354 12 0 α 0 3.650 13 0 0 - α 3.232 14 0 0 α 3.546 15 0 0 0 3.846 16 0 0 0 3.684 17 0 0 0 3.610 18 0 0 0 3.794 19 0 0 0 3.810 20 0 0 0 3.728 aX1, maltose; X2, yeast extract; X3, soytone.
where Y is the response value of the dependent variable (biomass of
Effects of Culture pH and Temperature
The effects of initial pH and incubation temperature were investigated to determine the optimal cultivation conditions for
Scale-Up Fermentation of Optimized Medium
Scale-up fermentation of
Statistical Analysis
All experiments were repeated in triplicate. Data are presented as the mean ± SD. One-way analysis of variance and Duncan’s multiple range test were used to determine the degree of significant differences. Values were considered significant at
Results and Discussions
Effects of Carbon and Nitrogen Sources on Biomass Production by L. plantarum 200655
Various carbon and nitrogen sources were evaluated to identify the factors that have profound effects on biomass production by
Maltose utilization is related to maltose-related genes such as MALS (encoding maltase gene), MALT (encoding maltose permease gene), and MALR (transcriptional activator gene of the MALS and MALT). These results are similar to those of a study by Yeo
Table 1 also shows the biomass according to the six nitrogen sources. It was obvious that different nitrogen sources in the medium had different effects on biomass production. When malt extract was added to the medium,
-
Table 1 . Effects of various carbon and nitrogen sources on biomass production of
L. plantarum 200655.Components Biomass (g/l) Carbon sources Glucose 1.703 Sucrose 1.753 Maltose 2.253 Fructose 1.459 Lactose 1.749 Galactose 1.401 Nitrogen sources Peptone 1.175 Soytone 1.480 Tryptone 1.371 Yeast extract 1.722 Beef extract 0.495 Malt extract 0.022
Plackett-Burman Design for Screening Variables
Because PBD is an efficient tool for estimating the main effect of each variable, it was applied to identify the most significant variables before performing RSM. Maltose, sucrose, lactose, yeast extract, soytone, and tryptone were selected as variables for PBD. They were set at two different levels (-1, 1), where the actual values were 10 to 30 g/l and 5 to 10 g/l respectively. The six variables with 12 experimental runs resulted in biomass production ranging from 2.036 to 2.561 g/l (Table 2). The experiments were conducted in triplicate in accordance with combinations of variables. The standardized effects of the variables are represented as a single column on the Pareto chart (Fig. 1). The vertical line through the column indicates whether the variables are statistically significant. Where the columns crosses over the line and extends to the right shows that the variables had a large effect on the biomass. Table 3 shows the effects of the variables and analysis of variance. The effects of maltose (X1), sucrose (X2), lactose (X3), yeast extract (X4), soytone (X5), and tryptone (X6) were 0.1884, 0.0236, -0.1318, 0.1667, 0.1456, and 0.1158, respectively. Only lactose had a negative effect on biomass production among the six variables, suggesting that biomass production was decreased by increasing the concentration of lactose in the medium in the tested concentration of 10 g/l to 30 g/l. In contrast, maltose had the highest effect, followed by the yeast extract and soytone. The analysis of variance showed that all components except sucrose were statistically significant (
-
Table 3 . Analysis of variables based on Plackett-Burman design.
Variables Effect Coefficient T -valueP -valueIntercept 2.2867 149.56 < 0.0001 X1 0.1884 0.0942 6.16 0.0016 X2 0.0236 0.0118 0.77 0.4759 X3 -0.1318 -0.0659 -4.31 0.0076 X4 0.1667 0.0833 5.45 0.0028 X5 0.1456 0.0728 4.76 0.0051 X6 0.1158 0.0579 3.79 0.0128
-
Fig. 1.
Pareto chart depicting the standard effects of six factors on the biomass production of L. plantarum 200655.
Optimization of Medium Components Using Response Surface Methodology
The CCD model was designed based on RSM to determine the optimal concentrations of maltose (X1, 8.18 g/l to 41.8 g/l), yeast extract (X2, 8.18 g/l to 41.8 g/l), and soytone (X3, 4.77 g/l to 55.2 g/l). The coded units and actual concentrations of the three independent variables are listed in Table 4. Table 5 shows the design matrix consisting of 20 runs and experimental responses, where biomass production varied from 2.348 to 3.892 g/l. The highest biomass was obtained in run 8 with a concentration of 35 g/l maltose, 35 g/l yeast extract, and 45 g/l soytone. The responses of biomass were analyzed by applying multiple regression analysis, and the second-order polynomial equation, which expressed the relationship between the predicted response and variables, was as follows:
where X1, X2, and X3 are maltose, yeast extract, and soytone, respectively.
The analysis of variance presented in Table 6 shows the significance of the quadratic regression model with linear, squared, and interaction terms. As shown in Table 6, a high
-
Table 6 . Analysis of variance of the response surface quadratic model.
Source DFa Adj SSb Adj MSc F -valueP -valueModel 9 3.99061 0.44340 33.22 < 0.0001 X1 1 1.95502 1.95502 146.46 < 0.0001 X2 1 0.13981 0.13981 10.47 0.0089 X3 1 0.42885 0.42885 32.13 0.0002 X1X2 1 0.00002 0.00002 0.00 0.9714 X1X3 1 0.00088 0.00088 0.07 0.8024 X2X3 1 0.00065 0.00065 0.05 0.8301 X12 1 1.23057 1.23057 92.18 < 0.0001 X22 1 0.14479 0.14479 10.85 0.0081 X32 1 0.28321 0.28321 21.22 < 0.0010 Residual 10 0.13349 0.01335 Lack of fit 5 0.09443 0.01889 2.42 0.1774 Pure error 5 0.03906 0.00781 Total 19 4.12410 R2 = 0.9676; adjusted R2 = 0.9385; R = 0.9837.
aDegree of freedom.
bSum of squares.
cMean square.
Three-dimensional response surface and contour plots were drawn to express the interactions between the two variables and derive the optimal concentration for maximal biomass production (Fig. 2). Graphical representations having a convex-shaped response surface were depicted based on the model equation. Each graph showed the infinite combinations of two independent variables, with the other one at a constant level. In the optimization of the medium for maximal actinorhodin production by
-
Fig. 2.
Response surface plots and contour plots for biomass production of (L. plantarum 200655.A, B ) Interaction between maltose (X1, g/l) and yeast extract (X2, g/l) with soytone (X3, g/l) at zero level. (C, D ) Interaction between maltose (X1, g/l) and soytone (X3, g/l) with yeast extract (X2, g/l) at zero level. (E, F ) Interaction between yeast extract (X2, g/l) and soytone (X3, g/l) with maltose (X1, g/l) at zero level.
The predicted model was validated under the optimized conditions by performing independent experiments in triplicate, and the biomass was found to be 3.845 g/l. The model showed good agreement with the predicted value, and the experimental value was within the 95% confidence interval range. The biomass in the optimized medium was 1.58-fold higher compared to that in the unoptimized medium (2.429 g/l). Manzoor
Effects of pH of Medium and Fermentation Temperature on Biomass
The effects of the initial pH and incubation temperature were investigated in the optimized medium since not only the medium formula but also the physicochemical parameters are important for bacterial growth. It seems that the pH and temperature affected the biomass production by
-
Fig. 3.
Effects of initial pH (A) and incubation temperature (B) on biomass production by Data are presented as the mean ± SD of independent experiments in triplicate. Different superscript letters of each figure are significantly different (L. plantarum 200655.p < 0.05).
The effects of temperature were also evaluated using the optimized medium. The maximal biomass (4.505 g/l) was attained at 30°C, which was statistically significant compared to that at other temperatures (
Comparison of Optimized and Unoptimized Media in a Bioreactor
Fig. 4 presents the time course of
-
Fig. 4.
The time course of biomass production and viable cells of L. plantarum 200655 in optimized medium (●, ■) and unoptimized medium (○, □).
Therefore, biomass production by
Acknowledgments
This work was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, and Forestry (IPET) through the Innovational Food Technology Development Program (#1190093), funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA).
Conflict of Interest
The authors have no financial conflicts of interest to declare.
References
- Tripathi MK, Giri SK. 2014. Probiotic functional foods: Survival of probiotics during processing and storage.
J. Funct. Foods 9 : 225-241. - Nataraj BH, Shivanna SK, Rao P, Nagpal R, Behare PV. 2020. Evolutionary concepts in the functional biotics arena: a mini-review.
Food Sci. Biotechnol. 30 : 487-496. - Lee NK, Kim SY, Han KJ, Eom SJ, Paik HD. 2014. Probiotic potential of
Lactobacillus strains with anti-allergic effects from kimchi for yogurt starters.LWT-Food Sci. Technol. 58 : 130-134. - Son SH, Yang SJ, Jeon HL, Yu HS, Lee NK, Park YS,
et al . 2018. Antioxidant and immunostimulatory effect of potential probioticLactobacillus paraplantarum SC61 isolated from Korean traditional fermented food, jangajji.Microb. Pathog. 125 : 486-492. - Asan-Ozusaglam M, Gunyakti A. 2019.
Lactobacillus fermentum strains from human breast milk with probiotic properties and cholesterol-lowering effects.Food Sci. Biotechnol. 28 : 501-509. - de Vries MC, Vaughan EE, Kleerebezem M, de Vos WM. 2006.
Lactobacillus plantarum - survival, functional and potential probiotic properties in the human intestinal tract.Int. Dairy J. 16 : 1018-1028. - Yan SJ, Lee JE, Lim SM, Kim YJ, Lee NK, Paik HD. 2019. Antioxidant and immune-enhancing effects of probiotic
Lactobacillus plantarum 200655 isolated from kimchi.Food Sci. Biotechnol. 28 : 491-499. - Cheon MJ, Lee NK, Paik HD. 2021. Neuroprotective effects of heat-killed
Lactobacillus plantarum 200655 isolated from kimchi against oxidative stress.Probiotics Antimicrob. Proteins . doi.org/10.1007/s12602-020-09740-w. - Yeo S, Shin HS, Lee HW, Hong D, Park H, Holzapfel W,
et al . 2018. Determination of optimized growth medium and cryoprotective additives to enhance the growth and survival ofLactobacillus salivarius .J. Microbiol. Biotechnol. 28 : 718-731. - Brinques GB, Peralba MC, Ayub MAZ. 2010. Optimization of probiotic and lactic acid production by
Lactobacillus plantarum in submerged bioreactor systems.J. Ind. Microbiol. 37 : 205-212. - Manzoor A, Qazi JI, Haq IU, Mukhtar H, Rasool A. 2017. Significantly enhanced biomass production of a novel bio-therapeutic strain
Lactobacillus plantarum (AS-14) by developing low-cost media cultivation strategy.J. Biol. Eng. 11 : 17. - Lee YM, Kim JS, Kim WJ. 2012. Optimization for the maximum bacteriocin production of
Lactobacillus brevis DF01 using response surface methodology.Food Sci. Biotechnol. 21 : 653-659. - Singh V, Haque S, Niwas R, Srivastava A, Pasupuleti M, Tripathi CKM. 2017. Strategies for fermentation medium optimization: an in-depth review.
Front. Microbiol. 7 : 2087. - Coelho LF, de Lima CJB, Rodovalho CM, Bernardo MP, Contiero J. 2011. Lactic acid production by new
Lactobacillus plantarum LMISM6 grown in molasses: optimization of medium composition.Braz. J. Chem. Eng. 28 : 27-36. - Liu B, Yang M, Qi B, Chen X, Su Z, Wan Y. 2010. Optimizing L-(+)-lactic acid production by thermophile
Lactobacillus plantarum As.1.3 using alternative nitrogen sources with response surface method.Biochem. Eng. J. 52 : 212-219. - Leal-Sánchez MV, Jiménez-Díaz R, Maldonado-Barragán A, Garrido-Fernández A, Ruiz-Barba JL. 2002. Optimization of bacteriocin production by batch fermentation of
Lactobacillus plantarum LPCO10.Appl. Environ. Microbiol. 68 : 4465-4471. - Desai KM, Akolkar SK, Badhe YP, Tambe SS, Lele SS. 2006. Optimization of fermentation media for exopolysaccharide production from
Lactobacillus plantarum using artificial intelligence-based techniques.Process Biochem. 41 : 1842-1848. - Othman NZ, Mohd Din ARJ, Azam Z, Rosli MA, Sarmidi MR. 2018. Statistical optimization of medium compositions for high cell mass and exopolysaccharide production by
Lactobacillus plantarum ATCC 8014.Appl. Food Biotechnol. 5 : 87-96. - Hwang CF, Chang JH, Houng JY, Tsai CC, Lin CK, Tsen HY. 2012. Optimization of medium composition for improving biomass production of
Lactobacillus plantarum Pi06 using the Taguchi array design and the Box-Behnken method.Biotechnol. Bioprocess Eng. 17 : 827-834. - Lim HS, Cha IT, Roh SW, Shin HH, Seo MJ. 2017. Enhanced production of gamma-aminobutyric acid by optimizing culture conditions of
Lactobacillus brevis HYE1 isolated from kimchi, a Korean fermented food.J. Miocrobiol. Biotechnol. 27 : 450-459. - Charalampopoulos D, Pandiella SS, Webb C. 2003. Evaluation of the effect of malt, wheat and barley extracts on the viability of potentially probiotic lactic acid bacteria under acidic conditions.
Int. J. Food Microbiol. 82 : 133-141. - Aasen IM, Møretrø T, Katla T, Axelsson L, Storrø I. 2000. Influence of complex nutrients, temperature and pH on bacteriocin production by
Lactobacillus sakei CCUG 42687.Appl. Microbiol. Biotechnol. 53 : 159-166. - Lee NK, Park YL, Choe GJ, Chang HI, Paik HD. 2010. Medium optimization for the production of probiotic
Lactobacillus acidophilus A12 using response surface methodology.Korean J. Food Sci. An. 30 : 359-364. - Selvamani S, Dailin DJ, Rostom M, Malek RA, Gupta VK, El-Enshasy HA. 2020. Optimizing medium components to enhance high cell mass production of biotherapeutic strain
Lactobacillus reuteri DSM 20016T by statistical method.J. Sci. Ind. Res. 79 : 798-803. - Elibol M. 2004. Optimization of medium composition for actinorhodin production by
Streptomyces coelicolor A3(2) with response surface methodology.Process Biochem. 39 : 1057-1062. - Hayek SA, Ibrahim SA. 2013. Current limitations and challenges with lactic acid bacteria: a review.
Food Nutr. Sci. 4 : 73-87. - Mbye M, Baig MA, AbuQamar SF, El-Tarabily KA, Obaid RS, Osaili TM,
et al . 2020. Updates on understanding of probiotic lactic acid bacteria responses to environmental stresses and highlights on proteomic analyses.Compr. Rev. Food Sci. Food Saf. 19 : 1110-1124. - Coghetto CC, Vasconcelos CB, Brinques GB, Ayub MAZ. 2016.
Lactobacillus plantarum BL011 cultivation in industrial isolated soybean protein acid residue.Braz. J. Microbiol. 47 : 941-948. - Fonteles TV, Costa MGM, de Jesus ALT, Rodrigues S. 2012. Optimization of the fermentation of cantaloupe juice by
Lactobacillus casei NRRL B-442.Food Bioprocess Tech. 5 : 2819-2826.
Related articles in JMB
Article
Research article
J. Microbiol. Biotechnol. 2021; 31(5): 717-725
Published online May 28, 2021 https://doi.org/10.4014/jmb.2103.03018
Copyright © The Korean Society for Microbiology and Biotechnology.
Optimization of Medium Composition for Biomass Production of Lactobacillus plantarum 200655 Using Response Surface Methodology
Ga-Hyun Choi, Na-Kyoung Lee, and Hyun-Dong Paik*
Department of Food Science and Biotechnology of Animal Resources, Konkuk University, Seoul 05029, Republic of Korea
Correspondence to:Hyun-Dong Paik, hdpaik@konkuk.ac.kr
Abstract
This study aimed to optimize medium composition and culture conditions for enhancing the biomass of Lactobacillus plantarum 200655 using statistical methods. The one-factor-at-a-time (OFAT) method was used to screen the six carbon sources (glucose, sucrose, maltose, fructose, lactose, and galactose) and six nitrogen sources (peptone, tryptone, soytone, yeast extract, beef extract, and malt extract). Based on the OFAT results, six factors were selected for the Plackett- Burman design (PBD) to evaluate whether the variables had significant effects on the biomass. Maltose, yeast extract, and soytone were assessed as critical factors and therefore applied to response surface methodology (RSM). The optimal medium composition by RSM was composed of 31.29 g/l maltose, 30.27 g/l yeast extract, 39.43 g/l soytone, 5 g/l sodium acetate, 2 g/l K2HPO4, 1 g/l Tween 80, 0.1 g/l MgSO4·7H2O, and 0.05 g/l MnSO4·H2O, and the maximum biomass was predicted to be 3.951 g/l. Under the optimized medium, the biomass of L. plantarum 200655 was 3.845 g/l, which was similar to the predicted value and 1.58-fold higher than that of the unoptimized medium (2.429 g/l). Furthermore, the biomass increased to 4.505 g/l under optimized cultivation conditions. For lab-scale bioreactor validation, batch fermentation was conducted with a 5-L bioreactor containing 3.5 L of optimized medium. As a result, the highest yield of biomass (5.866 g/l) was obtained after 18 h of incubation at 30°C, pH 6.5, and 200 rpm. In conclusion, mass production by L. plantarum 200655 could be enhanced to obtain higher yields than that in MRS medium
Keywords: Probiotics, Lactobacillus plantarum, medium optimization, plackett-burman design, response surface methodology
Introduction
The consumption of probiotics as functional foods or supplements has increased with the growing interest of consumers. Due to various health benefits such as balancing microbiota composition, improving immunity, protecting against intestinal pathogens, and controlling bowel diseases, the popularity of probiotics has grown remarkably in recent years [1]. As the importance of probiotics is consistently emphasized, research related to large-scale industrial production of probiotics has become necessary to satisfy the increasing demand.
Lactic acid bacteria (LAB), including the genera
Among the strains of lactobacilli, which consist of more than 50 species,
The commercial medium used for
To find effective ways to increase the yield of desired products, we optimized the medium using various strategies. One of the classical methods is one-factor-at-a-time (OFAT) in which only one factor is changed while the other factors are fixed constant. Although it is simple and convenient, OFAT involves a large number of experiments with an increasing number of factors. Moreover, the method disregards the interactions between factors [12, 13]. Thus, statistical methods such as Plackett-Burman design (PBD), factorial design, and response surface methodology (RSM) are preferred owing to their efficiency. Furthermore, RSM includes factorial designs and multiple regression analysis and considers the effects of factors, relationships between variables, and optimal conditions [14].
Although many studies have reported medium optimization of
Materials and Methods
Bacterial Strain and Medium
The basal medium was derived from MRS medium and consisted of 5 g/l sodium acetate (Sigma, USA), 2 g/l K2HPO4 (Samchun, Korea), 1 g/l Tween 80 (Yakuri Pure Chemicals Co. Ltd., Japan), 0.1 g/l MgSO4·7H2O (Shinyo Pure Chemicals Co., Japan), and 0.05 g/l MnSO4·H2O (Samchun). Basal medium was added to all media used in this study. The initial pH of all media was adjusted to 6.5 ± 0.05 with 1 M NaOH and 1 M HCl before sterilization. The MRS medium was used as a control (unoptimized medium) for comparison with the optimized medium.
Each fermentation condition for optimized medium was as follows: A colony was inoculated in MRS broth and incubated to an optical density (OD) of 0.5 ± 0.05 at 600 nm. One hundred milliliters of the medium were inoculated with 1% (v/v) of seed culture in a 250-ml Erlenmeyer flask and incubated at 37°C for 24 h without shaking.
Measurement of Biomass
Culture broth (50 ml) was collected and separated by centrifugation at 4,000 ×
Determination of Medium Compositions Using OFAT
To investigate the effects of carbon and nitrogen sources on bacterial cell mass, the OFAT method was used. This method can estimate the effects of each variable, but it does not consider the interactions between variables. Before proceeding with the OFAT method, the API 50 CHL Medium Kit (BioMerieux, France) was used to assess the availability of metabolizing carbohydrates. Glucose, sucrose, maltose, fructose, lactose, and galactose were selected, and six carbon sources were added individually to 100 ml of basal medium containing 10 g/l yeast extract at a concentration of 20 g/l. Likewise, six nitrogen sources (peptone, soytone, tryptone, yeast extract, beef extract, and malt extract) were added to 100 ml of basal medium containing 20 g/l glucose at a concentration of 10 g/l.
Plackett-Burman Design (PBD)
PBD was used to screen for significant factors affecting the biomass production by
-
Table 2 . Experimental design and response values of Plackett-Burman design..
Run Variablesa Biomass (g/l) X1 (g/l) X2 (g/l) X3 (g/l) X4 (g/l) X5 (g/l) X6 (g/l) 1 1 (30) -1 (10) 1 (30) -1 (5) -1 (5) -1 (5) 2.036 2 1 (30) 1 (30) -1 (10) 1 (10) -1 (5) -1 (5) 2.407 3 -1 (10) 1 (30) 1 (30) -1 (5) 1 (10) -1 (5) 2.109 4 1 (30) -1 (10) 1 (30) 1 (10) -1 (5) 1 (10) 2.407 5 1 (30) 1 (30) -1 (10) 1 (10) 1 (10) -1 (5) 2.561 6 1 (30) 1 (30) 1 (30) -1 (5) 1 (10) 1 (10) 2.360 7 -1 (10) 1 (30) 1 (30) 1 (10) -1 (5) 1 (10) 2.232 8 -1 (10) -1 (10) 1 (30) 1 (10) 1 (10) -1 (5) 2.181 9 -1 (10) -1 (10) -1 (10) 1 (10) 1 (10) 1 (10) 2.432 10 1 (30) -1 (10) -1 (10) -1 (5) 1 (10) 1 (10) 2.514 11 -1 (10) 1 (30) -1 (10) -1 (5) -1 (5) 1 (10) 2.122 12 -1 (10) -1 (10) -1 (10) -1 (5) -1 (5) -1 (5) 2.079 Actual values are presented in parentheses..
aX1, maltose; X2, sucrose; X3, lactose; X4, yeast extract; X5, soytone; X6, tryptone..
Central Composite Design and Response Surface Methodology
The RSM with CCD was conducted to optimize the concentration of medium components and estimate the effects of each variable and the interactions between variables. Maltose, yeast extract, and soytone were used as the independent variables. Variables were set at five different levels (-α, -1, 0, 1, α), and the range of actual values is presented in Table 4. The experimental design of the CCD was made up of a full factorial design with six center points (Table 5). The media were prepared according to a combination of variables in the experimental runs, and all experiments were conducted under static conditions. The obtained biomass was used to establish the regression model, and the quadratic polynomial equation for the variables was as follows:
-
Table 4 . Coded and real values of independent variables used in the central composite design..
Independent variables (g/l) Actual levels of coded values -α -1 0 1 α Maltose (X1) 8.18 15 25 35 41.8 Yeast extract (X2) 8.18 15 25 35 41.8 Soytone (X3) 4.77 15 30 45 55.2
-
Table 5 . Central composite design and response values..
Run Independent variablesa Biomass (g/l) X1 X2 X3 1 -1 -1 -1 2.484 2 1 -1 -1 3.216 3 -1 1 -1 2.620 4 1 1 -1 3.486 5 -1 -1 1 2.896 6 1 -1 1 3.714 7 -1 1 1 3.196 8 1 1 1 3.892 9 -α 0 0 2.348 10 α 0 0 3.570 11 0 - α 0 3.354 12 0 α 0 3.650 13 0 0 - α 3.232 14 0 0 α 3.546 15 0 0 0 3.846 16 0 0 0 3.684 17 0 0 0 3.610 18 0 0 0 3.794 19 0 0 0 3.810 20 0 0 0 3.728 aX1, maltose; X2, yeast extract; X3, soytone..
where Y is the response value of the dependent variable (biomass of
Effects of Culture pH and Temperature
The effects of initial pH and incubation temperature were investigated to determine the optimal cultivation conditions for
Scale-Up Fermentation of Optimized Medium
Scale-up fermentation of
Statistical Analysis
All experiments were repeated in triplicate. Data are presented as the mean ± SD. One-way analysis of variance and Duncan’s multiple range test were used to determine the degree of significant differences. Values were considered significant at
Results and Discussions
Effects of Carbon and Nitrogen Sources on Biomass Production by L. plantarum 200655
Various carbon and nitrogen sources were evaluated to identify the factors that have profound effects on biomass production by
Maltose utilization is related to maltose-related genes such as MALS (encoding maltase gene), MALT (encoding maltose permease gene), and MALR (transcriptional activator gene of the MALS and MALT). These results are similar to those of a study by Yeo
Table 1 also shows the biomass according to the six nitrogen sources. It was obvious that different nitrogen sources in the medium had different effects on biomass production. When malt extract was added to the medium,
-
Table 1 . Effects of various carbon and nitrogen sources on biomass production of
L. plantarum 200655..Components Biomass (g/l) Carbon sources Glucose 1.703 Sucrose 1.753 Maltose 2.253 Fructose 1.459 Lactose 1.749 Galactose 1.401 Nitrogen sources Peptone 1.175 Soytone 1.480 Tryptone 1.371 Yeast extract 1.722 Beef extract 0.495 Malt extract 0.022
Plackett-Burman Design for Screening Variables
Because PBD is an efficient tool for estimating the main effect of each variable, it was applied to identify the most significant variables before performing RSM. Maltose, sucrose, lactose, yeast extract, soytone, and tryptone were selected as variables for PBD. They were set at two different levels (-1, 1), where the actual values were 10 to 30 g/l and 5 to 10 g/l respectively. The six variables with 12 experimental runs resulted in biomass production ranging from 2.036 to 2.561 g/l (Table 2). The experiments were conducted in triplicate in accordance with combinations of variables. The standardized effects of the variables are represented as a single column on the Pareto chart (Fig. 1). The vertical line through the column indicates whether the variables are statistically significant. Where the columns crosses over the line and extends to the right shows that the variables had a large effect on the biomass. Table 3 shows the effects of the variables and analysis of variance. The effects of maltose (X1), sucrose (X2), lactose (X3), yeast extract (X4), soytone (X5), and tryptone (X6) were 0.1884, 0.0236, -0.1318, 0.1667, 0.1456, and 0.1158, respectively. Only lactose had a negative effect on biomass production among the six variables, suggesting that biomass production was decreased by increasing the concentration of lactose in the medium in the tested concentration of 10 g/l to 30 g/l. In contrast, maltose had the highest effect, followed by the yeast extract and soytone. The analysis of variance showed that all components except sucrose were statistically significant (
-
Table 3 . Analysis of variables based on Plackett-Burman design..
Variables Effect Coefficient T -valueP -valueIntercept 2.2867 149.56 < 0.0001 X1 0.1884 0.0942 6.16 0.0016 X2 0.0236 0.0118 0.77 0.4759 X3 -0.1318 -0.0659 -4.31 0.0076 X4 0.1667 0.0833 5.45 0.0028 X5 0.1456 0.0728 4.76 0.0051 X6 0.1158 0.0579 3.79 0.0128
-
Figure 1.
Pareto chart depicting the standard effects of six factors on the biomass production of L. plantarum 200655.
Optimization of Medium Components Using Response Surface Methodology
The CCD model was designed based on RSM to determine the optimal concentrations of maltose (X1, 8.18 g/l to 41.8 g/l), yeast extract (X2, 8.18 g/l to 41.8 g/l), and soytone (X3, 4.77 g/l to 55.2 g/l). The coded units and actual concentrations of the three independent variables are listed in Table 4. Table 5 shows the design matrix consisting of 20 runs and experimental responses, where biomass production varied from 2.348 to 3.892 g/l. The highest biomass was obtained in run 8 with a concentration of 35 g/l maltose, 35 g/l yeast extract, and 45 g/l soytone. The responses of biomass were analyzed by applying multiple regression analysis, and the second-order polynomial equation, which expressed the relationship between the predicted response and variables, was as follows:
where X1, X2, and X3 are maltose, yeast extract, and soytone, respectively.
The analysis of variance presented in Table 6 shows the significance of the quadratic regression model with linear, squared, and interaction terms. As shown in Table 6, a high
-
Table 6 . Analysis of variance of the response surface quadratic model..
Source DFa Adj SSb Adj MSc F -valueP -valueModel 9 3.99061 0.44340 33.22 < 0.0001 X1 1 1.95502 1.95502 146.46 < 0.0001 X2 1 0.13981 0.13981 10.47 0.0089 X3 1 0.42885 0.42885 32.13 0.0002 X1X2 1 0.00002 0.00002 0.00 0.9714 X1X3 1 0.00088 0.00088 0.07 0.8024 X2X3 1 0.00065 0.00065 0.05 0.8301 X12 1 1.23057 1.23057 92.18 < 0.0001 X22 1 0.14479 0.14479 10.85 0.0081 X32 1 0.28321 0.28321 21.22 < 0.0010 Residual 10 0.13349 0.01335 Lack of fit 5 0.09443 0.01889 2.42 0.1774 Pure error 5 0.03906 0.00781 Total 19 4.12410 R2 = 0.9676; adjusted R2 = 0.9385; R = 0.9837..
aDegree of freedom..
bSum of squares..
cMean square..
Three-dimensional response surface and contour plots were drawn to express the interactions between the two variables and derive the optimal concentration for maximal biomass production (Fig. 2). Graphical representations having a convex-shaped response surface were depicted based on the model equation. Each graph showed the infinite combinations of two independent variables, with the other one at a constant level. In the optimization of the medium for maximal actinorhodin production by
-
Figure 2.
Response surface plots and contour plots for biomass production of (L. plantarum 200655.A, B ) Interaction between maltose (X1, g/l) and yeast extract (X2, g/l) with soytone (X3, g/l) at zero level. (C, D ) Interaction between maltose (X1, g/l) and soytone (X3, g/l) with yeast extract (X2, g/l) at zero level. (E, F ) Interaction between yeast extract (X2, g/l) and soytone (X3, g/l) with maltose (X1, g/l) at zero level.
The predicted model was validated under the optimized conditions by performing independent experiments in triplicate, and the biomass was found to be 3.845 g/l. The model showed good agreement with the predicted value, and the experimental value was within the 95% confidence interval range. The biomass in the optimized medium was 1.58-fold higher compared to that in the unoptimized medium (2.429 g/l). Manzoor
Effects of pH of Medium and Fermentation Temperature on Biomass
The effects of the initial pH and incubation temperature were investigated in the optimized medium since not only the medium formula but also the physicochemical parameters are important for bacterial growth. It seems that the pH and temperature affected the biomass production by
-
Figure 3.
Effects of initial pH (A) and incubation temperature (B) on biomass production by Data are presented as the mean ± SD of independent experiments in triplicate. Different superscript letters of each figure are significantly different (L. plantarum 200655.p < 0.05).
The effects of temperature were also evaluated using the optimized medium. The maximal biomass (4.505 g/l) was attained at 30°C, which was statistically significant compared to that at other temperatures (
Comparison of Optimized and Unoptimized Media in a Bioreactor
Fig. 4 presents the time course of
-
Figure 4.
The time course of biomass production and viable cells of L. plantarum 200655 in optimized medium (●, ■) and unoptimized medium (○, □).
Therefore, biomass production by
Acknowledgments
This work was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, and Forestry (IPET) through the Innovational Food Technology Development Program (#1190093), funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA).
Conflict of Interest
The authors have no financial conflicts of interest to declare.
Fig 1.
Fig 2.
Fig 3.
Fig 4.
-
Table 1 . Effects of various carbon and nitrogen sources on biomass production of
L. plantarum 200655..Components Biomass (g/l) Carbon sources Glucose 1.703 Sucrose 1.753 Maltose 2.253 Fructose 1.459 Lactose 1.749 Galactose 1.401 Nitrogen sources Peptone 1.175 Soytone 1.480 Tryptone 1.371 Yeast extract 1.722 Beef extract 0.495 Malt extract 0.022
-
Table 2 . Experimental design and response values of Plackett-Burman design..
Run Variablesa Biomass (g/l) X1 (g/l) X2 (g/l) X3 (g/l) X4 (g/l) X5 (g/l) X6 (g/l) 1 1 (30) -1 (10) 1 (30) -1 (5) -1 (5) -1 (5) 2.036 2 1 (30) 1 (30) -1 (10) 1 (10) -1 (5) -1 (5) 2.407 3 -1 (10) 1 (30) 1 (30) -1 (5) 1 (10) -1 (5) 2.109 4 1 (30) -1 (10) 1 (30) 1 (10) -1 (5) 1 (10) 2.407 5 1 (30) 1 (30) -1 (10) 1 (10) 1 (10) -1 (5) 2.561 6 1 (30) 1 (30) 1 (30) -1 (5) 1 (10) 1 (10) 2.360 7 -1 (10) 1 (30) 1 (30) 1 (10) -1 (5) 1 (10) 2.232 8 -1 (10) -1 (10) 1 (30) 1 (10) 1 (10) -1 (5) 2.181 9 -1 (10) -1 (10) -1 (10) 1 (10) 1 (10) 1 (10) 2.432 10 1 (30) -1 (10) -1 (10) -1 (5) 1 (10) 1 (10) 2.514 11 -1 (10) 1 (30) -1 (10) -1 (5) -1 (5) 1 (10) 2.122 12 -1 (10) -1 (10) -1 (10) -1 (5) -1 (5) -1 (5) 2.079 Actual values are presented in parentheses..
aX1, maltose; X2, sucrose; X3, lactose; X4, yeast extract; X5, soytone; X6, tryptone..
-
Table 3 . Analysis of variables based on Plackett-Burman design..
Variables Effect Coefficient T -valueP -valueIntercept 2.2867 149.56 < 0.0001 X1 0.1884 0.0942 6.16 0.0016 X2 0.0236 0.0118 0.77 0.4759 X3 -0.1318 -0.0659 -4.31 0.0076 X4 0.1667 0.0833 5.45 0.0028 X5 0.1456 0.0728 4.76 0.0051 X6 0.1158 0.0579 3.79 0.0128
-
Table 4 . Coded and real values of independent variables used in the central composite design..
Independent variables (g/l) Actual levels of coded values -α -1 0 1 α Maltose (X1) 8.18 15 25 35 41.8 Yeast extract (X2) 8.18 15 25 35 41.8 Soytone (X3) 4.77 15 30 45 55.2
-
Table 5 . Central composite design and response values..
Run Independent variablesa Biomass (g/l) X1 X2 X3 1 -1 -1 -1 2.484 2 1 -1 -1 3.216 3 -1 1 -1 2.620 4 1 1 -1 3.486 5 -1 -1 1 2.896 6 1 -1 1 3.714 7 -1 1 1 3.196 8 1 1 1 3.892 9 -α 0 0 2.348 10 α 0 0 3.570 11 0 - α 0 3.354 12 0 α 0 3.650 13 0 0 - α 3.232 14 0 0 α 3.546 15 0 0 0 3.846 16 0 0 0 3.684 17 0 0 0 3.610 18 0 0 0 3.794 19 0 0 0 3.810 20 0 0 0 3.728 aX1, maltose; X2, yeast extract; X3, soytone..
-
Table 6 . Analysis of variance of the response surface quadratic model..
Source DFa Adj SSb Adj MSc F -valueP -valueModel 9 3.99061 0.44340 33.22 < 0.0001 X1 1 1.95502 1.95502 146.46 < 0.0001 X2 1 0.13981 0.13981 10.47 0.0089 X3 1 0.42885 0.42885 32.13 0.0002 X1X2 1 0.00002 0.00002 0.00 0.9714 X1X3 1 0.00088 0.00088 0.07 0.8024 X2X3 1 0.00065 0.00065 0.05 0.8301 X12 1 1.23057 1.23057 92.18 < 0.0001 X22 1 0.14479 0.14479 10.85 0.0081 X32 1 0.28321 0.28321 21.22 < 0.0010 Residual 10 0.13349 0.01335 Lack of fit 5 0.09443 0.01889 2.42 0.1774 Pure error 5 0.03906 0.00781 Total 19 4.12410 R2 = 0.9676; adjusted R2 = 0.9385; R = 0.9837..
aDegree of freedom..
bSum of squares..
cMean square..
References
- Tripathi MK, Giri SK. 2014. Probiotic functional foods: Survival of probiotics during processing and storage.
J. Funct. Foods 9 : 225-241. - Nataraj BH, Shivanna SK, Rao P, Nagpal R, Behare PV. 2020. Evolutionary concepts in the functional biotics arena: a mini-review.
Food Sci. Biotechnol. 30 : 487-496. - Lee NK, Kim SY, Han KJ, Eom SJ, Paik HD. 2014. Probiotic potential of
Lactobacillus strains with anti-allergic effects from kimchi for yogurt starters.LWT-Food Sci. Technol. 58 : 130-134. - Son SH, Yang SJ, Jeon HL, Yu HS, Lee NK, Park YS,
et al . 2018. Antioxidant and immunostimulatory effect of potential probioticLactobacillus paraplantarum SC61 isolated from Korean traditional fermented food, jangajji.Microb. Pathog. 125 : 486-492. - Asan-Ozusaglam M, Gunyakti A. 2019.
Lactobacillus fermentum strains from human breast milk with probiotic properties and cholesterol-lowering effects.Food Sci. Biotechnol. 28 : 501-509. - de Vries MC, Vaughan EE, Kleerebezem M, de Vos WM. 2006.
Lactobacillus plantarum - survival, functional and potential probiotic properties in the human intestinal tract.Int. Dairy J. 16 : 1018-1028. - Yan SJ, Lee JE, Lim SM, Kim YJ, Lee NK, Paik HD. 2019. Antioxidant and immune-enhancing effects of probiotic
Lactobacillus plantarum 200655 isolated from kimchi.Food Sci. Biotechnol. 28 : 491-499. - Cheon MJ, Lee NK, Paik HD. 2021. Neuroprotective effects of heat-killed
Lactobacillus plantarum 200655 isolated from kimchi against oxidative stress.Probiotics Antimicrob. Proteins . doi.org/10.1007/s12602-020-09740-w. - Yeo S, Shin HS, Lee HW, Hong D, Park H, Holzapfel W,
et al . 2018. Determination of optimized growth medium and cryoprotective additives to enhance the growth and survival ofLactobacillus salivarius .J. Microbiol. Biotechnol. 28 : 718-731. - Brinques GB, Peralba MC, Ayub MAZ. 2010. Optimization of probiotic and lactic acid production by
Lactobacillus plantarum in submerged bioreactor systems.J. Ind. Microbiol. 37 : 205-212. - Manzoor A, Qazi JI, Haq IU, Mukhtar H, Rasool A. 2017. Significantly enhanced biomass production of a novel bio-therapeutic strain
Lactobacillus plantarum (AS-14) by developing low-cost media cultivation strategy.J. Biol. Eng. 11 : 17. - Lee YM, Kim JS, Kim WJ. 2012. Optimization for the maximum bacteriocin production of
Lactobacillus brevis DF01 using response surface methodology.Food Sci. Biotechnol. 21 : 653-659. - Singh V, Haque S, Niwas R, Srivastava A, Pasupuleti M, Tripathi CKM. 2017. Strategies for fermentation medium optimization: an in-depth review.
Front. Microbiol. 7 : 2087. - Coelho LF, de Lima CJB, Rodovalho CM, Bernardo MP, Contiero J. 2011. Lactic acid production by new
Lactobacillus plantarum LMISM6 grown in molasses: optimization of medium composition.Braz. J. Chem. Eng. 28 : 27-36. - Liu B, Yang M, Qi B, Chen X, Su Z, Wan Y. 2010. Optimizing L-(+)-lactic acid production by thermophile
Lactobacillus plantarum As.1.3 using alternative nitrogen sources with response surface method.Biochem. Eng. J. 52 : 212-219. - Leal-Sánchez MV, Jiménez-Díaz R, Maldonado-Barragán A, Garrido-Fernández A, Ruiz-Barba JL. 2002. Optimization of bacteriocin production by batch fermentation of
Lactobacillus plantarum LPCO10.Appl. Environ. Microbiol. 68 : 4465-4471. - Desai KM, Akolkar SK, Badhe YP, Tambe SS, Lele SS. 2006. Optimization of fermentation media for exopolysaccharide production from
Lactobacillus plantarum using artificial intelligence-based techniques.Process Biochem. 41 : 1842-1848. - Othman NZ, Mohd Din ARJ, Azam Z, Rosli MA, Sarmidi MR. 2018. Statistical optimization of medium compositions for high cell mass and exopolysaccharide production by
Lactobacillus plantarum ATCC 8014.Appl. Food Biotechnol. 5 : 87-96. - Hwang CF, Chang JH, Houng JY, Tsai CC, Lin CK, Tsen HY. 2012. Optimization of medium composition for improving biomass production of
Lactobacillus plantarum Pi06 using the Taguchi array design and the Box-Behnken method.Biotechnol. Bioprocess Eng. 17 : 827-834. - Lim HS, Cha IT, Roh SW, Shin HH, Seo MJ. 2017. Enhanced production of gamma-aminobutyric acid by optimizing culture conditions of
Lactobacillus brevis HYE1 isolated from kimchi, a Korean fermented food.J. Miocrobiol. Biotechnol. 27 : 450-459. - Charalampopoulos D, Pandiella SS, Webb C. 2003. Evaluation of the effect of malt, wheat and barley extracts on the viability of potentially probiotic lactic acid bacteria under acidic conditions.
Int. J. Food Microbiol. 82 : 133-141. - Aasen IM, Møretrø T, Katla T, Axelsson L, Storrø I. 2000. Influence of complex nutrients, temperature and pH on bacteriocin production by
Lactobacillus sakei CCUG 42687.Appl. Microbiol. Biotechnol. 53 : 159-166. - Lee NK, Park YL, Choe GJ, Chang HI, Paik HD. 2010. Medium optimization for the production of probiotic
Lactobacillus acidophilus A12 using response surface methodology.Korean J. Food Sci. An. 30 : 359-364. - Selvamani S, Dailin DJ, Rostom M, Malek RA, Gupta VK, El-Enshasy HA. 2020. Optimizing medium components to enhance high cell mass production of biotherapeutic strain
Lactobacillus reuteri DSM 20016T by statistical method.J. Sci. Ind. Res. 79 : 798-803. - Elibol M. 2004. Optimization of medium composition for actinorhodin production by
Streptomyces coelicolor A3(2) with response surface methodology.Process Biochem. 39 : 1057-1062. - Hayek SA, Ibrahim SA. 2013. Current limitations and challenges with lactic acid bacteria: a review.
Food Nutr. Sci. 4 : 73-87. - Mbye M, Baig MA, AbuQamar SF, El-Tarabily KA, Obaid RS, Osaili TM,
et al . 2020. Updates on understanding of probiotic lactic acid bacteria responses to environmental stresses and highlights on proteomic analyses.Compr. Rev. Food Sci. Food Saf. 19 : 1110-1124. - Coghetto CC, Vasconcelos CB, Brinques GB, Ayub MAZ. 2016.
Lactobacillus plantarum BL011 cultivation in industrial isolated soybean protein acid residue.Braz. J. Microbiol. 47 : 941-948. - Fonteles TV, Costa MGM, de Jesus ALT, Rodrigues S. 2012. Optimization of the fermentation of cantaloupe juice by
Lactobacillus casei NRRL B-442.Food Bioprocess Tech. 5 : 2819-2826.