Skip to main content

Effect of Cyberlindnera jadinii yeast as a protein source on intestinal microbiota and butyrate levels in post-weaning piglets

Abstract

Background

Dietary yeast inclusions in a pig diet may drive changes both in gut bacterial composition and bacterial functional profile. This study investigated the effect of Cyberlindnera jadinii as a protein to replace 40% of the conventional proteins in a diet for weanling pigs on the microbiota in the small and large intestine, colonic short-chain fatty acid concentration, and colonic histopathology parameters. Seventy-two pigs weaned at 28 days of age were randomly assigned to either a control or a C. jadinii-based diet and followed for 2 weeks.

Results

Compared with the controls, higher numbers of cultivable lactic acid-producing bacteria in the small and large intestine were registered in the yeast group. Alpha and beta bacterial diversity were different between the diet groups with lower alpha-diversity and distinct bacterial composition in the large intestine in the yeast group compared with those of the controls. The large intestine microbiota in the yeast group had higher numbers of Prevotella, Mitsuokella and Selenomonas compared with those of the controls. The concentrations of colonic acetate and butyrate were higher in the controls compared with that of the yeast group. The colonic crypt depth was deeper in the control group. The gut histopathology of colonic tissues revealed no differences between the diets.

The colonic crypt depth tended to be deeper with higher relative abundance of an unclassified Spirochetes, higher colonic butyrate concentration, and higher bacterial richness. The concentration of colonic butyrate was positively associated with the relative abundance of the Faecalibacterium prausnitzii, Dialister, and an unclassified amplicon of the Spirochaetaceae family in the colon.

Conclusions

The replacement of the conventional proteins by proteins from Cyberlindnera jadinii in a weanling pig diet reshaped the large intestine microbiota structure. The novel yeast diet appeared to be selective for Lactobacillus spp., which may represent an added value resulting from using the sustainably produced yeast protein ingredient as an alternative to conventional protein ingredients in animal diets. The large intestine bacterial composition and their metabolites may be involved in an adaptive alteration of the colonic crypts without pathological consequences.

Introduction

Up to 70% of globally produced soy is used to maintain livestock production [1]. Sustainable protein alternatives are needed to reduce the dependency on soy and other conventional proteins as ingredients in the feed for animal husbandry. As early as in the 1940’s, researchers pursued the idea of replacing a substantial fraction of protein in animal feed with proteins from yeast derivatives [2, 3]. It soon became apparent that the lack of cost efficient methods for large scale yeast production would limit the use of yeast proteins [3]. Moreover, additional costs arose from the necessity of vitamin D and calcium supplementation to counteract the rachitogenic effect of yeast diets [4, 5]. Today, technology exists to produce yeast by industrial fermentation of Picea abies second-generation sugars as a carbon and energy source [6]. Our previous work has shown that up to 40% of conventional protein in a pig diet can be successfully replaced by proteins from a strain Cyberlindnera jadinii yeast [7]. The addition of yeast as a protein source supported high growth performance and improved gut health of the weanling pigs [7]. To date, a number of studies have investigated the effect of yeast supplementation on pig microbiota composition. Addition of live Saccharomyces cerevisiae to a diet promoted overgrowth of Mitsuokella bacterial genus in the large intestine microbiota in weanling piglets [8]. Inclusion of cider yeast probiotic to a diet shifted faecal microbiota towards higher numbers of Selenomonas and Prevotella in weanling piglets [9]. In addition, Upadrasta and co-workers reported reduction in Faecalibacterium, Roseburia, and Eubacterium in faeces of the yeast group. Impact of yeast derived components such as cell wall β-glucans and mannan-oligosaccharides on the gut microbiota in pigs has also been studied. Fouhse and co-workers reported high relative abundance of Mitsuokella and low relative abundance of Coprococcus and Roseburia in caecum of piglets supplemented with yeast derived mannan-rich fraction [10]. Nakashimada et al. studied changes in pig faecal bacterial composition using an in vitro intestinal model. These investigators found lower numbers of Faecalibacterium in the reactor system with addition of yeast cell wall components than without [11]. While supplementation of yeast ingredients does seem to promote distinct intestinal bacterial groups, the reduction in short-chain fatty acid (SCFA) producing bacteria may be another intrinsic feature of such diets. One of the major SCFAs produced by intestinal bacteria, butyrate, is an exogenous metabolite with a number of key functions related to gut homeostasis (reviewed in [12, 13]). While serving as fuel for colonocytes [14, 15], it is debatable whether high molarities of butyric acid are beneficial (reviewed in [16]). For instance, high colonic butyrate concentration is believed to modulate colonic crypt architecture [17], induce apoptosis in the stem cell compartment of crypts [18], and supress crypt stem cell proliferation [19]. Recent publications have been primarily focused on yeast as a feed additive and have investigated the effects of low levels of inclusion of yeast and its components on the gut microbiota in pigs [8, 9, 20]. However, little is known how inclusion of high levels of yeast affects microbial community of intestines. We hypothesize that the novel yeast diet can reshape intestinal microbiota composition in weanling piglets. The reason for featuring the post-weaning period in this study was because of the stress the animals experience during that period [21, 22], which may define the course of animal health development.

We used 16S rRNA bacterial gene sequencing and cultivation methods to compare the gut microbial consortia of yeast fed weanling piglets with that of the controls. Also, we investigated a possible role of individual bacterial groups in relation to the large intestine butyrate production and utilization.

Methods

Animals, housing, diet allocation

The trial was conducted at an experimental farm of the Norwegian University of Life Sciences (NMBU), Ås, Norway in the fall of 2017. A total of 72 crossbred [(Norwegian Landrace x Yorkshire z-line) x (Duroc) and (Norwegian Landrace) × (Duroc)] weanling piglets, selected from ten litters, was included in the experiment. The piglets were selected to enter the study based on their weight at the day of weaning, and after blocking by litter and body weight, the pigs were randomly allocated to either the control or the yeast diet. All animals were healthy during the nursery period and throughout the experiment. The animals were housed in environmentally controlled pens with a slatted floor at front and roofed resting area with a rubber mat. The animals were introduced to creep feed 2 weeks prior to weaning. The experiment was initiated when the piglets were weaned at 28 days of age (day 0 PW). Five to six piglets were grouped together in each pen and group-fed one of the allocated diets. All animals had ad libitum feeding and access to drinking water throughout the experiment. Diets were formulated to be isonitrogenous and isoenergetic based on the chemical composition of the ingredients and to meet, or exceed, the nutrient requirements of weanling pigs (Table 1). In the yeast diet, 40% of the crude protein derived from Cyberlindnera jadinii cells (LYCC 7549; Lallemand Yeast Culture Collection). The yeast cells were processed as described previously [6]. Briefly, after fermentation, the cells were washed, centrifuged, heat-inactivated, and dried. The diets were cereal-based (wheat, barley and oats), and the main protein ingredients in the control diet (soybean meal, potato protein concentrate, fish meal, and rapeseed meal) were partly replaced by yeast meal in the yeast diet (Table 1). At the days 2, 4, 7, and 14 PW, eight animals from each of the two feeding groups were sacrificed followed by sampling (Fig. 1). In addition, eight littermates were sampled at day zero to provide a baseline point for the day of weaning.

Table 1 Ingredients (g/kg as fed) and analysed chemical composition (g/kg DM, unless otherwise stated) of experimental diets
Fig. 1
figure 1

The outline of the study sampling design. The timeline is shown as circles connected by the dashed line from left to right. It starts from day 0 PW (d0 corresponds to 28th postnatal day) and continues till day 14 PW (d14). The sampling procedures are given above (16S rRNA gene sequencing, histology, and short-chain fatty acid analysis) and below (cultivation) the timeline

Bacterial cultivation / DNA extraction / 16S rRNA gene amplicon sequencing

