NC State
BioResources
Wu, J., Wang, Y., Zhang, Z., Lu, X., Wang, C., Zhong, Z., Qian, D., and Zhang, T. (2020). "Optimization of pretreating poplar wood shavings with cellulase to produce binderless fiberboard using response surface methodology," BioRes. 15(3), 6431-6441.

Abstract

Wood processing waste, poplar wood shavings, were used for fiberboard production, and the pretreatment conditions using cellulase were studied using response surface methodology (RSM). After single factors optimization, central level of temperature, dosage, and liquid/solid ratio (LSR) of cellulase pretreatment conditions were obtained. Further optimization to study the influence of the factors was carried out using Box-Behnken design of experiments. A second-order polynomial equation was obtained, and the low p-value (<0.007) implied that the model was highly significant by analysis of variance (ANOVA). The optimized cellulase pretreatment conditions for maximum bending strength (BS) of the fiberboard were determined by ridge analysis as 44.4 °C of temperature, 1.23 U/g of dosage, 4.2 of LSR, and 5 h of pretreatment time. Under the optimized conditions, the BS of the fiberboard reached 25.12 ± 0.35 MPa by validation experiment, which was twice that of the fiberboard without pretreatment. Thus, the cellulase pretreatment should be a good choice to produce high-strength binderless fiberboard.


Download PDF

Full Article

Optimization of Pretreating Poplar Wood Shavings with Cellulase to Produce Binderless Fiberboard Using Response Surface Methodology

Jianguo Wu,a,b,* Yuanyuan Wang,a Zhao Zhang,a Xuanxian Lu,a Chi Wang,a Zhiqing Zhong,a Dong Qian,a and Tong Zhang a,b

Wood processing waste, poplar wood shavings, were used for fiberboard production, and the pretreatment conditions using cellulase were studied using response surface methodology (RSM). After single factors optimization, central level of temperature, dosage, and liquid/solid ratio (LSR) of cellulase pretreatment conditions were obtained. Further optimization to study the influence of the factors was carried out using Box-Behnken design of experiments. A second-order polynomial equation was obtained, and the low p-value (<0.007) implied that the model was highly significant by analysis of variance (ANOVA). The optimized cellulase pretreatment conditions for maximum bending strength (BS) of the fiberboard were determined by ridge analysis as 44.4 °C of temperature, 1.23 U/g of dosage, 4.2 of LSR, and 5 h of pretreatment time. Under the optimized conditions, the BS of the fiberboard reached 25.12 ± 0.35 MPa by validation experiment, which was twice that of the fiberboard without pretreatment. Thus, the cellulase pretreatment should be a good choice to produce high-strength binderless fiberboard.

Keywords: Fiberboard; Binderless; Cellulase; Pretreatment; Poplar wood shavings; Response surface methodology; Optimization

Contact information: a: Jiangsu Key Laboratory for Eco-Agricultural Biotechnology around Hongze Lake, Jiangsu Key Construction Laboratory for Food Safe and Nutritional Function, School of Life Science, Huaiyin Normal University, Changjiang West Road 111, Huai’an 223300, China; b: Jiangsu Collaborative Innovation Center of Regional Modern Agriculture and Environmental Protection, Huaiyin Normal University, Changjiang West Road 111, Huai’an 223300, China;

* Corresponding author: jgwu_21@sina.cn

INTRODUCTION

As traditional fiberboard production widely employs formaldehyde- and other aldehyde-based adhesives, toxic and harmful gases that are released during the production seriously pollute the environment and endanger human health (Bertrand et al. 2007; Lithner et al. 2011; Piekarski et al. 2017; Pang et al. 2018). Consequently, eco-friendly binderless fiberboard production technology has been growing in recent years. Geng et al. (2006) found that the fiberboards made from alkaline-treated bark showed lighter color, higher internal bonding strength, higher modulus of rupture, and higher modulus of elasticity in comparison with the control panels. Nikhom and Okuma (2000) used oil palm frond as raw material to produce viscous substances by high temperature steam explosion, and they further studied the binderless fiberboard production technology. Binderless fiberboards using Triarrhena sacchariflora residue, corn stalk, and poplar wood shavings by bio-pretreatment with white-rot fungus Trametes hirsute or Trametes versicolor have been studied (Wu et al. 2011, 2016, 2020). However, physical and chemical methods had the disadvantages of high pollution and energy consumption, while bio-pretreatment methods had the disadvantages of long pretreatment time. Hence, more attention has been paid to enzyme pretreatment (Widsten and Kandelbauer 2008).

