Next Article in Journal
Anomaly Detection in Multi-Wavelength Photoplethysmography Using Lightweight Machine Learning Algorithms
Previous Article in Journal
A Review on Acoustic Emission Testing for Structural Health Monitoring of Polymer-Based Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Near Infrared Spectroscopy for Prediction of Yeast and Mould Counts in Black Soldier Fly Larvae, Feed and Frass: A Proof of Concept

by
Shanmugam Alagappan
1,2,
Anran Dong
3,
Deirdre Mikkelsen
1,3,
Louwrens C. Hoffman
1,2,4,
Sandra Milena Olarte Mantilla
1,
Peter James
5,
Olympia Yarger
6 and
Daniel Cozzolino
1,*
1
Centre for Nutrition and Food Sciences, Queensland Alliance for Agriculture and Food Innovation (QAAFI), The University of Queensland, Brisbane, QLD 4072, Australia
2
Fight Food Waste Cooperative Research Centre, Wine Innovation Central Building Level 1, Waite Campus, Urrbrae, SA 5064, Australia
3
School of Agriculture and Food Sustainability, Faculty of Science, University of Queensland, Brisbane, QLD 4072, Australia
4
Department of Animal Sciences, University of Stellenbosch, Private Bag X1, Matieland, Stellenbosch 7602, South Africa
5
Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation (QAAFI), The University of Queensland, Brisbane, QLD 4072, Australia
6
Goterra, 14 Arnott Street, Hume, Canberra, ACT 2620, Australia
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(15), 6946; https://doi.org/10.3390/s23156946
Submission received: 27 June 2023 / Revised: 1 August 2023 / Accepted: 2 August 2023 / Published: 4 August 2023
(This article belongs to the Section Optical Sensors)

Abstract

:
The use of black soldier fly larvae (BSFL) grown on different organic waste streams as a source of feed ingredient is becoming very popular in several regions across the globe. However, information about the easy-to-use methods to monitor the safety of BSFL is a major step limiting the commercialization of this source of protein. This study investigated the ability of near infrared (NIR) spectroscopy combined with chemometrics to predict yeast and mould counts (YMC) in the feed, larvae, and the residual frass. Partial least squares (PLS) regression was employed to predict the YMC in the feed, frass, and BSFL samples analyzed using NIR spectroscopy. The coefficient of determination in cross validation (R2CV) and the standard error in cross validation (SECV) obtained for the prediction of YMC for feed were (R2cv: 0.98 and SECV: 0.20), frass (R2cv: 0.81 and SECV: 0.90), larvae (R2cv: 0.91 and SECV: 0.27), and the combined set (R2cv: 0.74 and SECV: 0.82). However, the standard error of prediction (SEP) was considered moderate (range from 0.45 to 1.03). This study suggested that NIR spectroscopy could be utilized in commercial BSFL production facilities to monitor YMC in the feed and assist in the selection of suitable processing methods and control systems for either feed or larvae quality control.

1. Introduction

The world population is projected to grow up to 11 billion people by the end of this century, imposing a threat on global food security [1]. The growing population has also displayed changes in dietary patterns due to improved income [2,3]. These shifts in dietary habits and patterns have increased the global demand for protein, especially from animal sources. A recent report by the FAO has projected an increase (approx. 14%) in the demand and consumption of meat proteins by the year 2030 [4]. Patterns in consumption are also expected to increase up to 40% for poultry meat, followed by 34%, 20%, and 5% for pig meat, beef, and goat meat, respectively [4]. It has been recognized that livestock production by current practices accounts for about 18% of the total greenhouse gas emission derived from the agri-food production systems [5]. The production of meat protein using conventional feed ingredients is also challenged due to prevailing socio-environmental conditions, thereby suggesting the need for sustainable alternative feed sources [6,7,8].
The use of edible insects as a sustainable source of feed ingredient for livestock has been receiving increasing attention lately. The utilization of black solider fly larvae (BSFL), a saprophagous insect thriving in tropical and subtropical climatic conditions across the globe, has been of particular interest [9]. The BSFL undergoes six morphological stages or instars for a period of 10–52 days depending upon the feed source and rearing conditions, before they reach the pupal stage to become an adult fly [10].
Annually, 1.3 billion tons of the food generated for human consumption gets wasted [11]. Food waste constitutes for about 6% of the total anthropogenic greenhouse gas emissions [12]. The use of multi-stream food waste as substrate for rearing BSFL is becoming increasingly popular in commercial production facilities. BSFL reared on such food waste streams are rich in several nutrients including crude protein (range 32–46%), crude fat (range 15–38%), fibre (range 20.4–30.5%), and micronutrients [13,14,15,16]. Recent studies have also revealed that the replacement of conventional feed ingredients with BSFL meals in broilers, pigs, and aquaculture has positively influenced several growth performance parameters [17,18,19] making them a sustainable alternative to conventional feedstuff.
Feedstuffs and feed ingredients intended to be used as animal feed must be monitored for several microbial risks and regulated according to several legislations polices [20,21]. Yeast and mould counts (YMCs) are a combination of microbial organisms that are monitored in animal feeds in several regions across the globe [20,21]. The presence of these microorganisms can affect, among other things, the palatability of the feedstuff, growth, and the performance of domestic animals. It is also well established that several species of mould can produce mycotoxins, a secondary metabolite that can have adverse health effects in animals upon ingestion [22]. Previous studies have reported that BSFL reared on homogenous side stream wastes might vary in YMC between 2.0 to 7.0 log CFU/g depending upon the feed source and post-harvest conditions [23,24,25]. Therefore, it will be of importance for the feed industry to monitor YMC in BSFL when reared with different complex organic waste streams.
The detection and quantification of yeast and mould by conventional methods can be time consuming and laborious (e.g., culture of the microorganism in selective agars, microscopy, serology, histopathology, molecular markers, among others). On the other hand, the development of rapid and non-invasive methods will enable commercial BSFL production facilities to monitor microbial load as well as to assist in the implementation of a HACCP system or other feed safe protocols required by the industry [26]. Over the years, near infrared (NIR) spectroscopy combined with chemometric techniques has been showcased due to its potential to monitor product quality and ensure food safety in different fields [27,28,29,30]. NIR spectroscopy combined with chemometric techniques has been employed to determine the total viable counts and other microbial parameters in a diverse set of food samples (e.g., beef and pig meat, and fish) [31,32,33,34,35]. However, to date no research has been reported on the utilization of NIR spectroscopy to predict yeast or mould counts in BSFL.
The aim of this study is to evaluate the ability of NIR spectroscopy in combination with chemometric techniques to determine the YMC of feed, frass, untreated, and processed larvae reared with different food waste streams.

