Introduction
Prebiotics are typically fibrous compounds that pass undigested through the upper part of the gastrointestinal tract and stimulate the growth or activity of advantageous bacteria that colonize the bowel and contribute to the well-being of their host (Gibson and Roberfroid, 1995). Some of the most widely used prebiotics in the poultry industry are fructo-oligosaccharides, mannan-oligosaccharides, galacto-oligosaccharides and beta-glucans (Huyghebaert et al., 2011).
Selective fermentation of some prebiotics has been shown to induce changes in the composition and/or activity of the gastrointestinal microbiota, improving the health of the host (Gibson et al., 2004). Zhenping et al. (2013) showed increased growth performance, enhanced endocrine metabolism, and improved immune function in broiler chickens after in-feed supplementation with xylo-oligosaccharides prebiotics. Moreover, changes in enteric bacteria in the cecum (Spring et al., 2000) and improved intestinal morphology have been observed in broilers fed with dietary mannan-oligosaccharide (Baurhoo et al., 2009).
Because the microbiota can be modified (Khoruts et al., 2010; Borody and Khoruts, 2012) it constitutes an attractive target for therapeutic manipulation. However, successful outcome of such manipulations require a better understanding of the interactions between the host and its microbiota (Hamilton et al., 2013; Van Nood et al., 2013).
The composition and diversity of chicken intestinal microbiota were previously investigated using cultivation-based methodologies (Fernandez et al., 2002; Cross et al., 2007). However, the use of DNA-based molecular biology techniques, such as metagenomics and new generation DNA sequencing, allowed new opportunities to characterize uncultivable members of intestinal microbiota (Gong et al., 2002) shedding light on the composition and temporal spatial location of the microbial population in broiler's intestine. Nevertheless, the knowledge of the structure, interactions and functions of the intestinal microbiota is still limited and fragmented (Oakley and Kogut, 2016).
In this work the effect of dietary supplement with prebiotics derived from yeast wall (mannan-oligosaccharide and nucleotide) on the cecal microbiota of broilers was evaluated by the massive parallel sequencing of the 16S rRNA gene.
Materials and methods
Experimental design and sampling
Thirty female newborn chicks (Hubbart®) were placed on wood shavings litter for 35 days. Water and feed were given ad libitum. At arrival, chicks were randomly divided in groups for three different types of treatment (10 birds per treatment): treatment without additives in the feed (Neg), treatment with the prebiotics mannan-oligosaccharide (Mos), or nucleotides (Nuc). 200 g/ton of each prebiotic were incorporated in the ration from the first day until the 35th day. Briefly, the prebiotics used in this work were extracted from biomass derived from the sugar and beer industry. These by-products contain abundant in Saccharomyces cerevisiae yeast. The raw-matter was treated by industrial autolysis, obtaining two fractions. The first, composed mainly for yeast wall, was rich in mannan-oligosaccharides. The second fraction contained soluble yeast extract and was abundant in nucleotides (Chaud and Sgarbieri, 2006). All animal procedures were approved by the Animal Experimentation Ethics Committee of the Federal University of Paraná (authorization CEUA-Bio UFPR 898/2015). At 14 and 35 days of age, the chicks' body weights were recorded and samples were collected for DNA extraction, purification and sequencing. Immediately after euthanasia, the abdominal cavity was exposed and the cecum was dissected from the other intestinal sections. The cecum from each bird was cut open, and the contents were collected in a sterile 2-ml tube, stored on ice and later frozen and stored at −80°C until use. In total, 24 samples for cecal contents were collected (four chickens per treatment at 14 and 35 days of age).
DNA extraction, 16S rRNA gene amplification and sequencing
Genomic DNA from each sample was isolated from 200 mg of cecal luminal content using PowerFecal® DNA Isolation Kit (MO BIO laboratories, Inc.). The variable V4 region of 16S rRNA gene was amplified using the universal primers 515F and 806R (Caporaso et al., 2011) and KlenTaq Master Mix (Sigma). The PCR conditions used were 94°C for 3 min; 18 cycles of 94°C for 45 s, 50°C for 30 s and 68°C for 60 s; followed by 72°C for 10 min. The amplicons were quantified with Qubit using HS dsDNA kit (Invitrogen), diluted to 500pM and pooled. Then, 16pM of pooled DNA were sequenced using MiSeq reagent 500V2. Sequencing was performed using an Illumina MiSeq® sequencer (Illumina) obtaining paired-end reads of 250 bp as described (Caporaso et al., 2011). The dataset were submitted in NCBI site under the BioSample accession code SAMN07211773 (http://www.ncbi.nlm.nih.gov/biosample/7211773).
Diversity analysis
Sequencing data were analyzed with the QIIME pipeline (Caporaso et al., 2010). Since the Illumina output ranged from approximately 40,000 to 172,499 reads, the read number was re-sampled to 41,800 reads per sample, allowing for the diversity comparisons. Sequences were quality filtered and identified at the phyla and genus levels using the open-reference OTU method implemented in QIIME and the SILVA database (123 release) (Yilmaz et al., 2013; Rideout et al., 2014). Basic diversity analysis (OTU number, Unifrac-distances and cluster by Neighbor-Joining method) was conducted using QIIME, and the OTU table exported to R for further analyses.
Statistical analysis
Statistical analyses were performed using the Stats package included in RStudio software (RStudio-Team, 2015). First, data was explored using multivariate tools. Weighted unifrac-distances previously obtained from QIIME were used to sample ordination using PCoA and cluster analyzes. Since the samples showed a non-parametric distribution by the Shapiro Wilk test, the genera previously identified were compared according to the prebiotic treatments using the Kruskal Walli's test. Only statically significant results were reported (P < 0.05).
Results
The use of the prebiotics Mos and Nuc resulted in no changes in the body weight of the birds in 14 or 35 days (not shown). However, they did cause several changes in the gut microbiota. Using multivariate analysis, small changes in the microbiota were observed at 14 days, but differences were clearer after 35 days (Table 1). At 14 days, Firmicutes was the most predominant phylum in all treatments and, on average, accounted for 84.5% of identified sequences, followed by Proteobacteria (7.9%) and Tenericutes (4.3%).
Table 1
| Genera | Treatment | 14 day | K.W. | 35 day | K.W. | Time |
|---|---|---|---|---|---|---|
| Faecalibacterium | Neg | 21.5 ± 3.8 | A | 6.1 ± 1.0 | A | * |
| Mos | 25.7 ± 2.9 | A | 1.7 ± 0.5 | B | * | |
| Nuc | 13.6 ± 3.6 | B | 1.3 ± 0.1 | B | * | |
| Bacteroides | Neg | 0 | 10.9 ± 1.4 | B | * | |
| Mos | 0 | 17.7 ± 4.5 | B | * | ||
| Nuc | 0 | 35.9 ± 1.9 | A | * | ||
| Lactobacillus | Neg | 8.4 ± 2.2 | A | 3.4 ± 0.9 | B | * |
| Mos | 8.6 ± 2.0 | A | 3.5 ± 1.2 | B | * | |
| Nuc | 2.9 ± 0.7 | B | 6.8 ± 1.1 | A | * | |
| Streptococcus | Neg | 9.7 ± 1.2 | A | 3.3 ± 1.4 | B | * |
| Mos | 4.0 ± 3.8 | B | 5.9 ± 0.4 | A | * | |
| Nuc | 1.9 ± 1.7 | B | 5.6 ± 0.8 | A | * | |
| Lachnoclostridium | Neg | 5.4 ± 1.2 | 3.1 ± 0.5 | A | * | |
| Mos | 4.3 ± 0.3 | 1.6 ± 0.3 | B | * | ||
| Nuc | 5.2 ± 1.2 | 3.0 ± 0.1 | A | * | ||
| Anaerotruncus | Neg | 1.0 ± 0.2 | B | 6.3 ± 0.3 | A | * |
| Mos | 1.7 ± 0.3 | A | 3.5 ± 1.0 | B | * | |
| Nuc | 2.2 ± 0.4 | A | 5.9 ± 0.3 | A | * | |
| Ruminiclostridium 5 | Neg | 3.2 ± 0.3 | C | 2.4 ± 0.1 | A | |
| Mos | 5.1 ± 0.8 | B | 0.9 ± 0.3 | B | * | |
| Nuc | 6.5 ± 1.6 | A | 0.8 ± 0.1 | B | * | |
| Subdoligranulum | Neg | 7.2 ± 2.3 | A | 1.6 ± 0.2 | A | * |
| Mos | 1.5 ± 0.4 | B | 0.4 ± 0.1 | B | ||
| Nuc | 7.5 ± 2.9 | A | 0.1 ± 0.0 | B | * | |
| Parabacteroides | Neg | 0 | 0 | A | ||
| Mos | 0 | 11.4 ± 4.4 | B | * | ||
| Nuc | 0 | 0 | B | |||
| Escherichia-Shigella | Neg | 4.2 ± 1.5 | A | 0.4 ± 0.1 | A | * |
| Mos | 3.4 ± 2.0 | A | 0.2 ± 0.1 | B | * | |
| Nuc | 2.7 ± 2.5 | B | 0.1 ± 0.0 | B | * | |
| Alistipes | Neg | 0 | 4.9 ± 0.4 | A | * | |
| Mos | 0 | 4.1 ± 1.2 | A | * | ||
| Nuc | 0 | 0.9 ± 0.2 | B | * |
Most prevalent bacterial genera in broiler's cecum microbiota.
Data is shown by treatment (%, mean and standard error, n = 4). The K.W. columns show the result of statistical comparison between treatments (Neg, Mos, and Nuc) for each genus by the Kruskal Wallis test. Same uppercase letters indicate no statistical difference in genus abundance between treatments. The column Time shows statistical comparison between the two time points by the Tukey test. Asterisk (*) indicates significant difference between sampling times (P < 0.05). Control group (Neg), mannan-oligosaccharide (Mos) and nucleotide (Nuc).
After 35 days, Firmicutes remained the most predominant phylum in cecum in all treatments and, on average, accounted for 61.4% of all the bacterial sequences, followed by Bacteroidetes (29%), Proteobacteria (4.5%), and Tenericutes (4.2%). Interestingly, the phylum Bacteroidetes was only detected at 35 days (Supplementary Material 1). The prevalence of Firmicutes at 14 days might be explained by the higher need of butyrate by the young chicken growing intestine (Van der Wielen et al., 2000). Since in adult chickens maximal butyrate production is no longer required, digestion of complex polysaccharides by representatives of Bacteroidetes results in the production of both propionate and butyrate, which might be a more advantageous balance of nutrients (Polansky et al., 2016). In agreement with this assumption, Polansky et al. (2016) observed a gradual increase of Bacteroidetes counts in the cecum from the third week of life of chickens.
By further dissecting the 16S ribosomal sequences, a total of 665 OTUs corresponding to bacterial genera were detected, and were identified using the SILVA database (release 123) (Supplementary Material 2). Among all identified bacterial genera in the cecal microbiota, the most abundant in all treatments were Faecalibacterium, Bacteroides, Lactobacillus, Streptococcus, Lachnoclostridium, Anaerotruncus, Ruminiclostridium 5, Subdoligranulum, Parabacteroides, Escherichia-Shigella, and Alistipes (Supplementary Material 2). At 14 days, the most prevalent genera in all groups were Faecalibacterium accounting for, on average, 20.3% of all cecal bacterial sequences. The second most prevalent genera was Lactobacillus with 6.65% of bacterial sequences, followed by Subdoligranulum (5.42%) (Table 1). At 35 days, Bacteroides was the most predominant genera in all treatments, accounting for, on average, 21.5% of the cecal sequences. The second most prevalent was Anaerotruncus (5.59%), followed by Streptococcus (4.93%). Besides, at 35 days, the frequencies of 130 genera significantly differed between treatments. Over time an increase in the count of sequences was observed in Bacteroides, Anaerotruncus, and Alistipes genera. Interestingly, the Parabacteroides genus increased at 35 days only in the Mos treated group, while the genera Faecalibacterium, Lachnoclostridium, Ruminiclostridium 5, and Subdoligranulum decreased at 35 days (Supplementary Material 3).
Moreover, Principal Component Analysis of OTUs representation clearly showed the age explained most of the differences (57.8%), followed by treatment (17.4%) (Figure 1). This result was also observed in cluster analysis of microbial community by the Neighbor-Joining method (Supplementary Material 4).
Figure 1
The use of Nuc as prebiotics significantly increased the count of sequences of Lactobacillus and significantly decreased the count of sequences of Escherichia-Shigella at 35 day. On the other hand, the use of prebiotic Mos significantly increased the count of sequences of Butyricimonas and Roseburia genera at 35 day (Supplementary Material 6). The supplementation with both prebiotics significantly decreased the Chao1 index at 35 days, when compared to the negative group. This index is a species richness estimator, which estimates the total number of species present in a community (Supplementary Material 5).
Conclusion
In conclusion the composition of cecal microbiota is not constant and develops over time in chickens. Moreover, prebiotic supplementation significantly affected the microbial community structure in the cecum in broiler chickens with most significant shifts at 35 days. Although the prebiotics Mos and Nuc caused a decrease the richness and the diversity in the cecum, the prevalence of beneficial bacteria increased in both treatments. Nuc increased the count of sequences of Lactobacillus and decreased the count of sequences of Escherichia-Shigella, while Mos significantly increased the count of sequences of Butyricimonas and Roseburia. The Lactobacillus, Butyricimonas, and Roseburia genera are known for producing short-chain fatty acids (SCFAs) and have shown beneficial effects on the host's development and health (Pryde et al., 2002; Farnworth, 2005; Jhangi et al., 2014; Loman and Tappenden, 2016).
Statements
Author contributions
MN provided the prebiotics. DM and DL analyzed the data and performed the numerical and statistical analyses. DM and ES designed the work and wrote the manuscript. DL and LMC processed the sequence raw data. CS and LFC designed the protocol. DM performed DA extraction from samples. EB performed amplicons library preparation and sequencing. FP and ES supervised the work.
Acknowledgments
To CAPES/PROEX by the doctoral scholarship.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2017.00153/full#supplementary-material
References
1
BaurhooB.FerketP. R.ZhaoX. (2009). Effects of diets containing different concentrations of mannanoligosaccharide or antibiotics on growth performance, intestinal development, cecal and litter microbial populations, and carcass parameters of broilers. Poul. Sci.88, 2262–2272. 10.3382/ps.2008-00562
2
BorodyT. J.KhorutsA. (2012). Fecal microbiota transplantation and emerging applications. Nat. Rev. Gastroenterol. Hepatol.9, 88–96. 10.1038/nrgastro.2011.244
3
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. 10.1038/nmeth.f.303
4
CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.LozuponeC. A.TurnbaughP. J.et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A.108(Suppl. 1), 4516–4522. 10.1073/pnas.1000080107
5
ChaudS. G.SgarbieriV. C. (2006). Propriedades funcionais (tecnológicas) da parede celular de leveduras da fermentação alcoólica e das frações glicana, manana e glicoproteína. Food Sci. Technol.26, 369–379. 10.1590/S0101-20612006000200020
6
CrossD. E.McDevittR. M.HillmanK.AcamovicT. (2007). The effect of herbs and their associated essential oils on performance, dietary digestibility and gut microflora in chickens from 7 to 28 days of age. Br. Poul. Sci.48, 496–506. 10.1080/00071660701463221
7
FarnworthE. R. (2005). The Beneficial health effects of fermented foods-potential probiotics around the world. J. Nutraceut. Funct. Med. Foods4, 93–117. 10.1300/J133v04n03_07
8
FernandezF.HintonM.GilsB. V. (2002). Dietary mannan-oligosaccharides and their effect on chicken caecal microflora in relation to Salmonella Enteritidis colonization. Avian Pathol.31, 49–58. 10.1080/03079450120106000
9
GibsonG. R.ProbertH. M.LooJ. V.RastallR. A.RoberfroidM. B. (2004). Dietary modulation of the human colonic microbiota: updating the concept of prebiotics. Nut. Res. Rev.17, 259–275. 10.1079/NRR200479
10
GibsonG. R.RoberfroidM. B. (1995). Dietary modulation of the human colonic microbiota: introducing the concept of prebiotics. J. Nut.125, 1401–1412.
11
GongJ.ForsterR. J.YuH.ChambersJ. R.WheatcroftR.SabourP. M.et al. (2002). Molecular analysis of bacterial populations in the ileum of broiler chickens and comparison with bacteria in the cecum. FEMS Microbiol. Ecol.41, 171–179. 10.1111/j.1574-6941.2002.tb00978.x
12
HamiltonM. J.WeingardenA. R.UnnoT.KhorutsA.SadowskyM. J. (2013). High-throughput DNA sequence analysis reveals stable engraftment of gut microbiota following transplantation of previously frozen fecal bacteria. Gut Microbes4, 125–135. 10.4161/gmic.23571
13
HuyghebaertG.DucatelleR.ImmerseelF. V. (2011). An update on alternatives to antimicrobial growth promoters for broilers. Vet. J.187, 182–188. 10.1016/j.tvjl.2010.03.003
14
JhangiS.GandhiR.GlanzB.CookS.NejadP.WardD.et al. (2014). Increased archaea species and changes with therapy in gut microbiome of multiple sclerosis subjects (S24.001). Neurology82(10 Suppl.). Available online at: http://www.neurology.org/content/82/10_Supplement/S24.001.short
15
KhorutsA.DicksvedJ.JanssonJ. K.SadowskyM. J. (2010). Changes in the composition of the human fecal microbiome after bacteriotherapy for recurrent clostridium difficile-associated diarrhea. J. Clin. Gastroenterol.44, 354–360. 10.1097/MCG.0b013e3181c87e02
16
LomanB. R.TappendenK. A. (2016). Prebiotic short-chain fructooligosaccharides (scFOS) increases abundance of the butyrate producing microbial community differentially when administered with or without probiotic lactobacillus rhamnosus GG (LGG) in piglets with short-bowel syndrome (SBS). FASEB J.30(1 Suppl.), 683–685. Available online at: http://www.fasebj.org/content/30/1_Supplement/683.5.short
17
OakleyB. B.KogutM. H. (2016). Spatial and temporal changes in the broiler chicken cecal and fecal microbiomes and correlations of bacterial taxa with cytokine gene expression. Front. Vet. Sci.3:11. 10.3389/fvets.2016.00011
18
PolanskyO.SekelovaZ.FaldynovaM.SebkovaA.SisakF.RychlikI. (2016). Important metabolic pathways and biological processes expressed by chicken cecal microbiota. Appl. Environ. Microbiol.82, 1569–1576. 10.1128/AEM.03473-15
19
PrydeS. E.DuncanS. H.HoldG. L.StewartC. S.FlintH. J. (2002). The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett.217, 133–139. 10.1111/j.1574-6968.2002.tb11467.x
20
RStudio-Team (2015). RStudio: Integrated Development for R.Boston, MA: RStudio Team. Available online at: www.rstudio.com
21
RideoutJ. R.HeY.Navas-MolinaJ. A.WaltersW. A.UrsellL. K.GibbonsS. M.et al. (2014). Subsampled open-reference clustering creates consistent, comprehensive OTU definitions and scales to billions of sequences. PeerJ2:e545. 10.7717/peerj.545
22
SpringP.WenkC.DawsonK. A.NewmanK. E. (2000). The effects of dietary mannaoligosaccharides on cecal parameters and the concentrations of enteric bacteria in the ceca of salmonella-challenged broiler chicks. Poul. Sci.79, 205–211. 10.1093/ps/79.2.205
23
Van der WielenP. W. J. J.BiesterveldS.NotermansS.HofstraH.UrlingsB. A. P.van KnapenF. (2000). Role of volatile fatty acids in development of the cecal microflora in broiler chickens during growth. Appl. Environ. Microbiol.66, 2536–2540. 10.1128/AEM.66.6.2536-2540.2000
24
Van NoodE.VriezeA.NieuwdorpM.FuentesS.ZoetendalE. G.de VosW. M.et al. (2013). Duodenal infusion of donor feces for recurrent clostridium difficile. N. Eng. J. Med.368, 407–415. 10.1056/NEJMoa1205037
25
YilmazP.ParfreyL. W.YarzaP.GerkenJ.PruesseE.QuastC.et al. (2013). The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res.42, D643–D648. 10.1093/nar/gkt1209
26
ZhenpingS.WentingL.RuikuiY.JiaL.HonghongL.WeiS.et al. (2013). Effect of a straw-derived xylooligosaccharide on broiler growth performance, endocrine metabolism, and immune response. Can. J. Vet. Res.77, 105–109.
Summary
Keywords
mannan-oligosaccharides, chickens, gut microbiota, 16S rRNA, lactobacillus
Citation
Mesa D, Lammel DR, Balsanelli E, Sena C, Noseda MD, Caron LF, Cruz LM, Pedrosa FO and Souza EM (2017) Cecal Microbiota in Broilers Fed with Prebiotics. Front. Genet. 8:153. doi: 10.3389/fgene.2017.00153
Received
12 June 2017
Accepted
29 September 2017
Published
17 October 2017
Volume
8 - 2017
Edited by
Peter Dovc, University of Ljubljana, Slovenia
Reviewed by
Kieran G. Meade, Teagasc, The Irish Agriculture and Food Development Authority, Ireland; Shaojun Liu, Hunan Normal University, China
Updates
Copyright
© 2017 Mesa, Lammel, Balsanelli, Sena, Noseda, Caron, Cruz, Pedrosa and Souza.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dany Mesa dmesaf@unal.edu.co
This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.