Skip to main content
Download PDF

Effect of replacement of soybean oil by Hermetia illucens fat on performance, digestibility, cecal microbiome, liver transcriptome and liver and plasma lipidomes of broilers

Journal of Animal Science and Biotechnology volume 14, Article number: 20 (2023) Cite this article

Abstract

Background

In contrast to protein-rich insect meal, the feed potential of insect fat is generally less explored and knowledge about the suitability of insect fat as a fat source specifically in broiler diets is still limited. In view of this, the present study aimed to comprehensively investigate the effect of partial (50%) and complete replacement of soybean oil with insect fat from Hermetia illucens (HI) larvae in broiler diets on performance, fat digestibility, cecal microbiome, liver transcriptome and liver and plasma lipidomes. Thus, 100 male, 1-day-old Cobb 500 broilers were randomly assigned to three groups and fed three different diets with either 0 (group HI-0, n = 30), 2.5% (group HI-2.5, n = 35) or 5.0% (HI-5.0, n = 35) Hermetia illucens (HI) larvae fat for 35 d.

Results

Body weight gain, final body weight, feed intake, and feed:gain ratio during the whole period and apparent ileal digestibility coefficient for ether extract were not different between groups. Cecal microbial diversity did not differ between groups and taxonomic analysis revealed differences in the abundance of only four low-abundance bacterial taxa among groups; the abundances of phylum Actinobacteriota, class Coriobacteriia, order Coriobacteriales and family Eggerthellaceae were lower in group HI-5.0 compared to group HI-2.5 (P < 0.05). Concentrations of total and individual short-chain fatty acids in the cecal digesta were not different between the three groups. Liver transcriptomics revealed a total of 55 and 25 transcripts to be differentially expressed between groups HI-5.0 vs. HI-0 and groups HI-2.5 vs. HI-0, respectively (P < 0.05). The concentrations of most lipid classes, with the exception of phosphatidylethanolamine, phosphatidylglycerol and lysophosphatidylcholine in the liver and cholesterylester and ceramide in plasma (P < 0.05), and of the sum of all lipid classes were not different between groups.

Conclusions

Partial and complete replacement of soybean oil with HI larvae fat in broiler diets had no effect on growth performance and only modest, but no adverse effects on the cecal microbiome and the metabolic health of broilers. This suggests that HI larvae fat can be used as an alternative fat source in broiler diets, thereby, making broiler production more sustainable.

Introduction

In recent years, insect biomass has been increasingly recognized as an alternative and sustainable source of feed for monogastric livestock [1,2,3]. Consequently, the use of processed insect biomass as a feed for poultry and pigs has been authorized by the European commission in 2021 [4], with the aim of improving sustainability of food systems and securing food supply. Owing to its high protein content (40%–50% of dry matter (DM) depending on the insect species [5] and provided that the rearing substrate is suitable [6, 7]), research dealing with the feed potential of insect biomass has primarily focused on its role as a source of protein. In this context, a large number of studies in different monogastric livestock species demonstrated that insect larvae meal obtained from different edible insect species suitable for large-scale production, such as Tenebrio molitor (TM) and Hermetia illucens (HI), can replace conventional protein sources, like soybean meal, without impairing performance, metabolic health or product quality [8,9,10,11]. Apart from protein, insect larvae contain a significant amount of fat (up to 47% and 43% of DM in HI and TM, respectively, depending on the rearing substrate [12]), which can be obtained by a defatting process during insect meal production in order to increase the protein content and improve the storage stability of the insect meal. In contrast to insect meal, the feed potential of insect fat is far less explored. Despite that several studies have demonstrated that TM or HI fat has no negative impact on performance, gut morphology, selected blood parameters and product quality in broilers, laying hens and turkeys [13,14,15,16,17,18,19,20,21], knowledge about the suitability of insect fat as a fat source in broiler diets is still limited. In particular, in-depth analysis of the effects of insect fat on the gut microbiome and intermediary metabolism is lacking.

Despite that fatty acid composition of insect larvae is influenced by the fatty acids in the rearing substrate [5], the fatty acid composition of TM and HI fat are completely different. While TM fat consists mainly of unsaturated fatty acids (up to 80%), most of which are long-chain fatty acids (LCFA; C18:1 and C18:2; [11]), the majority of total fatty acids in HI fat are saturated fatty acids (SFA) [5, 16]. This explains that HI fat, but not TM fat has a hard consistency at room temperature. A further characteristic of HI fat is its high proportion of the medium chain-fatty acid (MCFA) lauric acid (C12:0) (≈ 40%; [16]). Interestingly, MCFA including lauric acid have been reported to exert antimicrobial effects (e.g., against Gram-positive cocci and Escherichia coli) in the gut of broilers [22, 23], and, thus, dietary inclusion of HI fat might alter the gut microbiota composition. Considering recent evidence that feed efficiency in broilers is affected by the gut microbiota composition through different mechanisms, such as a more complete digestion of substrates, an increased production of short-chain fatty acids (SCFA) or a decreased stimulation of the intestinal immune system [24, 25], an alteration of the gut microbiota might be of great relevance to broiler´s performance. In addition, convincing evidence exists that the gut microbiota has also a profound influence on energy metabolism and feeding behavior of the host due to the ability of the gut microbiota to communicate with the host along the gut-liver axis via different gut-derived compounds [26]. Moreover, MCFA-rich dietary fats might also exert a direct effect on energy metabolism, because triglyceride (TG)-bound MCFA are considered to be more efficiently absorbed than TG-bound LCFA due to easier emulsification and less dependence on pancreatic lipase [27, 28], thereby increasing digestible energy intake to support a higher growth performance. Against this background, the present study aimed to comprehensively investigate the effect of partial (50%) and complete replacement of soybean oil, the most commonly used fat source in commercial broiler diets, with HI larvae fat in broiler diets on performance, fat digestibility, cecal microbiome, liver transcriptome and liver and plasma lipidomes.

Methods

Animals and diets

The 35-d feeding trial was approved by the Animal Welfare Officer of the Justus Liebig University Giessen (approval no.: JLU 786_M). All experimental procedures described followed established guidelines for the care and handling of laboratory animals. The experiment included 100 male, 1-day-old broiler chickens (Cobb 500, Cobb, Wiedemar, Germany), which were randomly assigned to three groups (5 broilers/cage, group 1: 6 cages, group 2 and group 3: 7 cages). The mean initial body weight (BW) (44.9 ± 2.4 g; mean ± SD) was similar across the groups. The broilers were kept in 2.1 m2 cages equipped with nipple drinkers and feed automates and had free access to feed and water. At the floor of the cages there were cardboards, which were covered with litter to allow scratching, pecking and dustbathing. Cardboards and litter were exchanged 2 times per week during the first two weeks and every 2 d during the last three weeks of the trial. In addition, broilers were provided with perches in elevated position for resting and sleeping. Light intensity was constantly at 40 Lux and the light regime was 24 h:0 h, 23 h:1 h, 22 h:2 h, 21 h:3 h, 20 h:4 h, 19 h:5 h (light:dark) at d 1, 2, 3, 4, 5, 6, and 18 h:6 h from d 7 onward, as recommended by the breeder [29]. The room temperature decreased from 28–29 °C on d 1, measured at pen height, to 23–24 °C on d 35. During the first 6 d, infrared lamps (Albert Kerbl GmbH, Buchbach, Germany) were used as additional heat sources in order to adjust the temperature at the cage floor to 34 °C. Mean relative humidity was 60.0% ± 1.9%. The groups (group 1: HI-0, group 2: HI-2.5, group 3: HI-5.0) were fed three different diets, which varied only in the fat source (HI-0: 0 HI fat and 5% soybean oil, HI-2.5: 2.5% HI fat and 2.5% soybean oil, HI-5.0: 5.0% HI fat and 0 soybean oil), in a three-phase feeding system (starter diet form d 1 to 10, grower diet from d 11 to 21, finisher diet from d 22 to 35). The HI fat was purchased from Madebymade (Pegau, Germany) and stored at −20 °C until diet preparation. Prior to diet preparation, the stability of HI fat was assessed by determining the acid value and the percentage of polar compounds [30]. Low values for both parameters indicated no significant oxidation of the HI fat. Company’s details on the rearing conditions (substrate, duration) and further processing of the HI larvae are not available for reasons of confidentiality. The composition of the three diets are shown in Table 1. The finisher diets contained 0.5% titanium dioxide as indicator in order to calculate the apparent ileal digestibility (AID) coefficient for ether extract by the indicator method [31]. The diets met the broiler’s requirements of nutrients and energy according to the breeder’s recommendations [32]. For diet preparation, three different basal feed mixtures for the starter, grower and finisher period, respectively, were produced from all feed components except the fat source using mixing machines (V100 and V250, Diosna, Germany). The experimental diets were obtained by mixing the different basal feed mixtures with the intended amounts of fat sources (soybean oil, HI fat) using the mixing machine. Prior to adding the HI fat, the HI fat was melted in an electric oven at 50 °C. Subsequently, the experimental diets were pelleted using a pelleting device (V3/30 C, Simon-Heesen, Boxtel, Netherlands) and aliquots of all diets were stored at −20 °C for analysis of diet composition. Diets were fed in crumbled form during the first 3 d, and in pellet form (2 mm diameter) from d 3 until the end of the trial. Body weight (individually) and feed intake (per cage) were determined on d 1, 10, 21 and 35, and the feed:gain ratio was calculated from feed intake and body weight gains on cage basis.

