Optimization of Culture Conditions for Enrichment of Saccharomyces cerevisiae with Dl-α-Tocopherol by Response Surface Methodology

authors:

avatar Morteza Mohajeri Amiri a , avatar Mohammad Reza Fazeli b , avatar Mohsen Amini c , avatar Nasim Hayati Roodbari a , avatar Nasrin Samadi b , *

Department of Biology, Science and Research Branch, Islamic Azad University, Tehran, Iran.
Department of Drug and Food Control, Faculty of Pharmacy and Pharmaceutical Quality Assurance Research Center, Tehran University of Medical Sciences, Tehran, Iran.
Department of Medicinal Chemistry, Faculty of Pharmacy and Drug Design and Development Research Center, Tehran University of Medical Sciences, Tehran, Iran.

how to cite: Mohajeri Amiri M, Fazeli M R, Amini M, Hayati Roodbari N, Samadi N. Optimization of Culture Conditions for Enrichment of Saccharomyces cerevisiae with Dl-α-Tocopherol by Response Surface Methodology. Iran J Pharm Res. 2017;16(4):e125060. https://doi.org/10.22037/ijpr.2017.2146.

Abstract

Designing enriched probiotic supplements may have some advantages including protection of probiotic microorganism from oxidative destruction, improving enzyme activity of the gastrointestinal tract, and probably increasing half-life of micronutrient. In this study Saccharomyces cerevisiae enriched with dl-α-tocopherol was produced as an accumulator and transporter of a lipid soluble vitamin for the first time. By using one variable at the time screening studies, three independent variables were selected. Optimization of the level of dl-α-tocopherol entrapment in S. cerevisiae cells was performed by using Box-Behnken design via design expert software. A modified quadratic polynomial model appropriately fit the data. The convex shape of three-dimensional plots reveal that we could calculate the optimal point of the response in the range of parameters. The optimum points of independent parameters to maximize the response were dl-α-tocopherol initial concentration of 7625.82 µg/mL, sucrose concentration of 6.86 % w/v, and shaking speed of 137.70 rpm. Under these conditions, the maximum level of dl-α-tocopherol in dry cell weight of S. cerevisiae was 5.74 µg/g. The resemblance between the R-squared and adjusted R-squared and acceptable value of C.V% revealed acceptability and accuracy of the model.

Introduction

Yeast enrichment technique that is used in food microbial studies aims at increasing the content of essential micronutrients such as vitamins and minerals in food supplements to compensate micronutrients deficiencies and to provide a public health benefit (1). Saccharomyces cerevisiae, known as baker’s yeast, is a well-known probiotic which can grow in simple environmental conditions and minimal media (2-5). Because of its safety and unique fermentation properties, yeast cells are commonly used in the human and animal foods and biotechnology industries (2, 6). European food safety authority (EFSA) categorized S. cerevisiaeas as a microorganism that has QPS (Qualified Presumption of Safety) status (7).

α-tocopherol (vitamin E) is a member of eight fat soluble compounds with a chromanol ring and tocopherol chain (8). It may protect cell membranes, low density lipoprotein, and human plasma from oxidative stress by scavenging free-radicals (9). Additionally, α-tocopherol probably has beneficial effects on pregnancy (10), age-related macular degeneration (11), cancer (12), cataracts (13), glaucoma (14), Alzheimer’s, and Parkinson’s diseases (15). Furthermore α-tocopherol is important in human and animals as immune modulator (16, 17). Fully synthetic α-tocopherol "dl-α-tocopherol" and its acetate derivative, which is more stable, are usually used in food supplements (18).

There are a number of researches which have been used S. cerevisiae as a vehicle for some nutrients. Some studies reported uptake of essential trace metals such as zinc, iron, copper, and manganese by the yeast S. cerevisiae (1, 19 and 20). In addition to increasing the yeast biomass, the inorganic toxic metal elements with low bioavailability will be changed to safer and more bioactive types with the enhanced nutritional properties (21). Arshad et al. indicated that antioxidant activity of poultry meat was depended on the value of α-tocopherol concentration added to the diet (22). Furthermore α-tocopherol changed enzyme activity of the gastrointestinal epithelium in chicken which resulted in higher absorption of nutrients in digestive system (16). These studies indicated that using combination of α-tocopherol and S. cerevisiae may have some advantages including designing enriched probiotic supplements, protection of yeast as a probiotic microorganism from oxidative destruction, improving enzyme activity of the gastrointestinal tract and probably increasing half-life of α-tocopherol by protecting against destructor agents in digestive system.

