BioMed Research International

Volume 2014 (2014), Article ID 787989, 10 pages

http://dx.doi.org/10.1155/2014/787989

Research Article

Optimization of Milk-Based Medium for Efficient Cultivation of Bifidobacterium pseudocatenulatum G4 Using Face-Centered Central Composite-Response Surface Methodology

Khalilah Abdul Khalil,¹ Shuhaimi Mustafa,2,3 Rosfarizan Mohammad,3,4 Arbakariya Bin Ariff,³ Yamin Shaari,³ Yazid Abdul Manap,⁵ Siti Aqlima Ahmad,³ and Farrah Aini Dahalan⁶

¹Department of Biomolecular Sciences, Faculty of Applied Sciences, Universiti Teknologi MARA, Sec. 2, 40150 Shah Alam, Selangor, Malaysia

²Halal Products Research Institute, Universiti Putra Malaysia, Putra Infoport, 43400 Serdang, Selangor, 40150 Shah Alam, Selangor, Malaysia

³Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, 43400 13 Serdang, Selangor, Malaysia

⁴Institute of Tropical Forestry and Forest Products, Universiti Putra Malaysia, 43400 15 Serdang, Selangor, Malaysia

⁵Faculty of Food Science and Technology, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

⁶The School of Environmental Engineering, Universiti Malaysia Perlis, Kompleks Pengajian Jejawi 3, 02600 Arau, Perlis, Malaysia

Received 25 July 2013; Accepted 13 November 2013; Published 12 January 2014

Academic Editor: Filippo Canducci

Copyright © 2014 Khalilah Abdul Khalil et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study was undertaken to optimize skim milk and yeast extract concentration as a cultivation medium for optimal Bifidobacteria pseudocatenulatum G4 (G4) biomass and β-galactosidase production as well as lactose and free amino nitrogen (FAN) balance after cultivation period. Optimization process in this study involved four steps: screening for significant factors using 2³ full factorial design, steepest ascent, optimization using FCCD-RSM, and verification. From screening steps, skim milk and yeast extract showed significant influence on the biomass production and, based on the steepest ascent step, middle points of skim milk (6% wt/vol) and yeast extract (1.89% wt/vol) were obtained. A polynomial regression model in FCCD-RSM revealed that both factors were found significant and the strongest influence was given by skim milk concentration. Optimum concentrations of skim milk and yeast extract for maximum biomass G4 and β-galactosidase production meanwhile low in lactose and FAN balance after cultivation period were 5.89% (wt/vol) and 2.31% (wt/vol), respectively. The validation experiments showed that the predicted and experimental values are not significantly different, indicating that the FCCD-RSM model developed is sufficient to describe the cultivation process of G4 using skim-milk-based medium with the addition of yeast extract.

1. Introduction

Probiotic can be defined as “living microorganisms which when administered in adequate amounts confer a health benefit on the host” [1]. The beneficial effects of probiotic include regulation in immune system, anticarcinogenic, preventing the infection by exogenous microorganisms, vitamin synthesis, and enhancing digestion and absorption of nutrients [2].

Bifidobacteria, the major species of intestinal flora in healthy human, have been widely used in probiotics formulation [3]. Bifidobacteria are anaerobic rods, Gram-positive, and nonspore forming bacteria [4]. Bifidobacterium pseudocatenulatum G4, isolated from breast-fed infants feces, has been identified and characterized as a potential probiotic with many beneficial health effects [5–7]. The antimicrobial properties, adherence ability, antagonistic activity, and ability to survive in gut environment for Bifidobacterium pseudocatenulatum G4 have also been reported [8].

Proliferation of bifidobacteria requires appropriate medium which supplies nurturing substances for growth. The use of commercial media to support bacterial growth is limited by its cost [9] and off-flavor generated in food products [10]. Milk is a highly nutritious growth medium for many microorganisms as it is rich in carbohydrates, fat, casein, protein, vitamins, and minerals [11]. However, some microorganisms are unable to grow well in this medium when specific enzymes required for lactose metabolization process are lacking [12]. β-Galactosidase is an enzyme required to break lactose to galactose and glucose. Most of the bifidobacteria preferred these carbon sources to promote their growth. The production of β-galactosidase can be enhanced using yeast extract as a nitrogen source [13]. Various concentrations of yeast extract (ranging from 0.3% to 1.0% wt/vol) were added into milk (ranging from 10% to 12% wt/vol) to enhance growth of probiotics [14, 15]. Various concentrations of glucose (ranging from 0.5% to 2.0% wt/vol) have also been used as a carbon source in the cultivation of probiotic like bifidobacteria [16].