Luminal contents from ileum distalis, apex ceci, and apex coli spiralis were collected. Serial dilutions of 0.1 mg/ml of digesta in 0.9% saline were inoculated onto media. MacConkey, Tryptose Sulfite Cycloserine (TSC), de Man, Rogosa and Sharpe (MRS), and Slanetz and Bartley agar (Oxoid, Cambridge, UK) were used to recover and quantify coliforms, Clostridium perfringens, lactic acid bacteria (LAB), and enterococci, respectively. The dilution and incubation schemes were applied as described previously [23].

For 16S rRNA gene sequencing, digesta samples from ileum distalis, apex ceci, and apex coli spiralis were snap-frozen in liquid nitrogen and stored at -80 °C until DNA extraction. The DNA extraction was carried out on samples collected on days 0, 7 and 14 PW, according to a previously described protocol [24] with minor modifications. Briefly, 200 mg of thawed gut contents were added to 1 ml of InhibitEX Buffer (QIAGEN, GmbH, Hilden, Germany) following loading 500 mg of zirconia/silica beads ( = 0.1 mm, Carl Roth, Karlsruhe, Germany). The TissueLyser adaptors were cooled down at − 20 °C for 15 min prior to the bead-beating step. The bead-beating lasted for 1.5 min at 30 Hz in TissueLyser II (Qiagen, Retsch GmbH, Hannover, Germany). Proteins were digested with 30 μL of Proteinase K II (QIAGEN, GmbH, Hilden, Germany). DNA was bound to QIAamp spin column followed by washing with AW1 and AW2 buffers (QIAGEN, GmbH, Hilden, Germany). DNA was eluted with ATE buffer (QIAGEN, GmbH, Hilden, Germany). The yielded DNA purity was assessed by NanoDrop instrument (Thermo Fisher Scientific, Waltham, MA) with subsequent quantification by Qubit fluorometric broad range assay (Invitrogen, Eugene, OR, USA).

The library preparation and amplicon sequencing of V1-V3 hypervariable region of bacterial 16S rRNA gene were performed at GATC Biotech AG (Konstanz, Germany) using 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 534R (5′-ATTACCGCGGCTGCTGG-3′) primers. A microbiome standard, 20 Strain Staggered Mix Genomic Material (ATCC® MSA-1003™), was used as a positive control. The amplicon sequencing was run in three batches on an Illumina HiSeq 4000 sequencer. The resulting sequences were deposited in the SRA (PRJNA580284).

Illumina 16S rRNA gene amplicon data curation

Forward Illumina demultiplexed reads were taken to the analysis. The reads were analysed using DADA2 R package, version 1.8.0 [25]. The core DADA2 algorithm applied had the following setup: A) Quality filtering parameters: maxEE = 1, truncQ = 2 with forward primer clipping; B) Dereplication and denoising of the quality controlled reads; C) Resulting feature tables obtained from separate Illumina runs were merged with subsequent chimaera removal; D) Taxonomic assignment using RDP Naive Bayesian Classifier implemented in DADA2 R package (default settings). The GreenGenes database, version 13.8, [26] was used as a reference database for taxonomy assignment. The LULU post-clustering algorithm was applied to optimize diversity metrics [27].

Colonic SCFA, growth performance and liver index measurement

The colon digesta samples at day 7 PW and day 14 PW were the replicate samples of those used for the 16S rRNA gene sequencing analysis. The samples were thawed on ice. A mixture of 500 mg of gut contents sample and 500 μl of ice cold ddH2O was sonicated for 5 min in cold water. Next, after mixing and centrifugation (15 min, 4 °C, 15000 g), the supernatant was transferred to a spin column (45 kDa). After another centrifugation step (15 min, 4 °C, 15000 g), the samples were spiked with internal standards. Short-chain fatty acids (SCFA) were measured by “TRACE 1300 Gas Chromatograph” with autosampler, “AS 1310” (Thermo Fischer Scientific, Milan, Italy). Parameters of the capillary column were as following: model – “Stabilwax – DA”; length – 30 m; inner diameter – 0.25 mm; film thickness – 0.25 μm (Restek corporation, Bellefonte, PA, USA). The column operating protocol was as follows: starting temperature – 90 °C (2 min); temperature increase - 10 °C/min until 150 °C, 50 °C/min until 250 °C (1 min). The rate of Helium flow was 3 mL/min. Concentrations of acetic, propionic acids, as well as butyric, valeric acids and their isomers, were reported in μmol per gram of intestinal contents. The average daily gain (ADG) in this study was calculated as: (slaughter day body weight – body weight at weaning) / number of days PW. The liver index was calculated as used previously in [7]: liver index = liver weight (kg)/ live body weight (kg).

Histology

Colon tissue samples were collected within 20 min of euthanasia and fixed in 10% formalin. The gut contents were emptied, and mucosal surface was rinsed gently with cold water prior to formalin fixation. After 48 h of fixation, the tissues were routinely processed, embedded in paraffin and 4 μm sections were mounted on glass slides. The sections were subsequently deparaffinized in xylene and rehydrated in graded alcohol before routine staining with haematoxylin and eosin. The colonic tissues were evaluated histopathologically and scored semi-quantitatively where no pathology was scored 0, very mild tissue changes received the score 1, mild changes 2, moderate changes 3, and severe changes 4. Formalin-fixed, paraffin-embedded tissue sections were also stained with high iron diamine and alcian blue (HID-AB). Digital images of the intestinal sections were captured using NanoZoomer (Hamamatsu Photonics). Morphometric measurements were performed using the software Aperio Image Scope v12.3.3.5048 (Copyright Leica Biosystems Pathology Imaging, 2003–2016). For crypt depth (CD) measurements, the ten longest and well oriented crypts were selected, and micrographs were captured at 10× magnification. CD was measured from the crypt opening at the mucosal surface to the deepest portion of the crypt adjacent to the tunica muscularis mucosae.

Statistical analysis

The sample-size estimation was based on our pilot study with the similar design and feed composition (unpublished). To compare bacterial average CFUs recovered from selective plates, the non-parametric Mann-Whitney-Wilcoxon (MWW) test was applied. The linear regression model was used to predict variance in LAB colony-forming unit (CFU) with the diet as an explanatory variable.

The Shannon and Observed species alpha-diversity indices were calculated separately on the data at the ASV level and species level. To bin the ASVs to the species level, tax_glom() R function was applied [28]. Comparison of the resulting alpha-diversity figures between the diet groups was done using MWW test. The beta-diversity analysis was performed via principle coordinate analysis (PCoA) on Bray-Curtis dissimilarity matrix, and permutational multivariate analysis of variance (PERMANOVA) test for covariate significance using adonis() R function [29], 9999 permutations. The covariates included in the statistical model were the following: diet, sex, pen, and sow.

To screen for bacteria that appeared in higher numbers in one of the feeding groups compared with those of another group, or differentially abundant taxa, the analysis of composition of microbiomes (ANCOM) test was used (false discovery rate (FDR) = 0.05, multiple correction = 2). The test was performed at the phylum, family, and ASV levels [30].

To compare average concentrations of SCFA and to compare colonic CD in the colon between the diet groups, the MWW test was applied.

The loglinear analysis was applied to the histopathology parameter results for comparison between the diet groups. Multiple regression analysis was used to predict the colonic CD. The amplicon sequence variant (ASV) table was transformed to the relative abundances to derive individual bacterial ASV relative abundance figures.

To explore correlations between bacterial group relative abundance and metadata variables, the Pearson’s correlation coefficient (reported as r) was applied. To aid the graphical representation of the multiple regression modelling, the numeric variables of the model equation were subjected to principle component analysis (PCA) in R using prcomp() function to be further displayed on a biplot. The statistical significance was declared at p-values < 0.05 for all tests.

Results

Cultivation results

All the piglets appeared healthy throughout the experiment. There was no mortality, and no difference in feed intake and growth rate between the dietary treatments.

Lactobacillus spp.