Kharazipour and Mai (1998) first pretreated wood fiber with laccase without any traditional adhesive, and then they obtained fiberboard with certain mechanical strength by dry pressing and wet pressing, respectively. Felby et al. (2002) and Euring et al. (2011) pretreated raw materials by laccase, and they carried out small-scale factory production experiments of medium density fiberboard. Ibrahim et al. (2013) used laccase to modify kraft lignin to formulate adhesives combined with polyethyleneimine, chitosan, and soy protein, respectively. Wu et al. (2016) found that there was a positive correlation between the content of hydroxyl group and the mechanical strength of fiberboard. Cellulase could break the β-1,4-glycoside bond of cellulose into amorphous cellulose, cello-oligosaccharide, and small molecule sugars, which are rich in the hydroxyl group. Therefore, the cellulase pretreatment should be also a potential binderless fiberboard production method (Yu et al. 2003). However, still there are no relevant research reports available.

In China, the poplar wood shavings (PWS), the waste from poplar wood processing, are very easily available. PWS could be used for production of fiberboard, lipids, xylitol and bioenergy (Chen et al. 2016; Dalli et al. 2017; Li et al. 2019; Wu et al. 2020), but there have been no reports about binderless fiberboard using PWS by cellulase pretreatment. In this study, cellulase was firstly used to pretreat PWS for binderless fiberboard production. The cellulase pretreatment conditions were optimized using response surface methodology (RSM) (Gunny et al. 2013; Wang et al. 2018; Mohammed et al. 2020). The result would provide technical support for the production of binderless fiberboard by cellulase pretreatment.

EXPERIMENTAL

Materials and Agents

Well-distributed shavings of PWS collected from a wood processing factory in Huai’an (China), were dried at 60 °C in an oven to a constant weight. Cellulase was purchased from Huai’an (Bio-mass Green Bioenergy Co. Ltd, China). The filter paper activity (FPA), carboxymethylcellulose activity (CMC), β-glucosidase activity, and xylanase activity of the cellulase were about 350 U/g, 4000 U/g, 60 U/g, and 20000 U/g, respectively. Acetate buffer (pH 4.8) was prepared by mixing 12 mL of sodium acetate solution (1.0 M) and 8 mL of acetic acid solution (1.0 M) in 1 L purified water. Acetic acid and sodium acetate used were all of analytical grade and purchased from the Sinopharm Chemical Reagent Co. Ltd., China.

Cellulase Pretreatment

Approximately 60 g of PWS was immersed in 300 mL of acetate buffer with 60 U of cellulase (FPA). The pretreatment temperature was controlled at 50 °C by a thermostat (uncertainty of ± 0.1 °C) in the shaker. After 2 h of pretreatment, the PWS were loaded into six layers of gauze and squeezed by hands to about 60% moisture content to prepare fiberboard by hot pressing.

Single Factor Optimization

The initial single factor conditions used was 50 °C of temperature (T), 1.0 FPA U/g of PWS of enzyme dosage, liquid/solid ratio (LSR) of 5, and 2 h of pretreatment time (t). The single factor optimization method is shown in Table 1.

Table 1. Single Factor Optimization Method

Response Surface Methodology Optimization

Factor T, dosage, and LSR were optimized by response surface methodology (RSM) and the Box-Behnken design of experiments as shown in Table 2.

Table 2. Level of Variables for the Box-Behnken Design of Experiments

Preparation of Fiberboard

The well-distributed PWS were manually loaded into a stainless mold (100 mm × 100 mm × 3 mm), then pressed flat into a fiberboard by a R-3202 Hot-Press Model (Wuhan Qien Science and Technology Development Co. Ltd., Hubei, China). Hot pressing was carried out at 170 °C under 17 MPa of pressure for 10 min.

Property of Fiberboard

