DOI: https://doi.org/10.21203/rs.3.rs-61552/v1
Background: Weaning is associated with growth decline and morphological changes in the gastrointestinal tract of the pig. Insects, such as the Black soldier fly larvae (BSFL), are suggested as a sustainable novel protein source in diets for pigs. They contain chitin, medium chained fatty acids, and anti-microbial peptides, which could improve the gastrointestinal function and health in the post-weaning pig. The objective of this study was therefore to investigate the effect of increasing inclusion of full-fat BSFL in diets for post-weaning pigs on growth performance parameters, digestibility, gut morphology, and the microbial community in the colon. Eighty crossbred weanling pigs were weaned at approximately 32 days of age, with an average weaning weight of 10.6 ± 0.8 kg. For four weeks, the pigs were fed one of four dietary treatments: a control diet or one of three diets containing an increasing amount of full-fat BSFL meal: 4.76% (BSFL5), 9.52% (BSFL10), and 19.06% (BSFL20).
Results: Increased level of full-fat BSFL in the diet did not affect growth performance or fecal consistency. A reduction in apparent total tract digestibility (ATTD) of crude protein (P = 0.035) was found for increased inclusion of BSFL, whereas the ATTD of crude fat increased (P<0.001). Jejunal, ileal, or colonic morphometry was not affected by the BSFL inclusion.
Some changes in the colon microbiota were observed, but no differences in the short-chain fatty acid concentrations were detected between the dietary treatments. At the phylum level, dietary treatment affected the relative abundance of Bacteroidota, Firmicutes, Campilobacteria, and Thermoplasmatota, but there was no clear pattern relationship with the BSFL inclusion level. At the genus level, the inclusion of BSFL in the diet reduced the relative abundance of Lactobacillus (P = 0.015) compared to the control.
Conclusions: Collectively, the results indicate that up to 19.06% of full-fat BSFL meal could be included in a balanced diet for PW pigs without affecting growth performance, gut function, or health.
The world is facing a growing food demand by an ever-increasing number of people. To meet this challenge, new resources must be considered by using innovative solutions. Insects have been proposed as a high quality, efficient, and sustainable alternative protein source [1]. Black soldier fly larvae (BSFL), Hermetia illucens, is an easily reared species, capable of efficient conversion of a wide range of organic materials [2], and do not accumulate pesticides or mycotoxins [3, 4]. The BSFL has an overall favorable amino acid profile but is limiting in sulfur-containing amino acids and contains a high amount of minerals, especially calcium [5].
A well-functioning gastrointestinal tract (GIT) is of importance to the overall growth performance and health of pigs in all stages, and especially for the newly weaned pigs. If not defatted, BSFL is high in fat, especially in lauric acid (12:0) [5] which is categorized as a medium-chain fatty acid (MCFA) with antimicrobial effects, especially against gram-positive bacteria [6]. Insects and insect larvae also contain chitin, a dietary polysaccharide that can function as a prebiotic and an immunostimulant [7]. Finally, anti-microbial peptides, which are part of the insect immune system, are effective antimicrobial agents with low risk of development of bacteria resistance. Therefore, the anti-microbial peptides have good potential as health promoters in livestock, even though there is limited information about the in vivo effects [8].
Based on the content of these constituents, insects are interesting to investigate both as a feed ingredient and as a functional ingredient with potential health beneficial effects. Recently there has been a high focus devoted to the potential of insect meals in diets for monogastric farm animals [9–12], however, there is a scarcity of literature on the use of full-fat BSFL in diets for pigs. The objective of this study was to investigate the effect of increased inclusion of full-fat BSFL in diets for post-weaning (PW) pigs, focusing on growth performance parameters, digestibility, gut morphology, and microbial community in the colon.
A 27 days experiment was performed in February 2019 at the Center for livestock production (SHF, NMBU, Ås, Norway), which is an animal experimental unit approved by the National Animal Research Authority (permit no. 174). All pigs were handled under the applicable laws and regulations controlling experiments with live animals in Norway (regulated by the “Animal Welfare Act” and “The Norwegian Regulation on Animal Experimentation” derived from the “Directive 2010/63/EU on the protection of animals used for scientific purposes.
Eighty crossbred [Norwegian Landrace x Yorkshire z-line) x Duroc)] weanling pigs, selected from eleven litters (four or eight pigs from each litter) were included in the experiment. Pigs who had been under medical treatment in the nursing period were excluded. All piglets had access to the sows feed during the nursing period. The average weaning age was 32.8 ± 1.6 SD days and the average weaning weight 10.6 ± 0.8 SD kg. Pigs were blocked by litter, weight, and sex and assigned to one of four dietary treatments. Pigs were distributed in pens of four, and from a total of five pens per treatment, three of the pens per treatment were installed with rubber mats. The remaining pens had access to wood shavings as bedding material. The pen size was 1.6 m2. The room temperature was logged every morning, and the average temperature for the experimental period was 21.6 ± 1.4 °C.
BSFL meal was produced at HiProMine S.A., Poznan, Poland. The BSFL feed was normalized in terms of dry matter (DM) content by the addition of wheat middlings (17%) to fresh vegetables and fruit mix, consisting of apples (15%) carrots (50%) potatoes (15%) and cabbage (20%) and established at the level of 22% DM. Fresh vegetable pre-consumer waste was ground (2000 rpm/1 min, (HPM milling system, 55 kW, Poland) to pass a 2 mm screen and offered ad libitum to the BSFL. Substrates were not contaminated by any animal products in accordance with EC regulation (no 1069/09). At the prepupal stage (10th day of rearing), larvae were harvested, sieved through a 3 mm screen, and washed with water on drum separator at 90 ⁰C for 10 minutes (HPM cleaning system, Poland).
The dietary treatments included a control diet and three experimental diets with increased inclusion of BSFL; 4.76% (BSFL5), 9.52% (BSFL10), and 19.06% (BSFL20) (Table 1). Diets were formulated in collaboration with Felleskjøpet Fôrutvikling AS using their optimization least-cost program based on the Dutch energy evaluation system [13]. Diets were formulated on net energy and standardized ileal digestibility (SID) values to be isoenergetic, balanced for digestible amino acids, and to meet or exceed the nutritional requirements for this age pigs [14]. Literature values from Finke [5] for the amino acid content in the BSFL meal were used in the diet formulation. Digestibility coefficients for the amino acids in the BSFL meal were set to 83%. Yttrium oxide (Y2O3) was included as an inert marker in the diets for digestibility calculations. Pelleted diets were produced by the Center for Feed Technology (ForTek, NMBU, Ås, Norway). The feed mash was ground in a Münch hammer mill (HM 21.115, Wuppertal, Germany) fitted with a 3 mm screen before pelleting. The mash was steam conditioned at 82 °C in a double-pass pellet-press conditioner (Münch-Edelstahl, Germany) before pelleting (Münch-Edelstahl, Germany, 2 × 17 kW) through a 3.5 mm die with a production rate of 700 kg/h. The pigs had ad-libitum access to the experimental diets immediately PW through automatic feeders (FRH-2 Domino A/S, Tørring, Denmark) with 43 cm feeding space. The automatic feeders were checked daily and refilled when needed. Feed residues were registered weekly and average daily feed intake (ADFI) per pen was calculated. Clean drinking water was always available from a drinking nipple next to the feeder.
Dietary treatments |
||||
---|---|---|---|---|
Ingredients, g/kg as fed |
Control |
BSFL5 |
BSFL10 |
BSFL20 |
Wheata |
507.7 |
502.0 |
496.6 |
485.1 |
Barleyb |
200.0 |
200.0 |
200.0 |
200.0 |
Oatsc |
50.0 |
50.0 |
50.0 |
50.0 |
BSFL meald |
⎯⎯ |
47.6 |
95.2 |
190.6 |
Soybean meale |
70.8 |
59.7 |
48.2 |
25.5 |
Soy protein concentratef |
36.1 |
27.1 |
18.1 |
⎯⎯ |
Fish mealg |
34.1 |
25.5 |
17.1 |
⎯⎯ |
Rapeseed oil |
32.8 |
25.1 |
17.4 |
1.3 |
AkoFeed Gigant 60h |
11.0 |
8.1 |
5.1 |
⎯⎯ |
Rapeseed meali |
10.0 |
10.0 |
10.0 |
10.0 |
Monocalcium phosphate |
12.9 |
12.0 |
11.2 |
9.5 |
Limestone |
8.2 |
6.7 |
5.1 |
2.1 |
Sodium chloride |
5.3 |
5.4 |
5.5 |
5.8 |
Selenium premix |
0.9 |
0.9 |
0.8 |
0.8 |
Iron(II)fumarate |
0.4 |
0.4 |
0.4 |
0.4 |
Micro-mineral premixj |
2.0 |
2.0 |
2.0 |
2.0 |
Vitaminsk |
3.0 |
3.0 |
3.0 |
3.0 |
L-Lysin |
6.8 |
6.8 |
6.9 |
6.9 |
L-Methionine |
2.4 |
2.5 |
2.6 |
2.8 |
L-Threonine |
2.8 |
2.9 |
2.9 |
3.0 |
L-Valin |
1.1 |
0.8 |
0.5 |
0.0 |
L-Tryptophan |
0.9 |
0.8 |
0.7 |
0.6 |
Betaine |
0.7 |
0.7 |
0.7 |
0.7 |
Yttrium(III)oxide |
0.1 |
0.1 |
0.1 |
0.1 |
Calculated CP content |
181.0 |
181.0 |
181.0 |
181.0 |
Ratio CP from BSFL (% of total CP) |
0 |
10 |
20 |
39 |
Ratio CF from BSFL (% of total CF) |
0 |
19 |
39 |
69 |
aChemical composition per kg: 867 g DM, 14 g ash, 124 g CP, 488 g starch, 14 g CF, 116 g NDF, 45 g ADF, 16 MJ. | ||||
bChemical composition per kg: 858 g DM, 17 g ash, 92 g CP, 564 g starch, 11 g CF, 206 g NDF, 82 g ADF, 16 MJ. | ||||
cChemical composition per kg: 858 g DM, 23 g ash, 88 g CP, 408 g starch, 44 g CF, 270 g NDF, 150 g ADF, 17 MJ. | ||||
dHiProMine S.A., Poznanska Str, Poland. Chemical composition per kg: 905 g DM, 84 g ash, 380 g CP, 289 g CF, 17.4 g Ca, 2.4 g Mg, 7.8 g total P, 2.7 g. Amino acids per kg: 27.8 g alanine, 16.7 g arginine, 33.2 g aspartic acid, 2.7 g cysteine, 43.8 g glutamic acid, 18.2 g glycine, 9.0 g histidine, 16.4 g isoleucine, 30.4 g leucine, 21.9 g lysine, 6.5 g methionine, 15.0 g phenylalanine, 20.6 g proline, 14.2 g serine, 15.1 g threonine, 20.3 g tyrosine, 21.5 g valine. Lauric acid (C12:0): 86.7 g/kg. | ||||
eNon-GMO soybean meal, Denofa AS, Fredrikstad, Norway. Chemical composition per kg: 876 g DM, 53 g ash, 426 g CP, 11 g CF, 156 g NDF, 128 g ADF, 17 MJ. | ||||
fAX3 Gastric, TripleA a/s, Hornsyld, Denmark. | ||||
gNordsildmel AS, Egersund, Norway. Chemical composition per kg: 920 g DM, 147 g ash, 679 g CP, 101 g CF, 19 MJ. | ||||
hAAK AB, Malmö, Sweden. | ||||
iExpeller-pressed rapeseed cake, Mestilla, UAB, Klaipeda Lithuania. Chemical composition per kg: 928 g DM, 53 g ash, 342 g CP, 100 g CF, 223 g NDF, 179 g ADF, 20 MJ. | ||||
j“Mikro-Svin”; provided per kilogram of diet: 475 mg Ca; 3.4 mg Mg; 13.2 mg S; 120 mg Fe; 60 mg Mn; 120 mg Zn; 26 mg Cu; 0.6 mg I. | ||||
kProvided per kilogram of diet: 0.7 g Vitamin A; 1.2 g Vitamin E v5; 0.8 g Vitamin ADKB mix; 0.3 g Vitamin C (Stay C 35%). |
The farm had an ongoing issue with edema disease and symptoms resulted in antibiotic treatment of five pigs on days 8–10 PW. After one pig died without any registered symptoms on day 11 PW, all pigs were treated with intramuscular antibiotic injections (Borgal vet., Ceva Santé Animale, France) for three consecutive days (11–13 PW). After treatment, all pigs appeared healthy throughout the experiment.
Fecal consistency was assessed daily and registered for all pens based on the four category scale developed by Pedersen and Toft [15]. The daily fecal score was registered as a pen average with 0.25 intervals on the scale. All pigs were weighed weekly and average daily gain (ADG) and feed conversion ratio (FCR) was calculated per pen. Fecal samples were also collected every week for the determination of the fecal DM. An equal amount of feces from each pig was pooled for the pen before oven drying at 103 °C for 24 h. Individual fecal samples were collected on days 21, 22, 25, 26, and 27 for the determination of apparent total tract digestibility (ATTD). Individual samples from all days were pooled, freeze-dried, and ground at 0.5 and 1 mm (Retsch ZM 100 centrifugal mill, Haan, Germany) before chemical analysis. Apparent digestibility of nutrients was calculated as described by Maynard and Loosli [16].
Only pigs from the pens with rubber mats (n = 3 pens per treatment, a total of 47 pigs) were included in the terminal sampling. Pigs were fasting from the evening before but had access to feed three hours before euthanasia. Euthanasia was done using a captive bolt pistol followed by exsanguination. Immediately after exsanguination, the abdominal cavity was opened, and the GIT removed. Intestinal content was collected from the oral part of the jejunum and the top of the colon spiral. Tissue from oral jejunum, aboral ileum, and top of the colon spiral were collected for histological assessment. Intestinal content from the last two meters of the small intestine was collected and stored at -20 °C for the determination of apparent ileal digestibility (AID). This sample was unfortunately not collected from seven of the pigs (three control pigs, two pigs fed BSFL10, and two pigs fed BSFL20). The ileal contents were later freeze-dried, homogenized using a batch mill (A11 basic Analytical mill, IKA®, England) and chemically analyzed. Liver weight was recorded to calculate liver index (liver weight, kg / live body weight, kg * 100).
Dietary treatments |
||||
---|---|---|---|---|
Nutrients, g/kg DM |
Control |
BSFL5 |
BSFL10 |
BSFL20 |
Dry matter, g/kg |
890.3 |
888.9 |
887.7 |
893.0 |
Ash |
54.3 |
52.3 |
52.5 |
52.9 |
Crude protein |
196.8 |
196.1 |
199.8 |
206.8 |
Starch |
512.5 |
523.8 |
504.1 |
488.0 |
Crude fat |
69.3 |
79.7 |
78.6 |
89.1 |
NDF |
124.1 |
122.9 |
118.3 |
131.0 |
ADF |
52.9 |
52.0 |
50.2 |
53.8 |
Gross energy (MJ/kg) |
19.5 |
19.6 |
19.8 |
19.9 |
Phosphorus |
7.4 |
7.3 |
7.6 |
7.0 |
Pooled feed samples, collected during the experiment for each diet, were ground at 0.5 and 1 mm (Fritsch Pulverisette 19 cutting mill, Idar-Oberstein, Germany) before the chemical composition of nutrients were analyzed in triplicates. The analyzed chemical composition of the dietary treatments can be found in Table 2. Analyzed amino acid and fatty acid composition of the diets are shown in Table 3 and Table 4, respectively. The chemical analyses were performed by the LabTek group at the Department of Animal and Aquacultural Science, NMBU, Ås, Norway. In brief, DM was determined by drying to constant weight at 103 °C ± 2 °C [17], and ash was determined by complete combustion at 550 °C for at least 4 h [18]. Gross energy (GE) content was determined using a PARR 6400 Automatic Isoperibol Calorimeter (Parr Instruments, Moline, IL, USA) according to ISO 9831 [19]. Crude protein (CP) was analyzed with the Kjeldahl method (AOAC Official Method 2001.11), using a Digestor™ 2520 (FOSS Analytical, Hillerød, Denmark) and the Kjeltec™ 8400 analyzer (FOSS Analytical, Hillerød, Denmark). Crude fat (CF) was analyzed using Accelerated Solvent Extraction (ASE 350, Thermo Fisher Scientific Inc.). Extraction was conducted with 80% petroleum ether and 20% acetone at 125 °C. Starch was determined using an enzymatic-colorimetric method according to McCleary et al. [20], with some modifications. In brief, starch was degraded with heat-stable α-amylase and amyloglucosidase-enzymes to glucose. Glucose concentration was then determined using a spectrophotometer (RX Daytona +, Randox Laboratories Ltd., UK). Acid detergent fiber (ADF) and neutral detergent fiber (NDF) were analyzed according to the manufacturer’s methods using the Ankom200 Fiber Analyzer (ANKOM Technology, Macedon, New York, USA). Amino acids in diets and ileal samples were analyzed according to Commission Regulation (EC) No 152/2009. Amino acids were analyzed on a Biochrom 30 + Amino Acid Analyzer with an autosampler (Biochrom Ltd., Cambridge, UK). Tryptophan was analyzed in diets on a Dionex Ultimate 3000 HPLC system (Dionex Softron GmbH, Germering, Germany) with a Shimadzu RF-535 fluorescence detector (Shimadzu Corporation, Kyoto, Japan). The fatty acid composition of the diets was analyzed according to O'fallon et al. [21] by synthesizing the fatty acids to fatty acid methyl esters (FAME), in which concentrations were determined using a Trace GC Ultra gas chromatograph (Thermo Fisher Scientific, US). Total phosphorus was analyzed after combustion and acid digestion (Commission Regulation (EC) No 152/2009) using a commercial spectrophotometric kit (PH8328, Randox laboratories, County Antrim, UK). Yttrium concentration was determined after acid decomposition in a microwave digestion system (Start D, Milestone Srl, Italy), using a microwave plasma atomic emission spectrometer (MP-AES 4200, Agilent Technologies, US).
Dietary treatments |
||||
---|---|---|---|---|
Indispensable AA, g/16 g Na |
Control |
BSFL5 |
BSFL10 |
BSFL20 |
Arginine |
4.9 |
4.9 |
4.6 |
4.2 |
Histidine |
2.1 |
2.2 |
2.2 |
2.3 |
Isoleucine |
3.2 |
3.3 |
3.2 |
3.3 |
Leucine |
6.0 |
6.0 |
5.9 |
5.9 |
Lysine |
6.7 |
7.2 |
6.6 |
6.9 |
Methionine |
2.3 |
2.4 |
2.3 |
2.5 |
Phenylalanine |
3.8 |
3.7 |
3.6 |
3.8 |
Threonine |
4.1 |
4.3 |
4.0 |
4.3 |
Tryptophan |
1.4 |
1.2 |
1.4 |
1.4 |
Valine |
4.0 |
4.1 |
3.9 |
4.1 |
Dispensable AA, g/16 g N |
||||
Alanine |
3.3 |
3.6 |
3.7 |
4.1 |
Aspartic acid |
6.8 |
6.9 |
6.5 |
6.4 |
Cysteine |
1.3 |
1.3 |
1.2 |
1.1 |
Glutamic acid |
19.9 |
19.1 |
18.7 |
18.2 |
Glycine |
3.3 |
3.4 |
3.3 |
3.3 |
Proline |
5.9 |
5.9 |
6.0 |
6.3 |
Serine |
3.7 |
3.7 |
3.6 |
3.5 |
Tyrosine |
1.7 |
1.9 |
2.2 |
3.1 |
Total amino acids, g/16 g N |
83.1 |
83.8 |
81.6 |
83.5 |
aDetermined using water-corrected molecular weights |
Dietary treatments |
||||
---|---|---|---|---|
Fatty acids, g/kg DM |
Control |
BSFL5 |
BSFL10 |
BSFL20 |
C12:0 |
0.40 |
6.11 |
11.7 |
23.1 |
C14:0 |
0.42 |
1.35 |
2.33 |
4.25 |
C15:0 |
0.04 |
0.06 |
0.07 |
0.09 |
C16:0 |
13.3 |
12.4 |
11.4 |
9.96 |
C17:0 |
0.07 |
0.08 |
0.10 |
0.13 |
C18:0 |
4.31 |
3.49 |
2.76 |
1.54 |
C20:0 |
0.24 |
0.20 |
0.15 |
0.09 |
C21:0 |
0.01 |
0.09 |
0.28 |
0.32 |
C22:0 |
0.14 |
0.11 |
0.08 |
0.04 |
C23:0 |
0.02 |
0.01 |
0.01 |
- |
C24:0 |
0.05 |
0.04 |
0.03 |
0.02 |
Sum saturated fatty acids |
19.02 |
23.95 |
28.89 |
39.57 |
C14:1 |
- |
0.02 |
0.04 |
0.08 |
C16:1 |
0.40 |
0.58 |
0.76 |
1.13 |
C20:1 |
0.89 |
0.73 |
0.55 |
0.22 |
C24:1 |
0.11 |
0.08 |
0.04 |
0.03 |
Sum monounsaturated fatty acids |
1.40 |
1.41 |
1.39 |
1.46 |
C18:1n9t |
0.46 |
0.37 |
0.21 |
0.04 |
C18:1n9c |
25.9 |
22.8 |
19.0 |
12.2 |
C18:2n6c |
21.1 |
22.0 |
21.7 |
21.1 |
C18:3n3 |
4.56 |
4.15 |
3.56 |
2.39 |
C18:3n6 |
0.13 |
0.11 |
0.08 |
0.02 |
C20:2 |
0.13 |
0.11 |
0.09 |
0.05 |
C20:4n6 |
0.02 |
0.02 |
0.01 |
- |
C20:5n3 |
0.30 |
0.23 |
0.16 |
- |
C21:1n9 |
0.05 |
0.05 |
0.04 |
0.02 |
C22:2 |
0.04 |
0.02 |
0.02 |
0.02 |
C22:5n3 |
0.03 |
0.02 |
0.01 |
- |
C22:6n3 |
0.36 |
0.30 |
0.20 |
- |
Sum polyunsaturated fatty acids |
53.06 |
50.17 |
45.04 |
35.89 |
Trypsin and lipase activities were analyzed in jejunal content. 1.5 ml (600 µl for lipase analysis) ice-cold Milli-Q water were added to approximately 100 mg of the jejunal content, homogenized using a bead mill (TissueLyser, Qiagen Retsch, Haan, Germany) and sonicated in an ice-cold bath for three minutes (T 460/H, Elma Schmidbauer GmbH, Ransbach-Baumbach, Germany). After centrifugation at 21,100 × g for 10 min at 4 °C, the supernatant was collected, aliquoted, stored at -80 °C, and used for the analysis of lipase, trypsin, and total protein. Total protein concentration was determined in microtiter assay according to the Quick Start™ Bradford Protein Assay protocol (Bio-Rad Laboratories, Oslo, Norway). Absorbance was measured using a SpectraMax M2e Microplate Reader (Molecular Devices, LLC., San Jose, USA). Lipase and trypsin activity were analyzed using commercial kits (Lipase Activity Assay Kit III, MAK048-1KT, fluorometric, Sigma-Aldrich, Merck KGaA; Trypsin Activity Assay Kit, ab102531, colorimetric, Abcam), according to the manufacturer’s protocols.
Short-chain fatty acids (SCFA) were analyzed in colon content. Samples were thawed on ice and 500 mg of the colon content were mixed with 500 µl ice-cold internal standard solution (2-methyl valeric acid in 5% formic acid), before sonication for 5 min in cold water. After centrifugation for 15 min at 4 °C with 15000 × g, the supernatant was transferred to a spin column (45 kDa; VWR International, USA), and again centrifuged with the same parameters. SCFA concentration was determined by capillary gas chromatography on a stabilwax-DA, 30 m × 0.25 mm × 0.25 µm capillary column (Restek corporation, PA, USA) installed on a Trace 1300 gas chromatograph equipped with an AS 1310 autosampler, split injection, a flame ionization detector and Chromeleon software (Thermo Scientific, MA, USA). The initial oven temperature was 90 °C, held for 2 min, followed by a temperature increase of 10 °C/min to 150 °C and 50 °C/min to 250 °C and then held for 1 min. Helium was used as the carrier gas at a flow rate of 3 mL/min. The injector temperature was set at 260 °C and the detector temperature was set at 275 °C.
Jejunal, ileal, and colon morphology were blindly assessed and villus heights (VH) and crypt depths (CD) measured by Aquamedic AS, Oslo, Norway. Gross pathological observations were recorded before tissue sampling. Samples were fixated in 10% formalin up to 48 h before processing following standard histological methods for gut tissue. Sections were stained with hematoxylin and eosin (H&E) and evaluated by light microscopy, where the evaluation was done on morphological characteristics such as epithelial cell and barrier morphology and integrity, crypt changes such as hyperplasia, dilation or abscessation, degenerative and inflammatory mucosal changes including increased numbers of intraepithelial lymphocytes (IELs) and infiltration by leucocytes. Methodologies of the evaluation protocol are described by Day et al. [22] and Pérez de Nanclares et al. [23]. The morphological characteristics evaluated were graded using a semi-quantitative scoring system where score 0 is normal, score 1 represent mild changes, and score 2, 3, and 4 represent moderate, marked, and severe changes, respectively. VH and CD measurements were made on scanned whole-section images of the respective tissues captured using the PreciPoint M8 Microscope and Scanner (PreciPoint, Freising, Germany), and obtained using the ViewPoint software (PreciPoint, Freising, Germany). A minimum of three well-oriented crypts and villi were measured from each of the sections.
Total DNA was extracted from bacteria in approximately 190 mg of colon content using QIAamp Fast DNA Stolen Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s guidelines. The DNA concentration was determined using NanoDrop™ 8000 spectrophotometer (Thermo Fisher Scientific, Waltham, Massachusetts, USA). All samples were stored at -20 ºC until further analysis. Before preparing for 16S rRNA sequencing, samples were normalized to 10 ng/µL. The V3-V4 regions of the bacterial 16S rRNA gene were amplified using the primers Pro341f (5'-CCTACGGGNBGCASCAG-3') and Pro805r (5'-GACTACNVGGGTATCTAATCC-3'). The library preparation was conducted using the Miseq Reagent Kit V3 (Illumina, San Diego, California, USA) according to the Illumina 16 S Metagenomic Sequencing Library Preparation protocol (Illumina, San Diego, California, USA). For indexing, Nextera XT index kit V2 was used (Illumina, San Diego, California, USA). After the indexing reaction, all samples were measured by Qubit Fluorometer (Invitrogen, Carlsbad, California, USA) using Qubit™ 1X dsDNA HS Assay Kit (Invitrogen, Carlsbad, California, USA). An equal amount of each sample was pooled together, and spiked with 5% PhiX Control (Illumina, San Diego, Waltham, MA, USA). For sequencing, 10 pm of the pooled sample was loaded to a flow cell. The sequencing analysis was performed on the Miseq System (Illumina, San Diego, California, USA). The clustering density was 1256 k/mm2 and 88.3% of clusters were passing filter.
Raw sequence data were analyzed using DADA2 v. 1.12.1 [24] in R v. 3.5.0 [25]. Default parameters were used if other is not specified. In brief, primers were removed, and the forward reads were truncated at 280 bp and the reverse reads at 250 bp. Max expected errors were set to 7. This allowed 73% of the reads to pass the quality control. Error rates were estimated, and the core sample inference algorithm applied. The denoised read pairs were merged, and chimeras removed. The Silva v. 138 database [26, 27] was used as a reference database for the taxonomy assignment. A phyloseq object was built with the phyloseq v.1.26.1 for further analyses [28]. Figures were made with ggplot2 v.3.2.1 [29].
Outliers in the data were identified using the interquartile range (IQR) method. Outliers were defined as > Q3 + 3* IQR or < Q1 + 3*IQR. Statistical analysis was performed using R v.3.6.1 [25] in Rstudio v.1.2.5001 [30]. A mixed procedure was run, using the lme4 1.1–21 [31] and lmerTest 3.1-1 [32] packages. The following model was used for all individual data: Yijklmn = µ + dieti + sexj + beddingk + litterl + penm + εijklmn where Y is one observation on pig n; µ is the intercept; dieti is the fixed dietary treatment effect (i = 1:4); sexj is the fixed effect of the sex of the pig (j = 1,2); beddingk is the fixed effect of bedding material, rubber mat or wood shavings (k = 0,1); litterl is the random effect of the lth litter (l = 1:11); penm is the random effect of the pen (m = 1:20) and εijklmn is a random residual. The bedding was not included in the model when analyzing the terminal sampling parameter, as it was constant. For all pen-level parameters, the following model was used with the lm function in the stats package [25]: Yikn = µ + dieti + beddingk + εn where Y now is one observation on pen n. Results are presented as mean values and standard error of the mean (SEM). In addition, linear regression with the dietary inclusion level of BSFL as the only explanatory variable was run with the lm function for all dependent variables. For the effects of sex and bedding material, only significant results are presented. Effects are considered statistically significant when P ≤ 0.05 and tendencies are defined as P-values between 0.05 and 0.10.
To assess diversity in the microbial communities within pigs, Shannon alpha-diversity indices were calculated at the amplicon sequence variant (ASV) level. A comparison of the Shannon indices between dietary treatments was done using the Kruskal-Wallis Rank Sum Test in the stats package [25]. The beta-diversity (between pig variation in microbial communities) was assessed by principal coordinate analyses (PCoA) with the Bray-Curtis, unweighted, and weighted UniFrac distance matrices. Beta dispersion (variances) was calculated within each dietary group and a permutation-based test of multivariate homogeneity applied (vegan package v.2.6-6; Oksanen et al. [33]). Fulfilling the assumption of homogeneity in group variance, a PERMANOVA test was performed on the distance matrices with the dietary group as the independent variable. Pairwise PERMANOVA tests were performed to compare beta diversity between dietary treatments. P-values were adjusted for multiple testing with the method by Benjamini and Hochberg [34], also known as the false discovery rate.
The relative abundance was calculated, and the statistical difference in relative abundance between dietary treatments was analyzed with the Kruskal-Wallis test with dietary treatment as the explanatory variable.
To test for a statistical difference in relative abundance between dietary treatments, the Kruskal-Wallis test with dietary treatment as the explanatory variable was applied. Two-sample Wilcoxon tests, also known as the Mann-Whitney test, was used for pairwise comparison between the diets if the Kruskal Wallis test gave a significant effect of dietary treatment. P-values were corrected for multiple testing with the Benjamini and Hochberg [34] method.
To test for covariation between the microbiota and SCFA profile, dissimilarity indices with the Bray-Curtis indices were calculated separately for the microbiota and SCFA data with the vegdist function, and the covariation tested with the mantel function in the vegan package [33].
Inclusion of BSFL meal did not significantly affect ADG, ADFI, or FCR in the experimental period (Table 5). The ADFI ranged from 650 to 712 g, whereas the ADG ranged from 496 to 552 g for pigs fed increased BSFL during the total period (0–27 days). There was a tendency for lower FCR for day 0–14 PW with increasing inclusion levels of the BSFL meal (R2 = 0.194). Pigs in pens with wood shavings had significantly higher ADG (P = 0.040) for day 14–27 PW.
Dietary treatments | P-value | |||||||
---|---|---|---|---|---|---|---|---|
Day PW | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression | |
0–14 | ADFI | 485 | 448 | 459 | 464 | 9 | 0.644 | 0.714 |
ADG | 432 | 392 | 425 | 427 | 7 | 0.177 | 0.702 | |
FCR | 1.15 | 1.14 | 1.08 | 1.09 | 0.01 | 0.097 | 0.059 | |
14–27 | ADFI | 943 | 867 | 918 | 917 | 15 | 0.389 | 0.941 |
ADG | 683 | 608 | 661 | 636 | 15 | 0.262 | 0.559 | |
FCR | 1.39 | 1.44 | 1.39 | 1.45 | 0.01 | 0.538 | 0.426 | |
0–27 | ADFI | 712 | 650 | 680 | 682 | 11 | 0.324 | 0.775 |
ADG | 552 | 496 | 538 | 528 | 9 | 0.114 | 0.758 | |
FCR | 1.27 | 1.31 | 1.26 | 1.29 | 0.01 | 0.490 | 0.838 |
There were no significant differences for pH in stomach (N = 47; ANOVA P = 0.920; regression P = 0.658; SEM = 0.2) or pH in jejunum (N = 47; ANOVA P = 0.603; regression P = 0.612; SEM = 0.1) between the dietary treatments. Numerical means for pH measurements were; control: 4.5 and 4.7; BSFL5: 4.7 and 4.5; BSFL10: 4.3 and 4.8; BSFL20: 4.4 and 4.4, for stomach and jejunum, respectively.
There were no significant differences in fecal DM (Table 6) or fecal score (Table 7) between the dietary treatments. BSFL5 had the overall highest fecal score and BSFL10 the overall lowest fecal scores. The graph of the overall development in the fecal score during the experimental period (Supplementary Fig. 1) shows an increased score (> 2) in the first week PW, which seems to stabilize below score two after day 12, but again increasing to score > 2 on day 23–25. Pens with wood shavings had significant lower fecal scores in week 3 (P = 0.049), and numerical lower average fecal score all experimental weeks.
Dietary treatments | P-value | ||||||
---|---|---|---|---|---|---|---|
Fecal DM, % | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Day 7 | 26.6 | 25.7 | 25.3 | 24.3 | 0.7 | 0.725 | 0.235 |
Day 14 | 27.8 | 26.1 | 26.9 | 26.3 | 0.3 | 0.212 | 0.230 |
Day 21 | 26.3 | 26.2 | 26.1 | 24.7 | 0.4 | 0.414 | 0.105 |
Day 27 | 26.7 | 26.3 | 27.3 | 26.7 | 0.4 | 0.897 | 0.872 |
Total average | 26.8 | 26.1 | 26.4 | 25.5 | 0.3 | 0.500 | 0.159 |
Dietary treatments | P-value | ||||||
---|---|---|---|---|---|---|---|
Fecal Score | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Week 1 | 1.82 | 1.87 | 1.87 | 1.84 | 0.06 | 0.992 | 0.952 |
Week 2 | 1.83 | 1.99 | 1.88 | 1.91 | 0.06 | 0.800 | 0.855 |
Week 3 | 1.73 | 1.78 | 1.56 | 1.72 | 0.06 | 0.586 | 0.817 |
Week 4 | 1.86 | 2.00 | 1.86 | 1.95 | 0.05 | 0.760 | 0.775 |
Overall period | 1.82 | 1.92 | 1.80 | 1.86 | 0.05 | 0.811 | 0.960 |
There were no significant differences in jejunal trypsin (N = 42; ANOVA P = 0.856; regression P = 0.380; SEM = 0.14) or lipase (N = 45; ANOVA P = 0.988; regression P = 0.769; SEM = 0.01) activity between the dietary treatments. The highest numerical enzymatic activities were seen in BSFL10, and the lowest in BSFL20. Numerical means for trypsin activity were 1.04, 0.90, 1.06, 0.68 mU/µg protein, and for lipase activity 0.07, 0.07, 0.08 and 0.06 mU/mg protein, for the control, BSFL5, BSFL10 and BSFL20 diet, respectively.
No significant differences between dietary treatments were found in AID of the main nutrients given by the ANOVA (Table 8). The mean AID of CP ranged from 65.9 to 69.7 among the dietary groups, but no linear relationship with increased BSFL was seen. A significant linear regression was, however, found between increased inclusion of BSFL and increasing AID of CF (R2 = 0.108) and ash (R2 = 0.106). A significant effect of sex was found for AID of CP (P = 0.026), CF (P = 0.007) and ash (P = 0.025), which all where higher for gilts compared to barrows.
Dietary treatments | P-value | |||||||
---|---|---|---|---|---|---|---|---|
AID, % | N | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Dry matter | 39 | 71.1 | 67.4 | 68.0 | 68.0 | 1.0 | 0.536 | 0.507 |
Crude protein | 39 | 69.7 | 65.9 | 66.6 | 67.6 | 1.4 | 0.728 | 0.838 |
Starch | 40 | 97.4 | 94.9 | 95.3 | 95.8 | 0.5 | 0.434 | 0.586 |
Crude fat | 40 | 81.7 | 81.7 | 84.8 | 86.8 | 1.0 | 0.491 | 0.038 |
Ash | 37 | 30.9 | 28.1 | 26.6 | 40.7 | 2.0 | 0.147 | 0.049 |
There was a significant difference in AID of tyrosine between the dietary treatments. A tendency for increased AID with increasing inclusion of BSFL was found for alanine (R2 = 0.085) and proline (R2 = 0.079), whereas AID of lysine (R2 = 0.087) and tyrosine (R2 = 0.080) tended to decrease with increased inclusion of BSFL. For the AID of the other amino acids and total amino acids, no significant differences were found between the dietary treatments (Table 9). A significant effect of sex was found for alanine (P = 0.009), arginine (P = 0.004), aspartic acid (P = 0.004), cysteine (P = 0.046), histidine (P = 0.005), isoleucine (P = 0.003), leucine (P = 0.004), lysine (P < 001), methionine (P = 0.007), phenylalanine (P = 0.005), serine (P = 0.003), threonine (P = 0.005), valine (P = 0.007) and total amino acids (P = 0.020). All higher for gilts compared to barrows.
Dietary treatments | P-value | |||||||
---|---|---|---|---|---|---|---|---|
AID of AAa, % | N | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Indispensable AA | ||||||||
Arginine | 40 | 81.1 | 80.2 | 78.4 | 78.2 | 0.8 | 0.295 | 0.164 |
Histidine | 39 | 78.3 | 76.3 | 74.1 | 75.7 | 0.6 | 0.174 | 0.190 |
Isoleucine | 40 | 77.8 | 75.7 | 73.8 | 77.5 | 0.8 | 0.476 | 0.918 |
Leucine | 40 | 79.0 | 76.8 | 75.2 | 78.1 | 0.8 | 0.476 | 0.906 |
Lysine | 40 | 86.5 | 80.2 | 78.4 | 78.2 | 0.8 | 0.234 | 0.065 |
Methionine | 40 | 87.3 | 87.6 | 85.5 | 88.4 | 0.5 | 0.252 | 0.502 |
Phenylalanine | 40 | 79.5 | 76.0 | 73.6 | 77.8 | 0.8 | 0.185 | 0.740 |
Threonine | 40 | 75.8 | 75.2 | 73.8 | 76.4 | 0.9 | 0.784 | 0.775 |
Valine | 40 | 79.2 | 76.3 | 74.6 | 77.2 | 0.9 | 0.431 | 0.628 |
Dispensable AA | ||||||||
Alanine | 40 | 69.2 | 67.3 | 68.3 | 74.0 | 1.2 | 0.330 | 0.068 |
Aspartic acid | 40 | 67.5 | 65.5 | 66.7 | 69.9 | 1.3 | 0.687 | 0.331 |
Cysteine | 40 | 63.7 | 56.4 | 60.4 | 60.0 | 1.9 | 0.558 | 0.865 |
Glutamic acid | 40 | 81.1 | 77.0 | 79.7 | 81.3 | 1.1 | 0.480 | 0.519 |
Glycine | 40 | 42.5 | 39.7 | 43.2 | 53.1 | 3.2 | 0.453 | 0.152 |
Proline | 39 | 63.5 | 66.7 | 73.2 | 74.2 | 2.2 | 0.246 | 0.083 |
Serine | 40 | 73.5 | 70.9 | 70.6 | 71.5 | 1.0 | 0.681 | 0.692 |
Tyrosine | 40 | 57.6 | 50.3 | 41.1 | 48.1 | 1.7 | 0.008 | 0.078 |
Total amino acids | 40 | 75.0 | 73.3 | 73.4 | 75.6 | 1.0 | 0.745 | 0.655 |
aDetermined using water-corrected molecular weights |
The ATTD of CP, starch, CF, phosphorus, and ash significantly differed between the dietary treatments (Table 10). The ATTD of CP decreased with increasing inclusion of BSFL (R2 = 0.056), whereas by contrast, the ATTD of CF increased with increasing inclusion of BSFL (R2 = 0.147). The starch ATTD was overall high for all diets (> 99.5%) and pigs fed the BSFL20 obtained the highest ATTD of starch. The ATTD of phosphorus increased with increasing inclusion of BSFL (R2 = 0.153). A linear increase in ATTD of ash was also found for increasing inclusion levels of BSFL (R2 = 0.067).
Gilts had significantly higher ATTD of CP (P = 0.035), starch (P = 0.015) and phosphorus (P = 0.005). Using wood shavings as bedding material in the pens gave significantly higher ATTD of CP (P = 0.035), CF (P = 0.016), ADF (P = 0.016) and ash (P = 0.030).
Dietary treatments | P-value | |||||||
---|---|---|---|---|---|---|---|---|
ATTD, % | N | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Dry matter | 79 | 83.4 | 82.8 | 83.0 | 82.7 | 0.2 | 0.743 | 0.308 |
Crude protein | 79 | 79.6 | 78.4 | 77.5 | 77.3 | 0.4 | 0.045 | 0.035 |
Starcha | 79 | 99.7 | 99.7 | 99.8 | 99.8 | 0.02 | 0.045 | 0.079 |
Crude fat | 79 | 75.7 | 78.0 | 79.4 | 80.0 | 0.4 | 0.026 | < 0.001 |
ADF | 78 | 30.3 | 27.9 | 28.4 | 29.2 | 1.1 | 0.897 | 0.858 |
NDF | 77 | 39.9 | 36.9 | 37.5 | 39.4 | 0.7 | 0.565 | 0.870 |
Phosphorus | 79 | 51.3 | 51.5 | 55.3 | 56.4 | 0.6 | < 0.001 | < 0.001 |
Ash | 79 | 58.7 | 57.4 | 57.6 | 61.2 | 0.5 | 0.043 | 0.022 |
Energy | 79 | 82.6 | 82.0 | 82.3 | 81.8 | 0.2 | 0.676 | 0.274 |
a26 feces samples had starch levels below lower detection limit |
Macroscopic evaluation of the intestinal tissue during sampling revealed some changes, mostly in jejunum and colon, with reddening of the mucosa, and low muscle tone and no mucosal folding in the jejunum (Table 11, Fig. 1). The macroscopic appearance was distributed equally with no clear differences between the pigs fed the control diet and those fed increasing levels of BSFL. However, most macroscopic changes were found in the pigs fed BSFL5.
In the jejunum, morphological changes such as lamina propria lymphocyte infiltration, lamina propria edema, submucosal edema, intraepithelial lymphocyte numbers, and enterocyte vacuolization were present in several samples. The distribution of these main morphological changes showed no notable differences between the dietary treatments.
In the ileum, tissue sections were observed to have a normal and healthy mucosal appearance.
In the colon, mild to moderate inflammatory morphological changes were observed. The inflammation was largely characterized by an increased content of lymphocytes and plasma cells in the inter-crypt compartment (Fig. 2). Other notable changes observed were mild crypt dilatation as well as a few samples showing crypt abscessation characterized by an accumulation of sloughed-off epithelial cells and neutrophils. Pigs fed the control diet were all scored “mild” or “moderate” changes for inter-crypt area lymphocyte infiltration, whereas pigs fed the highest BSFL inclusion level had the least number affected with 6 out of 12 scored “normal”. For colon crypt dilatation, some pigs were scored “mild”, but no remarkable differences between the dietary treatments were found. Dietary treatment did not influence VH, CD, or VH:CD ratio in jejunum, ileum, or colon, as shown in Table 12.
Dietary treatments | ||||
---|---|---|---|---|
Macroscopic findingsa | Control | BSFL5 | BSFL10 | BSFL20 |
Reddening of the jejunal mucosa | 3 | 3 | 1 | 2 |
Jejunum with low muscle tone and no mucosal foldings | 2 | 4 | 2 | 1 |
Reddening of the ileal mucosa | 1 | 2 | 1 | |
Reddening of the colonic mucosa | 5 | 6 | 4 | 6 |
aNumbers represent total samples observed with specific changes (N = 47). |
Dietary treatments | P-value | ||||||
---|---|---|---|---|---|---|---|
Morphometry, µm | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Jejunum | |||||||
VH | 454 | 457 | 472 | 465 | 13 | 0.976 | 0.783 |
CD | 413 | 389 | 368 | 398 | 12 | 0.717 | 0.736 |
VH:CD ratio | 1.23 | 1.25 | 1.37 | 1.30 | 0.06 | 0.916 | 0.631 |
Ileum | |||||||
VH | 232 | 233 | 250 | 247 | 7 | 0.669 | 0.334 |
CD | 134 | 148 | 150 | 137 | 4 | 0.588 | 0.954 |
VH:CD ratio | 1.92 | 1.64 | 1.75 | 1.96 | 0.6 | 0.276 | 0.420 |
Colon | |||||||
CD | 486 | 479 | 480 | 482 | 9 | 0.990 | 0.927 |
There were no significant differences in liver index (N = 44; ANOVA P = 0.557; regression P = 0.221; SEM = 0.04) between the dietary treatments. Numerical means for liver index were 2.65, 2.58, 2.64 and 2.74 for the control, BSFL5, BSFL10 and BSFL20, respectively. Barrows had higher liver index compared to gilts (P = 0.007).
There was no difference in the SCFA profile in the colon between the dietary treatments (Table 13). However, the control diet had the numerical highest level of butyric acid and the highest numerical level of total SCFA. The SCFA profile differed between sexes. Barrows had higher levels of acetic- (P = 0.035) and propionic acid (P = 0.002), but lower levels of isobutyric acid (P = 0.045), isovaleric acid (P = 0.047) and of total SCFA (P = 0.030) compared to gilts.
Dietary treatments | P-value | |||||||
---|---|---|---|---|---|---|---|---|
SCFA, µmol/g | N | Control | BSFL5 | BSFL10 | BSFL20 | SEM | ANOVA | Regression |
Acetic acid | 47 | 56.5 | 51.9 | 52.7 | 52.1 | 1.63 | 0.708 | 0.474 |
Butyric acid | 47 | 13.5 | 10.4 | 12.4 | 12.6 | 0.62 | 0.221 | 0.978 |
Isobutyric acid | 47 | 0.87 | 0.69 | 0.87 | 0.78 | 0.06 | 0.643 | 0.839 |
Isovaleric acid | 47 | 1.05 | 0.79 | 1.11 | 0.92 | 0.08 | 0.400 | 0.860 |
Propionic acid | 47 | 25.6 | 21.8 | 24.2 | 24.5 | 0.84 | 0.438 | 0.921 |
Valeric acid | 46 | 1.96 | 1.75 | 1.70 | 1.77 | 0.11 | 0.889 | 0.622 |
Total SCFA | 47 | 99.8 | 87.4 | 93.0 | 92.3 | 2.68 | 0.437 | 0.602 |
Ammonia, mM | 47 | 3.00 | 2.72 | 2.87 | 2.70 | 0.08 | 0.574 | 0.327 |
There was no difference in alpha diversity indices in the colon of piglets fed the different dietary treatments (Fig. 3a; N = 45, P = 0.355) as measured by the Shannon diversity index. Beta-diversity was assessed by several distance methods (Bray-Curtis, unweighted and weighted UniFrac). There were no differences in colon microbiota variance (beta dispersion) between the dietary treatments with neither of the distance methods. The PERMANOVA test showed significant effect of dietary treatment for the Bray-Curtis (R2 = 0.094; P = 0.029) and unweighted UniFrac (R2 = 0.082; P = 0.041) distances, but not for the weighted UniFrac distances (R2 = 0.084; P = 0.203). However, when adjusted for multiple testing, the pairwise PERMANOVA tests gave no significant differences between any of the dietary treatments. The PCoA plot using Bray-Curtis distances (Panel b in Fig. 3) shows no clear grouping of dietary treatments.
At the phylum level, the colon microbiota was dominated by Bacteroidota and Firmicutes (Fig. 3, Panel c). The relative abundance of Bacteroidota significantly differed between the dietary treatments (P = 0.026), whereas a tendency was found for the Firmicutes (P = 0.055). The Bacteroidota and Firmicutes were contributing 46.3% and 44.2% to the control group, 51.0% and 39.4% to the BSFL5 group, 50.0% and 42.1% to the BSFL10 group, and 48.1% and 38.9% to the BSFL20 group, respectively. There was also a significant effect of dietary treatment on Campilobacteria (P = 0.042) and Thermoplasmatota (P = 0.002), but the relative abundance of these two phyla were low in all groups (Campilobacteria: 1.0% in Control, 1.5% in BSFL5, 0.6% in BSFL10 and 0.5% in BSFL20; Thermoplasmatota: <0.01% in all groups).
Figure 4 shows the relative abundance of the top 10 most abundant genera in the colon. Prevotella was the most abundant genus, contributing 15.2%, 18.2%, 14.0% and 13.6% to the colon microbiome in the pigs fed control, BSFL5, BSFL10, and BSFL20, respectively. A significant effect of dietary treatment was found for the relative abundance of Lactobacillus, where the control pigs had a significantly higher relative abundance of Lactobacillus compared to the pigs fed the BSFL20 diet. A significant effect of dietary treatment was also found for the Rikenellaceae RC9 gut group, where higher relative abundance was found in the colon of pigs fed the BSFL10 diet compared to pigs fed the BSFL5 diet. No covariation was found between the colonic SCFA concentrations and the microbiota composition (i.e. the overall ASV data) (P = 0.366).
The interest of using insects in animal feed is increasing [35, 36], and inclusion of BSFL have been investigated in diets for pigs [9], chicken [12] and fish [11, 37], but to this date, there is limited information about the effect of high dietary inclusion levels. The present study demonstrated that up to 19.06% full-fat BSFL meal can be included as an alternative protein and energy source in pelleted diets for weaned piglets without affecting growth performance.
Dietary inclusion of BSFL did not alter the ADG, ADFI, or FCR for a four-week period PW. There was, however, a tendency for decreased FCR with increased inclusion of BSFL the first 14 days PW. These results are partly consistent with Yu et al. [38], which reported a linear improvement in both ADG and FCR when feeding increasing BSFL inclusion from 0 to 4% in the two first weeks PW, whereas no differences were found for a four-week feeding period. Contrary, another study showed no differences in performance when feeding up to 8% full-fat BSFL for 15 days, after weaning at 21 days of age [9].
Insects contain varying levels of chitin, depending on both types of insects and life stage [5]. This polysaccharide is the major component of the insect’s cuticle and creates a strong skeleton together with minerals and proteins. The exoskeleton of insects is different from crustaceans by having more amino acids strongly bound to the chitin [39, 40] and the degree of chitin deacetylation is dependent on the insect life stage [41]. Inclusion of chito-oligosaccharide (derivates from chitosan) in a corn-soybean meal diet has been shown to reduce the incidences and score of diarrhea in piglets weaned at 16 days of age [42], and in piglets weaned at 21 days of age and challenged with Escherichia coli [43]. However, no difference in fecal score or fecal DM were found between the dietary treatments in this study and may be due to the complexity of the chitin in the insect cuticle and the low degree of acetalization and solubility compared to the chito-oligosaccharides [44]. The increased fecal score in the first week PW is in accordance with the critical period for PW diarrhea [45]. The increase in fecal score in the last experimental week could be related to increased stress from sampling, as individual fecal samples for digestibility were collected daily in this period.
There is scarce information regarding the nutrient digestibility of insects and only a few studies including BSFL. Newton et al. [46] found similar ATTD of CP in five weeks old barrows fed 33% BSFL meal compared to a soybean meal control. In the present study, there was reduced ATTD of CP with increased inclusion of BSFL, but no difference in trypsin activity was found. A decline in ATTD of CP was also reported by Yu et al. [38] when including low amounts of BSFL (< 4%). Reduced CP digestibility with increased dietary inclusion of BSFL has also been reported for broiler chickens [47], Atlantic salmon [37, 48], and rainbow trout [49]. The CP digestibility is found to negatively correlate with chitin content in BSFL, in vitro [50]. Because of the nitrogen-rich chitin, the Kjeldahl-N method overestimates the digestible protein content of insects [51]. This is supported by the fact that the AID of total amino acids were not affected by dietary treatment. Janssen et al. [52], suggested using a conversion factor of 5.6 to avoid the overestimation of protein content in BSFL.
Limited information is available about the amino acid digestibility of BSFL in pig diets. However, the amino acid digestibility of full-fat BSF prepupae was recently investigated by Tan et al. [53]. The AID coefficients for amino acids were found to be between 0.641 and 0.821, and the SID coefficients to be between 0.767 and 1.177. During the diet optimization in the present study, the SID coefficients of all amino acids were set to 0.83. AID of lysine tended to decrease with increased inclusion of BSFL. Tan et al. [53] reported a SID coefficient of 0.776 for lysine in BSFL. Overestimation of the SID coefficient in the diet formulation of lysine might be the cause of the tendency of reduced AID of lysine with increased BSFL. Moreover, AID of alanine and proline tended to increase with increased inclusion of BSFL. The SID values for BSFL published by Tan et al. [53] were close to the one used in the diet formulation, but Tan et al. [53] also reported a higher AID coefficient for alanine in BSFL than soybean meal. Proline is high in endogenous losses, attributed to low reabsorption of mucin [54], but no information about BSFL affecting endogenous losses was found in the literature. In general, the AID of amino acids in the control diets was slightly lower compared with the control diet in a previous experiment with this age pigs [55]. Especially the AID of tyrosine was in general low and differed between the dietary treatments. The reported SID coefficient for tyrosine in BSFL (82.4%; Tan et al. 2020) is only slightly lower than the coefficient used in the diet formulation (83%). In the method used for amino acid analysis, glucosamine, resulting from hydrolyzed insect chitin, give a ninhydrin-positive peak close to the tyrosine peak in the chromatogram, which may cause a higher uncertainty for accurate tyrosine analyses. It is important to remember that even with standardized feeding time before sampling, the AID calculations with the slaughtering method are snapshots compared to methods with continuous collections or pooled collections over several days. This might contribute to why the trend in AID and ATTD are the same, but only the ATTD results are significant.
There was a linear increase in AID and ATTD of CF with increased inclusion of BSFL, whereas lipase activity measured in the oral part of jejunum showed no difference. The fat from BSFL amounted to be 19–69% of the total dietary fat and was significantly affecting the fatty-acid profile of the diets. Saturated fatty acids increased with increased BSFL inclusion, whereas the total amount of polyunsaturated acids decreased. The increase in saturated fatty acids was dominated by lauric acid (C12:0). The MCFA lauric acid is a major constituent of the BSFLs fat. It is synthesized by the larvae, and therefore present at a high level independent of diet [56, 57]. Finke [5] reported that lauric acid constituted 42% of the total fatty acids in the BSFL. The MCFAs are rapidly absorbed in the mucosa and, unlike long-chain fatty acids, directly into the portal vein for transportation to the liver [6, 58]. The differences in fatty acid composition between the diets are therefore believed to be the cause of differences in ATTD of CF. Higher digestibility of fat was also reported with 8% [9] and 33% inclusion of BSFL [46]. Contrary, Yu et al. [38] found a lower ATTD of CF when including 4% full-fat BSFL. Their reported ATTD of CF was overall low (< 70%). Different feed ingredients causing a different fatty acid profile of the diets, and low inclusion level of BSFL limiting the differences in fatty-acid profiles between diets could explain the contradictory result by Yu et al. [38].
Gilts showed higher AID and ATTD of CP, CF, and ash, as well as for AID of several amino acids, compared to barrows. Higher ATTD of CP in gilts is previously reported by [59], including a nonsignificant increase in ATTD of ether extract. The use of wood shavings as bedding material also increased the numerical ATTD except for NDF. This is probably a result of the pigs eating the wood shavings. Dietary fiber is known to decrease ATTD [60], but when not included in the dietary composition it can disturb digestibility calculations by reducing the relative amount of the nutrient in the digesta.
The MCFA can directly supply the enterocytes with energy and thereby improve morphological changes in periods with nutrient deficiency such as weaning [6]. It is well known that weaning causes villus atrophy. However, both in this study and the study by Spranghers et al. [9], no effect of BSFL inclusion was found on intestinal morphometry. The effect of MCFA on morphology could be more evident the first days PW as the VH in the small intestine seems to recover within two weeks PW [61]. Biasato et al. [10] also reported no difference in intestinal morphometric indices after feeding pigs diets with up to 10% defatted BSFL. In accordance with our results, they also observed greater morphometric indices in the jejunum than the ileum.
The gut health was mildly suboptimal due to the observation of inflammation and edema in the jejunum and colon, but the morphological disorders did not appear to have a strong association with any of the dietary treatments. The morphological changes in edema and inflammation could be related to the edema disease outbreak. The results are consistent with Biasato et al. [10] which replaced defatted BSFL with soybean meal, concluding that up to 10% BSFL inclusion does not negatively affect gut morphology or histological features.
The microbiota of the gut is involved in the nutritional, physiological, and immunological functions of the pig [62]. Variation in the diet composition is one of the most important factors affecting the GIT microbial ecosystem of the pig [63]. The microbial community in the colon was dominated by Bacteroidota and Firmicutes, also reported by Yu et al. [64]. Yu and coworkers found no differences in phyla abundances when including 4 or 8% full-fat BSFL in diets for finishing-pigs. By contrast, the relative abundance of Bacteroidota, Firmicutes, Campilobacteria, and Thermoplasmatota differed between the dietary treatments in the present study, but there was no clear dose-effect of the dietary BSFL inclusion. However, the highest relative abundance of Firmicutes and lowest relative abundance of Bacteroidota was found in the pigs fed the control diet. In accordance with Yu et al. [64], no significant differences in Shannon indices were observed.
Prevotella was the most abundant genus in the colon of the pigs regardless of dietary treatment. This is in correspondence with several authors reporting Prevotella to be the most abundant genus in the fecal microbiome in the PW period [65–67]. Lactobacillus was more abundant in the colon of the control piglets. This result is contrary to what recently was reported by Yu et al. [64], where inclusion of 4% BSFL increased the abundance of Lactobacillus. In the study by Yu et al. [64], finishing-pigs were fed corn-based diets contrary to our PW pigs, which were fed wheat and barley-based diets. Also, when they increased the inclusion to 8% BSFL, Yu et al. [64] did not find the same beneficial effects on the colonic microbiota as when feeding the 4% inclusion of BSFL. Lactobacillus is a genus of gram-positive lactic acid-producing bacteria in general considered favorable in the gut microbiome. It is known that lauric acid, abundantly found in BSFL fat, has antimicrobial effects, especially against gram-positive bacteria [6]. Spranghers et al. [9], observed inhibition of the growth of lactobacilli in an in vitro study with BSFL fat but did not observe differences in vivo in the small intestine. They also observed that a high amount of the lauric acid was absorbed already in the stomach and therefore had little opportunity to modulate the colon microbiome. In rainbow trout, it is shown that the BSFL life stage and lipid content of the BSFL are important factors influencing the gut microbiome [68].
It is important to mention the treatment with antibiotics of all pigs in the second week of the experiment. The effect of the antibiotic on the gut microbiota is antibiotic specific [62], and it is therefore difficult to know what impact the antibiotic treatment of the pig in this experiment had on their gut microbiota. Indeed, the antibiotic treatment could have affected microbiota results and be the reason for the small differences in the colon microbiome between the dietary treatments. Because of the antibiotic treatment in the second experimental week, it might be that the feeding period was too short to cause significant changes in the gut microbiome. Also, therapeutic doses of antibiotics might suppress the innate immune defenses and contribute to increased disease susceptibility [62]. Almost half of the pigs, independent of dietary treatment, had reddening of the colonic mucosa, indicating some inflammation in the colonic intestinal wall. However, the overall Shannon diversity indices reported were similar to what was reported by Yu et al. [64] in finishing pigs.
To conclude, up to 19.06% of full-fat BSFL meal could be included in a balanced diet for PW pigs without affecting performance results. Increased inclusion of BSFL decreased ATTD of CP and increased ATTD of CF but did not affect the general gut function, as assessed by enzyme activity, morphometry, and histological evaluation. Some changes in the colon microbiota composition were observed, such as decreased relative abundance of Lactobacillus with dietary inclusion of BSFL. However, there are limited information and contradictive results on the effect of BSFL on the gut microbiome in pigs, thus further investigations are needed.
ADF: Acid detergent fiber
ASV: Amplicon sequence variant
AID: Apparent ileal digestibility
ATTD: Apparent total tract digestibility
ADF: Average daily gain
ADFI: Average daily feed intake
BSFL: Black soldier fly larvae
CF: Crude fat
CP: Crude protein
CD: Crypt depth
DM: Dry matter
FCR: Feed conversion ratio
GIT: Gastrointestinal tract
MCFA: Medium-chain fatty acid
NDF: Neutral detergent fiber
PW: Post-weaning
PCoA: Principal coordinate analyses
SCFA: Short-chain fatty acid
SEM: Standard error of the mean
SID: Standardized ileal digestibility
VH: Villus height
The experiment was conducted at an animal experimental unit approved by the National Animal Research Authority (permit no. 174). The authors state that all piglets were handled under the applicable laws and regulations controlling experiments with live animals in Norway, and in accordance with the EU Directive 2010/63/EU for animal experiments.
Not applicable
The datasets analyzed in the present study are available from the corresponding author on reasonable request.
The authors declare that they have no competing interests.
This work was supported by Foods of Norway, a Center for Research-based Innovation (The Research Council of Norway; grant number 237841/030) and FeedMileage, Efficient use of Feed Resources for a Sustainable Norwegian Food Production (The Research Council of Norway; grant no. 310 233685/E50).
IMH, JØH, MØ, and LTM contributed to the conceptualization and study design. IMH, JØH, LTM, GHG, and RMÅ contributed to sample collections. RMÅ conducted enzyme activity analyses, 16s rRNA extraction, and sequencing. IMH and GHG performed statistical analysis on growth performance, liver index, fecal consistency, and digestibility data. IMH performed statistical analysis on pH, enzyme activities, short-chain fatty acids, and raw sequence data. IMH wrote the original draft. All authors critically reviewed the original draft and approved the final manuscript.
The authors would like to thank the staff at the Center for Livestock Production (SHF), NMBU, Ås, for taking care of the daily animal management during the experiment, and members of the Foods of Norway team for their help during sampling. The authors also want to thank Elvis, Amir, and Thijs for contributing to sample collection and histological evaluation of intestinal tissue.