In the jejunum, LAB were found in higher numbers on average in the yeast group (9.57 logCFU/g) compared with the controls (7.30 logCFU/g) at day 4 PW (p < 0.001) (Additional file 1). The same pattern was observed in the ileum at days 4 and 7 PW (yeast = 9.48 vs control = 8.44 logCFU/g, and yeast = 10.0 vs control = 8.61 logCFU/g, respectively) (p < 0.001 for both tests). The variance in the LAB counts in the ileum at days 4, 7 and 14 PW was explained by diet (R2adj = 0.45, p < 0.05). Interestingly, when only day 7 was considered, the same linear model could explain 65% of the diet-related variance in the LAB counts. In the cecum and colon, LAB counts were also higher in the yeast group at day 7 PW, and day 14 PW (p < 0.05 for both), except in the colon at day 14 PW.

Enterococcus spp.

The counts of enterococci (8.98 logCFU/g) were found to be higher in the ileum of the yeast group at day 4 PW compared with those of the controls (8.09 logCFU/g) (p < 0.001). At day 7 PW, enterococci in the colon of the yeast group were higher than those of the control group (9.03 logCFU/g vs 8.20 logCFU/g) (p < 0.001).

Coliforms

Coliforms were at higher numbers (9.72 logCFU/g) in the cecum of the yeast group at day 7 PW compared with those of the controls (8.47 logCFU/g) (p < 0.001).

C. perfringens

No statistically significant difference was observed between the two feeding groups.

16S rRNA gene sequencing results

There were on average 449,177 (SD = 57,148) reads per sample after filtering, denoising, and chimeric amplicon removal. There were 2100, 3301, and 3485 ASVs detected in the ileum, cecum, and colon samples, respectively. At day 0 PW, there were 2645 ASVs identified for all sampled gut locations. Similarly, 3378 ASVs and 2994 ASVs were found at day 7 PW and day 14 PW, respectively. The results of sequencing of the positive controls (mock communities) are given in Additional file 2.

Alpha microbial diversity

The microbial communities in the large intestine of the yeast group were less diverse in comparison to those of the control pigs at day 7 and 14 PW (Fig. 2a, Additional file 3: Table S3) as measured by the Shannon diversity index at the ASV level. There were more distinct ASVs identified in the cecum of the controls than those of the yeast group at day 7 PW and day 14 PW (p < 0.05 for both). There was no difference in alpha microbial diversity when comparing ileum microbiotas between the two diets (Fig. 2a, Additional file 3). Interestingly, when compared at the species level, the Shannon diversity index was higher in cecum microbiota of the controls at day 7 PW only. Otherwise, alpha diversity analysis at the species level showed no difference between the diets (Additional file 3).

Fig. 2
figure 2

Gut microbiota diversity and composition in postweaning pigs fed with control and yeast diets. Panel a: Shannon diversity index in the pig gut microbiota denoted by day, gut site, and diet (0 day 0 PW, 7 day 7 PW, 14 day 14 PW; IL ileum, CE caecum, CO colon; B baseline diet, C control diet, Y yeast diet) and coloured by diet (baseline orange; control pink; yeast dark cyan). The p-values derive from MWW test comparing the averages of the Shannon diversity index between the diets at each sampling day (n = 8). The box size corresponds to IQRs with the median value represented as the lines inside the box. The whiskers represent upper and lower quartiles of the diversity estimates. Panel b: Principal coordinate analysis plot of pig gut microbiota coloured by diet and gut segment (B.CE baseline caecum, orange; B.CO baseline colon, dark orange; B.IL baseline ileum, yellow; C.CE controls caecum, dark pink; C.CO controls colon, indian red; C.IL controls ileum, light pink; Y.CE yeast caecum, dark cyan; Y.CO yeast colon, teal; Y.IL yeast ileum, cyan) and shaped by the day PW ( day 0; ▲ day 7; ■ day 14). Panel c: Stacked bar plot showing group average relative abundance of six top abundant bacterial populations at the phylum level in the pig gut denoted by day, gut site, and diet (0 day 0 PW, 7 day 7 PW, 14 day 14 PW; IL ileum, CE caecum, CO colon; B baseline diet, C control diet, Y yeast diet) and coloured by diet (baseline orange; control pink; yeast dark cyan). The x-axis shows the relative proportions of the bacterial groups coloured with distinct colours. Panel d: Stacked bar plot showing group average relative abundance of ten top abundant bacterial populations at the family level in the pig gut denoted as in the panel c

Beta microbial diversity

The PCoA was conducted to compare the gut microbial compositions of individuals by visualizing Bray-Curtis dissimilarity matrices on the plot followed by PERMANOVA statistical test for significance of study covariates. The large intestine microbiotas tended to cluster together according to the diet type (Fig. 2b). At day 7 PW, the diet accounted for 23% of the variance in the large intestine microbial composition (p < 0.001). Even more of the explained variance in microbiota (26%) was attributed to the diet when the large intestine data at day 14 PW were analysed (p < 0.001). Notably, 12% of the variance in the large intestine microbiota was attributed to the sex of the animals at day 7 PW (p < 0.05).

The variance in the ileal microbiota composition was not explained by the dietary treatment, nor due to any other tested covariates (i.e. sex, pen, and sow).

Distribution of major taxa and differentially abundant taxa

Ileum

The ileal microbial consortia primarily consisted of Firmicutes (78%), Proteobacteria (9%), Tenericutes (7%), and Bacteroidetes (1%) on average irrespective of the diet (Fig. 2c).

Phylum level

There were no differentially abundant phyla identified when comparing the diet groups at any of the sampling days.

Family level

The Porphyromonadaceae family was more abundant in the control group compared with that of the yeast group (Additional file 4B).

ASV level

Clostridiaceae 02d06 ASV and family S24 unclassified ASV of Bacteroidales order were found in higher numbers in the control group at day 7 PW compared with those of the yeast group (Additional file 4A). While Prevotella ASV9, Prevotella copri ASV23 and ASV33, and unclassified ASV of Lactobacillus genus were overrepresented in the control ileal microbiota day 14 PW, Clostridium perfringens ASV2 and ASV7, and Lactobacillus salivarius ASV5 were differentially abundant in the yeast group (Additional file 4B).

Caecum

The most abundant phyla in the cecum were the following: Bacteroidetes (64%), Firmicutes (30%), Proteobacteria (3%), and Cyanobacteria (1%) (Fig. 2c).

Phylum level

There was a higher proportion of Bacteroidetes in the yeast group at day 7 PW compared with those of the control group (Fig. 3c, Additional file 4D). Low abundant Spirochaetes phylum was overrepresented in the control group at day 14 PW compared with that of the yeast group (Additional file 4E).

Fig. 3
figure 3

Gut health parameters in postweaning pigs fed with control and yeast diets. Panel a Boxplot showing comparison of colonic crypt depth (CD) measurements between the control and yeast group at day 7 PW and day 14 PW coloured and denoted by the diet and day (7.C controls day 7 PW pink, 7.Y yeast day 7 PW dark cyan, 14.C controls day 14 PW pink,14.Y yeast day 14 PW dark cyan). MWW test p-values are provided above the boxes. Panel b Representative section of colon mucosal crypts from the control group at day 14 PW (HID-AB stain, scale bar 250 μm). Panel c Representative section of colon mucosal crypts from the yeast group at day 14 PW (HID-AB stain, scale bar 250 μm). Panel d: Principal coordinate analysis illustrating contribution of Spirochetes relative abundance, colonic butyrate concentration, and 16S rRNA gene richness to the colonic crypt depth. The individual observations and the correlation circles are coloured by the diet (controls pink, yeast dark cyan). The adjusted R2 is given at the right top corner. Panel e: Histopathological assessment of the colon mucosa, evaluating epithelial damage (ED), number of intraepithelial lymphocytes (IELs), presence of crypt abscess (CA), and infiltration of leukocytes neutrophils (N), macrophages (M), lymphocytes (L), eosinophils (E), in addition the number of piglets that were diagnosed with a very mild or mild colitis (DI). The horizontal stacked bar plot shows the number of animals with none (0), very mild (1), mild (2), or moderate histopathological changes