The length, width, and thickness of fiberboard were measured by a micrometer, and the mass of fiberboard was measured by electronic balance (Model BSA124S, Sartorius, Gottingen, Germany), so the density of fiberboard was calculated by dividing the fiberboard mass (kg) by its volume (m3). The fiberboard was cut into specimens of 100 mm in length and 20 mm in width. The bending strength (BS) of fiberboard specimens was determined by performing the three point flex test with 60 mm of span length at 3 mm/min crosshead speed using a micro-computer controlled electron universal testing machine (Model ETM104B, Wance Group, Shenzhen, China). The BS could be calculated as Eq. (1). The water swelling ratio was measured from the different thickness (h) before and after immersing in water for 24 h, and the water swelling ratio could be calculated as Eq. (2). Detailed procedures were carried out as shown in National Standard (GB/T 11718-2009) (Standardization Administration of China 2009).

 (1)

 (2)

where F is the maximum force to break the fiberboard specimens (N); L is the span length, = 60 mm; b is the width of fiberboard specimens (mm); h is the thickness of fiberboard specimens (mm); h2 is the final thickness of fiberboard after immersing in water (mm); and h1 is the initial thickness of fiberboard before immersing in water (mm).

Statistical Analysis

Test Pilot software (Wance Group, Shenzhen, China) was used to calculate BS, ORIGIN PRO 8.0 (OriginLab Corp., Wellesley, MA, USA) was used to analyze the data and draw bar graphs, and Statistical Analysis System 9.4 (SAS 9.4) (SAS Institute Inc., North Carolina, USA) was used to design experiment and analyze the data of RSM. The error bars in bar graphs correspond to standard errors of three tests.

RESULTS AND DISCUSSION

Single Factor Optimization of Cellulase Pretreatment Conditions

The results of enzyme dosage, LSR, temperature, and pretreatment time for BS of fiberboard are shown in Fig. 1. After PWS was pretreated by cellulase, it was made into fiberboard without adhesive. With the increase of dosage, LSR, and temperature, the BS of fiberboard increased first and then began to decrease. At a dosage of 1.25 U/g, LSR of 4, and at 45 °C, the BS of fiberboard indicated a peak value. The pretreatment time of cellulase had a positive correlation with the BS of fiberboard within 5 hours of pretreatment. The BS of fiberboard reached the maximum at 5 h of pretreatment time, and it decreased slightly after prolonging the pretreatment time.

Fig. 1. Sing factor optimization of cellulase pretreatment conditions

The interaction effect of dosage, LSR, and temperature might appear in the process of cellulase pretreatment. It was necessary to optimize these three factors by RSM. To ensure the pretreatment effect, the pretreatment time was kept for 5 h.

Optimization by Response Surface Methodology

Response surface methodology using Box-Behnken design (BBD) (Li et al. 2018) was made to determine the optimal levels of the three selected factors (T for x1, dosage for x2, LSR for x3). The experimental design and results are shown in Table 3.

Table 3. Box-Behnken Design with BS as Response

The relationships between BS and the tested factors were obtained by application of RSM. By employing multiple regression analysis on the experimental data, BS and the tested factors were related by the following second-order polynomial equation as Eq. (3). The 3D response surfaces and 2D contour plots could be generated by SAS 9.4 based on the Eq. (3).

 (3)

where Y was BS of the fiberboard, x1 was temperature (T), x2 was dosage, and x3 was LSR.

Table 4. Analysis of Variance (ANOVA) for the Second-Order Polynomial Model

The analysis of variance (ANOVA) data for the selected quadratic polynomial model is listed in Table 4. The high model F-value (11.406) and low p-value (<0.007) implied that the model was highly significant. The coefficient of determination R2 was calculated as 0.954, which implied that the model was reliable for the BS of fiberboard.

Table 5. Regression Results of the Box-Behnken Design

Fig. 2. 3D Response surfaces showing the effects of temperature (T), Dosage, and LSR on the BS of fiberboard

The coefficient estimates of model equation, along with the corresponding p-values, are presented in Table 5. The p-value was employed to check the significance of each coefficient, which also indicated the interactions between the variables (Liu et al. 2003; Sakr et al. 2019). The corresponding coefficient was more significant as the p-value was smaller. Temperature (x1) was highly significant with p-values less than 0.05, indicating that temperature was the most effective factor (Souri et al. 2018).

