Next Article in Journal
A Novel Biofilm Model System to Visualise Conjugal Transfer of Vancomycin Resistance by Environmental Enterococci
Next Article in Special Issue
Encapsulated Mixture of Methyl Salicylate and Tributyrin Modulates Intestinal Microbiota and Improves Growth Performance of Weaned Piglets
Previous Article in Journal
Molecular Detection of Drug-Resistance Genes of blaOXA-23-blaOXA-51 and mcr-1 in Clinical Isolates of Pseudomonas aeruginosa
Previous Article in Special Issue
Effects of the Inclusion of Fermented Mulberry Leaves and Branches in the Gestational Diet on the Performance and Gut Microbiota of Sows and Their Offspring
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 4
10.3390/microorganisms9040787
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gut Microbiota Dynamics, Growth Performance, and Gut Morphology in Broiler Chickens Fed Diets Varying in Energy Density with or without Bacitracin Methylene Disalicylate (BMD)

by
Deborah Adewole
* and
Fisayo Akinyemi
Department of Animal Science and Aquaculture, Faculty of Agriculture, Dalhousie University, Truro, NS B2N 5E3, Canada
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(4), 787; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040787
Submission received: 17 March 2021 / Accepted: 7 April 2021 / Published: 9 April 2021
(This article belongs to the Special Issue Gut Microbiota Development in Farm Animals)
Browse Figures
Versions Notes

Abstract

:
High-energy-density diet could increase body weight at the expense of the intestinal health of the animals. In order to optimize production without negatively influencing the gut health of chickens, dietary supplementation with bacitracin methylene disalicylate (BMD) is a common feeding strategy adopted to enhance production performance and intestinal health. Studies have suggested that BMD could improve chicken growth performance and gut health through modulation of the gut microbiota. The current study investigated the effect of BMD supplementation in a normal-energy (NE) or high-energy (HE) diet on growth performance, organ weights, jejunal morphology, and gut microbiota of broiler chickens at different growth stages. Birds were allocated to four treatments: normal-energy basal diet (NE-BAS), normal-energy BMD diet (NE-BMD), high-energy basal diet (HE-BAS), and high-energy BMD diet (HE-BMD). In the starter phase, body weight and body weight gain were reduced significantly (p < 0.05) in chickens fed HE diets compared to those fed NE diets. The FCR was significantly higher (p < 0.05) in birds fed HE-BMD diets in the starter phase but lower (p < 0.05) during the grower phase when compared to other treatments. Moreover, the relative bursa weight increased significantly (p = 0.0220) among birds that received HE diets. Birds fed HE-BMD had greater villus height (p = 0.054) than NE-BMD group. Among the chickens fed the HE diets, those that received BMD treatment had a significantly increased (p = 0.003) villus width (13.3% increase) compared to those that received the basal diet. Improved population of Firmicutes was observed in chickens fed HE-BMD diet when compared to HE-BAS. Our results imply that BMD may be more effective in improving intestinal health when supplemented in a high-energy diet for broiler chickens.

1. Introduction

Dietary energy density can be referred to as the amount of available energy per unit weight [1]. In the poultry industry, dietary energy and nutrient density have been reported to have a significant impact on gut health and growth performance [2,3,4]. The impact of high-energy-density diets in poultry could be beneficial, and at the same time detrimental. High-energy-density diets have been shown to increase body weight gain [5,6] and improve feed conversion efficiency [7,8] in poultry. Lamot et al. [9] reported an increase in gain-to-feed ratio during the first week of life for broiler chickens exposed to higher diet densities. However, aside from weight gain, the intestinal health of poultry birds is equally important as this can affect the farmer’s cost of maintaining a healthy flock. Studies have shown that changes in dietary energy density induce rapid changes in the composition of bacteria that colonizes the intestinal tract of mammals [10]. Other studies with mice and rats have shown that high-energy diets can trigger microbiota dysbiosis due to an imbalance between energy intake and expenditure [11,12,13]. A recent study with Pekin ducklings revealed that dietary energy content altered microbiota composition and diversity in the cecum [14].
In order to optimize production without negatively influencing the gut health of chickens, dietary supplementation with bacitracin methylene disalicylate (BMD) is a common feeding strategy adopted to enhance production performance and intestinal health. In Canada, BMD is regarded as an antibiotic of medium importance (that is, it has alternatives available) [15]. Although there are plans to eliminate the preventive use of BMD because of the worldwide concern about antibiotic resistance and the quest to preserve the potency of antibiotics for human and animal health [16], BMD is still currently being used in chicken production as a preventive measure against diseases [15]. BMD increases body weight gain and feed efficiency and promotes the cecal microbiota diversity in several poultry species [17,18,19]. Broiler chickens fed BMD showed a reduction in the concentrations of coliforms and Lactobacillus salivarius in the lower part of the intestine, along with significant growth-promoting effects [20].
It is plausible that bacitracin coupled with a diet of high or low energy density may elicit beneficial effects on the intestinal health of chickens. While various studies have focused on the effects of BMD on intestinal microbiota and growth performance of chickens [20,21,22,23], research on the effect of energy density and its interaction with the use of antibiotics in broiler chickens is limited. Therefore, it is important to give attention to dietary energy, how it affects gut microbiota in relation to antibiotic usage, and its influence on the growth performance of chickens. We hypothesized that BMD would complement the shortcomings resulting from feeding chickens with either low or high energy alone. In the present study, we determined the effects of BMD, two levels of dietary energy, and their interactions on the growth performance, organ weights, jejunal morphology, and gut microbiota of broiler chickens at different growth stages.

2. Materials and Methods

2.1. Management and Housing

Animal care approval for the experimental methods was granted by Dalhousie University Animal Use and Care Committee, and all chickens were handled and cared for according to the recommendations of the Canadian Council of Animal Care [24]. Broiler chicks (Ross 308) obtained at 1 day old from a commercial source arrived at the Atlantic Poultry Research Center, Dalhousie University Faculty of Agriculture, Truro, NS, and chicks were weighed in groups of 25 birds. Each group was assigned to a floor pen (0.93 m × 2.14 m), at a stocking density of 0.076 m2/bird. Room temperature was monitored daily and was gradually reduced from 31 to 22.6 °C from day 0 to day 42. The lighting program was set to produce 18 h of light and 6 h of darkness throughout the experimental period, and illumination was gradually reduced from 20 lx on day 0 to 5 lx.

2.2. Experimental Diets

Chickens received feed and water ad libitum via a phase-feeding program, which consisted of starter (0–14 days), grower (15–24 days), and finisher (25–42 days). Diets were fed as crumbled pellets during the starter phase and as pellets during the grower and finisher phases. Birds were randomly sorted into four dietary treatments, consisting of eight replicate pens per treatment. The experiment was designed as a 2 × 2 factorial design consisting of two levels of energy: (1) normal energy (NE) density and (2) high energy (HE) density and two levels of BMD inclusion: (1) a basal diet containing 0% of BMD (Basal) and (2) Basal + 0.05% BMD. The HE diets were formulated to contain an additional 100 kcal/kg ME above the NE diets, similar to the study of Kindlein et al. [25]. The NE and HE diets were formulated to have a similar ratio of metabolizable energy (ME) to crude protein and digestible amino acids, as presented in Table 1.

2.3. Data and Sample Collection

Feed intake (FI), body weight (BW), and mortality were recorded weekly, and body weight gain (BWG) and feed conversion ratio (FCR) were calculated. On days 21, 36, and 42 of age, one chicken per pen (eight chickens per treatment) was randomly selected and euthanized by electrical stunning and exsanguination. Birds were bled for 3–5 min after carotid and jugular veins were cut. After slaughter, digesta from the pair of ceca were mixed, sampled, and stored in plastic RNAse- and DNAse-free tubes, placed in liquid nitrogen, and afterward kept at −80 °C until ready for DNA extraction and sequencing. On day 42, visceral organs, including spleen, ceca, liver, bursa of Fabricius, and heart, were harvested and weighed by trained personnel. Samples from the jejunum (1.5 cm length midway between the point of entry of the bile ducts and Meckel’s diverticulum) were collected and fixed in 10% neutral buffered saline for histomorphological processing.

2.4. Analysis of Diets and Histological Samples

2.4.1. Diet Analysis

Nitrogen content in the diets was determined using the combustion method [26], Method 990.03, with a nitrogen analyzer (Model Leco CN828 Carbon Nitrogen Determinator, St. Joseph, MO, USA) and CP was calculated as N × 6.25. The ether extract in samples was determined after hexane extraction [26] using Method 920.39 in an Ankom XT10 Fat Extractor system (Macedon, NY, USA). The metabolizable energy content of diets was determined at the Central Testing Laboratory, Winnipeg, MB, using their standard laboratory procedures.

2.4.2. Histological Measurements