Family level

An unclassified family of the Tremblayales order was found in higher numbers in the control group day 7 when compared with that of the yeast group (Additional file 5C). At day 14 PW, the family s24.7 of Bacteroides was more predominant in the control group, whereas an unclassified family of Alphaproteobacteria phylum (RF32 order) was overrepresented in the yeast group (Additional file 5D).

ASV level

At day 7 PW, three Prevotella affiliated amplicons (ASV28, ASV33, ASV50), as well as Mitsuokella genus ASV17, were differentially abundant in the yeast group when compared with those of the control group. Conversely, Faecalibacterium prausnitzii ASV35 was found to be more abundant in the control group at the same time point (Additional file 5A). Seven different Veillonellaceae family ASVs, including Selenomonas ruminantium, Bulleidia p.1630.c5, Parabacteroides, and four other unclassified taxa were overrepresented in the yeast group at day 14 PW compared with those of the control group. Three variants of unclassified Prevotella ASV, Selenomonas, Mitsuokella, Mucispirillum schaedleri, and unclassified ASV of the Porphyromonadaceae family were overrepresented in the yeast group at day 14 PW when compared with those of the control group. An ASV classified as Selenomonas ruminantium and Prevotella, Lactobacillus, Campylobacter, and Anaerovibrio ASV33 genera were differentially abundant in the control group at day 14 PW in comparison with those of the yeast group (Additional file 5B).

Colon

The colon microbiota structure resembled that of the cecum with Bacteroidetes (73%), Firmicutes (19%), Proteobacteria (3%), Spirochaetes (1.8%) and Cyanobacteria (1.6%) representing the most dominant phyla (Fig. 2c).

Phylum level

There was no statistically significant difference in bacterial phyla abundances between the diet groups at any of the sampling days.

Family level

At day 7 PW, more Victivallaceae family ASVs were detected in the control group than in the yeast group (Additional file 6C). At day 14 PW, Succinivibrionaceae and Bacteroidales p.2534.18B5 families were found in higher numbers in the control group, whereas Oxalobacteraceae bacterial family was more abundant in the yeast group (Additional file 6D).

ASV level

One Mitsuokella amplicon variant was more abundant in the yeast group, while an unclassified species of the Bacteroidales order (p.2534.18B5 family) was less abundant in the same group at day 7 PW. Two distinct Prevotella amplicon variants (ASV2 and ASV50) were differentially abundant in, respectively, control and yeast group. Notably, ASV50, also previously identified in the cecum of the same time point, ranked first on relative abundance, when all samples were considered (Additional file 6A). At day 14 PW, four Prevotella ASVs (ASV17, 50, 35, 41), two distinct ASVs classified as Selenomonas ruminantium (ASV10, 2), Mitsuokella ASV17, Parabacteroides ASV22 and ASV23, Bulleidia p.1630.c5 and an unclassified ASV of the Veillonellaceae family were found in higher amounts in the yeast group than those in the control group. Four different Prevotella amplicons (ASV18, ASV9, 67, 96), Selenomonas ruminantium ASV9, Anaerovibrio ASV16, and two unclassified ASVs of Bacteroidales and YS2 bacterial orders were more abundant in the control group than in the yeast group (Additional file 6B).

SCFA in the colon

The total colonic SCFA concentration did not differ between the diet groups when measured at day 7 PW (p = 0.32). However, at day 14 PW, the levels of total SCFA tended to be higher in the control group compared with those of the yeast group (p = 0.065). For day 14 PW, butyrate and acetate were found at higher concentration in the control group compared with those of the yeast group (p < 0.05) (Table 2). The concentrations of propionate, valerate, iso-butyrate, and iso-valerate did not differ between the two groups at the statistically significant level.

Table 2 Colonic SCFA concentration comparison between the control and the yeast group. The concentration values are reported as μmol/g of colon digesta. Comparison pairs that correspond to p-values less than 0.05 (MWW test) are given in bold

Histology

The colonic CD in the control group was on average deeper than that of the yeast group at both day 7 PW and day 14 PW (p-value = 0.007 and p-value = 0.01, respectively) (Fig. 3a, b, c). The colonic butyrate concentration positively correlated with the crypt depth irrespective of the diet (r = 0.55, p-value = 0.001). The prediction model of the colonic CD showed that the depth tended to be deeper with higher numbers of an unclassified Spirochaetaceae ASV (β = 25,500, SE = 5920, t-value = 4.3, p = 0.0002), higher colonic butyrate concentration (β = 6.33, SE = 2.11, t-value = 2.9, p = 0.006), and richer colon microbiota as calculated at ASV level (β = 0.198, SE = 0.094, t-value = 2.1, p = 0.044). Overall the statistical model could be predictive of 65% of variance in colonic CD (Fig. 3d). No statistically significant difference was found in histopathological parameters in the colon comparing the two feeding groups (Fig. 3e).

Colonic bacteria – colonic butyrate association

The concentration of colonic butyrate positively correlated with the relative abundance of F. prausnitzii in the colon (r = 0.73, p < 0.0001) (Additional file 7A). The latter accounted for 52% variance in butyrate concentration as estimated by the linear model equation. Inclusion of an unclassified Spirochaetaceae family ASV and Dialister genus relative abundance to the existing linear model improved the model with 72% of variance in colonic butyrate production explained Additional file 7B). Notably, the Oxalobacter genus, a member of the differentially abundant Oxalobacteriaceae family in the yeast group, negatively correlated with the colon butyrate concentration (r = − 0.71, p = 0.002).

Butyrate – liver – colonic bacteria association

The linear model for predicting the liver index from individual bacterial groups, revealed that F. prausnitzii ASV and an unclassified Spirochaetaceae family ASV could explain 47% of variance in the liver index (β = 43.4, SE = 13.3, t-value = 3.2, p = 0.006 and β = 49.7, SE = 22.0, t-value = 2.26, p = 0.04, respectively) (Additional file 8A).

The colon concentration of butyrate positively correlated with the liver index when the two sampling days, day 7 and 14 PW, were considered (r = 0.65, p < 0.0001). However, when stratified by diet, the strength of association was different for the control group (r = 0.9, p < 0.0001) and the yeast group (r = 0.66, p < 0.01) (Additional file 8B, C).

ADG-butyrate association

At day 7 PW, the relative abundance of the Proteobacteria phylum in colon was negatively correlated with ADG (r = − 0.65, p = 0.009). The same trend but of lesser magnitude was observed in the ileum and colon samples at both day 7 PW, and day 14 PW (r = − 0.37, p = 0.04 and r = − 0.41, p = 0.02, respectively). Notably, the relative abundance of the Succinivibrionaceae family was positively correlated with ADG (r = 0.47, p = 0.006). A positive correlation was found between the relative abundance of the Prevotellaceae family in the colon of the yeast group and ADG (r = 0.53, p = 0.04).

Discussion

We investigated the effect of a high level Cyberlindnera jadinii yeast diet on the gut bacterial compositions in weanling piglets. Protein from C. jadinii yeast was used to replace 40% of crude protein in a conventional Norwegian piglet diet. The growth performance and histopathology analysis of the gut tissues indicated that the data obtained was from equally healthy animals. To explore the bacterial composition, 16S rRNA gene sequencing and bacterial cultivation were used. The sensitivity of the V1-V3 16S rRNA gene amplification assay varied in relation to different bacterial groups. For instance, there was low sensitivity to Escherichia coli detection as revealed by the sequencing of a known community standard. With help of cultivation techniques, a fair comparison of coliform numbers between the control and yeast group was obtained. To compare the microbial composition between the diet groups, a phylogeny-agnostic method, permutational multivariate analysis of variance on Bray-Curtis dissimilarity matrix, was used. This decision was based on our observations that the single end sequencing data had provided a limited phylogenetic signal for phylogeny-informed beta diversity method estimates. In contrast, the use of Bray-Curtis dissimilarity matrix was discriminative enough for the comparison of the microbiotas that were presumably similar in their compositions. As expected, the microbiota of the ileum was structurally different from that of the large intestine. Conversely, there was similarity of microbiota compositions between the cecum and colon gut segments. The diets did shape the composition of the large intestine microbiota. The ileum microbiota composition, in contrast, was less affected by diet with only few differences in microbiota at the family and ASV level. This may due to several factors that the host exerts on the bacterial succession in the small intestine. These factors, such as peristalsis, bile acids, pancreatic enzymes, hydrogen ion concentration, and local immunity seem to limit bacterial colonization of the small intestine to those bacterial species recognized by the host immune system [31, 32].