The 3D response surfaces and 2D contour plots (Figs. 2 and 3) generated by SAS 9.4 were the graphical representations of the regression equation (Eq. 3). These plots visualized the relationship between the response and each variable, and the interactions between two tested variables. The 3D response surfaces and 2D contour plots could also locate the optimum values of the factors for the maximum of the response. The maximum predicted response was indicated by the surface confined in the smallest ellipse in the contour diagram (Muralidhar et al. 2001; He et al. 2013; Zuorro et al. 2019). Peaks were found in 3D response surfaces (Fig. 2), and the smallest ellipse was also found in 2D contour plots (Fig. 3), which indicated that there was indeed an optimal condition to get the maximum BS of the fiberboard.

Fig. 3. 2D Contour plots showing the effects of temperature (T), Dosage, and LSR on the BS of fiberboard

Optimization Results and Verification

The optimum conditions for cellulase pretreatment could be determined by ridge analysis (Kim and Rhee 2015), and the results are shown in Table 6. Based on ridge analysis, the maximum BS of fiberboard could reach to 24.586 ± 0.275 MPa, when temperature was 44.4 °C, dosage was 1.23 U/g, and LSR was 4.20. To verify the predicted results, the validation experiments were performed in triplicate. Under the optimized condition, the BS of fiberboard reached 25.12 ± 0.35 MPa with implying that experimental and predicted values (24.586 MPa) were in good agreement. The BS of the fiberboard using PWS pretreated under the optimized cellulase pretreatment conditions was twice that of the fiberboard without pretreatment (12.35 MPa).

The density of fiberboard was 952 kg/m3, which indicated that the fiberboard belonged to high density fiberboard based on the National Standard of the People’s Republic of China (GB/T 31765-2015) (Standardization Administration of China 2015). However, BS of the fiberboard was lower than that of national standard (37 MPa). Water swelling ratio of the fiberboard was only 11.8%, which reached the national standard (<12%). Comparing the national standard of medium density fiberboard (GB/T 11718-2009) (Standardization Administration of China 2009), water swelling ratio of the fiberboard was better than that of the national standard of medium density fiberboard (35%), and BS of the fiberboard reached the national standard (25 MPa). In addition, the BS of the transversely cutting specimen, longitudinally cutting specimen, and 45° obliquely cutting specimen were 25.3 MPa, 25.0 MPa, and 24.8 MPa, respectively. The fiberboard could be approximated as an in-plane isotropic material. Consequently, the cellulase pretreatment should be a good choice to produce high-strength binderless fiberboard, and the fiberboard could be used as high quality medium density fiberboard. Hydrolysis of cellulose by cellulase could increase the content of hydroxyl group in material, which was beneficial to the strength of fiberboard. Steam explosion or acid hydrolysis could also increase the hydroxyl group availability in materials. Whether it can also be used to improve the strength of fiberboard? That will be an interesting research topic.

Table 6. Ridge Analysis for BS as Response

CONCLUSIONS

The optimization condition of cellulase pretreatment was obtained as 44.4 C of temperature, 1.23 U/g of dosage, 4.2 of LSR, and 5 h of pretreatment time. Under the optimized condition, the BS of fiberboard was reached 25.12±0.35 MPa, which was twice that of the fiberboard without pretreatment. Therefore, cellulase pretreatment should be a good choice to produce high strength binderless fiberboard.

ACKNOWLEDGMENTS

The authors are grateful for the support of the National Natural Science Foundation of China, Grant No. 31870543; and Qing Lan Project of Jiangsu Province of China; and Doctor Jianlong He for the suggestions on the modification of this paper.

REFERENCES CITED

Bertrand, J . P., Simon, V., and Chau, N. (2007). “Associations of symptoms related to isocyanate, ureaformol, and formophenolic exposures with respiratory symptoms and lung function in coal miners,” Int. J. Occup. Environ. Health 13(2), 181-187. DOI: 10.1179/oeh.2007.13.2.181

Chen, D., Li, Y., Cen, K., Luo, M., Li, H., and Lu, B. (2016). “Pyrolysis polygeneration of poplar wood: Effect of heating rate and pyrolysis temperature,” Bioresour. Technol. 218, 780-788. DOI: 10.1016/j.biortech.2016.07.049

Dalli, S. S., da Silva, S. S., Uprety, B. K., and Rakshit, S. K. (2017). “Enhanced production of xylitol from poplar wood hydrolysates through a sustainable process using immobilized new strain Candida tropicalis UFMG BX 12-a,” Appl. Biochem. Biotechnol. 182(3), 1053-1064. DOI: 10.1007/s12010-016-2381-4

