1. Introduction
Current predictions estimate a world population of 9.7 billion by 2050 [
1]. In addition to this increase, prosperity in current developing countries is expected to grow, typified by a change in diet whereby more animal proteins are consumed. In order to produce animal proteins, high feed inputs are necessary. For example, currently 36% of all cultivated crops are used as animal feed [
2]. Proteins are an important component in these feeds and are mainly derived from oilseed meal (in particular soy), by-products of biodiesel production and fish [
3]. However, since 1970, overfishing has been a major issue and current fisheries exceed the maximum capacity for ensuring a sustainable biological system [
4]. This means a further depletion is detrimental for all fish stocks. Moreover, the sustainability of soy production can be questioned regarding the deforestation of rainforests [
5]. High demand for soy protein for both human food and livestock feed leads to excessively high production of soy in South-American countries. To meet the soy demands, deforestation (not only by cutting, but also by burning) happens rapidly to increase the area for soy cultivation. Strikingly, since 2021, this has led to the Amazon rainforests, previously known as the lungs of the earth, to become CO
2 emitters [
6,
7].
In order to thrive on an overcrowded planet in the long-term, it is important to relieve pressure on the most conventionally used protein sources such as fishmeal and soy and produce proteins sustainably [
8]. Insect farming, specifically the production of the larvae of
Hermetia illucens (Linnaeus), commonly known as the black soldier fly (BSF), has the potential to serve as a crucial component in the agro and feed industry [
9]. This is attributed not only to their exceptional bioconversion efficiency but also their remarkable digestion plasticity, which enables their growth using diverse organic side-streams [
10,
11]. It has already been shown that BSF larvae have a favourable feed conversion ratio, nutritional composition, and amino acid composition [
12,
13,
14]. The dry weight can contain up to 35% fat and 50% protein, with an amino acid profile similar to that of fishmeal [
15]. The side-streams that allow BSF rearing include various sources, including manure, waste-streams with high microbial loads, supermarket and domestic food waste, as well as agricultural side-streams that are unsuitable for human and animal consumption and therefore treated as organic waste [
10,
11]. Studies have shown that BSFs are extremely efficient in converting diverse sources of organic materials into larval biomass [
16,
17], showing reductions in organic waste up to 84.8% [
17,
18]. Organic waste constitutes 47% of the global total waste [
19]; moreover, roughly one-third of the global food production for human consumption (ca. 1.3 billion tons per year) is lost or wasted [
20]. Since BSFs can be used for the reduction of waste streams, as well as the production of high-quality protein, they contribute to increased circularity in our agro and feed industry, relieving the pressure on conventional protein sources [
11,
15].
Although the BSF industry has made significant progress in recent years, the use of waste or side-streams as a substrate for the production of BSFs remains largely unconsolidated. Moreover, the waste streams used are often already valorised as feed for conventional livestock, thus creating direct competition with this sector [
21]. However, the limited knowledge regarding the nutritional requirements of BSFs and the lack of a stable supply chain for waste streams as feed for BSFs, often makes producers reluctant to use them.
This is in sharp contrast to conventional livestock. For pig, cattle, and poultry, we possess detailed knowledge regarding the feed requirements for each developmental stage. This knowledge allows feed companies to successfully valorise side-streams as tailor-made feed for livestock. Even though the insect sector has made big advancements in the last decade, there are several studies that have evaluated growth of BSF larvae on specific waste streams [
11,
12,
22]. Unfortunately, due to large variations in experimental setup and composition of the side-streams (even with similar streams sampled at different time spans or locations), often varying results are obtained [
11]. Some of these studies have also attempted to model the nutritional requirements of BSF larvae, but modelling showed varying results due to the complexity of these streams [
23,
24,
25].
In order to optimise BSF production, it is essential to unravel how the substrate macronutrient composition influences larval growth and the larval nutritional composition [
23,
26]. However, modelling larval growth based on side-streams is demanding due to the variation in nutritional composition, including challenging-to-measure components such as different types of mono- and polysaccharides, and variations in amino acid and fatty acid profiles [
11]. Thus, to address these challenges, several studies used artificial substrates to overcome this issue [
23,
27].
Variations in amino acid profiles have a significant influence on larval growth, while the types of carbohydrates, including short sugars, long sugars, and different starch types, play a pivotal role in altering the digestive processes [
28]. Furthermore, diverse fat sources are reflected in larval fatty acid compositions [
28]. Notably, the fibre profile, encompassing lignin, cellulose, and hemicellulose contents, as well as the types of hemicelluloses, also significantly impact larval digestive processes [
29]. Some fibres can serve as an energy source, whereas others can negatively affect larval growth [
29].
Previous research has demonstrated the potential of artificial substrates in modelling the nutritional requirements for larvae. However, in previous studies, little attention was given to the physical composition of the substrates (e.g., adapted substrate dry matter contents based on their water holding capacity (WHC)) [
30]. This study takes a step forward from previous research by incorporating WHC methods developed by Frooninckx et al. (2024), utilising artificial substrates and employing a more intricate experimental design. This study aimed to develop a reliable model for accurately predicting larval growth parameters, specifically focusing on the larval mass produced when reared on various side-streams. By examining the effects of feed macronutrient contents on the growth of black soldier fly larvae, we aimed to provide valuable insights into optimising substrate formulations. This approach offers practical benefits for the industry by enabling the blending of nutritionally diverse side-streams to create ideal feed substrates, ultimately enhancing the efficiency and sustainability of BSF production systems.
2. Materials and Methods
2.1. Preparation of Feed Substrates
For composing artificial substrates, sunflower oil (Vandemoortele, Ghent, Belgium), wheat starch (Tarwezetmeel, Van Beekum specerijen, Harderwijk, The Netherlands), casein (organic casein protein, Ekopura, Hoofddorp, The Netherlands) and cellulose (Alphacel, MP Biomedicals, Santa Ana, CA, USA) were used. Composing artificial substrates, solely on cellulose, starch, casein, and sunflower oil, can potentially lead to deficits in essential micronutrients or substances such as sterols [
31]. Therefore, it was decided to use a base amount of chicken start mash (Chicken Start Mash 259, AVEVE, Geel, Belgium) of 8% dry matter (DM) for each different substrate composition. For the optimisation of the substrate dry matter content, potato starch was used instead of wheat starch. However, due to poor digestibility, for later experiments it was replaced with wheat starch. The macronutrient compositions of the ingredients are displayed in
Table 1 and were provided by the suppliers of the feed ingredient producers.
The experimental setup was divided into three phases:
Optimisation of substrate dry matter content (1): This phase aimed to determine the optimal dry matter content for the artificial substrates, ensuring larval growth.
Therefore, different hydration levels for a range of artificial substrate compositions were evaluated to identify the hydration condition that maximises performance or suitability for subsequent experiments.
Development of artificial substrates and prediction models (2): This phase involved the creation and testing of artificial substrates. Data from these experiments were used to develop predictive models that can forecast the performance of the substrates under various conditions.
Validation of the experimental design (3): This phase was dedicated to validating the accuracy and reliability of the developed models. This involved testing the models under different conditions with real side-streams to ensure their predictive validity and robustness.
2.1.1. Optimisation of Substrate Dry Matter Content
In the experiment to
optimise the substrate dry matter content, the design displayed in
Table 2 was used. Substrates with the displayed macronutrient contents were hydrated up to (1) a dry matter content of 30% and (2) their maximal water holding capacity (WHC). The amounts of each ingredient added to formulate the diet are displayed based on the dry matter content of the substrate.
WHC was determined by adding 3 g of fresh substrate to a 50 mL falcon and adding an excess amount of water (25 mL). The substrates underwent vortexing for homogenisation and were subsequently left at room temperature for 1 h. Following this, the falcons were subjected to centrifugation at 10,000×
g for a duration of 30 min. After centrifugation, the water layer was decanted, and the residual liquid was allowed to drain for 30 min. The substrate, along with the bound water, was then reweighed to determine the water holding capacity using the formula provided below [
30].
The WHC indicates the amount of water required to fully hydrate the substrate without any free water remaining. For example, a WHC of 150% means that 20 g of substrate would need 30 g of water to achieve optimal hydration.
2.1.2. Development of Artificial Substrates and Prediction Models
To formulate diverse feed substrates, a combination of a central composite and a Box–Behnken design methodology [
32] was used.
Figure 1 illustrates the arrangement of the Central Composite design (on the left) and the Box–Behnken design (on the right). Our aim was to assess the impact of substrate carbohydrate, protein, and fat content on larval bioconversion efficiency; thus, we set boundaries based on established substrates commonly used in black soldier fly (BSF) breeding [
11].
The macronutrient contents were categorised as (−α, −1, 0, +1, +α), wherein −1 represented the lower boundary, +1 the upper boundary, and 0 denoted the midpoint. Additionally, axial points (depicted as blue dots in
Figure 1) were introduced to the model to capture curvature, with values of −α and +α. We employed the Central Composite Circumscribed (CCC) methodology with an α value of 1.682 (
, with n as the number of factors (3), being substrate protein, fat, and carbohydrate content).
Table 3 displays the actual values corresponding to the theoretical points for each macronutrient content.
This approach allows for a comprehensive exploration of the effects of various substrate compositions on larval growth, contributing to our understanding of optimal feeding conditions for BSF larvae.
2.1.3. Preparation of Feed Substrates for Artificial Substrates
Table 4 shows both the macronutrient composition, and the formulation used to create the substrates, as used in the artificial substrate experiment. Conditions 1–33 were tested in 5-fold and condition 34 (the centrepoint) was tested in 10-fold, resulting in a total of 175 experimental units.
2.1.4. Validation of the Experimental Design
For the validation of the model, a mix of organic side-streams was used. A recent study focused on the growth of
H. illucens on 12 organic side-streams, which were chemically characterised, as presented in
Table 5 [
11]. The substrates consisted of pulp (pulp left behind after apple-juice production), beer draff (side-stream from beer-brewing), industrial food waste (supermarket and restaurant waste), chicken manure (mix of wood-pulp bedding and chicken manure), corn meal (ground corn), forced chicory roots (matured roots of chicory), fruit puree (fruit overproduction, mixed into a slurry), grain middlings (side-stream from wheat-industry), household food waste (organic food waste picked up from household containers, mixed into a slurry), hydrolysed feather meal (ground-up feathers from the chicken-industry, that underwent hydrolysation), and vegetable overproduction (overproduction from auctions, mixed into a slurry). The side-streams were selected on the basis of their availability in the Flanders region, Belgium. The side-streams were ground (Robot Coupe blixer 23, robot-coupe, Utrecht, The Netherlands) and frozen at −20 °C until used [
11].
For the validation experiment, the substrate mixtures were composed as presented in
Figure 2. Substrates were composed of 12%, 15.5%, and 30% crude protein content, allowing maximising inclusion of the different side-streams.
2.2. Black Soldier Fly Rearing and Maintenance
BSF larvae are continuously maintained by the Centre of Expertise in Sustainable Biomass and Chemistry at Thomas More University of Applied Sciences, Belgium at the Insect Pilot Plant. An egg quantity of 1 g was harvested (eggs deposited over the course of 48 h) and collected in a plastic weighing dish. This dish was placed upon a mixture of 100 g chicken feed (Chicken Start Mash 259, AVEVE, Belgium) and 100 mL tap water (total dry matter content of 45%). The container was incubated in a climate chamber at 27 ± 1 °C at 60% RH (relative humidity). On day 3 after collecting the eggs, the weighing dish was removed and the substrate including neonates was gently mixed using a table spoon. At day 4, the young larvae including substrate/frass mixture were transferred to a larger container containing 240 g of chicken feed and 360 g of water. For this experiment, 8-day-old larvae, calculated as days after harvest, were used at the start of the experiment (day 0).
2.3. Experimental Setup
For the optimisation of substrate dry matter content, 500 larvae (mean larval mass of 3–6 mg) were counted and separated from the nursing container using tweezers. In triplicate, plastic containers (17.5 × 11.9 × 5.9 cm), closed with a mesh lid, were filled with 100 g dry matter of artificial substrate to which the appropriate amount of water was added (targeting 30% DM or maximal WHC for artificial substrates) and 500 larvae were added to each of these containers. The amount of added substrate corresponds to 0.2 g DM per larva. A sheet of aluminium foil was placed on the substrate directly contacting approximately 80% of the substrate surface area, to prevent moisture loss through evaporation. At day 9 of the experiment, the larvae were separated from their frass, counted, weighed, and dry matter content was determined, by oven drying at 105 °C for 24 h.
For the development of artificial substrates and prediction models, 100 larvae (3–6 mg) were counted and separated from the nursing container using tweezers. In five-fold, plastic containers (6 cm diameter) were filled with 20 g dry matter of artificial substrate to which the appropriate amount of water was added (maximal WHC). To each of these containers, 100 larvae were added. The amount of added substrate corresponds to 0.2 g DM per larva. Also, a sheet of aluminium foil was placed on the substrate directly contacting approximately 80% of the substrate surface area, to prevent moisture loss through evaporation. At day 9 of the experiment, the larvae were separated from their frass, larvae were counted, weighed, and dry matter content was determined. Additionally, the total mass and dry matter content of the residues was determined. Larval survival was determined by counting larvae at the start and end of the feed experiments. At the end of the experiment, the total larval yield and residue were measured, and subsequently, the determination of dry matter content was performed (by oven drying at 105 °C for 24 h). This sequential approach allowed the calculation of larval growth parameters, such as bioconversion efficiency.
For the validation of the experimental design, 500 larvae (3–6 mg) were counted and separated from the nursing container using tweezers. In triplicate, plastic containers (17.5 × 11.9 × 5.9 cm) were filled with 100 g dry matter of various side-streams-mixes, and brought to maximal water holding capacity. The amount of added substrate corresponds to 0.2 g DM per larva. Again, aluminium foil was added on top of the substrate to prevent moisture loss through evaporation. At day 9 of the experiment, the larvae were separated from their frass, larvae were counted, weighed, and dry matter content was determined. Additionally, the total mass and dry matter content of the residues was determined. Larval survival was determined by counting larvae at the start and end of the feed experiments. At the end of the experiment, the total larval yield and residue were measured, and subsequently, the determination of dry matter content was performed (by oven drying at 105 °C for 24 h). This sequential approach allowed the calculation of larval growth parameters, such as bioconversion efficiency.
2.4. Calculations
The calculations displayed below were executed according to Broeckx et al., 2021 [
11]. The larval survival rates were determined through the division of the larval count at the end of the experiment (day 9) by the initial larval population count at the start of the experiment.
Bioconversion efficiency (BE) was calculated as following:
with L
end being the larval biomass at the end of the experiment and L
start being the larval biomass at the start of the experiment (both expressed in dry matter). D corresponds to the amount of added substrate (expressed in g dry matter).
2.5. Statistical Analysis
All statistical analysis and graphical illustrations were drafted using the JMP Pro17.0.0 software package from SAS (Buckinghamshire, UK). Data were tested for normality using a Shapiro–Wilk′s test and Levene′s test to examine the homogeneity of variance.
To determine significant differences in larval survival ratios, a t-test was applied with a significance threshold of p = 0.05. The prediction model was constructed in JMP. Least square linear regression was performed to investigate the relation between substrate macronutrient contents (independent variable) and bioconversion efficiency (dependent parameters). A full factorial approach was applied, including both main and interaction effects of the factors. Additionally, the quadratic effects of the factors (e.g., protein content2) were incorporated into the model.