In optimization of milk-based-medium for bifidobacteria growth, it is important to consider β-galactosidase production by the cells in order to facilitate lactose breakdown into more readily available sugar. In addition, the concomitant use of nitrogen sources is also important to avoid nutrients wastage. Most of the previous studies on the medium optimization for bifidobacteria focused mainly on the maximum biomass production [8, 17]. To the best of our knowledge, justification involving other responses of growth activities of bifidobacteria in milk medium in order to ensure that the medium has been used competently is not available in the literature.

The optimization of medium concentration or components by conventional approach involves large numbers of experiments which are time consuming and expensive and lack in representing the effects of individual factors. Thus, a statistical approach known as response surface methodology (RSM) has been applied as a tool for process improvement. The main advantage of RSM is the reduced number of experiments or trials required to evaluate multiple parameters and their interactions [18, 19]. The most commonly design approached is the central composite design. The design allows estimation of all the regression parameters required to fit a second-order model of given responses. Rotability character is the most preferred in any central composite design. This is because this characteristic provides constant variance of the estimated response corresponding to all new observation points that are at the same distance from the center point of the design (in terms of the coded variable).

Therefore, the main objective of this study was to optimize milk-based medium for the improvement of B. pseudocatenulatum G4 cultivation. In this optimization process, there were four responses; namely, (i) cells production, (ii) accumulation of lactose degrading enzyme (β-galactosidase), and (iii) residual concentration of lactose and free amino nitrogen remaining in the culture at the end of cultivation, were considered.

2. Materials and Methods

2.1. Bacterium and Inoculum Preparation

Probiotic bacterium, B. pseudocatenulatum G4 (G4), was used throughout this study and this bacterium was obtained from the Probiotic Laboratory, Faculty of Food Science and Technology, Universiti Putra Malaysia. G4 was previously isolated from breast-fed infant feces [5, 6] and was stored at −20°C in a mixture of glycerol and Trypticase Phytone Yeast (TPY) extract broth (Scharlau-Chemie, Barcelona, Spain), at a ratio of 80 : 20. TPY medium (Scharlau-Chemie, Barcelona, Spain) was used to maintain and to propagate the bacterium [4]. A single colony of the bacterium was transferred from TPY agar to TPY broth and incubated anaerobically at 37°C for 24 h. After two successive transfers in TPY broth at 37°C for 24 h under anaerobic condition using Anaerocult A gas packs (Merck, Darmstadt, Germany), the activated culture was then properly diluted and served as the standard inoculum for all cultivations.

2.2. Medium Preparation

Cultivation was carried out in 250 mL SCHOTT DURAN bottle (Schott Duran, Mainz, Germany). The medium components consist of skim milk (NZMP medium heat skim milk powder, Auckland, New Zealand), yeast extract (Bio Springer, Maisons-Alfort Cedex, France), and glucose (Merck, Darmstadt, Germany). All components were prepared in separate bottle and sterilized at 121°C for 15 min. Once cooled at room temperature, the components were mixed prior to inoculation with 10⁴ cfu/mL of G4. The bottle was then incubated at 37°C for 21 h in an anaerobic condition. The culture pH was not controlled in all experiments. During the cultivation, 10 mL of samples was taken at 3 h intervals for analysis.

2.3. Lactose Determination

The amount of lactose in the samples was determined using HPLC method as described by Hou et al. with some modifications [20]. One mL of sample was centrifuged at 13000 rpm for 10 min. The clear fraction was filtered through a 0.2 μm nylon membrane filter and injected onto HPLC system (Alliance 2690/5: Waters Corporation, Milford, CA), equipped with a 4.6 mm × 150 mm Agilent Zorbax Carbohydrate Analysis column (Agilent Technologies Inc., USA). The mobile phase used is 75% (v/v) acetonitrile (Fisher, HPLC grade). The flow rate was set at 1.4 mL min⁻¹ and the analysis was carried out at 30°C using a refractive index detector (RI-1371, Waters Corporation, Milford, CA).

2.4. Free Amino Nitrogen Analyses

Free amino nitrogen (FAN) presence in the sample was determined using ninhydrin analysis described by Hwang and Ederer with some modifications [21]. A sample (1 mL) was centrifuged at 13000 rpm for 10 min. The clear fraction was collected and diluted into 50 mL distilled water and 2 mL of the diluted sample was then transferred into test tubes. Ninhydrin solution (0.35 g ninhydrin into 100 mL of a 1 : 1 mixture of acetone and butanol) was then added and the tubes were covered with a piece of parafilm to avoid the loss of solvent due to evaporation. The tubes were then heated in a boiling water bath (80°C–100°C) for 15 min before being transferred to a cold water bath. The dilution reagent (5 mL) was then added and mixed and the absorbance was read at 570 nm using spectrophotometer (UV-Vis 1601 spectrophotometer, Shimadzu, Japan). Distilled water was used as a blank and glycine (Merck, Darmstadt, Germany) was used as a standard.

