Next Article in Journal
Economic Growth, Income Inequality and Food Safety Risk
Previous Article in Journal
Research Advances of Modification and Nutraceutical Properties of Polysaccharide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of the Heat-Drying Conditions of Drone Pupae by Response Surface Methodology (RSM)

1
Research Center for Environmentally Friendly and Quality Livestock Production Technology, Hankyong National University, Anseong-si 17579, Gyeonggi-do, Republic of Korea
2
School of Animal Life Convergence Science, Hankyong National University, Anseong-si 17579, Gyeonggi-do, Republic of Korea
*
Author to whom correspondence should be addressed.
Foods 2023, 12(16), 3062; https://doi.org/10.3390/foods12163062
Submission received: 19 June 2023 / Revised: 6 August 2023 / Accepted: 14 August 2023 / Published: 15 August 2023
(This article belongs to the Topic Sustainable Food Processing)

Abstract

:
Recent research has been conducted on various types of pre-processing methods for insects, including freeze-drying, microwave drying, hot air heat drying, and non-heat drying. This study aimed to identify the factors that have the greatest impact on heat drying conditions and establish the optimal heat drying conditions for drone pupae (Apis melifera L.) using response surface methodology (RSM) to minimize quality changes. Drone pupae were treated under various conditions, including blanching time (53–187 s) (X1), drying temperatures (41.6–58.4 °C) (X2), and drying time (266–434 min) (X3). The effect of these treatments on response variables, including the color parameter (WI, YI, BI, △E, and BD), AV, and TB of the dried drone pupae, was evaluated using a central composite design. The whole design consisted of 20 experimental points carried out in random order, which included eight factorial points, six center points, and six axial points. The optimal drying conditions for drone pupae were determined to be a blanching time of 58 s, a drying temperature of 56.7 °C, and a drying time of 298 min. The response variables were most affected by drying temperature and drying time and to a lesser extent by blanching time. The processed drone pupae using the optimized drying conditions resulted in the color parameters (WI, BI, YI, ΔE, and BD) being found to be 66.67, 21.33, 26.27, 31.27 and 0.13, respectively. And TB (log CFU/g) and AV (mg/g) values were found to be 3.12 and 4.33, respectively. The estimated and actual values for dried drone pupae showed no significant difference (p < 0.05). Comparing the physicochemical and microbiological properties of freeze-dried and optimal heat-dried drone pupae, the L and b value as well as PV were significantly lower in the heat-dried samples, while no significant difference was observed in the a value and AV (p < 0.05). Our study suggests that the model we developed can be applied to the large-scale production of drying conditions for use in the pharmaceutical and food industries.

1. Introduction

With the increased population and demand for higher protein sources, edible insects are gaining recognition as an alternative protein source [1]. The United Nations Food and Agriculture Organization (FAO) has proposed insect consumption as a viable substitute for meat [2]. It is projected that the edible insect market in the United States and Europe will reach approximately $3.3 billion by 2025, and in Asia, it is expected to grow to approximately $270 million by 2024. The global insect industry is anticipated to experience significant growth, with the market size expanding to tens of billions of dollars, particularly in sectors that utilize insects for applications such as new materials and pharmaceutical development [3].
Insects have advantages as a protein substitute food, but the challenge of overcoming consumer perception due to disgust for insects remains. To easily apply insects to food, insects can be processed into food ingredients such as insect powder or protein extract [4], but pre-processing is required for this. Pre-processing methods have a significant impact on the final product’s lipid oxidation, coloration, and microbial load. Insects undergo various reactions that cause color changes, including enzymatic polyphenol oxidation, complex formation between iron and polyphenols, and the formation of Maillard reaction products [5,6,7]. Traditional methods of pre-processing insects include roasting and sun-drying [8], and further research has been conducted recently on various types of pre-processing methods such as freeze-drying, microwave drying, hot-air heat drying, and non-heat drying. Browning is one of the main issues in several edible insect species [9]. Several enzymes, including phenol oxidase, lactase, tyrosine hydroxylase, decarboxylase, and peroxidase, play a role in the browning process. In insects, most phenolic compounds are produced from L-tyrosine in the shikimic acid pathway and are often modified by enzymes. Phenolic compounds can also react with proteins and amino acids non-enzymatically, resulting in the formation of dark compounds. Blanching is a widely used food pre-processing method for preparing insects in the food industry, and it is a common treatment aimed at deactivating enzymes, which can lead to undesirable changes in the food product. During the blanching process, samples are exposed to high temperatures ranging from 70 to 150 °C for a short period of time using hot water or steam [10,11]. Blanching is a method of immersing the material in boiling water, and when it is applied to sacrificing insects, hot water quickly penetrates into the insects, inhibiting denaturation during sacrifice, and it can also achieve enzyme inactivation and some sterilization effects, making it a sacrificial method with many advantages.
Bees are used as food in many countries around the world, such as Thailand and India [8,12,13,14,15], and in the case of bee larvae, there is evidence that it was used as a medicinal material in traditional East Asian medicine, such as ancient traditional Chinese therapies and Korean traditional medicine. Recently, research has been conducted on the drone pupae of the Apis melifera L. as a food ingredient in South Korea [16,17,18,19,20,21,22]. Drone pupae contain large amounts of major nutrients, amino acids, minerals, etc., and they have reported antioxidant, blood sugar enhancing, and hair loss improvement effects. However, after harvesting, the rich nutrients of drone pupae can be affected by lipid oxidation and microbial contamination depending on the storage method. After drone pupae are harvested, they should be stored in a freezer or freeze-dried, but most bee farms in Korea are small-scale farms with difficulties in facility investment. Therefore, it is necessary to develop a low-cost pre-processing method to reduce the production cost of drone pupae and to create high-value-added materials.
Response surface methodology (RSM) is defined as a statistical analysis method for the surface of a reaction where a dependent variable is affected by multiple independent variables with complex interactions [23]. These independent variables are factors that affect the response, and the dependent variable is a response variable that represents a reaction affected by the factors. The statistical approach involves setting up an appropriate experimental design based on independent variables (X1, X2, X3…) to obtain measured values of dependent variables (Y1, Y2, Y3…) and analyzing these variables to achieve the purpose of the study. This experimental design analysis is widely used in various aspects such as cost reduction, process development, product development, method development, and quality control. To optimize it, a one-factor-at-a-time method was traditionally applied in which all variables except one were fixed and only one variable was changed at a time to observe its effect. However, in cases where various independent variables interact with each other, such as in food research, a method of changing only one variable to study the interaction requires a lot of time and expense. Therefore, RSM, which can observe interactions between factors and save time and expenses, is widely used.
The purpose of this study was to investigate the factors that have the greatest influence on heat drying and to establish the optimal heat drying conditions for drone pupae using response surface methodology (RSM) to minimize quality changes such as discoloration and acidity according to these factors. To achieve this, we conducted experiments with blanching time, drying temperature, and drying time as independent variables for heat drying conditions, with color parameters, TB, and AV set as dependent variables. In addition, microbial and acidity changes were analyzed in drone pupae dried using optimized heat drying conditions and in those that were freeze-dried.