We detected a difference in the alpha bacterial diversity of the large intestine between the yeast and the control groups based on the Shannon index. The control group appeared more diverse at the amplicon level. However, when the bacterial richness and bacterial diversity were examined at the species level, the difference between the diets became less apparent. Only in cecum samples day 7 PW did the diversity figures differ at a statistically significant level. To explain this discrepancy between the analyses done at different resolution levels, namely, ASV and species level, we adopted the following logic; the calculation of bacterial community richness is affected by the number of ASVs inferred from DADA2 pipeline. The more 16S rRNA gene amplicon variants detected, the richer the community. On average, the number of the 16S rRNA gene copies per genome of the Firmicutes (F) tend to be twice as high as that of the Bacteroidetes (B) [33]. This estimate prompted us to revisit the difference in F:B ratio at day 7 PW. Indeed, the results of the ANCOM analysis at the phylum level supported that the F:B ratio was higher in the cecum samples of the control group at day 7 PW compared with that of the yeast group. Further, we attempted to verify that it was the intragenomic variability in the number of the 16S rRNA gene copies that contributed most to the Shannon diversity calculation. To achieve this, we compared the bacterial richness at the species level between the feeding groups. As expected, there was no difference in species richness between the feeding groups when ASVs were binned at the species level. Therefore, we conclude that the overrepresentation of Bacteroidetes-affiliated ASVs in the yeast large intestine resulted in a Shannon diversity that was lower compared with that of the control group. Firmicutes may appear in higher numbers in the microbiota of the control pigs because of differences in the diet formulation. In the yeast diet, yeast proteins mostly replaced the conventional protein sources used in the control diet, i.e. soybean meal, potato protein concentrate, fish meal, and rapeseed meal. The presence of various dietary fibres in a diet influences the gut microbiota composition by promoting the growth of Firmicutes (reviewed in [34]). It is known that up to 10% of soybean meal, and a considerable portion of rapeseed meal (> 15%) is neutral detergent fibre [7]. Thus, the higher level of soybean- and rapeseed meal in the control diet could account for the high presence of fibre-degrading Firmicutes in the large intestine of the control group.

According to our findings, a Prevotellaceae family-related amplicon, ASV50, was predominant in the large intestine of pigs in the yeast group. Even though the analysis of the microbiota composition was confined to sequencing of the 16S rRNA bacterial gene only, the overrepresentation of the Prevotellaceae family might be related to the availability of the non-digested parts of the yeast cells in the diet. It is conceivable that the method of yeast processing partially precluded its digestibility in the small intestine thus making the yeast cells available for microbial fermentation in the large intestine. The overgrowth of Prevotella in the yeast driven microbiota might also be attributed to the microbial peptidase and proteinase activities of this bacterial groups [35]. Previous studies by Mach et al. [36] and by Ramayo-Caldas et al. [37], showed that the enterotype dominated by Prevotella and Mitsuokella species is associated with lowered alpha diversity and improved growth performance. These findings are in line with our results on alpha bacterial diversity and ADG in the yeast group. The lower levels of butyrate and acetate in the colon of the yeast fed piglets may be due to the predominance of the Prevotella-Mitsuokella-affiliated groups and hence a suppression of certain SCFA-producers [37]. Higher abundance of Mitsuokella in the large intestine of yeast fed animals is consistent with the studies where yeast was supplemented [8]. Furthermore, in a study using similar dietary formulations as the present study, Cruz and co-workers [7] showed that the total tract digestibility of phosphorous was higher in the yeast group than in the control group. As Mitsuokella and Selenomonas genera are reported to release phosphorous from phytate [38, 39], it is tempting to ascribe this metabolic activity to these bacteria. The Mitsuokella and Selenomonas genera were found in higher numbers in the yeast group compared with the numbers in the control group. However, the resolution of the 16S rRNA gene method does not always provide enough confidence in assigning the PCR amplicons to the species level. To learn about the functional potential and contribution to host metabolism of the mentioned Selenomonadaceae and Prevotellaceae, use of anaerobic cultivation techniques may be necessary. It has previously been reported that in gut eco-systems supplied with low levels of yeast-derived components, the outgrowth of Prevotella, Selenomonas, and Mitsuokella commonly co-occurs with the reduction in SCFA producing bacteria [9,10,11]. The interpretation of the ANCOM analysis revealed more F. prausnitzii (97.5% identity to F. prausnitzii strain ATCC 27768, GenBank accession: NR_028961) in the caecum of the control group than that of the yeast group. To our knowledge, the ANCOM test we used for identification of differentially abundant taxa between the two dietary groups performs better than other tests with respect to false discovery rate control. However, when applied to groups with less than twenty samples per group, the sensitivity of ANCOM decreases [40]. In our 16S rRNA gene sequencing setup, F. prausnitzii on average represented 0.9% of the caecum microbiota population and 0.3% of the colon microbiota population. It is possible that eight samples per group were not enough for ANCOM to detect differences in rare colonic F. prausnitzii between the two groups. However, irrespective of the dietary interventions, we found that F. prausnitzii was positively correlated with the colonic butyrate concentration. The opposite relationship was found for the Oxalobacter genus. Oxalobacter strain OxB, an oxalate degrader, was studied by Allison et al. [41]. Allison and co-workers showed that acetate is an essential nutrient for growth of the bacterium. Their findings suggest that a competition for the nutrient between F. prausnitzii and Oxalobacter is conceivable [41, 42]. There was a mutual exclusion relationship between the two bacteria at a statistically significant level when examined with the CoNet co-occurrence network analysis [43] (data not shown). In addition to the role of F. prausnitzii in butyrate production, our results suggest that Dialister and an unclassified member of the Spirochaetaceae family may contribute to the colonic butyrate pool.

Next, we have found an association between the concentration of colonic butyrate and the liver index. Similarly, the number of F. prausnitzii in the colon correlated with the liver index. The portal vein concentration of butyrate is known to reflect the production levels of butyrate in the colon [44]. Thus, it is likely that the liver index was related to the uptake and metabolism of butyrate in the liver in our study. It is intriguing to think that gut microbiota members may be involved in the butyrate metabolism to the extent where the size of the liver is affected. Reduced level of butyric acid has been shown to be associated with pathologic conditions in man [45]. Butyrate has been implicated to play a role in the integrity of the intestinal wall, serving as an energy source for colonocytes and as a regulatory molecule [46, 47]. However, it is unclear what concentration of butyrate is optimal to maintain gut integrity and homeostasis in weanling piglets. Furthermore, butyrate has been shown to impact actively on the colonic crypt stem cells [18, 19]. Wang and co-workers demonstrated that butyrate diminished the crypt cell proliferation in a dose-response manner in an in vitro human colon crypt array [17]. In the present study, we observed a difference in the colonic CD between the two diet groups that possibly could be attributed to the altered abundance of butyrate-producing bacteria. Similar findings but of a lesser magnitude were reported by Mentschel and Claus in a study where piglets were fed with resistant potato starch [18]. In the light of our findings and previous publications [17,18,19], there is a good reason to believe that crypt elongation is a compensatory change to protect the crypt stem cell compartment from butyrate toxicity. The correlation between the colonic butyrate and the liver size suggests that the colonocytes received butyrate levels exceeding their metabolic capacity, with the butyrate surplus being transported to the liver. Histopathological examination of intestinal tissues did not reveal any difference in gut health parameters between the two feeding groups. Thus, it is tempting to speculate that there is a saturation point in the butyrate microbial production beyond which butyrate is not required as a fuel for colonocytes.