Table 1 Composition of the broiler diets, g/kg
Full size table

Analysis of diet composition

Concentrations of DM, crude protein, crude ash, ether extract, crude fiber and amino acids in the main diet components (wheat, maize, soybean extraction meal) and the experimental diets were determined by official methods [33]. Concentrations of sugar and starch were analyzed in the experimental diets by official methods (sugar according to method 7.2.1 and starch according to method 7.1.3 [33]). Total lipid fatty acid composition of HI fat, soybean oil and diets were analyzed as described below. Apart from total lipid fatty acid composition, the detailed composition (all major lipid classes and individual lipid species) of the HI fat was comprehensively analyzed using lipidomics, as described for liver and plasma (see below). The apparent N-corrected metabolizable energy (AMEN) content of the diets was calculated according to the formula of the World’s Poultry Science Association (WPSA) for poultry compound feed [34]:

$${\mathrm{AME}}_{\mathrm N}\;(\mathrm{MJ}/\mathrm{kg})\:=\left[\left(0.01551\:\times\:\mathrm{crude}\;\mathrm{protein}\right)\:+\left(0.03431\:\times\:\mathrm{ether}\;\mathrm{extract}\right)\:+\;\left(0.01669\:\times\:\mathrm{starch}\right)\:+\:\left(0.01301\:\times\:\mathrm{sugar}\right)\right]$$

Considering that the AME and AMEN content of HI larvae fat for broiler chickens was shown to be similar to that of soybean oil [35], the formula is appropriate to calculate the ME contents of the experimental diets.

Sample collection

A total of 12 (group HI-0) and 14 (groups HI-2.5 and HI-5.0) broilers per group (two broilers from each cage), whose body weights represented the mean body weight of the whole group, were selected for sample collection in order to avoid that effects were biased by random selection of broilers with very low or very high body weights. All analyses described below were carried out in these animals. The animals were killed by bleeding (opening of Vena jugularis and Arteria carotis) under electrical anesthesia using a BTG-40A stunning device (Westerhoff Geflügeltechnik, Hoogstede, Germany) in accordance with the European legislation for euthanasia of animals [36]. Whole blood was collected into ethylenediaminetetraacetic acid-coated polyethylene tubes (9 mL S-Monovette, Sarstedt, Nümbrecht, Germany). Plasma was prepared by centrifugation (1100 × g, 10 min) at 4 °C and stored at −20 °C. The liver was excised, washed in ice-cold NaCl solution (0.9%), weighted and small aliquots were snap-frozen in liquid nitrogen and stored at −80 °C. The gastrointestinal tract was removed and digesta from the ileum (segment between Meckel´s diverticulum and the ileo-cecal junction) and the cecum was collected. Tissue and digesta samples were snap-frozen in liquid nitrogen and stored at −80 °C pending analysis.

Determination of AID coefficient for ether extract

The AID coefficient for ether extract was determined at the end of the experiment by the indicator method using titanium dioxide as an inert marker [31]. Prior to analysis, ileal digesta samples were freeze-dried and ground using a centrifugal mill (Retsch, Haan, Germany). Ileal digesta concentration of the indigestible indicator TiO2 was determined by the method of Brandt and Allam with slight modifications [37]. Concentration of ether extract in the ileal digesta was determined by official methods as described above. Based on the ileal concentrations of indicator, the AID coefficient for ether extract was calculated according to the following formula:

$$\mathrm{AID}\;\mathrm{coefficient}\;(\%)\:=\:100-\lbrack({\mathrm{TiO}}_{2_\_\mathrm{Diet}}/{{\mathrm{TiO}}_2}_{{}_\_\mathrm{Di}\mathrm g\mathrm e\mathrm s\mathrm t\mathrm a})\:\times\:({{\mathrm{EE}}_\_}_{\mathrm{Digesta}}/{\mathrm{EE}}_{{}_\_\mathrm{Di}\mathrm e\mathrm t})\:\times\:100\rbrack,$$

in which TiO2_Diet is the TiO2 concentration in the diet (% DM), TiO2_Digesta is the TiO2 concentration in the ileal digesta (% DM), EE_Digesta is the ether extract concentration in ileal digesta (% DM), and EE_Diet is the ether extract concentration in the diet (% DM).

Determination of microbiota composition and diversity in the cecal digesta

Metagenomic DNA was isolated from approximately 180–200 mg of cecal digesta using genomic DNA columns (Macherey‐Nagel, Düren, Germany) according to Lagkouvardos et al. [38]. V3-V4 regions of the 16S rRNA genes were amplified using bacteria‐specific primers following a two-step procedure according to the Illumina sequencing protocol as described [38]. Amplicons were sequenced using a MiSeq system (Illumina, Inc., San Diego, CA, USA). Further processing of raw sequences was carried out as described recently [39]. Finally, sequences with a relative abundance > 0.1% in at least one sample were sorted, merged and operational taxonomic units (OTU) were picked at a threshold of 97% similarity. Taxonomic classification to the OTU was assigned using the SILVA database [40]. Further downstream analyses were done using Rhea (https://lagkouvardos.github.io/Rhea/). The differential abundance analysis of taxa was performed on the aggregated data at the different taxonomic levels as described [38]. For estimation of diversity within samples (α-diversity), the Shannon and Simpson indices, the most common indices to compare diversity, were calculated and transformed to the corresponding effective number of species according to Jost [41], because they are better suited at indicating the true diversity between samples and are minimally affected by the number of rare species. To measure the similarity between different microbial profiles, the β-diversity was determined by calculating generalized UniFrac distances with PERMANOVA statistical test as described previously [38]. Visualization of bacterial profiles among different groups was done by computation of non-metric multidimension distance scaling (NMDS) [42].

Determination of SCFA concentrations in the cecal digesta

Cecal digesta SCFA concentrations were determined as described previously [43]. In brief, 50 mg aliquots of cecal digesta were mixed with 0.5 mL 5% o-phosphoric acid containing internal standard (0.15 mg/mL crotonic acid). Extraction was carried out by vortexing for 3 min and subsequent centrifugation at 21,100 × g at 4 °C for 10 min. Prior to injection, the supernatant was centrifuged again at 21,100 × g at 4 °C for 5 min. 1 μL of the extract was injected into a gas chromatograph (Clarus 580 GC system, Perkin Elmer, Waltham, USA) equipped with a polar capillary column (10 m free fatty acid phase, 0.32 mm internal diameter, 0.25 μm film thickness; Macherey and Nagel, Düren, Germany) and a flame ionisation detector.

RNA extraction and hepatic transcript profiling

Total RNA from liver aliquots (20 mg) were isolated using TRIzol reagent (Invitrogen, Karlsruhe, Germany) according to the manufacturer’s protocol. RNA quantity and quality were assessed spectrophotometrically using an Infinite 200 M microplate reader equipped with a NanoQuant plate (both from Tecan, Mainz, Germany). The average RNA concentration and the A260/A280 ratio of all total RNA samples (n = 40, means ± SD) were 421 ± 48 ng/μL and 1.90 ± 0.02. For hepatic transcript profiling, total RNA samples from six randomly selected broilers/group were sent on dry-ice to the Genomics Core Facility “KFB—Center of Excellence for Fluorescent Bioanalytics” (Regensburg, Germany). Following a further RNA quality check using an Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany), which revealed an average RNA integrity number (RIN) value of 8.23 ± 0.19 for all samples (n = 18, means ± SD), total RNA samples were processed using an Affymetrix GeneChip Array (Chicken Gene 1.0 ST), which covers 18,214 genes represented by 439,582 probes, according to the Applied Biosystems™ GeneChip™ Whole Transcript (WT) PLUS Reagent Kit User Guide (Thermo Fisher Scientific, Waltham, MA, USA). Following scanning of the processed GeneChips, cell intensity files, which provided a single intensity value for each probe cell, were generated from the image data using the Command Console software (Affymetrix). The compressed array image files (CEL files) were imported into the Applied Biosystems™ Transcriptome Analysis Console (TAC) (v. 4.0.2) software (Thermo Fisher Scientific) for calculation of summarized probe set signals (in log2 scale) using the Robust Multichip Analysis algorithm, comparison fold changes (FC) and significance P-values (ANOVA). Annotation of the microarrays was performed with the “ChiGene-1_0-st-v1.na36.galgal3.transcript.csv” annotation file. The microarray data of this study have been deposited in MIAME compliant format in the NCBI´s Gene Expression Omnibus public repository [44]. Owing to the rather moderate differences in the hepatic transcriptomes between groups HI-5.0 vs. HI-0 and groups HI-2.5 vs. HI-0, the differentially expressed transcripts were filtered based on a fold change > 1.5 or < −1.5 and a P-value < 0.05. Identical or similar filter criteria were also applied in several recent studies [45, 46]. Filtering of differentially expressed transcripts using the Benjamini & Hochberg false discovery rate adjustment method could not be applied, because the adjusted P-values for all transcripts were > 0.05.

Gene set enrichment analysis (GSEA) was performed with the identified differentially expressed transcripts in order to identify enriched Gene Ontology (GO) terms within GO category biological process using the Database for Annotation, Visualization and Integrated Discovery (DAVID) 6.8 bioinformatic resource [47, 48]. Biological process and molecular function terms were considered as enriched if P < 0.05.

