Preparation of fermented rapeseed-seaweed feed
The FRS feed was provided by FermentationExperts (Denmark), which was a blend of rapeseed meal (Brassica napus), wheat bran (Triticum eastivum) and two types of brown seaweed (Saccharina latissima and Ascophylum nodosum) prepared via a controlled two-step solid state fermentation. The inoculum consisted of three lactic acid bacteria: Pediococcus acidilactici (DSM 16243), Pediococcus pentosaceus (DSM 12834) and Lactobacillus plantarum (DSM 12837). The addition of the inoculant controlled the process by acidifying the blend within the first 24 hours, and assuring an almost entirely anaerobic process. The process continued for 11 days at 38 ˚C. The fermented material was then dried in a spin flash dryer, with a temperature setting and pass-through-speed that preserved the viable bacteria and the microbial thermolabile metabolites.
Animal feeding and performance recording
The feeding trial was carried out on a commercial pig farm (Kawiks Farm, Patoki 23. 98–170 Widawa. Province. Lodz city, Poland) in 2018, where groups of piglets were weaned one day a week over a 5-week period. The trial procedure and sample collection were approved by the Local Ethical Commission of Olsztyn University of Life Sciences (Olsztyn, Poland) with regards to experimentation and animal care. A total of 690 piglets were tested under three different feeding regimens (230 piglets per feeding treatment) from 28 days of age (10 days before weaning) until 85 days of age when the piglets exited the nursing unit. One group was a control group fed a basal feed according to Danish nutritional recommendations[17] (0% FRS), and the other two groups received supplementation of 2.5% or 5% FRS to the basal feed (feed dry matter basis in Table 1). Piglets on each dietary regime were housed in nursing pens holding an average of 48 animals per pen. Each dietary treatment was repeated 5 times (1 repetition per experimental week and 1 pen representing a repetition) and the control was repeated 4 times. None of the diets included growth promoters, prescription antibiotics or zinc oxide. Piglets that experienced diarrhea or any other serious health conditions were removed from the experiment and treated elsewhere and counted as piglets that did not complete the experiment. Feed and fresh water were supplied ad libitum throughout the experiment. Litter weights in the nursing period were recorded every week, and feed intake was recorded daily. Performance indicators such as body weight, ADFI (average daily feed intake), average daily weight gain (ADG), feed conversion ratio (FCR) and the completion rate were calculated by pen as previously outlined[5].
Table 1
Feed formulations used for the experiment.
| Pre-starter diet | Starter diet |
| 0% FRS | 2.5% FRS | 5% FRS | 0% FRS | 2.5% FRS | 5% FRS |
Ingredients | (g/kg) |
Wheat (11.2%) | 613.26 | 597.88 | 583.86 | 539.88 | 525.11 | 508.24 |
Barley (10.6%) | 100.00 | 100.00 | 100.00 | 200.00 | 200.00 | 200.00 |
Soybean meal (46.0%) | 0.00 | 0.00 | 0.00 | 170.00 | 170.00 | 169.69 |
Digestible soy | 79.08 | 69.35 | 59.32 | 10.63 | 0.77 | 0.00 |
Fermented rapeseed-seaweed meal | 0.00 | 25.00 | 50.00 | 0.00 | 25.00 | 50.00 |
Potato protein | 40.00 | 40.00 | 40.00 | 0.00 | 0.00 | 0.00 |
Fish meal (70%) | 40.00 | 40.00 | 40.00 | 14.00 | 14.00 | 7.00 |
Whey protein | 50.00 | 50.00 | 50.00 | 0.00 | 0.00 | 0.00 |
Soy bean oil | 37.28 | 38.93 | 40.17 | 23.74 | 25.20 | 26.50 |
Limestone (Ca 38.5%) | 0.00 | 0.00 | 0.00 | 5.71 | 5.57 | 5.66 |
Calcium formate | 5.00 | 5.00 | 5.00 | 0.00 | 0.00 | 0.00 |
Calcium phosphate | 9.12 | 8.68 | 8.23 | 8.83 | 8.38 | 8.29 |
Sodium chloride | 3.83 | 3.70 | 3.58 | 4.98 | 4.85 | 4.89 |
Summer fruit | 2.00 | 2.00 | 2.00 | 2.00 | 2.00 | 2.00 |
Tretracid liquid | 5.82 | 5.83 | 5.84 | 5.60 | 5.61 | 5.75 |
Lysine HCl (98%) | 5.00 | 4.00 | 2.50 | 5.00 | 4.00 | 2.50 |
Methionine DL (99%) | 0.77 | 0.73 | 0.69 | 1.02 | 0.97 | 1.01 |
Threonine L (99%) | 1.88 | 1.83 | 1.77 | 1.94 | 1.88 | 1.85 |
Valine (98%) | 1.11 | 1.09 | 1.06 | 1.00 | 0.97 | 0.96 |
Tryptophan (99%) | 0.65 | 0.65 | 0.65 | 0.37 | 0.37 | 0.36 |
Microbial phytase | 0.15 | 0.15 | 0.15 | 0.15 | 0.15 | 0.15 |
Microbial xylanase, beta-glucanase | 0.15 | 0.15 | 0.15 | 0.15 | 0.15 | 0.15 |
Vitamin E (50%) | 0.03 | 0.03 | 0.30 | 0.10 | 0.01 | 0.01 |
Vitamin-mineral premixa | 5.00 | 5.00 | 5.00 | 5.00 | 5.00 | 5.00 |
Calculated nutritive value (%) | |
Dry weight | 88.50 | 88.50 | 88.60 | 87.50 | 87.50 | 87.60 |
Metabolizable energy (MJ) | 14.30 | 14.30 | 14.30 | 13.50 | 13.50 | 13.50 |
Crude protein | 19.70 | 19.70 | 19.70 | 18.40 | 18.40 | 18.40 |
Crude fat | 5.56 | 5.75 | 5.90 | 4.03 | 4.20 | 4.31 |
Crude fiber | 2.20 | 2.38 | 2.57 | 3.11 | 3.29 | 3.51 |
Ashes | 5.37 | 5.42 | 5.45 | 5.69 | 5.71 | 5.87 |
Starch (g) | 471.40 | 408.30 | 399.90 | 425.60 | 416.90 | 406.80 |
Lactose (g) | 36.50 | 36.50 | 36.50 | 0.00 | 0.00 | 0.00 |
Calcium | 0.83 | 0.83 | 0.84 | 0.82 | 0.82 | 0.82 |
Total phosphorous | 0.65 | 0.65 | 0.65 | 0.58 | 0.58 | 0.58 |
Digestible phosphorous | 0.59 | 0.59 | 0.59 | 0.51 | 0.51 | 0.51 |
Sodium | 0.23 | 0.23 | 0.23 | 0.22 | 0.22 | 0.22 |
Chlorine | 0.57 | 0.56 | 0.56 | 0.51 | 0.50 | 0.50 |
Potassium | 0.65 | 0.65 | 0.66 | 0.67 | 0.68 | 0.69 |
Lysine | 1.46 | 1.46 | 1.46 | 1.28 | 1.28 | 1.28 |
Methionine | 0.45 | 0.45 | 0.45 | 0.41 | 0.41 | 0.41 |
Met + Cystb | 0.77 | 0.77 | 0.78 | 0.72 | 0.73 | 0.73 |
Threonine | 0.92 | 0.92 | 0.92 | 0.81 | 0.81 | 0.81 |
Tryptophane | 0.30 | 0.31 | 0.31 | 0.26 | 0.26 | 0.26 |
Valine | 1.03 | 1.03 | 1.03 | 0.90 | 0.90 | 0.90 |
Isoleucine | 0.78 | 0.78 | 0.78 | 0.70 | 0.69 | 0.69 |
aProvided the following per kilogram of feed: Vitamin A 13,000 IU; vitamin D3 2000 IU; Vitamin E 165 mg; vitamin B1 2.5 mg; vitamin B2 7.0 mg; biotin 200 mcg; vitamin B6 4 mg; vitamin B12 50 mcg; vitamin K 3 mg; Niacin 35 mg; folic acid 1.5 mg; pantotenic acid 21.7 mg; vitamin C 100 mg; choline 0 mg; Fe 180 mg; Zn 150 mg; Cu 0 mg; Mn 55 mg; Se 0.4 mg; I 0.6 mg; Mg 0 mg. bmethionine + cysteine. |
Biological sample collection
A total of 10 piglets from each treatment (5 in each of two experimental weeks) were randomly selected and euthanized 3 weeks after weaning. The animals were euthanized by stun gunning with a captive bolt immediately followed by de-bleeding at the farm slaughtering facilities under strict sanitary regulations. Whole blood samples and serum for clinical analysis, the digesta from the colon for microbiome analysis, and jejunum and colon tissues for histopathological analyses, were collected in that order immediately after slaughtering.
A blood sample from each piglet was deposited in a tube with the anti-coagulant EDTA and preserved on ice until taken to the laboratory, where it was stored at 2–8 ˚C until analysis. Another blood sample was collected in a tube without anticoagulant, and serum separated by centrifugation, which was then stored at -20 ˚C until analysis. Gut tissues and colon contents were sampled after opening of the abdominal wall, and the stomach, small and large intestines were occluded at both ends and removed. Approximately 2 cm3 of colon content was collected from the apex of the ascending spiral of the colon with a sterile spatula and deposited in cryotubes with RNAlater™ (Sigma-Aldrich, Munich, Germany). Tubes with colon contents were kept at room temperature for less than 24 h, followed by cryopreservation in the laboratory. Tissue samples (approximate 2 cm long) of the whole transection of the jejunum and colon were excised and carefully rinsed from gut contents by flushing with saline (0.9% NaCl). For each tissue a sterilized blade was used. Tissues were preserved in 10% formaldehyde and kept at room temperature for no longer than 24 hours until further processing[18].
Blood hematology, blood biochemistry and serum immunoglobulin analysis
Full blood counts (erythrocyte, hemoglobin, hematocrit, mean corpuscular volume, mean corpuscular hemoglobin, mean corpuscular hemoglobin concentration, red cell distribution width) and differential white blood cells count (platelets, leucocytes, lymphocytes, monocytes, neutrophils, eosinophils, basophils) were performed using a Sysmex XT 2000i analyzer (Sysmex Corporation, Kobe, Japan). Serum analysis measured concentrations of the following, using standardized quantification methods: alanine aminotransferase, glutamic pyruvic transaminase, aspartate aminotransferase, glutamic-oxaloacetic transaminase, lactate dehydrogenase, lysozyme, glucose, total protein, blood urea nitrogen, uric acid, phosphorous, total cholesterol, triglycerides, low density lipoprotein, high density lipoprotein and immunoglobulin G (IgG) according to previously described procedures[5].
Histological morphometric analysis of intestinal tissues
The histological analysis of mid-jejunal and colonic tissues was conducted by a commercial analytical laboratory (ALAB Weterynaria, Warsaw, Poland) according to previous procedures[5]. In short, tissue sections fixed in 10% formaldehyde were dehydrated by means of graded ethanol and xylene baths and embedded in paraffin wax, and 3–4 µm section were then stained with haematoxylin and eosin. Histopathological evaluations (at different lens magnifications) measured gut-associated lymphoid tissue (GALT), intraepithelial lymphocytes (IELs) and lymphatic infiltration of the stromal mucosa (stromal lymphocytes, SL) counts. For GALTs, the numbers of lymphoid follicles per millimeter square were counted. For IEL scoring, the following scale was used: 0-normal (0–10 IELs/100 enterocytes), 1-low (10–15 IELs/100 enterocytes), 2-moderate (15–20 IELs/100 enterocytes; this level suggests chronic subclinical inflammation, where the intestinal-blood barrier may be damaged; weak lymphocytic inflammation), 3-severe (> 20 IELs/100 enterocytes; this level indicates chronic inflammation with infiltration damaging the epithelium and intestinal-blood barrier; moderate lymphocytic inflammation). For SL, the visual scoring scale was: 0-normal (single lymphocytes in stromal connective tissues of villus and crypts), 1-low (increased number of lymphocytes, but no damage to the stroma structures), 2-moderate (abundant infiltration of lymphocytes in stroma, damaging blood vessel walls, connective tissue fiber, reducing visibility of stroma structures), 3-severe (lymphocyte infiltration completely disrupts and conceals the stroma). In a blinded fashion, 10 fields of view per piglet at 4 × magnification were used for evaluation of GALT structures and numbers of lymphoid follicles. IEL and SL were evaluated at 40 × magnification. The analysis used a standard light microscope Olympus BX41 and Cell Sens software (Olympus Corporation, Tokyo, Japan). The gut tissues samples which could not reach the requirements for histological analysis were discarded, resulting in n = 9, 8, 10 for the 0%, 2.5% and 5% FRS group, respectively.
16S rRNA gene amplicon sequencing of colon content
Collected colon contents was stored at -60 ˚C prior to the analysis. Two types of 16S rRNA gene amplicon sequencing strategies were adopted to characterize the prokaryotic community: Illumina, NextSeq (Illumina, CA, USA) and MinION (Oxford Nanopore Technologies, Oxford, UK). The genomic DNA was extracted using Bead-Beat Micro AX Gravity Kit (A&A Biotechnology, Gdynia, Poland) according to the manufacturer’s instruction. DNA concentration and purity were measured using NanoDrop ND-1000 spectrophotometer (Saveen and Werner AB, Sweden).
Extracted DNA was diluted to 10 ng/µL prior to library preparation. The V3 hypervariable region of 16S rRNA gene was amplified and sequenced with Illumina technology as previously described[19]. Near full-length 16S rRNA gene amplicons were amplified and sequenced with ONT targeting V1-V8 hypervariable region using following primers: ONT_27Fa: GTCTCGTGGG CTCGGAGATG TGTATATAGA TCGCAGAGTT TGATYMTGGCTCAG; ONT_27Fb: GTCTCGTGGG CTCGGAGATG TGTATATAGA TCGCAGAGTT TGATCCTGGCTTAG and ONT_1540_R: GTCTCGTGGG CTCGGAGATG TGTATACTCT CTATTACGGY TACCTTGTTACGACT. Custom designed barcoding system was developed to tag encode up to 96 samples during the second round of PCR, and the PCR primer sequence is given in Table S1 (Additional file). The PCR1 reaction mix contained 5 µl of PCRBIO buffer and 0.25 µL PCRBIO HiFi polymerase (PCR Biosystems Ltd, London, United Kingdom), 1 µL of primers mix (5 µM of ONT_27Fa and ONT_27Fb, and 10 µM of ONT_1540_R, see above), 5 µL of genomic DNA (~ 10 ng/µL) and nuclease-free water to a total value of 25 µL. The PCR thermal conditions were as follows: denaturation at 95 °C for 5 min; 33 cycles of 95 °C for 20 s, 55 °C for 20 s and 72 °C for 45 s; followed by final elongation at 72 °C for 4 min.
PCR products were verified by agarose gel electrophoresis and then subjected for barcoding (PCR2). The PCR2 mix composed of 5 µL PCRBIO buffer, 0.25 µL PCRBIO HiFi polymerase (PCR Biosystems Ltd, London, United Kingdom), 2 µL of barcode primers (5 µM), 1 µL of PCR1 template and DEPC water up to 25 µL. The PCR2 thermal conditions were as follows: denaturation at 95 °C for 2 min; 13 cycles of 95 °C for 20 s, 55 °C for 20 s, 72 °C for 40 s; final elongation at 72 °C for 4 min. The final PCR products were purified using AMPure XP beads (Beckman Coulter Genomic, CA, USA) and pooled in equimolar concentrations. The pooled barcoded amplicons were subjected to 1D genomic DNA by ligation protocol (SQK-LSK109) to complete library preparation for MinION sequencing. Approximate 0.2 µg of amplicons were used for the initial step of end-prep. And 40 ng of prepared amplicon library was loaded on a R9.4.1 flow cell.
Sequencing data analysis
The raw Illumina derived dataset containing pair-ended reads with corresponding quality scores were merged and trimmed using fastq_mergepairs and fastq_filter scripts implemented in the USEARCH pipeline as described previously[19]. Purging the dataset from chimeric reads and constructing zero radius Operational Taxonomic Units (zOTUs) was conducted using the UNOISE[20]. The Greengenes (13.8) 16S rRNA gene collection was used as a reference database[21].
Data generated by MinION were collected using MinKnow software v19.06.8 (https://nanoporetech.com). The Guppy v3.2.2 basecalling toolkit was used to base call raw fast5 to fastq (https://nanoporetech.com). Porechop v0.2.2 was used for adapter trimming and sample demultiplexing (https://github.com/rrwick/Porechop). The Porechop adapter list was (adapters.py) edited accordingly and is given in Table S1 (Additional file). Sequences containing quality scores (fastq files) were quality corrected using NanoFilt (q ≥ 10; read length > 1Kb). Taxonomy assignment of quality corrected reads against Greengenes (13.8) database was conducted using uclast method implemented in parallel_assign_taxonomy_uclust.py (QIIME v1.9.1). The uclust settings were tuned on mock communities (--similarity 0.8; min_consensus_fraction 0.51) assuring annotations to the lowest taxonomic level with no false positive annotations. The settings allowed it to treat individual amplicon sequence variants as individual “seeds”. Reads classified to at least phylum level were subjected for further analysis.
Statistics
All the statistical analysis concerning phenotypic data was performed with R (v3.6.2). The difference of piglet production performance was evaluated using linear mixed model as previously outlined[5] and orthogonal polynomial contrast was used to appreciate the effect on increasing dose of the FRS (0%, 2.5%, 5%). The blood hematology and biochemistry data was analyzed by R package compareGroups[22] (v4.0) using “comparaGroups” command and the descriptive table was generated by “createTable” command. In R package compareGroups, the significant differences among groups were determined by anova and Tukey’s procedure for post hoc tests. Wilcoxon rank-sum test was used to evaluate the histological difference between groups.
For microbiome analysis, QIIME 2[23] (v2018.11) combined with R packages (ggplot2, vegan, corrplot, Rhea, rstatix, vennDiagram) were used. Three samples were removed due to inadequate library size (< 1000 counts), resulting in n = 9, 8, 10 for 0%, 2.5% and 5% FRS group, respectively. For both sequencing strategies, all the samples were summarized at the L7 levels (species) and rarefied to the same sequencing depth (11000 reads/sample) for alpha and beta diversity calculations. Rarefaction on the zOTU table (Illumina data) was adopted as comparison for rarefaction on the species-level summarized table. Principal coordinate analysis (PCoA) plots were generated using binary Jaccard and Bray Curtis distance metrics, and PERMANOVA was performed to determine differences between groups and p values were adjusted by Benjamini-Hochberg correction. ANCOM[24] was adopted to identify differentially abundant taxa between groups at summarized L7 level. For taxa identified by ANCOM, Wilcoxon rank-sum test was adopted for pairwise comparison. Phenotypic data were integrated with species-level bacterial abundances by Pearson’s correlation analysis using R package Rhea[25]. Rare microbial features were removed with a cutoff of mean relative abundance > 0.1% and minimal presence among 30% of samples. Zeros were regarded as NA. Centered log-ratio transformation was conducted in both the microbial relative abundance and phenotypic data.