The cultivation results demonstrated that overall the differences between the feeding groups in the counts of LAB, enterococci, and coliforms were detected from day 7 PW. Much of the inter-individual variation before day 7 PW may be attributed to the weaning event. The bacterial succession of the gut is governed by, but not limited to, the substrate availability, gut physiology, and immune status. Feed intake during the time of weaning, when the piglets shift from milk to solid feed, is a key factor for immune system maturation [22], and luminal wall development [21]. The weaning event entails an irregular and variable timing in the acceptance of the new type of diet. This, consequently, leads to a transient starvation in some animals. According to our observations, albeit non-systematic, this was the case in our experiment. To this end, it is to be expected that the major variability in bacterial succession occurs during the first two weeks PW. LAB were consistently found in higher numbers in the ileum and large intestine of the yeast group compared with intestines of the control group. This bacterial group has a range of bioactive properties known to benefit mammals (reviewed in [48]). Attempts have been made to graft LAB into GI tract to improve health or ameliorate disease [49]. The ileum has a very dynamic gut environment, where bacteria must overcome multiple factors (e.g. digesta flow, peristalsis, microbe-host interaction, and microbe-microbe interaction) if they are to colonize and persist in the intestinal segment [32]. It has been reported that LAB are capable of adhesion to the intestinal cell wall [48]. Russo et al. showed that adhesion to human enterocytes of some LAB strains in vitro was inducible by β-D-glucan extracted from Pediococcus parvulus [50]. Therefore, the higher abundance of intestinal LAB in the yeast group in the present study may be attributable to the presence of a β-glucan fraction from the yeast cell wall in the feed. It is also tempting to speculate that the presence of yeast cell wall glucans in the feed affects digesta viscosity in the lumen, which is a factor that would favour LAB colonization. Snart and colleagues demonstrated that high-viscosity dietary fibre β-glucans supplementation was positively associated with higher numbers of lactobacilli in the caecum of rats [51]. Supplementation with the yeast cell wall was implicated in an increase in lactobacilli numbers in the ileal digesta of broilers in the studies by Liy et al. and Ghosh et al. [52, 53]. Their findings suggest that the yeast cell wall, or its components, may have selective properties towards LAB in a range of hosts. It is notable that L. salivarius 16S rRNA gene amplicon relative abundance was found in higher numbers in the ileum of the pigs fed yeast at day 14 PW compared with that of the control pigs. The concordance between the results obtained from culture-dependent and -independent methods strengthens the validity of our findings. The augmentation of intestinal LAB is a promising aspect of the yeast-derived diet in GI tract of pigs. However, further research is needed to elucidate whether it is the yeast wall β-glucans or other ingredients of the diet that favour the LAB increase.

Conclusions

The replacement of 40% of the crude protein from the main protein sources traditionally used in Norway with proteins from Cyberlindnera jadinii in a weanling piglet diet reshaped the large intestine microbiota structure. The microbiota of yeast fed piglets showed a dominance of Prevotella-, Mitsuokella- and Selenomonas-related taxa along with the decreased alpha-diversity. Larger numbers of viable LAB cells were recovered from both small and large intestines of the yeast fed piglets compared with the control piglets. Owing to the functional capacity of the above bacterial groups, we believe that Cyberlindnera jadinii yeast, in addition of being a high-quality protein source, promote growth of beneficial gut microbes.

Availability of data and materials

The raw sequencing reads are deposited in the SRA archive: PRJNA580284.

Abbreviations

16S rRNA :

16S ribosomal ribonucleic acid

ADG:

Average daily gain

ANCOM:

Analysis of compositions of microbiomes

ASV:

Amplicon sequence variant

CD:

Crypt depth

CFU:

Colony forming unit

DNA :

Deoxyribonucleic acid

HID-AB:

High iron diamine-alcian blue

LAB:

Lactic acid-producing bacteria

MWW:

Mann-Whitney-Wilcoxon

PERMANOVA:

Permutational multivariate analysis of variance

PW:

Post-weaning

PCoA:

Principle coordinate analysis

PCA:

Principle component analysis

SCFA:

Short-chain fatty acid

SE:

Standard error

SRA:

Sequence Read Archive