2. Materials and Methods

2.1. Sample Preparation

Five different feeds (2 homogenous/single food waste streams; 3 complex/heterogenous food waste streams) were used for the rearing of BSFL. The experiment was carried out using a homogeneous waste stream comprised of soy waste (SW) and customized bread-vegetable (BV) and the rearing conditions defined by the commercial site “A” were followed. The experiments carried out using the complex food waste streams included supermarket waste (food waste primarily comprising of bread, packed pastries, bottled milk, and chips) childcare center waste (primary constituents include salads, sandwiches, and cooked food) and a food waste mixture (mixture of waste from supermarkets, retail stores, childcare centres, and fast-food restaurants). The packaging material from these food wastes were removed and was used as feed for BSFL where the rearing conditions defined by the commercial site “B” were followed.

2.1.1. Experiment 1—Rearing of BSFL Using Homogenous Waste

The homogenous SW obtained from a custard manufacturer was utilized as feed for rearing BSFL. Customized BV waste was prepared by creating a mixture containing bread 60% (w/w), and a combination of peaches, potatoes, carrots, apples, and bananas (8% each w/w) purchased from a local supermarket. Two kg of SW was added to 6 trays (dimensions 60 × 40 × 12 cm) whilst the remaining feed required throughout the experiment was stored at 4 °C. A feed subsample was collected and stored at −20 °C prior to NIR analysis. A fraction of this feed was suspended in 40% (v/v) glycerol at a 1:1 (w/v) ratio and stored at −80 °C for YMC analysis.
Two kilograms of the starter feed from the BV experiment were prepared by chopping 1.2 kg of bread and 160 g of each of the fruits and vegetables. The chopped bread and vegetables were mixed and added to 6 trays (60 × 40 × 12 cm). Unlike SW, the feed utilized in this experiment was prepared fresh, where a subsample (approx. 20 g) was collected each time the BSFL in the trays were being fed and stored at −20 °C. The stored BV feed samples were pooled together before NIR analysis and stored at −20 °C. A sub-sample was suspended in 40% (v/v) glycerol at a 1:1 (w/v) ratio and stored at −80 °C for YMC analysis.
Approximately 10,000 5-day old (chronological age) BSFL obtained from hatched eggs were added into each of the 6 trays containing the starter feed (SW and BV feed). The trays were then subjected to the rearing conditions established by the industrial facility and fed ad libitum. The trays for both experiments were prepared for harvesting the 5th instar BSFL upon sighting the 1st pre-pupae. The larvae were picked using sterile tweezers from random spots in each of the trays (n = 6 trays for each waste type). The frass available upon harvesting the 5th instar larvae was collected from random spots in each of these trays. The 5th instar larvae and 5th instar frass samples were stored at −20 °C prior to NIR analysis where a fraction of these samples was suspended in 40% (v/v) glycerol at a 1:1 (w/v) ratio and stored at −80 °C for YMC analysis. The 6th instar larvae and 6th instar frass were harvested following the same protocol described above upon sighting of the first pupae and, stored at −20 °C prior to NIR analysis where a small sample fraction was suspended in 40% glycerol as stated above and stored at −80 °C for YMC analysis.

2.1.2. Experiment 2—Rearing of BSFL Using Heterogenous Waste

The three heterogenous food waste streams: supermarket waste, childcare centre waste, and a food waste mixture as described in Section 2.1 were received from local suppliers, de-packaged, and grinded using an industrial grinder. Three kg of the grinded waste from each of the waste streams were placed into trays (n = 3 for each waste, 60 × 40 × 12 cm) where approximately 15,000 of 5-day old BSFL (chronological age) were added to each tray. The trays were then subjected to the rearing conditions established by the industrial facility B and fed ad libitum. The feed for these experiments was de-packaged and grinded on the day of feeding where a subsample (approx. 20 g) from each of the waste streams was collected, pooled, and stored at −20 °C prior to NIR analysis. A fraction of the stored feed samples from the three different waste streams was suspended in 40% (v/v) glycerol at a 1:1 (w/v) ratio and stored at −80 °C for yeast and mould analysis. The 5th instar BSFL and frass from these trials were harvested and stored for analysis as mentioned above in Section 2.1.1. It is to be noted that 6th instar BSFL and frass were not collected in these experiments.
The temperature and other production conditions were monitored by each of the facilities. The feeding, waste management system, and rearing conditions are protected by industrial rights. Overall, this study has used industrial operating procedures and standards that can be used by similar types of industries worldwide.

2.1.3. Post-Harvest Processing of BSFL

The BSFL after harvesting was subjected to three different processing treatments. The whole live larvae samples were homogenized using a food processor (Cuisinart Food Processor, Sydney, NSW, Australia). The second treatment involved blanching the live larvae at 100 °C for 5 min and homogenizing them using a food processor (Cuisinart Food Processor, Australia). The final treatment involved blanching live larvae at the above-mentioned conditions, accompanied with drying using a commercial dehydrator at 70 °C for 12 h followed by homogenization with the food processor (Cuisinart Food Processor, Australia). After processing, the samples were stored at −20 °C prior to NIR analysis. A sub-sample of this was suspended in 40% (v/v) glycerol at a 1:1 (w/v) ratio and stored at −80 °C for YMC analysis. It is to be noted that blanching and drying treatments were carried out only on BSFL samples from Facility B.

2.2. Determination of Yeast and Mould