Validation of microarray data using qPCR analysis

Microarray data of 16 differentially expressed transcripts were validated by qPCR. For qPCR analysis, total RNA from all broilers (n = 12–14/group) was used. Synthesis of cDNA and qPCR analysis was performed with a Rotor-Gene Q system (Qiagen, Hilden, Germany) as described recently in detail [49]. Gene-specific primers were synthesized by Eurofins MWG Operon (Ebersberg, Germany). Characteristics of primers are listed in Additional file 1: Table S1. Normalization was carried out using multiple reference genes as described recently [50].

Tissue homogenization and lipid extraction

Frozen liver tissue was homogenized in methanol/water (50/50, v/v) with addition of 1% sodium laurylsulfate using bead-based homogenization at a concentration of 0.05 mg/µL [51]. Lipid class specific, non-endogenous internal standards were added prior to lipid extraction. An amount of 2 mg liver (wet weight) or a volume of 10 µL plasma was subjected to lipid extraction according to the protocol by Bligh and Dyer [52]. A volume of 0.5 mL of the chloroform phase was recovered by a pipetting robot and vacuum dried. The residue was dissolved in 1.2 mL chloroform/methanol/2-propanol (1:2:4, v/v/v) with 7.5 mmol/L ammonium formate.

Lipidomic analysis of major lipid classes by mass spectrometry

The analysis of lipids was performed by direct flow injection analysis (FIA) using a high-resolution Fourier Transform (FT) hybrid quadrupole-Orbitrap mass spectrometer (FIA-FTMS) [53]. TG, diglycerides (DG) and cholesteryl esters (CE) were recorded in positive ion mode as [M + NH4]+ in m/z range 500–1000 and a target resolution of 140,000 (at m/z 200). CE species were corrected for their species-specific response [54]. Ceramides (Cer), phosphatidylcholines (PC), ether PC (PC O), phosphatidylethanolamines (PE), ether PE (PE O), phosphatidylglycerols (PG), phosphatidylinositols (PI), and sphingomyelins (SM) were analyzed in negative ion mode in m/z range 520–960; lysophosphatidylcholines (LPC) and lysophosphatidylethanolamine (LPE) in m/z range 400–650. Multiplexed acquisition (MSX) was applied for free cholesterol (FC) and the internal standard FC[D7] [54]. Lipid annotation is based on the latest update of the shorthand notation [55].

The datasets from liver and plasma lipidomes were subjected to principal component analysis (PCA) using the MetaboAnalystR 3.2 package for R version 4.2.1. For the PCA, the relative metabolite composition of individual lipid species within the different lipid classes were used. Prior to the PCA, variables with missing values were either excluded from the analyzes if more than 50% of the samples were missing or the missing values were replaced by the limit of detection (1/5 of the minimum positive value of each variable). After normalization by log transformation and autoscaling the remaining values were used for the PCA.

Determination of fatty acid composition of total lipids of the diets and the liver

Fatty acid composition of total lipids of the liver and the diets was determined by gas chromatography-flame ionization detection (GC-FID). Briefly, total lipids were extracted from 75 mg liver aliquots with a 3:2 (v/v)-mixture of n-hexane and isopropanol containing C19:0 (50 mg/mL) as internal standard. After extraction, samples were centrifuged (1200 × g, 10 min) and an aliquot of the supernatant was evaporated under a stream of N2 at 37 °C. Lipids were subsequently transmethylated using trimethylsulfonium hydroxide solution (Sigma-Aldrich) and the resulting fatty acid methyl esters (FAME) were separated by a GC-FID system described in detail recently [56].

Statistical analysis

Statistical analysis was conducted with SPSS 27 software (IBM, Armonk, NY, USA). The cage served as the experimental unit for feed intake and feed:gain ratio and the individual animal for all other data. All parameters were tested with the Shapiro–Wilk test for normal distribution and with the Levene’s test for homoscedasticity. When the normal distribution was followed only after a log transformation, the log transformed data were used for statistical analysis. Differences between the three groups were analyzed by one-way analysis of variance (one-way ANOVA) followed by a Tukey’s post-hoc test when the data were normally distributed and the variances were homogeneous. If the data showed heterogeneity of variance, the means of the three groups were analyzed using Welch's ANOVA in conjunction with the post-hoc Games-Howell test. If the normal distribution was not followed, a Kruskal–Wallis one-way ANOVA was performed using the Mann–Whitney U test with Bonferroni correction as post-hoc test. For all tests, a P-value < 0.05 was considered statistically significant.

Results

Lipid composition of the HI fat

Analysis of the composition of major lipid classes of the HI fat revealed that HI total lipids consisted almost completely of TG (99.2% of total fat). All other lipid classes detected (PG, PC, DG, PE and SM) made up < 0.3% of total lipids (Table 2). Analysis of individual TG species demonstrated that TG 36:0, TG 38:0, TG 40:0 and TG 42:2 were the most abundant TG species, whereas all other TG species contributed < 5% of all TG species (Table 2). The majority of TG species of the HI fat was saturated (56.7%), while TG species with one, two and three or more double bonds made up 14.0%, 17.3% and 12.0%, respectively. Due to the low proportions of PG, PC, DG, PE, PC-O, SM and LPC in the HI fat, the individual lipid composition of these non-TG lipid classes is not reported. Analysis of the fatty acid composition of HI fat by GC-FID revealed that SFA were the dominating fatty acids, with C12:0 (57.2%), C14:0 (8.7%) and C16:0 (10.9%) contributing to 76.8% of total fatty acids (data not shown). The essential fatty acids C18:2 n-6 and C18:3 n-3 made up 10.6% and 0.8%, of total fatty acids, respectively.

Table 2 Lipid composition of the Hermetia illucens larvae fat
Full size table

Composition of the experimental diets

The three experimental diets within each feeding phase had similar concentrations of crude nutrients, sugar, starch and energy (Table 3) and amino acids (Additional file 1: Table S2), but substantially differed in the fatty acid composition of dietary total lipids (Table 3). With increasing replacement of soybean oil by HI fat, the proportions of C10:0, C12:0 and C14:0 markedly increased, while those of C18:0, C18:1, C18:2 n-6 and C18:3 n-3 decreased. As a consequence, the dominating fatty acids in the HI-5.0 diet were in decreasing order: C12:0, C18:2 n-6, C18:1, C16:0 and C14:0; in the HI-2.5 diet the dominating fatty acids were in decreasing order: C18:2 n-6, C18:1, C12:0, C16:0 and C18:3 n-3; in the HI-0 diet the dominating fatty acids were in decreasing order: C18:2 n-6, C18:1, C16:0, C18:3 n-3 and C18:0. Monounsaturated fatty acids and polyunsaturated fatty acids in the HI fat total lipids made up 9.4% and 11.4%, respectively.

Table 3 Concentrations of nutrients and energy in the broiler diets
Full size table

Performance and ether extract digestibility

Performance data (body weight gain, final body weight, feed intake, and feed:gain ratio) during the whole period were not different between groups. However, body weight gain and feed intake during the starter period were higher in groups HI-2.5 and HI-5.0 than in group HI-0 (P < 0.05, Table 4). No differences of performance data between groups were found during the grower and the finisher period. AID coefficient for ether extract also did not differ between groups.

Table 4 Performance data and apparent ileal digestibility (AID) coefficient for ether extract of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d
Full size table

Cecal microbiota diversity and composition

In order to identify alterations of the cecal microbiota structure of the broilers, 16S rRNA-based high-throughput sequencing was applied. Following quality check, chimera check and filtering, the high-quality sequences obtained from the cecum digesta samples of the 40 broilers were delineated into 90 OTUs at 97% sequence identity (Additional file 1: Table S3). Treatment effect on microbial diversity was evaluated by the use of different diversity metrics. None of the metrics used to describe α-diversity (effective species richness, shannon effective, simpson effective, evenness) differed between groups (Fig. 1a). In addition, β-diversity of cecal bacterial community calculated based on generalized UniFrac distances did not differ among groups. The MetaNMDS plot generated to visualize the difference in β-diversity of cecal bacterial community among groups shows that no clustering was most visible among groups (Fig. 1b).

Fig. 1
figure 1

Analysis of the cecal microbiome. Indicators of α-diversity (effective richness, shannon effective, simpson effective, evenness) of the cecal bacterial community (A), visualization of the difference in β-diversity of cecal bacterial community between groups by non-metric multidimension distance scaling plot (B), and distribution of cecal bacteria at different taxonomic levels (phylum, order, family) (C) of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. A: Box and whisker plots for n = 12–14 broilers/group; C: Data are means for n = 12–14 broilers/group. a,bMeans without a common letter differ across the groups, P < 0.05

Full size image