Experimental design and statistical analysis are important tools for optimization (23). These techniques can be used for assessment of multivariable system with minimal experimental trials (23, 24). One of the most popular approach in experimental design is response surface methodology (RSM). RSM is a group of techniques for designing experiments, evaluating the interaction of variables, fitting appropriate mathematical models, and optimization of arbitrary responses (24).

In this study S. cerevisiae was used as a probiotic model for dl-α-tocopherol accumulation and transport for the first time. Three independent variables were selected for determination of optimal conditions for enrichment of S. cerevisiae biomass with dl-α-tocopherol by using statistical analysis and response surface methodology.

Experimental

Microorganism, culture condition and chemicals

S. cerevisiae ATCC 9763 was maintained on Sabouraud Dextrose agar (SDA, Merck Co. Germany) and stored at 4 °C. The dl-α-tocopherol was obtained from DSM Co. Basel, Switzerland. All other chemicals were obtained from Sigma-Aldrich Co. USA.

Fermentation conditions

Shake flask experiments by using Sabouraud dextrose broth (SDB) were conducted to optimize cultivation conditions. The 500-mL Erlenmeyer flasks containing 200 mL of SDB medium were inoculated with 5% (v/v) of a 24-h preculture to give an initial yeast cell concentration of approximately equivalent to 0.1 g/L. Vitamin E (dl-α-tocopheryl acetate acetate) was added to medium at an initial concentrations of 2500-20000 µg/mL in preliminary studies. The cultures were incubated on a rotary shaker (J Labtech, DaihanLabtech Co.) at 100-200 rpm in dark condition. The samples were withdrawn at regular intervals and analyzed for determination of dry cell weight (g/L), and dl-α-tocopherol uptake (µg/g).

Dry cell weight measurement

For measurement of dry cell weight, 15 mL of SDB culture at each time interval was centrifuged at 8000×g for 10 min and washed twice with deionized water. The biomass was dried at 50 °C under reduced pressure and weighed for determination of dry cell weight.

Determination of dl-α-tocopherol uptake

Extraction procedure

Extraction of dl-α-tocopheryl acetate from the yeast cells was performed according to the method proposed by Popović et al. (25). The culture was centrifuged for 10 min at 1500×g and washed twice with phosphate-buffered saline (pH 7.4). Then, approximately 0.5 g of biomass was added to 30 mL of methanol kept at room temperature and dark condition for 16 h. The suspension was sonicated in an ultrasonic bath (Q Sonica, Q700) for 15-20 min. The solvent was evaporated at 40 °C under nitrogen steam. Then, 10 mL of n-hexan was added to the upper phase and sonicated for 15 min. Subsequently, n-hexan was evaporated. The residue was dissolved in 1 mL of n-hexan and filtered (0.22 μm pore size) prior to HPLC analysis.

HPLC conditions

Analysis was performed with Chemstation Software (Agilent Technologies). Chromatographic separation was performed on an Agilent column (Zorbax eclipse plus-C18 column), 10 cm, 2.1 mm, and 1.8 μm. Column oven temperature ranged from 25 to 40 °C. Mobile phase consisted of 100% n-hexan with a flow-rate of 0.3 mL/min and a pressure of 105 bars. The injection volume was 5 μL. Total run time was 20 min and experiment wavelength was 285 nm according to Popović et al. (25). The amount of uptaken dl-α-tocopherol was expressed as µg per gram of dry S. cerevisiae cells.

Experimental Design and Statistical Analysis

Screening Study

One-variable-at-a-time experiments were carried out to find the range of the parameters used for screening study. The levels of pH (4.5 to 7.5), temperature (25-35 °C), carbon source (glucose, sucrose, molasses, glycerol, maltose, starch, and fructose), nitrogen source (ammonium sulfate, meat extract, yeast extract, peptone and tryptone), primary vitamin E (dl-α-tocopheryl acetate) concentration, shaking speed, and incubation time were used in the optimization studies.

Our preliminary studies showed that among the selected factors, initial dl-α-tocopherol concentration, shaking speed and addition of sucrose to SDB medium play important roles in dl-α-tocopherol enrichment inside the yeast cells.

