Next Article in Journal
Ustiloxin G, a New Cyclopeptide Mycotoxin from Rice False Smut Balls
Next Article in Special Issue
Comparative Transcriptome Analysis of Penicillium citrinum Cultured with Different Carbon Sources Identifies Genes Involved in Citrinin Biosynthesis
Previous Article in Journal
Multi‐Mycotoxin Analysis in Durum Wheat Pasta by Liquid Chromatography Coupled to Quadrupole Orbitrap Mass Spectrometry
Previous Article in Special Issue
Biodegradation Mechanisms of Patulin in Candida guilliermondii: An iTRAQ-Based Proteomic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization for the Production of Deoxynivalenoland Zearalenone by Fusarium graminearum UsingResponse Surface Methodology

1
Institute of Crop Science, Chinese Academy of Agricultural Sciences, No.12 Zhongguancun South St., Haidian District, Beijing 100081, China
2
Laboratory of Quality and Safety Risk Assessment for Cereal Products (Beijing), Ministry of Agriculture, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Submission received: 22 November 2016 / Revised: 26 January 2017 / Accepted: 4 February 2017 / Published: 10 February 2017

Abstract

:
Fusarium mycotoxins deoxynivalenol (DON) and zearalenone (ZEN) are the most common contaminants in cereals worldwide, causing a wide range of adverse health effects on animals and humans. Many environmental factors can affect the production of these mycotoxins. Here, we have used response surface methodology (RSM) to optimize the Fusarium graminearum strain 29 culture conditions for maximal toxin production. Three factors, medium pH, incubation temperature and time, were optimized using a Box-Behnken design (BBD). The optimized conditions for DON production were pH 4.91 and an incubation temperature of 23.75 °C for 28 days, while maximal ZEN production required pH 9.00 and an incubation temperature of 15.05 °C for 28 days. The maximum levels of DON and ZEN production were 2811.17 ng/mL and 23789.70 ng/mL, respectively. Considering the total level of DON and ZEN, desirable yields of the mycotoxins were still obtained with medium pH of 6.86, an incubation temperature of 17.76 °C and a time of 28 days. The corresponding experimental values, from the validation experiments, fitted well with these predictions. This suggests that RSM could be used to optimize Fusarium mycotoxin levels, which are further purified for use as potential mycotoxin standards. Furthermore, it shows that acidic pH is a determinant for DON production, while an alkaline environment and lower temperature (approximately 15 °C) are favorable for ZEN accumulation. After extraction, separation and purification processes, the isolated mycotoxins were obtained through a simple purification process, with desirable yields, and acceptable purity. The mycotoxins could be used as potential analytical standards or chemical reagents for routine analysis.

1. Introduction