To analyze the effect on the microbiota composition, the microbial community was analyzed at different taxonomic levels (phylum, class, order, family, genus). In total, only the relative abundances of four different bacterial taxa differed among groups. At the phylum level, only the relative abundance of the least abundant bacterial phylum, Actinobacteriota, contributing less than 0.15% to total bacterial phyla was different among groups. While the relative abundance of this phylum was lower in group HI-5.0 than in group HI-2.5 (P < 0.05), it did not differ between group HI-5.0 and HI-0 and between group HI-2.5 and HI-0 (Fig. 1c). Amongst five different bacterial classes identified in the cecal digesta of the broilers, only the relative abundance of Coriobacteriia, which belong to the phylum Actinobacteriota, was affected by treatment; while the relative abundance of this phylum was lower in group HI-5.0 than in group HI-2.5 (P < 0.05), it did not differ between group HI-5.0 and HI-0 and between group HI-2.5 and HI-0. At the order level, the relative abundance of Coriobacteriales (class Coriobacteriia) was identified as the only bacterial order differing between groups; its relative abundance was lower in group HI-5.0 than in group HI-2.5 (P < 0.05), but was not different between groups HI-5.0 and HI-2.5 (Fig. 1c). At the family level, the only taxon identified to be altered was the Eggerthellaceae (order Coriobacteriales), whose abundance was lower in group HI-5.0 than in group HI-2.5 (P < 0.05), but was not different between group HI-5.0 and HI-0 and between group HI-2.5 and HI-0 (Fig. 1c). No differences between groups were found regarding the relative abundances of bacterial genera.

Cecal concentrations of short-chain fatty acids

The concentrations of total and individual SCFA [acetic acid (C2:0), propionic acid (C3:0), isobutyric acid (iC4:0), butyric acid (C4:0), isovaleric acid (iC5:0) and valeric acid (C5:0)] in the cecal digesta were not different between the three groups (Fig. 2).

Fig. 2
figure 2

Analysis of microbial fermentation profile in the gut. Concentrations of individual [acetic acid (C2:0), propionic acid (C3:0), butyric acid (C4:0), isobutyric acid (iC4:0), valeric acid (C5:0), isovaleric acid (iC5:0)] and total (= sum of individual) short-chain fatty acids (SCFA) in the cecal digesta of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Data are means ± SD, n = 12–14 broilers/group

Full size image

Liver transcriptome

According to the filter criteria applied (P < 0.05; fold change > 1.5 and < −1.5), a total of 55 transcripts were identified as differentially expressed (upregulated: 31, downregulated: 24) in the liver between groups HI-5.0 and HI-0. In Fig. 3a, the differentially expressed transcripts between groups HI-5.0 and HI-0 are illustrated as red dots in the Volcano plot. Amongst the upregulated genes, only three genes (BDKRB1, XDH, IGJ) exhibited a regulation > 2.0-fold. The top 10 upregulated transcripts were in decreasing order of their fold change (in brackets): BDKRB1 (2.57), XDH (2.36), IGJ (2.10), DPP4 (1.90), ENDOUL (1.83), UPP2 (1.82), EMB (1.81), BASP1 (1.76), SIK1 (1.71) and B3GALT2. Amongst the downregulated genes, five genes (IL22RA2, ANKRD22, LOC101749538, ATP2B2, SLC6A13) were regulated < −2.0-fold. The top 10 downregulated transcripts were in increasing order of their fold change: IL22RA2 (−5.09), ANKRD22 (−3.76), LOC101749538 (−2.55), ATP2B2 (−2.15), SLC6A13 (−2.12), DDO (−1.89), SCAP (−1.88), SPATA4 (−1.84), CA4 (−1.82) and ULK1 (−1.81). The fold change and P-value of all differentially expressed transcripts between group HI-5.0 vs. HI-0 are shown in Additional file 1: Table S4.

Fig. 3
figure 3

Differential transcriptome analysis of the liver. Volcano plot illustrating the differentially expressed transcripts in the liver of broilers between group HI-5.0 vs. HI-0 (A) and group HI-2.5 vs. HI-0 (B). The double filtering criteria are indicated by horizontal (P-value < 0.05) and vertical (fold change > 1.5 or < −1.5) dashed lines. Red dots in the upper left and the upper right corner represent the downregulated and the upregulated transcripts, respectively. C Most enriched gene ontology (GO) biological process terms associated with the transcripts differentially expressed between group HI-5.0 vs. HI-0. GO terms are sorted by their enrichment P-values (EASE score) (top: lowest P-value, bottom: highest P-value)

Full size image

Considering the same filter criteria as for the comparison of groups HI-5.0 vs. HI-0, a total of 25 transcripts were identified as differentially expressed (upregulated: 19, downregulated: 6) in the liver between groups HI-2.5 and HI-0 (Fig. 3b). Only three genes (NSMF, CHIA, BDKRB1) were regulated > 2.0-fold and no genes were regulated < −2.0-fold. The 10 most strongly upregulated transcripts were in decreasing order of their fold change (in brackets): NSMF (2.29), CHIA (2.16), BDKRB1 (2.04), ENPEP (1.80), PFKFB3 (1.80), MIR1790 (1.72), SIK1 (1.71), DPP4 (1.71), LOC395159 (1.69) and CP (1.66). The 6 downregulated transcripts were in increasing order of their fold change: GUCY2C (−1.53), ATP2B2 (−1.53), INDOL1 (−1.51), MIR17 (−1.51), DDO (−1.51) and PROM1 (−1.50). The fold change and P-value of all differentially expressed transcripts between groups HI-2.5 vs. HI-0 are shown in Additional file 1: Table S5.

Microarray data of 16 differentially expressed transcripts between groups HI-5.0 and HI-0 were validated by qPCR. As shown in Additional file 1: Table S6, the effect direction (positive or negative fold change) was the same between microarray and qPCR for all validated transcripts, whereas the effect size (value of fold change) differed to some extent for the validated transcripts between microarray and qPCR. Statistical analysis of qPCR data revealed that twelve transcripts were regulated either significantly (CA4, DDO, IL22RA2, HMGCS1, SCAP, SH2D4A, ULK1, XDH) or at a P-value < 0.1 (AKR1D1, ANKRD22, MTHFS, OGN), whereas four transcripts (DPP4, FANCL, SIK1, SPATA4) were not regulated.

In order to extract biological meaning from the transcripts differentially expressed between groups HI-5.0 and HI-0, GSEA was performed using GO category biological process. Due to the low number of differentially expressed transcripts, GSEA was not conducted for the comparison of groups HI-2.5 and HI-0. Within GO category biological process, the most enriched biological process terms assigned to the transcripts regulated between groups HI-5.0 and HI-0 were (in increasing order of their P-values): small molecule biosynthetic process, cholesterol metabolic process, secondary alcohol metabolic process, sterol metabolic process, lipid biosynthetic process, steroid biosynthetic process, cholesterol biosynthetic process and secondary alcohol biosynthetic process (Fig. 3c).

Liver and plasma lipidomes

The major lipid classes in the liver were in decreasing order: TG, PC, PE, FC, PI, DG, SM, CE, PE O, Cer, LPC, LPE, PC O and PG. Despite that the sum of total lipid classes in the liver were not different between groups (HI-0: 47.6 ± 8.9 µmol/g, HI-2.5: 42.4 ± 8.9 µmol/g, HI-5.0: 47.7 ± 6.5 µmol/g, mean + SD, n = 12–14/group, P = 0.170), hepatic concentrations of PE, LPC, LPE and PG differed between groups (P < 0.05), even though only slightly (Fig. 4a). Hepatic concentrations of LPE and PG were higher in group HI-5.0 than in groups HI-0 and HI-2.5 (P < 0.05). Concentrations of PE and LPC in the liver were higher in group HI-5.0 than in group HI-2.5 (P < 0.05), but not compared to group HI-0.

Fig. 4
figure 4

Analysis of the liver lipidome. Concentrations of the major lipid classes in the liver (A) and lipid species composition within the three most abundant lipid classes in the liver, TG (B), PC (C) and PE (D), of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Lipid species were grouped according to the number of double bonds (DB). Data are means ± SD, n = 12–14 broilers/group. Abbreviations: CE, cholesteryl ester; Cer, ceramide; DG, diglycerides; FC, free cholesterol; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; PC, phosphatidylcholine; PC O, PC-ether; PE, phosphatidylethanolamine; PE O, PE-ether; PG, phosphatidylglycerol; PI, phosphatidylinositol; SM, sphingomyelin; TG, triglycerides

Full size image

In contrast to the slight or absent differences in the concentrations of the lipid classes in the liver between groups, pronounced differences between groups were seen with regard to the composition of individual lipid species within the different lipid classes. Regarding TG, the dominating lipid class in the liver, relative proportions of species with zero and one double bonds were higher in group HI-5.0 than in group HI-0, while those with five and six double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 4b). With regard to PC, the second most abundant hepatic lipid class, relative proportions of species with one and three double bonds were higher and those with four, five and six double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 4c). Within PE, the third most abundant lipid class in the liver, relative proportions of species with two and three double bonds were higher and those with four, six and eight double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 4d). The proportions of all individual TG, PC and PE species are shown in Additional file 1: Table S7.