2. Materials and Methods

2.1. Materials

Ethanol, phenolphthalein solution, ethyl ether, and chloroform were obtained from Daejung Chemical & Metals Co., Ltd. (Seoul, South Korea). Potassium hydroxide solution (in ethanol) and acetic acid were obtained from Duksan Science (Seoul, South Korea). Potassium iodide and sodium thiosulfate pentahydrate were acquired from Sigma-Aldrich (St. Louis, MO, USA).
The drone pupae of Apis mellifera L. used in the experiment were raised in bee hives specifically at a beekeeping farm in Jeonbuk, South Korea in May 2022. The bee hives were removed after 20–23 days and immediately transferred to a freezer for storage before use in the experiment. The drone pupae in the bee hives were manually crushed, and they were separated and used. The drone pupae were blanched for 53–187 s in boiling water with a water:pupae ratio of 100:1 (v:w). After blanching, drone pupae were cooled in cold (4 °C) water to interrupt heat processes. And then, they were dried a laboratory oven for 266–434 min at 41.6–58.4 °C and used for experimental sample. For freeze-drying, the drone pupae blanched for 120 s in boiling water was first frozen at −70 °C for 48 h. The frozen sample was then freeze-dried using a freeze dryer (FD8512P, Ilshin-bio Co., Ltd., South Korea) at a temperature of −30 °C and pressure of 120 Pa for 48 h.

2.2. Color Parameters

2.2.1. Color Indices

The color of each experimental sample was measured using a colorimeter (Color i7, X-rite, Grand Rapids, MI, USA) by taking five repeated measurements of the L* (lightness), a* (redness), and b* (yellowness) value. The standard white plate had an L value of 95.75, an a value of 0.15, and a b value of 2.82. The total color change (△E) was calculated using Equation (1). The L value difference (△L), a value difference (△a), and b value difference (△b) were analyzed using the initial value color difference. The whiteness (WI), yellowness (YI) and browning (BI) indices were calculated according to Tellez-Morales et al. [23] following Equations (2)–(4):
E = ( L * ) 2 + ( a * ) 2 + ( b * ) 2
WI = ( 100 L * 2 ) + a * 2 + b * 2  
YI = 142.86 × b * 2   L * 2  
BI = [ 100 × ( X 0.31 ) ] 0.17     where   X = ( a * + 1.75 × L * ) ( 5.645 × L *   + a *   3.012 × b * )  

2.2.2. Browning Degree (BD)

BD was measured according to the method by Chang et al. [24]. Experimental samples were extracted by adding 5 g of the sample to 50 mL of 50% ethanol and extracting for 24 h at room temperature. The extracted solution was filtered through filter paper (Whatman No. 2), and the absorbance was measured at 420 nm using a UV spectrophotometer (Optizen 2120UV, Mecasys Co., Ltd., Daejeon, Republic of Korea).

2.3. Acid Value (AV)

AV was measured according to the method by Kimetal [18]. A five-gram sample was mixed with a 100 mL mixture of ethyl ether and ethanol (2:1, v/v), shaken, and then 100 μL of a 1% phenolphthalein solution was added as an indicator. The amount of 0.1 N potassium hydroxide solution required to titrate the sample until a light pink color was maintained for 30 s and used to calculate the acidity value (Equation (5)).
AV ( mg / g ) = 5.611 × ( S B ) × N   S a m p l e   w t   ( g )  
where S = sample titration (mL); B = blank titration; N = normality of KOH–ethanol solution.

2.4. Peroxide Value (PV)

PV assays was carried out following the method by Kim et al. [18] with slight modifications. Briefly, 20 mL of acetic acid/chloroform (3:2 v/v) was used to dissolve a 1 g sample of oil, and then, a burette was used to add 1 mL of a saturated potassium iodide solution. Iodine was created after the reaction had occurred. Deionized water (50 mL) was added. Using a burette, iodine was titrated with 0.01 N sodium thiosulfate pentahydrate. The following Equation (6) was used to compute the peroxide value (expressed as milliequivalents of peroxide per kilogram of sample):
PV ( meq / kg ) = ( S B ) × N   S a m p l e   w t   ( g )   ×   1000
where S = sample titration (uL); B = blank titration; N = normality of Na2S2O3

2.5. Microbial Analysis

Total aerobic bacteria and total coliform in the sample were determined using the methods described by Lee et al. [25].

2.5.1. Total Aerobic Bacteria (TB)

Each 25 g sample was blended with 225 mL of sterile saline for 2 min using a stomacher (BagMixer® 400, Interscience Inc., St. Nom, France). Then, a series of 10-fold serial dilutions was prepared using sterile saline, and each dilution (1 mL) was inoculated on Petrifilm aerobic count plates (Petrifilm™ 3M, St. Paul, MN, USA) that were incubated at 30 °C for 48 h.

2.5.2. Total Coliform (TC)

To count total coliform, Petrifilm coliform count plates (Petrifilm™ 3M, St. Paul, MN, USA) were used, and the samples were incubated at 37 °C for 24 h. The identified colonies were counted and expressed as the log of the number of colony-forming units (CFUs)/g.

2.6. Experimental Design and Statistical Analysis