Euring, M., Ruhl, M., Ritter, N., Kues, U., and Kharazipour, A. (2011). “Laccase mediator systems for eco-friendly production of medium-density fiberboard (MDF) on a pilot scale: Physicochemical analysis of the reaction mechanism,” Biotechnol. J. 6(10), 1253-1261. DOI: 10.1002/biot.201100119

Felby, C., Hassingboe, J., and Lund, M. (2002). “Pilot-scale production of fiberboards made by laccase oxidized wood fibers: board properties and evidence for cross-linking of lignin,” Enzyme Microb. Tech. 31(6), 736-741. DOI: 10.1016/S0141-0229(02)00111-4

Geng, X., Zhang, S. Y., and Deng, J. (2006). “Alkaline treatment of Black Spruce Bark for the manufacture of binderless fiberboard,” J. Wood Chem. Technol. 26(4), 313-324. DOI: 10.1080/02773810601076857

Gunny, A. A., Arbain, D., and Sithamparam, L. (2013). “Optimization of complex fermentation media for glucose oxidase production using statistical approach,” Pak. J. Biol. Sci. 16(18), 960-964. DOI: 10.3923/pjbs.2013.960.964

He, S., Wang, H., Wu, B., Zhou, H., Zhu, P., Yang, R., and Yan, X. (2013). “Response surface methodology optimization of fermentation conditions for rapid and efficient accumulation of macrolactin A by marine Bacillus amyloliquefaciens ESB-2,” Molecules 18, 408-417. DOI: 10.3390/molecules18010408

Ibrahim, V., Mamo, G., Gustafsson, P. J., and Hatti-Kaul, R. (2013). “Production and properties of adhesives formulated from laccase modified Kraft lignin,” Ind. Crop. Prod. 45, 343-348. DOI: 10.1016/j.indcrop.2012.12.051

Kharazipour, A., and Mai, C. A. H. (1998). “Polyphenoles for compounded materials,” Polym. Degrad. Stabil. 59, 237-243. DOI: 10.1016/S0141-3910(97)00157-2

Kim, S. A. and Rhee, M. S. (2015). “Predictive model and optimization of a combined treatment of caprylic acid and citric acid for the reduction of Escherichia coli O157:H7 using the response surface methodology,” Int. J. Food Microbiol. 197, 9-14. DOI: 10.1016/j.ijfoodmicro.2014.12.006

Li, H., van den Driesche, S., Bunge, F., Yang, B., and Vellekoop, M. J. (2018). “Optimization of on-chip bacterial culture conditions using the Box-Behnken design response surface methodology for faster drug susceptibility screening,” Talanta 194, 627-633. DOI: 10.1016/j.talanta.2018.10.048

Li, X., He, Y., Zhang, L., Xu, Z., Ben, H., Gaffrey, M.J., Yang, Y., Yang, S., Yuan, J.S., Qian, W.J., and Yang, B. (2019). “Discovery of potential pathways for biological conversion of poplar wood into lipids by co-fermentation of Rhodococci strains,” Biotechnol. Biofuels 12, 60. DOI: 10.1186/s13068-019-1395-x

Lithner, D., Larsson, A., and Dave, G. (2011). “Environmental and health hazard ranking and assessment of plastic polymers based on chemical composition,” Sci. Total Environ. 409(18), 3309-3324. DOI: 10.1016/j.scitotenv.2011.04.038

Liu, J. Z., Weng, L. P., Zhang, Q. L., Xu, H., and Ji, L. N. (2003). “Optimization of glucose oxidase production by Aspergillus niger in a benchtop bioreactor using response surface methodology,” World J. Microb. Biot.19, 317-323. DOI: 10.1023/A:1023622925933

Mohammed, A. H., Mohammad, O. H., and Mehmet, F. O. (2020). “Optimization of a confined jet geometry to improve film cooling performance using response surface methodology (RSM),” Processes 8(2), 232. DOI: 10.3390/pr8020232

Muralidhar, R. V., Chirumamila, R. R., Marchant, R., and Nigam, P. (2001). “A response surface approach for the comparison of lipase production by Candida cylindracea using two different carbon sources,” Biochem. Eng. J. 9, 17-23. DOI: 10.1016/S1369-703X(01)00117-6

Nikhom, L., and Okuma, M. (2000). “Development of boards made from oil palm frond II: Properties of binderless boards from steam-exploded fibers of oil palm frond,” J. Wood Sci. 46, 322-326. DOI: 10.1007/BF00766224