2.5. β-Galactosidase Analyses

β-Galactosidase was assayed on nonpermeabilized and permeabilized cells. The culture was centrifuged at 11000 rpm for 10 min at 4°C. Subsequently, the resulting supernatant was discarded by aspiration. The cell pellet was washed twice with Z buffer (0.1 M phosphate—pH 7; 10 mM MgSO₄·7H₂O; 1 mM CaCl₂) and concentrated by 10-fold. Permeabilization was accomplished by the addition of 0.5 mL of Triton X-100 (5% vol/vol in Z buffer) into 0.5 mL cell suspension with mixing using vortex mixer. The mixture was incubated at 37°C for 10 min and then centrifuged at 14000 rpm for 15 min and the resulting supernatant served as the enzyme source. The activity of β-galactosidase was then assayed essentially according to the method described by Nagy et al. (2001) [22]. The reaction mixture was composed of 0.5 mL of enzyme sample and 0.5 mL of 15 mM o-nitrophenyl β-d-galactopyranoside (Sigma Chemical Co., St. Louis, MO) in 0.03 M sodium phosphate buffer (pH 6.8). Once a faint yellow tint appeared, the reaction was terminated by adding 0.1 M sodium carbonate. The absorbance was measured using a spectrophotometer at 420 nm wavelength (UV-Vis 1601 spectrophotometer, Shimadzu, Japan). One unit of β-galactosidase activity was defined as the amount of enzyme catalyzing the formation of 1 μmol of o-nitrophenyl per min under the assay condition.

2.6. Microbiological Analyses

For viable cells enumeration, samples were serially diluted using 0.1% (wt/vol) sterile peptone water (Merck, Darmstadt, Germany) and plated in duplicate onto TPY agar. The plates were incubated in anaerobic jars containing Anaerocult A (Merck, Darmstadt, Germany) at 37°C for 48 h. All plates with 30 to 300 colonies were counted and the viability was expressed as log₁₀ cfu/mL.

2.7. Experimental Design and Statistical Analyses

Design-Expert software (Stat-Ease Inc., Minneapolis, MN, USA, version 6.0.6) was used to design the experiment in this study. First, initial screening experiments were performed to evaluate the significant effects of skim milk, yeast extract, and glucose on biomass production using a 2³ full factorial design. Each factor was coded at three levels between −1 and +1, with 5 times replication of center point (0), where the variables of factors were changed in the ranges. The factorial design resulting in 26 experimental runs (including duplicates and 5 middle points run). In the second step, a first-order empirical equation was applied to exclude insignificant factors and to generate the steepest ascent path to facilitate the maximum increase in response. The steepest ascent design was based on results of the equation from the screening step of which only factors with significant effects on response were used and fixed at their zero coded value. The insignificant factors were eliminated from the design. The path begins at the centre point of design and stretches out linearly in order to determine the suitable ranges within significant factors which tend to approach the optimal condition. The third step involved further optimization of significant factors using face-centered central composite design (FCCD). FCCD is an effective design that is used for sequential experimentation and provides reasonable amount of information for testing the goodness of fit and does not require large number of design points thereby reducing the overall cost associated with the experiment [23]. In this design, four responses, namely, biomass (, log₁₀ cfu/mL) and β-galactosidase production, (, U/mL) as well as residual lactose (, g/L) and FAN (, mg/L) remaining in the culture were determined. There were three coded factor levels, −1, 0, and +1 (where −1 corresponded to the low level of each factor, 0 to the middle level, and +1 to the high level). The coded values were determined using

The center point was repeated five times in order to evaluate the curvature and the experiment replication facilitated the pure error estimation, so that the significant lack of fit of the models could be predicted. All the 26 experiments were carried out in duplicate for 21 h.

3. Results and Discussion

3.1. Initial Screening of Significant Medium Components and the Steepest Ascent

The effect of skim milk (), yeast extract (), and glucose () on biomass production () during the cultivation of G4 using 2³ full factorial design is shown in Table 1. The maximum biomass (ranging from 4.495 to 7.638 log₁₀ cfu/mL) was obtained at 21 h of cultivation. The analysis of variance (ANOVA) of the first-order model is shown in Table 2 while the regression analysis of the model is shown in Table 3. The model was significant and only 1.12% of the total variation was observed due to noise. The model was linear with insignificant () curvature. The regression analysis of the variables showed that skim milk () was the only significant factor while yeast extract and glucose were insignificant with values higher than 0.05 (Table 3). Since yeast extract () was the sole nitrogen, it must be considered in the medium optimization. Improvement of growth of bifidobacteria could be achieved using yeast extract at concentration higher than 1.0% without the addition of glucose [17]. As a result, final-order equation (coded term) was generated based on the first-order model to determine the biomass production response () to the medium composition containing skim milk () and yeast extract () factors, which give

