Marine killer yeast Metschnikowia saccharicola active against pathogenic yeast in crab and an optimization of the toxin production

1 College of Food Science and Engineering, Ocean University of China, Qingdao 266003, China. 2 South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, National Research and Development Center for Aquatic Product Processing, Key Laboratory of Aquatic Product Processing, Ministry of Agriculture, Guangzhou 510300, China. 3 Beihang-goertec Microelectronics Institute, Beihang Qingdao Research Institute, Qingdao 266041, China. 4 Department of Quality and Regulatory Affairs, Air Liquide Medical Systems, 92182 Antony CEDEX, France.

Author(s) agree that this article remain permanently open access under the terms of the Creative Commons Attribution License 4.0 International License bicuspidata has been detected in crab Portunus trituberculatus cultured in 2001 in Zhoushan, Zhejiang Province, China (Xu et al., 2003).However, anti-yeast compounds such as nystatin, benzalkonium bromide, and extracts of gold thread root and garlic are toxic to the crab, so killer toxins produced by some yeast species against pathogenic yeasts may be a good alternative (Philliskink and Young, 1975).Studies have shown that killer yeasts and the killer toxins they produce have antimicrobial activity and can be used to control the growth of pathogenic yeasts in humans, animals and plants (Comitini et al., 2004;Magliani et al., 2008;Chi et al., 2010;Wang et al., 2013).
Yeast-derived killer toxins have been widely used to control harmful yeast growth in industrial fermentation systems, food production, animal breeding, agriculture, and medicine as antifungal drug formulations (Magliani et al., 2008;Wang et al., 2013).In the food industry, killer toxins inhibit the growth of wild yeast and prevent contamination, thus providing a good system environment and improving product quality (Schmitt and Breinig, 2002).In aquaculture, some diseases of marine animals, including P. trituberculatus can be controlled by certain marine killer toxin-producing yeast (Wang et al., 2007;Chi et al., 2010).
Optimizing the process of yeast fermentation is important to obtain maximum yield of killer toxins (Bandeira et al., 2006;Çorbacı and Uçar, 2017).However, the reaction system cannot be comprehensively analyzed using a traditional onevariable-at-a-time technique if more than one variable is present (Bezerra et al., 2008).Therefore, it is necessary to employ an appropriate statistical approach to evaluate the relationship among measurable variables.Response surface methodology (RSM) comprises mathematical and statistical techniques aimed at optimizing system response influenced by several independent variables (Bezerra et al., 2008), such as fermentation temperature and the amount of starter culture to obtain the best performance during fermentation (Yaakob et al., 2011).The objective of this study was to screen marine yeast species for the production of killer toxins against pathogenic yeast and determine the optimal conditions for killer toxin production in culture by investigating the interactions between process variables providing a basis for further research of disease prevention and control in marine aquaculture.

Yeast strains and media
M. bicuspidata WCY (collection number 2E00088 at the Marine Microorganisms Culture Collection of China) was confirmed pathogenic yeast in P. trituberculatus (Wang et al., 2007).Candida tropicalis and Candida albicans isolated from different marine environments were used as susceptible yeast strains.Yeast strains were grown in yeast extract peptone dextrose (YPD) medium containing 1.0% yeast extract, 2.0% glucose, 2.0% peptone, and 3% agar.The medium for killer toxin production was prepared aseptically by adding 2% NaCl and 15% glycerol to YPD, and medium pH was initially adjusted to 4.5 with 0.05 mol/L citric acidhydrogen phosphate disodium buffer.Killer toxin activity was assayed in YPD agar supplemented with 3 mg/mL methylene blue dissolved in ethanol (final concentration 0.003%) (Guo et al., 2013a).

Screening of marine killer yeast
Each yeast strain from the slants was grown in YPD liquid medium at 28°C, 140 rpm for 24 h; then, 2 ml of culture was centrifuged at 4°C for 5 min (5,000 × g) and washed three times with sterile water.Cell suspension was adjusted to 1 × 10 7 cells/ml, and 0.2 ml was inoculated into assay YPD agar plates seeded with susceptible pathogenic yeast strains.After 2 to 3 days of incubation at 28°C, a clear killing zone was observed around the colonies of killer yeast, which were then selected based on killing activity (determined as the ratio of the inhibition-zone diameter to colony diameter) (Wang et al., 2007;Peng, 2010).

Phylogenetic analysis
Yeast phylogenetic analysis was performed by the neighbor joining method using MEGA 5.1 (Tamura et al., 2011).Evolutionary distances were calculated using the p-distance model of MEGA 5.1, and bootstrap analysis was performed on 1,000 random resamplings.Reference sequences were retrieved from GenBank (accession numbers are indicated on the tree).