Optimization Design

Optimum level of each parameter was determined by a statistical method termed response surface methodology (RSM) by using Design-Expert (version 10.0.1, Stat-Ease, Inc., Minneapolis, MN, USA). Optimization of the level of dl-α-tocopherol entrapment in S. cerevisiae cells was performed by Box-Behnken (26) as a RSM method. Variables that were considered significant in screening study, including initial dl-α-tocopherol concentration (A), shaking speed (B), and concentration of sucrose added to SDB medium (C), were defined at three levels (low, basal and high), coded as (-1, 0 and +1) (Table 1). In Box-Behnken technique designed by Design-Expert® software, a set of seventeen experiments with five central point replicates was conducted.

The equation of the system was described by the modified quadratic polynomial model:

Y = α0 +∑αiXi + ∑αiiXi2 + ∑ αijXiXj + ε

where Y is the response, α0 is a constant, αi and αii reveal linear and quadratic effects respectively, αij is the quadratic effects of the interactions, Xi and Xj are coded values of the parameters and ε is the error of model. In this study Y reveals the expected dl-α-tocopherol amount in dry weight of yeast biomass.

The significance of the model on the response was checked by the analysis of variance (ANOVA) and was determined by a p-value below 0.05. Goodness of fit for modified quadratic polynomial model equation was determined by the multiple correlation coefficient (R2) and adjusted R2. In this study, response surface plot indicating the effect of coded variable (dl-α-tocopherol concentration, shaking speed and sucrose) interactions and relations on dl-α-tocopherol amount in dry cell weight were main responses.

Plot of experimental versus predicted response values and plot of studentized residual versus predicted response values were checked.

Results and Discussion

Box-Behnken Design and experimental responses

Our preliminary study showed that among selected factors, the dl-α-tocopherol concentration, shaking speed, and addition of sucrose to fermentation medium play important roles.

The data obtained from experiments were analyzed statistically to determine which interaction was significant. According to Box-Behnken design, each parameter was changed at three levels while the other parameters were kept constant. The results of Box- Behnken design are shown in Table 2.

Response Surface Methodology Analysis

The information obtained from the design showed that the range of dl-α-tocopherol in dry weight was 1.48-5.64 (µg/g). To select the appropriate model to fit the data, the analysis of variance by calculating F-value was used (Table 3). According to the results, we found that a modified quadratic polynomial equation appropriately fit the data. The lack of fit F-value of 2.93 suggested that the lack of fit is not significant that presents acceptability of the model.

The R-squared (R2) and adjusted R-squared (R2 adj) coefficients of this model are 0.993 and 0.986, respectively, indicating that the noise of the system is less than one percent. Furthermore, the resemblance between the R2 and adjusted R2 reveals the acceptability of the model to anticipate the results in the optimization. The predicted R2 (0.933) shows a very appropriate harmony between the value predicted by the model and the actual data (Figure 1).

Furthermore, the relative standard deviation (RSD) or C.V.% is 4.42 that indicates the accuracy and repeatability of the model. Adeq precision ratio of 32.833 indicates an adequate signal. Adeq precision measures the signal to noise ratio and a ratio of greater than 4 is desirable. It means that the model can be used to navigate the design space.

Absence of trends in the plot of studentized residual versus the values predicted by the model shows that the variances in the data are acceptable and no outliers are present in the experiments (Figure 2).

Table 1

Factors used in the Box-Behnken design and their levels

ParametersUnitSymbolLevel -1Level 0Level +1
α-tocopherol concentrationµg/mLA2500.06250.010000.0
Shaking speedrpmB100.0150.0200.0
Sucrose% w/vC2.05.08.0
Table 2

The Box-Behnken experiments and the dl-α-tocopherol amounts in dry cell weight of yeast.

Response: dl-α-tocopherol per dry cell weight (µg/g)Factor 3C: sucrose (% w/v)Factor 2B: shaking speed (rpm)Factor 1A: Initial dl-α-tocopherol con. (µg/mL)BlockRunStd
2.26785.00100.002500.00Block121
4.18875.00100.0010000.00Block112
1.48055.00200.002500.00Block193
2.05535.00200.0010000.00Block1174
2.46032.00150.002500.00Block185
3.46572.00150.0010000.00Block1166
2.89038.00150.002500.00Block147
5.20198.00150.0010000.00Block1128
3.16822.00100.006250.00Block159
2.08742.00200.006250.00Block11010
4.58008.00100.006250.00Block1611
3.34728.00200.006250.00Block1312
5.50765.00150.006250.00Block11313
5.46265.00150.006250.00Block11514
5.64175.00150.006250.00Block11415
5.37445.00150.006250.00Block1716
5.33785.00150.006250.00Block11117
Table 3