For every unit increase in , an increase of 1.09 units in is predicted (2). On the other hand, for every unit increase in , 0.30 units of increase in is expected. Between these two factors, skim milk () gave more effect than yeast extract () on biomass production (). This equation was further used as the fundamental scale in the subsequent step which was the steepest ascent. The path of the steepest ascent was determined based on the increase in 0.10% (wt/vol) of and the movement was generated along the path until no improvement occurred. Five design units were developed based on 0.10/0.10 = 1. Hence, the movement of was 0.28 design units (). As for natural factor, movement was based on coded factor × lower point of (). Thus, movement of natural factor involved . The path coordination of the steepest ascent was generated and is shown in Table 4. The highest biomass production can be observed in the third step of the steepest ascent coordinates with the value of log₁₀ cfu/mL from the combination of skim milk (6% (wt/vol)) and yeast extract (1.89% (wt/vol)). After the third step of coordinates, declining in biomass production was observed. Consequently, this combination was selected as the middle point for further optimization.

Table 1: The 2³ full factorial design and responses for screening experiments.

Run
Factors
Response

Skim milk (% wt/vol)
Yeast extract (% wt/vol)
Glucose (% wt/vol)
Maximum biomass (log₁₀ CFU per milliliter)

1
2 (−1)
0.5 (−1)
0.5 (−1)
4.495

2
2 (−1)
0.5 (−1)
1.5 (+1)
4.789

3
2 (−1)
3 (+1)
0.5 (−1)
5.335

4
2 (−1)
3 (+1)
1.5 (+1)
5.800

5
4 (0)
1.75 (0)
1 (0)
6.417

6
4 (0)
1.75 (0)
1 (0)
6.671

7
4 (0)
1.75 (0)
1 (0)
6.075

8
4 (0)
1.75 (0)
1 (0)
6.824

9
4 (0)
1.75 (0)
1 (0)
6.874

10
6 (+1)
0.5 (−1)
0.5 (−1)
7.158

11
6 (+1)
0.5 (−1)
1.5 (+1)
7.173

12
6 (+1)
3 (+1)
0.5 (−1)
7.204

13
6 (+1)
3 (+1)
1.5 (+1)
7.638

Table 2: ANOVA results of the first-order model for 2³ full factorial design.

Source of variation or factors
DFa
Sum of squares
Mean square
or value
value

Modelb
7
10.75
1.54
14.09
0.0112

Curvature
7
0.43
0.43
0.071
0.1182

Pure error
4
0.44
0.11
3.94

Correlation total
12
11.61

DF: degree of freedom.

b: 96.10%.

Significant at .

Table 3: Regression analysis of the 2³ full factorial design with maximum biomass (log₁₀ cfu/mL) as the response.

Variable
Value
value

87.94
0.0007

6.40
0.0647

1.67
0.2658

2.06
0.2246

0.11
0.7572

0.40
0.5623

0.07
0.8034

Model
14.09
0.0112

Curvature
3.94
0.1182

C.V = 2.78%

Adjusted

  represents skim milk concentration in % (wt/vol).

represents yeast extract concentration % (wt/vol).

represents glucose concentration % (wt/vol).

Table 4: The steepest ascent coordination path for all chosen factors at coded and natural levels.

Step
Coded factora
Natural factorb (% wt/vol)
Maximum biomass (log₁₀ cfu/mL)e

Basec
0
0
4
1.75

Δd
1
0.28
2
0.14

Base + Δ
1
0.28
6
1.89

Base + 2Δ
2
0.56
8
2.03

Base + 3Δ
3
0.84
10
2.17

Base + 4Δ
4
1.12
12
2.31

Base + 5Δ
5
1.40
14
2.45

Coded factor : skim milk, : yeast extract.

bNatural factor : skim milk, : yeast extract.

cBased: 0 for coded factor and middle point for natural factor.

dMovement units (coded factor : 0.1/0.1, : 0.3/1.09) (natural factor: coded factor lower point).

eMaximum biomass of B. pseudocatenulatum G4 (log₁₀ cfu/mL) at 21 h of cultivation period. Values are mean of triplicates.

3.2. Optimization of Medium Component

The experimental responses for the optimization of skim milk and yeast extract are shown in Table 5. Skim milk concentration () was in the range of 4% (wt/vol) to 8% (wt/vol) with 6% (wt/vol) as a center point whereas, for yeast extract (), the range of 1% (wt/vol) to 2.8% (wt/vol) with 1.9% (wt/vol) as a center point was fixed further optimization. Center points with a coded (0,0) were repeated five times. The importance of medium components for cultivation process can be considered by their effect on biomass production () and β-galactosidase production () as well as the residual lactose () and free amino nitrogen () remaining in the culture. The responses (, , , and ) were fitted with quadratic polynomial model and subsequently produced the response surface as expressed in (3). Consider