RSM was used to establish the optimization of the heat-drying conditions [26], and the optimization was performed using Minitab Statistical software Cloud-Based Version 20 (Minitab LLC., State College, PA, USA). A central composite design was employed to design the experimental data. The independent variables were set as blanching time (X1), drying temperature (X2), and drying time (X3) with a total 20 samples. Table 1 presents the values of α for the three independent variables, which were encoded into 5 levels, and the value of α was set to 1.68 [27]. Based on the results of preliminary experiments, the outline of experimental design with the coded and actual levels is presented in Table 2. The dependent variables (Y) were color parameters (WI, YI, BI, △E, and BD), TB and AV. These values were related to the coded variables by a quadratic polynomial using the equation below (7):
Y = β0 + β1X1 + β2X2 + β3X3 + β11X12 + β22X22 + β33X32 + β12X1X2 + β13X1X3 + β23X2X3
The coefficients of the polynomial model were represented by β0 (constant term), β1, β2 and β3 (linear effects), β11, β22 and β33 (quadratic effects), and β12, β13 and β23 (interaction effects). Analysis of variance (ANOVA) was used to determine the significance of the regression coefficient, and the developed regression models were validated using the conducted statistical analysis. The adequacy of the model was evaluated by the coefficient of multiple determination (R2). The terms in the polynomial were considered statistically different when p < 0.05.
A t-test was used to determine the significance of differences between optimal heat-dried and freeze-dried drone pupae. Statistical analyses were performed using SPSS 22.0 for Windows (SPSS Inc., Chicago, IL, USA).

3. Results and Discussion

For the sake of simplicity, some of the three-dimensional response surface plots are not depicted in Figure 1.

3.1. Model Fitting

The corresponding fitting of the explanatory models and the variation of the color parameters (WI, BI, YI, △E, and BD), TB, and AV were analyzed by the sum of squares of the sequential model. The estimated regression coefficients of the quadratic polynomial models for the response variables, along with the corresponding coefficients of determination (R2), are presented in Table 3. The coefficient of determination, R2, is the proportion of variation in the response attributed to the model rather than random error. It is recommended that an R2 value of at least 80% is necessary for a well-fitted model, and as the R2 value approaches unity, it indicates the empirical model’s suitability to the actual data. Conversely, a lower value of R2 indicates the unsuitability of the model to explain the relationship between variables [28]. Our results indicate that the R2 values for response variables were higher than 0.80. Specifically, the R2 values for WI, BI, YI, △E, BD, TB, and AV were found to be 0.8155, 0.9557, 0.9254, 0.9132, 0.9778, 0.9126, and 0.8693, respectively. These values suggest that the models accurately predicted the properties of these variables.

3.2. Color Parameters

The results in Table 3 indicate that the WI was statistically significant (p < 0.05) in the linear term of drying temperature (X2) and interaction term of blanching time and drying time (X1X3). The regression model fitted to the experimental results of the WI showed a correlation coefficient of R2 of 0.8155. The lack of fit (Table 4) was not significant, showing that these models are reliable predictors of responses. The negative coefficient of the linear term of drying temperature (X2) indicated that the WI decreased with an increase in drying temperature, which is to be expected by the dark color of the dried drone pupae. On the other hand, the positive coefficient of the interaction term (X1X3) of blanching time and drying time resulted in an increase in the WI. This is evident in Figure 1a; it was observed that the WI gradually increased with longer blanching and drying times. The results of this study are consistent with the findings of Cacchiarelli et al. [9], who reported higher WI values in larvae subjected to longer blanching times and higher blanching temperatures. For the BI (Table 3), it was observed that the significant (p < 0.001) effects were linear drying temperature (X2) and linear drying time (X3), interaction of blanching time and drying temperature (X1X2), quadratic blanching time (X1) and quadratic drying time (X3). Based on the sum of squares (Table 4), the importance of the independent variables on the BI could be ranked in the order of drying time (X3) and drying temperature (X2). Figure 1b describes the interaction between blanching time (X1) and drying temperature (X2). It was observed that the BI increased as the blanching time and drying temperature increased. Similarly, for the YI (Table 3), it was observed that the significant (p < 0.001) effects were linear drying temperature (X2) and linear drying time (X3), interaction of blanching time and drying temperature (X1X2), quadratic blanching time (X1) and quadratic drying time (X3). For blanching, although the linear effect of blanching time on the BI and YI was not significant, the coefficient of quadratic effects was significant. Based on the sum of squares (Table 4), the importance of the independent variables on the YI could be ranked in the order of drying time (X3) and drying temperature (X2).
In Figure 1d, samples dried under different drying conditions showed an increase in color change, with a significant increase in the △E value observed with increasing drying time. Based on the sum of squares, the importance of the independent variables could be ranked in the order of drying time (X3) and drying temperature (X2). Figure 1e,f show that the significant (p < 0.001; Table 3) interactions in the response variable for BD were all linear, quadratic, and interaction terms except for the interaction of blanching time and drying temperature (X1X2). The BD is defined as the degree of purity of the brown color and is one of the most commonly used indicators of browning in various food products [29]. These data are consistent with the WI, where higher drying temperatures trigger the Maillard reaction and caramelization, resulting in a darker hue and increased intensity of color in the drone pupae. After blanching, BD gradually increased with increasing drying time. This becomes more evident with the three-dimensional plot, where the linear blanching time and drying temperature term was significant (p < 0.001; Table 3) with a negative effect. On the other hand, BD gradually increases as the drying time increases. Therefore, the longer drying time causes the increase in BD, which is to be expected by the dark color of the dried drone pupae.

3.3. Total Aerobic Bacteria

The regression equation for expansion as TB is shown in Table 3. It was observed that the drying temperature and time had very significant positive linear effects (p < 0.001). Based on the sum of squares (Table 4), the importance of the independent variables could be ranked in the order of drying temperature (X2) and drying time (X3). It can be observed that there is an increasing trend in TB with drying temperature and time. During the drying of drone pupae, an increase in microbial counts is observed, which is expected due to the rich nutrients in them.

3.4. Acid Value

The regression equation for expansion as AV is shown in Table 3. It was observed that there is a decrease in AD (mg/kg) as the blanching time increases with the evaluation of drying temperature and drying time. However, when the changes in AD values were measured according to the drying time, it was observed that the AD value increased as the drying temperature increased. These findings are supported in Table 3, where significant effects (p < 0.05) were observed for the terms of linear blanching time (X1) (negative effect), drying temperature (X2) (positive effect), drying time (X3) (positive effect) and quadratic blanching time (X1) (negative effect). Therefore, it can be expected that the acid value of dried drone pupae would increase as a result of the drying treatment after blanching, since this treatment affects the oxidation of the product.

3.5. Optimal Drying Conditions

The optimal parameters for the development of drying conditions of drone pupae to improve pre-processing quality were found to be the minimum values of BI, YI, △E, BD, TB, and AV, as well as the maximum value of WI. The optimal drying conditions were determined to be a blanching time of 58 s, a drying temperature of 56.7 °C, and a drying time of 298 min. The predicted and actual values for the optimal drying conditions are given in Table 5. The difference of process parameters between the actual data and the statistically predicted data in the processed drone pupae with optimized drying conditions showed a minimum of 1.41% difference and a maximum of 29.38% difference. A 95% confidence level t-test (IBM SPSS Statistics 22) was used to compare the predicted and actual values for dried drone pupae in optimum drying condition circumstances, and no significant difference was detected.

3.6. Physicochemical and Microbiological Properties of Freeze Dried and Heat-Dried Drone Pupae

3.6.1. Acid Value and Peroxide Value of Freeze-Fried and Heat-Dried Drone Pupae

To compare the AV and PV of freeze-dried drone pupae, measurements were taken of the heat-dried drone pupae optimized in this study using a commonly used preprocessing method for edible insects (Table 6). The AV did not show a significant difference between the two samples. However, the PV was higher in the heat-dried sample. The PV, as an indicator of the initial stage of oxidation spoilage, is represented by the amount of peroxide, which is a primary oxidation product formed during the initial oxidation process of lipid components. The PV of the freeze-dried sample was 5.74 ± 1.90 meq/kg, while that of the heat-dried sample was significantly higher at 10.36 ± 1.02 meq/kg. The processing of raw materials can significantly affect quality, nutritional composition, digestibility, and palatability in food manufacturing. This is particularly true for animal-based raw materials such as edible insects, which require a slaughtering stage and subsequent microbiological and chemical stabilization by drying to reduce water activity [14]. Defatting is a common step in the use of edible insects to produce high-protein products, which is the primary purpose for using them. Lipids can dilute protein concentration, hinder manufacturing technology, and affect palatability, as well as the productive yields and lipid composition of animals in the case of feed. The three main processes involved in processing insect larvae, which are slaughter, drying, and defatting, can significantly affect lipid oxidation in fats and meals. Blanching and freezing are the most popular methods for slaughtering insects. After testing different killing methods, blanching was preferred by Larouche et al. [30] for all quality parameters of black soldier fly (Hermetia illucens) larvae meals, including low lipid oxidation. Thermal methods, pro-oxidant conditions, and the degradation of minor antioxidant compounds due to processing can contribute to lipid oxidation. Both thermal methods (such as oven or microwave) and non-thermal methods (such as freeze-drying) have been tested for drying insect larvae, but conclusive results regarding the impact of drying methods on lipid oxidation have not been reached, which is likely due to dependence on the slaughtering method used [31,32]. In the present study, there was no significant difference in acid value between freeze-dried and heat-dried drone pupae.

3.6.2. Color Indices of Freeze-Fried and Heat-Dried Drone Pupae

The results of the investigation into the color of freeze-dried and heat-dried drone pupae showed significant lowering of L and b values (p < 0.05) in the heat-dried drone pupae, which is consistent with a study conducted by Brishti et al. [33]. They suggested that the formation of brown pigments, such as melanoidins, through the Maillard reaction under heated conditions during oven drying may be the reason for the observed result. This finding is supported by the lowest amount of lysine content in oven dried, as lysine is actively involved in the Maillard reaction, and lower lysine content confirms the occurrence of the Maillard reaction. Color is an essential quality parameter associated with product or raw material quality attributes such as freshness, sensory and nutritional properties, and the presence of visual and non-visual defects. Additionally, color can indirectly control the Maillard reaction [34]. In the present study, it was observed that the color of the drone pupae became darker brown depending on the blanching time and drying conditions. The results of this study are similar to those found in other research. Cacchiarelli et al. [9] reported that insects that have undergone blanching treatment experience a decrease in water-holding capacity. They reported that when color parameters were measured in different pH solutions after boiling treatments (60 °C for 5 min and 90 °C for 1 min), WI and BI decreased more slowly than in the untreated control, indicating a browning effect. They also reported that the main effect on color differences (ΔE), which was calculated as the difference in color between time 0 and subsequent times, was due to the blanching treatment, which was followed by the effect of pH. Furthermore, these findings are consistent with the report of Saucier et al. [35] of a significant effect of blanching on black soldier fly larvae. They reported that pre-treatment of larvae with puncturing, blanching, and scalding in boiling water resulted in reduced drying times, which was presumably due to the impact on the wax-coated cuticle that shields the larvae from desiccation. Khatun et al. [36] investigated the color changes of freeze-dried, oven-dried, and blanched house crickets, and it was found that freeze-drying resulted in significantly higher levels of lightness than blanching and a lower browning index for cricket species. Another study observed the color of four insect species treated with the blanching method during low-temperature storage and reported that while the change in the L value was not significant over the storage period, there was a tendency for the a and b values to decrease due to the destruction of pigments and the formation of browning reaction products [37]. Kim et al. [38] conducted a study to measure the color of cricket powder that was freeze-dried and then oven-dried at various storage temperatures. The study found that at 25 °C, the L value increased slightly, but it decreased at 35 °C and 40 °C, with a trend of decreasing as the storage temperature increased over a long period of 6 months. Additionally, the a and b values decreased at all storage temperatures with a significantly greater reduction observed at higher temperatures compared to lower temperatures. Several studies have investigated the impact of freezing and blanching methods on the protein and lipid quality, as well as the color stability, of black soldier fly larvae. Leni et al. [39] reported that blanching, a method of killing insects, inhibits the browning reaction and other enzymatic changes that occur during slow killing by freezing. This leads to an increase in the extractability of proteins in aqueous solutions, prevents essential amino acid loss, and improves enzymatic digestibility. Caligiani et al. [40] reported that prepupae killed by freezing showed a significant reduction in acylglycerols during storage and the release of free fatty acids, which was likely due to the activation of lipases. In contrast, prepupae killed by blanching have a stable lipid fraction consisting mainly of triacylglycerols. Thus, the killing method has a significant impact on black soldier fly oil composition and potential applications. Larouche et al. [30] conducted a study to optimize larval killing methods and evaluate their impact on the nutritional and microbiological quality of black soldier fly larvae. Ten different methods were tested, resulting in varying the coloration of the freeze-dried and granulated larvae. Asphyxiation and cold-killing methods produced similar colors, while mechanical disruption and heating methods resulted in significantly different colors. Desiccation resulted in larvae with significantly higher color intensity than blanched or frozen larvae. High-pressure processing resulted in a product with significantly higher lightness and color intensity than grinding. Asphyxiation resulted in larvae with higher lightness and hue angle, resulting in a color closer to yellow compared to heat methods.

3.6.3. Microbial Analysis of Freeze-Fried and Heat-Dried Drone Pupae

The total aerobic bacteria count was significantly higher in heat-dried drone pupae (2.99 Log CFUs/g) compared to freeze-dried (1.99 Log CFUs/g). Regarding total coliform, no colonies were detected after any treatment. In mulberry silkworm powder, the total microorganism count was 6.54 log CFUs/g. It was reported that frozen silkworms (Bombyx Mori), bamboo caterpillars (Omphisa Fuscidentalis), and field crickets (Gryllus Bimaculatus) had total aerobic counts of 7.9–8.3 log CFUs/g, yeast and mold counts of 5.5–6.2 log CFUs/g, and Enterobacteriaceae counts of 4.7–6.5 log CFUs/g [31,41]. Insects bred in contact with soil or sawdust are prone to contamination with various microorganisms, resulting in generally high microbial counts in insects. However, since the drone pupae grow in beehives and consume honey and royal jelly, their internal microbial count is lower compared to other edible insects.

4. Conclusions

To optimize the drying conditions of drone pupae, the response surface methodology (RSM) was used, and optimization was performed using Minitab software Cloud-Based Version 20. Our results showed that the R2 values for these response variables were higher than 0.80, indicating that the regression models were suitable to explain the behavior. The optimal drying conditions were determined to be a blanching time of 58 s, a drying temperature of 56.7 °C, and a drying time of 298 min. The difference in process parameters between the actual data and the statistically predicted data in the processed drone pupae with optimized drying conditions showed almost identical results, indicating that the pre-treatment conditions for edible insects using RSM are a suitable method. When comparing the physicochemical and microbiological characteristics of drone pupae processed by freeze-drying and drone pupae processed by optimized heat treatment, our results showed that the a value in color and acidity were not significantly different between the two groups. Our study suggests that the model we obtained can be applied to the large-scale production of drying conditions for further use in the pharmacy/food industries. Optimized conditions may be used in research involving test foods for the starting procedure of insect pre-processing for analysis of the samples. RSM is an effective tool for accounting for surface factors for an outcome and was used here to save labor and cost.

Author Contributions

Conceptualization, I.N.; Methodology, S.B.; Software, S.B.; Validation, A.S.M.; Formal analysis, A.S.M.; Investigation, A.S.M.; Data curation, S.B.; Writing—original draft, S.B.; Writing—review and editing, I.N.; Project administration, I.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Rural Development Administration of Korea (PJ016177).

Data Availability Statement