ANOVA for response surface modified quadratic model analysis of variance.

SourceSum of SquaresdfMean SquareF-Valuep-value Prob > F
Model33.4684.18148.35< 0.0001 significant
A-Initialdl-α-tocopherol concentration4.2214.22149.82< 0.0001
B-shaking speed3.4213.42121.49< 0.0001
C-sucrose2.9312.93103.78< 0.0001
AB0.4510.4516.070.0039
AC0.4310.4315.130.0046
A28.0118.01284.01< 0.0001
B210.62110.62376.55< 0.0001
C21.4211.4250.470.0001
Residual0.2380.028
Lack of Fit0.1740.0422.930.1615 not significant
Pure Error0.05740.014
Cor Total33.6816
R20.9933
R2 adj0.9866
C.V.%4.42
Adeq Precision32.833
Table 4

Validity of optimization design

ResponseSoftware predictionValue (µg/g)Validity experiment value (µg/g)95% CI low95% CI high
dl-α-tocopherol per dry cell weight5.875.745.696.04
Plot of experimental versus predicted values of dl-α-tocopherol amounts per dry cell weight of S. cerevisiae
Plot of studentized residual versus predicted values of α-tocopherol amounts per dry cell weight of S. cerevisia
Response surface plot indicating the effect of shaking speed and initial dl-α-tocopherol concentration interaction on dl-α-tocopherol amount per dry cell weight of S. cerevisia
Response surface plot indicating the effect of sucrose (% w/v) added to culture medium and initial dl-α-tocopherol concentration interaction on dl-α-tocopherol amount per dry cell weight of S. cerevisia
Response surface plot indicating the effect of sucrose (% w/v) added to culture medium and shaking speed (rpm) interaction on dl-α-tocopherol amount per dry cell weight of S. cerevisia

Eventually, in this study, the modified quadratic polynomial equation is an appropriate model to explain the level of dl-α-tocopherol amounts entrapped in S. cerevisiae cells. The regression analysis report for modified quadratic model showed that linear, squared, and interaction coefficients are significant (p < 0.05) except one linear interaction (Table 3).

The equation that software suggested to this experiment is as the follows:

Y = 5.46 + 0.73 A - 0.65 B + 0.60C - 0.34 AB + 0.33 AC- 1.38 A2 - 1.59 B2 - 0.58 C2

where Y is the dl-α-tocopherol amount in dry weight of yeast cells. A, B, and C are initial concentration of dl-α-tocopherol in culture medium, shaking speed, and concentration of sucrose added to SDB medium, respectively. The coefficients of A and C are positive that means positive effects of dl-α-tocopherol and sucrose concentrations on dl-α-tocopherol entrapment in S. cerevisiae cells. Negative coefficient of B means that the dl-α-tocopherol amount in dry cell weight of S. cerevisiae increases when shaking speed decreases.

The linear and quadratic coefficients of AB, A2, B2 and C2 are negative with p < 0.05, showing the significant influence of these factors. Interaction between dl-α-tocopherol concentration and shaking speed (AB) is more significant than interaction between dl-α-tocopherol concentration and sucrose (AC). The interaction between shaking speed and sucrose is insignificant (p > 0.05). Therefore, the modified quadratic polynomial model was used for optimization.

Three-dimentional rosponce plots

One of the most efficient tools in analysis of interactions is using graphical plots, especially three-dimensional response surface plot. In Figure 3, the relationship between shaking speed and initial dl-α-tocopherol concentrations investigated while the sucrose concentration was constant at the middle value. This plot indicates the quadratic coefficients of shaking speed and initial dl-α-tocopherol concentration is significant. The convex shape of three-dimentional plot reveals that we could calculate the optimal point of the response in the range of parameters. Figure 4 shows that when the sucrose and initial dl-α-tocopherol concentrations raised, the entrapment of dl-α-tocopherol in yeast cells increased and reached to an optimum level. The shape of plot reveals that the interaction between dl-α-tocopherol concentration and shaking speed is also significant. By raising the shaking speed and concentrations of sucrose and dl-α-tocopherol, the enriched cells of S. cerevisiae increased to an optimum point and then decreased. In Figure 5, the relationship between shaking speed and sucrose added to medium was investigated while the initial dl-α-tocopherol concentration was constant at the middle value. This plot shows that interaction between shaking speed and sucrose added to media is significant too. When the sucrose and shaking speed raised, the entrapment of dl-α-tocopherol in yeast cells increased and reached to an optimum level.