In plasma, the major lipid classes were in decreasing order: PC, CE, FC, TG, PE, SM, LPC, PC O, PE O, DG, LPE and Cer (Fig. 5a). Like the concentrations of most lipid classes (PC, FC, TG, PE, SM, LPC, PC O, PE O, DG, LPE), the sum of total lipid classes in plasma did not differ among groups (HI-0: 7992 ± 1423 µmol/L, HI-2.5: 7683 ± 1742 µmol/L, HI-5.0: 8550 ± 752 µmol/L, mean + SD, n = 12–14/group, P = 0.358). In contrast, the plasma concentrations of CE and Cer differed among groups (P < 0.05), but only slightly; plasma concentrations of CE and Cer were higher in group HI-5.0 than in group HI-0. Plasma concentration of Cer in group HI-5.0 was also higher than in group HI-2.5 (P < 0.05), whereas plasma concentration of CE did not differ between group HI-5.0 and group HI-2.5. Similar as in the liver, pronounced differences among groups were observed regarding the composition of individual lipid species within the different lipid classes in plasma. With regard to CE, the dominating lipid class in the plasma, relative proportions of species with zero and one double bonds were higher in group HI-5.0 than in group HI-0, while those of species with two, three, four, five and six double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 5b). Regarding PC, the second most abundant plasma lipid class, the relative proportions of PC species with zero, one and three double bonds were higher, whereas those with four, five, six and seven double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 5c). With regard to TG, relative proportions of species with zero, one, and two double bonds were higher and those with four, five, six, seven and eight double bonds were lower in group HI-5.0 than in group HI-0 (P < 0.05, Fig. 5d). The proportions of all individual CE, PC and TG species in plasma are shown in Additional file 1: Table S8.

Fig. 5
figure 5

Analysis of the plasma lipidome. Concentrations of the major lipid classes in the plasma (A) and lipid species composition within the three most abundant lipid classes in plasma, CE (B), PC (C) and TG (D), of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Lipid species were grouped according to the number of double bonds (DB). Data are means ± SD, n = 12–14 broilers/group. Abbreviations: CE, cholesteryl ester; Cer, ceramide; DG, diglycerides; FC, free cholesterol; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; PC, phosphatidylcholine; PC O, PC-ether; PE, phosphatidylethanolamine; PE O, PE-ether; SM, sphingomyelin; TG, triglycerides

Full size image

The pronounced differences among groups with regard to the composition of individual lipid species within the different lipid classes in liver and plasma were reflected by the results from PCA analysis which was carried out with the relative proportions of individual lipid species of all lipid classes. The dimensional reduction of the lipidome datasets from the liver (Fig. 6a) and the plasma (Fig. 7a) revealed a clear separation (rightward shift) between the three feeding groups with increasing dietary levels of HI fat. For the liver lipidome the cumulative proportion is 54.6% with principal component 1 accounting for more than 40% of the variance of the dataset. In the plasma, the cumulative proportion adds up to 51.8% with principal component 1 accounting for 38.8%. The loading plot for the liver lipidome dataset shows that the left shift of the HI-0 group was mainly caused by individual TG and PC lipid species with three or more double bonds and the right shift of the HI-2.5 and HI-5.0 groups was mainly caused by individual TG, PC and PE lipid species with less than 3 double bounds (Fig. 6b). In the plasma, the right shift of the HI-2.5 and HI-5.0 groups was largely driven by individual TG, PC and PE species with less than three double bounds, while the left shift of the HI-0 group is driven by individual TG, PC and CE species with more than 3 double bounds (Fig. 7b).

Fig. 6
figure 6

Principal component analysis of liver lipidome. Scores plot with plotted 5% confidence interval (A) and associated loading plot (B) of principal component analysis (PCA) of the liver lipidome of the broilers. Broilers were fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Data are principal components (PC 1 or PC 2) and their loadings calculated based on the relative abundances of the lipid species within their lipid classes, n = 12–14 broilers/group. Abbreviations: Cer, ceramide; DG, diglycerides; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; PC, phosphatidylcholine; PC O, PC-ether; PE, phosphatidylethanolamine; PE O, PE-ether; PG, phosphatidylglycerol; PI, phosphatidylinositol; TG, triglycerides

Full size image
Fig. 7
figure 7

Principal component analysis of plasma lipidome. Scores plot with plotted 5% confidence interval (A) and associated loading plot (B) of principal component analysis (PCA) of the plasma lipidome of the broilers. Broilers were fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Data are principal components (PC 1 or PC 2) and their loadings calculated based on the relative abundances of the lipid species within their lipid classes, n = 12–14 broilers/group. Abbreviations: CE, cholesteryl ester; Cer, ceramide; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; PC, phosphatidylcholine; PC O, PC-ether; PE, phosphatidylethanolamine; PE O, PE-ether; SM, sphingomyelin; TG, triglycerides

Full size image

Fatty acid composition of liver total lipids

As expected, fatty acid composition of hepatic total lipids was affected by dietary inclusion of HI fat. Proportions of C12:0, C14:0 and C16:1 were higher in group HI-5.0 than in groups HI-2.5 and HI-0, whereas proportion of C18:3 n-3 was lower in groups HI-5.0 and HI-2.5 than in group HI-0 (P < 0.05, Table 5).