Measurement of killer toxin activity
Killer toxin activity was determined by a diffusion test using Oxford cups (6.0 × 10.0 mm) placed on the surface of assay agar plates (Santos et al., 2000) seeded with strain WCY.Killer toxin supernatant (250 μL) was added to each cup, and plates were incubated at 28°C for 2 days; the diameter of the inhibition zone was used as a measure of killer toxin activity (Peng, 2010).And each experiment was repeated 3 times.

Killer toxin production
Before designing the RSM experiment, selecting the inoculum  concentration, pH of the fermentation broth, fermentation time, and temperature as four independent variables for single-factor screening experiments were done.M. saccharicola was seeded at different inoculum concentrations (2 to 12%) in 500-mL Erlenmeyer flasks containing 100 mL of toxin production medium (pH 4.5) in a rotary bed shaker (140 rpm) at 28°C for 48 h.Based on the results, medium pH was varied from 3.5 to 6.5 using 0.05 M citric acidhydrogen phosphate disodium buffer, and yeast was cultured at different temperatures (20 to 35°C) for 12 to 72 h.

Experimental design and statistical analysis
According to the results of single factor tests, a three-level-fourfactor Box-Behnken Design (BBD) was applied to optimize fermentation conditions using the Design Expert (V 8.0.5) software.
The diameter of the inhibition zone (Y) was taken as a response of the design experiment.Full quadratic models (Bezerra et al., 2008) for killer toxin activity were established using the following mathematical equation: where Y is the response (diameter of inhibition zone, cm), β0 is the intercept, βi, βii, and βij are regression coefficients of the linear, quadratic, and interactive terms, respectively; and Xi and Xj represent coded independent variables.The fitted polynomial equation can be expressed as surface plots to visualize the relationships between responses and investigated parameters (Lwa et al., 2015).ANOVA was employed to evaluate the empirical mathematical model at 5% significance level and to test the significance of the difference between two or more sample mean differences.The statistical significance was considered by Pvalue，where the calculated P-value should be greater than the tabulated P-value to reject the null hypothesis and all the regression coefficients were 0 (Devore and Farnum, 2004).

Screening of marine killer yeast against pathogenic yeast strains
The yeast strain DD21-2 isolated from river bed sediments could secrete killer toxins onto culture medium and kill pathogenic yeast M. bicuspidata WCY, C. tropicalis Ct, and C. albicans YTS-03 (Figure 1).The killer activity of strain DD21-2 against M. bicuspidata WCY was higher than that against Ct and YTS-03, as evidenced by the ratio of diameter of inhibition zone to diameter of the colony been ≥1.5 (Figure 1); therefore, M. bicuspidata WCY was used as a sensitive strain in further experiments.

Identification of yeast strain DD21-2
Based on colony and cell morphology (Figure 2), fermentation and carbon source assimilation (Table 2),  and strain types of marine yeasts (Kurtzman and Fell, 2011), the yeast strain DD21-2 was related to M. saccharicola.Phylogenetic analysis based on ITS and the D1/D2 domain of 26S rDNA showed that many phylogenetically related yeast species were similar to the marine yeast strain and confirmed that DD21-2 was closely related to M. saccharicola (Figures 3 and 4).
Although the phylogenetic analysis of 26S rDNA showed that Metschnikowia is similar to Candida, they differ greatly in morphology and physiology (Kurtzman and Fell, 2011;Kaewwichian et al., 2012).Moreover, the nucleotide blast on NCBI also showed that only Metschnikowia saccharicola has 100% similarity to reference sequence.Phylogenetic analysis of ITS sequencing also fully proved this point (Figure 4).Therefore, the topology of the phylogenetic in Figures 3  and 4 confirms that strain DD21-2 was M. saccharicola (NCBI accession number KY849584 and MF115995).

Model fitting for RSM
The complete design matrix of the Box-Behnken design and the interaction effect of variables on the response are given in Table 3.A total of 29 experiments were conducted to determine the optimum fermentation conditions, and the data indicated that different levels of variables resulted in different response.
Multiple regression fitting was performed using Design Expert (V 8.0.5) (Shen et al., 2010), and full quadratic models were established using Equation 1. Through regression analysis, the relationship between inoculum concentration (X 1 ), pH of the fermentation broth (X 2 ), fermentation time (X 3 ), fermentation temperature (X 4 ), and the diameter of the inhibition zone (Y) was described by quadratic regression Equation 2: ANOVA showed that the P value of the model was 0.0228 (Table 4), indicating that regression was significant and that the experimental data and model were well fitted.
The P test showed that the quadratic effect or interaction of independent variables significantly influenced the diameter of the inhibition zone.The coefficient of determination (R

2
) for the quadratic model was 0.903, indicating that it was a well-fitting model that could adequately describe the relationships between the response (inhibition-zone diameter) and variables within the studied ranges.These results suggested that the empirical models for killer toxin activity provided good predictions for optimization.