Optimization value

The optimum points of independent parameters to maximize the response were dl-α-tocopherol initial concentration of 7625.82 µg/mL, sucrose concentration of 6.86% w/v and shaking speed of 137.70 rpm. When the unstudied parameters were considered constant at the middle point, the maximum level of dl-α-tocopherol in biomass dry cell weight was predicted to be 5.866 µg/g.

Confirmation of the model

For validation of the model, an experiment was performed with the variable values suggested by the software. By applying the optimum conditions, 5.74 µg of dl-α-tocopherol was extracted from each gram of dried yeast cells; approximately this was 97.8% of the value predicted by the modified quadratic model (Table 4). The observed values were in good agreement with the predicted value. These data confirm adequacy and significance of the model.

Conclusions

This study is the first report about using S. cerevisiae as an accumulator and transporter of a lipid-soluble vitamin, dl-α-tocopherol. Our results showed that dl-α-tocopherol absorption in S. cerevisiae cells was affected significantly by shaking speed, sucrose, and dl-α-tocopherol concentrations. The modified quadratic model obtained from Design-Expert® software was adequate (p < 0.001). The resemblance between the R-squared and adjusted R-squared and acceptable value of C.V% revealed acceptability and accuracy of model. Mechanisms that govern movement of this lipid-soluble substance along the plasma membrane of eukaryotic cells need further studies.