where and represent coded independent factors of skim milk and yeast extract, respectively. Meanwhile, , , , , and are coefficients and is the random error.

Table 5: Experimental design and results using face-centered full factorial design (FCCD).

Run
Skim milk,
Yeast extract,
Responsesa

Coded value
Actual value
Coded value
Actual value
Max. biomass, (log₁₀ cfu/mL)
β-Galactosidase production (U/mL)
Lactose balance (g/L)
FAN balance (mg/L)

1
−1
4.0
−1
1.0
5.15
0.19
4.59
198.24

2
−1
4.0
0
1.9
5.18
3.03
4.51
209.58

3
−1
4.0
+1
2.8
5.22
3.13
4.40
213.24

4
0
6.0
−1
1.0
7.09
7.06
10.83
229.15

5
0
6.0
0
1.9
7.15
8.37
11.53
237.06

6
0
6.0
0
1.9
7.11
9.40
11.58
231.04

7
0
6.0
0
1.9
7.12
9.32
11.81
234.65

8
0
6.0
0
1.9
7.14
8.35
11.91
238.98

9
0
6.0
0
1.9
7.16
9.38
11.84
233.14

10
0
6.0
+1
2.8
7.13
8.41
12.90
359.76

11
+1
8.0
−1
1.0
7.26
7.02
22.41
287.53

12
+1
8.0
0
1.9
7.23
6.27
19.10
290.59

13
+1
8.0
+1
2.8
7.11
6.18
25.25
318.70

All factorial and axial points are means of duplicate.

3.2.1. Optimization of Biomass Production ()

By applying regression analysis on the experimental data, biomass production can be described by the second-order equation

Skim milk () was shown to be highly significant () whereas yeast extract () presented less influence towards biomass production. As a result, the quadratic model was further reduced with insignificant model term excluded. However, the addition of yeast extract (), as nitrogen source, in the milk medium must be considered to enhance growth of bifidobacteria [13]. Therefore, yeast extract is included in (4) to (5). Consider

Table 6 shows ANOVA result of quadratic model. Biomass production is significant as indicated by value (). The implies the sample variation of 99.93% for biomass production () and this indicates that only about 0.07% of total variation is not explained by the model, indicating good agreement between the experimental and predicted values for biomass production. The lack of fit measures the failure of the model to represent data in the experimental domain at points, which are not included in the regression [24]. The value of lack of fit for regression is insignificant (), suggesting that the model fitted well to the data in the experimental region. Moreover, second-order terms were found sufficient and higher order terms were not necessary.

Table 6: ANOVA and regression analysis for the response of biomass production (): optimization of skim milk (% wt/vol) and yeast extract (% wt/vol).

Source of variation
Sum of squire or coefficient estimate
DFa
MS or Std error
value
value

Modelb
8.95
4
2.24
2893.35
<0.001

Residual

8

Lack of fit

4

2.60
0.1890

Pure error

4

Total
8.95
12

Factorc

Intercept
7.14
1

1.01
1
0.011
7889.68
<0.001d

1
0.011
0.34
0.5732

−0.94
1
0.015
3667.73
<0.001d

−0.055
1
0.014
15.65
0.0042d

DF: degree of freedom.

b = 99.93%.

c, skim milk (% wt/vol), , yeast extract (% wt/vol).

dSignificant at .

The response surface model shown in Figure 1 indicated that skim milk concentration (% wt/vol) was the most important factor in the medium. Lactose in milk was able to supply readily utilizable carbon for growth of G4. Increased milk concentration, with small addition of yeast extract as growth enhancer, significantly increased biomass production. Lactose in milk induced β-galactosidase biosynthesis in Bifidobacterium strain, which break down lactose into glucose and galactose more efficiently during the cultivation process [13]. The response surface obtained is a stationary edge system, whereby the maximum value is a plane platform rather than point form. Thus, there is flexibility in choosing the appropriate optimum points [25].

Figure 1: Response surface methodology for biomass production (log₁₀ cfu/mL) from the adjusted quadratic mathematical model.

Based on the optimum point acquired from response surface methodology, the combination of 7.17% (wt/vol) skim milk, , and 1.02% (wt/vol) yeast extract, , was predicted to produce 7.44 log₁₀ cfu/mL biomass. This prediction was verified by a validation experiment. A maximum biomass of 7.35 log₁₀ cfu/mL was obtained from replication experiments. Even though the experimental value was lower than the predicted value, no statistical difference () was observed. In order to observe other growth activities of G4 in optimized medium, analysis of β-galactosidase production as well as determination of lactose and free amino nitrogen residual in the culture at 21 h of cultivation was performed. As shown in Table 7, the optimized medium (7.17% (wt/vol) skim milk and 1.02% (wt/vol) yeast extract) was able to support G4 growth up to 7 log₁₀ cfu/mL with substantially high β-galactosidase production at 21 h of cultivation. However, high amount of lactose and free amino nitrogen residual was still observed at the end of cultivation. To avoid wastage of nutrients remaining in the medium, further optimization step was carried out by considering residual nutrients remaining in the culture at the end of cultivation. This step was conducted and presented in the following section.