Pang, B., Li, M. K., Yang, S., Yuan, T. Q., Du, G. B., and Sun, R. C. (2018). “Eco-friendly phenol-urea-formaldehyde co-condensed resin adhesives accelerated by resorcinol for plywood manufacturing,” ACS Omega 3(8), 8521-8528. DOI: 10.1021/acsomega.8b01286

Piekarski, C. M., de Francisco, A. C., da Luz, L. M., Kovaleski, J. L., and Silva, D. A. L. (2017). “Life cycle assessment of medium-density fiberboard (MDF) manufacturing process in Brazil,” Sci. Total Environ. 575, 103-111. DOI: 10.1016/j.scitotenv.2016.10.007

Sakr, M., Hanafi, R., Fouad, M., Al-Easa, H., and El-Moghazy, S. (2019). “Design and optimization of a luminescent Samarium complex of isoprenaline: A chemometric approach based on Factorial design and Box-Behnken response surface methodology,” Spectrochim. Acta A 208, 114-123. DOI: 10.1016/j.saa.2018.09.061

Souri, M., Hoseinpour, V., Shakeri, A., and Ghaemi, N. (2018). “Optimisation of green synthesis of MnO nanoparticles via utilising response surface methodology,” IET Nanobiotechnol. 12(6), 822-827. DOI: 10.1049/iet-nbt.2017.0145

Standardization Administration of China. (2015) “National standard of the people’s Republic of China: High density fiberboard (GB/T 31765-2015),” Beijing: Standards Press of China.

Standardization Administration of China. (2009) “National standard of the people’s Republic of China: Medium density fiberboard (GB/T 11718-2009),” Beijing: Standards Press of China.

Wang, W., Cheng, Y., and Tan, G. (2018). “Design optimization of SBS-modified asphalt mixture reinforced with eco-friendly basalt fiber based on response surface methodology,” Materials (Basel) 11(8), 1311. DOI: 10.3390/ma11081311

Widsten, P., and Kandelbauer, A. (2008). “Adhesion improvement of lignocellulosic products by enzymatic pre-treatment,” Biotechnol. Adv. 26(4), 379-386. DOI: 10.1016/j.biotechadv.2008.04.003

Wu, J., Chen, C., Zhang, H., Xia, L., Huang, Y., Huang, H., Wang, Y., Qian, D., Wang, J., Wang, X., and Zhang, T. (2020). “Eco-friendly fiberboard production without binder using poplar wood shavings bio-pretreated by white rot fungi Coriolus versicolor,” Constr. Build. Mater. 236, 117620. DOI: 10.1016/j.conbuildmat.2019.117620

Wu, J., Zhang, X., Liu, J., Xiong, M., Lu, X., Fan, H., Wang, X., and Zhang, X. (2016). “Medium density fibreboard production by hot pressing without adhesive using Triarrhena sacchariflora residue bio-pretreated by white-rot fungus Coriolus versicolor,” J. Appl. Microbiol. 121(2), 415-421. DOI: 10.1111/jam.13148

Wu, J., Zhang, X., Wan, J., Ma, F., Tang, Y., and Zhang, X. (2011). “Production of fiberboard using corn stalk pretreated with white-rot fungus Trametes hirsute by hot pressing without adhesive,” Bioresource Technol. 102(24), 11258-11261. DOI: 10.1016/j.biortech.2011.09.097

Yu, P., McKinnon, J. J., Maenz, D. D., Olkowski, A. A., Racz, V. J., and Christensen, D. A. (2003). “Enzymic release of reducing sugars from oat hulls by cellulase, as influenced by Aspergillus ferulic acid esterase and Trichoderma xylanase,” J. Agric. Food Chem. 51(1), 218-223. DOI: 10.1021/jf020476x

Zuorro, A., Levecchia, R., González-Delgado, Á. D., García-Martinez, J. B., and L’Abbate, P. (2019). “Optimization of enzyme-assisted extraction of flavonoids from corn husks,” Processes 7(11), 804. DOI: 10.3390/pr7110804

Article submitted: March 14, 2020; Peer review completed: May 24, 2020; Revised version received and accepted: May 30, 2020; Published: July 6, 2020.

DOI: 10.15376/biores.15.3.6431-6441