References

  • 1.

    Stehlik-Tomas V, Zetić VG, Stanzer D, Grba S, Vahčić N. Zinc, copper and manganese enrichment in yeast Saccharomyces cerevisiae. Food Technol. Biotechnol. 2004;42:115-20.

  • 2.

    Volkov V. Quantitative description of ion transport via plasma membrane of yeast and small cells. Front. Plant Sci. 2015;6:425. [PubMed ID: 26113853].

  • 3.

    Lodolo EJ, Kock JLF, Axcell BC, Brooks M. The yeast Saccharomyces cerevisiae: The main character in beer brewing. FEMS Yeast Res. 2008;8:1018-36. [PubMed ID: 18795959].

  • 4.

    Fazeli MR, Bahmani S, Jamalifar H, Samadi N. Effect of probiotication on antioxidant and antibacterial activities of pomegranate juices from sour and sweet cultivars. Nat. Prod. Res. . 2011;25:288-97. [PubMed ID: 21294041].

  • 5.

    Moslehi-Jenabian S, Pedersen LL, Jespersen L. Beneficial effects of probiotic and food borne yeasts on human health. Nutrients. 2010;2:449-73. [PubMed ID: 22254033].

  • 6.

    European Food Safety Authority (EFSA). Scientific opinion on the safety and efficacy of MycoCell (Saccharomyces cerevisiae) as a feed additive for dairy cows. EFSA J. 2014;12:3830.

  • 7.

    European Food Safety Authority (EFSA). Scientific opinion on the maintenance of the list of QPS biological agents intentionally added to food and feed (2013 update). EFSA J. 2013;11:3449.

  • 8.

    Ernst IMA, Pallauf K, Bendall JK, Paulsen L, Nikolai S, Huebbe P, Roeder T, Rimbach G. Vitamin E supplementation and lifespan in model organisms. Ageing Res. Rev. 2013;12:365-75. [PubMed ID: 23099151].

  • 9.

    Traber MG, Atkinson J. Vitamin E, antioxidant and nothing more. Free Radic. Biol. Med. 2007;43:4-15. [PubMed ID: 17561088].

  • 10.

    Bashir S, Hasan SS, Ahmed SN, Basar T. Role of alpha tocopherol as an adjuvant therapy in pregnancy induced hypertension. Pak. J. Pharmacol. 2011;28:59-64.

  • 11.

    Ueda K, Zhao J, Kim HJ, Sparrow JR. Photodegradation of retinal bisretinoids in mouse models and implications for macular degeneration. Proc. Natl. Acad. Sci. U. S. A. 2016;113:6904-6909. [PubMed ID: 27274068].

  • 12.

    Das Gupta S, Suh N. Tocopherols in cancer: An update. Mol. Nutr. Food Res. . 2016;60:1354-63. [PubMed ID: 26751721].

  • 13.

    Engin KN. Alpha-tocopherol: looking beyond an antioxidant. Mol. Vis. 2009;15:855-60. [PubMed ID: 19390643].

  • 14.

    Engin KN, Engin G, Kucuksahin H, Oncu M, Engin G, Guvener B. Clinical evaluation of the neuroprotective effect of alpha tocopherol against glaucomatous damage. Eur. J. Ophthalmol. 2007;17:528-33. [PubMed ID: 17671926].

  • 15.

    Ricciarelli R, Argellati F, Pronzato MA, Domenicotti C. Vitamin E and neurodegenerative diseases. Mol. Aspects Med. 2007;28:591-606. [PubMed ID: 17306357].

  • 16.

    Farrokhifar SH, Ali Jafari R, Erfani Majd N, Fatemi Tabatabaee SR, Mayahi M. Effects of dietary vitamin E on mucosal maltase and alkaline phosphatase enzyme activities and on the amount of mucosal malonyldialdehyde in broiler chickens. Vet. Res. Forum . 2013;4:221-5. [PubMed ID: 25568675].

  • 17.

    Castellani ML, Shaik-Dasthagirisaheb YB, Tripodi D, Anogeianaki A, Felaco P, Toniato E, De Lutiis MA, Fulcheri M, Tetè S, Galzio R, Salini V, Caraffa A, Antinolfi P, Frydas I, Sabatino G, Kempuraj D. Interrelationship between vitamins and cytokines in immunity. J. Biol. Regul. Homeost. Agents . 2010;24:385-90. [PubMed ID: 21122276].

  • 18.

    Ajandouz el H, Castan S, Jakob S, Puigserver A. A fast, sensitive HPLC method for the determination of esterase activity on α-tocopheryl acetate. J. Chromatogr. Sci. . 2006;44:631-3. [PubMed ID: 17254374].

  • 19.

    Kaur T, Bansal MP. Selenium enrichment and anti-oxidant status in baker’s yeast, Saccharomyces cerevisiaeat different sodium selenite concentrations. Nutr. Hosp. . 2006;21:704-8. [PubMed ID: 17147070].

  • 20.

    Gaensly F, Wille GMFC, Brand D, Bonfim TMB. Iron enriched Saccharomyces cerevisiae maintains its fermenting power and bakery properties. Ciênc. Tecnol. Aliment. 2011;31:980-983.

  • 21.

    Esmaeili S, Khosravani-Darani K, Pourahmad R, Komeili R. An experimental design for production of selenium-enriched yeast. World Appl. Sci. J. 2012;19:31-7.

  • 22.

    Arshad MS, Anjum FM, Asghar A, Khan MI, Yasin M, Shahid M, El-Ghorab AH. Lipid stability and antioxidant profile of microsomal fraction of broiler meat enriched with α-lipoic acid and α-tocopherol acetate. J. Agric. Food Chem. . 2011;59:7346-52. [PubMed ID: 21557618].

  • 23.

    Yin H, Fan G, Gu Z. Optimization of culture parameters of selenium-enriched yeast (Saccharomyces cerevisiae) by response surface methodology (RSM). Food Sci. Tech. 2010;43:666-9.

  • 24.

    Khoshayand F, Goodarzi S, Shahverdi AR, Khoshayand MR. Optimization of culture conditions for fermentation of soymilk using Lactobacillus casei by response surface methodology. Probiotics Antimicrob. Proteins . 2011;3:159-67. [PubMed ID: 26781676].

  • 25.

    Popović S, Kostadinović LM, Brkljača JS, Krulj JA, Manojlović MS, Solarov MIB. The development and validation of HPLC method for quantification of DL-α-tocopherol in quinoa seeds (Chenopodium quinoa Willd). Food Feed Res. 2014;41:147-52.

  • 26.

    Box GEP, Behnken DW. Some new three level designs for the study of quantitative variables. Technometrics. 1960;2:455-75.