Table 7: Comparison experiments using optimized skim milk and yeast extract based on biomass production, β-galactosidase production, lactose, and free amino nitrogen residual at 21 h of cultivation period.

Skim milk (% wt/vol)*
Biomass production (log₁₀ cfu/mL)
β-Galactosidase production (U/mL)
Lactose (g/L)
Free amino nitrogen (mg/L)

0

4

7.17 (optimized concentrations)

10

12

TPY

Cultivation was conducted in different concentrations of skim milk with 1.02% yeast extract supplementation.

**Values in the same column with different letters were significantly different ().

3.2.2. β-Galactosidase Production (), Lactose (), and Free Amino Nitrogen () Presence as Responses

Further optimization of G4 medium cultivation was carried out by including β-galactosidase production, lactose, and free amino nitrogen balance at 21 h of cultivation as responses. The responses (, , and ) were fitted with the second-order polynomial equations

The statistical significance of the model was evaluated by the value of the analysis of variance (ANOVA). ANOVA statistics for the three responses, (β-galactosidase), (lactose residual), and (free amino nitrogen residual), at 21 h of the cultivation period are shown in Table 8. Quadratic models for and as well as linear model for are highly significant, as shown by low probability value (). The models fitted well to the experimental design as lack of fit for all three responses is insignificant (). Moreover, the coefficient of determinations () is close to 1 (: 0.97, : 0.99, and : 0.98), indicating high degree of correlation between observed and predicted values.

Table 8: ANOVA analysis for responses (β-galactosidase: U/mL), (lactose: g/L), and (FAN: mg/L).

Source
Sum squares
DF
Mean square
value
value

    Model
96.50
5
19.30
52.91
<0.0001

    Residual
2.55
7
0.36

    Lack of fit
1.33
3
0.44
1.46
0.3521

    Pure error
1.22
4
0.30

    Model
151.33
2
75.66
615.28
<0.0001

    Residual
1.23
10
0.12

    Lack of fit
1.05
6
0.17
3.83
0.1075

    Pure error
0.18
4
0.046

    Model
14609.65
5
2921.93
123.58
0.0001

    Residual
165.51
7
23.64

    Lack of fit
126.17
3
42.06
4.28
0.0971

    Pure error
39.34
4
9.84

The response surface plot of β-galactosidase production () is shown in Figure 2(a). High production of β-galactosidase was observed at the middle range of skim milk (6.00 to 7.00% (wt/vol)). On the other hand, reduced production of β-galactosidase was detected at high and low point of skim milk concentration. The use of high skim milk concentration (8.00% (wt/vol)) might contribute to the presence of high lactose concentration and this could suppress β-galactosidase production. Medium containing high lactose concentration (>5% [wt/vo]), may attributed to increase in the concentration of internally glucose. It has been reported that the increased of internal glucose would suppressed the biosynthesis of –galactosidase by tested organisms [26]. The presence of low lactose concentration in the medium using low skim milk concentration (4.00% (wt/vol)) may not induce β-galactosidase production.

(a)

(b)

(c)

Figure 2: Response surface plot of the effect of skim milk and yeast extract concentration (% wt/vol) for the responses of (a) β-galactosidase production (U/mL), (b) lactose residual (g/L), and (c) free amino nitrogen residual (mg/L) at 21 h of cultivation period.

A linear response surface plot of lactose residual () is shown in Figure 2(b). At high skim milk concentration (8.00% (wt/vol)) residual lactose concentration at 21 h of cultivation reached a maximum value with minor changes in yeast extract (1% to 2.8% (wt/vol)). In contrast, very low residual glucose concentration was detected at low skim milk concentration (4.00% (wt/vol)). This response indicated that skim milk concentration (4% to 8% (wt/vol)) was closely associated with residual lactose concentration. On the other hand, yeast extract ranging from 1% to 2.8% (wt/vol)) had less significant effect on the response. However, the presence of yeast extract was required to support β-galactosidase production (Figure 2(a)). Minimum residual lactose concentration at 21 h of cultivation was considered in this optimization method to avoid wastage of carbon source. From this point of view, low skim milk concentration is preferred to be used as cultivation medium for G4.