Fungi in the genus Fusarium are important pathogens of small-grain cereals, including wheat, maize, barley, and oats, especially in the temperate regions of the world [1,2]. Toxigenic fungi produce a variety of toxic metabolites that contaminate cereal grains and cereal-based food products, resulting in economic losses and potentially threatening the health of humans and animals [3,4,5]. Fusarium head blight (FHB), a common fungal disease of cereals, is caused by several Fusarium species worldwide, including F. graminearum, F. asiaticum, F. avenaceum, F. culmorum, F. tricinctum, F. langsethiae, F. sporotrichioides and F. poae [6,7,8]. Among these species, F. graminearum is ubiquitous and the most prevalent species in temperate regions, such as China, United States and other countries [9,10]. However, F. asiaticum is also a dominant species in China, especially in the southern rainy regions [11].
F. graminearum produces zearalenone (ZEN) and trichothecene mycotoxins, such as deoxynivalenol (DON), 3-acetyl-deoxynivalenol (3-ADON), 15-acetyl-deoxynivalenol (15-ADON), nivalenol (NIV) and 4-acetylnivalenol (4-ANIV) [12,13]. Based on the different type B trichothecene mycotoxins produced, F. graminearum can be classified into different chemotypes: chemotype I, which produces DON and/or its acetylated derivatives, and chemotype II, which produces NIV and/or 4-ANIV [14]. The DON chemotype can be further broken down into chemotype IA (producing DON and 3-ADON) and IB (producing DON and 15-ADON [15,16]. Additionally, some Fusarium isolates that produce both NIV and DON (NIV/DON chemotype) have been described as “unknown” chemotypes [17]. The most well-known and commonly found mycotoxins produced by different isolates of F. graminearum are DON and ZEN, although 3-ADON and 15-ADON can also be found frequently [18,19], accounting for 0%–16% of DON content [9].
Mycotoxin biosynthesis is a complex process that is regulated by genetic mechanisms, which can be affected by various environmental stimuli [20]. It has been reported that both the growth of fungi and their toxigenic potential can be affected by several environmental factors, including temperature, water activity, pH, and nutrient composition and availability [20,21,22,23]. Published data suggest that the same fungus can produce a different range of mycotoxins under different conditions, and even the biosynthesis of mycotoxins for toxin-producing strains is not necessary under certain conditions [24,25]. Merhej et al. [26] demonstrated that the pH regulatory factor Pac1 can regulate expression of the Tri genes, which are associated with trichothecene biosynthesis in F. graminearum and suggested that the production of trichothecene was induced only under acidic pH conditions.
A published report concerning the impact of environmental factors and fungicides on the growth and deoxinivalenol production of F. graminearum isolates from Argentinian wheat, showed that DON production in the presence of fungicides was influenced by complex interactions between water activity, temperature, fungicide concentration and incubation time [27]. However, there are limited data on the influence of environmental factors and their effects on mycotoxin production in Fusarium species. A better understanding of the role of culture conditions, including medium pH, incubation temperature and time, as well as their interactions on the production of DON and ZEN, is required, especially for F. graminearum strains. On the other hand, with increased concerns on the potential threat of Fusarium mycotoxin to human health, there has been a growing demand for mycotoxin standards to use in the relevant studies. In this way, maximizing the production of targeted mycotoxins may be possible by regulating the incubation conditions of F. graminearum. Mycotoxins with acceptable purity could be obtained by separation and purification processes and used as reference materials.
Response surface methodology (RSM), one of the most popular optimization techniques, has been widely used to estimate the relationships between a set of controllable experimental factors and obtained data [28,29]. Based on the model, independent variables can be controlled and optimal responses for the production of prime products can be obtained [30]. This methodology can greatly reduce the number of necessary experiments and provide a set of mathematical equations for theoretical process optimization [31]. Box-Behnken design (BBD), consisting of central and middle points on the edges of the cube circumscribed on the sphere, has been widely used for multivariate optimization of RSM. BBD produces a statistical model with fewer design points than central composite designs (CCD) [32].
This work was conducted to investigate the influence of cultivation conditions on the synthesis of mycotoxins and to maximize the production of DON and ZEN from F. graminearum using response surface methodology (RSM). After separation and purification, the isolated mycotoxin samples could have potential use as mycotoxin standards or chemical reagents for routine analysis.

2. Results and Discussion

The most well-known and commonly encountered mycotoxins of F. graminearum are ZEN, DON, and its acetylated derivatives (3-ADON, 15-ADON). Based on the reported literature, mycotoxin production can be affected by medium composition and environmental conditions, such as incubation temperature, time, and water activity [11,22,23,33]. In this study, the effect of incubation conditions on DON and ZEN production was evaluated, and the optimization of independent variables was performed by RSM.

2.1. Model Fitting and Statistical Analysis

Preliminary experiments showed that all of the variables examined in this study had an effect on DON and ZEN production. Therefore, the effects of three variables, including medium pH, incubation temperature and time on the responses (levels of DON and ZEN) were examined using Box-Behnken design (BBD). The complete design matrix, together with the response values are shown in Table 1. By applying multiple regression analysis, the predicted response Y for the yields of DON and ZEN could be obtained by the following second-order polynomial Equations (1) and (2):
YDON = −10314.19 + 2303.58X1 + 322.65X2 + 220.03X3 − 8.87X1X2 − 16.88X1X3 + 4.72X2X3 − 165.15X12 − 8.66X22 − 3.73X32
YZEN = 18106.60 − 2969.35X1 − 260.82X2 − 1090.25X3 − 93.59X1X2 + 56.38X1X3 − 10.59X2X3 + 521.77X12 + 19.47X22 + 29.44X32
where Y is the predicted response (levels of DON or ZEN), and X1, X2 and X3 are three independent variables, namely, media pH, incubation temperature, and time.
The results of the analysis of variance by Fisher’s F test, goodness-of-fit and the adequacy of the models are summarized in Table 2. For DON, a highly significant quadratic regression model was obtained with a high F value and a very low p-value (p = 0.0001 < 0.01), indicating that the combined effects of all the independent variables contributed significantly to maximizing the response [34]. The lack of fit (p = 0.1110 > 0.05) also suggested that the obtained data was a good fit with the model. The value of the determination coefficient (R2 = 0.9983) implied that the sample variation of 99.83% for DON production was attributable to the independent variables. The adjusted correlation coefficient (R2adj = 0.9962) and predicted determination coefficient (R2pre = 0.9793) for Equation (1) were also satisfactory to confirm the significance of the model. The higher of the adjusted correlation coefficients showed a better degree of correlation between the actual and predicted values [35].
For Equation (2), the p-value (p < 0.0001) and F-value (F = 2884.71) of the model suggested that it was significant. The value of the determination coefficient (R2 = 0.9997) indicated that approximately 0.03% of the total variance was not explained by the response. The adequacy of the model was also justified by the obtained adjusted determination coefficient (R2adj = 0.9994) and predicted determination coefficient (R2pre = 0.9958). The aforementioned data suggested that the mathematical model was reliable and could be used for ZEN production.
The model coefficients obtained by regression analysis for each variable are shown in Table 2. The significance of each variable and the strength of any interaction between each independent variable were tested using the p-value. Table 2 shows that the regression coefficients of all linear, quadratic and interaction terms were significant at the 1% level. The corresponding variables became more significant when the F-value was larger and the p-value was smaller [36]. From the values of the coefficients in the regression model, the order in which the independent variables affected DON production was in line with that of ZEN, which was pH (X1) > Time (X3) > Temperature (X2).

2.2. Response Surface Analysis

Using RSM, the effects of the independent variables (medium pH, incubation temperature and time) and their interactions on the yield of DON and ZEN could be graphically described using three-dimensional response surface plots and two-dimensional contour plots. Based on the regression model described in Equations (1) and (2), the responses were predicted and the optimum values for the production of DON and ZEN were determined [37,38].
The interaction effect of initial medium pH (X1) and incubation temperature (X2) on of DON and ZEN production is shown in Figure 1 (a1,a2) and Figure 2 (a1,a2). The pH value significantly affected the DON production, and was the most important factor which affected the response (Figure 1 (a1,a2)). This showed that DON levels tended to increase with increasing medium pH values within the range of 3.00–4.80, irrespective of the incubation temperature used, whereas DON concentration decreased with increasing pH in the ranged of 5.00–9.00. When pH was set to 9.00 (alkaline environment) irrespective of other factors, the model predicted that no DON would be produced by F. graminearum. This was in accordance with published reports, where Merhej et al. showed that the pH regulatory factor Pac1 regulated Tri gene expression and trichothecene production in F. graminearum, and the synthesis of trichothecene is induced only under acidic conditions [26]. Additionally, Gardiner et al. suggested that low extracellular pH both promoted and was required for DON production in F. graminearum [39]. Although low pH induces DON accumulation, it is not suitable for the growth and development of fungi [40], therefore the optimized pH was 4.70–5.10 in this study. For the effect of temperature, in the range of 15.0 °C–24.0 °C, the contour plot showed an increase in DON production with rising temperature, although DON levels tended to decrease when the temperature exceed 24 °C. Figure 2 (a1,a2) shows that ZEN production increased with increasing pH (4.80–9.00) and reached a maximum level at pH 9.00. However, the yield of ZEN seemed to be reduced with increasing temperature (15.5 °C–25 °C), meaning that temperature negatively affected the production of ZEN when pH and incubation time were fixed. Additionally, there were significant interaction effects between pH and temperature on ZEN production (Table 2, Figure 2 (a1,a2).
Figure 1 (b1,b2) shows that there was an increasing trend of DON level with increasing incubation time (14–28 days) when other factors were fixed. The same trend could be seen from three dimensional (3D) surface plot (Figure 2 (b1,b2)), which indicated that there was positive correlation between ZEN levels and incubation time. Thus, the maximum content of DON and ZEN was obtained when the fungi were cultivated for 28 days. As for DON, based on the interaction between pH and incubation time, the optimized plot was determined at pH of 4.50–5.20 and an incubation time of 28 days.
The interactions of incubation temperature and time on DON and ZEN levels are represented in Figure 1 (c1,c2) and Figure 2 (c1,c2). As shown in Figure 1c (c1,c2), both an increase in incubation temperature and time within the designed range resulted in increased DON yield, whereas no significant increase was found when the temperature reached 24 °C–25 °C. Thus, the maximum DON level was obtained when the fungi were cultivated for 28 days at 24 °C. On the contrary, temperature negatively affected ZEN production. The optimized condition for ZEN production was cultivation for 28 days at 15 °C. According to Garcia et al., the optimal temperature for mycotoxin production from F. graminearum was 15 °C for zearalenone and 20 °C for deoxynivalenol [22]. However, Llorens et al. showed that temperature significantly affected mycotoxin production and the optimal values were 28 °C, 20 °C and 15 °C for DON, NIV and 3-ADON production, respectively [23]. It was found that there were some differences in the optimal temperature for mycotoxin production among published reports. This may be attributed to the studied strains, which were from different localities and countries, and the different compositions of the culture medium. In this study, the optimal temperature for DON production was 23 °C–24 °C and for ZEN production was 15 °C–15.5 °C.

2.3. Optimization of Independent Variables and Validation of the Model

In order to obtain the desired response goal, optimization was performed by applying a Box-Behnken design of RSM to predict the optimum levels of the independent variables (medium pH, incubation temperature and time). Numerical and graphical optimization procedures were conducted to find the optimal process. Considering the experimental purpose, two different goals were determined: (1) optimize the conditions to achieve maximum levels of DON and ZEN, separately; and (2) find a combined condition to balance the production of DON and ZEN and obtain a satisfactory total yield for DON and ZEN. Based on the first target, the maximum DON level of 2811.17 ng/mL, under the condition pH 4.91, incubation temperature of 23.75 °C and 28 days of incubation, was predicted. However, the predicted maximum level of ZEN (23789.70 ng/mL) was produced at a pH of 9.0, an incubation temperature of 15.05 °C and incubation time of 28 days.
As the most favorable conditions for DON and ZEN production were inconsistent, we attempted to balance the two levels. The program combined the individual desirability into a single number, and then attempted to maximize this function [38]. The optimum working conditions and respective responses (DON and ZEN) were estimated, and the desirable solutions were found. Under these conditions, the model predicted DON and ZEN levels of 1718.74 ng/mL and 11797.9 ng/mL, respectively. The optimum conditions were pH of 6.86, incubation temperature of 17.75 °C, and an incubation time of 28 days.
To verify the accuracy of the model and equations, an optimization experiment was conducted under the following conditions: (i) pH value of 6.80, incubation temperature of 17.7 °C, and incubation time of 28 days; (ii) pH value of 4.90, incubation temperature of 23.80 °C, and incubation time of 28 days; and (iii) pH value of 9.00, incubation temperature of 15.00 °C, and incubation time of 28 days. The obtained DON and ZEN levels were 1658.05 ng/mL and 118103.40 ng/mL, respectively, under condition (i), while 2771.35 ng/mL of DON was obtained under condition (ii), and 23218.88 ng/mL of ZEN was obtained under condition (iii). These data showed that the actual values were all close to the predicted results. This suggested that the optimum conditions determined by RSM could be used to optimize the culture conditions for DON and ZEN production.
Several reports have shown that acid and alkaline environments can markedly affect the production of trichothecene mycotoxins, especially for DON, and suggested that DON could be only produced under acidic conditions from F. graminearum [23,36,37]. However, few publications are available concerning the acid-base conditions required for ZEN production. Our data showed that a higher pH value (weakly alkaline) was favorable for ZEN accumulation. This knowledge can provide further information to control and prevent mycotoxin contamination in the planting, harvesting and storage of cereals.

2.4. Purification and Purity of Mycotoxin Samples for Potential Analytical Standards

Based on our early experiments, we obtained acceptable recoveries (>80%) of mycotoxins, including DON and ZEN, following solvent extraction and hexane partitioning (data not shown). According to Wu et al., approximately 80% of fusaproliferin from cultures of Fusarium subglutinans was recovered after hexane partitioning. This demonstrated that hexane partitioning was an effective way to remove impurities from crude fusaproliferin extracts prior to HPLC [41].
After extraction of mycotoxins from potato dextrose broth (PDB) medium, DON and ZEN were isolated and purified by preparative HPLC. As shown in Figure 3, the targeted fractions were collected by an automatic fraction collector. The collected fractions were concentrated using a rotary evaporator and a freeze–dryer to a dry powder. Moisture content was then determined by UPLC-MS/MS and analytical HPLC based on comparing their retention times and MS data with the analytical standards and published data [42]. DON and ZEN purity was measured by areas of peak normalization method of HPLC [43], which showed that DON was >91% pure, while ZEN was >87% pure, suggesting that these isolated mycotoxins could potentially be used as analytical standards. In this study, a common nutrient broth of PDB with simple constituents was used. However, the yield of DON and ZEN might have been lower than some other F. graminearum isolates that were incubated with cereals, such as wheat and maize. Molto et al. showed that twenty 27 of 27 F. graminearum isolates produced deoxynivalenol (384–5745 µg/kg) and 13/27 produced zearalenone (200–35,045 µg/kg) when they were cultured with Argentinian maize [44]. Due to the simple purification steps, desirable yields, and acceptable purity of the isolated mycotoxins, their production is feasible, and the purified mycotoxin samples may be used as potential analytical standards or chemical reagents for routine analysis.

3. Conclusions

The proposed model equations illustrated both the quantitative effect and the interactions of the variables on mycotoxin production. By applying RSM, the optimized conditions for F. graminearum strain 29 were determined and desirable yields of the targeted mycotoxins were obtained, which could be further supplied as mycotoxin standards after purification and separation. Under the optimized conditions, maximal levels of DON (2811.17 ng/mL) and ZEN (23789.70 ng/mL) were obtained. A combined solution with a desirable total yield of both DON and ZEN was obtained, using a medium pH of 6.86, incubation temperature of 17.76 °C, and time of 28 days. The corresponding experimental values, from the validation experiments, fitted well with these predictions, suggesting that RSM could be applied to optimize Fusarium mycotoxin production. After extraction, separation and purification processes, the isolated mycotoxins were obtained with a desirable yield and acceptable purity. The mycotoxins could be used as potential analytical standards or chemical reagents for routine analysis.
This work also shows that the ambient pH value is a crucial factor for F. graminearum toxin production. Acidic pH is a determinant for DON production, while an alkaline environment is favorable for ZEN accumulation. As for incubation temperature, higher temperatures (23–25 °C) are suitable for DON accumulation, although a lower temperature (approximately 15 °C) favors a higher yield of ZEN.

4. Materials and Methods

4.1. Chemicals and Reagents

Mycotoxin standards, including DON and ZEN were purchased from Sigma-Aldrich (Sigma-Aldrich, St. Louis, MO, USA). Mycotoxin solid standards were dissolved in methanol (100 μg/mL) and stored at −20 °C in a sealed vial until use. All organic solvents, including methanol, acetonitrile and formic acid were of HPLC grade and were supplied by Fisher Scientific (Shanghai, China). Sodium hydroxide and concentrated hydrochloric acid were all analytical reagents and purchased from Beijing Chemical Reagent (Beijing, China). Pure water was obtained from a milli-Q system (Millipore, Billerica, MA, USA).

4.2. Fungal and Inoculum Preparation

Three fungal strains were kindly supplied by Professor Xiaoming Wang, Institute of Crop Science, Chinese Academy of Agricultural Sciences (CAAS, Beijing, China). Strains 001, 029 and 075 were obtained by single spore isolation from diseased corn spikes in China, which originated in China. The identification of the F. graminearum isolates was confirmed at the molecular level by a multiplex-PCR with primers Tri7F340/Tri7R965, 3551H/4056H, Tri3F971/Tri3R1679 and Tri3F1325/Tri3R1679 [45]. All strains belonged to the F. graminearum DON chemotype and can be further defined as the 15-ADON chemotype, which mainly produce ZEN, DON, and some of its acetyl-derivatives, namely, 15-ADON and 3-ADON. These isolates were incubated for 4–5 days at 25 °C on potato dextrose agar (PDA; potato extract, 4.0 g/L; dextrose, 20.0 g/L; agar, 15.0 g/L). Colonies were subcultured on new PDA plates and stored on PDA slants at 4 °C for further use. Following preliminary experiment to compare toxin production in the strains, strain 029 was chosen due to its high level of toxin production. Plugs (3 mm diameter) were collected with a cork borer from the margin of colonies grown on PDA medium in 9 cm Petri dishes for seven days at 25 °C. All the plugs were prepared and incubated in the fluid nutrient medium of potato dextrose broth (PDB).

4.3. Culture Conditions for Mycotoxin Production of F.graminearum Strain 29

Potato dextrose broth was used (PDB; potato extracts, 6.0 g/L; dextrose, 20.0 g/L) as the culture medium for fungal growth and mycotoxin production. The fungi was incubated in 200 mL of PDB for each run with different environmental factors. Medium pH (3–9), incubation temperature (15–25 °C), and time (14–28 days) were evaluated to find the most suitable culture condition for mycotoxin accumulation.

4.4. Experimental Design for Mycotoxin Production of F. graminearum Strain 29

To optimize the toxin-producing conditions for F. graminearum, the different culture conditions, including medium pH, incubation temperature and time, were evaluated using RSM. RSM was performed to assess the relationship between responses (mycotoxin production) and three independent variables, as well as to optimize the relevant conditions of variables aimed to predict the best value of responses [38].
A three-factor, three-level design was developed by Box-Behnken design (BBD; Design Expert software, Trial Version 8.0.6, Stat-Ease Inc., Minneapolis, MN, USA) to obtain a second-order polynomial model. Thereby, a suitable combination of variables could be determined to give the highest production of the targeted mycotoxins. Three variables, namely, medium pH (X1), incubation temperature (X2), and time (X3) and their appropriate ranges were determined based on single-factor experiments.
The coded and actual levels of the independent variables are shown in Table 3. A total of 17 experimental runs were carried out, and five replicates at the center point of the design were used for estimating the experimental error sum of squares. Triplicate analyses were performed at all design points in a random order.
The experimental data was analyzed using multiple regression and the second order polynomial model fitted for predicting optimal levels is expressed in Equation (3):
Y = β 0 + β 1 X 1 + β 2 X 2 + β 3 X 3 + β 12 X 1 X 2 + β 13 X 1 X 3 + β 23 X 2 X 3 + β 11 X 1 2 + β 22 X 2 2 + β 33 X 3 2
where Y is the predicted response; X1, X2, and X3 are independent variables; β0 is a constant coefficient; β1, β2, β3 are linear regression coefficients; β12, β13, β23 are interactive regression coefficients; and β11, β22, β33 are quadratic regression coefficients.

4.5. Mycotoxin Determination

UPLC–MS/MS was used for the analysis and determination of DON and ZEN. The methods for extraction, clean-up and determination of DON and ZEN from the cultures was based on previously published work [9] with slight modifications. Briefly, 5 mL culture filtrate was extracted with 20 mL of acetonitrile for 30 min using an automatic shaker. After centrifugation, 8 mL of the supernatant was passed through a MycoSep 226 Aflazon + multifunctional column (Romer Labs, Inc. Union, MO, USA) and 4 mL of the purified extract was evaporated to dryness under a stream of nitrogen. The residue was then dissolved in 1 mL of methanol–water (50/50, v/v), followed by filtering, and was subsequently used for analysis. Data acquisition and processing were performed using MassLynx v4.1 and Quanlynx (Waters, Milford, MA, USA).

4.6. Mycotoxins Extraction, Purification and Analysis for Potential Analyte Standards

4.6.1. Mycotoxin Extraction from Growth Medium

The extraction of DON and ZEN was carried out after 28 days of incubation in PDB based on the optimized culture conditions. The cultures were filtered with filter paper, and the mycelia were extracted with methanol, then the filtrate and extract were combined and concentrated under reduced pressure using a rotary evaporator at 65 °C. Subsequently, methanol/water (80/20, v/v) was added and the extracts were washed with hexane (60 mL) to defat. The aqueous and methanol layers were collected, concentated and filtered through a 0.45 µm nylon filter, and stored in brown glass vials at −20 °C until further purification and isolation [46,47].

4.6.2. Mycotoxin Seperation and Purification by Preparative Column HPLC

The mycotoxin-containing extracts dissolved in methanol were purified by preparative HPLC. The LC system was equipped with a LC-20APpump, SIL-10AP injector, FRC-10A fraction collector and a diode array detector (DAD) (Shimadzu, Kyoto, Japan), and a reverse-phase SHIM-pack PRC-ODS C18 column (20 mm × 250 mm, 5 μm, Shimadzu, Kyoto, Japan ) was used for mycotoxin seperation. The mobile phases were acetonitrile/methanol (50:50, v/v, solvent A) and H2O (solvent B) at a flow rate of 8 mL/min. The detector was set at λ1 = 220 nm for DON and its acetylated derivatives, and λ2 = 236 nm for ZEN. The linear gradient program started from 20% A for 8 min and increased to 80% A within 1 min, then held for 11 min. Finally, the initial composition of 20% A was re-established followed by equilibration for 5 min. LC retention times and UV absorbance profiles of the purified mycotoxins were compared to those of standard solutions. The targeted fractions were collected by a fraction collector. The LC-purified mycotoxin fractions were concentrated using a rotary evaporator at 45 °C, frozen overnight at −20 °C and then dried in a freeze–dryer for 48 h.

4.6.3. Identification and Purity of the Isolated Mycotoxin Samples

The identfication of the isloated mycotoxins was performed by UPLC-MS/MS with ESI (+) mode and analytical HPLC equipped with a DAD detector by comparing their retention times and MS data with the analytical standards and published data [9,42]. As for MS detection, the precursor ion (m/z) of 297.28, quantitative ion (Q) of 249.18 and qualitative ion (q) of 203.16 were used for DON analysis, while the precursor ion (m/z) of 319.05, quantitative ion (Q) of 283.03 and qualitative ion (q) of 187.36 were used for ZEN analysis. The purity of the isolated samples was determined by analytical HPLC using areas of peak normalization method [43].

4.7. Statistical Analysis

Stat-Ease software (Design-Expert 8.0.5 version, Stat-Ease Corporation, Minneapolis, MN, USA) was used for the regression analysis of the data, and to plot the response surface graphs. Analysis of variance (ANOVA) was performed and the values were considered significant when p < 0.05. The variability and accuracy of the model were determined according to the regression coefficient (R2) and lack of fit, respectively. The response surface and contour plots of the predicted responses of the model were used to assess interactions between the significant factors. Additionally, numerical optimization was carried out by performing three-dimensional response surface analysis of the independent and dependent variables.

Acknowledgments

We appreciate Xiaoming Wang of Institute of Crop Science, Chinese Academy of Agricultural Sciences (CAAS) for kindly supplying the F. graminearum strains used in this study. This work was financially supported by the National Key Program on Quality and Safety Risk Assessment for Agro-products and the Agricultural Science and Technology Program for Innovation Team on Quality and Safety Risk Assessment of Cereal Products, CAAS. We are grateful for the reviewers for their valuable suggestions on our paper.

Author Contributions

Li Wu conceived, performed the experiment and wrote the manuscript, Lijuan Qiu revised the language of manuscript, Huijie Zhang, Xuexu Hu and Juan Sun carried out the literature research, and Bujun Wang provides useful suggestions on the revision of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, S.; Hartman, G.L.; Domier, L.L.; Boykin, D. Quantification of Fusarium solani f. sp. glycines isolates in soybean roots by colony-forming unit assays and real-time quantitative PCR. Theor. Appl. Genet. 2008, 117, 343–352. [Google Scholar] [CrossRef] [PubMed]
  2. Youssef, S.A.; Maymon, M.; Zveibil, A.; Klein-Gueta, D.; Sztejnberg, A.; Shalaby, A.A.; Freeman, S. Epidemiological aspects of mango malformation disease caused by Fusarium mangiferae and source of infection in seedlings cultivated in orchards in Egypt. Plant Pathol. 2007, 56, 257–263. [Google Scholar] [CrossRef]
  3. Oliveira, P.M.; Zannini, E.; Arendt, E.K. Cereal fungal infection, mycotoxins, and lactic acid bacteria mediated bioprotection: From crop farming to cereal products. Food Microbiol. 2014, 37, 78–95. [Google Scholar] [CrossRef] [PubMed]
  4. Przybylska-Gornowicz, B.; Tarasiuk, M.; Lewczuk, B.; Prusik, M.; Ziółkowska, N.; Zielonka, Ł.; Gajęcki, M.; Gajęcka, M. The effects of low doses of two Fusarium toxins, zearalenone and deoxynivalenol, on the pig jejunum. A light and electron microscopic study. Toxins 2015, 7, 4684–4705. [Google Scholar] [CrossRef] [PubMed]
  5. Antonissen, G.; Martel, A.; Pasmans, F.; Ducatelle, R.; Verbrugghe, E.; Vandenbroucke, V.; Li, S.; Haesebrouck, F.; Van Immerseel, F.; Croubels, S. The impact of Fusarium mycotoxins on human and animal host susceptibility to infectious diseases. Toxins 2014, 6, 430–452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Yli-Mattila, T.; Mach, R.L.; Alekhina, I.A.; Bulat, S.A.; Koskinen, S.; Kullnig-Gradinger, C.M.; Kubicek, C.P.; Klemsdal, S.S. Phylogenetic relationship of Fusarium langsethiae to Fusarium poae and Fusarium sporotrichioides as inferred by IGS, ITS, β-tubulin sequences and UP-PCR hybridization analysis. Int. J. Food Microbiol. 2004, 95, 267–285. [Google Scholar] [CrossRef] [PubMed]
  7. Bottalico, A.; Perrone, G. Toxigenic Fusarium species and mycotoxins associated with head blight in small-grain cereals in Europe. In Mycotoxins in Plant Disease: Under the aegis of COST Action 835 ‘Agriculturally Important Toxigenic Fungi 1998-2003’, EU project (QLK 1-CT-1998-01380), and ISPP ‘Fusarium Committee’; Logrieco, A., Bailey, J.A., Corazza, L., Cooke, B.M., Eds.; Springer: Dordrecht, The Netherlands, 2002; pp. 611–624. [Google Scholar]
  8. Chilaka, C.; De Boevre, M.; Atanda, O.; De Saeger, S. Occurrence of Fusarium mycotoxins in cereal crops and processed products (Ogi) from Nigeria. Toxins 2016, 8, 342. [Google Scholar] [CrossRef] [PubMed]
  9. Wu, L.; Wang, B. Evaluation on levels and conversion profiles of DON, 3-ADON, and 15-ADON during bread making process. Food Chem. 2015, 185, 509–516. [Google Scholar] [CrossRef] [PubMed]
  10. Simsek, S.; Burgess, K.; Whitney, K.L.; Gu, Y.; Qian, S.Y. Analysis of deoxynivalenol and deoxynivalenol-3-glucoside in wheat. Food Control 2012, 26, 287–292. [Google Scholar] [CrossRef]
  11. Dong, F.; Qiu, J.; Xu, J.; Yu, M.; Wang, S.; Sun, Y.; Zhang, G.; Shi, J. Effect of environmental factors on Fusarium population and associated trichothecenes in wheat grain grown in Jiangsu province, China. Int. J. Food Microbiol. 2016, 230, 58–63. [Google Scholar] [CrossRef] [PubMed]
  12. Tralamazza, S.M.; Bemvenuti, R.H.; Zorzete, P.; de Souza Garcia, F.; Corrêa, B. Fungal diversity and natural occurrence of deoxynivalenol and zearalenone in freshly harvested wheat grains from Brazil. Food Chem. 2016, 196, 445–450. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, T.; Lee, S.-H.; Shin, J.Y.; Kim, H.-K.; Yun, S.-H.; Kim, H.-Y.; Lee, S.; Ryu, J.-G. Comparison of trichothecene biosynthetic gene expression between Fusarium graminearum and Fusarium asiaticum. Plant Pathol. J. 2014, 30, 33–42. [Google Scholar] [CrossRef] [PubMed]
  14. Sydenham, E.W.; Marasas, W.F.O.; Thiel, P.G.; Shephard, G.S.; Nieuwenhuis, J.J. Production of mycotoxins by selected Fusarium graminearum and F. crookwellense isolates. Food Addit. Contam. 1991, 8, 31–41. [Google Scholar] [CrossRef] [PubMed]
  15. Miller, J.D.; Greenhalgh, R.; Wang, Y.; Lu, M. Trichothecene chemotypes of three Fusarium species. Mycologia 1991, 83, 121–130. [Google Scholar] [CrossRef]
  16. Haratian, M.; Sharifnabi, B.; Alizadeh, A.; Safaie, N. PCR analysis of the Tri13 gene to determine the genetic potential of Fusarium graminearum isolates from Iran to produce nivalenol and deoxynivalenol. Mycopathologia 2008, 166, 109–116. [Google Scholar] [CrossRef] [PubMed]
  17. Quarta, A.; Mita, G.; Haidukowski, M.; Logrieco, A.; Mule, G.; Visconti, A. Multiplex PCR assay for the identification of nivalenol, 3- and 15-acetyl-deoxynivalenol chemotypes in Fusarium. FEMS Microbiol. Lett. 2006, 259, 7–13. [Google Scholar] [CrossRef] [PubMed]
  18. Edwards, S.G. Fusarium mycotoxin content of UK organic and conventional wheat. Food Addit. Contam. A 2009, 26, 496–506. [Google Scholar] [CrossRef] [PubMed]
  19. Ok, H.E.; Kim, H.J.; Cho, T.Y.; Oh, K.S.; Chun, H.S. Determination of deoxynivalenol in cereal-based foods and estimation of dietary exposure. J. Toxicol. Environ. Health A 2009, 72, 1424–1430. [Google Scholar] [CrossRef] [PubMed]
  20. Shwab, E.K.; Keller, N.P. Regulation of secondary metabolite production in filamentous ascomycetes. Mycol. Res. 2008, 112, 225–230. [Google Scholar] [CrossRef] [PubMed]
  21. Xue, H.; Bi, Y.; Tang, Y.; Zhao, Y.; Wang, Y. Effect of cultivars, Fusarium strains and storage temperature on trichothecenes production in inoculated potato tubers. Food Chem. 2014, 151, 236–242. [Google Scholar] [CrossRef] [PubMed]
  22. Garcia, D.; Barros, G.; Chulze, S.; Ramos, A.J.; Sanchis, V.; Marin, S. Impact of cycling temperatures on Fusarium verticillioides and Fusarium graminearum growth and mycotoxins production in soybean. J. Sci. Food Agric. 2012, 92, 2952–2959. [Google Scholar] [CrossRef] [PubMed]
  23. Llorens, A.; Mateo, R.; Hinojo, M.J.; Valle-Algarra, F.M.; Jimenez, M. Influence of environmental factors on the biosynthesis of type B trichothecenes by isolates of Fusarium spp. from Spanish crops. Int. J. Food Microbiol. 2004, 94, 43–54. [Google Scholar] [CrossRef] [PubMed]
  24. Kokkonen, M.; Ojala, L.; Parikka, P.; Jestoi, M. Mycotoxin production of selected Fusarium species at different culture conditions. Int. J. Food Microbiol. 2010, 143, 17–25. [Google Scholar] [CrossRef] [PubMed]
  25. Vogelgsang, S.; Sulyok, M.; Banziger, I.; Krska, R.; Schuhmacher, R.; Forrer, H.R. Effect of fungal strain and cereal substrate on in vitro mycotoxin production by Fusarium poae and Fusarium avenaceum. Food Addit. Contam. A 2008, 25, 745–757. [Google Scholar] [CrossRef] [PubMed]
  26. Merhej, J.; Richard-Forget, F.; Barreau, C. The pH regulatory factor Pac1 regulates Tri gene expression and trichothecene production in Fusarium graminearum. Fungal Genet. Biol. 2011, 48, 275–284. [Google Scholar] [CrossRef] [PubMed]
  27. Ramirez, M.L.; Chulze, S.; Magan, N. Impact of environmental factors and fungicides on growth and deoxinivalenol production by Fusarium graminearum isolates from Argentinian wheat. Crop Prot. 2004, 23, 117–125. [Google Scholar] [CrossRef]
  28. Gupta, S.; Shah, J.; Balasubramanian, B.A. Strategies for reducing colorectal cancer among blacks. JAMA Intern. Med. 2012, 172, 182–184. [Google Scholar] [CrossRef] [PubMed]
  29. Prasad, K.N.; Hassan, F.A.; Yang, B.; Kong, K.W.; Ramanan, R.N.; Azlan, A.; Ismail, A. Response surface optimisation for the extraction of phenolic compounds and antioxidant capacities of underutilised Mangifera pajang Kosterm. peels. Food Chem. 2011, 128, 1121–1127. [Google Scholar] [CrossRef] [Green Version]
  30. Nwabueze, T.U. Review article: Basic steps in adapting response surface methodology as mathematical modelling for bioprocess optimisation in the food systems. Int. J. Food Sci. Technol. 2010, 45, 1768–1776. [Google Scholar] [CrossRef]
  31. Box, G.E.P.; Wilson, K.B. On the Experimental Attainment of Optimum Conditions. In Breakthroughs in Statistics: Methodology and Distribution; Kotz, S., Johnson, N.L., Eds.; Springer: New York, NY, USA, 1992; pp. 270–310. [Google Scholar]
  32. Ferreira, S.L.C.; Bruns, R.E.; Ferreira, H.S.; Matos, G.D.; David, J.M.; Brandão, G.C.; da Silva, E.G.P.; Portugal, L.A.; dos Reis, P.S.; Souza, A.S.; et al. Box-Behnken design: An alternative for the optimization of analytical methods. Anal. Chim. Acta 2007, 597, 179–186. [Google Scholar] [CrossRef] [PubMed]
  33. West, J.S.; Holdgate, S.; Townsend, J.A.; Edwards, S.G.; Jennings, P.; Fitt, B.D.L. Impacts of changing climate and agronomic factors on Fusarium ear blight of wheat in the UK. Fungal Ecol. 2012, 5, 53–61. [Google Scholar] [CrossRef] [Green Version]
  34. Singh, P.; Shera, S.S.; Banik, J.; Banik, R.M. Optimization of cultural conditions using response surface methodology versus artificial neural network and modeling of l-glutaminase production by Bacillus cereus MTCC 1305. Bioresour. Technol. 2013, 137, 261–269. [Google Scholar] [CrossRef] [PubMed]
  35. Tian, Y.; Zeng, H.; Xu, Z.; Zheng, B.; Lin, Y.; Gan, C.; Lo, Y.M. Ultrasonic-assisted extraction and antioxidant activity of polysaccharides recovered from white button mushroom (Agaricus bisporus). Carbohyd. Polym. 2012, 88, 522–529. [Google Scholar] [CrossRef]
  36. Muralidhar, R.V.; Chirumamilla, R.R.; Ramachandran, V.N.; Marchant, R.; Nigam, P. Racemic resolution of RS-baclofen using lipase from Candida cylindracea. Mededelingen (Rijksuniv. Gent. Fak. Landbouwkund. Toegep. Biol. Wet.) 2001, 66, 227–232. [Google Scholar] [PubMed]
  37. Sahoo, C.; Gupta, A.K. Optimization of photocatalytic degradation of methyl blue using silver ion doped titanium dioxide by combination of experimental design and response surface approach. J. Hazard. Mater. 2012, 215–216, 302–310. [Google Scholar] [CrossRef] [PubMed]
  38. Kontogiannopoulos, K.N.; Patsios, S.I.; Karabelas, A.J. Tartaric acid recovery from winery lees using, cation exchange resin: Optimization by response surface methodology. Sep. Purif. Technol. 2016, 165, 32–41. [Google Scholar] [CrossRef]
  39. Gardiner, D.M.; Osborne, S.; Kazan, K.; Manners, J.M. Low pH regulates the production of deoxynivalenol by Fusarium graminearum. Microbiology 2009, 155, 3149–3156. [Google Scholar] [CrossRef] [PubMed]
  40. Merhej, J.; Boutigny, A.L.; Pinson-Gadais, L.; Richard-Forget, F.; Barreau, C. Acidic pH as a determinant of TRI gene expression and trichothecene B biosynthesis in Fusarium graminearum. Food Addit. Contam. A 2010, 27, 710–717. [Google Scholar] [CrossRef] [PubMed]
  41. Wu, X.; Leslie, J.F.; Thakur, R.A.; Smith, J.S. Purification of fusaproliferin from cultures of Fusarium subglutinans by preparative high-performance liquid chromatography. J. Agric. Food Chem. 2003, 51, 383–388. [Google Scholar] [CrossRef] [PubMed]
  42. Sun, J.; Li, W.; Zhang, Y.; Hu, X.; Wu, L.; Wang, B. QuEChERS purification combined with ultrahigh-performance liquid chromatography tandem mass spectrometry for simultaneous quantification of 25 mycotoxins in cereals. Toxins (Basel) 2016, 8, 375. [Google Scholar] [CrossRef] [PubMed]
  43. Papadoyannis, I.N.; Gika, H.G. Peak purity determination with a diode array detector. J. Liq. Chromatogr. Relat. Technol. 2004, 27, 1083–1092. [Google Scholar] [CrossRef]
  44. Molto, G.A.; Gonzalez, H.H.; Resnik, S.L.; Pereyra Gonzalez, A. Production of trichothecenes and zearalenone by isolates of Fusarium spp. from Argentinian maize. Food Addit. Contam. 1997, 14, 263–268. [Google Scholar] [CrossRef] [PubMed]
  45. Castañares, E.; Albuquerque, D.R.; Dinolfo, M.I.; Pinto, V.F.; Patriarca, A.; Stenglein, S.A. Trichothecene genotypes and production profiles of Fusarium graminearum isolates obtained from barley cultivated in Argentina. Int. J. Food Microbiol. 2014, 179, 57–63. [Google Scholar] [CrossRef] [PubMed]
  46. Nilanonta, C.; Isaka, M.; Kittakoop, P.; Trakulnaleamsai, S.; Tanticharoen, M.; Thebtaranonth, Y. Precursor-directed biosynthesis of beauvericin analogs by the insect pathogenic fungus Paecilomyces tenuipes BCC 1614. Tetrahedron 2002, 58, 3355–3360. [Google Scholar] [CrossRef]
  47. Meca, G.; Sospedra, I.; Soriano, J.M.; Ritieni, A.; Valero, M.A.; Manes, J. Isolation, purification and antibacterial effects of fusaproliferin produced by Fusarium subglutinans in submerged culture. Food Chem. Toxicol. 2009, 47, 2539–2543. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Three-dimensional response surface plots and corresponding contour plots of variables for DON production of F.graminearum strain 29. The interactions between medium pH and incubation temperature (a1,a2), pH and incubation time (b1,b2), and incubation temperature and time (c1,c2) are shown.
Figure 1. Three-dimensional response surface plots and corresponding contour plots of variables for DON production of F.graminearum strain 29. The interactions between medium pH and incubation temperature (a1,a2), pH and incubation time (b1,b2), and incubation temperature and time (c1,c2) are shown.
Toxins 09 00057 g001aToxins 09 00057 g001bToxins 09 00057 g001c
Figure 2. Three-dimensional response surface plots and corresponding contour plots of variables for ZEN production of F.graminearum strain 29. The interactions between medium pH and incubation temperature (a1,a2), pH and incubation time (b1,b2), and incubation temperature and time (c1,c2) are shown.
Figure 2. Three-dimensional response surface plots and corresponding contour plots of variables for ZEN production of F.graminearum strain 29. The interactions between medium pH and incubation temperature (a1,a2), pH and incubation time (b1,b2), and incubation temperature and time (c1,c2) are shown.
Toxins 09 00057 g002aToxins 09 00057 g002bToxins 09 00057 g002c
Figure 3. UV dual-wavelength spectrophotometry chromatogram of mycotoxins analyzed by preparative HPLC.
Figure 3. UV dual-wavelength spectrophotometry chromatogram of mycotoxins analyzed by preparative HPLC.
Toxins 09 00057 g003
Table 1. Experimental design used in the response surface methodology studies based on three independent variables and the observed responses of F. graminearum strain 29.
Table 1. Experimental design used in the response surface methodology studies based on three independent variables and the observed responses of F. graminearum strain 29.
No.X1X2X3Mycotoxin Levles (ng/mL)
PHTemperature (°C)Time (day)DONZEA
1−1 (3) −1 (15) 0 (21) 961.02 ± 35.78419.97 ± 30.25
21 (9) −1 (15) 0 (21) <LOD19084.20 ± 1102.56
3−1 (3) 1 (25) 0 (21) 1493.05 ± 91.23609.66 ± 41.35
41 (9) 1 (25) 0 (21) <LOD13457.90 ± 978.25
5−1 (3) 0 (20) −1 (14) 585.59 ± 39.75818.34 ± 54.25
61 (9) 0 (20) −1 (14) <LOD13941.4 ± 886.14
7−1 (3) 0 (20) 1 (28) 2003.54 ± 123.472288.42 ± 231.52
81 (9) 0 (20) 1 (28) <LOD20047.30 ± 1203.59
90 (6) −1 (15) −1 (14) 1559.09 ± 101.274118.31 ± 334.59
100 (6) 1 (25) −1 (14) 1405.5 ± 89.582259.09 ± 189.67
110 (6) −1 (15) 1 (28) 2098.03 ± 128.718661.47 ± 678.15
120 (6) 1 (25) 1 (28) 2605.26 ± 134.555319.75 ± 408.14
130 (6) 0 (20) 0 (21) 2276.33 ± 121.303198.34 ± 278.27
140 (6) 0 (20) 0 (21) 2313.33 ± 108.243133.31 ± 256.77
150 (6) 0 (20) 0 (21) 2332.08 ± 107.573126.45 ± 281.73
160 (6) 0 (20) 0 (21) 2283.66 ± 119.353156.87 ± 312.25
170 (6) 0 (20) 0 (21) 2376.34 ± 123.833185.81 ± 243.51
Table 2. ANOVA for response surface quadratic models for targeted mycotoxin production.
Table 2. ANOVA for response surface quadratic models for targeted mycotoxin production.
SourceSum of SquaresdfMean SquaresF-Valuep-Value
DON
Model1.520 × 107 91.689 × 106461.940.0001
X13.179 × 10613.179 × 106869.41<0.0001
X298051.42198051.4226.820.0013
X31.246 × 10611.246 × 106340.770.0001
X1 X270763.98170763.9819.350.0032
X1 X35.030 × 10515.030 × 105137.57<0.0001
X2 X31.092 × 10511.092 × 10529.860.0009
X129.302 × 10619.302 × 1062544.13<0.0001
X221.974 × 10511.974 × 10553.990.0002
X321.408 × 10511.408 × 10538.510.0005
Residual25592.7973656.11
Lack of Fit19068.0036356.003.900.1110
Pure Error6524.7941631.20
Cor Total1.523 × 10716
ZEN
Model6.534 × 10897.260 × 1072884.71<0.0001
X14.882 × 10814.882 × 10819398.730.0001
X21.414 × 10711.414 × 107562.050.0001
X32.918 × 10712.918 × 1071159.650.0001
X1 X27.885 × 10617.885 × 106313.310.0001
X1 X35.607 × 10615.607 × 106222.800.0001
X2 X35.495 × 10515.495 × 10521.830.0023
X129.285 × 10719.285 × 1073689.540.0001
X229.977 × 10519.977 × 10539.650.0004
X328.764 × 10618.764 × 106348.240.0001
Residual1.762 × 105725166.28
Lack of Fit175500357393.405.760.0710
Pure Error3983.754995.94
Cor Total6.536 × 10816
Table 3. Experimental range and levels of the independent variables.
Table 3. Experimental range and levels of the independent variables.
VariablesSymbolsRange and Levels
Low (−1)Medium (0)High (+1)
Temperature (°C)X1152025
pHX2369
Time (d)X3142128

Share and Cite

MDPI and ACS Style

Wu, L.; Qiu, L.; Zhang, H.; Sun, J.; Hu, X.; Wang, B. Optimization for the Production of Deoxynivalenoland Zearalenone by Fusarium graminearum UsingResponse Surface Methodology. Toxins 2017, 9, 57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins9020057

AMA Style

Wu L, Qiu L, Zhang H, Sun J, Hu X, Wang B. Optimization for the Production of Deoxynivalenoland Zearalenone by Fusarium graminearum UsingResponse Surface Methodology. Toxins. 2017; 9(2):57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins9020057

Chicago/Turabian Style

Wu, Li, Lijuan Qiu, Huijie Zhang, Juan Sun, Xuexu Hu, and Bujun Wang. 2017. "Optimization for the Production of Deoxynivalenoland Zearalenone by Fusarium graminearum UsingResponse Surface Methodology" Toxins 9, no. 2: 57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins9020057

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

Article Metrics

Back to TopTop