Samples (feed, larvae, and frass) collected from the five different waste streams described in the previous sections were analyzed using routine mould and yeast reference methods. The Australian standard: AS 5013.29-2009 was followed with some modifications [36]. In brief, 0.1 g of sample was suspended in 40% glycerol broth and homogenized with 0.9 mL peptone water (0.1%, w/v, Merck, Rahway, NJ, USA) and serial diluted to appropriated concentrations with the same diluent. Samples with water activity (aw) > 0.95 were plated onto Dichloran Rose Bengal Chloramphenicol (DRBC, Merck) agar plates where samples with aw < 0.95 were plated onto Dichloran 18% Glycerol (DG18, Merck) agar plates. The Aw of each sample was measured using a Novasina LabSwift-AW water activity meter (Novasina AG, Lachen, Switzerland). All plates were aerobically incubated with the lid in uppermost position at 25 °C for 5 days before enumeration. The Penicillium commune (laboratory strain) and Saccharomyces cerevisiae (laboratory strain) were streaked on both DRBC and DG18 agar plates and used as positive controls. Both plates were incubated under the same incubation conditions as described above. All the analysis was carried out in biological triplicates.

2.3. NIR Analysis

Stored samples were thawed at room temperature and analyzed using a Fourier Transform-NIR instrument (Tango-R, Bruker Optics GmbH, Ettlingen, Germany). Samples were placed in a glass cuvette (10 mm diameter) where the spectra were collected as the average of 64 interferograms at a resolution of 4 cm−1 in the wavenumber range of 11,550 to 3950 cm−1 (OPUS software, version 8.5, Bruker Optics GmbH, Ettlingen, Germany). The cuvettes were cleaned with 70% ethanol and wiped dry using paper towels between samples. The NIR raw spectra, baseline corrected spectra, first (1, 10, 10, 1), and second (2, 10, 10, 2) derivatives are shown in the Supplementary File.

2.4. Data Analysis

The Vektor Direktor (version 1.0, KAX Group, Sydney, NSW, Australia) was used for multivariate data analysis. The NIR spectral data were smoothed and pre-processed using the Savitzky–Golay second derivative (second order polynomial and a smoothing window size of 10 points) before analysis. Trends, patterns, and outliers in the data set were visualized by performing principal component analysis (PCA) [37]. Partial least squares regression (PLS) [38] was utilized to develop models for the prediction of YMC in the feed, frass, and BSFL samples obtained from the different experiments [39]. A summary of the different data sets (number of samples) used to develop the PLS models is depicted in Table 1. To evaluate the performance of the PLS models, the original dataset was split into two subsets, namely calibration and validation, using the Kennard-Stone algorithm provided by the software (Vektor Direktor, version 1.0, KAX Group, Sydney, NSW, Australia). In this study, approximately 70% of the samples were selected and used to develop the cross-validation models, while 30% of the samples were used for validation. By performing data partitioning, knowledge of the training dataset did not affect the test dataset, and the predictive power of the created model subsequently increased. A leave-one-out cross validation was applied during the development of the PCA and PLS models. The cross-validation models were evaluated using the coefficient of determination (R2CV), the standard error in cross validation (SECV), bias, slope, and the residual predictive deviation (RPD) (SD/SECV) [40,41]. The validation models were evaluated using the standard error in prediction [40,41].

3. Results and Discussion

3.1. Yeast and Mould Counts

The average and standard deviation of the YMC determined in each of the samples is shown in Figure 1. It can be observed that the moulds were distributed at higher concentrations in frass than in larvae and feed. It has been reported that BSFL gut is characterized with a diverse set of fungal communities including Candida, Dipodascaceae, Trichosporon, Pichia, species where their occurrence depends upon the substrate used for rearing the larvae [42,43,44]. Recent studies have also reported the passage of different genera of bacteria present in the guts of BSFL through their digestive tracts as consequence of the digestion of different waste streams or feeds used to grow the larvae [45]. Therefore, it can be suggested that fungal communities could also occur in the residual frass after the digestion process. However, the fungal count in the supermarket waste frass was found to be lower than that in the waste source. This can possibly be accounted for due to the relatively low water activity (0.85) observed in the supermarket frass compared to the other frass samples (>0.91).
In this study, the YMC in the frass obtained from the 6th instars samples from site A was higher than the frass obtained from the 5th instar (Figure 1). It has been reported that BSFL upon reaching 6th instar (pre-pupae) will stop feeding on the substrate, empty their digestive tract and crawl off to a dry clean site to commence pupation [46]. The emptying of the digestive tract during the pre-pupal stage could explain the high counts in the frass obtained from the 6th instar frass compared to the 5th instar. Furthermore, the 6th instar BSFL, despite having emptied their guts, was found to have higher YMC when compared to the 5th instar larvae when reared with both SW and BV waste. It is to be noted that the larvae samples in these experiments were not surface sterilized upon harvest. Therefore, it is suggested that the high counts observed in 6th instar larvae is probably due to cross-contamination of the exterior of the larvae that is in contact with the frass.
The YMC for the larvae reared on heterogenous waste streams was observed to be low when compared to that of BSFL reared on homogenous waste streams (SW and BV waste). It is well established that the differences in the rearing substrate and rearing conditions modulate the nutritional composition of the larvae and microbial counts [47,48,49]. The BSFL has also been reported to exert antimicrobial action against pathogenic microbes based on the feed source and rearing conditions employed [50,51,52].
It can be also seen from Figure 1 that the YMC of BSFL larvae subjected to processing, such as blanching and drying, were reduced to 1.3 log CFU/g. Campbell and collaborators [23] reported similar levels of reduction in YMC upon subjecting the BSFL to two different thermal processing treatments. The ability of thermal processing to reduce the microbial load of several pathogenic microorganisms in BSFL has been demonstrated in previous studies [53,54].

3.2. Principal Component Analysis

A principal component analysis (PCA) was performed using the NIR spectra collected from all the experiments (n = 75). The first two principal components (PC) explain 88.34% of variance in the samples. It can be observed from the PCA scores plot (Figure 2) that the frass and larvae samples from the different rearing experiments were well separated and formed separate clusters unlike the feed samples which exhibited limited separation. It is also noteworthy that the BSFL samples subjected to the drying process were found to be clustered and separated from the rest of the larvae samples. The PCA loading plot showed that the highest loadings in PC1 had characteristic bands in the region around 7200 cm−1 and 5300 cm−1 which is associated with the water content [55,56] (Figure 3). The bands observed in the region around 5800 cm−1 and between 4900 cm−1 to 4600 cm−1 could be attributed to lipid and protein content, respectively [57,58,59]. The highest loadings in PC2 exhibited positional shifts in the bands around 7100 cm−1 and 5300 cm−1 which could be associated with changes in water activity of the sample. Similar positional shifts were observed in regions around 4900 cm−1 and 4600 cm−1 (Figure 3). These shifts could be attributed to the difference in physical characteristics of the samples and explained by processing and other factors (e.g., site, waste stream, etc.) [26].

3.3. Descriptive and Cross Validation Statistics

The descriptive statistics (mean, standard deviation, range, and CV) for the YMC measured in the different data sets used to develop the PLS models are shown in Table 2. The mean YMC values range from 2.1 log CFU/g in the larvae samples to 4.6 log CFU/g in the frass samples. The CV for the YMC values showed that the complete data set has the highest variability where the lowest was found to be in the frass samples. Overall, the descriptive statistics indicated that there was good dispersion in the samples to attempt PLS calibration for YMC using NIR spectroscopy.
The cross-validation and prediction statistics for the measurement of YMC in the different data sets (larvae, feed, and frass) from the different waste streams and sites are reported in Table 3. The coefficient of determination in cross validation (R2CV) and the standard error in cross validation (SECV) obtained for the prediction of YMC were for feed (R2cv: 0.98 and SECV: 0.20), frass (R2cv: 0.81 and SECV: 0.90), larvae (R2cv: 0.91 and SECV: 0.27). The RPD value (4) obtained for the prediction of YMC in the feed suggested that this model could be used to predict YM counts in the feed sources intended to be used for feeding the BSFL. The RPD values for the prediction of YMC in the frass and larvae were 3.1 and 4, respectively. The cross-validation statistics using samples from facility A were R2cv: 0.90, SECV: 0.45, and RPD: 3.5 while for facility B they were R2cv: 0.80, SECV: 0.67, and RPD: 2. These cross-validation models were considered moderately accurate, suggesting that they could be used for screening the samples for the occurrence of YMC in the samples. The cross-validation statistics for the models developed using all samples (combined set) were R2cv: 0.74, SECV: 0.82, and RPD: 2.1. The predictive ability of the cross-validation models was evaluated using the validation set created by applying the Kennard-Stone algorithm explained in the Section 2. Moderate predictions of YMC in all the samples were obtained as indicated by the moderate standard error of prediction (SEP) obtained (range 0.45 to 1.03). The predicted vs. reference values plot for the different cross-validation models are displayed in Figure 4A–F. The number of latent variables (LV) used to develop the cross-validation models were 2 for feed and frass, 5 for facility A and the combined set, 6 LV were used to develop the models using the samples from facility B, and 8 for larvae.
Previous results from our group indicated that BSFL can be cluster based on it its morphological stage and the feed source used to rear the larvae [26]. It was discussed in the previous section that the feed, frass, and the larvae samples employed in this study were clustered in groups, suggesting spectral differences among the samples. It can also be suggested that the poor cross-validation statistics obtained for some of the models developed could be related to some of the limitations of the reference method. Furthermore, some of the data sets have shown a very low concentration of YMC.
The loadings plot for the PLS models used to predict YMC in the feed, larvae, and frass samples are reported in Figure 5A–C. The loadings derived from the PLS regression models were found to be different with few similarities between the feed and frass samples analyzed. The highest and most similar loadings for the prediction of YMC in the feed and frass samples were observed in the region around 5100 cm−1 which is associated with O-H stretching, HOH, and N-H bending combinations that can be associated with water and proteins, whilst that around 5800 cm−1 is associated with asymmetric and symmetric stretching of C-H in CH2 in the first overtone region, which is mainly related to lipids. The bands around 7100 cm−1 are associated with O-H stretching and bending vibrations in the first overtone attributed to the water content in the samples [55,56] and, the band around 8700 cm−1 is due to the C-H vibrations in the second overtone region that can be linked to either lipids or aromatic groups [57,60].
In addition, the feed samples showed the highest loadings around 5300 cm−1 which is associated with moisture content of the samples, and around 5968 cm−1 which is linked with C-H aromatic groups while around 6928 cm−1 with N-H asymmetric vibrations in the 1st overtone region—which are attributed to proteins [55,56]. The bands around 5968 cm−1 and 6928 cm−1 are also associated with C-H (2υ) vibrations in aromatic groups [58,60]. The highest loadings for the prediction of YMC using the frass samples were observed around 4352 cm−1 and 4784 cm−1, associated with C-H bending (aromatic groups) and O-H, respectively [58,59,61] whilst the band around 6032 cm−1 corresponds with C-H vibrations in the 1st overtone region attributed to lipids [59,60]. The highest loadings used to predict YMC using the larvae samples were observed around the regions 4848 cm−1, 5088 cm−1, 5952 cm−1, and 6800 cm−1. The band at 4848 cm−1 is associated with C-H bending (aromatic groups), around 5952 cm−1 with C-H (2υ) vibrations in the aromatic groups, while around 6800 cm−1 with protein [57,58,60]. The band at 5088 cm−1 is associated with O-H and C-H combinations related to polysaccharides [60].

4. Conclusions

The findings of this study indicated that NIR spectroscopy combined with PLS regression was able to predict YMC with good and moderate accuracy in the feed samples and residual frass. The cross-validation statistics for the prediction of YMC in the larvae, despite exhibiting different spectral characteristics, did not yield prediction models suitable for screening fungal contamination. The reference data imposed certain restraints that influenced the accuracy of the PLS model developed for the larvae samples. The use of NIR spectroscopy in combination with chemometrics can be used as an instantaneous tool to monitor the quality of the feed samples intended to be used as substrate for rearing BSFL.
Overall, this study suggested that NIR spectroscopy could be utilized in commercial BSFL production facilities to monitor YMC in the feed, assist in the selection of suitable processing methods to be applied based on the extent of contamination observed in the substrate, and implement control systems either for feed or larvae quality control.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s23156946/s1, Figure S1: Near infrared raw spectra of all samples (feed, frass and larvae) analysed.; Figure S2. Near infrared baseline corrected spectra of all samples (feed, frass and larvae) analysed. Figure S3. First derivative of the NIR spectra of all samples (feed, frass and larvae) analysed. Figure S4. Second derivative of the NIR spectra of all samples (feed, frass and larvae) analysed.

Author Contributions

S.A. and A.D., sampling, data collection, analysis, data interpretation, draft preparation; D.C., sample preparation, analysis, draft preparation; L.C.H., D.M., S.M.O.M., P.J., O.Y., D.C., A.D. and S.A., writing—review and editing of manuscript; L.C.H., D.M., S.M.O.M., P.J., O.Y. and D.C., project administration; L.C.H.: funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

The work has been supported by the Fight Food Waste Cooperative Research Centre (Project 2.4.1), whose activities are funded by the Australian Government’s Cooperative Research Centre Program.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

The authors acknowledge Goterra Pvt Ltd. and Flyfarm Pvt Ltd. for their support with sample preparation and collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Erhard, A.L.; Águas Silva, M.; Damsbo-Svendsen, M.; Menadeva Karpantschof, B.-E.; Sørensen, H.; Bom Frøst, M. Acceptance of insect foods among Danish children: Effects of information provision, food neophobia, disgust sensitivity, and species on willingness to try. Food Qual. Prefer. 2023, 104, 104713. [Google Scholar] [CrossRef]
  2. Ismail, B.P.; Senaratne-Lenagala, L.; Stube, A.; Brackenridge, A. Protein demand: Review of plant and animal proteins used in alternative protein product development and production. Anim. Front. 2020, 10, 53–63. [Google Scholar] [CrossRef] [PubMed]
  3. Kim, S.W.; Less, J.F.; Wang, L.; Yan, T.; Kiron, V.; Kaushik, S.J.; Lei, X.G. Meeting Global Feed Protein Demand: Challenge, Opportunity, and Strategy. Annu. Rev. Anim. Biosci. 2019, 7, 221–243. [Google Scholar] [CrossRef] [PubMed]
  4. OECD; Food and Agriculture Organization of the United Nations. OECD-FAO Agricultural Outlook 2021–2030; OECD Publishing: Paris, France, 2021. [Google Scholar]
  5. FAO. Greenhouse Gas Emissions from Agrifood Systems: Global, Regional and Country Trends, 2000–2020; FAOSTAT Analytical Brief No. 50; FAO: Rome, Italy, 2022; p. 12. [Google Scholar]
  6. Kim, T.K.; Yong, H.I.; Kim, Y.-B.; Kim, H.-W.; Choi, Y.-S. Edible Insects as a Protein Source: A Review of Public Perception, Processing Technology, and Research Trends. Food Sci. Anim. Resour. 2019, 39, 521–540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kolasa-Więcek, A. Modeling greenhouse gas emissions from livestock farming in Poland with the use of stepwise multiple regression. J. Res. Appl. Agric. Eng. 2013, 58, 78–85. [Google Scholar]
  8. van Huis, A. Insects as Food in sub-Saharan Africa. Int. J. Trop. Insect Sci. 2003, 23, 163–185. [Google Scholar] [CrossRef]
  9. Kaya, C.; Generalovic, T.N.; Ståhls, G.; Hauser, M.; Samayoa, A.C.; Nunes-Silva, C.G.; Roxburgh, H.; Wohlfahrt, J.; Ewusie, E.A.; Kenis, M.; et al. Global population genetic structure and demographic trajectories of the black soldier fly, Hermetia illucens. BMC Biol. 2021, 19, 94. [Google Scholar] [CrossRef]
  10. Liew, C.S.; Yunus, N.M.; Chidi, B.S.; Lam, M.K.; Goh, P.S.; Mohamad, M.; Sin, J.C.; Lam, S.M.; Lim, J.W.; Lam, S.S. A review on recent disposal of hazardous sewage sludge via anaerobic digestion and novel composting. J. Hazard. Mater. 2022, 423, 126995. [Google Scholar] [CrossRef]
  11. Fowles, T.M.; Nansen, C. Insect-based bioconversion: Value from food waste. In Food Waste Management; Springer: Berlin/Heidelberg, Germany, 2020; pp. 321–346. [Google Scholar]
  12. Amicarelli, V.; Lagioia, G.; Bux, C. Global warming potential of food waste through the life cycle assessment: An analytical review. Environ. Impact Assess. Rev. 2021, 91, 106677. [Google Scholar] [CrossRef]
  13. Shumo, M.; Osuga, I.M.; Khamis, F.M.; Tanga, C.M.; Fiaboe, K.K.M.; Subramanian, S.; Ekesi, S.; van Huis, A.; Borgemeister, C. The nutritive value of black soldier fly larvae reared on common organic waste streams in Kenya. Sci. Rep. 2019, 9, 10110. [Google Scholar] [CrossRef] [Green Version]
  14. Nyakeri, E.; Ogola, H.; Ayieko, M.; Amimo, F. Valorisation of organic waste material: Growth performance of wild black soldier fly larvae (Hermetia illucens) reared on different organic wastes. J. Insects Food Feed. 2017, 3, 193–202. [Google Scholar] [CrossRef]
  15. Gao, Z.; Wang, W.; Lu, X.; Zhu, F.; Liu, W.; Wang, X.; Lei, C. Bioconversion performance and life table of black soldier fly (Hermetia illucens) on fermented maize straw. J. Clean. Prod. 2019, 230, 974–980. [Google Scholar] [CrossRef]
  16. Gold, M.; Cassar, C.M.; Zurbrügg, C.; Kreuzer, M.; Boulos, S.; Diener, S.; Mathys, A. Biowaste treatment with black soldier fly larvae: Increasing performance through the formulation of biowastes based on protein and carbohydrates. Waste Manag. 2020, 102, 319–329. [Google Scholar] [CrossRef]
  17. Phaengphairee, P.; Boontiam, W.; Wealleans, A.; Hong, J.; Kim, Y.Y. Dietary supplementation with full-fat Hermetia illucens larvae and multi-probiotics, as a substitute for antibiotics, improves the growth performance, gut health, and antioxidative capacity of weaned pigs. BMC Vet. Res. 2023, 19, 7. [Google Scholar] [CrossRef] [PubMed]
  18. Nampijja, Z.; Kiggundu, M.; Kigozi, A.; Lugya, A.; Magala, H.; Ssepuuya, G.; Nakimbugwe, D.; Walusimbi, S.S.; Mugerwa, S. Optimal substitution of black soldier fly larvae for fish in broiler chicken diets. Sci. Afr. 2023, 20, e01636. [Google Scholar] [CrossRef]
  19. Maranga, B.; Kagali, R.; Mbogo, K.; Orina, P.; Munguti, J.; Ogello, E. Growth Performance of African Catfish (Clarias gariepinus) Fed on Diets Containing Black Soldier Fly (Hermetia illucens) Larvae Under Aquaponic System. Aquac. Stud. 2022, 23. Available online: https://www.aquast.org/uploads/pdf_565.pdf (accessed on 1 June 2023). [CrossRef]
  20. FAO; IFIF. Good Practices for the Feed Sector—Implementing the Codex Alimentarius Code of Practice on Good Animal Feeding; FAO: Rome, Italy, 2020; p. 140. [Google Scholar]
  21. Dairy Australia. Mould and Mycotoxin Risks in Feed; Dairy Australia: Melbourne, VIC, Australia, 2022. [Google Scholar]
  22. Zain, M.E. Impact of mycotoxins on humans and animals. J. Saudi Chem. Soc. 2011, 15, 129–144. [Google Scholar] [CrossRef] [Green Version]
  23. Campbell, M.; Ortuño, J.; Stratakos, A.C.; Linton, M.; Corcionivoschi, N.; Elliott, T.; Koidis, A.; Theodoridou, K. Impact of Thermal and High-Pressure Treatments on the Microbiological Quality and In Vitro Digestibility of Black Soldier Fly (Hermetia illucens) Larvae. Animals 2020, 10, 682. [Google Scholar] [CrossRef] [Green Version]
  24. Kamau, E.; Mutungi, C.; Kinyuru, J.; Imathiu, S.; Affognon, H.; Ekesi, S.; Nakimbugwe, D.; Fiaboe, K.K.M. Changes in chemical and microbiological quality of semi-processed black soldier fly (Hermetia illucens L.) larval meal during storage. J. Insects Food Feed. 2020, 6, 417–428. [Google Scholar] [CrossRef]
  25. Were, G.J.; Irungu, F.G.; Ngoda, P.N.; Affognon, H.; Ekesi, S.; Nakimbugwe, D.; Fiaboe, K.K.M.; Mutungi, C.M. Nutritional and microbial quality of extruded fish feeds containing black soldier fly (Hermetia illucens L.) larvae meal as a replacement for fish meal for Nile Tilapia (Oreochromis niloticus) and African sharptooth catfish (Clarius gariepinus). J. Appl. Aquac. 2022, 34, 1036–1052. [Google Scholar] [CrossRef]
  26. Alagappan, S.; Hoffman, L.C.; Mantilla, S.M.O.; Mikkelsen, D.; James, P.; Yarger, O.; Cozzolino, D. Near Infrared spectroscopy as a traceability tool to monitor black soldier fly larvae (Hermetia illucens) intended as animal feed. Appl. Sci. 2022, 12, 8168. [Google Scholar] [CrossRef]
  27. Wang, X. 7—Near-infrared spectroscopy for food quality evaluation. In Evaluation Technologies for Food Quality; Zhong, J., Wang, X., Eds.; Woodhead Publishing: Sawston, UK, 2019; pp. 105–118. [Google Scholar]
  28. Cortés, V.; Blasco, J.; Aleixos, N.; Cubero, S.; Talens, P. Monitoring strategies for quality control of agricultural products using visible and near-infrared spectroscopy: A review. Trends Food Sci. Technol. 2019, 85, 138–148. [Google Scholar] [CrossRef]
  29. Pasquini, C. Near infrared spectroscopy: A mature analytical technique with new perspectives–A review. Anal. Chim. Acta 2018, 1026, 8–36. [Google Scholar] [PubMed]
  30. Cozzolino, D. 7—Near Infrared Spectroscopy and Food Authenticity. In Advances in Food Traceability Techniques and Technologies; Espiñeira, M., Santaclara, F.J., Eds.; Woodhead Publishing: Sawston, UK, 2016; pp. 119–136. [Google Scholar]
  31. Achata, E.M.; Oliveira, M.; Esquerre, C.A.; Tiwari, B.K.; O’Donnell, C.P. Visible and NIR hyperspectral imaging and chemometrics for prediction of microbial quality of beef Longissimus dorsi muscle under simulated normal and abuse storage conditions. LWT 2020, 128, 109463. [Google Scholar] [CrossRef]
  32. Barbin, D.F.; Elmasry, G.; Sun, D.-W.; Allen, P.; Morsy, N. Non-destructive assessment of microbial contamination in porcine meat using NIR hyperspectral imaging. Innov. Food Sci. Emerg. Technol. 2013, 17, 180–191. [Google Scholar] [CrossRef]
  33. Cheng, J.-H.; Sun, D.-W. Rapid and non-invasive detection of fish microbial spoilage by visible and near infrared hyperspectral imaging and multivariate analysis. LWT-Food Sci. Technol. 2015, 62, 1060–1068. [Google Scholar] [CrossRef]
  34. Liu, Y.; Dixit, Y.; Reis, M.M.; Prabakar, S. Towards the non-invasive assessment of staling in bovine hides with hyperspectral imaging. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2023, 289, 122220. [Google Scholar] [CrossRef]
  35. Tito, N.B.; Rodemann, T.; Powell, S.M. Use of near infrared spectroscopy to predict microbial numbers on Atlantic salmon. Food Microbiol. 2012, 32, 431–436. [Google Scholar] [CrossRef]
  36. AS 5013.29-2009; Standards Australia Committee FT-024 Food Products, Subcommittee FT-024-01 Food Microbiology. Method 29: Examination for Specific Organisms- Colony Count of Yeasts and Moulds. Standards Australia: Sydney, NSW, Australia, 2009.
  37. Cozzolino, D.; Power, A.; Chapman, J. Interpreting and Reporting Principal Component Analysis in Food Science Analysis and Beyond. Food Anal. Methods 2019, 12, 2469–2473. [Google Scholar] [CrossRef]
  38. Wold, S.; Sjöström, M.; Eriksson, L. PLS-regression: A basic tool of chemometrics. Chemom. Intell. Lab. Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  39. Cozzolino, D. The Ability of Near Infrared (NIR) Spectroscopy to Predict Functional Properties in Foods: Challenges and Opportunities. Molecules 2021, 26, 6981. [Google Scholar] [CrossRef] [PubMed]
  40. Williams, P.; Dardenne, P.; Flinn, P. Tutorial: Items to be included in a report on a near infrared spectroscopy project. J. Near Infrared Spectrosc. 2017, 25, 85–90. [Google Scholar] [CrossRef]
  41. Fearn, T. Assessing calibrations: SEP, RPD, RER and R2. NIR News 2002, 13, 12–13. [Google Scholar] [CrossRef]
  42. Vitenberg, T.; Opatovsky, I. Assessing Fungal Diversity and Abundance in the Black Soldier Fly and its Environment. J. Insect Sci. 2022, 22, 3. [Google Scholar] [CrossRef]
  43. Varotto Boccazzi, I.; Ottoboni, M.; Martin, E.; Comandatore, F.; Vallone, L.; Spranghers, T.; Eeckhout, M.; Mereghetti, V.; Pinotti, L.; Epis, S. A survey of the mycobiota associated with larvae of the black soldier fly (Hermetia illucens) reared for feed production. PLoS ONE 2017, 12, e0182533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Tanga, C.M.; Waweru, J.W.; Tola, Y.H.; Onyoni, A.A.; Khamis, F.M.; Ekesi, S.; Paredes, J.C. Organic Waste Substrates Induce Important Shifts in Gut Microbiota of Black Soldier Fly (Hermetia illucens L.): Coexistence of Conserved, Variable, and Potential Pathogenic Microbes. Front. Microbiol. 2021, 12, 635881. [Google Scholar] [CrossRef]
  45. Gold, M.; von Allmen, F.; Zurbrügg, C.; Zhang, J.; Mathys, A. Identification of Bacteria in Two Food Waste Black Soldier Fly Larvae Rearing Residues. Front. Microbiol. 2020, 11, 582867. [Google Scholar] [CrossRef] [PubMed]
  46. Soetemans, L.; Uyttebroek, M.; Bastiaens, L. Characteristics of chitin extracted from black soldier fly in different life stages. Int. J. Biol. Macromol. 2020, 165, 3206–3214. [Google Scholar] [CrossRef]
  47. Alagappan, S.; Rowland, D.; Barwell, R.; Mantilla, S.; Mikkelsen, D.; James, P.; Yarger, O.; Hoffman, L. Legislative landscape of black soldier fly (Hermetia illucens) as feed. J. Insects Food Feed. 2021, 8, 343–355. [Google Scholar] [CrossRef]
  48. Purkayastha, D.; Sarkar, S. Sustainable waste management using black soldier fly larva: A review. Int. J. Environ. Sci. Technol. 2021, 19, 12701–12726. [Google Scholar] [CrossRef]
  49. Seyedalmoosavi, M.M.; Mielenz, M.; Veldkamp, T.; Daş, G.; Metges, C.C. Growth efficiency, intestinal biology, and nutrient utilization and requirements of black soldier fly (Hermetia illucens) larvae compared to monogastric livestock species: A review. J. Anim. Sci. Biotechnol. 2022, 13, 31. [Google Scholar] [CrossRef] [PubMed]
  50. Lalander, C.H.; Fidjeland, J.; Diener, S.; Eriksson, S.; Vinnerås, B. High waste-to-biomass conversion and efficient Salmonella spp. reduction using black soldier fly for waste recycling. Agron. Sustain. Dev. 2015, 35, 261–271. [Google Scholar] [CrossRef] [Green Version]
  51. De Smet, J.; Vandeweyer, D.; Van Moll, L.; Lachi, D.; Van Campenhout, L. Dynamics of Salmonella inoculated during rearing of black soldier fly larvae (Hermetia illucens). Food Res. Int. 2021, 149, 110692. [Google Scholar] [CrossRef] [PubMed]
  52. Erickson, M.C.; Islam, M.; Sheppard, C.; Liao, J.; Doyle, M.P. Reduction of Escherichia coli O157: H7 and Salmonella enterica serovar enteritidis in chicken manure by larvae of the black soldier fly. J. Food Prot. 2004, 67, 685–690. [Google Scholar] [CrossRef] [Green Version]
  53. Bessa, L.W.; Pieterse, E.; Marais, J.; Dhanani, K.; Hoffman, L.C. Food Safety of Consuming Black Soldier Fly (Hermetia illucens) Larvae: Microbial, Heavy Metal and Cross-Reactive Allergen Risks. Foods 2021, 10, 1934. [Google Scholar] [CrossRef]
  54. Larouche, J.; Deschamps, M.-H.; Saucier, L.; Lebeuf, Y.; Doyen, A.; Vandenberg, G.W. Effects of Killing Methods on Lipid Oxidation, Colour and Microbial Load of Black Soldier Fly (Hermetia illucens) Larvae. Animals 2019, 9, 182. [Google Scholar] [CrossRef] [Green Version]
  55. Santos, P.M.; Simeone, M.L.F.; Pimentel, M.A.G.; Sena, M.M. Non-destructive screening method for detecting the presence of insects in sorghum grains using near infrared spectroscopy and discriminant analysis. Microchem. J. 2019, 149, 104057. [Google Scholar] [CrossRef]
  56. Wenz, J.J. Examining water in model membranes by near infrared spectroscopy and multivariate analysis. Biochim. Biophys. Acta (BBA)-Biomembr. 2018, 1860, 673–682. [Google Scholar] [CrossRef]
  57. Kröncke, N.; Böschen, V.; Woyzichovski, J.; Demtröder, S.; Benning, R. Comparison of suitable drying processes for mealworms (Tenebrio molitor). Innov. Food Sci. Emerg. Technol. 2018, 50, 20–25. [Google Scholar] [CrossRef]
  58. Benes, E.; Biró, B.; Fodor, M.; Gere, A. Analysis of wheat flour-insect powder mixtures based on their near infrared spectra. Food Chem. X 2022, 13, 100266. [Google Scholar] [CrossRef]
  59. Cruz-Tirado, J.P.; Amigo, J.M.; Barbin, D.F. Determination of protein content in single black fly soldier (Hermetia illucens L.) larvae by near infrared hyperspectral imaging (NIR-HSI) and chemometrics. Food Control 2023, 143, 109266. [Google Scholar] [CrossRef]
  60. Weyer, L. Practical Guide to Interpretive Near-Infrared Spectroscopy; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  61. Riu, J.; Vega, A.; Boqué, R.; Giussani, B. Exploring the Analytical Complexities in Insect Powder Analysis Using Miniaturized NIR Spectroscopy. Foods 2022, 11, 3524. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Yeast and mould count of BSFL, feed, and frass samples employed in the study.
Figure 1. Yeast and mould count of BSFL, feed, and frass samples employed in the study.
Sensors 23 06946 g001
Figure 2. Principal component score plot of all samples (data sets) analyzed using near infrared reflectance spectroscopy.
Figure 2. Principal component score plot of all samples (data sets) analyzed using near infrared reflectance spectroscopy.
Sensors 23 06946 g002
Figure 3. Principal component loadings plot of all samples (data sets) analyzed using near infrared reflectance spectroscopy.
Figure 3. Principal component loadings plot of all samples (data sets) analyzed using near infrared reflectance spectroscopy.
Sensors 23 06946 g003
Figure 4. Predicted vs. reference plot for the prediction of yeast and mould counts in the different data sets analyzed using near infrared reflectance spectroscopy.
Figure 4. Predicted vs. reference plot for the prediction of yeast and mould counts in the different data sets analyzed using near infrared reflectance spectroscopy.
Sensors 23 06946 g004
Figure 5. Partial least squares loadings for the optimal models used to predict yeast and mould counts in the feed, larvae, and frass samples analyzed using near infrared reflectance spectroscopy.
Figure 5. Partial least squares loadings for the optimal models used to predict yeast and mould counts in the feed, larvae, and frass samples analyzed using near infrared reflectance spectroscopy.
Sensors 23 06946 g005
Table 1. Summary of different datasets used to develop the partial least square regression models.
Table 1. Summary of different datasets used to develop the partial least square regression models.
Data SetDetails of Samples
Feed (n = 15)Substrate used for rearing BSFL: SW, customised BV, supermarket waste, childcare waste, waste mixture.
Larvae (n = 39)5th and 6th instar BSFL from SW and BV waste; unprocessed, blanched, and dried larvae from supermarket waste, childcare waste, waste mixture.
Frass (n = 21)5 and 6th instar frass from soy and BV waste trials. 5th instar frass from supermarket waste, childcare waste, waste mixture
Facility A (n = 30)Feed, 5th and 6th instar BSFL and frass from soy waste and BV waste. Homogenous data set.
Facility B (n = 45)Feed, unprocessed, blanched, dried BSFL, and frass from supermarket waste, childcare waste, waste mixture.
Both facilities together (n = 75)Facility A + B
Table 2. Descriptive statistics for the yeast and mould count in the different datasets used to develop the partial least squares regression analyzed using near infrared spectroscopy.
Table 2. Descriptive statistics for the yeast and mould count in the different datasets used to develop the partial least squares regression analyzed using near infrared spectroscopy.
Samples SetNMeanSDMaxMinCV%
Feed152.91.44.61.448.2
Larvae392.11.14.41.452.3
Frass214.61.66.5234.7
Facility A304.21.66.51.438.1
Facility B 452.21.25.71.454.5
All samples (feed, frass and facilities) 7531.76.51.456.6
N: number of samples, SD: standard deviation, max: maximum, min: minimum, CV: coefficient of variation (SD/mean × 100).
Table 3. Cross validation and prediction statistics for the measurement of yeast and mould count in the different data sets analyzed using near infrared reflectance spectroscopy. N: number of samples, R2cv: coefficient of determination in cross validation, SECV: standard error in cross validation, RPD: residual predictive deviation (SD/SECV), LV: latent variables, SEP: standard error of prediction.
Table 3. Cross validation and prediction statistics for the measurement of yeast and mould count in the different data sets analyzed using near infrared reflectance spectroscopy. N: number of samples, R2cv: coefficient of determination in cross validation, SECV: standard error in cross validation, RPD: residual predictive deviation (SD/SECV), LV: latent variables, SEP: standard error of prediction.
nR2cvSECVBiasSlopeLVSEP
Feed100.980.200.0420.9320.45
Frass180.810.900.130.8020.98
Larvae300.910.27−0.0290.9580.65
Facility A200.900.45−0.0060.9550.67
Facility B300.800.62−0.0070.9360.97
All samples (feed, frass, and facilities)650.740.82−0.0250.8051.03
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

Alagappan, S.; Dong, A.; Mikkelsen, D.; Hoffman, L.C.; Mantilla, S.M.O.; James, P.; Yarger, O.; Cozzolino, D. Near Infrared Spectroscopy for Prediction of Yeast and Mould Counts in Black Soldier Fly Larvae, Feed and Frass: A Proof of Concept. Sensors 2023, 23, 6946. https://doi.org/10.3390/s23156946

AMA Style

Alagappan S, Dong A, Mikkelsen D, Hoffman LC, Mantilla SMO, James P, Yarger O, Cozzolino D. Near Infrared Spectroscopy for Prediction of Yeast and Mould Counts in Black Soldier Fly Larvae, Feed and Frass: A Proof of Concept. Sensors. 2023; 23(15):6946. https://doi.org/10.3390/s23156946

Chicago/Turabian Style

Alagappan, Shanmugam, Anran Dong, Deirdre Mikkelsen, Louwrens C. Hoffman, Sandra Milena Olarte Mantilla, Peter James, Olympia Yarger, and Daniel Cozzolino. 2023. "Near Infrared Spectroscopy for Prediction of Yeast and Mould Counts in Black Soldier Fly Larvae, Feed and Frass: A Proof of Concept" Sensors 23, no. 15: 6946. https://doi.org/10.3390/s23156946

APA Style

Alagappan, S., Dong, A., Mikkelsen, D., Hoffman, L. C., Mantilla, S. M. O., James, P., Yarger, O., & Cozzolino, D. (2023). Near Infrared Spectroscopy for Prediction of Yeast and Mould Counts in Black Soldier Fly Larvae, Feed and Frass: A Proof of Concept. Sensors, 23(15), 6946. https://doi.org/10.3390/s23156946

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