Front. Vet. Sci.Frontiers in Veterinary ScienceFront. Vet. Sci.2297-1769Frontiers Media S.A.10.3389/fvets.2019.00114Veterinary ScienceOriginal ResearchDifferential Effects of Bacitracin Methylene Disalicylate (BMD) on the Distal Colon and Cecal Microbiota of Young Broiler ChickensProctorAlexandraPhillipsGregory J.*Department of Veterinary Microbiology and Preventative Medicine, College of Veterinary Medicine, Iowa State University, Ames, IA, United States
Edited by: Timothy J. Johnson, University of Minnesota Twin Cities, United States
Reviewed by: Guillermo Tellez, University of Arkansas, United States; Sherry Layton, Vetanco, Argentina
*Correspondence: Gregory J. Phillips gregory@iastate.edu
This article was submitted to Veterinary Infectious Diseases, a section of the journal Frontiers in Veterinary Science
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) and the copyright owner(s) 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.
Antibiotics have been used extensively for growth promotion in poultry, along with other food production animals, as well as therapeutically to treat infectious diseases. However, with concerns over selection for drug antibiotic resistant bacteria the practice of using subtherapeutic doses of antibiotics is under increased scrutiny. Consequently, we assessed the impact of the commonly used antibiotic bacitracin methylene disalicylate (BMD) on the gastrointestinal microbiota of chickens. For this we administered therapeutic doses of BMD as a feed additive and 16s rRNA gene amplicon sequencing to measure changes in taxonomic abundance on the distal colon and cecal microbiota of young broiler chickens. While BMD treatment was found to impact the abundance of selected taxa and overall beta diversity, significant changes were, in general, limited to the colon of the treated birds. Selected taxa at the phylum, class, and genus levels that were most impacted were identified. The composition of the cecum remained relatively stable in BMD-treated animals. As poultry production practices seek alternatives to growth promoting antibiotic feed additives, manipulation of the gastrointestinal microbiota holds promise. These results suggest that targeting the cecum may offer a means to promote changes to the microbiota that maximize the benefits for the hosts.
Along with therapeutic use, antibiotics are well-established for their ability to promote growth through improved weight gain and feed efficiency in livestock (1), including in broiler chickens (2). In poultry, bacitracin methylene disalicylate (BMD) is commonly used for growth promotion (3). Compared to other growth promoting antibiotics, such as virginiamycin, BMD is a relatively narrow spectrum antibiotic that targets primarily gram-positive bacteria, including Streptococci, Staphylococci, Clostridia, Fusobacterium, and Actinomyces. BMD interferes with protein synthesis and cell wall production and induces cell lysis in these microorganisms (4, 5). The antibiotic is not well-absorbed by the intestine and therefore primarily acts on bacteria in the gastrointestinal (GI) tract of the animals through delivery as a feed additive (6).
Feeding chickens low doses of BMD benefited the birds, including increased villus height throughout the small intestine and improved digestion of dietary components that correlated with increased body weight and feed consumption (7).
BMD is also used as to treat and prevent necrotic enteritis caused by Clostridium perfringens (8, 9), which is a cause of significant economic loss in the poultry industry (10). While antibiotic growth promoters make important contributions to the overall efficiency of livestock production, they are also not without their concern as sub-therapeutic doses used are also associated with selection and spread of drug resistant bacterial pathogens (11, 12). Concerns over widespread use of antibiotics in agriculture has prompted a ban on their use in the European Union with increased scrutiny for their use in the United States (13).
Given the importance of poultry for human nutrition and the food animal industry world-wide, emphasis has been placed on characterizing the chicken microbiome as a means to improve our understanding of antibiotic growth promotion and to identify alternative strategies that do not select for drug resistant bacteria (14–18). Toward this, numerous studies, representing a variety of methods, have assessed the impact of antibiotic treatment on the microbiota of poultry (7, 9, 19–32). In general, these studies have shown that growth-promoting antibiotics can have significant effects on the structure and function of the microbiota colonizing the GI tract. As to be expected with a list of wide-ranging studies, there are few bacterial taxa that are consistently altered by antibiotics that can explain their growth promoting activities since the composition and activity of the chicken microbiome is highly dependent on environmental conditions, feed composition and method of assessment of the microbioal communities. Interestingly, however, chicken microbiota studies have revealed that BMD, along with other growth-promoting antibiotics, can deplete species of Lactobacillus, as well as other probiotic species (32, 33). This observation has led to the suggestion that a reduction in bile-salt hydrolase activity encoded by many of these bacteria may contribute to growth promotion by reversing the negative effects on fat metabolism of these enzymes (34). Clearly additional studies are needed to better understand how changes to the microbiota by low-dose antibiotics contribute to animal growth enhancement.
To further out understanding of how the growth promoting feed additive BMD impacts the chicken microbiota, we have focused on distinguishing between the effects of the antibiotic on cecal vs. colon bacterial populations. These two compartments of the chicken digestive tract are colonized with distinct microbial communities (35, 36). Also, while metabolism and adsorption of macronutrients occurs primarily in the colon, fermentation of complex polysaccharides occurs primarily in the cecum (14, 15, 36). Because of these spatial and functional differences, we sought to determine the extent to which BMD impact the microbiota of the distinct compartments of the GI tract. For this, we conducted 16s rRNA gene amplicon taxonomic profiling of the microbiota of the distal colon and cecum from young broiler birds using therapeutic doses (8, 37) of the antibiotic to accentuate differences in microbial composition in the GI tract.
Materials and MethodsAnimal Model and Housing
This study was carried out with the approval of the Iowa State University Institutional Animal Care and Usage Committee under protocol number 6-11-7167-G. The design followed a necrotic enteritis model, however, no pathogens were administered to the chickens and only antimicrobial feed additives were added to the experimental group.
Approximately 30 day old jumbo Cornish/Rock broiler chicks were obtained from Welp Hatchery (Bancroft, IA) and housed in pens created by tying two 32″ × 8′ × 1/8″ (81 cm × 2.45 m × 0.3 cm) pegboards together to form a circle. This circular pen was divided into three equally sized areas with similar pegboard material. Each pen was bedded using ~3″ (7.5 cm) of wood shavings. Heat lamps were made available for each pen. One two-gallon, galvanized waterer and one galvanized metal feeder were supplied to each pen.
Groups of 15 birds were housed in the pens described. On days 1–7, each group received 1 kg of a low-protein chick starter (LPF) once a day. On days 8–10, each group was given 1 kg of a high protein feed (HPF) once a day. On days 11–18, the control group remained on the same HPF feed while the challenge group received the HPF supplemented with BMD (200 g/ton). On these days, each group was fed 1 kg HPF with or without BMD twice a day. On day 19, the chicks were euthanized and samples were collected. Distal colon and cecal contents were collected and stored at −80°C until total DNA was isolated.
DNA Isolation
Total genomic DNA was extracted using the PowerSoil DNA Isolation Kit (MoBio, Carlsbad, CA). The manufacturer's protocol was followed with the exception of the initial vortexing step, which was extended to 20 min to thoroughly homogenize the samples. Purified genomic DNA extracts were quantified using a Quibt 2.0 Fluorometer (Life Technologies, Carlsbad, CA), and DNA stored at −20°C in the provided 10 mM Tris buffer.
Sequencing and Analysis
PCR amplification of the variable regions 3–5 of the 16S rRNA gene was done using region specific primers 357F (CTCCTACGGGAGGCAGCAG)−926R (CCGTCAATTCMTTTRAGTTT). Amplicon sequencing was performed at the Research and Testing Laboratory (Lubbock, TX) using the Roche 454 Titanium platform.
The resulting DNA sequences were analyzed using QIIME (Quantitative Insights into Microbial Ecology) (38). The reads were first demultiplexed and binned per sample. The reads were also quality filtered during this step to remove poor quality sequences using default quality filtering values. Denoising of the 454 reads was performed using Denoiser (39). Chimeric sequences were identified using USEARCH and removed (40, 41). The remaining sequences were clustered into OTUs at 97% similarity using USEARCH and the open reference OTU picking strategy in QIIME. Sequences were aligned to the Greengenes (13_5) rRNA sequence core reference database using PyNAST (42, 43). Taxonomic assignments were made using the RDP Classifier 2.2 at a 97% similarity to the Greengenes reference database (44). Phylogenetic trees were built using FastTree 2.1.3 (45). Alpha and beta diversity, analysis of similarity (ANOSIM) and Adonis tests were generated using QIIME. PCoA plots were generated using Emperor in QIIME (46). Mann–Whitney U-tests were performed on taxonomic summaries using a custom R script (R Project) developed at the Institute of Genome Sciences at the University of Maryland-Baltimore. Sequence reads have been submitted to NCBI's short read archives.
ResultsAlpha and Beta Diversity
A total of 79,670 sequences were analyzed using QIIME. After filtering based on quality scores, 73,529 sequences corresponded to 619 OTUs with an average of 2,298 ± 1,335 sequences per sample.
Figure 1 summarizes the alpha diversity measurements used to determine the extent to which BMD altered the composition of the microbiota. This included measurements of observed OTUs (1A), Faith's phylogenetic diversity (1B), Shannon (1C), and Simpson (1D) alpha diversity indices to compare treatment groups. Non-parametric two-sample t-tests were used to identify significant differences among the groups. While significant differences in OTUs were not observed among the groups by the Simpson and Shannon diversity metrics (data not shown), differences were noted by Faith's phylogenetic diversity (Table 1).
Comparison of the composition of the microbiota in treatment groups. (A) Observed OTUs. (B) Faith's phylogenetic diversity. (C) Shannon alpha diversity index. (D) Simpson alpha diversity index. Treatment groups include: BMD treated, cecum (red line); BMD treated, distal colon (blue line); Control, cecum (orange line); Control, distal colon (green line).
Significant differences among the treatment and control groups as assessed by Faith's phylogenetic diversity non-parametric two-sample t-tests.
Group 1
Group 2
Group 1 mean
Group 1 std
Group 2 mean
Group 2 std
t statistic
p-value
Control-cecum
BMD-cecum
5.48
1.13
7.32
2.41
−1.80
0.66
BMD-cecum
BMD-colon
7.33
2.41
11.86
1.76
−3.89
0.04
Control-cecum
Control-colon
5.48
1.13
5.77
1.32
−0.43
1
BMD-colon
Control-colon
11.86
1.76
5.76
1.32
7.30
0.01
Within the BMD treated group, the distal colon showed higher phylogenetic diversity compared to the cecal samples (p = 0.036). The phylogenetic diversity of the colon of the BMD treated birds was also greater than the same samples from the control group (p = 0.006). This could indicate the BMD treatment is causing a few OTUs to become depleted hence allowing for other unique OTUs to establish or increase in their proportion within the community.
Figure 2 shows the weighted (abundance considered) Unifrac PcoA beta diversity plots of the treatment groups (Figure 2A), while Figure 2B shows the unweighted (abundance independent) plots. In the weighted Unifrac PCoA analysis, cecal samples from both treated and control groups showed greater similarity than samples from the colon. Conversely, the colon samples of the control and BMD groups showed greater variability, with less distinct clustering. Also, the control samples from the colon clustered closer to each other and also closer to the cecal samples than the BMD colon samples.
The unweighted Unifrac PCoA plots showed that treatment influenced the dissimilarity of the samples more than GI tract location as samples in each group were not clustered as tightly as the weighted PCoA. Analysis of similarity (ANOSIM) and Adonis tests were performed on the weighted and unweighted Unifrac distances obtained from the beta diversity workflow in QIIME. The ANOSIM test based on both treatment and GI location resulted in a p-value of 0.001 and a test statistic of 0.379 for weighted and a p-value of 0.001 and test statistic of 0.556 for unweighted Unifrac analysis. These metrics indicated that the grouping based on the variables of treatment and GI location is weak (i.e., the differences can be explained by randomness). Adonis tests also based on both treatment and GI location resulted in a p-value of 0.001 and an R2 value of 0.441 for weighted and a p-value of 0.001 and R2 value of 0.404 for unweigthed Unifrac analysis.
Relative Abundance
Differences in relative abundance among all the treatment groups at different taxonomic levels were assessed using Mann–Whitney U-tests. Figure 3 shows the taxonomic summary for each group at the phylum level and Table 2 shows the p-values for the Mann–Whitney U-tests for phylum level differences. As evident, the dominant phylum for each treatment group in both the distal colon and cecum was Firmicutes (96.5–98.7%), with other phyla including Proteobacteria (0.9–3.0%), Actinobacteria (0.1–0.5%), and Bacteroidetes (~0.1%).
Relative abundance at the phylum level for each of the treatment groups. Specific phyla are shown in the key on the right.
Bacterial taxa with significant differences in abundance in pairwise comparisons between treatment groups as determined by Mann–Whitney U-tests.
Taxonomic rank
Treatment Groupsa
Taxa with significant differences
Mann–Whitney p value
Control
BMD
Co
Ce
Co
Ce
Phylum
Proteobacteria
0.043
Class
Bacilli
0.019
Class
Clostridia
0.007
Genus
Lachnospiraceae
0.008
Genus
Oscillospira
0.008
Genus
Erysipelotrichaceae cc_115
0.034
Genus
Enterococcus
0.009
Genus
Peptostreptococcaceae
0.035
Genus
Lachnospiraceae
0.049
Co, colon; Ce, cecum. Filled cells show pairwise comparisons associated with significant differences in taxonomic abundance.
At the class level (Table 2), Bacilli were depleted in the cecum of the control fed group compared to the distal colon site (p = 0.01865) and Clostridia were enriched in the cecum of the control birds compared to the distal colon (p = 0.007). Figure 4 shows the relative abundance of the classes from each treatment group. The Firmicutes Clostridia and Bacilli were the dominant class (77.7–94.5 and 2.0–15.6%, respectively), with the remaining classes including Erysipelotrichi (1.5–3.2%), Gammaproteobacteria (0.9–2.9%), and Actinobacteria (0.1–0.5%). The Clostridia appeared to comprise a greater relative abundance in the cecum (84.3–94.5%) of both treatment groups compared to the distal colon (60.2–87.2%). There were no significant differences between the cecal samples of the control and BMD treated groups or the cecal samples and the distal colon samples of the BMD group.
Relative abundance at the class level for each of the treatment groups. Specific classes are shown in the key on the right.
Selected genera were also significantly altered by the BMD treatment (Table 2). Specifically, two genera, Oscillospira and an unnamed member within the Erysipelotrichaceae family (cc_115), were depleted in the distal colon of the BMD supplemented group compared to the distal colon of the control group (p = 0.008 and 0.034, respectively). Conversely, an unknown genus in the family Lachnospiraceae was enriched in the distal colon of the BMD treated group compared to the distal colon of the control group (p = 0.008). The same microorganism was also depleted in the cecum compared to the colon of the BMD treated birds (p = 0.049). Only an unknown genus in the family Peptostreptococcaceae was depleted in the cecum of the BMD treated group compared to the cecum of the control group (p = 0.034). Figure 5 summarizes the genera that comprised each of the treatment groups and reveals that no one genus dominated in abundance.
Relative abundance at the genus level for each of the treatment groups. Specific genera are shown in the key on the right.
Discussion
Consistent with previous studies, the distal GI tracts of the birds surveyed here were dominated by Firmicutes (Figure 3), which included the classes Clostridia and Bacilli. Members of the phyla Proteobacteria, Actinobacteria, Bacteroidetes, and Tenericutes were also detected (47–49). The differences in the distribution of the microbiota were more evident at the class and genus level with the greatest changes observed in the distal colon in the BMD treated birds. Faith's phylogenetic diversity revealed an increase in diversity in the colon of the BMD treated birds. One explanation for this is that the antibiotic treatment depleted some gram-positive species that allowed others to expand in their place without significantly altering the total number of OTUs observed within the samples. In general, the composition of cecal samples showed more similarity than samples recovered from the distal colon for both the control and BMD treated birds. There were also fewer changes between the cecal samples of the control and BMD groups compared to the distal colon samples. The observed differences in the microbial populations and relative stability of the cecum is consistent with previous studies comparing cecal populations and feces in broiler chickens (36), and indicates that the cecal microbiota is buffered to some extent from the antimicrobial effects of BMD. This likely holds true for other feed additives in chickens (31, 50).
Individual variation was more also evident in the distal colon samples than the cecal samples. Bird to bird variation is not uncommon in chickens and may be explained by the immediate environment having significant impact on how domestic fowl acquire their microbiota over maternal sources as is observed in most mammals. Chicks hatched in commercial settings are typically not exposed to the hen's microbiota post hatching, therefore colonization depends on environmental factors and could be affected by the surroundings, litter management practices, and contact with other chicks (50, 51).
As cited in the introduction, there have been several studies showing that sub-therapeutic doses of antibiotics can alter the microbiota. To enhance these effects, we used a BMD dose considered to be therapeutic designed to reduce potential pathogens during outbreaks of GI diseases such as necrotic enteritis (8, 37). At this dosage, the Phylum Proteobacteria was reduced (p = 0.04) in the distal colon of the BMD treated birds when compared to the control. BMD treatment also decreased Oscillospira, Peptostreptococcaceae, and an unknown Erysipelotrichaceae in the colon of the birds, while only Peptostreptococcaceae was depleted in the cecum of the BMD group.
Between treatment groups, the control samples clustered closer together in the Unifrac PCoA plots compared to the BMD groups. The individual variation among birds, as well as variability in feed consummation may have contributed to the lack of tight clustering in the BMD groups. “Pecking order” among the birds may also influence feed (and BMD) consumption within the groups (52).
In general, few OTUs were significantly different between the control and BMD treated group in the distal colon. While the more proximal GI has greater susceptibility to antibiotics than the distal GI, this may also indicate the bacteria in the distal colon of the chicks have a higher proportion of bacteria that are resistant to bacitracin (22, 53, 54).
These results have implications for development of new strategies as use of antibiotics for growth promotion is being phased out of commercial use. Specifically, selective manipulation of the microbiome through alternative supplementation is growing in interest as an alternative to antibiotics (55, 56). This can include the use of beneficial bacteria as probiotics, prebiotic supplementation, phytobiotics, or enzymes (30). Antibiotic alternatives can confer resistance to colonization of pathogens through competitive inhibition, decreasing pH of the GI tract, or by contributing to overall animal health through immune modulation (55, 57, 58).
Probiotics can consist of one or more Gram-positive bacteria, such as Lactobacillus, Enterococcus, or Bacillus, as well as multiple strains of the same species. For example, Nisbet found administration of 29 microorganisms from ceca of older, Salmonella free chicks protected 1-day old chicks from various enteric pathogens (59, 60). Probiotics also have potential to replace growth-promoting antibiotics as evidenced by studies showing broilers fed Bacillus subtilis daily had increased weight gain and improved feed conversion ratios than the control animals (58). Additional benefits of probiotics include enhanced production of short chained fatty acids (butyrate, acetate, and proprionate) in broilers fed non-digestible carbohydrates or oligosaccharides that are fermented by members of the microbiota, as well as enhanced protection and antimicrobial production (30). Prebiotics can increase adaptive immune responses when administered in ovo and aided in intestinal development in newly hatched chicks (61).
The relative stability of the chicken microbial community in the cecum compared to the distal colon may prove to be beneficial to production practices that seek to exploit the microbiome to enhance production. For example, modern poultry production typically prevents contact between chicks and older birds. This means the chicks are exposed to environmental bacteria rather than those associated with a healthy bird. Cecal transplants, or beneficial bacteria sourced from the cecum, could be used as “seeds” for post hatched chicks. While more studies are needed, the cecal microbiota presents itself as a potential source for colonizing newly hatched chicks with a healthy, chicken specific community that can speed GI development and help prevent diseases such as necrotic enteritis.
Ethics Statement
This study was carried out with the approval of the Iowa State University Institutional Animal Care and Usage Committee under protocol number 6-11-7167-G.
Author Contributions
AP designed and conducted the experiments, analyzed the data, and wrote the manuscript. GP assisted with data analysis and helped write the manuscript.
Conflict of Interest Statement
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.
The authors thank Jerry Snyder and Glenn Songer for their help with the animal experiments.
ReferencesFrostAJWoolcockJB. Antibiotics and animal production. In: Woolcock JB, editor. Microbiology of Animals and Animal Products. New York, NY: Elsevier (1991). p. 181–94.DibnerJJRichardsJD. Antibiotic growth promoters in agriculture: history and mode of action. Poult Sci. (2005) 84:634–43. 10.1093/ps/84.4.63415844822SinghMChauhanSSKumarP. Effect of supplementation of diets with BMD and virginiamycin on the growth performance, carcass characteristics and bacterial population in broiler chickens. Vet World. (2008) 1:141–3.SmithJLWeinbergED. Mechanisms of antibacterial action of bacitracin. J Gen Microbiol. (1962) 28:559–69. 10.1099/00221287-28-3-55913914300ButayePDevrieseLAHaesebrouckF. Antimicrobial growth promoters used in animal feed: effects of less well known antibiotics on gram-positive bacteria. Clin Microbiol Rev. (2003) 16:175–88. 10.1128/CMR.16.2.175-188.200312692092DroumevD. Review of antimicrobial growth promoting agents available. Vet Res Commun. (1983) 7:85–99. 10.1007/BF022286016364557MurugesanGRSyedBHaldarSPenderC. Phytogenic feed additives as an alternative to antibiotic growth promoters in broiler chickens. Front Vet Sci. (2015) 2:21. 10.3389/fvets.2015.00021PrescottJFSivendraRBarnumDA. The use of bacitracin in the prevention and treatment of experimentally-induced necrotic enteritis in the chicken. Can Vet J. (1978) 19:181–3. 698898FasinaYONewmanMMStoughJMLilesMR. Effect of Clostridium perfringens infection and antibiotic administration on microbiota in the small intestine of broiler chickens. Poult Sci. (2015) 95:247–60. 10.3382/ps/pev32926567176KaldhusdalMBenestadSLLøvlandA. Epidemiologic aspects of necrotic enteritis in broiler chickens - disease occurrence and production performance. Avian Pathol. (2016) 45:271–4. 10.1080/03079457.201626956946MarshallBMLevySB. Food animals and antimicrobials: impacts on human health. Clin Microbiol Rev. (2011) 24:718–33. 10.1128/CMR.00002-1121976606LinJHunkapillerAALaytonACChangYJRobbinsKR. Response of intestinal microbiota to antibiotic growth promoters in chickens. Foodborne Pathog Dis. (2013) 10:331–7. 10.1089/fpd.2012.134823461609MaronDFSmithTJSNachmanKE. Restrictions on antimicrobial use in food animal production: an international regulatory and economic survey. Global Health. (2013) 2013:48. 10.1186/1744-8603-9-48YeomanCJChiaNJeraldoPSiposMGoldenfeldNDWhiteBA. The microbiome of the chicken gastrointestinal tract. Anim Health Res Rev. (2012) 13:89–99. 10.1017/S146625231200013822853945OakleyBBLillehojHSKogutMHKimWKMaurerJJPedrosoA. The chicken gastrointestinal microbiome. FEMS Microbiol Lett. (2014) 360:100–12. 10.1111/1574-6968.1260825263745Borda-MolinaDSeifertJCamarinha-SilvaA. Current perspectives of the chicken gastrointestinal tract and Its microbiome. Comput Struct Biotechnol J. (2018) 16:131–9. 10.1016/j.csbj.2018.03.00230026889ClavijoVFlórezMJV. The gastrointestinal microbiome and its association with the control of pathogens in broiler chicken production: a review. Poult Sci. (2018) 97:1006–21. 10.3382/ps/pex35929253263KersJGVelkersFCFischerEAJHermesGDAStegemanJASmidtH. Host and environmental factors affecting the intestinal microbiota in chickens. Front Microbiol. (2018) 9:235. 10.3389/fmicb.2018.0023529503637StutzMWJohnsonSLJudithFR. Effects of diet, bacitracin, and body weight restrictions on the intestine of broiler chicks. Poult Sci. (1983) 62:1626–32. 6634597EngbergRMHedemannMSLeserTDJensenBB. Effect of zinc bacitracin and salinomycin on intestinal microflora and performance of broilers. Poult Sci. (2000) 79:1311–9. 10.1093/ps/79.9.131111020077KnarreborgASimonMAEngbergRMJensenBBTannockGW. Effect so fdietary fat source and subtherapeutic levels of antibiotic on the bacterial community in the ileum of broiler chickens at variou sages. Appl Environ Microbiol. (2002) 68:5918–24. 10.1128/AEM.68.12.5918-5924.2002DumonceauxTJHillJEHemmingsenSMVan KesselAG. Characterization of intestinal microbiota and response to dietary virginiamycin supplementation in the broiler chicken. Appl Environ Microbiol. (2006) 72:2815–23. 10.1128/AEM.72.4.281516597987PedrosoAAMentenJFLambaisMRRacanicciAMLongoFASorbaraJO. Intestinal bacterial community and growth performance of chickens fed diets containing antibiotics. Poult Sci. (2006) 85:747–52. 10.1093/ps/85.4.74716615359GongJYuHLiuTGillJJChambersJRWheatcroftR. 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–82. 10.1111/j.1365-2672.2007.03699.x18201175LuJHofacreCSmithFLeeMD. Effects of feed additives on the development on the ileal bacterial community of the broiler chicken. Animal. (2008) 2:669–76. 10.1017/S175173110800189422443592DanzeisenJLKimHBIsaacsonRETuZJJohnsonTJ. Modulations of the chicken cecal microbiome and metagenome in response to anticoccidial and growth promoter treatment. PLoS ONE. (2011) 6:e27949. 10.1371/journal.pone.002794922114729TorokVAHughesRJMikkelsenLLPerez-MaldonadoRBaldingKMacAlpineR. Identification and characterization of potential performance-related gut microbiotas in broiler chickens across various feeding trials. Appl Environ Microbiol. (2011) 77:5868–78. 10.1128/AEM.00165-1121742925LaVorgnaMSchaefferJLBadeDDicksonJCooksonKDavisSW. Performance of broilers fed a broader spectrum antibiotic (virginiamycin) or a narrower spectrum antibiotic (bacitracin methylene disalicylate) over 3 consecutive grow-out cycles. J Appl Poult Res. (2013) 22:574–82. 10.3382/japr.2012-00703NeumannAPSuenG. Differences in major bacterial populations in the intestines of mature broilers after feeding virginiamycin or bacitracin methylene disalicylate. J Appl Microbiol. (2015) 119:1515–26. 10.1111/jam.12960PourabedinMGuanLZhaoX. Xylo- oligosaccharides and virginiamycin differentially modulate gut microbial composition in chickens. Microbiome. (2015) 3:15. 10.1186/s40168-015-0079-425874109CostaMCBessegattoJAAlfieriAAWeeseJSFilhoJABObaA. Different antibiotic growth promoters induce specific changes in the cecal microbiota membership of broiler chicken. PLoS ONE. (2017) 12:e0171642. 10.1371/journal.pone.017164228222110Díaz CarrascoJMRedondoEAPin VisoNDRedondoLMFarberMDFernández MiyakawaME. Tannins and bacitracin differentially modulate gut microbiota of broiler chickens. Biomed Res Int. (2018) 2018:1879168. 10.1155/2018/187916829682522GubanJKorverDRAllisonGETannockGW. Relationship of dietary antimicrobial drug administration with broiler performance, decreased population levels of Lactobacillus salivarius, and reduced bile salt deconjugation in the ileum of broiler chickens. Poult Sci. (2006) 85:2186–94. 10.1093/ps/85.12.218617135676LinJ. Antibiotic growth promoters enhance animal production by targeting intestinal bile salt hydrolase and its producers. Front Microbiol. (2014) 5:33. 10.3389/fmicb.2014.0003324575079ChoiJHKimGBChaCJ. Spatial heterogeneity and stability of bacterial community in the gastrointestinal tracts of broiler chickens. Poult Sci. (2014) 93:1942–50. 10.3382/ps.2014-0397424931967OakleyBBKogutMH. Spatial and temporal changes in the broiler chicken cecal and fecal microbiomes and correlations of bacterial taxa with cytokine gene expression. Front Vet Sci. (2016) 3:11. 10.3389/fvets.2016.0001126925404HofacreCLFroymanRGeorgeBGoodwinMABrownJ. Use of Aviguard, virginiamycin, or bacitracin MD against Clostridium perfringens-associated necrotizing enteritis. J Appl Poultry Res. (1998) 7:412–8. 10.1093/japr/7.4.412CaporasoJGKuczynskiJStombaughJBittingerKBushmanFDCostelloEK. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods. (2010) 7:335–6. 10.1038/nmeth.f.30320383131ReederJKnightR. Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Methods. (2010) 7:668–9. 10.1038/nmeth0910-668b20805793EdgarRC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. (2010) 26:2460–1. 10.1093/bioinformatics/btq46120709691EdgarRCHaasBJClementeJCQuinceCKnightR. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. (2011) 27:2194–200. 10.1093/bioinformatics/btr38121700674CaporasoJGBittingerKBushmanFDDesantisTZAndersenGLKnightR. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. (2010) 26:266–7. 10.1093/bioinformatics/btp63619914921DeSantisTZHugenholtzPLarsenNRojasMBrodieELKellerK. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. (2006) 72:5069–72. 10.1128/AEM.03006-0516820507McDonaldDPriceMNGoodrichJNawrockiEPDesantisTZProbstA. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. (2011) 6:610–8. 10.1038/ismej.2011.13922134646PriceMNDehalPSArkinAP. FastTree 2 - approximately maximum-likelihood trees for large alignments. PLoS ONE. (2010) 5:e9490. 10.1371/journal.pone.000949020224823Vázquez-BaezaYPirrungMGonzalezAKnightR. EMPeror : a tool for visualizing high - throughput microbial community data. Giga Sci. (2013) 2:1–4. 10.1186/2047-217X-2-1624280061SekeljaMRudIKnutsenSHDenstadliVWesterengBNæsT. Abrupt temporal fluctuations in the chicken fecal microbiota are explained by its gastrointestinal origin. Appl Environ Microbiol. (2012) 78:2941–8. 10.1128/AEM.05391-1122307311WeiSMorrisonMYuZ. Bacterial census of poultry intestinal microbiome. Poult Sci. (2013) 92:671–83. 10.3382/ps.2012-0282223436518MancabelliLFerrarioCMilaniCMangifestaMTurroniFDurantiS. Insights into the biodiversity of the gut microbiota of broiler chickens. Environ Microbiol. (2016) 18:4727–38. 10.1111/1462-2920.1336327129897TorokVAllisonGEPercyNJOphel-KellerKHughesRJ. Influence of antimicrobial feed additives on broiler commensal posthatch gut microbiota development and performance. Appl Environ Microbiol. (2011) 77:3380–90. 10.1128/AEM.02300-1021441326CressmanMDYuZNelsonMCMoellerSJLilburnMSZerbyHN. Interrelations between the microbiotas in the litter and in the intestines of commercial broiler chickens. Appl. Environ. Microbiol. (2010) 76:6572–82. 10.1128/AEM.00180-1020693454PagelMDawkinsMS. Peck orders and group size in laying hens: ‘futures contracts’ for non-aggression. Behav Process. (1997) 40:13–25. 24897609AndreaFGoldenbergerDTeuberM. Antimicrobial susceptibility of intestinal bacteria from Swiss poultry flocks before the ban of antimicrobial growth promoters. Syst Appl Microbiol. (2001) 24:116–21. 10.1078/0723-2020-00004WiseMGSiragusaGR. Quantitative analysis of the intestinal bacterial community in one- to three-week-old commercially reared broiler chickens fed conventional or antibiotic-free vegetable-based diets. J Appl Microbiol. (2007) 102:1138–49. 10.1111/j.1365-2672.2006.03153.x17381758Montalban-ArquesADe SchryverPBossierPGorkiewiczGMuleroVGatlinDMR. Selective manipulation of the gut microbiota Improves immune status in vertebrates. Front Immunol. (2015) 6:512. 10.3389/fimmu.2015.0051226500650FooJLLingHLeeYSChangMW. Microbiome engineering: current applications and Its future. Biotechnol J. (2017) 12:1–11. 10.1002/biot.20160009928133942WalshMCGardinerGEHartOMLawlorPGDalyMLynchB. Predominance of a bacteriocin-producing Lactobacillus salivarius component of a five-strain probiotic in the porcine ileum and effects on host immune phenotype. FEMS Microbiol Lett. (2008) 64:317–27. 10.1111/j.1574-6941.2008.00454.x18373687SenSIngaleSLKimYWKimJSKimKHLohakareJD. Effect of supplementation of Bacillus subtilis LS 1-2 to broiler diets on growth performance, nutrient retention, caecal microbiology and small intestinal morphology. Res Vet Sci. (2012) 93:264–8. 10.1016/J.RVSC.2011.05.02121757212NisbetDJTellezGILowryVKAndersonRCGarciaGNavaG. Effect of a commercial competitive exclusion culture (Preempt) on mortality and horizontal transmission of Salmonella gallinarum in broiler chickens. Avian Dis. (1998) 42:651–6. 9876831NisbetD. Refined competitive exclusion cultures in the prevention of enteropathogen colonisation in poultry and swine. Antonie Van Leeuwenhoek. (2002) 81:481–6. 10.1023/A:1020541603877SlawinskaAPlowiecASiwekMJaroszewskiMBednarczykM. Long-term transcriptomic effects of prebiotics and synbiotics delivered in ovo in broiler chickens. PLoS ONE. (2016) 11:e0168899. 10.1371/journal.pone.016889928002487