Table 5 Fatty acid composition of hepatic total lipids of broilers fed diets with either 0 (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d
Full size table

Discussion

In the present study, the effect of replacement of soybean oil with HI larvae fat in broiler diets on the cecal microbiome and metabolic health was comprehensively investigated using omics technologies. Since no study existed in the literature reporting the lipid composition of HI fat in detail, with the exception of total lipid fatty acid composition, lipidomics was applied for the first time to fully characterize the HI fat. Our lipidomic analysis showed that the HI larvae fat used consisted almost completely (> 99%) of TG, whereas PG, PC, DG, PE and SM were only minor lipid classes. Amongst the TG species, saturated ones (34:0, 36:0, 38:0, 40:0, 42:0) made up approximately 60% of all TG species. TG species with one or two double bonds contributed to nearly 30% of all TG species, while those with three or more double bonds contributed less than 10% of all TG species. In line with the high proportion of saturated TG species and TG species with a low number of double bonds, additional fatty acid analysis of HI larvae fat total lipids revealed that C12:0, C16:0, C14:0 and C10:0 in decreasing order of their proportions made up nearly 80% of total fatty acids. Our observation that saturated TG species made up 60% of all TG species, while saturated fatty acids made up 80% of all fatty acids in HI fat total lipids, however, is not a discrepancy. This is likely attributed to the fact that saturated fatty acids are not exclusively found in TG species with zero double bonds, but are also present in TG species with 1 or more double bonds. Regardless of this, our data from lipidomics and fatty acid analysis clearly show that the HI fat used in the present study contained large amounts of saturated fatty acids, which largely explains the hard consistency of the HI fat at room temperature. Despite the marked difference in the fatty acid profile between soybean oil, the commonly used fat source in broiler diets, and HI larvae fat, partial or complete replacement of soybean oil with HI larvae fat had no effect on performance parameters (feed intake, body weight gain, feed:gain ratio) of the broilers over the whole 35 d-period. In line with the unaltered growth performance, ileal digestibility of ether extract determined during the finishing period did not differ among groups. This indicates that soybean oil replacement with HI larvae fat has neither beneficial nor adverse effects on fat digestibility, feed utilization and growth performance of broilers. However, our study revealed that feed intake and body weight gain, but not feed efficiency were increased by dietary inclusion of HI larvae fat during the starter period, whereas no effect was found during the grower and the finisher period. This suggests that feed acceptance was obviously improved in broilers during the early d of life by the use of HI larvae fat. The reason underlying this observation is unclear, but we exclude the possibility that the taste of the soybean oil used in the HI-0 diet was impaired due to fat deterioration during storage or diet preparation, because the same batch of soybean oil was also used for the diet of group HI-2.5, despite feed intake of this group was similar as in group HI-5.0. Unfortunately, comparable data from other studies are not available to explain our observation. While in two other studies dealing with the use of HI larvae fat as an alternative fat source in broiler diets performance data were reported either for the whole period [16] or only the finishing period [13], a further study did not provide any performance data [17]. However, the latter study is not suitable due to inadequate control of nutrient composition of the diets as evident from the marked difference of crude fat content between the starter diets (3.6% and 5.9% in the control diet and HI larva fat diet, respectively).

Despite studies demonstrating antimicrobial effects of MCFA in the intestine of broilers at dosages between 0.35% and 1.4% in the feed [22, 23], the present study, in which the dietary concentration of MCFA (sum of C10:0 and C12:0) was about 1.25% and 2.5% in groups HI-2.5 and HI-5.0, respectively, showed only a very little impact of dietary inclusion of HI larvae fat on the cecal microbiome. This was evident from the observation that cecal microbial diversity (α and β) did not differ between groups and taxonomic analysis revealed differences in the abundance of only four low-abundance bacterial taxa among groups; namely, the relative abundances of the phylum Actinobacteriota, the class Coriobacteriia, the order Coriobacteriales and the family Eggerthellaceae (all belonging to Actinobacteriota) were lower in group HI-5.0 compared to group HI-2.5, but the abundances of these bacterial taxa did not differ between groups HI-5.0 and HI-0 and between groups HI-2.5 and HI-0. No differences at all were found between groups in the abundances of bacterial families belonging to the two main phyla, Bacteroidota and Firmicutes, together accounting for approximately 99% of all bacteria in the three groups. Considering this and the rather small reduction in the abundances of the abovementioned bacterial taxa suggests that the impact of HI larvae fat on the cecal microbiota structure was very low. One important reason might be that the amount of MCFA reaching the cecum was not sufficient to exert pronounced antimicrobial effects. This assumption is supported by our observation that the ileal digestibility of crude fat was very high (approximately 90% in all groups) and the fact that TG-bound MCFA are hydrolysed and absorbed more rapidly than TG-bound long-chain fatty acids due to easier emulsification and less dependence on pancreatic lipase activity [27, 28]. In agreement with this, an increased ileal fat digestibility has been observed in a recent study with broilers by dietary substitution of soybean oil by HI fat [20]. Nevertheless, reports exist demonstrating that short-term administration of specific MCFA, such as caprylic acid (C8:0) and 1-monoglyceride of capric acid (C10:0), via the feed or the drinking water reduces cecal colonization of pathogenic bacteria, such as Campylobacter jejuni, after artificial infection [22, 23, 57]. However, contradictory results have been reported from another study, in which addition of either non-coated or coated specific MCFA (C6:0, C8:0, C10:0) to the feed at comparable concentrations as in the studies from Solis de los Santos et al. [22, 23] did not reduce cecal Campylobacter colonization in broilers [58]. The reason for these contradictory results is not clear, but differences between studies in the formulation of the MCFA in the feed may explain that the MCFA were less absorbed in the small intestine, thereby, reaching the ceca at higher concentrations in the one study than in the other. Thus, it is possible that the formulation of the diets containing HI larvae fat in the present study allowed an efficient absorption of MCFA in the small intestine, but a low passage of MCFA into the ceca of the broilers.

Owing to differences between bacterial families of different taxa with regard to the use of fermentation substrates and the metabolic pathways engaged in substrate utilization, a shift in the gut microbial community is typically accompanied by an altered profile of microbial fermentation products, such as SCFA, in the gut. Since dietary inclusion of HI larvae fat into the broiler diets did not cause a substantial shift in the bacterial community of the cecum, the unaltered concentrations of total and individual SCFA in the cecal digesta among the three groups are not surprising and are a further indication of the limited impact of HI larvae fat on the cecum microbiota of the broilers. In agreement with our results, a recent study demonstrated that complete replacement of 5% corn oil in the diet with HI larvae fat did not affect the concentrations of total and individual SCFA and viable counts of Clostridium perfringens in the ileum and cecum of broilers during a 30-d-feeding period [16]. In addition, feeding a diet containing about 1% MCFA as lauric acid had no influence on the DNA copy numbers of some beneficial bacteria (Lactobacillus and Bifidobacterium spp.), opportunistic pathogens (Enterobacteriaceae, Escherichia coli) and the pathogen Campylobacter jejuni in the jejunal digesta of broilers [59].

Consistent with our finding that the dietary fats including the HI larvae fat were readily absorbed in the small intestine, the present study revealed some effects of dietary inclusion of HI larvae fat on the hepatic and plasma lipidomes of the broilers. However, the effects of dietary inclusion of HI larvae fat were mainly related to the individual lipid species composition of the lipid classes, whereas the concentrations of the different lipid classes in liver and plasma were either not or only marginally (liver: PE, PG, LPC and LPE; plasma: CE and Cer) affected. In line with this, the sum of total lipids in liver and plasma were not influenced by dietary treatment. With regard to the individual lipid species composition of the different lipid classes, the most striking observation was that the proportions of lipid species with zero or few double bonds (e.g., in the liver: TG species with 0 and 1 double bonds, PC species with 1 and 3 double bonds, PE species with 2 and 3; in plasma: TG species with 0, 1 and 2 double bonds) were increased but those with four or more double bonds (e.g., in the liver: TG species with 5 and 6 double bonds, PC species with 4, 5 and 6 double bonds, PC species with 4, 6 and 8 double bonds; in plasma: TG species with 4, 5, 6, 7 and 8 double bonds) were decreased in group HI-5.0 compared to group HI-0. These alterations in the individual lipid species composition of lipid classes likely reflected the higher percentage of saturated fatty acids in the HI-5.0 broiler diets (59%–60%) than in the HI-2.5 (37%–38%) and the HI-0 (16%–17%) broiler diets, thereby, resulting in a higher intake of saturated fatty acids from HI larvae fat. Although most of the alterations in the plasma and liver lipidomes among the treatment groups were probably attributed to the marked difference in fatty acid composition between soybean oil and HI larvae fat, certain changes in the liver or plasma lipidomes might also reflect specific metabolic effects of the characteristic fatty acids, such as MCFA, contained in the HI larvae fat. For instance, the reason for the slight increase of CE by HI larvae fat might be based on the fact that MCFA, such as capric acid and lauric acid, and the long-chain saturated fatty acid myristic acid and palmitic acid exert hypercholesterolemic effects and increase blood levels of low-density lipoprotein (LDL) [60,61,62], in which CE makes up approximately 50%, at least in humans. Considering this, it appears possible that the raise in plasma CE levels in group HI-5.0 resulted from elevated LDL levels. The latter might be explained by an effect of HI larvae fat on different aspects of lipoprotein metabolism, such as the catabolism of very-LDL (VLDL) particles by lipoprotein lipase in extrahepatic tissues (white adipose tissue, skeletal muscle), hepatic clearance of LDL particles from plasma or the exchange of CE between lipoproteins which is governed by CE transfer protein.

In order to provide a deeper insight into the metabolic effects of HI larvae fat in broilers, differential transcriptome analysis of the liver was carried out. The finding that only 55 annotated genes out of more than 19,000 genes screened were found to be differentially expressed between group HI-5.0 and group HI-0 and the number of genes regulated greater twofold was very low, suggests that the impact of HI larvae fat inclusion on the hepatic transcriptome was modest. As expected, the impact of HI larvae fat inclusion on the hepatic transcriptome was less at the lower dietary inclusion level of HI larvae fat; i.e., only 25 annotated genes were identified as differentially expressed between groups HI-2.5 and HI-0, from which 7 genes were also differentially expressed between groups HI-5.0 and HI-0. Bioinformatic GSEA of these 55 differentially expressed hepatic genes revealed a particular involvement of the encoded proteins in biological processes dealing with sterol metabolism, such as cholesterol metabolic process, cholesterol biosynthetic process, sterol metabolic process, steroid biosynthetic process and lipid biosynthetic process. This observation was likely explained by the finding that several genes involved in sterol synthesis, such as HMGCS1 (encoding 3-hydroxy-3-methylglutaryl-CoA synthase 1), SCAP (encoding SREBF chaperone), AKR1D1 (encoding aldo–keto reductase family 1 member D1), FDFT1 (encoding farnesyl-diphosphate farnesyltransferase 1) and HSD17B7 (encoding hydroxysteroid 17-beta dehydrogenase 7), were amongst the genes differentially expressed between groups HI-5.0 or HI-2.5 and HI-0. Despite that the proteins encoded by HMGCS1, FDFT1 and SCAP are involved in cholesterol synthesis [63], the unaltered levels of free cholesterol and CE in the liver between groups suggests that the regulation of these genes by HI larvae fat inclusion had no impact on hepatic cholesterol levels. Plausible reasons might be that the regulation of these genes occurred either only at the transcriptional level or at a too weak level to substantially alter the abundance of the encoded proteins. In any case, our findings from hepatic transcriptome analysis indicate that the effect of soybean oil replacement with HI larvae fat on intermediary metabolism of broilers is relatively weak.

Conclusion

Comprehensive analysis of the effect of partial and complete replacement of soybean oil with HI larvae fat in broiler diets on growth performance, cecal microbiome, liver transcriptome and liver and plasma lipidomes revealed only very modest, but not any adverse effects of dietary HI larvae fat inclusion. Interestingly, growth performance of the broilers during the starter period was even improved by dietary inclusion of HI larvae fat. The findings of this study suggest that HI larvae fat can be used as an alternative fat source in broiler diets, which makes broiler production more sustainable through the exclusion of soybean oil and the utilization of HI larvae fat from regional production.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The microarray data of this study have been deposited in MIAME compliant format in the NCBI´s Gene Expression Omnibus public repository.

Abbreviations

AID:

Apparent ileal digestibility

AKR1D1:

Aldo–keto reductase family 1 member D1

CE:

Cholesteryl esters

Cer:

Ceramides

DG:

Diglycerides

DM:

Dry matter

FAME:

Fatty acid methyl esters

FC:

Free cholesterol

FDFT1:

Farnesyl-diphosphate farnesyltransferase 1

GO:

Gene ontology

GSEA:

Gene set enrichment analysis

HI:

Hermetia illucens

HMGCS1:

3-Hydroxy-3-methylglutaryl-CoA synthase 1

HSD17B7:

Hydroxysteroid 17-beta dehydrogenase 7

LCFA:

Long-chain fatty acids

LDL:

Low-density lipoprotein

LPC:

Lysophosphatidylcholines

LPE:

Lysophosphatidylethanolamines

MCFA:

Medium-chain fatty acids

NMDS:

Non-metric multidimension distance scaling

OUT:

Operational taxonomic units

PC:

Phosphatidylcholines

PC O:

Ether PC

PCA:

Principal component analysis

PE:

Phosphatidylethanolamines

PE O:

Ether PE

PG:

Phosphatidylglycerols

PI:

Phosphatidylinositols

SCAP:

SREBF chaperone

SCFA:

Short-chain fatty acids

SFA:

Saturated fatty acids

SM:

Sphingomyelins

TG:

Triglycerides

TM:

Tenebrio molitor

VLDL:

Very low-density lipoprotein

WPSA:

World´s Poultry Science Association

References

  1. van Huis A. Potential of insects as food and feed in assuring food security. Annu Rev Entomol. 2013;58:563–83.

    Article  Google Scholar 

  2. Gasco L, Biasato I, Dabbou S, Schiavone A, Gai F. Animals fed insect-based diets: state-of-the-art on digestibility, performance and product quality. Animals (Basel). 2019;9:170.

    Article  Google Scholar 

  3. Sogari G, Amato M, Biasato I, Chiesa S, Gasco L. The potential role of insects as feed: a multi-perspective review. Animals (Basel). 2019;9:119.

    Article  Google Scholar 

  4. EUR-Lex - 32021R1372. Commission Regulation (EU) 2021/1372 of 17 August 2021 amending Annex IV to Regulation (EC) No 999/2001 of the European Parliament and of the Council as regards the prohibition to feed non-ruminant farmed animals, other than fur animals, with protein derived from animals. URL: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:32021R1372. Site accessed at August 10 2022.

  5. Makkar HPS, Tran G, Heuzé V, Ankers P. State-of-the-art on use of insects as animal feed. Anim Feed Sci Technol. 2014;197:1–33.

    Article  CAS  Google Scholar 

  6. Pinotti L, Giromini C, Ottoboni M, Tretola M, Marchis D. Review: Insects and former foodstuffs for upgrading food waste biomasses/streams to feed ingredients for farm animals. Animal. 2019;13:1365–75.

    Article  CAS  Google Scholar 

  7. Pinotti L, Ottoboni M. Substrate as insect feed for bio-mass production. J Insects Food Feed. 2021;7:585–96.

    Article  Google Scholar 

  8. Bovera F, Piccolo G, Gasco L, Marono S, Loponte R, Vassalotti G, et al. Yellow mealworm larvae (Tenebrio molitor, L.) as a possible alternative to soybean meal in broiler diets. Br Poult Sci. 2015;56:569–75.

    CAS  Google Scholar 

  9. Bovera F, Loponte R, Marono S, Piccolo G, Parisi G, Iaconisi V, et al. Use of larvae meal as protein source in broiler diet: Effect on growth performance, nutrient digestibility, and carcass and meat traits. J Anim Sci. 2016;94:639–47.

    Article  CAS  Google Scholar 

  10. Meyer S, Gessner DK, Braune MS, Friedhoff T, Most E, Höring M, et al. Comprehensive evaluation of the metabolic effects of insect meal from Tenebrio molitor L. in growing pigs by transcriptomics, metabolomics and lipidomics. J Anim Sci Biotechnol. 2020;11:20.

    Article  CAS  Google Scholar 

  11. Hartinger K, Greinix J, Thaler N, Ebbing MA, Yacoubi N, Schedle K, et al. Effect of graded substitution of soybean meal by Hermetia illucens larvae meal on animal performance, apparent ileal digestibility, gut histology and microbial metabolites of broilers. Animals (Basel). 2021;11:1628.

    Article  Google Scholar 

  12. Benzertiha A, Kierończyk B, Rawski M, Mikołajczak Z, Urbański A, Nogowski L, et al. Insect fat in animal nutrition – a review. Ann Anim Sci. 2020;20:1217–40.

    Article  Google Scholar 

  13. Schiavone A, Dabbou S, De Marco M, Cullere M, Biasato I, Biasibetti E, et al. Black soldier fly larva fat inclusion in finisher broiler chicken diet as an alternative fat source. Animal. 2018;12:2032–9.

    Article  CAS  Google Scholar 

  14. Benzertiha A, Kierończyk B, Rawski M, Kołodziejski P, Bryszak M, Józefiak D. Insect oil as an alternative to palm oil and poultry fat in broiler chicken nutrition. Animals (Basel). 2019;9:116.

    Article  Google Scholar 

  15. Cullere M, Schiavone A, Dabbou S, Gasco L, Dalle Zotte A. Meat quality and sensory traits of finisher broiler chickens fed with black soldier fly (Hermetia illucens l.) larvae fat as alternative fat source. Animals (Basel). 2019;9:140.

    Article  Google Scholar 

  16. Kim YB, Kim DH, Jeong SB, Lee JW, Kim TH, Lee HG, et al. Black soldier fly larvae oil as an alternative fat source in broiler nutrition. Poult Sci. 2020;99:3133–43.

    Article  CAS  Google Scholar 

  17. Kim B, Bang HAT, Jeong JY, Kim M, Kim KH, Chun JL, et al. Effects of dietary supplementation of black soldier fly (Hermetia illucens) larvae oil on broiler health. J Poult Sci. 2021;58:222–9.

    Article  CAS  Google Scholar 

  18. Sypniewski J, Kierończyk B, Benzertiha A, Mikołajczak Z, Pruszyńska-Oszmałek E, Kołodziejski P, et al. Replacement of soybean oil by Hermetia illucens fat in turkey nutrition: effect on performance, digestibility, microbial community, immune and physiological status and final product quality. Br Poult Sci. 2020;61:294–302.

    Article  CAS  Google Scholar 

  19. Heuel M, Sandrock C, Leiber F, Mathys A, Gold M, Zurbrügg C, et al. Black soldier fly larvae meal and fat can completely replace soybean cake and oil in diets for laying hens. Poult Sci. 2021;100:101034.

    Article  CAS  Google Scholar 

  20. Hartinger K, Fröschl K, Ebbing MA, Bruschek-Pfleger B, Schedle K, Schwarz C, et al. Suitability of Hermetia illucens larvae meal and fat in broiler diets: effects on animal performance, apparent ileal digestibility, gut histology, and microbial metabolites. J Anim Sci Biotechnol. 2022;13:50.

    Article  CAS  Google Scholar 

  21. Kim B, Kim M, Jeong JY, Kim HR, Ji SY, Jung H, et al. Black soldier fly (Hermetia illucens) larvae oil as an alternative fat ingredient to soybean oil in laying hen diets. Anim Biosci. 2022;35:1408–17.

    Article  CAS  Google Scholar 

  22. Solis de los Santos F, Donoghue AM, Venkitanarayanan K, Metcalf JH, Reyes-Herrera I, Dirain ML, et al. The natural feed additive caprylic acid decreases Campylobacter jejuni colonization in market-aged broiler chickens. Poult Sci. 2009;88:61–4.

    Article  CAS  Google Scholar 

  23. Solis de los Santos F, Donoghue AM, Venkitanarayanan K, Dirain ML, Reyes-Herrera I, Blore PJ, et al. Caprylic acid supplemented in feed reduces enteric Campylobacter jejuni colonization in ten-day-old broiler chickens. Poult Sci. 2008;87:800–4.

    Article  CAS  Google Scholar 

  24. Huang Y, Lv H, Song Y, Sun C, Zhang Z, Chen S. Community composition of cecal microbiota in commercial yellow broilers with high and low feed efficiencies. Poult Sci. 2021;100:100996.

    Article  CAS  Google Scholar 

  25. Lv H, Huang Y, Wang T, Zhai S, Hou Z, Chen S. Microbial composition in the duodenum and ileum of yellow broilers with high and low feed efficiency. Front Microbiol. 2021;12: 689653.

    Article  Google Scholar 

  26. Ringseis R, Gessner DK, Eder K. The gut-liver axis in the control of energy metabolism and food intake in animals. Annu Rev Anim Biosci. 2020;8:295–319.

    Article  CAS  Google Scholar 

  27. Hanczakowska E. The use of medium-chain fatty acids in piglet feeding–a review. Ann Anim Sci. 2017;17:967–77.

    Article  CAS  Google Scholar 

  28. Shah ND, Limketkai BN. The use of medium-chain triglycerides in gastrointestinal disorders. Practical Gastroenterol. 2017;160:20–5.

    Google Scholar 

  29. Cobb-Vantress.com. Broiler management guide. URL: https://www.cobb-vantress.com/assets/5c7576a214/Broiler-guide-R1.pdf. Site accessed at August 15 2021.

  30. Keller U, Brandsch C, Eder K. Supplementation of vitamins C and E increases the vitamin E status but does not prevent the formation of oxysterols in the liver of guinea pigs fed an oxidised fat. Eur J Nutr. 2004;43:353–9.

    Article  CAS  Google Scholar 

  31. Short FJ, Gorton P, Wiseman J, Boorman KN. Determination of titanium dioxide added as an inert marker in chicken digestibility studies. Anim Feed Sci Technol. 1996;59:215–21.

    Article  CAS  Google Scholar 

  32. Cobb-Vantress.com. Cobb500 broiler performance & nutrition supplement. URL: https://www.cobb-vantress.com/assets/5a88f2e793/Broiler-Performance-Nutrition-Supplement.pdf. Site accessed at August 15 2021.

  33. VDLUFA (Verband Deutscher Landwirtschaftlicher Untersuchungs- und Forschungsanstalten), Die chemische Untersuchung von Futtermitteln. VDLUFA-Methodenbuch. Band III, Ergänzungslieferungen von 1983, 1988, 1992, 1997, 2004, 2006, 2007, VDLUFA-Verlag, Darmstadt, Germany; 2007.

  34. WPSA (World’s Poultry Science Association). The prediction of apparent metabolizable energy values for poultry in compound feeds. World’s Poult Sci J. 1984;40:181–2.

    Google Scholar 

  35. Kierończyk B, Rawski M, Stuper-Szablewska K, Józefiak D. First report of the apparent metabolisable energy value of black soldier fly larvae fat used in broiler chicken diets. Animal. 2022;16: 100656.

    Article  Google Scholar 

  36. EUR-Lex - 32009R1099. Council Regulation (EC) No 1099/2009 of 24 September 2009 on the protection of animals at the time of killing. URL: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=celex%3A32009R1099. Site accessed at October 11 2022.

  37. Brandt M, Allam SM. Analytik von TiO2 im Darminhalt und Kot nach Kjeldahlaufschluss. Arch Anim Nutr. 1987;37:453–4.

    Google Scholar 

  38. Lagkouvardos I, Kläring K, Heinzmann SS, Platz S, Scholz B, Engel KH, et al. Gut metabolites and bacterial community networks during a pilot intervention study with flaxseeds in healthy adult men. Mol Nutr Food Res. 2015;59:1614–28.

    Article  CAS  Google Scholar 

  39. Meyer S, Gessner DK, Maheshwari G, Röhrig J, Friedhoff T, Most E, et al. Tenebrio molitor larvae meal affects the cecal microbiota of growing pigs. Animals (Basel). 2020;10:1151.

    Article  Google Scholar 

  40. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.

    Article  CAS  Google Scholar 

  41. Jost L. Partitioning diversity into independent alpha and beta components. Ecology. 2007;88:2427–39.

    Article  Google Scholar 

  42. Minchin PR. Simulation of multidimensional community patterns: towards a comprehensive model. Vegetatio. 1987;71:145–56.

    Article  Google Scholar 

  43. Fiesel A, Gessner DK, Most E, Eder K. Effects of dietary polyphenol-rich plant products from grape or hop on pro-inflammatory gene expression in the intestine, nutrient digestibility and faecal microbiota of weaned pigs. BMC Vet Res. 2014;10:196.

    Article  Google Scholar 

  44. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.

    Article  CAS  Google Scholar 

  45. Gessner DK, Schwarz A, Meyer S, Wen G, Most E, Zorn H, et al. Insect meal as alternative protein source exerts pronounced lipid-lowering effects in hyperlipidemic obese zucker rats. J Nutr. 2019;149:566–77.

    Article  Google Scholar 

  46. Ringseis R, Zeitz JO, Weber A, Koch C, Eder K. Hepatic transcript profiling in early-lactation dairy cows fed rumen-protected niacin during the transition from late pregnancy to lactation. J Dairy Sci. 2019;102:365–76.

    Article  CAS  Google Scholar 

  47. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.

    Article  Google Scholar 

  48. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  49. Keller J, Ringseis R, Koc A, Lukas I, Kluge H, Eder K. Supplementation with l-carnitine downregulates genes of the ubiquitin proteasome system in the skeletal muscle and liver of piglets. Animal. 2012;6:70–8.

    Article  CAS  Google Scholar 

  50. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:research0034.1.

    Article  Google Scholar 

  51. Höring M, Krautbauer S, Hiltl L, Babl V, Sigruener A, Burkhardt R, et al. Accurate lipid quantification of tissue homogenates requires suitable sample concentration, solvent composition, and homogenization procedure-a case study in murine liver. Metabolites. 2021;11:365.

    Article  Google Scholar 

  52. Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37:911.

    Article  CAS  Google Scholar 

  53. Höring M, Ejsing CS, Krautbauer S, Ertl VM, Burkhardt R, Liebisch G. Accurate quantification of lipid species affected by isobaric overlap in Fourier-transform mass spectrometry. J Lipid Res. 2021;62:100050.

    Article  Google Scholar 

  54. Höring M, Ejsing CS, Hermansson M, Liebisch G. Quantification of cholesterol and cholesteryl ester by direct flow injection high-resolution fourier transform mass spectrometry utilizing species-specific response factors. Anal Chem. 2019;91:3459–66.

    Article  Google Scholar 

  55. Liebisch G, Fahy E, Aoki J, Dennis EA, Durand T, Ejsing CS, et al. Update on LIPID MAPS classification, nomenclature, and shorthand notation for MS-derived lipid structures. J Lipid Res. 2020;61:1539–55.

    Article  CAS  Google Scholar 

  56. Schlegel G, Ringseis R, Shibani M, Most E, Schuster M, Schwarz FJ, et al. Influence of a rumen-protected conjugated linoleic acid mixture on carcass traits and meat quality in young Simmental heifers. J Anim Sci. 2012;90:1532–40.

    Article  CAS  Google Scholar 

  57. Hilmarsson H, Thormar H, Thráinsson JH, Gunnarsson E, Dadadóttir S. Effect of glycerol monocaprate (monocaprin) on broiler chickens: an attempt at reducing intestinal Campylobacter infection. Poult Sci. 2006;85:588–92.

    Article  CAS  Google Scholar 

  58. Hermans D, Martel A, Van Deun K, Verlinden M, Van Immerseel F, Garmyn A, et al. Intestinal mucus protects Campylobacter jejuni in the ceca of colonized broiler chickens against the bactericidal effects of medium-chain fatty acids. Poult Sci. 2010;89:1144–55.

    Article  CAS  Google Scholar 

  59. Zeitz JO, Fennhoff J, Kluge H, Stangl GI, Eder K. Effects of dietary fats rich in lauric and myristic acid on performance, intestinal morphology, gut microbes, and meat quality in broilers. Poult Sci. 2015;94:2404–13.

    Article  CAS  Google Scholar 

  60. Denke MA, Grundy SM. Comparison of effects of lauric acid and palmitic acid on plasma lipids and lipoproteins. Am J Clin Nutr. 1992;56:895–8.

    Article  CAS  Google Scholar 

  61. Dreon DM, Fernstrom HA, Campos H, Blanche P, Williams PT, Krauss RM. Change in dietary saturated fat intake is correlated with change in mass of large low-density-lipoprotein particles in men. Am J Clin Nutr. 1998;67:828–36.

    Article  CAS  Google Scholar 

  62. Tholstrup T, Ehnholm C, Jauhiainen M, Petersen M, Høy CE, Lund P, et al. Effects of medium-chain fatty acids and oleic acid on blood lipids, lipoproteins, glucose, insulin, and lipid transfer protein activities. Am J Clin Nutr. 2004;79:564–9.

    Article  CAS  Google Scholar 

  63. Shimano H. Sterol regulatory element-binding proteins (SREBPs): transcriptional regulators of lipid synthetic genes. Prog Lipid Res. 2001;40:439–52.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Open Access funding enabled and organized by Projekt DEAL. Lea Schäfer was financially supported by H. Wilhelm Schaumann foundation.

Author information

Authors and Affiliations

  1. Institute of Animal Nutrition and Nutrition Physiology, Justus Liebig University Giessen, Heinrich-Buff-Ring 26-32, 35392, Giessen, Germany

    Lea Schäfer, Sarah M. Grundmann, Erika Most, Klaus Eder & Robert Ringseis

  2. Institute of Nutritional and Food Sciences, Molecular Food Technology, University of Bonn, Friedrich-Hirzebruch-Allee 7, 53115, Bonn, Germany

    Garima Maheshwari

  3. Institute of Clinical Chemistry and Laboratory Medicine, University Hospital of Regensburg, Franz-Josef-Strauss-Allee 11, 93053, Regensburg, Germany

    Marcus Höring & Gerhard Liebisch

Authors
  1. Lea Schäfer
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sarah M. Grundmann
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Garima Maheshwari
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Marcus Höring
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Gerhard Liebisch
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Erika Most
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Klaus Eder
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Robert Ringseis
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

RR and KE designed research and coordinated the study. LS performed the experiment, the laboratory analysis and statistical analysis of all data. SMG performed principal component analysis. GM performed analysis of the microbiome data. MH and GL performed lipidomics analysis. RR performed bioinformatic analysis of transcriptomic data. LS, KE and RR interpreted the data. LS wrote the paper. KE and RR reviewed and edited the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Robert Ringseis.

Ethics declarations

Ethics approval and consent to participate

The animal study trial was approved by the Animal Welfare Officer of the Justus Liebig University Giessen (approval no.: JLU 786_M).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1. Table S1

Characteristics of Gallus gallus gene-specific primers used for qPCR analysis. Table S2 Concentrations of amino acids in the broiler diets. Table S3 Operational taxonomic units (OTU) identified in cecum digesta of the broilers. Table S4 List of differentially expressed transcripts in the liver of broilers between group HI-5.0 vs. HI-0. Table S5 List of differentially expressed transcripts in the liver of broilers between group HI-2.5 vs. HI-0. Table S6 qPCR validation of microarray data. Table S7 Individual lipid species composition of triacylglycerols (TG), phosphatidylcholine (PC) and phosphatidylethanolamine (PE) in the liver of broilers fed diets with either 0% (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d. Table S8 Individual lipid species composition of cholesteryl esters (CE), phosphatidylcholine (PC) and triacylglycerols (TG) in plasma of broilers fed diets with either 0% (HI-0), 2.5% (HI-2.5) or 5.0% (HI-5.0) Hermetia illucens (HI) larvae fat for 35 d.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schäfer, L., Grundmann, S.M., Maheshwari, G. et al. Effect of replacement of soybean oil by Hermetia illucens fat on performance, digestibility, cecal microbiome, liver transcriptome and liver and plasma lipidomes of broilers. J Animal Sci Biotechnol 14, 20 (2023). https://doi.org/10.1186/s40104-023-00831-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-023-00831-6

Keywords

Journal of Animal Science and Biotechnology

ISSN: 2049-1891