Fixed jejunum samples were processed using the same procedure as Oladokun et al. [24]. Briefly, villus height (from the base of the intestinal mucosa to the tip of the villus excluding the intestinal crypt), villus width (halfway between the base and the tip), and crypt depth (from the base upward to the region of transition between the crypt and villi) were determined by taking approximately 10 measurements of each component per slide using an image processing and analysis system (Leica DC480; Leica Microsystems Imaging Solutions Ltd., Concord, ON, Canada).

2.5. DNA Extraction and 16S rRNA Gene Sequencing

Specimens were placed into a MoBio PowerMag Soil DNA Isolation aBead Plate. DNA was extracted following MoBio’s instructions on a KingFisher robot. Bacterial 16S rRNA genes were PCR-amplified with dual-barcoded primers targeting the V4 region (515F 5′-GTGCCAGCMGCCGCGGTAA-3′, and 806R 5′-GGACTACHVGGGTWTCTAAT-3′), as per the protocol of Kozich et al. [27]. Amplicons were sequenced with an Illumina MiSeq using the 300-bp paired-end kit (v.3). Sequences were denoised, taxonomically classified using Silva (v. 138) as the reference database, and clustered into 97%-similarity operational taxonomic units (OTUs) with the mothur software package (v. 1.44.1) [28], following the recommended procedure (https://www.mothur.org/wiki/MiSeq_SOP; accessed on 11 October 2020). The potential for contamination was addressed by co-sequencing DNA amplified from specimens and from template-free controls (negative control) and processing the extraction kit reagents the same way as the specimens. A positive control from “S00Z1-” samples consisting of cloned SUP05 DNA was also included. Operational taxonomic units were considered putative contaminants (and were removed) if their mean abundance in controls reached or exceeded 25% of their mean abundance in specimens.

2.6. Bioinformatics and Statistical Analyses

We sequenced 16Sv4 amplicons generated from cecal digesta samples on a MiSeq. MiSeq-generated Fastq files were quality-filtered and clustered into 97% similarity operational taxonomic unit (OTUs) using the mothur software package (http://www.mothur.org; accessed on 11 October 2020). The resulting dataset had 94,509 OTUs (including singletons). An average of 42,742 quality-filtered reads was generated per sample. Sequencing quality for R1 and R2 was determined using FastQC 0.11.5 (Supplementary Figures S1 and S2). The Shannon index was used to estimate alpha diversity [29] on raw OTU abundance tables after filtering out contaminants. The significance of diversity differences was tested with ANOVA or linear mixed model depending on the study design. To estimate beta diversity across samples, we excluded OTUs occurring with a count of less than 3 in at least 10% of the samples and then computed Bray–Curtis indices. Abundance-weighted sample pairwise differences were calculated using the Bray–Curtis dissimilarity. Bray–Curtis dissimilarity is calculated by the ratio of the summed absolute differences in counts to the sum of abundances in the two samples [30]. Beta diversity was visualized using principal coordinate analysis (PCoA) ordination, emphasizing differences across samples. Permutational multivariate analyses of variance (PERMANOVA) was used to assess variation in community structure, with treatment group as the main fixed factor, and 9999 permutations were used for significance testing [31]. We conducted all analyses in the R environment. DESeq2 R package was used to identify differentially abundant taxa among diet variables. The 2-sided Welch’s t-test [32] and Benjamini–Hochberg false discovery rate correction were adopted for the group analysis. We used STAMP software [33] to analyze the significant microbes among the treatments and different age groups. The principal component analysis (PCA) analysis of the microbial community relationship was also performed with STAMP software using ANOVA as statistical test (p-value < 0.05) and Tukey–Kramer as post hoc test. Further statistical analyses were performed using GraphPad prism Program (version 7.0.1). p-values < 0.05 were regarded as significant values. Results are presented as means ± standard deviation (SD). LEfSe was performed (http://huttenhower.sph.harvard.edu/galaxy; accessed on 20 November 2021) by using the relative abundance of genus to identify different taxa microbes between the time points and different diet groups. We only showed the differential features using the default parameters (LDA Score > 2, p < 0.05) [34].
The growth performance, jejunal morphology, and organ weights data were subjected to analysis of variance using the general linear models (GLM) procedure of SAS with energy level, BMD, and energy × BMD interaction as factors and the following parameters as variables: feed intake, body weight, body weight gain, feed conversion ratio, organ weights, villus heights, villus width, crypt depth, and ratio of villus height to crypt depth. p-values less than 0.05 were considered significant.

3. Results

3.1. Growth Performance, Organ Weight, and Jejunal Morphology

The effects of BMD and energy level on the growth performance of broiler chickens are shown in Table 2. There were no significant effects of dietary treatments on overall FI, BW, BWG, and FCR. However, during the starter phase, chickens fed BMD in the high-energy group had lower BW and BWG and higher FCR compared to those fed the basal diet in the normal-energy group, while other treatments were intermediate. During the starter phase, chickens in the high-energy BMD group had a higher FCR compared to those fed the basal diet in the normal-energy group. However, this was compensated for by a lower FCR during the grower phase compared to the other treatments. As illustrated in Table 3, bursa’s relative weight significantly increased (p = 0.0220) among the birds that received HE diets compared to those that received NE diets. There were no effects of BMD and energy level on the relative weights of the remaining organs. As presented in Table 4, birds fed HE-BMD had greater villus height (p = 0.054) than those provided with NE-BMD. Among the chicken birds fed the HE diets, those that received BMD treatment had a significantly increased (p = 0.003) villus width (13.3% increase) compared to those that received the basal diet. There were no effects of dietary treatment on crypt depth and VH/CD.

3.2. Composition of Chicken Gut Microbiota of Chickens Fed with or without BMD at Varying Energy Levels

A total of 94,509 OTUs were obtained from all samples. An average of 42,742 quality-filtered reads was generated per sample (Figure 1).
In the current study, we examined the effects of BMD on the chicken gut bacterial community at different energy levels. Alpha diversity analysis using the Shannon index, which accounts for richness and evenness, showed no differences among the four treatments, indicating that BMD did not influence microbiome evenness at different energy levels (Figure 2a). The PCA analysis revealed a weak influence of the treatments on gut microbial composition in the sampled chickens. A high level of similarity was observed among the treatments (Figure 2b). Additionally, we assessed the relative abundance among all the treatment groups at different taxonomic levels. The top classes were Clostridia, Bacteriodia, and Bacilli. Clostridia emerged as the most abundant in chickens fed NE, with a proportion of 58.9 and 58.0% for basal and BMD diets, respectively. Bacteriodia emerged as the most abundant in chickens fed high-energy diets, with a proportion of 34.8 and 32.0% for basal and BMD diets, respectively. Bacilli were most dominant in chickens fed NE (8.99 and 8.09% for basal and BMD groups, respectively) (Figure 2c). At the phylum level, the top two phyla among all treatments were Firmicutes and Bacteroidota. Firmicutes were more dominant in chickens fed NE diets, followed by those fed HE diets. The relative abundance was 69.6, 67.7, 63.3, and 65.2% for NE-BAS, NE-BMD, HE-BAS, and HE-BMD, respectively. Bacteroidota was seen as the most abundant in chickens fed HE diets, which include HE-BAS (34.77%) and HE-BMD (32.04%). The lowest proportion of Bacteroidota (27.95%) was observed among the chickens fed NE-BAS, while 29.72% was seen in NE-BMD (Figure 2d).
At the genus level, the dominant genera among the chickens fed with NE-BAS, NE-BMD, HE-BAS, and HE-BMD were Lachnospiraceae_unclassified, Bacteroides, Faecalibacterium, Clostridia_unclassified, Alistipes, Clostridia_UCG_014_ge, Lactobacillus, Bacterodales, Peptostreptococcacea_unclassified, and Blautia. The most dominant genus among chickens fed NE-BMD was Lachnospiraceae_unclassified (24.08%), while Bacteroides was the most abundant among chickens fed HE-BAS (20.11%). Irrespective of dietary energy density, we observed almost the same proportion of Faecalibacterium among chickens that received BMD or BAS (17.1 and 17.0%, respectively). Alistipes was most dominant in chickens fed BMD diets, irrespective of energy density. Chickens fed NE-BAS had the lowest proportion of Alistipes and the highest proportion of Lactobacillus, while those fed NE-BMD had the lowest abundance of Lactobacillus (Figure 2e). Two out of 83 genera were significantly different in abundance between the chickens fed BAS and BMD. Streptococcus was significantly abundant in chickens fed with BAS, while BMD inclusion triggered the drastic reduction of this microbe. Oscillospirales_ge was found to be differentially abundant in the chickens that received BMD treatment (Figure 2f). Moreover, we observed a reduction in the relative abundance of some intestinal bacteria in the BMD treatment group.
We identified 15 OTUs that were differentially abundant amongst treatments (NE-BAS, NE-BMD, HE-BAS, and HE-BMD) (FDR < 0.05), namely Otu00013, Otu00008, Otu00210, Otu00016, Otu00081, Otu00004, Otu00309, Otu00004, Otu00309, Otu00209, Otu00001, Otu00039, Otu00280, Otu00071, Otu00079, Otu00229, and Otu00150, which belong to the two phyla Firmicutes and Bacteroidota (Table 5).
Six phyla and 52 genera were observed in the broiler chickens. Further analysis using t-test was performed to evaluate the significantly abundant genera. Among the NE group, Clostridia_UCG_014 was significantly different in abundance between chickens fed with basal and BMD diets. Streptococcus, Oscillospirales, Fusicatenibacter, and Ruminococcacea_unclassified were significantly different in abundance between chickens fed NE-BMD and HE-BAS. Among chickens fed HE, Streptococcus, Oscillospirales, and Fusicatenibacter were differential between basal and BMD diets. A higher abundance of Escherichia–Shigella was found to be significant in chickens fed NE-BMD diet. Bacilli_unclassified was significantly enriched (p < 0.05) in chickens fed NE-BAS (Figure 3).

3.3. Dynamic Changes in Microbial Taxa across Growth Stages

The principal component analysis (PCA) analysis showed significant differences across growth stages. We observed a high level of similarities between days 36 and 43 compared to day 21 (Figure 4a). Alpha diversity index using Shannon revealed a significant difference between days 21 and 43 and between days 36 and 43 (p < 0.05), although the difference between days 21 and 36 was not significant (Figure 4b). Moreover, the PCoA plot obtained from beta analysis showed that the interaction between diets and age influenced the microbial community of the broiler chickens. Greater similarity among diets on days 36 and 43 compared to day 21 was observed. Moreover, samples from different diets on day 21 clustered closer than those on days 36 and 43 (Figure 4c).
From the permutational analysis of variance (PERMANOVA) that was conducted to determine significant differences in beta diversity among sampling factors (Supplementary Table S1), we observed there were no significant differences in bacteria composition along with the treatments. Further post hoc pairwise testing indicated that there was a significant difference (p-value < 0.05) in gut microbial flora across the treatments at different days (Supplementary Table S2), which can be visualized by the above PCoA plot (Figure 4c). Pairwise contrasts are presented in Supplementary Table S2, and the FDR method was used to correct p-values for multiple comparisons.
LEfSe analysis result revealed the effect size of some significant taxa on different days (Figure 4d). Additionally, Welch’s t-test using Stamp software displayed the relative abundance of significant microbes among days. Specifically, 17 genera changed significantly between days 21 and 36, while 3 changed between days 36 and 43 (q-value < 0.05). Some of the genera that changed significantly were Alistipes, Odoribacter, Bacillales_unclassified, Hungateiclostridiaceae_unclassified, Corynebacterium, and Parasutterella (Figure 4e).
We further compared the abundance of gut microbiota in the broiler chickens at different growth stages. The most abundant phyla in all three time points were Firmicutes and Bacteroidota. On day 21, the proportion of Firmicutes was 86.3%, and the proportion of this phylum was reduced drastically to 52.8% and then increased to 61.0% on day 43. The highest proportion of Bacteroidota (44.49%) was seen on day 36, and the proportion of this phylum was later reduced to 36.6% on day 43. The lowest proportion of Bacteroidota (11.7%) was observed on day 21 (Figure 4f). A similar occurrence was seen at the genus level. The most abundant genera were Faecalibacterium, Lachnospiraceae_unclassified, Bacteroides, Alistipes, and Clostridia_unclassified. On day 21, the proportion of Faecalibacterium was 35.4%, and this genus also emerged as the most abundant. A minimal proportion of Faecalibacterium was observed on days 36 and 43 (8.45 and 6.75%, respectively). On day 36, Bacteroides emerged as the most abundant with a proportion of 26.6%, followed by Lachnospiraceae_unclassified (21.8%). On day 43, Lachnospiraceae_unclassified was the most dominant (26.6%). Notably, the relative abundance of Lachnospiraceae_unclassified increased with growth in broiler chickens (18.1, 21.8, and 26.5% on days 21, 36, and 46, respectively). Moreover, the relative abundance of Alistipes increased as broiler chickens increased with age. A low proportion of Alistipes (1.77%) was seen on day 21, and the proportion of this genus increased to 8.67% on day 36 and 13.50% on day 43. On the contrary, the relative abundance of Faecalibacterium decreased with growth changes (Figure 4g).

4. Discussion

4.1. Effect of BMD and Energy Density on Growth Performance, Organ Weight, and Jejunal Morphology

High-energy diet has been reported to improve BW, BWG [35,36], and FCR [37] of broiler chickens. In the current study, there was no significant effect of dietary energy on overall growth performance. There was an increase in FCR among the chickens fed high-energy diets during the starter phase, but this was later compensated for during the grower and finisher phases. Our result supports a previous study by Coon et al. [38] who found that chick weight gain for 0–28 days was not significantly different for chicks fed high-energy diet. Moreover, FCR of chicks fed high-energy starter diet was superior to those fed low-energy diet [38]. Wang et al. (2014) [37] noted that high nutrient density improves FCR during the grower feeding phase. Additional findings in our study showed that BW and BWG were not significantly influenced by BMD, which is in agreement with a previous study by Damron et al. [39]. On the contrary, BMD improved the BW in broilers [40]. In the current study, we found that BMD significantly improved FCR in chickens on days 25–35. Our result is supported by a previous study that also established that FCR was significantly improved in birds treated with bacitracin [41]. Sims et al. [17] affirmed that at week 18, FCR was significantly lower in chickens fed with BMD than in the control group. Surprisingly, there was no interaction between BMD and energy level on growth performance parameters in our study.
The relative weight of immune organs is one of the primary determining factors of poultry birds’ immune status [42]. Our findings indicated no significant differences among the treatments in the relative weight of organs except for bursa weight, which was significantly increased (p = 0.0220) among the chicken birds that received high-energy diets compared to control. The bursa of Fabricius (BF) has been studied and revealed as a primary lymphoid organ in birds; it plays a vital role in the differentiation of B lymphocytes for immunoglobulin production during an immune challenge [43] and the development of the adaptive immune response of birds [44]. Conflicting with our findings, Cho et al. [45] concluded that emulsifier and multienzyme in low-density diets could increase the relative bursa weight. A recent study reported that a low-energy-density diet could increase chicken bursa weight [46]. High-energy diet could also be associated with the roles of the bursa in poultry birds. Chickens fed with a high-energy diet with BMD had significantly increased villus height and width in the current study. This observation indicates that BMD may be more effective in improving intestinal health when a high-energy diet is fed to chickens compared to when a normal-energy diet is fed. Adebowale et al. [47] reported that high dietary energy density induced stress in piglets and reduced villus height/crypt depth ratio in the ileum and duodenum. This indicates that a high-energy diet could be stressful to the chickens’ intestinal functionality. In the current study, the negative effect of feeding high-energy diet on the intestinal structure was ameliorated by the administration of BMD to the chickens. This agrees with previous studies that reported that low doses of BMD increased villus height through the small intestine in chickens [48]. The increase in villus height and surface area indicates an increase in cell proliferation, which may contribute to intestinal epithelial integrity maintenance [49].

4.2. Effect of BMD and Energy Density on Gut Microbiota Dynamics

Diets and antibiotics have been reported to shape the composition of gut microbiota [50]. Antibiotics sometimes have continuing effects on microbial flora in the gut, leading to a decrease in the abundance of beneficial microbes and promoting harmful bacteria [51]. Contrarily, a study revealed that the addition of BMD or Bacillus subtilis significantly reduced the population of harmful bacteria in the intestine [40]. Dietary components are vital for poultry gut health [52]. We investigated the effects of BMD on chicken gut microbiota at different dietary energy levels. The most abundant phyla in broiler chickens were Firmicutes and Bacteroidota, which is consistent with the results of our previous study [7]. Firmicutes produce molecules directly absorbed by the host gut wall as an energy source, and their abundance has been linked to weight gain in chickens [53,54]. It is worth noting that the inclusion of BMD into HE diets improved the population of Firmicutes slightly, compared to the control (HE-BAS), in this current study. Additionally, the abundance of Bacteroidota in broiler chickens that we observed in chickens fed a high-energy diet could result from this microbe’s involvement in the fermentation of dietary carbohydrates to produce short-chain fatty acids (SCFAs) [55]. SCFAs are the main end products of fermentation of nondigestible carbohydrates that colonize the human gut; they serve as sources of energy and have lasting impacts on host physiology [56].
We observed a slight reduction in the relative abundance of intestinal bacteria in the BMD treatment group, which is consistent with our previous findings [7]. The reduction in the relative abundance of Bacteroidetes observed among chickens fed with BMD agrees with a recent study that reported that antibiotic treatment reduced Bacteroidetes diversity in infants [57]. Bacteroidetes members are essential in developing a stable and healthy gut microbiota [58]. Moreover, the BMD effect was seen to significantly decrease Streptococcus abundance in our study. Streptococcus regulates the production of bacteriocin through quorum sensing [59]. Recently, the leading causative agent of dental caries in humans has been linked with Streptococcus mutants [60]. Further research is required to establish the exact role of Streptococcus in chickens.
We detected significant differences in microbial abundance among dietary treatment groups using STAMP software. Oscillopirales_ge abundance was significantly increased in chickens fed with high- and normal-energy diets containing BMD compared to those fed with high-energy basal diet. This suggests that BMD may contribute to indigestion in chickens. A high abundance of Oscillopirales may play an increasing role in constipation [61]. A significant increase in Streptococcus in chickens fed a normal-energy basal diet inhibited the growth of Oscillospirales. In contrast, a high-energy diet with BMD drastically reduced the abundance of Streptococcus. This phenomenon supports a previous study that established negative correlations between Streptococcus and some bacteria. Streptococcus has been found to produce metabolites that work against or kill some other bacteria [62,63], such as Oscillospirales. Meanwhile, the abundance of some microbes, such as Faecalibacterium, was not affected by diet but changed dynamically with time.
One of the essential factors that influences gut bacterial composition is the birds’ age [64]. We examined how the gut microbiota changes dynamically with growth stages and found that Faecalibacterium decreased as chickens’ age increased in our study, which is inconsistent with the findings of Donaldson et al. [65], which recorded an increase in the relative abundance of Faecalibacterium with time, though in slow proportion. The enrichment of Faecalibacterium is associated with a healthy human adult [65] and might influence chicken immune development [66]. The variation in the two studies might be due to diet or antibiotic treatment, which has been described to alter gut microbiota’s composition despite its stability [65]. The proportion of Alistipes increased as broiler chickens increased with age in the present study. Alistipes are commensal bacteria that belong to the phylum Bacteroidata and are assumed to benefit the host gut [67]. We further observed significant microbes at different time points of our study (days 21, 36, and 43). Ruminococcus was significantly abundant on day 36 of our study. Previous studies have reported that Ruminococcus might play a role in increasing mucin production by goblet cells [68] and improving host resistance to harmful bacterial invasion in calves [69]. The high abundance of Bacteroides and Ruminococcus suggests a healthy gut microbiota [70]. Lactobacillus was significantly abundant on day 43 based on the LEfSe result in our study. Increased Lactobacillus abundance could help inhibit some pathogens due to the production of vitamins and organic acids by this beneficial probiotic [71,72]. In the current study, the decrease in Lactobacillus abundance resulting from BMD inclusion signifies an adverse effect of antibiotics on beneficial microbes. A decline in Lactobacillus abundance is consistent with an earlier study that reported that antimicrobial growth promoters were associated with a reduction in the abundance of Lactobacillus species [73]. An upsurge of lactobacilli may be related to broiler growth depression in relation to competing for nutrient uptake or fat absorption in broiler chickens [10].
Our results showed that there were some interactions between diets and age, based on post hoc pairwise test, that influence the microbial community of the broiler chickens. This finding is consistent with a previous study that has established that the chicken gut microbiota is mostly shaped by the age and diet of the birds [64].

5. Conclusions

In conclusion, dietary supplementation of BMD in a high-energy diet improved FCR during the grow-out phase, increased villus height and width, and increased the relative abundance of beneficial gut microbiota members such as Firmicutes in broiler chickens. This study is novel because, to our knowledge, it is the first to present the effects of BMD in a high-energy-density diet on growth performance, organ weights, jejunal morphology, and gut microbiota of broiler chickens. Our results imply that BMD may be more effective in improving intestinal health when supplemented in a high-energy diet for broiler chickens. It also affirmed that chicken gut microbiota is mostly shaped by age and diet.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms9040787/s1, Figure S1: Sequences quality for R1; Figure S2: Sequences quality for R2; Table S1: Beta-diversity among sampling factors results from PERMANOVA; Table S2: Post-Hoc Pairwise test contrasts for multiple comparisons.

Author Contributions

Data curation, F.A.; funding acquisition, D.A.; investigation, D.A.; methodology, D.A.; project administration, D.A.; resources, D.A.; software, D.A. and F.A.; supervision, D.A.; writing—original draft, D.A. and F.A.; writing—review and editing, D.A. and F.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NSERC Discovery Grant, grant number 34288, and Agriculture and Agri-Food Canada (Pan Atlantic Program), grant number 53631. The APC was funded by Dalhousie University Start-Up grant (34741).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Canadian Council of Animal Care, and approved by Dalhousie University Animal Care Committee (protocol code 2019-008; March 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge Erin Maxwell for animal care, Janice MacIsaac for help with diet formulation, and Samson Oladokun for helping with gut morphology analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Drewnowski, A.; Almiron-Roig, E.; Marmonier, C.; Lluch, A. Dietary energy density and body weight: Is there a relationship? Nutr. Rev. 2004, 11, 403–413. [Google Scholar] [CrossRef]
  2. Abdollahi, M.R.; Zaefarian, F.; Ravindran, V.; Selle, P.H. The Interactive Influence of Dietary Nutrient Density and Feed form on the Performance of Broiler Chickens. Anim. Feed Sci. Technol. 2018, 239, 33–43. [Google Scholar] [CrossRef]
  3. Brickett, K.E.; Dahiya, J.P.; Classen, H.L.; Gomis, S. Influence of Dietary Nutrient Density, Feed form, and Lighting on Growth and Meat Yield of Broiler Chickens. Poult. Sci. 2007, 86, 2172–2181. [Google Scholar] [CrossRef]
  4. Zeng, Q.F.; Cherry, P.; Doster, A.; Murdoch, R.; Adeola, O.; Applegate, T.J. Effect of Dietary Energy and Protein Content on Growth and Carcass Traits of Pekin Ducks. Poult. Sci. 2015, 94, 384–394. [Google Scholar] [CrossRef] [PubMed]
  5. Fan, H.P.; Xie, M.; Wang, W.W.; Hou, S.S.; Huang, W. Effects of Dietary Energy on Growth Performance and Carcass Quality of White Growing Pekin Ducks from Two to Six Weeks of Age. Poult. Sci. 2008, 87, 1162–1164. [Google Scholar] [CrossRef] [PubMed]
  6. Grover, R.M.; Anderson, D.L.; Damon, R.A.; Ruggles, L.H. The Effects of Bird Density, Dietary Energy, Light Intensity, and Cage Level on the Reproductive Performance of Heavy Type Chickens in Wire Cages. Poult. Sci. 1972, 51, 565–575. [Google Scholar] [CrossRef]
  7. Adewole, D. Effect of Dietary Supplementation with Coarse or Extruded Oat Hulls on Growth Performance, Blood Biochemical Parameters, Ceca Microbiota and Short Chain Fatty Acids in Broiler Chickens. Animals 2020, 9, 1160. [Google Scholar] [CrossRef] [PubMed]
  8. Kuttappan, V.A.; Brewer, V.B.; Apple, J.K.; Waldroup, P.W.; Owens, C.M. Influence of Growth Rate on the Occurrence of White Striping in Broiler Breast Fillets. Poult. Sci. 2012, 91, 2677–2685. [Google Scholar] [CrossRef]
  9. Lamot, D.M.; Sapkota, D.; Wijtten, P.J.A.; Van Den Anker, I.; Heetkamp, M.J.W.; Kemp, B.; Van Den Brand, H. Diet Density during the First Week of Life: Effects on Growth Performance, Digestive Organ Weight, and Nutrient Digestion of Broiler Chickens. Poult. Sci. 2019, 98, 789–795. [Google Scholar] [CrossRef]
  10. Daniel, H.; Gholami, A.M.; Berry, D.; Desmarchelier, C.; Hahne, H.; Loh, G.; Mondot, S.; Lepage, P.; Rothballer, M.; Walker, A.; et al. High-Fat Diet Alters Gut Microbiota Physiology in Mice. ISME J. 2014, 8, 295–308. [Google Scholar] [CrossRef]
  11. Hildebrandt, M.A.; Hoffmann, C.; Sherrill-Mix, S.A.; Keilbaugh, S.A.; Hamady, M.; Chen, Y.Y.; Knight, R.; Ahima, R.S.; Bushman, F.; Wu, G.D. High-Fat Diet Determines the Composition of the Murine Gut Microbiome Independently of Obesity. Gastroenterology 2009, 137, 1716–1724.e2. [Google Scholar] [CrossRef] [Green Version]
  12. Minaya, D.M.; Turlej, A.; Joshi, A.; Nagy, T.; Weinstein, N.; DiLorenzo, P.; Hajnal, A.; Czaja, K. Consumption of a High Energy Density Diet Triggers Microbiota Dysbiosis, Hepatic Lipidosis, and Microglia Activation in the Nucleus of the Solitary Tract in Rats. Nutr. Diabetes 2020, 10, 20. [Google Scholar] [CrossRef]
  13. De La Serre, C.B.; Ellis, C.L.; Lee, J.; Hartman, A.L.; Rutledge, J.C.; Raybould, H.E. Propensity to High-Fat Diet-Induced Obesity in Rats Is Associated with Changes in the Gut Microbiota and Gut Inflammation. Am. J. Physiol. Gastrointest. Liver Physiol. 2010, 299, G440–G448. [Google Scholar] [CrossRef]
  14. Liu, J.B.; Yan, H.L.; Zhang, Y.; Hu, Y.D.; Zhang, H.F. Effects of Dietary Energy and Protein Content and Lipid Source on Growth Performance and Carcass Traits in Pekin Ducks. Poult. Sci. 2019, 98, 4829–4837. [Google Scholar] [CrossRef] [PubMed]
  15. Government of Canada. Categorization of Antimicrobial Drugs Based on Importance in Human Medicine. 2009. Available online: https://www.canada.ca/en/health-canada/services/drugs-health-products/veterinary-drugs/antimicrobial-resistance/categorization-antimicrobial-drugs-based-importance-human-medicine.html#shr-pg0 (accessed on 20 December 2020).
  16. Chicken Farmers of Canada Annual Report. 2019. Available online: https://www.chickenfarmers.ca/wp-content/uploads/2020/03/Annual-report-2019-%E2%80%93-E-Web.pdf (accessed on 20 December 2020).
  17. Sims, M.D.; Dawson, K.A.; Newman, K.E.; Spring, P.; Hooge, D.M. Effects of Dietary Mannan Oligosaccharide, Bacitracin Methylene Disalicylate, or Both on the Live Performance and Intestinal Microbiology of Turkeys. Poult. Sci. 2004, 83, 1148–1153. [Google Scholar] [CrossRef] [PubMed]
  18. Huyghebaert, G.; Ducatelle, R.; Immerseel, F. Van. An Update on Alternatives to Antimicrobial Growth Promoters for Broilers. Vet. J. 2011, 187, 182–188. [Google Scholar] [CrossRef] [Green Version]
  19. Crisol-Martínez, E.; Stanley, D.; Geier, M.S.; Hughes, R.J.; Moore, R.J. Understanding the Mechanisms of Zinc Bacitracin and Avilamycin on Animal Production: Linking Gut Microbiota and Growth Performance in Chickens. Appl. Microbiol. Biotechnol. 2017, 101, 4547–4559. [Google Scholar] [CrossRef] [PubMed]
  20. Engberg, R.M.; Hedemann, M.S.; Leser, T.D.; Jensen, B.B. Effect of Zinc Bacitracin and Salinomycin on Intestinal Microflora and Performance of Broilers. Poult. Sci. 2000, 79, 1311–1319. [Google Scholar] [CrossRef]
  21. Proctor, A.; Phillips, G.J. Differential Effects of Bacitracin Methylene Disalicylate (BMD) on the Distal Colon and Cecal Microbiota of Young Broiler Chickens. Front. Vet. Sci. 2019, 6, 114. [Google Scholar] [CrossRef]
  22. Gong, J.; Yu, H.; Liu, T.; Gill, J.J.; Chambers, J.R.; Wheatcroft, R.; Sabour, P.M. Effects of Zinc Bacitracin, Bird Age and Access to Range on Bacterial Microbiota in the Ileum and Caeca of Broiler Chickens. J. Appl. Microbiol. 2008, 104, 1372–1382. [Google Scholar] [CrossRef]
  23. Oladokun, S.; Koehler, A.; MacIsaac, J.; Ibeagha-Awemu, E.M.; Adewole, D.I. Bacillus Subtilis Delivery Route: Effect on Growth Performance, Intestinal Morphology, Cecal Short-Chain Fatty Acid Concentration and Cecal Microbiota in Broiler Chickens. Poult. Sci. 2021, 100, 100809. [Google Scholar] [CrossRef] [PubMed]
  24. CCAC. The Care and Use of Farm Animals in Research, Teaching and Testing; CCAC: Ottawa, ON, Canada, 2009.
  25. Kindlein, L.; Ferreira, T.Z.; Driemeier, D.; Nascimento, V.P.D.; Vieira, S.L.; Moraes, L.E.F.D.D.; King, A.; Sainz, R.D. Occurrence and Severity of White Striping in Broilers Until 50d of Age Fed with High and Low-Energy Diets: Body Weight, Histopathological Changes and Meat Quality. J. Vet. Sci. Technol. 2017, 8, 478. [Google Scholar]
  26. AOAC. Official Methods of Analysis, 15th ed.; TrAC Trends Analytical Chemist: Washington, DC, USA, 1990. [Google Scholar]
  27. Kozich, J.J.; Westcott, S.L.; Baxter, N.T.; Highlander, S.K.; Schloss, P.D. Development of a Dual-Index Sequencing Strategy and Curation Pipeline for Analyzing Amplicon Sequence Data on the MiSeq Illu- mina Sequencing Platform. Appl. Environ. Microbiol. 2013, 79, 5112–5120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  29. Weaver, W. Recent Contributions to the Mathematical Theory of Communication. In The Mathematical Theory of Communication; Institute of General Semantics: Forest Hills, NY, USA, 1949. [Google Scholar]
  30. Bray, J.R.; Curtis, J.T. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 1957, 27, 325–349. [Google Scholar] [CrossRef]
  31. Anderson, M.J.; Walsh, D.C.I. PERMANOVA, ANOSIM, and the Mantel Test in the Face of Heterogeneous Dispersions: What Null Hypothesis Are You Testing? Ecol. Monogr. 2013, 83, 557–574. [Google Scholar] [CrossRef]
  32. Welch, B.L. The Generalisation of Student’s Problems When Several Different Population Variances Are Involved. Biometrika 1947, 34, 28–35. [Google Scholar]
  33. Parks, D.H.; Beiko, R.G. Identifying Biologically Relevant Differences between Metagenomic Communities. Bioinformatics 2010, 26, 715–721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic Biomarker Discovery and Explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhao, J.P.; Chen, J.L.; Zhao, G.P.; Zheng, M.Q.; Jiang, R.R.; Wen, J. Live Performance, Carcass Composition, and Blood Metabolite Responses to Dietary Nutrient Density in Two Distinct Broiler Breeds of Male Chickens. Poult. Sci. 2009, 88, 2575–2584. [Google Scholar] [CrossRef]
  36. Perween, S.; Kumar, K.; Chandramoni; Kumar, S.; Singh, P.K.; Kumar, M.; Dey, A. Effect of Feeding Different Dietary Levels of Energy and Protein on Growth Performance and Immune Status of Vanaraja Chicken in the Tropic. Vet. World 2016, 9, 893–899. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, X.; Peebles, E.D.; Zhai, W. Effects of Protein Source and Nutrient Density in the Diets of Male Broilers from 8 to 21 Days of Age on Their Subsequent Growth, Blood Constituents, and Carcass Compositions. Poult. Sci. 2014, 93, 1463–1474. [Google Scholar] [CrossRef]
  38. Coon, C.N.; Becker, W.A.; Spencer, J.V. The Effect of Feeding High Energy Diets Containing Supplemental Fat on Broiler Weight Gain, Feed Efficiency, and Carcass Composition. Poult. Sci. 1981, 60, 1264–1271. [Google Scholar] [CrossRef]
  39. Damron, B.L.; Wilson, H.R.; Fell, R.V. Growth and Performance of Broiler Breeders Fed Bacitracin Methylene Disalicylate and Zinc Bacitracin. Poult. Sci. 1991, 70, 1487–1492. [Google Scholar] [CrossRef] [PubMed]
  40. Manafi, M.; Khalaji, S.; Hedayati, M.; Pirany, N. Efficacy of Bacillus Subtilis and Bacitracin Methylene Disalicylate on Growth Performance, Digestibility, Blood Metabolites, Immunity, and Intestinal Microbiota after Intramuscular Inoculation with Escherichia Coli in Broilers. Poult. Sci. 2017, 96, 1174–1183. [Google Scholar] [CrossRef] [PubMed]
  41. Yang, C.; Kennes, Y.M.; Lepp, D.; Yin, X.; Wang, Q.; Yu, H.; Yang, C.; Gong, J.; Diarra, M.S. Effects of Encapsulated Cinnamaldehyde and Citral on the Performance and Cecal Microbiota of Broilers Vaccinated or Not Vaccinated against Coccidiosis. Poult. Sci. 2020, 99, 936–948. [Google Scholar] [CrossRef]
  42. Bai, W.Q.; Zhang, K.Y.; Ding, X.M.; Bai, S.P.; Wang, J.P.; Peng, H.W.; Zeng, Q.F. High Dietary Energy Content Increases Inflammatory Markers after Lipopolysaccharide Challenge in Meat Ducks. Poult. Sci. 2019, 98, 164–171. [Google Scholar] [CrossRef] [PubMed]
  43. Schat, K.A.; Skinner, M.A. Avian Immunosuppressive Diseases and Immunoevasion. In Avian Immunology, 2nd ed.; Academic Press: Cambridge, MA, USA, 2013; pp. 275–297. [Google Scholar]
  44. Grond, K.; Sandercock, B.K.; Jumpponen, A.; Zeglin, L.H. The Avian Gut Microbiota: Community, Physiology and Function in Wild Birds. J. Avian Biol. 2018, 49, e01788. [Google Scholar] [CrossRef] [Green Version]
  45. Cho, H.J.; Zhao, P.; Kim, I.H. Effects of Emulsifier and Multi-Enzyme in Different Energy Densitydiet on Growth Performance, Blood Profiles, and Relative Organ Weight in Broiler Chickens. J. Agric. Sci. 2012, 4, 161. [Google Scholar] [CrossRef] [Green Version]
  46. Hussein, E.O.S.; Suliman, G.M.; Alowaimer, A.N.; Ahmed, S.H.; Abd El-Hack, M.E.; Taha, A.E.; Swelum, A.A. Growth, Carcass Characteristics, and Meat Quality of Broilers Fed a Low-Energy Diet Supplemented with a Multienzyme Preparation. Poult. Sci. 2020, 99, 1988–1994. [Google Scholar] [CrossRef] [PubMed]
  47. Adebowale, T. The Effect of Dietary High Energy Density and Carbohydrate Energy Ratio on Digestive Enzymes Activity, Nutrient Digestibility, Amino Acid Utilization and Intestinal Morphology of Weaned Piglets. J. Anim. Physiol. Anim. Nutr. 2019, 103, 1492–1502. [Google Scholar] [CrossRef] [PubMed]
  48. Murugesan, G.R.; Syed, B.; Haldar, S.; Pender, C. Phytogenic Feed Additives as an Alternative to Antibiotic Growth Promoters in Broiler Chickens. Front. Vet. Sci. 2015, 2, 21. [Google Scholar]
  49. Abdelqader, A.; Al-Fataftah, A.R. Effect of Dietary Butyric Acid on Performance, Intestinal Morphology, Microflora Composition and Intestinal Recovery of Heat-Stressed Broilers. Livest. Sci. 2016, 183, 78–83. [Google Scholar] [CrossRef]
  50. Broom, L.J.; Kogut, M.H. The Role of the Gut Microbiome in Shaping the Immune System of Chickens. Vet. Immunol. Immunopathol. 2018, 204, 44–51. [Google Scholar] [CrossRef]
  51. Zimmermann, P.; Curtis, N. The Effect of Antibiotics on the Composition of the Intestinal Microbiota—A Systematic Review. J Infect. 2019, 79, 471–489. [Google Scholar] [CrossRef]
  52. Dal Pont, G.C.; Farnell, M.; Farnell, Y.; Kogut, M.H. Dietary Factors as Triggers of Low-Grade Chronic Intestinal Inflammation in Poultry. Microorganisms 2020, 8, 139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Den Besten, G.; Van Eunen, K.; Groen, A.K.; Venema, K.; Reijngoud, D.J.; Bakker, B.M. The Role of Short-Chain Fatty Acids in the Interplay between Diet, Gut Microbiota, and Host Energy Metabolism. J. Lipid Res. 2013, 54, 2325–2340. [Google Scholar] [CrossRef] [Green Version]
  54. Angelakis, E.; Raoult, D. The Increase of Lactobacillus Species in the Gut Flora of Newborn Broiler Chicks and Ducks Is Associated with Weight Gain. PLoS ONE 2010, 5, e10463. [Google Scholar] [CrossRef]
  55. Louis, P.; Flint, H.J. Formation of Propionate and Butyrate by the Human Colonic Microbiota. Environ. Microbiol. 2017, 19, 29–41. [Google Scholar] [CrossRef] [Green Version]
  56. Morrison, D.J.; Preston, T. Formation of Short Chain Fatty Acids by the Gut Microbiota and Their Impact on Human Metabolism. Gut Microbes 2016, 7, 189–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Eck, A.; Rutten, N.B.M.M.; Singendonk, M.M.J.; Rijkers, G.T.; Savelkoul, P.H.M.; Meijssen, C.B.; Crijns, C.E.; Oudshoorn, J.H.; Budding, A.E.; Vlieger, A.M. Neonatal Microbiota Development and the Effect of Early Life Antibiotics Are Determined by Two Distinct Settler Types. PLoS ONE 2020, 15, e0228133. [Google Scholar] [CrossRef] [Green Version]
  58. Jandhyala, S.M.; Talukdar, R.; Subramanyam, C.; Vuyyuru, H.; Sasikala, M.; Reddy, D.N. Role of the Normal Gut Microbiota. World J. Gastroenterol. 2015, 21, 8787–8803. [Google Scholar] [CrossRef]
  59. Shanker, E.; Federle, M.J. Quorum Sensing Regulation of Competence and Bacteriocins in Streptococcus Pneumoniae and Mutans. Genes 2017, 8, 15. [Google Scholar] [CrossRef]
  60. Matsumoto-Nakano, M. Role of Streptococcus Mutans Surface Proteins for Biofilm Formation. Jpn. Dent. Sci. Rev. 2018, 54, 22–29. [Google Scholar] [CrossRef] [PubMed]
  61. Chen, Y.R.; Zheng, H.M.; Zhang, G.X.; Chen, F.L.; Chen, L.D.; Yang, Z.C. High Oscillospira Abundance Indicates Constipation and Low BMI in the Guangdong Gut Microbiome Project. Sci. Rep. 2020, 10, 9364. [Google Scholar] [CrossRef] [PubMed]
  62. García-Bayona, L.; Comstock, L.E. Bacterial Antagonism in Host-Associated Microbial Communities. Science 2018, 361, eaat2456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wang, C.Y.; Dawid, S. Mobilization of Bacteriocins during Competence in Streptococci. Trends Microbiol. 2018, 26, 389–391. [Google Scholar] [CrossRef]
  64. Shang, Y.; Kumar, S.; Oakley, B.; Kim, W.K. Chicken Gut Microbiota: Importance and Detection Technology. Front. Vet. Sci. 2018, 5, 254. [Google Scholar] [CrossRef]
  65. Rodríguez, J.M.; Avershina, E.; Ross, R.P.; Murphy, K.; Stanton, C.; Rudi, K.; Kober, O.I.; Marchesi, J.R.; Narbad, A.; Juge, N.; et al. The Composition of the Gut Microbiota throughout Life, with an Emphasis on Early Life. Microb. Ecol. Health Dis. 2015, 26, 26050. [Google Scholar] [CrossRef]
  66. Donaldson, E.E.; Stanley, D.; Hughes, R.J.; Moore, R.J. The Time-Course of Broiler Intestinal Microbiota Development after Administration of Cecal Contents to Incubating Eggs. PeerJ 2017, 5, e3587. [Google Scholar] [CrossRef] [Green Version]
  67. Jiang, W.; Wu, N.; Wang, X.; Chi, Y.; Zhang, Y.; Qiu, X.; Hu, Y.; Li, J.; Liu, Y. Dysbiosis Gut Microbiota Associated with Inflammation and Impaired Mucosal Immune Function in Intestine of Humans with Non-Alcoholic Fatty Liver Disease. Sci. Rep. 2015, 5, 8096. [Google Scholar] [CrossRef]
  68. Graziani, F.; Pujol, A.; Nicoletti, C.; Dou, S.; Maresca, M.; Giardina, T.; Fons, M.; Perrier, J. Ruminococcus Gnavus E1 Modulates Mucin Expression and Intestinal Glycosylation. J. Appl. Microbiol. 2016, 120, 1403–1417. [Google Scholar] [CrossRef] [Green Version]
  69. Song, Y.; Malmuthuge, N.; Steele, M.A.; Guan, L.L. Shift of Hindgut Microbiota and Microbial Short Chain Fatty Acids Profiles in Dairy Calves from Birth to Pre-Weaning. FEMS Microbiol. Ecol. 2018, 94, fix179. [Google Scholar] [CrossRef]
  70. Hollister, E.B.; Gao, C.; Versalovic, J. Compositional and Functional Features of the Gastrointestinal Microbiome and Their Effects on Human Health. Gastroenterology 2014, 146, 1449–1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Neal-McKinney, J.M.; Lu, X.; Duong, T.; Larson, C.L.; Call, D.R.; Shah, D.H.; Konkel, M.E. Production of Organic Acids by Probiotic Lactobacilli Can Be Used to Reduce Pathogen Load in Poultry. PLoS ONE 2012, 7, e43928. [Google Scholar] [CrossRef]
  72. Huang, P.; Zhang, Y.; Xiao, K.; Jiang, F.; Wang, H.; Tang, D.; Liu, D.; Liu, B.; Liu, Y.; He, X.; et al. The Chicken Gut Metagenome and the Modulatory Effects of Plant-Derived Benzylisoquinoline Alkaloids 06 Biological Sciences 0605 Microbiology. Microbiome 2018, 6, 211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Robinson, K.; Becker, S.; Xiao, Y.; Lyu, W.; Yang, Q.; Zhu, H.; Yang, H.; Zhao, J.; Zhang, G. Differential Impact of Subtherapeutic Antibiotics and Ionophores on Intestinal Microbiota of Broilers. Microorganisms 2019, 7, 282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Microorganisms 09 00787 g001 550
Figure 1. The total number of quality-filtered reads per sample. These reads reflect the total number of high-quality sequences that align with 16Sv4, clustered into OTUs, and were assigned taxonomic classification. Diets 1, 2, 6, and 7 represent normal-energy basal diet (NE-BAS), normal-energy diet with bacitracin methylene disalicylate (NE-BMD), high-energy basal diet (HE-BAS), and high-energy diet with bacitracin methylene disalicylate (HE-BMD), respectively.
Figure 1. The total number of quality-filtered reads per sample. These reads reflect the total number of high-quality sequences that align with 16Sv4, clustered into OTUs, and were assigned taxonomic classification. Diets 1, 2, 6, and 7 represent normal-energy basal diet (NE-BAS), normal-energy diet with bacitracin methylene disalicylate (NE-BMD), high-energy basal diet (HE-BAS), and high-energy diet with bacitracin methylene disalicylate (HE-BMD), respectively.
Microorganisms 09 00787 g001
Microorganisms 09 00787 g002a 550Microorganisms 09 00787 g002b 550
Figure 2. (a) Alpha analysis using the Shannon index showing the diversity of gut microbiota among the treatments (NE-BAS, NE-BMD, HE-BAS, and HE-BMD); the significance was tested with ANOVA and plot was generated using GraphPad prism. (b) The principal component analysis (PCA) plot, created using STAMP software, compares genus-level taxonomic profiles among different diets. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet with bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet with bacitracin methylene disalicylate. (c) Relative abundance of top-class in broiler chickens fed different diets. (d) Relative abundance of top phyla in broiler chickens fed different diets. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet containing bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet containing bacitracin methylene disalicylate. (e) Relative abundance of top genera in broiler chickens fed different diets (NE-BAS, NE-BMD, HE-BAS, and HE-BMD). Different colors represent different genera. (f) Significant differences in genera between the chickens fed basal (BAS) and bacitracin methylene disalicylate (BMD), assessed using STAMP software package (p < 0.05).
Figure 2. (a) Alpha analysis using the Shannon index showing the diversity of gut microbiota among the treatments (NE-BAS, NE-BMD, HE-BAS, and HE-BMD); the significance was tested with ANOVA and plot was generated using GraphPad prism. (b) The principal component analysis (PCA) plot, created using STAMP software, compares genus-level taxonomic profiles among different diets. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet with bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet with bacitracin methylene disalicylate. (c) Relative abundance of top-class in broiler chickens fed different diets. (d) Relative abundance of top phyla in broiler chickens fed different diets. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet containing bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet containing bacitracin methylene disalicylate. (e) Relative abundance of top genera in broiler chickens fed different diets (NE-BAS, NE-BMD, HE-BAS, and HE-BMD). Different colors represent different genera. (f) Significant differences in genera between the chickens fed basal (BAS) and bacitracin methylene disalicylate (BMD), assessed using STAMP software package (p < 0.05).
Microorganisms 09 00787 g002aMicroorganisms 09 00787 g002b
Microorganisms 09 00787 g003 550
Figure 3. Significant differences in microbes among different diets using STAMP software package (p-value < 0.05). (a) Significant differences in microbes between chickens fed NE-BAS and NE-BMD. (b) Significant differences in microbes between chickens fed NE-BAS and HE-BAS. (c) Significant differences in microbes between chickens fed NE-BAS and HE-BMD. (d) Significant differences in microbes between chickens fed NE-BMD and HE-BMD. (e) Significant differences in microbes between chickens fed NE-BMD and HE-BAS. (f) Significant differences in microbes between chickens fed HE-BAS and HE-BMD. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet containing bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet containing bacitracin methylene disalicylate.
Figure 3. Significant differences in microbes among different diets using STAMP software package (p-value < 0.05). (a) Significant differences in microbes between chickens fed NE-BAS and NE-BMD. (b) Significant differences in microbes between chickens fed NE-BAS and HE-BAS. (c) Significant differences in microbes between chickens fed NE-BAS and HE-BMD. (d) Significant differences in microbes between chickens fed NE-BMD and HE-BMD. (e) Significant differences in microbes between chickens fed NE-BMD and HE-BAS. (f) Significant differences in microbes between chickens fed HE-BAS and HE-BMD. NE-BAS = normal-energy basal diet; NE-BMD = normal-energy diet containing bacitracin methylene disalicylate; HE-BAS = high-energy basal diet; HE-BMD = high-energy diet containing bacitracin methylene disalicylate.
Microorganisms 09 00787 g003
Microorganisms 09 00787 g004a 550Microorganisms 09 00787 g004b 550
Figure 4. (a) The principal component analysis (PCA) plot, created using STAMP software, compares genus-level taxonomic profiles across growth stages (day 21 (D21), day 36 (D36), and day 43 (D41)). (b) Alpha analysis using the Shannon index, showing the diversity of gut microbiota at different growth stages (D21, D36, and D41). (c) Principal coordinate analysis (PCoA) plot obtained from beta analysis shows differences in bacterial community regardless of time. (d) LEfSe analysis result revealed the effect size of some significant taxa at different stages of growth. The plot was generated using the online LEfSe project. Red, green, and blue represent the enriched taxa in chickens on days 21, 36, and 43, respectively. (e) Significant differences in microbes between days 21 and 36 and between days 36 and 43 (q-value < 0.05), assessed using STAMP software package. (f) Relative abundance of top phyla during different stages of growth (days). (g) Relative abundance of genera during different stages of growth (days). Bar plots were generated using GraphPad Prism.
Figure 4. (a) The principal component analysis (PCA) plot, created using STAMP software, compares genus-level taxonomic profiles across growth stages (day 21 (D21), day 36 (D36), and day 43 (D41)). (b) Alpha analysis using the Shannon index, showing the diversity of gut microbiota at different growth stages (D21, D36, and D41). (c) Principal coordinate analysis (PCoA) plot obtained from beta analysis shows differences in bacterial community regardless of time. (d) LEfSe analysis result revealed the effect size of some significant taxa at different stages of growth. The plot was generated using the online LEfSe project. Red, green, and blue represent the enriched taxa in chickens on days 21, 36, and 43, respectively. (e) Significant differences in microbes between days 21 and 36 and between days 36 and 43 (q-value < 0.05), assessed using STAMP software package. (f) Relative abundance of top phyla during different stages of growth (days). (g) Relative abundance of genera during different stages of growth (days). Bar plots were generated using GraphPad Prism.
Microorganisms 09 00787 g004aMicroorganisms 09 00787 g004b
Table
Table 1. Ingredient, calculated, and analyzed compositions of the basal diets used in the study.
Table 1. Ingredient, calculated, and analyzed compositions of the basal diets used in the study.
StarterGrowerFinisher
Normal EnergyHigh EnergyNormal EnergyHigh EnergyNormal EnergyHigh Energy
Ingredient composition
Corn41.3342.6644.3239.6848.4844.11
Soybean meal (46.5% CP)40.1736.4336.4838.7031.5233.52
Wheat 10.0010.0010.0010.0010.0010.00
Vegetable Oil3.435.874.597.005.678.04
Limestone1.801.801.651.641.521.51
Dicalcium phosphate1.231.211.061.050.930.92
Pellet Binding Agent Y0.500.500.500.500.500.50
DL-Methionine Premix X0.610.630.530.560.490.52
Vitamin–Mineral Premix W,V0.500.500.500.500.500.50
Salt0.400.400.370.370.380.38
HCL Lys0.030.000.000.000.010.00
Calculated composition
Folic acid (ppm)2.202.202.202.202.202.20
Digestible Trp0.250.260.230.240.210.22
Digestible Thr0.870.900.820.840.740.76
Digestible Met + Cys0.950.980.870.900.800.83
Digestible Lys1.281.321.161.211.031.07
ME, kcal/kg300031003100320032003300
Crude protein23.023.821.522.219.520.1
Calcium0.960.960.870.870.790.79
Available P0.480.480.440.440.400.40
Sodium0.190.190.180.180.180.18
ME/CP130130144144164164
ME/Digestible Lys234423492672264531073084
ME/Met + Cys315831633563355640003976
ME/Digestible Thr344834443781381043244342
ME/Trp12,00011,92313,47813,33315,23815,000
Analyzed composition (%, except where otherwise stated)
ME, kcal/kg302231453109320332203317
Crude protein24.225.221.722.920.220.3
Crude fat6.118.036.418.697.046.70
Dry matter88.889.186.287.987.786.7
X Supplied/kg premix: DL-methionine, 0.5 kg; wheat middlings, 0.5 kg. Y Pel-stik: Uniscope, Inc., Johnstown, CO, USA. W Starter vitamin–mineral premix contained the following per kg of diet: 9750 IU vitamin A, 2000 IU vitamin D3, 25 IU vitamin E, 2.97 mg vitamin K, 7.6 mg riboflavin, 13.5 mg Dl Ca-pantothenate, 0.012 mg vitamin B12, 29.7 mg niacin, 1.0 mg folic acid, 801 mg choline, 0.3 mg biotin, 4.9 mg pyridoxine, 2.9 mg thiamine, 70.2 mg manganese, 80.0 mg zinc, 25 mg copper, 0.15 mg selenium, 50 mg ethoxyquin, 1543 mg wheat middlings, and 500 mg ground limestone. V Grower and Finisher vitamin–mineral premix contained the following per kg of diet: 9750 IU vitamin A, 2000 IU vitamin D3, 25 IU vitamin E, 2.97 mg vitamin K, 7.6 mg riboflavin, 13.5 mg Dl Ca-pantothenate, 0.012 mg vitamin B12, 29.7 mg niacin, 1.0 mg folic acid, 801 mg choline, 0.3 mg biotin, 4.9 mg pyridoxine, 2.9 mg thiamine, 70.2 mg manganese, 80.0 mg zinc, 25 mg copper, 0.15 mg selenium, 50 mg ethoxyquin, 1543 mg wheat middlings, and 500 mg ground limestone.
Table
Table 2. Effect of dietary energy density and bacitracin withdrawal on growth performance of broiler chickens.
Table 2. Effect of dietary energy density and bacitracin withdrawal on growth performance of broiler chickens.
Normal EnergyHigh EnergySEM 1p-Value
Basal BMD 2Basal BMD BMDEnergyBMD × Energy
Feed intake, g/bird
Days 0–71641641671715.070.6700.3630.724
Days 8–143543623503519.910.6360.4480.727
Days 15–24112910681288105997.80.1490.4470.398
Days 25–35183717951824175325.40.0350.2790.574
Days 35–42118614021198118783.80.2320.2350.186
Days 0–4246704791482645201370.5070.6790.130
Body weight, g
Day 043.643.343.243.20.290.6580.3900.658
Day 7186 a177 ab177 ab171 b3.700.1310.0100.892
Day 1450149250250613.50.8300.5730.611
Day 24131012901284128713.60.5170.2890.404
Day 35247124532457247723.00.9660.8290.433
Day 42323332473246329138.20.4470.4640.681
Body weight gain, g/bird
Days 0–7142 a138 ab134 ab127 b3.730.1510.0120.842
Days 8–1431731432633213.20.9290.3320.733
Days 15–2480779578178515.90.7890.2760.639
Days 25–35116111641173119015.90.5450.2280.675
Days 25–4276279478981422.00.2040.2970.902
Days 0–42319032043203324838.20.4460.4550.687
Feed conversion ratio
Days 0–71.15 b1.19 ab1.25 ab1.38 a0.060.1940.0170.609
Days 8–141.131.161.091.070.050.8710.1730.663
Days 15–241.411.341.671.350.130.1350.3600.414
Days 25–351.59 a1.54 a1.56 a1.47 b0.020.0030.0090.286
Days 25–421.551.801.521.460.120.6670.1170.220
Days 0–421.471.501.511.390.010.3350.3780.099
In a row, means assigned different lowercase letters are significantly different, p < 0.05 (Tukey’s procedure). 1 Standard error of the mean. 2 Bacitracin methylene disalicylate.
Table
Table 3. Effect of dietary energy density and BMD on organ weights (g/kg body weight).
Table 3. Effect of dietary energy density and BMD on organ weights (g/kg body weight).
ParameterNormal EnergyHigh EnergySEM 1p-Value
Basal BMD 2Basal BMDBMDEnergyBMD × Energy
Spleen0.740.860.830.760.070.7380.8900.148
Ceca5.704.554.614.600.510.2600.3140.264
Liver15.815.515.415.00.530.4840.4170.883
Bursa1.481.281.651.720.140.6080.0220.281
Heart4.805.204.804.990.210.1750.6090.622
1 Standard error of the mean. 2 Bacitracin methylene disalicylate.
Table
Table 4. Effect of dietary energy level and BMD on jejunal morphology (mm) of broiler chickens.
Table 4. Effect of dietary energy level and BMD on jejunal morphology (mm) of broiler chickens.
ParameterNormal EnergyHigh EnergySEM 1p-Value
Basal BMD 2Basal BMDBMDEnergyBMD × Energy
Villus height1.51 ab1.33 b1.51 ab1.63 a0.040.6300.0190.051
Villus width0.18 a0.19 a0.12 b0.19 a0.010.0070.5190.017
Crypt depth0.170.180.150.170.010.1790.5520.666
VH/CD 39.727.789.639.360.390.1210.0620.453
In a row, means assigned different lowercase letters are significantly different, p < 0.05 (Tukey’s procedure). 1 Standard error of the mean; 2 Bacitracin methylene disalicylate. 3 Villus height/crypt depth ratio.
Table
Table 5. Fifteen OTUs that were differentially abundant amongst treatments (NE-BAS, NE-BMD, HE-BAS, and HE-BMD) (FDR < 0.05).
Table 5. Fifteen OTUs that were differentially abundant amongst treatments (NE-BAS, NE-BMD, HE-BAS, and HE-BMD) (FDR < 0.05).
log2 Fold Changep-ValuePadjPhylumGenus
Otu0001310.276283.77 × 10−1043.76 × 10−101BacteroidotaBacteroidia_unclassified
Otu0000812.164691.10 × 10−525.50 × 10−50BacteroidotaAlistipes
Otu002106.0157422.53 × 10−348.40 × 10−32BacteroidotaBacteroidia_unclassified
Otu0001612.062191.26 × 10−293.15 × 10−27FirmicutesPeptostreptococcaceae_unclassified
Otu00081−5.505315.47 × 10−231.09 × 10−20FirmicutesFaecalibacterium
Otu0000412.521446.05 × 10−211.00 × 10−18BacteroidotaBacteroides
Otu003096.3011781.10 × 10−201.57 × 10−18FirmicutesRuminococcus
Otu002095.8071561.99 × 10−192.47 × 10−17BacteroidotaAlistipes
Otu00001−3.825544.49 × 10−194.96 × 10−17FirmicutesFaecalibacterium
Otu000396.4944052.03 × 10−162.02 × 10−14BacteroidotaBacteroides
Otu00280−5.76752.99 × 10−162.70 × 10−14FirmicutesFaecalibacterium
Otu0007110.031661.45 × 10−151.21 × 10−13FirmicutesBacillales_unclassified
Otu000798.4994491.78 × 10−151.36 × 10−13FirmicutesLachnospiraceae_unclassified
Otu002296.9410291.93 × 10−151.37 × 10−13FirmicutesFirmicutes_unclassified
Otu001503.5742642.30 × 10−151.53 × 10−13FirmicutesLachnospiraceae_unclassified
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adewole, D.; Akinyemi, F. Gut Microbiota Dynamics, Growth Performance, and Gut Morphology in Broiler Chickens Fed Diets Varying in Energy Density with or without Bacitracin Methylene Disalicylate (BMD). Microorganisms 2021, 9, 787. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040787

AMA Style

Adewole D, Akinyemi F. Gut Microbiota Dynamics, Growth Performance, and Gut Morphology in Broiler Chickens Fed Diets Varying in Energy Density with or without Bacitracin Methylene Disalicylate (BMD). Microorganisms. 2021; 9(4):787. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040787

Chicago/Turabian Style

Adewole, Deborah, and Fisayo Akinyemi. 2021. "Gut Microbiota Dynamics, Growth Performance, and Gut Morphology in Broiler Chickens Fed Diets Varying in Energy Density with or without Bacitracin Methylene Disalicylate (BMD)" Microorganisms 9, no. 4: 787. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040787

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop