Food Science and Preservation
The Korean Society of Food Preservation
Research Article

Effects of pre-fasting diets on nutritional composition and volatile profiles of Protaetia brevitarsis larvae

Min Ji Park1https://orcid.org/0009-0006-5065-8742, Ji-Youn Kim2https://orcid.org/0009-0002-8371-4457, Ji Yeong Park2https://orcid.org/0000-0002-9927-3513, Sun Young Kim1,*https://orcid.org/0000-0002-9953-0110
1Industrial Entomology Division, National Institute of Agricultural Sciences, Rural Development Administration, Wanju 55365, Korea
2Fermented & Processed Food Research Division, National Institute of Crop and Food Science, Rural Development Administration, Wanju 55365, Korea
*Corresponding author Sun Young Kim, Tel: +82-63-238-2978, E-mail: carp0120@korea.kr

Citation: Park MJ, Kim JY, Park JY, Kim SY. Effects of pre-fasting diets on nutritional composition and volatile profiles of Protaetia brevitarsis larvae. Food Sci. Preserv., 33(3), 442-453 (2026)

Copyright © The Korean Society of Food Preservation. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Mar 31, 2026; Revised: May 19, 2026; Accepted: Jun 01, 2026

Published Online: Jun 30, 2026

Abstract

In this study, we evaluated whether formulated depuration diets can modulate the nutritional quality, chemical safety, and volatile profile of Protaetia brevitarsis larvae. Late third-instar larvae were provided with glutinous rice powder (GRP) or one of three formulated diets (FD1-FD3) based on wheat bran, eggshell powder, and brewer’s spent grain for 48 h, followed by 48 h of fasting. Freeze-dried larval powders were analyzed for proximate composition, dietary fiber, fatty acids, minerals, and heavy metals, and volatile profiles were assessed using an electronic nose. Proximate composition, dietary fiber, and carbohydrate showed numerical treatment-dependent variation but no significant differences among groups by Tukey’s HSD test. Oleic acid was the predominant fatty acid, with FD3 forming the highest Tukey group. Calcium and iron were significantly higher in FD1-FD3. Pb was slightly but significantly higher in FD1-FD3 than in GRP, whereas Cd and As did not differ significantly; all measured heavy metals remained low and were below the Korean Food Code limits. Electronic-nose analysis revealed clear diet-dependent differences in volatile patterns: GRP was characterized mainly by butanal, 2,3-pentanedione, heptanal, and butane-2,3-dione, whereas FD2 was more closely associated with propyl propanoate, 2-heptanone, and anethole. These results indicate that diet formulation during depuration can serve as a practical pre-processing strategy to modulate selected nutritional and odor-related quality attributes of edible P. brevitarsis larvae without appreciable heavy-metal contamination.

Keywords: Protaetia brevitarsis; pre-fasting diet; nutritional composition; heavy metals; volatile profile

1. Introduction

Edible insects are increasingly explored as ingredients in food systems, yet their broader adoption depends on robust evidence of safety and consistent quality. Risk assessments and supply-chain analyses emphasize that insect-based materials may pose biological and chemical hazards originating from production substrates and processing conditions, and that standardization is essential to reduce variability and strengthen consumer trust (Fuso et al., 2024; Liceaga et al., 2022; Traynor et al., 2024). For food-oriented applications, pre-processing steps that influence hygiene, chemical safety, and product quality warrant systematic evaluation (Acosta-Estrada et al., 2021; Fraqueza and Patarata, 2017; Kooh et al., 2020).

Larvae of the white-spotted flower chafer, Protaetia brevitarsis, have been investigated as an insect resource for food applications in Korea (Kim et al., 2024; Park and Kim, 2025; Park et al., 2024; Park et al., 2025b). Previous research has proposed pre-treatment conditions for P. brevitarsis larvae intended for use as food ingredients, highlighting the importance of defining processing workflows before further conversion into powders or other forms (Kim et al., 2022; Kim et al., 2025). Despite this progress, key pre-processing control points that can affect both yield and product quality require additional evidence, particularly under practical farming and processing contexts.

A central pre-processing step for many edible insects is gut evacuation, often termed depuration or fasting, which is applied to reduce intestinal contents prior to downstream processing (Ahn et al., 2015; Chung et al., 2013; Mancini et al., 2020). Because insects are frequently processed and consumed as whole-body materials, residual gut contents can contribute to hygienic concerns and sensory issues (Baigts-Allende and Stathopoulos, 2023; Cappelli et al., 2020; Kooh et al., 2019; Zhang et al., 2024). Feed withdrawal has been used to promote fecal clearance in farmed insects; for example, in the edible cricket Gryllus bimaculatus, fecal output decreased substantially after a short fasting period, supporting depuration as a practical measure to reduce gut contents (Ahn et al., 2015; Kwon et al., 2013; Park et al., 2025a). At the same time, prolonged fasting can impose physiological stress and reduce harvest yield through weight loss or mortality, indicating a trade-off between gut clearance and production efficiency (Egonyu et al., 2025; Han et al., 2024). Therefore, depuration strategies that preserve gut-evacuation benefits while minimizing production losses are relevant to food-grade insect processing.

Diet and processing conditions are well-recognized sources of variability in insect composition (Cacchiarelli et al., 2022; Megido et al., 2017; Selaledi and Mabelebele, 2021). Reviews on edible insects have shown that macronutrient composition, lipid composition, and overall acceptability can vary with species, developmental stage, diet, and processing method (Kouřimská and Adámková, 2016; Meyer-Rochow et al., 2021; Oonincx and Finke, 2021). For P. brevitarsis, changes in the type of feed provided prior to processing have been associated with measurable shifts in nutritional composition, indicating that the pre-harvest period can influence final material characteristics (Ham et al., 2021; Noh et al., 2015). From a safety perspective, chemical contaminants such as heavy metals are relevant hazards for insect-based foods, and empirical studies in farmed insects indicate that diet can influence heavy-metal levels in larvae (Gori et al., 2025; Meyer et al., 2021). Heavy-metal accumulation has also been discussed in relation to insect substrates and feeding conditions, supporting the need for monitoring when alternative ingredients are introduced (EFSA Scientific Committee, 2015; Schrögel and Wätjen, 2019).

Sensory quality is another practical barrier to consumer acceptance of insect-based foods. Studies on beetle larvae have shown that processing interventions can reduce characteristic off-odors associated with specific volatile compounds, underscoring the value of volatile profiling for evaluating strategies designed to improve acceptability (Kim et al., 2021). Because odor-related attributes are often driven by complex mixtures of compounds, multivariate approaches can help interpret treatment-dependent shifts in overall volatile patterns (Baldovini and Chaintreau, 2020; Lubes and Goodarzi, 2017).

In this context, we evaluated whether formulated depuration diets provided during the gut-evacuation step could serve as a controllable pre-processing strategy to modulate the nutritional quality, chemical safety, and volatile profiles of P. brevitarsis larvae. Three formulated depuration diets composed of wheat bran, eggshell powder, and brewer’s spent grain at graded proportions (FD1-FD3) were compared with a glutinous rice flour-based control reflecting common practice. We assessed treatment-dependent differences in nutritional composition (proximate composition and fatty acid/mineral profiles), measured selected heavy metals relevant to food safety, and summarized multivariate variation in volatile profiles using principal component analysis (PCA).

2. Materials and methods

2.1. Developmental characteristics of P. brevitarsis larvae
2.1.1. Preparation of the formulated diet

Eggshells used to prepare the formulated diets were immersed in water and washed to remove foreign materials. After washing, excess surface water was removed at room temperature until no visible moisture remained. The eggshells were then dried in a hot-air dryer at 105°C for 2 h to remove residual moisture and sterilize the material. The dried eggshells were ground into a powder and used as an ingredient. Wheat bran and brewer’s spent grain were also dried to reduce moisture. All ingredients were weighed according to the composition shown in Table 1 and mixed thoroughly to obtain a uniform formulated diet for the fasting process.

Table 1. Composition (%) of the control diet and formulated diets provided to Protaetia brevitarsis larvae during the gut evacuation (depuration) step
Ingredient (%) Control Formulated Diet
GRP1),2) FD13) FD2 FD3
Glutinous rice powder 1004) - - -
Wheat bran -5) 47 73 90
Eggshell powder - 3 2 1
Brewer’s spent grain - 50 25 9

GRP, larvae provided with the control diet (glutinous rice powder) during depuration; FD1-FD3, larvae provided with formulated diets 1-3 during depuration.

The control diet consisted of glutinous rice powder only (100%, w/w).

FD1-FD3 indicate formulated diets prepared by blending wheat bran, eggshell powder, and brewer’s spent grain at the listed proportions.

Values are expressed as percentage by weight in the diet mixture (%, w/w) on a dry-matter basis.

“−” indicates that the ingredient was not included. Diet codes (GRP, FD1-3) correspond to the larval treatment groups receiving each diet during depuration.

Download Excel Table
2.1.2. Experimental insects and rearing

P. brevitarsis larvae were reared in the insect-rearing facility of the Industrial Insect Division, National Institute of Agricultural Sciences, Rural Development Administration (RDA, Korea). Rearing conditions were maintained at 25±1°C, 50-60% relative humidity, and a 16L:8D photoperiod. Fermented oak sawdust (Healing Bugs, Gimje, Korea) was used as the rearing substrate/feed. Only late third-instar larvae meeting the harvest body-mass criterion (approximately 2.5 g per individual) were selected (Park and Kim, 2025).

2.1.3. Fasting treatment and sample preparation

For gut-content removal as a pretreatment, larvae were provided with the control diet (glutinous rice powder, GRP) or one of three formulated diets for 48 h. Each diet was spread to a thickness of 1-2 cm in plastic rearing containers (54×40.5×18.5 cm; length×width×height). Two hundred larvae were placed in each container, and two independent replicate containers were prepared per treatment group under the same rearing conditions. The rearing container was defined as the biological experimental unit (n=2 containers per treatment). After feeding, the diet was removed, and larvae were fasted for 2 days (48 h) with adequate ventilation under a dry environment.

After fasting, larvae were washed three times with water to remove surface contaminants and gently blotted on filter paper to remove excess surface water. Washed larvae were rapidly frozen at -70°C for at least 24 h using an ultra-low temperature freezer (NIHON Freezer, Tokyo, Japan). Frozen samples were then freeze-dried for 65±5 h using a freeze dryer (Illshin Biobase, Dongducheon, Korea). The dried larvae were ground using a 100-mesh grinder (Korea Industry, Siheung, Korea), and the resulting powders were stored at -80°C until analysis.

2.2. Nutritional analysis of P. brevitarsis
2.2.1. Proximate composition

The proximate composition of P. brevitarsis larvae was determined according to the Korean Food Code of the Ministry of Food and Drug Safety (MFDS, Korea) and Association of Official Analytical Chemists (AOAC) methods. Moisture content was measured by atmospheric-pressure oven drying at 105°C to constant weight. Crude ash content was determined by dry ashing in a muffle furnace at 550°C. Crude protein content was determined using the micro-Kjeldahl method after digestion with a catalyst mixture (protein digestion accelerator) and concentrated sulfuric acid. Total nitrogen was quantified using an automated Kjeldahl analyzer (Kjeltec 8400, GRPS Analytical, Hillerød, Denmark), and crude protein was calculated using a nitrogen-to-protein conversion factor of 6.25. Crude fat content was determined using a Soxhlet apparatus with diethyl ether as the solvent, and the extracted residue was quantified gravimetrically after solvent removal. Crude fiber was determined using a filter bag filtration method. Carbohydrate content was calculated as 100 − (moisture + crude protein + crude fat + crude ash). Results were expressed as weight percentages (% w/w), as specified.

2.2.2. Fatty acid composition

Fatty acid composition was determined according to the Korean Food Code using a Folch-type solvent extraction method, conversion to fatty acid methyl esters (FAMEs), and subsequent GC analysis. Briefly, 50 g of homogenized larval sample was extracted with 250 mL chloroform:methanol (2:1, v/v) using a homogenizer operated at 3,000 rpm. The extract was dehydrated over anhydrous sodium sulfate and concentrated in a water bath. The lipid fraction was saponified with 1 mL of 0.5 N NaOH at 100°C for 20 min and cooled at room temperature for 30 min, then methylated by adding 2 mL boron trifluoride reagent, followed by an additional 20-min heating step and a 30-min cooling period. Phase separation was induced by adding 1 mL n-heptane and 8 mL saturated NaCl solution, and the upper organic layer containing FAMEs was collected. FAMEs were analyzed using GC (US/HP 6890, Agilent Technologies, Seoul, Korea) equipped with an SP-2560 silica capillary column (100 m×0.25 mm i.d., 0.20 μm). The inlet temperature was set to 225°C, and the sample injection volume was 1 μL.

2.2.3. Mineral composition

Mineral contents were analyzed according to the Korean Food Code. Dried larval powder (50 mg) was subjected to dry ashing via sequential drying and carbonization, followed by complete ashing at 450-550°C in a muffle furnace. The resulting ash was dissolved in hydrochloric acid and filtered through a glass fiber filter. Elemental concentrations were quantified using inductively coupled plasma-optical emission spectrometry (ICP-OES; Horiba, Kyoto, Japan). Molybdenum (Mo) and selenium (Se) were analyzed after dry-digestion pretreatment. The chlorine (Cl) content was determined by ashing 1 g of a salt-containing sample; the ash was dissolved in water, filtered, and titrated with silver nitrate using potassium chromate as an indicator.

2.2.4. Heavy metal composition

Heavy metal contents were determined according to the Korean Food Code. Samples were dry-ashed at 600°C for >2 h in an electric furnace. After cooling to room temperature, the ash was dissolved in 1:1 (v/v) hydrochloric acid and allowed to stand for 18 h. The solution was filtered through Whatman No. 6 filter paper (Whatman International Co., Maidstone, UK). Lead (Pb), cadmium (Cd), and arsenic (As) were quantified using ICP-OES (Horiba, Kyoto, Japan) at 220.35, 226.50, and 193.70 nm, respectively.

2.3. Electronic nose analysis

Electronic nose (E-nose) analysis was performed following the methods described by Kang et al. (2020) and Oh et al. (2025) with some modifications. For multicomponent odor profiling, 2.5 g of finely ground sample was mixed with 5 mL of triple-distilled water and stirred. The mixture was transferred to a 20-mL vial (Ls-Phs-Psck GmbH, Langerwehe, Germany), and headspace volatile compounds were generated by incubating the vial at 40°C for 20 min with agitation at 500 rpm. Volatile profiles were analyzed using an E-nose system (Alpha MOS, Toulouse, France) equipped with two columns (MTX-5 and MTX-1701) and a flame ionization detector (FID). The column oven temperature program was as follows: 50°C for 2 s, increased at 1°C/s to 80°C, then increased at 3°C/s to 250°C and held for 21 s. Injection was performed using a syringe-type injector at 25°C with a column head pressure of 1.0 psi. During analysis, injector pressure and detector pressure were set to 1.0 and 39.0 psi, respectively. Volatile compounds were tentatively identified based on retention time and the instrument library, and peak-area data were used for subsequent analyses. PCA was performed using Alpha MOS software to visualize treatment-dependent differences in overall volatile profiles, and discrimination power (DP) values were used to identify marker compounds contributing to sample differentiation.

2.4. Statistical analysis

Results are presented as mean±standard deviation. Statistical comparisons among the four treatment groups (GRP, FD1, FD2, and FD3) were performed using one-way analysis of variance (ANOVA). When significant treatment effects were detected, Tukey’s honestly significant difference (HSD) test was applied for multiple comparisons. Statistical analyses were performed using PASW Statistics 22.0 for Windows (IBM Corp., Chicago, IL, USA), and differences were considered significant at p>0.05. Electronic-nose data were further analyzed using Alpha MOS software (Alpha MOS, Toulouse, France). PCA was applied to the peak-area data to compare overall volatile patterns among treatments, and score and loading plots were used for visualization. Compounds with DP values >0.90 were regarded as major discriminant volatiles.

3. Results and discussion

3.1. Nutritional evaluation of P. brevitarsis larvae following a gut-emptying pretreatment
3.1.1. Comparative analysis of proximate composition

Proximate composition was determined in whole-larvae powders collected after a gut-emptying pretreatment, comparing the control group (GRP) with larvae subjected to formulated-diet treatments (FD1-FD3) (Table 2). Crude protein content was 53.2±2.9% in GRP and ranged from 55.6±0.7% to 57.4± 0.7% in the FD groups. Although FD1 and FD2 showed numerically higher crude protein than GRP, Tukey’s HSD test did not detect significant differences among treatment groups for crude protein or any other variables listed in Table 2 (all NS).

Table 2. Proximate composition and dietary fiber of Protaetia brevitarsis larvae in the control (GRP) and formulated diet (FD1-FD3) treatment groups (%)
Proximate composition (%) Control Formulated Diet
GRP1) FD1 FD2 FD3
Crude protein 53.2±2.92) 56.9±0.3 57.4±0.7 55.6±0.7
Crude fat 12.5±1.0 10.8±0.3 11.4±0.1 11.9±0.0
Crude ash 7.5±1.2 9.7±0.4 8.6±0.6 9.9±0.2
Dietary fiber 18.7±0.0 17.4±0.3 17.6±1.2 16.5±0.1
Carbohydrate3) 26.3±1.4 22.6±0.4 20.9±2.3 22.7±0.5NS4)

GRP, larvae provided with the control diet (glutinous rice powder) during depuration (control group); FD1-FD3, larvae provided with formulated diets 1-3 during depuration (formulated-diet groups).

All values are mean±SD (n=2, % w/w).

All values are expressed on a dry matter (DM) basis. Carbohydrate was calculated by difference [100 - (crude protein + crude fat + crude ash)].

NS, not significant differences were observed among treatment groups for any variable in this table by one-way ANOVA (p>0.05).

Download Excel Table

Crude fat was numerically lower in FD1-FD3 (10.8-11.9%) than in GRP (12.5±1.0%), whereas crude ash was numerically higher in FD1-FD3 (8.6-9.9%) compared with GRP (7.5±1.2%). Dietary fiber ranged from 16.5±0.1% to 18.7±0.0%, and carbohydrate ranged from 20.9±2.3% to 26.3±1.4% across the treatment groups. These numerical differences were not statistically significant under the present experimental design.

Collectively, the formulated-diet pretreatment should be interpreted as producing numerical tendencies in proximate composition, dietary fiber, and carbohydrate rather than statistically supported compositional shifts. Such diet-dependent variability in insect proximate composition is well documented; species, developmental stage, feeding substrate, and processing conditions—particularly the feeding substrate—can influence crude protein and crude fat levels (Janssen et al., 2017; Oonincx et al., 2015; Rumpold and Schl, 2013). However, interpretation of short-term pretreatment effects requires caution because the current analysis was performed on whole larvae, including the gastrointestinal tract. Consequently, observed numerical shifts may reflect not only changes in larval tissues but also differences in residual gut contents (e.g., feed residues and frass) remaining at sampling.

To disentangle these effects, future work should (1) quantify the residual mass of gut contents before and after pretreatment and (2) evaluate gut-emptying efficiency metrics (e.g., frass output and gut-content weight) alongside proximate composition, thereby clarifying the pretreatment’s role as a controllable process variable affecting analytical outcomes (Kipkoech, 2023; Yoon et al., 2020).

3.1.2. Comparative analysis of fatty acid composition

Thirteen individual fatty acids were quantified (Fig. 1), and oleic acid was consistently the predominant fatty acid across all groups. Oleic acid was highest in FD3 larvae (11.275 g/100 g), followed by FD2 (7.500 g/100 g), GRP (7.050 g/100 g), and FD1 (6.227 g/100 g). Tukey’s HSD grouping indicated that FD3 was significantly higher than FD2, GRP, and FD1; FD2 was significantly higher than FD1; and GRP did not differ significantly from FD1 or FD2.

kjfp-33-3-442-g1
Fig. 1. Fatty acid concentrations of Protaetia brevitarsis larvae in the control (GRP) and formulated-diet (FD1-FD3) treatment groups after gut-evacuation pre-treatment. All values are mean± SD of two biologically independent replicate containers (n=2 per treatment). Fatty acids are grouped as saturated and unsaturated, as indicated in the figure. Different lowercase letters above bars indicate significant differences among treatment groups for the same fatty acid at p>0.05 by one-way ANOVA followed by Tukey’s HSD test; NS indicates no significant difference among treatment groups. GRP, larvae provided with the control diet (glutinous rice powder); FD1-FD3, larvae provided with formulated diets 1-3 during the depuration step.
Download Original Figure

Among the remaining major fatty acids, palmitic acid (2.035-2.524 g/100 g), palmitoleic acid (0.960-2.246 g/100 g), and linoleic acid (0.726-1.486 g/100 g) were consistently abundant. Significant Tukey groupings were also observed for palmitoleic acid, linoleic acid, α-linolenic acid, eicosenoic acid, EPA, and lauric acid, whereas the remaining fatty acids were annotated as NS. Total quantified fatty acids were 11.959 g/100 g in GRP, 10.824 g/100 g in FD1, 13.237 g/100 g in FD2, and 16.819 g/100 g in FD3 larvae, with FD3 showing the highest total fatty acid value. When fatty acids were grouped by saturation class, the sum of saturated fatty acids (SFAs) ranged from 2.481 to 3.080 g/100 g, whereas the sum of unsaturated fatty acids (UFAs) ranged from 8.334 to 14.338 g/100 g, indicating that UFAs predominated across treatments (values calculated from Fig. 1).

Insects can exhibit diet-responsive fatty acid profiles, as dietary fatty acids and carbon sources can influence the relative abundance of specific fatty acids, including unsaturated fatty acids (Oonincx et al., 2015; Yoon et al., 2020). In P. brevitarsis larvae, variability in the distribution of major fatty acids—particularly oleic, palmitic, palmitoleic, and linoleic acids—has likewise been reported across rearing conditions and feeding substrates (Eom et al., 2024; Ham et al., 2021). Collectively, these results support the interpretation that formulated-diet pretreatment could serve as a controllable process variable to modulate selected fatty acids, as reflected by the Tukey letter groupings for oleic acid and other diet-responsive fatty acids.

Nevertheless, it cannot be determined conclusively whether the values in Fig. 1 primarily reflect physiological remodeling of lipid synthesis and/or storage during the pretreatment period, or whether they partly represent carryover effects from residual gut contents derived from the pretreatment diets. Although a 48-h starvation step was applied, the completeness of gut-content removal was not independently verified using quantitative indicators. Therefore, follow-up studies should integrate (1) the fatty acid composition of each diet, (2) quantitative measures of gut-emptying efficiency (e.g., frass output and residual gut-content mass), and (3) larval fatty acid composition to distinguish absorption and metabolic changes from residual-diet effects (Oonincx et al., 2015; Stanley‐Samuelson et al., 1988).

3.1.3. Comparative analysis of mineral composition

In the macromineral analysis (mg/100 g), Calcium (Ca) was significantly higher in all formulated-diet treatment groups than in GRP according to one-way ANOVA followed by Tukey’s HSD test. Ca increased from 243.7±6.2 mg/100 g in GRP to 694.1±45.3 mg/100 g in FD1, 644.1±47.1 mg/100 g in FD2, and 744.2±66.2 mg/100 g in FD3. Phosphorus (P) ranged from 993.5 to 1,191.5 mg/100 g in FD1-FD3 larvae and did not differ significantly from GRP larvae (1,015.0± 44.9 mg/100 g). Potassium (K) was significantly higher in FD1 larvae (3,407.6±81.4 mg/100 g) than in GRP larvae (3,089.5±79.9 mg/100 g) (p>0.05), whereas magnesium (Mg) did not differ significantly among treatments (388.0-428.5 mg/100 g). Regarding microminerals, zinc (Zn) was significantly higher in FD1 (11.5±0.4 mg/100 g) and FD2 (11.0±0.9 mg/100 g) than in GRP (8.8±0.0 mg/100 g), whereas FD3 was intermediate and did not differ significantly from either GRP or FD1-FD2 (Fig. 2B). Iron (Fe) was significantly higher in all FD groups than in GRP; among the FD groups, Fe was highest in FD1, lower in FD2, and intermediate in FD3.

kjfp-33-3-442-g2
Fig. 2. Mineral concentrations of Protaetia brevitarsis larvae in the control (GRP) and formulated-diet (FD1-FD3) treatment groups after gut-evacuation pre-treatment. (A), Macrominerals (Ca, P, K, and Mg); (B), Microminerals (Zn and Fe). Mineral contents were determined after dry ashing using ICP-OES. All values are mean±SD of two biologically independent replicate containers (n=2 per treatment). Different lowercase letters above bars indicate significant differences among treatment groups for the same mineral at p>0.05 by one-way ANOVA followed by Tukey’s HSD test; NS indicates no significant difference among treatment groups. GRP, larvae provided with the control diet (glutinous rice powder) during the depuration step; FD1-FD3, larvae provided with formulated diets 1-3 during the depuration step.
Download Original Figure

Overall, the FD pretreatments produced statistically supported increases in Ca and Fe (approximately 2.6-3.1-fold and 2.2-2.8-fold relative to GRP, respectively), whereas no significant treatment effects were detected for P, K, or Mg. Collectively, these results suggest that the short-term formulated diets used during the pretreatment step can alter—or potentially confound—the apparent mineral profile of P. brevitarsis larvae, particularly for Ca and Fe. Previous studies have reported that larval composition varies with feeding sources and rearing conditions (Ham et al., 2021; Yoon et al., 2020), and the present data further suggest that detectable shifts may occur even during brief gut-clearing treatments (Ham et al., 2021; Park et al., 2025a).

Notably, Ca was highest in FD3 despite FD3 having the lowest proportion of eggshell powder, indicating that the observed Ca increase cannot be explained solely by the eggshell inclusion rate. Several non-mutually exclusive factors may contribute, including baseline mineral contributions from other ingredients (e.g., wheat bran and brewer’s spent grain), variation in the extent of gut-content removal, and/or short-term redistribution following intake. To disentangle these effects, future studies should quantify mineral concentrations in individual ingredients (especially Ca, Fe, Zn, and K) and incorporate sampling strategies that more effectively separate larval tissue from residual gastrointestinal contents after pretreatment (Jankauskienė et al., 2024; Noyens et al., 2023).

3.1.4. Comparative analysis of hazardous substances

Concentrations of Pb, Cd, and total As (mg/kg, dry weight) were low across all treatments (Table 3). Pb showed a significant Tukey grouping, with FD1-FD3 assigned to a higher group than GRP; however, the absolute difference was small (approximately 0.03 mg/kg). Cd and total As were annotated as NS, indicating no significant differences among treatment groups. Overall, the maximum observed concentrations were 0.069 mg/kg for Pb, 0.020 mg/kg for Cd, and 0.031 mg/kg for total As.

Table 3. Heavy metal contents of Protaetia brevitarsis larvae in the control (GRP) and formulated-diet (FD1-FD3) treatment groups after gut-evacuation pre-treatment
Hazardous substance Content Control Formulated diets
GRP1) FD1 FD2 FD3
Heavy metals (mg/kg) Lead (Pb) 0.040±0.0022)b3) 0.068±0.002a 0.067±0.010a 0.069±0.006a
Cadmium (Cd) 0.011±0.001NS4) 0.020±0.002NS 0.020±0.004NS 0.019±0.004NS
Arsenic (As) 0.018±0.005NS 0.031±0.002NS 0.031±0.003NS 0.028±0.004NS

Abbreviations: GRP, larvae provided with the control diet (glutinous rice powder) during depuration (control group); FD1-FD3, larvae provided with formulated diets 1-3 during depuration (formulated-diet groups).

All values are mean±SD (n=2, % w/w).

Different superscript lowercase letters, within the same row indicate significant differences among treatment groups at p<0.05 by one-way ANOVA followed by Tukey’s HSD test.

NS indicates no significant differences among treatment groups within the same row (p>0.05).

Download Excel Table

Because chemical hazards in insect-based foods—particularly toxic elements such as heavy metals—can be strongly influenced by the contamination profile of the rearing feed/substrate, monitoring is widely recognized as a key control point when feed ingredients or substrate sources change (EFSA Scientific Committee, 2015; van der Fels-Klerx et al., 2016). From a regulatory perspective, the Korean Food Code was revised under MFDS Notification No. 2021-54 to harmonize heavy-metal criteria across edible insect species and to regulate arsenic based on inorganic arsenic rather than total arsenic. Under these provisions for dried edible insects (dry-weight basis), the maximum levels are ≤0.1 mg/kg for Pb, Cd, and inorganic As, with a higher Pb limit (≤0.3 mg/kg) applied to P. brevitarsis larvae and Allomyrina dichotoma larvae (EFSA Scientific Committee, 2015; van der Fels-Klerx et al., 2018).

In Korea, the Food Code sets maximum levels for dried edible insects of 0.1 mg/kg for Pb (with a higher Pb limit of 0.3 mg/kg applied to P. brevitarsis larvae and Allomyrina dichotoma larvae), 0.1 mg/kg for Cd, and 0.1 mg/kg for As (MFDS, 2021). In this study, the maximum observed concentrations of Pb and Cd were 0.07 and 0.02 mg/kg, respectively, both well below these limits. Arsenic concentrations were also low, with total As not exceeding 0.03 mg/kg in any treatment and remaining below the 0.1 mg/kg regulatory threshold.

These results indicate that, despite the small statistically significant increase in Pb, the short-term formulated-diet gut-clearing pretreatment (FD1-FD3) did not introduce appreciable Pb, Cd, or As contamination under the present experimental conditions. Nevertheless, chemical hazards in edible insects depend strongly on the contamination profile of the rearing substrate and can exhibit element- and species-specific transfer and bioaccumulation (Malematja et al., 2023; Meyer et al., 2021; Truzzi et al., 2019; van der Fels-Klerx et al., 2016). Accordingly, for industrial application, changes in feed formulation or raw-material sourcing should be accompanied by routine screening of feed ingredients and lot-based monitoring of Pb, Cd, and As in both substrates and larvae.

3.2. E-nose analysis of volatile profiles

The selected volatile compounds detected by the E-nose and their peak areas are summarized in Table 4, and the overall treatment-dependent variation is presented in the PCA score/loading plot in Fig. 3. In the loading plot, none of the selected discriminant volatiles were located in the first quadrant. Instead, 2,3-pentanedione and butanal clustered in the second quadrant; heptanal and butane-2,3-dione in the third quadrant; and anethole, 2-heptanone, and propyl propanoate in the fourth quadrant. All seven compounds showed DP values >0.90, indicating strong contributions to sample differentiation. PC1 and PC2 explained 99.992% and 0.007428% of the total variance, respectively, together accounting for 99.999428% of the overall variation. In the score plot, FD1 and FD3 were positioned close to each other, indicating similar overall volatile profiles. Conversely, the control group (GRP) was clearly separated from the formulated-diet treatments, and FD2 was also distinctly separated from the other groups and positioned closer to the fourth-quadrant marker compounds. Because odor-related attributes arise from complex mixtures of volatile compounds, this PCA pattern suggests that the depuration diets altered the overall odor-related characteristics of P. brevitarsis larvae rather than changing only a single volatile compound (Cai et al., 2024).

Table 4. Selected volatile compounds, odor descriptors, and peak areas detected by electronic nose in Protaetia brevitarsis larvae from the control (GRP) and formulated-diet (FD1-FD3) treatment groups after gut-evacuation pretreatment
Column Compound1) Rt (min)2) Odor descriptor3) Peak area
Control (GRP4)) FD1 FD2 FD3
MTX-5 Butane-2,3-dione 19.1 Butter, caramelized, chlorine, creamy, fruity 1,298.1±76.75)a6) 732.0±26.7b 872.0±73.8b 293.4±199.4c
Butanal 22.3 Chocolate, cocoa, green, malty, musty 239,932.6±47,551.3a 14,072.7±648.5b 19,854.1±1428.8b 17,802.7±435.2b
2,3-Pentanedione 32.9 Almond, apple, butter, burnt, butterscotch 4,445.7±1,032.3a ND7) ND ND
Propyl propanoate 44.8 Apple, chemical, fruity, fruity (sweet), pineapple 1,175.6±116.5c 1,770.1±219.2b 3,310.5±194.6a 1,542.7±207.1bc
2-Heptanone 52.0 Banana, cheese, coconut, fruity, gaseous 209.3±82.3c 447.2±30.0b 1,286.9±96.9a 414.6±57.6b
Heptanal 52.9 Aldehydic, citrus, fatty, fish (dry), fresh 2,391.3±443.7a 570.4±81.4c 1,164.3±97.0b 649.7±67.1bc
MTX-1701 Anethole 83.1 Anise, herbaceous, licorice, oil, spicy ND ND 91.6±11.9a ND

Volatile compounds were tentatively identified by electronic-nose analysis based on retention time and the instrument library.

Rt, retention time.

Odor descriptors are library-based terms associated with the tentatively identified volatile compounds and do not represent sensory-panel data.

GRP, larvae provided with the control diet (glutinous rice powder) during the depuration step (control group); FD1-FD3, larvae provided with formulated diets 1-3 during the depuration step (formulated-diet groups).

All values are mean±SD (n=3).

Different superscript letters within the same row indicate significant differences among treatment groups at p<0.05 by one-way ANOVA followed by Tukey’s HSD test.

ND, not detected.

Download Excel Table

Table 4 further supports the PCA-based separation. Butanal exhibited the largest peak area in GRP and was significantly lower in all formulated-diet groups, without a clear formulation-dependent trend. Butane-2,3-dione and heptanal were also highest in GRP and decreased in the formulated-diet groups, whereas 2,3-pentanedione was detected only in GRP. Conversely, propyl propanoate and 2-heptanone showed their highest peak areas in FD2, and anethole was detected exclusively in FD2. Taken together with the library-assigned odor descriptors in Table 4, these results indicate that GRP was more closely associated with butanal, 2,3-pentanedione, heptanal, and butane-2,3-dione, whereas FD2 was more closely associated with propyl propanoate, 2-heptanone, and anethole. This shift suggests that the conventional glutinous-rice-based control and the formulated depuration diets generated different volatile signatures, with FD2 showing the most distinctive treatment-specific pattern among the formulated-diet groups. The close placement of FD1 and FD3 in the score plot, despite differences in formulation, further suggests that their overall volatile patterns were more similar to each other than to GRP or FD2. Similar treatment-dependent changes in volatile composition have been reported for P. brevitarsis under different fermentation or rearing conditions, supporting the interpretation that the depuration diets used in the present study contributed to the observed differences in volatile profiles (Cha et al., 2024a; Cha et al., 2024b).

kjfp-33-3-442-g3
Fig. 3. PCA score/loading plot of electronic-nose-derived volatile profiles of Protaetia brevitarsis larvae from the control (GRP) and formulated-diet (FD1-FD3) treatment groups after gut-evacuation pretreatment. PCA was performed using peak-area data obtained from electronic-nose analysis. Sample positions represent treatment-group scores, and labeled compounds indicate the loadings of the selected volatile compounds contributing to group discrimination. Samples located closer together exhibit more similar overall volatile profiles. PC1 and PC2 explained 99.992% and 0.007428% of the total variance, respectively. Only compounds with DP values greater than 0.90 are shown. GRP, larvae provided with the control diet (glutinous rice powder) during the depuration step (control group); FD1-FD3, larvae provided with formulated diets 1-3 during the depuration step (formulated-diet groups).
Download Original Figure

These findings are important because sensory quality remains a practical barrier to broader acceptance of insect-based foods, and odor profile is closely associated with consumer acceptance of P. brevitarsis larvae products (Cha et al., 2024a). However, the current E-nose approach is more suitable for rapid VOC fingerprinting and comparative pattern recognition than for definitive identification of key odorants. Accordingly, the odor descriptors in Table 4 should be interpreted as library-based tentative assignments rather than direct sensory-panel results (Cai et al., 2024).

In the present study, the formulated depuration diets did not simply reduce the abundance of selected GRP-associated compounds; rather, they shifted the dominant volatile markers toward a different profile. However, the odor descriptors listed in Table 4 are library-based terms for tentatively identified compounds and should not be interpreted as direct sensory-panel results. In addition, because the larvae were processed as whole-body materials and the analyzed samples were whole-larva powders obtained after 48 h of diet provision followed by 48 h of fasting, the observed differences may reflect both changes in larval tissue composition and differences in residual gut contents. Nevertheless, these data indicate that diet formulation during depuration can influence the volatile profile of P. brevitarsis larvae and may therefore serve as a practical pre-processing factor for modulating odor-related quality.

4. Conclusions

Formulated depuration diets altered the nutritional composition and volatile characteristics of P. brevitarsis larvae while maintaining Pb, Cd, and As contamination at negligible levels under the present conditions. Electronic-nose analysis further revealed clear diet-dependent separation of volatile profiles, with FD2 showing the most distinct pattern among the formulated-diet groups. These findings suggest that depuration diet formulation can serve as a practical pre-processing factor for modulating the quality of P. brevitarsis larvae, rather than functioning only as a gut-emptying step. However, because this study used a limited number of independent replicate containers, further validation with larger biological replication is required before broad industrial generalization.

Acknowledgements

None.

Conflict of interests

The authors declare no potential conflicts of interest.

Author contributions

Conceptualization: Park MJ, Kim SY. Methodology: Park MJ, Kim JY, Park JY. Formal analysis: Park MJ, Kim JY, Park JY. Validation: Kim SY. Data curation: Park MJ. Writing - original draft: Park MJ. Writing - review & editing: Kim JY, Park JY, Kim SY.

Ethics approval

This article does not require IRB/IACUC approval because there are no human and animal participants.

Funding

This work was supported by the Cooperative Research Program for Agriculture Science and Technology Development, Rural Development Administration, Republic of Korea (grant no. PJ01749501). Min Ji Park was additionally supported by the 2025 RDA Fellowship Program, National Institute of Agricultural Sciences, Rural Development Administration, Republic of Korea.

References

1.

Acosta-Estrada BA, Reyes A, Rosell CM, Rodrigo D, Ibarra-Herrera CC. Benefits and challenges in the incorporation of insects in food products. Front Nutr. 8:687712 2021;

2.

Ahn MY, Hwang JS, Yoon HJ, Park K, Kim S, Kim EM. Fasting conditions and dietary phenamena of edible cricket (Gryllus bimaculatus). J Seric Entomol Sci. 53:78-81 2015;

3.

Baigts-Allende DK, Stathopoulos C. Overcoming obstacles in insect utilization. Eur Food Res Technol. 249:849-860 2023;

4.

Baldovini N, Chaintreau A. Identification of key odorants in complex mixtures occurring in nature. Nat Prod Rep. 37:1589-1626 2020;

5.

Cacchiarelli C, Fratini F, Puccini M, Vitolo S, Paci G, Mancini S. Effects of different blanching treatments on colour and microbiological profile of Tenebrio molitor and Zophobas morio larvae. LWT. 157:113112 2022;

6.

Cai W, Zhou W, Liu J, Wang J, Kuang D, Wang J, Long Q, Huang D. An exploratory study on the rapid detection of volatile organic compounds in Gardenia fruit using the heracles neo ultra-fast gas phase electronic nose. Metabolites. 14:445 2024;

7.

Cappelli A, Cini E, Lorini C, Oliva N, Bonaccorsi G. Insects as food: A review on risks assessments of Tenebrionidae and Gryllidae in relation to a first machines and plants development. Food Control. 108:106877 2020;

8.

Cha JY, Han J, Heo JA, Yu HH, Kim YJ, Jang HW, Kim MR, Choi YS. Variation of volatile compounds and sensory profile for Protaetia brevitarsis larvae fermented with lactic acid bacteria and yeast. Food Chem. 452:139480 2024a;

9.

Cha JY, Kim TK, Kim YJ, Lee JH, Kang MC, Jang HW, Choi YS. Comparison of plastic and glass cages on volatile compounds in protein extracted from Protaetia brevitarsis larvae. Future Foods. 10:100429 2024b;

10.

Chung MY, Kwon EY, Hwang JS, Goo TW, Yun EY. Pre-treatment conditions on the powder of Tenebrio molitor for using as a novel food ingredient. J Seric Entomol Sci. 51:9-14 2013;

11.

EFSA Scientific Committee. Risk profile related to production and consumption of insects as food and feed. EFSA J. 13:4257 2015;

12.

Egonyu JP, Labu S, Nyangena DN, Khamis F, Cheseto X, Mbi C, Sevgan T. Pre‑harvest starvation effects on microbial load, weight loss and proximate composition of edible field cricket (Gryllus bimaculatus) and desert locust (Schistocerca gregaria). Int J Trop Insect Sci. 45:2539-2547 2025;

13.

Eom H, Kim H, Hahn D, Kwon O. Development of feed material and its effect on the nutritional composition of Protaetia brevitarsis larvae. Ento Res. 54:e12711 2024;

14.

Fraqueza MJR, Patarata LAdaSC. Constraints of HACCP application on edible insect for food and feed. In: Future Foods. In: Mikkola H, editor.IntechOpen. Rijeka, Croatia: p p. 89-113 2017;

15.

Fuso A, Leni G, Prandi B, Lolli V, Caligiani A. Novel foods/feeds and novel frauds: The case of edible insects. Trends Food Sci Technol. 147:104457 2024;

16.

Gori A, Armani A, Pedonese F, Benini O, Mancini S, Nuvoloni R. Heavy metals (Pb, Cd, Ni) in insect-based products for human consumption sold by e-commerce in the EU market: Occurrence and potential health risk associated with dietary exposure. Food Control. 167:110781 2025;

17.

Ham YK, Kim SW, Song DH, Kim HW, Kim IS. Nutritional composition of white-spotted flower chafer (Protaetia Brevitarsis) larvae produced from commercial insect farms in Korea. Food Sci Anim Resour. 41:416-427 2021;

18.

Han G, Choi M, Jo H, Nam M, Lee J, Choi Y, Lee G, Kim B. Identification condition of according to pre-treatment conditions of Protaetia brevitarsis defecation-induction feeding using processing. Korean J Appl Entomol. 63:323-333 2024; .

19.

Jankauskienė A, Falkauskas R, Aleknavičius D, Kiseliovienė S, Antanaitis Š, Šumskienė M, Juknienė I, Kabašinskienė A. The influence of different sustainable substrates on the nutritional value of Tenebrio molitor larvae. Foods. 13:365 2024; .

20.

Janssen RH, Vincken JP, van den Broek LAM, Fogliano V, Lakemond CMM. Nitrogen-to-protein conversion factors for three edible insects: Tenebrio molitor, Alphitobius diaperinus, and Hermetia illucens. J Agric Food Chem. 65:2275-2278 2017;

21.

Kang JE, Kim HE, Lim BR, Choi HS, Jeong ST. Analysis of physicochemical characteristics in yakju with variety nuruk. Korean J Community Living Sci. 31:15-24 2020;

22.

Kim HY, Kim D, Chun J. Quality evaluation of hot air-dried edible insect ‘Protaetia brevitarsis larvae’ including sacrificing method, storage temperature, and defatting condition. Food Eng Prog. 26:69-77 2022;

23.

Kim J, Lee HE, Kim Y, Yang J, Lee SJ, Jung YH. Development of a post-processing method to reduce the unique off-flavor of Allomyrina dichotoma: Yeast fermentation. LWT. 150:111940 2021;

24.

Kim S, Park MJ, Lee JH, Kim SY. Comparative analysis of nutritional components and harmful substances of Protaetia brevitarsis larvae processed by different drying methods. Korean J Appl Entomol. 64:395-403 2025; .

25.

Kim SY, Park MJ, Lee H, Koo B. Analysis of nutritional components and harmful substances of protaetia brevitarsis larvae fed fermented oyster mushroom cultivation by-products. Korean J Appl Entomol. 63:373-382 2024; .

26.

Kipkoech C. Beyond proteins: Edible insects as a source of dietary fiber. Polysaccharides. 4:116-128 2023;

27.

Kooh P, Jury V, Laurent S, Audiat-Perrin F, Sanaa M, Tesson V, Federighi M, Boué G. Control of biological hazards in insect processing: Application of HACCP method for yellow mealworm (Tenebrio molitor) powders. Foods. 9:1528 2020;

28.

Kooh P, Ververis E, Tesson V, Boué G, Federighi M. Entomophagy and public health: A review of microbiological hazards. Health. 11:1272-1290 2019;

29.

Kouřimská L, Adámková A. Nutritional and sensory quality of edible insects. NFS J. 4:22-26 2016;

30.

Kwon EY, Yoo J, Yoon YI, Hwang JS, Goo TW, Kim MA, Choi YC, Yun EY. Pre-treatment of the white-spotted flower chafer (protaetia brevitarsis) as an ingredient for novel foods. J Korean Soc Food Sci Nutr. 42:397-402 2013;

31.

Liceaga AM, Aguilar-Toalá JE, Vallejo-Cordoba B, González-Córdova AF, Hernández-Mendoza A. Insects as an alternative protein source. Annu Rev Food Sci Technol. 13:19-34 2022;

32.

Lubes G, Goodarzi M. Analysis of volatile compounds by advanced analytical techniques and multivariate chemometrics. Chem Rev. 117:6399-6422 2017;

33.

Malematja E, Manyelo TG, Sebola NA, Kolobe SD, Mabelebele M. The accumulation of heavy metals in feeder insects and their impact on animal production. Sci Total Environ. 885:163716 2023;

34.

Mancini S, Fratini F, Tuccinardi T, Degl’Innocenti C, Paci G. Tenebrio molitor reared on different substrates: Is it gluten free?. Food Control. 110:107014 2020;

35.

Megido RC, Desmedt S, Blecker C, Béra F, Haubruge É, Alabi T, Francis F. Microbiological load of edible insects found in Belgium. Insects. 8:12 2017;

36.

Meyer AM, Meijer N, Hoek-van den Hil EFH, van der fels-Klerx HJ. Chemical food safety hazards of insects reared for food and feed. J Insects Food Feed. 7:823-831 2021;

37.

Meyer-Rochow VB, Gahukar RT, Ghosh S, Jung C. Chemical composition, nutrient quality and acceptability of edible insects are affected by species, developmental stage, gender, diet, and processing method. Foods. 10:1036 2021;

38.

Noh CW, Jeon SH, Son D, Cho YS, Lee BJ. Changes of nutritive component with before processing feeding type for larva of Protaetia brevitarsis. J Korean Soc Int Agric. 27:675-681 2015;

39.

Noyens I, Schoeters F, Peer MV, Berrens S, Goossens S. The nutritional profile, mineral content and heavy metal uptake of yellow mealworm reared with supplementation of agricultural sidestreams. Sci Rep. 13:11604 2023;

40.

Oh BY, Shin SJ, Oh SW. Quality changes in wine with rose petals (Magic Red) during alcohol fermentation and aging. Food Sci Preserv. 32:553-564 2025;

41.

Oonincx DGAB, Finke MD. Nutritional value of insects and ways to manipulate their composition. J Insects Food Feed. 7:639-659 2021;

42.

Oonincx DGAB, Van Broekhoven S, Van Huis A, Van Loon JJA. Feed conversion, survival and development, and composition of four insect species on diets composed of food by-products. PLoS One. 10:e0144601 2015;

43.

Park KM, Kim SY, Lee YS. Behavior of inoculated Staphylococcus aureus during the rearing of Protaetia brevitarsis larvae on fermented oak sawdust as a food product. J Food Prot. 88:100588 2025a;

44.

Park MJ, Kim SY. Persimmon leaf as a functional additive for enhancing productivity and nutritional value in Protaetia brevitarsis larvae. Int J Ind Entomol Biomater. 49:43-54 2025; .

45.

Park MJ, Park K, Lee Y, Lee KY, Kim SY. By-products of king oyster mushroom cultivation in feed enhance growth and nutrition in Protaetia brevitarsis larvae. J Environ Sci Int. 34:645-651 2025b;

46.

Park MJ, Yoon HJ, Kim SY. Analyzing the nutritional value of Protaetia brevitarsis larvae feeding on coffee and oyster mushroom cultivation by-products. Int J Ind Entomol Biomater. 48:147-155 2024; .

47.

Rumpold BA, Schlüter OK. Nutritional composition and safety aspects of edible insects. Mol Nutr Food Res. 57:802-823 2013;

48.

Schrögel P, Wätjen W. Insects for food and feed-safety aspects related to mycotoxins and metals. Foods. 8:288 2019;

49.

Selaledi L, Mabelebele M. The influence of drying methods on the chemical composition and body color of yellow mealworm (Tenebrio molitor L.). Insects. 12:333 2021;

50.

Stanley‐Samuelson DW, Jurenka RA, Cripps C, Blomquist GJ, de Renobales M. Fatty acids in insects: Composition, metabolism, and biological significance. Arch Insect Biochem Physiol. 9:1-33 1988;

51.

Traynor A, Burns DT, Wu D, Karoonuthaisiri N, Petchkongkaew A, Elliott CT. An analysis of emerging food safety and fraud risks of novel insect proteins within complex supply chains. Npj Sci Food. 8:7 2024;

52.

Truzzi C, Illuminati S, Girolametti F, Antonucci M, Scarponi G, Ruschioni S, Riolo P, Annibaldi A. Influence of feeding substrates on the presence of toxic metals (Cd, Pb, Ni, As, Hg) in larvae of Tenebrio molitor: Risk assessment for human consumption. Int J Environ Res Public Health. 16:4815 2019;

53.

van der Fels-Klerx HJ, Camenzuli L, Belluco S, Meijer N, Ricci A. Food safety issues related to uses of insects for feeds and foods. Compr Rev Food Sci Food Saf. 17:1172-1183 2018;

54.

van der Fels-Klerx HJ, Camenzuli L, van der Lee MK, Oonincx DGAB. Uptake of cadmium, lead and arsenic by Tenebrio molitor and Hermetia illucens from contaminated substrates. PLoS One. 11:e0166186 2016;

55.

Yoon CH, Jeon SH, Ha YJ, Kim SW, Bang WY, Bang KH, Gal SW, Kim IS, Cho YS. Functional chemical components in Protaetia brevitarsis Larvae: Impact of supplementary feeds. Food Sci Anim Resour. 40:461-473 2020;

56.

Zhang ZQ, Chen SC, Xiao JH, Huang DW. State-of-the-art review of edible insect: From bioactives, pretreatment to enrichment. Food Biosci. 59:103879 2024;