Figure 2(c) shows response surface plot representing free amino nitrogen residual concentration in the culture at 21 h of cultivation. The presence of nitrogen in medium was attributed to yeast extract and skim milk as well as dead cells and was believed to sustain bacterial growth. Similar to lactose, minimum residual concentration of nitrogen remaining in the culture after 21 h was considered in the optimization to avoid wastage of nitrogen source. The response surface generated based on the second-order coefficient shows that skim milk and yeast extract presented a quadratic effect (Figure 2(c)). The use of skim milk at high concentration with increasing yeast extract concentration (ranging from 1% to 2.8% (wt/vol)) significantly influenced the residual concentration of free amino nitrogen in the culture. On the other hand, reduced residual concentration of free amino nitrogen in the culture was observed at low skim milk concentration (4.00% (wt/vol)) with decreasing yeast extract concentration (ranging from 2.8% to 1% (wt/vol)). This result indicated that interactions between skim milk and yeast extract might have a stronger influence on residual concentration of free amino nitrogen remaining in the culture. Consequently, the use of low skim milk and yeast extract concentration significantly improved the nitrogen uptake efficiency with low residual concentration remaining in the culture at the end of cultivation.

3.3. Validation of Optimized Medium

The optimum concentrations for skim milk and yeast extract were 5.89% (wt/vol) and 2.31% (wt/vol), respectively. The responses criteria for the optimization process were maximum for biomass and β-galactosidase production. On the other hand, the minimum values were set for residual concentration of lactose and free amino nitrogen. An experiment was performed under the predicted optimal conditions in order to validate the optimized medium. The experimental values fitted well to the predicted results with no significant difference (). Therefore, this result was encountered in the process of validation of response surface methodology optimization.

4. Conclusion

Results from this study have demonstrated that the optimization of milk-based medium using response surface methodology (RSM) greatly improved G4 cultivation performance, in terms of final cell concentration and β-galactosidase as well as residual concentration of lactose and amino nitrogen remaining in the culture. The quadratic model (biomass, β-galactosidase production, and free amino nitrogen residual) and linear model (lactose residual) were found sufficient for the optimization of medium for G4. The optimal medium that consists of 5.89% (wt/vol) skim milk and 2.31% (wt/vol) yeast extract gave the final biomass count of 10⁷ cfu/mL, β-galactosidase activity of 8.74 U/mL, residual lactose concentration of 9.15 g/L, and residual free amino nitrogen concentration of 226.07 mg/mL. In addition, the use of this optimal medium gave comparable biomass count and β-galactosidase activity as those obtained in cultivation using high skim milk concentrations (10% to 13% (wt/vol)). However, residual lactose and free amino nitrogen concentration remaining in the culture at the end of cultivation were reduced by about two times lower as compared to those observed in cultivation using high skim milk concentration.

Conflict of Interests

The authors declare no conflict of interests.

Acknowledgments

This study was funded by the Ministry of Higher Education, Malaysia, under Project no. 600-RMI/FRGS 5/13 15/2013. The first author is also grateful to Universiti Teknologi MARA (UiTM), Malaysia, and Universiti Putra Malaysia (UPM) for the support throughout the conduct of this study.

References

FAO/WHO, “Health and nutritional properties of probiotics in food including powder milk with live lactic: report of a joint FAO/WHO expert consultation on acid bacteria expert consultation on evaluation of health and nutritional properties of probitics in food including powder milk with live lactic acid bacteria,” Tech. Rep., 2001.

Y. Ohashi and K. Ushida, “Health-beneficial effects of probiotics: its mode of action,” Animal Science Journal, vol. 80, no. 4, pp. 361–371, 2009.

T. Mitsouka, “Taxonomy and ecology of bifidobacteria,” Bifidobacteria Microflora, no. 3, pp. 11–28, 1984.

V. Scardovi, “Genus Bifidobacterium,” in Bergey’s Manual of Systemic Bacteriology, N. N. S. Mair, Ed., pp. 1418–1434, William and Wilkins, New York, NY, USA, 1986.

M. Shuhaimi, A. M. Ali, N. M. Saleh, and A. M. Yazid, “Differentiation of Bifidobacterium isolates from faeces of infant by RAPD,” Bioscience Microflora, vol. 20, pp. 89–94, 2001.

M. Shuhaimi, A. M. Ali, N. M. Saleh, and A. M. Yazid, “Utilisation of Enterobacterial Repetitive Intergenic Consensus (ERIC) sequence-based PCR to fingerprint the genomes of Bifidobacterium isolates and other probiotic bacteria,” Biotechnology Letters, vol. 23, no. 9, pp. 731–736, 2001.

M. Shuhaimi, A. M. Ali, M. S. Norihan, and A. M. Yazid, “Classification of Bifidoacteriumisolates from infant faeces using PCR-based and 16S rdnapartial sequence analysis methods,” Bioscience Microflora, vol. 21, pp. 155–161, 2002.

S. Wong, Y. K. Wendy, B. M. Kabeir, M. Shuhaimi, M. Rosfarizan, and A. M. Yazid, “Survival of Bifidobacterium pseudocatenulatum strain isolated from breast-fed infants to simulated gastric pH environment,” Malaysian Applied Biology, vol. 35, pp. 57–62, 2006.

L. R. Rodrigues, J. A. Teixeira, and R. Oliveira, “Low-cost fermentative medium for biosurfactant production by probiotic bacteria,” Biochemical Engineering Journal, vol. 32, no. 3, pp. 135–142, 2006.

J. Mark, “Successful probiotic bacteria,” in Handbook of Functional Dairy Products, C. Shortt and J. O’Brien, Eds., pp. 36–61, CRC Press, Boca Raton, Fla, USA, 2004.

V. M. Marshall and A. Y. Tamime, “Starter cultures employed in the manufacture of biofermented milks,” International Journal of Dairy Technology, vol. 50, no. 1, pp. 35–41, 1997.

H. M. Østlie, M. H. Helland, and J. A. Narvhus, “Growth and metabolism of selected strains of probiotic bacteria in milk,” International Journal of Food Microbiology, vol. 87, no. 1-2, pp. 17–27, 2003.

C. A. Hsu, R. C. Yu, and C. C. Chou, “Production of β-galactosidase by Bifidobacteria as influenced by various culture conditions,” International Journal of Food Microbiology, vol. 104, no. 2, pp. 197–206, 2005.

L. Avonts, E. V. Uytven, and L. D. Vuyst, “Cell growth and bacteriocin production of probiotic Lactobacillus strains in different media,” International Dairy Journal, vol. 14, no. 11, pp. 947–955, 2004.

A. Shihata and N. P. Shah, “Proteolytic profiles of yogurt and probiotic bacteria,” International Dairy Journal, vol. 10, no. 5-6, pp. 401–408, 2000.

R. I. Dave and N. P. Shah, “Evaluation of media for selective enumeration of Streptococcus thermophilus, Lactobacillus delbrueckii ssp. bulgaricus, Lactobacillus acidophilus, and Bifidobacteria,” Journal of Dairy Science, vol. 79, no. 9, pp. 1529–1536, 1996.

W. Stephanie, B. M. Kabier, M. Shuhaimi, M. Rosfarizan, and A. M. Yazid, “Growth optimization of a probiotic candidate, Bifidobacterium pseudocatenulatum G4, in milk medium using response surface methodology,” Biotechnology and Bioprocess Engineering, vol. 12, no. 2, pp. 106–113, 2007.

M.-J. Chen, K.-N. Chen, and C.-W. Lin, “Optimization on response surface models for the optimal manufacturing conditions of dairy tofu,” Journal of Food Engineering, vol. 68, no. 4, pp. 471–480, 2005.

F. Karacan, U. Ozden, and S. Karacan, “Optimization of manufacturing conditions for activated carbon from Turkish lignite by chemical activation using response surface methodology,” Applied Thermal Engineering, vol. 27, no. 7, pp. 1212–1218, 2007.

J.-W. Hou, R.-C. Yu, and C.-C. Chou, “Changes in some components of soymilk during fermentation with bifidobacteria,” Food Research International, vol. 33, no. 5, pp. 393–397, 2000.

M. N. Hwang and G. M. Ederer, “Rapid hippurate hydrolysis method for presumptive identification of group B streptococci,” Journal of Clinical Microbiology, vol. 1, no. 1, pp. 114–115, 1975.

Z. Nagy, T. Kiss, A. Szentirmai, and S. Biro, “Carbon source regulation of β-galactosidase biosynthesis in Penicillium chrysogenum,” Journal of Basic Microbiology, vol. 41, pp. 351–362, 2001.

B. K. Körbahti, N. Aktaş, and A. Tanyolaç, “Optimization of electrochemical treatment of industrial paint wastewater with response surface methodology,” Journal of Hazardous Materials, vol. 148, no. 1-2, pp. 83–90, 2007.

H. Gao and W.-Y. Gu, “Optimization of polysaccharide and ergosterol production from Agaricus brasiliensis by fermentation process,” Biochemical Engineering Journal, vol. 33, no. 3, pp. 202–210, 2007.

R. H. Myers and D. C. Montgomery, Response Surface Methodology: Process and Product Optimization Using Designed Experiments, John Wiley & Sons, New York, NY, USA, 2nd edition, 2002.

E. Demirhan, D. K. Apar, and B. Özbek, “Product inhibition of whey lactose hydrolysis,” Chemical Engineering Communications, vol. 195, no. 3, pp. 293–304, 2008.