References

  1. Brack D, Glover A, Wellesley L. Agricultural Commodity Supply Chains: Trade, Consumption and Deforestation. London: Royal Institute of International Affairs (Chatham House); 2016;1–80.

  2. Braude R, Henry KM, Kon SK. Further studies of the Rachitogenic effect of dried yeast in pig diets. Br J Nutr. Cambridge University Press; 1948;2:66–75.

  3. Macrae TF, El-Sadr MM, Sellers KC. The nutritive value of yeast protein: comparison of the supplementary values of yeast protein and casein for maize protein in the nutrition of the pig. Biochem J. 1942;36:460–77 Available from: https://www.ncbi.nlm.nih.gov/pubmed/16747548.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Hoff-Jørgensen E. Yeast and rickets. Nature. 1947;159:99–100. https://doi.org/10.1038/159099b0.

    Article  PubMed  Google Scholar 

  5. Braude R, Kon SK, White EG. Yeast as a protein supplement for pigs; further observations on its rachitogenic effect. J Comp Pathol Ther. 1944;54:88–96. https://doi.org/10.1016/S0368-1742(44)80009-1.

    Article  CAS  Google Scholar 

  6. Sharma S, Hansen LD, Hansen JO, Mydland LT, Horn SJ, Øverland M, et al. Microbial Protein Produced from Brown Seaweed and Spruce Wood as a Feed Ingredient. J Agric Food Chem. Am Chem Soc. 2018;66:8328–35.

    Article  CAS  Google Scholar 

  7. Cruz A, Håkenåsen IM, Skugor A, Mydland LT, Åkesson CP, Hellestveit SS, et al. Candida utilis yeast as a protein source for weaned piglets: Effects on growth performance and digestive function. Livest Sci. 2019;226:31–9 Available from: https://linkinghub.elsevier.com/retrieve/pii/S1871141319301003.

    Article  Google Scholar 

  8. Kiros TG, Derakhshani H, Pinloche E, D’Inca R, Marshall J, Auclair E, et al. Effect of live yeast Saccharomyces cerevisiae (Actisaf Sc 47) supplementation on the performance and hindgut microbiota composition of weanling pigs. Sci Rep. 2018;8:1–13.

    Article  CAS  Google Scholar 

  9. Upadrasta A, O’Sullivan L, O’Sullivan O, Sexton N, Lawlor PG, Hill C, et al. The effect of dietary supplementation with spent cider yeast on the swine distal gut microbiome. PLoS One [Internet]. Public Library of Science; 2013;8:e75714. Available from: https://doi.org/10.1371/journal.pone.0075714.

  10. Fouhse JM, Dawson K, Graugnard D, Dyck M, Willing BP. Dietary supplementation of weaned piglets with a yeast-derived mannan-rich fraction modulates cecal microbial profiles, jejunal morphology and gene expression. Animal. 2019;13:1591–8.

  11. Nakashimada Y, Michinaka A, Watanabe K, Nishio N, Fujii T. Brewer’s yeast cell wall affects microbiota composition and decreases Bacteroides fragilis populations in an anaerobic gut intestinal model. J Biosci Bioeng. 2011;111:178–84. https://doi.org/10.1016/j.jbiosc.2010.09.005.

    Article  CAS  PubMed  Google Scholar 

  12. Furusawa Y, Obata Y, Fukuda S, Endo TA, Nakato G, Takahashi D, et al. Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature. 2013;504:446–50. https://doi.org/10.1038/nature12721.

    Article  CAS  PubMed  Google Scholar 

  13. Koh A, De Vadder F, Kovatcheva-Datchary P, Bäckhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. 2016;165:1332–45.

    Article  CAS  PubMed  Google Scholar 

  14. Roediger WE. Role of anaerobic bacteria in the metabolic welfare of the colonic mucosa in man. Gut. 1980;21:793–8 Available from: https://www.ncbi.nlm.nih.gov/pubmed/7429343.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Luciano L, Groos S, Busche R, Von Engelhardt W, Reale E. Massive apoptosis of colonocytes induced by butyrate deprivation overloads resident macrophages and promotes the recruitment of circulating monocytes. Cell Tissue Res. 2002;309:393–407.

    Article  CAS  PubMed  Google Scholar 

  16. Lupton JR. Diet induced changes in the colonic environment and colorectal cancer. J Nutr. 2004;134:479–82.

    Article  CAS  PubMed  Google Scholar 

  17. Wang Y, Kim R, Gunasekara DB, Reed MI, DiSalvo M, Nguyen DL, et al. Formation of human colonic crypt Array by application of chemical gradients across a shaped epithelial monolayer. Cmgh. 2018;5:113–30. https://doi.org/10.1016/j.jcmgh.2017.10.007.

    Article  PubMed  Google Scholar 

  18. Mentschel J, Claus R. Increased butyrate formation in the pig colon by feeding raw potato starch leads to a reduction of colonocyte apoptosis and a shift to the stem cell compartment. Metabolism. 2003;52:1400–5.

    Article  CAS  PubMed  Google Scholar 

  19. Kaiko GE, Ryu SH, Koues OI, Collins PL, Solnica-Krezel L, Pearce EJ, et al. The Colonic Crypt Protects Stem Cells from Microbiota-Derived Metabolites. Cell. 2016;165:1708–20. https://doi.org/10.1016/j.cell.2016.05.018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Murphy P, Dal Bello F, O’Doherty J, Arendt EK, Sweeney T, Coffey A. Analysis of bacterial community shifts in the gastrointestinal tract of pigs fed diets supplemented with β-glucan from Laminaria digitata , Laminaria hyperborea and Saccharomyces cerevisiae. Animal. 2013;7:1079–87.

    Article  CAS  PubMed  Google Scholar 

  21. van Beers-Schreurs HMG, Nabuurs MJA, Vellenga L, der Valk HJK, Wensing T, Breukink HJ. Weaning and the weanling diet influence the villous height and crypt depth in the small intestine of pigs and Alter the concentrations of short-chain fatty acids in the large intestine and blood. J Nutr. 1998;128:947–53 Available from: https://academic.oup.com/jn/article/128/6/947/4722367.

    Article  PubMed  Google Scholar 

  22. Horn N, Ruch F, Little CR, Miller G, Ajuwon KM, Adeola O. Impact of acute feed and water deprivation at weaning and subsequent heat stress on growth performance and ileal morphology in nursery pigs. J Anim Sci. 2016;94:289–93.

    Article  CAS  Google Scholar 

  23. Umu ÖCO, Fauske AK, Åkesson CP, De Nanclares MP, Sørby R, Press CML, et al. Gut microbiota profiling in Norwegian weaner pigs reveals potentially beneficial effects of a high-fiber rapeseed diet. PLoS One. 2018;13:e0209439.

    Article  Google Scholar 

  24. Knudsen BE, Bergmark L, Munk P, Lukjancenko O, Priemé A, Aarestrup FM, et al. Impact of Sample Type and DNA Isolation Procedure on Genomic Inference of Microbiome Composition. mSystems [Internet]. 2016;1:e00095-16. Available from: http://msystems.asm.org/content/1/5/e00095-16.abstract.

  25. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3 Available from: https://www.ncbi.nlm.nih.gov/pubmed/27214047.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72 Available from: https://www.ncbi.nlm.nih.gov/pubmed/16820507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Frøslev TG, Kjøller R, Bruun HH, Ejrnæs R, Brunbjerg AK, Pietroni C, et al. Algorithm for post-clustering curation of DNA amplicon data yields reliable biodiversity estimates. Nat Commun. 2017;8:1188 Available from: https://www.ncbi.nlm.nih.gov/pubmed/29084957.

    Article  PubMed  PubMed Central  Google Scholar 

  28. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217. https://doi.org/10.1371/journal.pone.0061217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, et al. Package vegan: Community Ecology Package [Internet]. R Package. version 2.3–1. 2013. p. 264. Available from: http://cran.r-project.org/package=vegan.

  30. Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663 Available from: https://www.ncbi.nlm.nih.gov/pubmed/26028277.

    PubMed  Google Scholar 

  31. Savage DC. Microbial ecology of the gastrointestinal tract. Annu Rev Microbiol. 1977;31:107–33.

    Article  CAS  PubMed  Google Scholar 

  32. Rolfe RD. Colonization resistance. In: Roderick M, Bryan W, Isaacson RE, editors. Gastrointest Microbiol Vol 2 Gastrointest microbes host Interact. 1st ed. US: Springer; 1996. p. 510–3.

    Google Scholar 

  33. Sun DL, Jiang X, Wu QL, Zhou NY. Intragenomic heterogeneity of 16S rRNA genes causes overestimation of prokaryotic diversity. Appl Environ Microbiol. 2013;79:5962–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Scott KP, Duncan SH, Flint HJ. Dietary fibre and the gut microbiota. Nutr Bull. 2008;33:201–11.

    Article  Google Scholar 

  35. Avgustin G, Wallace RJ, Flint HJ. Phenotypic diversity among Ruminal isolates of Prevotella ruminicola: proposal of Prevotella brevis sp. nov., Prevotella bryantii sp. nov., and Prevotella albensis sp. nov. and redefinition of Prevotella ruminicola. Int J Syst Bacteriol. 2009;47:284–8.

    Article  Google Scholar 

  36. Mach N, Berri M, Estellé J, Levenez F, Lemonnier G, Denis C, et al. Early-life establishment of the swine gut microbiome and impact on host phenotypes. Environ Microbiol Rep. 2015;7:554–69.

    Article  CAS  PubMed  Google Scholar 

  37. Ramayo-Caldas Y, Mach N, Lepage P, Levenez F, Denis C, Lemonnier G, et al. Phylogenetic network analysis applied to pig gut microbiota identifies an ecosystem structure linked with growth traits. ISME J. 2016;10:2973–7.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Tang HC, Sieo CC, Abdullah N, Mohamad R, Omar AR, Chong CW, et al. Production of phytase by mitsuokella jalaludinii in semi-solid state fermentation of agricultural by-products. Sains Malaysiana. 2018;47:277–86.

    CAS  Google Scholar 

  39. Yanke LJ, Bae HD, Selinger LB, Cheng KJ. Phytase activity of anaerobic ruminal bacteria. Microbiology. 1998;144:1565–73.

    Article  CAS  PubMed  Google Scholar 

  40. Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome Microbiome. 2017;5:1–18.

    Article  Google Scholar 

  41. Allison MJ, Dawson KA, Mayberry WR, Foss JG. Oxalate-degrading anaerobes that inhabit the gastrointestinal tract. Arch Microbiol. 1985;141:1–7 Available from: https://link.springer.com/content/pdf/10.1007%2FBF00446731.pdf.

    Article  CAS  PubMed  Google Scholar 

  42. Duncan SH, Hold GL, Harmsen HJM, Stewart CS, Flint HJ. Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. Nov., comb. nov. Int J Syst Evol Microbiol. 2002;52:2141–6.

    CAS  PubMed  Google Scholar 

  43. Faust K, Raes J. CoNet app: inference of biological association networks using Cytoscape. F1000Research. 2016;5:1519 Available from: https://www.ncbi.nlm.nih.gov/pubmed/27853510.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Knudsen KEB, Serena A, Canibe N, Juntunen KS. New insight into butyrate metabolism. Proc Nutr Soc. 2003;62:81–6.

    Article  CAS  Google Scholar 

  45. Imhann F, Vich Vila A, Bonder MJ, Fu J, Gevers D, Visschedijk MC, et al. Interplay of host genetics and gut microbiota underlying the onset and clinical presentation of inflammatory bowel disease. Gut. 2018;67:108–19 Available from: https://www.ncbi.nlm.nih.gov/pubmed/27802154.

    Article  CAS  PubMed  Google Scholar 

  46. Donohoe DR, Garge N, Zhang X, Sun W, O’Connell TM, Bunger MK, et al. Microbiome, the regulate, butyrate metabolism. Energy Cell Metab. 2011;13:517–26.

    Article  CAS  PubMed  Google Scholar 

  47. Feng W, Wu Y, Chen G, Fu S, Li B, Huang B, et al. Sodium butyrate attenuates diarrhea in weaned piglets and promotes tight junction protein expression in Colon in a GPR109A-dependent manner. Cell Physiol Biochem. 2018;47:1617–29.

    Article  CAS  PubMed  Google Scholar 

  48. Pessione E. Lactic acid bacteria contribution to gut microbiota complexity: lights and shadows. Front Cell Infect Microbiol. 2012;2:1–15.

    Article  Google Scholar 

  49. Preston K, Krumian R, Hattner J, Demontigny D, Stewart M, Gaddam S. Lactobacillus acidophilus CL1285, lactobacillus casei LBC80R and lactobacillus rhamnosus CLR2 improve quality-of-life and IBS symptoms: a double-blind, randomised, placebo-controlled study. Benef Microbes. 2018;9:697–706.

  50. Russo P, López P, Capozzi V, de Palencia PF, Dueñas MT, Spano G, et al. Beta-glucans improve growth, viability and colonization of probiotic microorganisms. Int J Mol Sci. 2012;13:6026–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Snart J, Bibiloni R, Grayson T, Lay C, Zhang H, Allison GE, et al. Supplementation of the diet with high-viscosity beta-glucan results in enrichment for lactobacilli in the rat cecum. Appl Environ Microbiol. 2006;72:1925–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ghosh TK, Haldar S, Bedford MR, Muthusami N, Samanta I. Assessment of yeast cell wall as replacements for antibiotic growth promoters in broiler diets: effects on performance, intestinal histo-morphology and humoral immune responses. J Anim Physiol Anim Nutr (Berl). 2012;96:275–84.

    Article  CAS  Google Scholar 

  53. Liu N, Wang JQ, Jia SC, Chen YK, Wang JP. Effect of yeast cell wall on the growth performance and gut health of broilers challenged with aflatoxin B 1 and necrotic enteritis. Poult Sci. 2018;97:477–84.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the Bacteriology group at the Department of Food Safety and Infection Biology at Norwegian University of Life Sciences for their support during sampling and in the culturing study. We are grateful to the Histology group at the Department of Basic Sciences and Aquatic Medicine Biology at Norwegian University of Life Sciences for their contribution during sampling and histopathology analysis. We also thank Mathieu Castex and Caroline Achard for critical comments on this manuscript.

Funding

SI was supported by a PhD grant from Department of Food Safety and Infection Biology at Norwegian University of Life Sciences. The authors thank Foods of Norway Centre for Research-based Innovation (The Research Council of Norway, Lysaker, Norway, grant number 237841/030) and the Centre’s industrial partners for financial support.

Author information

Authors and Affiliations

Authors

Contributions

CMP, LTM, and MØ contributed to the study conceptualization; CPA, CMP, HS, IMH, LTM, MØ, and OCOU contributed to the study design; CPA, CMP, OCOU, IMH, and SI contributed to data acquisition and data analysis; CPA, HS, IMH, LTM, MØ, OCOU, and SI were involved in data interpretation; SI wrote the original draft; all authors reviewed and edited the manuscript prior to submission; all authors read and approved the final version of the manuscript and approved publication.

Corresponding authors

Correspondence to Özgün C. O. Umu or Margareth Øverland.

Ethics declarations

Ethics approval and consent to participate

All animal management procedures were in compliance with the Norwegian Animal Welfare Act 10/07/2009 and the EU Directive on the Protection of Animals used for Scientific Purposes (2010/63/EU).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Comparison of bacterial CFUs on the selective agar plates. The values are the group medians (IQR) of logCFUs per gram of lumen contents. The bold font indicates statistically significant level (p < 0.05) of MWW test.

Additional file 2.

Relative abundance of bacterial genera in the sequenced mock community standards along with their expected abundance.

Additional file 3.

Alpha bacterial diversity as measured by a) Observed species and b) Shannon diversity index. The average values between the control and yeast group are represented as mean values with standard deviation (SD) along with medians with inter-quartile ranges (IQR). Comparison pairs that correspond to p-values less than 0.05 (MWW test) are given in bold.

Additional file 4.

Differentially abundant ASVs between the yeast and control diets (ileum, caecum). A ileum, d 7 PW, B ileum, d 14 PW, C ileum, d 14 PW (family level), D caecum, d 7 PW (phylum level), E caecum, d 14 PW (phylum level). All taxonomic entities appeared as differentially abundant at FDR = 0.05.

Additional file 5.

Differentially abundant ASVs between the yeast and control diets (caecum). A, caecum, d 7 PW, B caecum, d 14 PW, C caecum, d 7 PW (family level), D caecum, d 14 PW (family level). All taxonomic entities appeared as differentially abundant at FDR = 0.05.

Additional file 6.

Differentially abundant ASVs between the yeast and control diets (colon). A colon, d 7 PW, B colon, d 14 PW, C colon, d 7 PW (family level), D colon, d 14 PW (family level). All taxonomic entities appeared as differentially abundant at FDR = 0.05.

Additional file 7.

Association of colonic butyrate concentration with individual bacterial groups. Panel A: Correlation plot of colonic butyrate concentration (measured in μM per gram of intestinal contents) against F. prausnitzii relative abundance measured at days 7 and 14 PW (n = 32). The dots are coloured by the diet (control pink; yeast dark cyan). Pearson’s rho is reported above the regression line. Panel B: Principal component analysis performed on the relative abundance of Spirochaetaceae, Faecalibacterium, Dialister and molarities of butyrate in the colon of pigs measured at days 7 and 14 PW (n = 32 but 3 dots are not shown). The dots are coloured by the diet (control pink; yeast dark cyan). The vectors represent the degree of correlation between the bacterial groups data and the butyrate concentration data.

Additional file 8.

Association of liver index with individual bacterial groups, and colonic butyrate concentration. Panel A: Principal component analysis performed on the relative abundance of Spirochaetaceae, Faecalibacterium and the liver index of pigs measured at days 7 and 14 PW (n = 32 but 2 dots are not shown). The dots are coloured by the diet (control pink; yeast dark cyan). The vectors represent the degree of correlation between the bacterial groups data and the liver index data. Panel B: Correlation plot of colonic butyrate concentration (measured in μM per gram of intestinal contents) against liver index measured in the control group pigs at days 7 and 14 PW (n = 16). The dots are coloured by the diet (control pink). Pearson’s rho is reported above the regression line. Panel C: Correlation plot of colonic butyrate concentration (μM/gram of intestinal contents) against liver index measured in the yeast group pigs at days 7 and 14 PW (n = 16). The dots are coloured by the diet (yeast dark cyan). Pearson’s rho is reported above the regression line.

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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Iakhno, S., Umu, Ö.C.O., Håkenåsen, I.M. et al. Effect of Cyberlindnera jadinii yeast as a protein source on intestinal microbiota and butyrate levels in post-weaning piglets. anim microbiome 2, 13 (2020). https://doi.org/10.1186/s42523-020-00031-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-020-00031-x

Keywords