Interaction effect of process variables on the diameter of the inhibition zone
The three-dimensional (3D) surfaces and the corresponding contour plots generated by the Design Expert software were used to study the effects of the parameters and their interactive effects on the killer toxin activity (Figure 6).According to the equation, the coefficients of the quadratic terms were negative, and the parabolic opening in 3D surfaces was downward, which indicated a maximum value point.Contour plots revealed significant interactive effects of the four factors on the diameter of the inhibition zone.The increment of the inhibition zone diameter with the increase of inoculum concentration and pH up to a critical point (Figure 6a) is in agreement with the fact-finding results of killer toxin' s optimal condition of Liu et al. (2012), and the ellipsoid contour plot of the inhibition-zone diameter indicated that the interaction between inoculum concentration and pH was considerable (Figure 6b).
The interaction effects of inoculum concentration and fermentation time on the inhibition zone diameter revealed that longer fermentation time produced a larger inhibition zone (Figure 6c).However, as the inoculum concentration and fermentation time increased, killer toxin activity decreased (Figure 6c), which could be    attributed to inactivation of the killer toxin during longer fermentation (Jia, 2012).In the interaction of the inoculum concentration and fermentation temperature, the diameter of the inhibition zone increased with the increase of the inoculum concentration and temperature up to the optimum level (Figure 6e).Similarly, the  interaction of pH and fermentation time significantly influenced the diameter of the inhibition zone (Figure 6g).The inhibition zone diameter was found to decrease linearly with the increase in the fermentation temperature above 27.5°C and pH above 5.5 (Figure 6i), consistent with the results obtained for Kloeckera apiculata KY-13c with optimal growth condition at pH 5.5 (Wu et al., 2015).Thus, high pH and fermentation temperature inhibited killer toxin activity.However, the DD21-2 strain could exert killer effects in a wider range of fermentation conditions compared with some other marine yeast, such as Williopsis saturnus WC91-2, which showed optimum toxin production at pH 3 to 3.5 and 16°C (Wang et al., 2012).Furthermore, fermentation time and temperature also had a significant effect on DD21-2 killer activity; the activity first increased and then decreased (Figure 6k), which was similar to that observed for Saccharomyces cerevisiae (Zhou et al., 2014).Thus, the interaction among the analysed factors significantly affected the inhibition-zone diameter, that is, killer toxin production.

Optimization of killer activity
Killer toxin production by M. saccharicola DD21-2 was optimized with RSM based on the four-factor, three-level BBD method.The results indicate that the four factors were significant in affecting killer activity, and the best fermentation conditions were determined for killer toxin production: Inoculum concentration 4.1% (v/v), pH 5.5, fermentation temperature 28°C, and fermentation time 2.7 days.These conditions were verified experimentally, and the diameters of the obtained triplicate inhibition zones were 1.45, 1.42 and 1.45 cm, which were very close to the predicted value of 1.41 cm.The error between the predicted and experimental values (2.84%) was within the 5% level of significance, indicating that the model is acceptable and the response surface optimization of fermentation conditions is effective.Thus, higher killer toxin yields can be achieved by the aforementioned methodology and the study suggests practical directions for further research in purification and characterization, and finally provides a basis for disease prevention and control in marine aquaculture.

Figure 3 .
Figure 3. Phylogenetic tree of yeast strain DD21-2 and 21 of its closest relatives (type strains) constructed using neighbor parsimony analysis of the D1/D2 domain of 26S rDNA.Bootstrap support values were calculated based on 1,000 pseudoreplications; values ≥50% are shown above the branches.Strain and sequence accession numbers are shown.

Figure 4 .
Figure 4. Phylogenetic tree of yeast strain DD21-2 and 23 of its closest relatives (type strains) based on the neighbor parsimony analysis of ITS sequences.Bootstrap support values were calculated based on 1,000 pseudoreplications; values ≥50% are shown above the branches.Strain and sequence accession numbers are shown.

Figure 5 .
Figure 5. Effects of inoculum concentration (a), pH (b), and fermentation time (c) and temperature (d) on the activity of the killer toxin.

Figure 6 .
Figure 6.The three-dimensional surface map and contour map of interactions between two variables on killer toxin production by yeast strains.

Table 1 .
Factors and levels of RSM test.

Table 3 .
Design and results of central composite tests.