The data are available from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tilman, D.; Balzer, C.; Hill, J.; Befort, B.L. Global food demand and the sustainable intensification of agriculture. Proc. Natl. Acad. Sci. USA 2011, 108, 20260–20264. [Google Scholar] [CrossRef]
  2. Arnold, V.H. Potential of Insects as Food and Feed in Assuring Food Security. Annu. Rev. Entomol. 2013, 58, 563–583. [Google Scholar]
  3. Arnold, V.H. Edible insects are the future. Proc. Nutr. Soc. 2013, 75, 294–305. [Google Scholar]
  4. Van Huis, A.; Van Itterbeeck, J.; Klunder, H.; Mertens, E.; Halloran, A.; Muir, G. Edible Insects: Future Prospects for Food and Feed Security; FAO: Rome, Italy, 2013; p. 171. [Google Scholar]
  5. Janssen, R.H.; Vincken, J.-P.; Arts, N.J.; Fogliano, V.; Lakemond, C.M. Effect of endogenous phenoloxidase on protein solubility and digestibility after processing of Tenebrio molitor, Alphitobius diaperinus and Hermetia illucens. Food Res. Int. 2018, 121, 684–690. [Google Scholar] [CrossRef] [PubMed]
  6. David-Birman, T.; Raften, G.; Lesmes, U. Effects of thermal treatments on the colloidal properties, antioxidant capacity and in-vitro proteolytic degradation of cricket flour. Food Hydrocoll. 2018, 79, 48–54. [Google Scholar] [CrossRef]
  7. Janssen, R.H.; Canelli, G.; Sanders, M.G.; Bakx, E.J.; Lakemond, C.M.M.; Fogliano, V.; Vincken, J.-P. Iron-polyphenol complexes cause blackening upon grinding Hermetia illucens (black soldier fly) larvae. Sci. Rep. 2019, 9, 1–11. [Google Scholar] [CrossRef]
  8. Mozhui, L.; Kakati, L.; Kiewhuo, P.; Changkija, S. Traditional Knowledge of the Utilization of Edible Insects in Nagaland, North-East India. Foods 2020, 9, 852. [Google Scholar] [CrossRef]
  9. Cacchiarelli, C.; Fratini, F.; Puccini, M.; Vitolo, S.; Paci, G. Effects of different blanching treatments on colour and microbio-logical profile of Tenebrio molitor and Zophobas morio larvae. LWT-Food Sci. Technol. 2022, 157, 13112. [Google Scholar] [CrossRef]
  10. Andrea, M.L. Processing insects for use in the food and feed industry. Insect Sci. 2021, 48, 32–36. [Google Scholar]
  11. Mancini, S.; Fratini, F.; Tuccinardi, T.; Turchi, B.; Nuvoloni, R.; Paci, G. Effects of different blanching treatments on microbiological profile and quality of the mealworm (Tenebrio molitor). J. Insects Food Feed 2019, 5, 225–234. [Google Scholar] [CrossRef]
  12. Hider, R.C. Honeybee venom: A rich source of pharmacologically active peptides. Endeavour 1988, 12, 60–65. [Google Scholar] [CrossRef] [PubMed]
  13. Banerjee, P.; Sohoo, K.N.; Biswas, T.K.; Basu, S.K.; Chatterjee, J.; Hui, A.K.; Chakraborty, N.C.; Debnath, P.K. Bees make medicine for mankind. Ind. J. Trad. Knowl. 2003, 2, 22–26. [Google Scholar]
  14. Al-Waili, N.S. Mixture of Honey, Beeswax and Olive Oil Inhibits Growth of Staphylococcus aureus and Candida albicans. Arch. Med. Res. 2005, 36, 10–13. [Google Scholar] [CrossRef] [PubMed]
  15. Thomas, S.; Andrews, A.M.; Hay, N.P.; Bourgoise, S. The anti-microbial activity of maggot secretions: Results of a preliminary study. J. Tissue Viability 1999, 9, 127–132. [Google Scholar] [CrossRef]
  16. Ghosh, S.; Sohn, H.-Y.; Pyo, S.-J.; Jensen, A.B.; Meyer-Rochow, V.B.; Jung, C. Nutritional Composition of Apis mellifera Drones from Korea and Denmark as a Potential Sustainable Alternative Food Source: Comparison Between Developmental Stages. Foods 2020, 9, 389. [Google Scholar] [CrossRef] [Green Version]
  17. Kim, S.G.; Woo, S.O.; Bang, K.W.; Jang, H.R.; Han, S.M. Chemical composition of drone pupa of Apis mellifera and its nu-tritional evaluation. J. Apic. 2018, 33, 103086076. [Google Scholar]
  18. Kim, H.Y.; Woo, S.O.; Kim, S.G.; Bang, K.W.; Choi, H.M.; Moon, H.J.; Han, S.M. Analysis of Oxidative Stability in Drone Pupae (Apis mellifera L.). J. Apic. 2019, 34, 63–66. [Google Scholar] [CrossRef]
  19. Kim, H.Y.; Woo, S.O.; Kim, S.G.; Choi, H.M.; Moon, H.J.; Han, S.M. Antioxidant and Antihyperglycemic Effects of Honeybee Drone Pupae (Apis mellifera L.) Extracts. J. Apic. 2020, 35, 33–39. [Google Scholar] [CrossRef]
  20. Kim, J.-E.; Kim, D.-I.; Koo, H.-Y.; Kim, H.-J.; Kim, S.-Y.; Lee, Y.-B.; Moon, J.-H.; Choi, Y.-S. Evaluation of Honey Bee (Apis mellifera L.) Drone Pupa Extracts on the Improvement of Hair Loss. J. Apic. 2020, 35, 179–188. [Google Scholar] [CrossRef]
  21. Kim, H.Y.; Woo, S.O.; Kim, S.G.; Bang, K.W.; Choi, H.M.; Moon, H.J.; Han, S.M. Anti-inflammatory Activities of Drone Pupae (Apis mellifera L.) in Macrophages. J. Apic. 2019, 34, 255–259. [Google Scholar] [CrossRef]
  22. Kim, J.-E.; Kim, D.-I.; Kim, H.-J.; Kim, S.-Y.; Lee, Y.-B.; Moon, J.-H.; Park, H.-G.; Choi, Y.-S. Characteristics of Hydrolysis of Protein in Drone Pupa (Apis mellifera L.). J. Apic. 2020, 35, 169–177. [Google Scholar] [CrossRef]
  23. Téllez-Morales, J.A.; Hernández-Santos, B.; Navarro-Cortez, R.O.; Rodríguez-Miranda, J. Impact of the addition of cricket flour (Sphenarium purpurascens) on the physicochemical properties, optimization and extrusion conditions of extruded nixtamalized corn flour. Appl. Food Res. 2022, 2, 100149. [Google Scholar] [CrossRef]
  24. Chang, M.-J.; Han, M.-R.; Kim, M.-H. Effects of Salt Addition in Sugar Based Osmotic Dehydration on Mass Transfer and Browning Reaction of Carrots. Prev. Nutr. Food Sci. 2003, 8, 230–234. [Google Scholar] [CrossRef]
  25. Lee, S.; Choi, Y.S.; Jo, K.; Kim, T.K.; Yong, H.I.; Jung, S. Quality characteristics and protein digestibility of protaetia brevi-tarsis larvae. J. Anim. Sci. Technol. 2020, 62, 741–752. [Google Scholar] [CrossRef] [PubMed]
  26. Yolmeh, M.; Jafari, S.M. Applications of Response Surface Methodology in the Food Industry Processes. Food Bioprocess Technol. 2017, 10, 413–433. [Google Scholar] [CrossRef]
  27. Deaton, M.L.; Khuri, A.I.; Cornell, J.A. Response Surfaces: Designs and Analysis. J. Am. Stat. Assoc. 1989, 84, 332. [Google Scholar] [CrossRef]
  28. Koocheki, A.; Taherian, A.R.; Razavi, S.M.A.; Bostan, A. Response surface methodology for optimization of extraction yield, viscosity, hue and emulsion stability of mucilage extracted from Lepidium perfoliatum seeds. Food Hydrocoll. 2009, 23, 2369–2379. [Google Scholar] [CrossRef]
  29. Pathare, P.B.; Opara, U.L.; Al-Said, F.A.-J. Colour Measurement and Analysis in Fresh and Processed Foods: A Review. Food Bioprocess Technol. 2013, 6, 36–60. [Google Scholar] [CrossRef]
  30. Larouche, J.; Deschamps, M.H.; Saucier, L.; Lebeuf, Y.; Doyen, A.; Vandenberg, G.W. Effects of killing methods in lipid ox-idation, colour and microbial load of Black Soldier Fly (Hermetia illucens) larvae. Animals 2019, 9, 182. [Google Scholar] [CrossRef] [Green Version]
  31. Kurdi, P.; Chaowiwat, P.; Weston, J.; Hansawasdi, C. Studies on microbial quality, protein yield, and antioxidant of some frozen edible insects. Int. J. Food Sci. 2021, 5580976, 1–7. [Google Scholar] [CrossRef]
  32. Janssen, R.H.; Lakemond, C.M.M.; Fogliano, V.; Renzone, G.; Scaloni, A.; Vincken, J.P. Involvement of phenoloxidase in browning during grinding of Tenebrio molitor larvae. PLoS ONE 2017, 12, e0189685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Brishti, F.H.; Chay, S.Y.; Muhammad, K.; Ismail-Fitry, M.R.; Zarei, M.; Karthikeyan, S.; Saari, N. Effects of drying techniques on the physicochemical, functional, thermal, structural and rheological properties of mung bean (Vigna radiata) protein isolate powder. Food Res. Int. 2020, 138, 109783. [Google Scholar] [CrossRef]
  34. Cristina, R. Hot air and freeze-drying of high-value foods: A review. J. Food Eng. 2001, 49, 311–319. [Google Scholar]
  35. Saucier, L.; M’ballou, C.; Ratti, C.; Deschamps, M.-H.; Lebeuf, Y.; Vandenberg, G. Comparison of black soldier fly larvae pre-treatments and drying techniques on the microbial load and physico-chemical characteristics. J. Insects Food Feed. 2021, 8, 45–64. [Google Scholar] [CrossRef]
  36. Khatun, H.; Claes, J.; Smets, R.; De Winne, A.; Akhtaruzzaman, M.; Van Der Borght, M. Characterization of freeze-dried, oven-dried and blanched house crickets (Acheta domesticus) and Jamaican field crickets (Gryllus assimilis) by means of their physicochemical properties and volatile compounds. Eur. Food Res. Technol. 2021, 247, 1291–1305. [Google Scholar] [CrossRef]
  37. Son, Y.-J.; Ahn, W.; Kim, S.-H.; Park, H.-N.; Choi, S.-Y.; Lee, D.-G.; Kim, A.-N.; Hwang, I.-K. Study on the Oxidative and Microbial Stabilities of Four Edible Insects during Cold Storage after Sacrificing with Blanching Methods. Korean J. Food Nutr. 2016, 29, 849–859. [Google Scholar] [CrossRef] [Green Version]
  38. Kim, D.H.; Kim, E.M.; Chang, Y.J.; Ahn, M.Y.; Lee, Y.H.; Park, J.J.; Lim, J.H. Determination of the shelf life of cricket powder and effects of storage on its quality characteristics. Korean J. Food Preserv. 2016, 23, 211–217. [Google Scholar] [CrossRef] [Green Version]
  39. Leni, G.; Caligiani, A.; Sforza, S. Killing method affects the browning and the quality of the protein fraction of Black Soldier Fly (Hermetia illucens) prepupae: A metabolomics and proteomic insight. Food Res. Int. 2018, 115, 116–125. [Google Scholar] [CrossRef]
  40. Caligiani, A.; Marseglia, A.; Sorci, A.; Bonzanini, F.; Lolli, V.; Maistrello, L.; Sforza, S. Influence of the killing method of the black soldier fly on its lipid composition. Food Res. Int. 2019, 116, 276–282. [Google Scholar] [CrossRef]
  41. Anuduang, A.; Loo, Y.Y.; Jomduang, S.; Lim, S.J.; Mustapha, W.A.W. Effect of Thermal Processing on Physico-Chemical and Antioxidant Properties in Mulberry Silkworm (Bombyx mori L.) Powder. Foods 2020, 9, 871. [Google Scholar] [CrossRef]
Figure 1. The three-dimensional response surface plots showing interactive effects of blanching time (X1), drying temperature (X2), and drying time (X3) on WI (a), BI (b), YI (c), △E (d), BD (e,f). Each net surface represents the response surface predicted with the quadratic model as a function of each variable and described by the equations given in Table 3. *: significant at p < 0.05. **: significant at p < 0.01.
Figure 1. The three-dimensional response surface plots showing interactive effects of blanching time (X1), drying temperature (X2), and drying time (X3) on WI (a), BI (b), YI (c), △E (d), BD (e,f). Each net surface represents the response surface predicted with the quadratic model as a function of each variable and described by the equations given in Table 3. *: significant at p < 0.05. **: significant at p < 0.01.
Foods 12 03062 g001
Table 1. Experimental range and levels of the independent variables.
Table 1. Experimental range and levels of the independent variables.
Independent VariableXiLevel
−1.68−1011.68
Blanching Time (s)X15380120160182
Drying Temperature (°C)X241.645505558.4
Drying Time (min)X3266300350400434
Table 2. Central composite experimental design for drying treatment with coded and actual levels of the independent variable.
Table 2. Central composite experimental design for drying treatment with coded and actual levels of the independent variable.
RunPoint TypeCoded VariablesActual Variables
X1X2X3X1X2X3
1Factorial−1−1−18045300
2Factorial1−1−116055300
3Factorial−11−18045300
4Factorial11−116055300
5Factorial−1−118045400
6Factorial1−1116055400
7Factorial−1118045400
8Factorial11116055400
9Axial−1.68005850350
10Axial1.680018250350
11Axial0−1.68012041.6350
12Axial01.68012058.4350
13Axial00−1.6812050266
14Axial001.6812050434
15Central00012050350
16Central00012050350
17Central00012050350
18Central00012050350
19Central00012050350
20Central00012050350
Table 3. Polynomial regression equations for color parameters (WI, BI, YI, △E, and BD), TB, and AV by independent variables.
Table 3. Polynomial regression equations for color parameters (WI, BI, YI, △E, and BD), TB, and AV by independent variables.
Response Surface ModelR2p Value
WI=67.842 − 0.305 X2 − 0.399 X32 + 0.521 X1X30.81550.010
BI=24.187 − 2.792 X2 − 3.338 X3 + 1.456 X12 − 1.132 X32 + 1.798 X1X20.9557<0.001
YI=28.316 − 2.447 X2 − 3.250 X3 + 1.623 X12 − 1.363 X32 + 1.898 X1X20.9254<0.001
△E=27.62 + 8.94 X2 + 12.52 X3 − 4.22 X12 + 4.71 X22 + 7.02 X32 − 6.74 X1X20.9132<0.001
BD=0.21796 − 0.00782 X1 − 0.03660 X2 + 0.04627 X3 − 0.01752 X12 − 0.01155 X22 + 0.00985 X32− 0.01046 X1X3 − 0.04085 X2X30.9778<0.001
TB=3.345 + 0.741 X2 + 0.563 X3 + 0.3209 X120.9126<0.001
AV=0.12702 − 0.01017 X1 + 0.00643 X2 + 0.01034 X3 − 0.00629 X120.8693<0.001
Table 4. Analysis of variance results and determination coefficients for RSM models.
Table 4. Analysis of variance results and determination coefficients for RSM models.
SourceWIBIYI△E
Sum of SquaresF Valuep ValueSum of SquaresF Valuep ValueSum of SquaresF Valuep ValueSum of SquaresF Valuep Value
Model10.87814.910.010351.42423.99<0.001342.85813.79<0.0015280.8411.69<0.001
Linear2.24193.040.080264.15854.09<0.001228.51627.570.0013230.0521.46<0.001
X10.91113.700.0835.7803.550.0892.7230.990.3440.470.010.925
X21.27155.170.046106.33165.32<0.00181.69229.57<0.0011090.8921.740.001
X30.05930.240.634152.04793.40<0.001144.10152.17<0.0012138.6942.63<0.001
Crossed4.41755.980.01331.4426.440.01136.6014.420.0321349.238.960.003
X1X21.13174.600.05825.85315.880.00328.8210.430.009256.055.10.047
X1X32.17288.830.0142.3011.410.2625.2941.920.196318.386.350.030
X2X31.11294.520.0593.2882.020.1862.4870.900.365708.7114.130.004
Quadratic4.21875.710.01555.82411.430.00177.7419.380.003701.564.660.028
X120.70612.870.12130.43518.700.00237.84813.700.004363.487.240.023
X220.81093.300.1001.4780.910.3635.7242.070.181135.52.70.131
X322.29129.310.01218.42311.320.00726.6889.660.011202.584.040.072
Lack-of-fit2.08335.020.05513.4854.830.05523.7456.120.034455.149.770.013
Pure error0.3777 2.794 3.879 46.59
Residual2.4609 16.279 27.624 501.73
Total13.3391 367.703 370.482 5782.57
SourceBDTBAV
Sum of SquaresF Valuep ValueSum of SquaresF Valuep ValueSum of SquaresF Valuep Value
Model0.0710049.04<0.00114.909011.60<0.0010.004457.390.002
Linear0.04832100.14<0.00111.967427.92<0.0010.0034417.12<0.001
X10.000845.190.0460.15041.0500.3290.0014121.09<0.001
X20.01828113.65<0.0017.492452.45<0.0010.000578.440.016
X30.02921181.58<0.0014.324630.27<0.0010.0014621.82<0.001
Crossed0.0145930.24<0.0010.35610.830.5070.000341.690.232
X1X20.000372.300.1600.00230.020.9020.000142.140.174
X1X30.000885.440.0420.14301.000.3410.000050.710.418
X2X30.0133582.97<0.0010.21071.480.2520.000152.210.168
Quadratic0.0080816.75<0.0012.58566.030.0130.000673.360.063
X120.0044127.42<0.0011.479610.360.0090.000578.490.015
X220.0019211.920.0060.50113.510.0910.000142.020.185
X320.001408.670.0150.46733.270.1010.000000.020.892
Lack-of-fit0.0013354.870.0541.01662.470.1720.0002560.620.694
re error0.000274 0.4119 0.000413
Residual0.001609 1.4286 0.000669
Total0.072607 16.3376 0.005116
Table 5. Predicted and experimental values of the responses at optimum conditions.
Table 5. Predicted and experimental values of the responses at optimum conditions.
Process ParameterTargetPredicted DataActual Data% Difference
WIMaximize69.502 a66.67 ± 0.56 a4.16 ± 0.84
BIMinimize22.898 a21.33 ± 0.43 a7.09 ± 2.04
YIMinimize25.909 a26.27 ± 0.48 a1.41 ± 1.79
△EMinimize25.55 a31.27 ± 9.79 a28.43 ± 15.46
BDMinimize0.135 a0.13 ± 0.01 a4.54 ± 3.59
TB (log CFU/g)Minimize4.193 a3.12 ± 0.12 a29.38 ± 3.61
AV (mg/g)Minimize4.114 a4.33 ± 0.18 a4.98 ± 4.19
Actual data and % difference indicate mean ± SD of triplicate measurements. Means within a row with same letters are not significantly different (p < 0.05).
Table 6. The comparison of the physicochemical and microbiological properties of the drone pupae treated with optimal heat drying and freeze drying.
Table 6. The comparison of the physicochemical and microbiological properties of the drone pupae treated with optimal heat drying and freeze drying.
ParameterFreeze DriedHeat Driedp Value
L value69.03 ± 0.49 a63.23 ±0.60 b0.008
a value2.60 ± 0.22 a1.66 ±0.17 a0.093
b value20.27 ± 0.82 a12.73 ±0.29 b0.003
AV (mg/g)3.83 ± 0.06 a4.20 ± 0.26 a0.062
PV (meq/kg)5.74 ± 1.90 b10.36 ± 1.02 a0.033
TB (log CFU/g)1.99 ± 0.35 b2.99 ± 0.09 a0.032
All values are mean ± SD of triplicate measurements. Means within a row with different letters are significantly different (p < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Baek, S.; Mae, A.S.; Nam, I. Optimization of the Heat-Drying Conditions of Drone Pupae by Response Surface Methodology (RSM). Foods 2023, 12, 3062. https://doi.org/10.3390/foods12163062

AMA Style

Baek S, Mae AS, Nam I. Optimization of the Heat-Drying Conditions of Drone Pupae by Response Surface Methodology (RSM). Foods. 2023; 12(16):3062. https://doi.org/10.3390/foods12163062

Chicago/Turabian Style

Baek, SeungHee, Agapito Sheryl Mae, and InSik Nam. 2023. "Optimization of the Heat-Drying Conditions of Drone Pupae by Response Surface Methodology (RSM)" Foods 12, no. 16: 3062. https://doi.org/10.3390/foods12163062

APA Style

Baek, S., Mae, A. S., & Nam, I. (2023). Optimization of the Heat-Drying Conditions of Drone Pupae by Response Surface Methodology (RSM). Foods, 12(16), 3062. https://doi.org/10.3390/foods12163062

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