Bacitracin Methylene Disalicylate Improves Intestinal Health by Modulating Its Development and Microbiota in Weaned Rabbits

Intestinal infections are a major cause of morbidity and mortality in humans and agricultural animals, especially newborns and weaned animals. Preventive treatments that help weaned animals maintain homeostasis and balance the hindgut microbial populations are desirable. The present study aimed to explore the impact of bacitracin methylene disalicylate (BMD) on the intestinal health by analyzing the intestinal environment, morphology, expression of peptidoglycan recognition proteins (PGRPs), and flora of weaned rabbits. A total of 300 New Zealand weaned rabbits were randomly divided into the following five treatment groups for a 35-day feed trial: control group (basal diet), bacitracin zinc (BZ) group (50 mg/kg BZ), BMDa group (100 mg/kg BMD), BMDb group (50 mg/kg BMD), and BMDc group (rabbits fed a basal diet supplemented with 25 mg/kg BMD). In each treatment group, 28 rabbits were slaughtered for experimental analysis. The results showed that the supplementation of BMD increased the environmental acidity of the cecum of the weaned rabbits and reduced the ammonia-nitrogen concentration, which was beneficial to the survival of useful bacteria in the intestine. The morphology analysis of the duodenum using hematoxylin and eosin staining revealed that the villus length, villus/crypt ratio, and intestinal wall thickness increased in the BMD group, thereby improving the structure of the duodenum and the absorption capacity of the small intestine. Moreover, real-time polymerase chain reaction test showed that PGRPs (especially PGLYRP-1 and PGLYRP-2) in the intestinal had an antagonistic effect with BMD in the process of inhibiting pathogenic bacteria, resulting in their decreased expression (P < 0.05). Furthermore, through 16S rRNA sequencing in the cecal content, the abundance of the predominant phyla in the BMDa and BZ groups was found to be the closest. The abundance of the genera Lachnospira, Erysipelotrichaceae (p-75-a5), Paraprevotellaceae (YRC22), Mogibacterium, Peptococcaceae (rc4-4), Anaerovibrio, Succinivibrio, and Sphaerochaeta increased in the BMDa and BZ groups (P < 0.05). The relative abundance of Alistipes, Sedimentibacter, and Dorea significantly increased only in the BMDa group (P < 0.05). Conclusively, BMD, as well as microbes, improved the intestinal environment and structure to maintain the intestinal health of weaned rabbits.

Introduction

Intestinal infections are a major cause of morbidity and mortality worldwide in humans and agricultural animals, especially newborns and weaned animals. The neonatal phase in humans and animals is one of the critical periods in life, during which nutrition and management have a long-term effect on the overall performance (Videhult and West, 2016). In this stage, they grow fast and have a strong metabolism; however, the digestive system is incompletely developed, and the intestinal flora balance is easily affected. The intestine contains a microbial ecosystem, and a balanced microbiome is important in the development of a healthy intestinal mucosal immune system (Laukoetter et al., 2008). Under normal circumstances, the intestinal flora and intestinal immunity are in a dynamic balance through precise control mechanisms. Once this mechanism is destroyed, a series of diseases ensue. Therefore, preventive treatments that help weaned rabbits maintain homeostasis and balance the gut microbial populations are desirable.

Bacitracin methylene disalicylate (BMD) is the reaction product of Bacillus licheniformis secondary metabolites, bacitracin and methylene salicylic acid. BMD has a strong inhibitory effect on Gram-positive bacteria. On the contrary, BMD has no inhibitory effect on intestinal probiotics, such as Lactobacilli, Bifidobacteria, and Bacillus. It has a promoting effect on the growth and intestinal health of livestock, and has been used in the breeding of broilers, pigs, and cattle (Damron and Wilson, 1985; Song et al., 2011; Xie et al., 2013). The role of BMD in intestinal flora balance in rabbits needs further investigation. Moreover, BMD can restore intestinal flora balance and control gastrointestinal infections in rabbits.

Besides the intestinal mucosal epithelial barrier and the intestinal flora, the intestinal health also depends on the immune barrier formed by intestinal immune cells and their secretions. Peptidoglycan recognition protein (PGRP), as a bactericidal protein in innate immunity, can maintain normal intestinal flora and is associated with inflammatory bowel disease (IBD). PGRP-Iα and Nod2 can cooperate to maintain the intestinal flora balance and resist colitis (Jing et al., 2014). Several non-catalytic PGRPs can act as selective peptidoglycan receptors in the upstream of Toll and IMD pathways in Drosophila (Zaidman-Remy et al., 2011). A previous study found that diarrhea in weaned rabbits led to the upregulation of PGRPs in the cecum and duodenum (Chen et al., 2017). Therefore, the effect of the interaction between BMD and PGRP on the homeostasis of intestinal microorganisms and intestinal immunity needs further analysis.

Further, a previous study found that the addition of BMD to the diet of weaned rabbits improved their growth performance and significantly reduced diarrhea (16.67 vs. 28.33%) and mortality (5 vs. 16.67%) rates (Hao et al., 2016). Based on these responses, this study aimed to further compare the effects of BMD in the diet on intestinal health by analyzing the intestinal environment, morphology, expression of intestinal PGRPs, and intestinal flora of weaned rabbits. The findings revealed the feasibility of using BMD to restore intestinal flora balance and control gastrointestinal infections.

Materials and Methods

Ethics Statement

This study was strictly performed following the regulations for experimental animals of the China Department of Agriculture and approved by the Institutional Animal Care and Use Committee of Yangzhou University (Jiangsu, China).

Experimental Design and Animal Management

A total of 300 New Zealand weaned rabbits (35 days of age, 633.9 ± 32.62 g) were selected as experimental animals. The animals were randomly divided into the following five treatment groups (n = 60, including 30 male and 30 female) for a 35-day feed trial: (1) control group (rabbits fed a basal diet), (2) bacitracin zinc (BZ) group (rabbits fed a basal diet supplemented with 50 mg/kg BZ), (3) BMDa group (rabbits fed a basal diet supplemented with 100 mg/kg BMD), (4) BMDb group (rabbits fed a basal diet supplemented with 50 mg/kg BMD), and (5) BMDc group (rabbits fed a basal diet supplemented with 25 mg/kg BMD). The basal diet was formulated according to the nutrient requirements of weaned rabbits. The ingredient composition and nutrient levels are shown in Supplementary Table 1. This study used BZ, a polypeptide antibiotic, as the positive control. The feeding and management methods of rabbits in each group were the same. The feed was processed into pellets for free feeding and drinking. No vaccines or antibiotics were administered to these rabbits throughout the experiment.

Sampling Procedure

The rabbits were euthanized, and the pH value of the cecal contents was immediately measured. Afterward, the cecal contents were taken in a 2-mL cryopreservation tube and stored in liquid nitrogen for detecting of volatile fatty acids (VFAs) and ammonia-nitrogen (NH₃-N). Subsequently, the parts of the digestive tract were carefully separated to identify different segments of the intestines, and a 2- to 3-cm section of the intestine was cut from the same part of the duodenum. After cleaning the intestinal canal of the duodenum with 0.9% physiological saline, the canal was placed in a preformulated paraformaldehyde solution for fixation and used for preparing intestinal tissue sections. A 1- to 2-cm section of the intestinal canal was taken from the duodenum and cecum of the rabbits into a 2-mL cryopreservation tube and stored in liquid nitrogen for total RNA extraction. The samples were collected from 28 rabbits in each group, and 8 were used for microbial analysis.

Determination of Intestinal Environmental Indicators

The contents of the cecum were immediately taken out, and the pH value was quickly determined using a DELTA 320 pH meter (Mettler-Toledo, Shanghai, China). Furthermore, 2 mL of distilled water was added to 2 mL of cecal contents and centrifuged at 10,000 rpm for 10 min. Then, 1 mL of supernatant was taken, 0.2 mL of 20% metaphosphoric acid solution containing 60 mM crotonic acid was added, and 0.4 μL of supernatant was taken for injection analysis. The VFA was measured using a GC-9A Gas Chromatograph (Shimadzu, Kyoto, Japan), with nitrogen as the carrier gas and a flow rate of 30 mL/min. The standard solutions were 3.353 mg/mL acetic acid, 1.189 mg/mL propionic acid, and 0.793 mg/mL butyric acid. In addition, after centrifuging 4 mL of cecal content diluent for 10 min, 50 μL of the supernatant was taken in a 10-mL test tube. Subsequently, 3 mL each of phenol and sodium hypochlorite was added to the test tube. The test tube was cooled in a 60°C water bath for 10 min and then immediately with cold water. Then, the optical density (OD) at 546 nm was determined using a UV756 spectrophotometer (UNICO, CA, United States) to calculate the NH₃-N concentration in the cecum. The ammonia standard stock solution contained 32 mg/100 mL NH₃-N.

Determination of Intestinal Villus Height and Crypt Depth

The duodenal intestinal contents were cleaned with ice-cold physiological saline, and the intestinal canal of about 2–3 cm at the same location was collected and fixed in 4% paraformaldehyde solution. Subsequently, the fixed tissue was taken out, rinsed with running water, and then put into gradient alcohol for dehydration. The samples were cleared with xylene for about 15–20 min, dipped in wax, and embedded at 60°C. Next, after sectioning, patching, and dewaxing, hematoxylin–eosin was used for staining. The slices were observed using an optical microscope, photographed, and stored for later analysis. Each section of the intestinal tissue was cut continuously into five pieces, and six typical fields of view (clear pictures and complete villi) were selected from each tissue section. For each field of view, the longest villi and the corresponding crypt depth were selected for measurement. The average values of the measurements were recorded, and the villus height and crypt depth were calculated.

Analysis of Intestinal Gene Expression of PGRPs

The total RNA of the intestinal tissue was extracted with TRIzol. The reverse transcription (RT) reaction system was prepared using a SuperRT cDNA Kit following the manufacturer’s protocols. The system included 2 μL of the RNA template, 4 μL of dNTP Mix (2.5 mmol/L), 2 μL of Prime Mix, 4 μL of 5 × RT buffer, 1 μL of SuperRT 200 U/μL, and RNase-free water to a volume of 20 μL. Real-time quantitative RT–polymerase chain reaction (qRT-PCR) was carried out using an Applied Biosystems 7500 Real-Time RT-PCR System (Applied Biosystems, CA, United States). The following program was used: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s. According to the mRNA sequences of rabbit PGRPs published in GenBank, qRT-PCR primers of peptidoglycan recognition protein (PGLYRP)-1, PGLYRP-2, and PGLYRP-3 were designed, and are shown in Supplementary Table 2. Each sample was measured in triplicate. The 2^–ΔΔCt method was used to quantitate and analyze the expression of target genes, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the reference gene.

Microbial Analysis

The total DNA of the cecal contents was extracted using hexadecyl trimethyl ammonium bromide. Genomic DNA was amplified using the specific primers 515F 5′-GTGCCAGCMGCCGCGGTAA-3′ and 806R 5′-GGACTACHV GGGTWTCTAAT-3′ with barcode-specific primers in the V4 region of the bacterial 16S rRNA gene. The PCR products were evaluated using a Qubit R 2.0 Fluorometer (Life Technologies, CA, United States) and an Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, United States). The amplicon sequencing library was subsequently constructed, and Illumina HiSeq was used for 16S rRNA gene sequencing.

According to the barcode and primer sequence, data of each sample were obtained from the off-machine sequence. After truncating the barcode and primer sequences, the reads of each sample were spliced using FLASH (version 1.2.7¹) (Magoc and Salzberg, 2011). The Quantitative Insights Into Microbial Ecology (QIIME, version 1.7.0²) process was used for quality control and data processing (Caporaso et al., 2010). The effective tags of the samples were clustered using the UPARSE software (version 7.0.1001³), and the sequences were clustered into operational taxonomic units (OTUs) with >97% identity (Edgar, 2013). Meanwhile, the sequences with the highest frequency were selected as the representative sequences of OTUs according to their algorithm principles. Species annotations were analyzed using Ribosomal Database Project (RDP, version 2.2⁴) Classifier and Greengene databases⁵ (with a threshold of 0.8–1), and the community composition of each sample was counted at the kingdom, phylum, class, order, family, genus, and species levels (DeSantis et al., 2006; Wang et al., 2007). The data were normalized with the least amount of data as the standard. The alpha- and beta-diversity analyses were assessed on the data after normalization by in-house Perl scripts. The QIIME software was used to calculate the observed species, Chao1, Shannon, Simpson, and abundance-based coverage estimator (ACE) indices. The vegan package of R software (version 2.15.3) was used for non-metric multidimensional scaling (NMDS) analysis to reflect the differences between and within samples (Noval Rivas et al., 2013). The raw data for this project were deposited with the national center for biotechnology information (NCBI, accession number: PRJNA530911).

Statistical Analysis

For microbial analysis, the multiresponse permutation procedure (MRPP) was used to analyze whether significant differences existed in microbial community structure. Metastat analysis was carried out using R software (version 2.15.3) to find the species with significant differences (White et al., 2009). At each classification level, permutation tests were performed to obtain the P-values, and the Benjamini and Hochberg false discovery rate method was used to further correct the P-values. Other indices except microbial analysis were used to evaluate the significance using the one-way analysis of variance of SPSS 21. The results were expressed as mean ± standard deviation in tables and figures at the level of significance (P < 0.05).

Results

Intestinal Morphology in Weaned Rabbits

As shown in Figure 1 and Table 1, duodenal histological analysis revealed that the villi were significantly longer in the BMDa, BMDb, and BZ groups than in the BMDc and control groups (P < 0.05). The crypt depth was significantly smaller in the BZ group than in other groups (P < 0.05), while it was not significantly different in the BMD groups and control group (P > 0.05). The villus/crypt ratio was significantly larger in the BMDa group than in the control group (P < 0.05). The duodenal wall thickness was significantly greater in the BMD and BZ groups than in the control group (P < 0.05). These findings indicated that BMD could affect the morphology to improve the absorption capacity of the intestine.

FIGURE 1

Figure 1. Effect of BMD on the morphology of duodenum in weaned rabbits.

TABLE 1

Table 1. Comparison of BMD based on the morphology of the duodenum in weaned rabbits.

Intestinal Environment in Weaned Rabbits

As shown in Table 2, the concentration of acetic acid in the cecal contents was significantly higher in the BMD and BZ groups than in the control group due to BMD supplementation (P < 0.05). BMD, BZ, and control groups had no significant difference in propionic acid concentration, butyric acid concentration, and acetic acid/(propionic acid + butyric acid) (P > 0.05). In the cecal contents, the pH value was significantly lower in the BMDb group than in the other groups (P < 0.05); no significant difference was found in the pH values of the other groups (P > 0.05). The NH₃-N content decreased gradually with the increase in BMD addition and was significantly lower in the BMD group than in the control group (P < 0.05). The NH₃-N content was significantly lower in the BMDa group than in the BMDb and BMDc groups (P < 0.05), and significantly lower in the BZ group compared with the BMDc and control groups (P < 0.05). Therefore, supplementation with 50–100 mg/kg BMD increased the environmental acidity of the cecum in the weaned rabbits and reduced the NH₃-N content to benefit the survival of useful bacteria in the intestines.

TABLE 2

Table 2. Comparison of BMD based on cecal fermentation in weaned rabbits.

Gene Expression of Intestinal PGRPs in Weaned Rabbits

The expression levels of PGRPs in different dose groups were tested to further verify the beneficial effects of BMD on the intestinal flora. The expression levels of PGLYRP-1, PGLYRP-2, and PGLYRP-3 in the duodenum and cecum were detected. The expression levels of PGRPs in the duodenum and cecum followed the same trend, that is, the BMD and BZ groups had decreased expression levels compared with the control group (Figure 2). The expression levels of PGLYRP-1 and PGLYRP-2 were significantly lower in the BMD and BZ groups than in the control group (P < 0.05). The expression level of PGLYRP-3 was not significantly different among BMD, BZ, and control groups (P > 0.05).

FIGURE 2

Figure 2. Expression of PGRPs in the cecum and duodenum. (A) Expression of PGRPs in the cecum. (B) Expression of PGRPs in the duodenum. a,b indicate significant differences within a row (P < 0.05).

Diversity, Richness, and Composition of Bacterial Communities in the Cecal Content

In terms of alpha diversity (Table 3), no differences in observed species, Shannon and Simpson indices, were noted in the BMDa, BMDb, BZ, and control groups. BMD and BZ supplementation increased the ACE and Chao1 richness compared with the control group, especially in the BMDa and BZ groups (P < 0.05). The NMDS analysis based on OTUs was conducted for comparison between the groups (Supplementary Figure 1). The stress was less than 0.2, indicating that the NMDS analysis had a certain reliability. The samples in these groups clustered due to different treatments. The findings combined with the MRPP analysis results showed that the differences in the microbial community structure among the groups were highly significant (P < 0.05), except the BMDa vs. BZ groups (Supplementary Figure 1 and Supplementary Table 3; P > 0.05).

TABLE 3

Table 3. Number of observed species, richness, and diversity indices in the cecal content of each treatment.

All the reads of the samples in the cecum were processed together. In total, bacterial sequencing generated 2,324,160 raw reads. Based on the > 97% sequence similarity, 2010 OTUs were identified, all of which belonged to the bacterial domain by the Greengenes classification. Moreover, 1,154 ± 91 OTUs per sample, on average, were identified. Within the bacterial population, the top 10 predominant phyla in the cecal content were identified across all samples at the class level (Figure 3). Firmicutes, Bacteroidetes, Verrucomicrobia, and Proteobacteria were the four dominant phyla, representing 66.29 ± 4.77%, 15.65 ± 4.89%, 4.50 ± 2.14%, and 3.74 ± 0.85% of the total sequences, respectively. Synergistetes and Tenericutes represented average percentages of 1.41 ± 1.19% and 2.77 ± 0.75%, respectively, of the total sequences. The proportion of phyla Actinobacteria, Cyanobacteria, Spirochaetes, and Fusobacteria was less than 1% of the total microbial community. Besides, the abundance of Firmicutes, Bacteroidetes, Proteobacteria, and Synergistetes presented a correlated change trend with the addition of BMD. In this aspect, the proportion of these phyla in the BMDa group most closely matched with that in the BZ group, followed by BMDb and control groups. Comparing the BMDa and BZ (positive control) groups with the control group revealed that the proportion of Bacteroidia, Gammaproteobacteriac, Deltaproteobacteria, Betaproteobacteria, and Bacillus increased while the proportion of Verrucomicrobiae and Mollicutes decreased due to the addition of high-dose BMD and BZ.

FIGURE 3

Figure 3. Percentage composition of the top 10 predominant phyla in the cecal content.

Effect of BMD Supplementation on the Relative Abundance of Bacterial Communities

The similarities and differences between the communities in different groups were compared and presented as a Venn diagram (Figure 4). The BMDa, BMDb, BMDc, BZ, and control groups had 1171 OTUs in common, with common OTUs comprising 67.84, 78.64, 79.88, 68.36, and 80.54% of the sequences in these communities, respectively.

FIGURE 4

Figure 4. Shared operational taxonomic unit (OTU) analysis of the different groups. Each circle in the figure indicates a group, and the numbers in the circle and circle overlap represent the number of OTUs between the groups. The number not in overlap indicates the number of unique OTUs in the group.

The effect of BMD supplementation on the relative abundance of bacterial genera in the cecal content using 16S rRNA sequencing is illustrated in Table 4. At the genus level, 113 taxa were detected, and 30 genera were affected (P < 0.05) by treatments. Compared with the control group, BMD and BZ supplementation increased (P < 0.05) the proportions of Turicibacter, Epulopiscium, and others. Specifically, the abundance of Lachnospira, Erysipelotrichaceae (p-75-a5), Paraprevotellaceae (YRC22), Mogibacterium, Peptococcaceae (rc4-4), Anaerovibrio, Succinivibrio, and Sphaerochaeta increased in the BMDa and BZ groups (P < 0.05). The relative abundance of Alistipes, Sedimentibacter, and Dorea significantly increased only in the BMDa group (P < 0.05).

TABLE 4

Table 4. Effect of BMD supplementation on the relative abundance of bacterial genera in the cecal content detected using 16S rRNA sequencing.

Discussion

The intestinal mucosa serves as a filter between the intestinal lumen environment and homeostasis. It not only regulates the passage of nutrients and molecules but also prevents the penetration of bacteria, toxins, and dietary antigens into the submucosal tissues or circulatory system of the body (Groschwitz and Hogan, 2009). An effective intestinal barrier contains many constituent elements, such as the mucus layer, immune effect factor, and microbiome.

Changes in the intestinal tissue morphology can reflect the health of animals and their ability to digest and absorb nutrients to a certain extent (DeLegge, 2008). Weaned rabbits were susceptible to weaning stress and other changes in conditions, resulting in changes in tissue structure, such as intestinal villous atrophy, increased crypt depth, and mucosal damage, reducing the ability to absorb nutrients. The addition of BMD significantly increased the thickness of the intestinal wall, and also had a significant effect on the villus length and the villus/crypt ratio. The villi were longer in the BMD treatment groups compared with the control group in turkeys and pigs, which was typically equated with excellent gut health and high absorptive efficiency (Sims et al., 2004; Song et al., 2011). Moreover, the villus length in the BMDa and BMDb groups increased, indicating that the number of intestinal epithelial cells and the absorptive efficiency increased. The villus/crypt ratio was larger in the BMD group than in the control group, indicating that the addition of BMD changed the morphological structure of the duodenum, which was beneficial to the absorption of nutrients and maintained the steady state of the intestine.

The cecum is vital in the digestion and absorption of nutrients in rabbit diets, especially crude fiber. An appropriate cecal internal environment helped synthesize microbial proteins and prevent rabbit intestinal inflammation (Gidenne and Bellier, 2000). The main relevant indicators of the cecal internal environment included pH, VFA, and NH₃-N. In this study, the addition of BMD increased the content of acetic acid, and the content of VFA stimulated the growth of intestinal mucosa, reduced alkaline damage, and enabled beneficial microorganisms in the intestine to grow and reproduce in a suitable environment. The content of NH₃-N was significantly lower in the BMD group than in the control group. The content of NH₃-N showed a significantly decreasing trend with the increase in the amount of BMD added, suggesting that the ability of microorganisms to use NH₃-N was enhanced to facilitate the synthesis of more bacterial proteins. In summary, the addition of BMD caused pH, VFA, and NH₃-N to interact with each other, resulting in a decrease in the alkalinity of the intestinal environment, which helped maintain a healthy state of the rabbit intestine.

It has been reported that BMD can regulate the innate immune responses in broilers (Brennan et al., 2013). However, the effect of BMD on intestinal immunity of rabbits is unclear. PGRPs were known innate immune-activating molecules highly conserved in insects and mammals (Liu et al., 2001; Dziarski, 2004). Up to now, 19 kinds of PGRPs have been found in insects. They activated Toll and IMD signaling pathways, induced protein lyase expression and phagocytosis, and participated in the hydrolysis of peptidoglycan and immunity against infection (Werner et al., 2000). Mammals contain four kinds of PGRPs, named PGLYRP-1, PGLYRP-2, PGLYRP-3, and PGLYRP-4 (formerly known as PGRP-S, PGRP-L, PGRP-Iα, and PGRP-Iβ, respectively), all possessing very high antibacterial activity (Kang et al., 1998; Liu et al., 2001). In previous studies, the expression levels of PGLYRP-1, PGLYRP-2, and PGLYRP-3 in the duodenum and cecum of rabbits with diarrhea increased, which were related to the increase in pathogenic bacteria in the intestines (Chen et al., 2017). After adding BMD, the expression levels of PGLYRP-1, PGLYRP-2, and PGLYRP-3 reduced compared with those of the control group. Both BMD and PGRPs have antibacterial effects, and their antagonism may cause the downregulation of PGRPs. However, the regulation mechanism in animals is very complicated. A more in-depth study is needed to prove the relationship between BMD, intestinal harmful flora, and PGRPs.

Since colonization, intestinal flora has become an important part of human and animal life. During the growth of the organisms, the flora and host interact to form a complex ecosystem that constantly changes and is relatively stable for a long time. Studies have shown that the imbalances caused by the changes in the composition, ratio, and diversity of early intestinal flora in infants and young children might be a potential factor for the occurrence of diseases, such as allergies, obesity, IBD, and intestinal colic, and largely affect the growth and development of diseases in the future (Videhult et al., 2016). Among monogastric livestock, the relative volume of rabbit cecum is larger. Microbes digested food residues in the cecum, and the cecum provided suitable conditions for the activities of microorganisms. In this study, the abundance of Clostridia, Bacteroidia, Synergistia, Deltaproteobacteria, and Betaproteobacteria displayed a correlated change trend with the addition of BMD, suggesting that BMD affected the diversity, richness, and composition of bacterial communities in the intestinal contents of weaned rabbits. Further, the analysis of the relative abundance of bacterial communities revealed that the proportion of Turicibacter, Epulopiscium, Lachnospira, Erysipelotrichaceae (p-75-a5), Paraprevotellaceae (YRC22), Mogibacterium, Peptococcaceae (rc4-4), Anaerovibrio, Succinivibrio, and Sphaerochaeta increased due to the addition of BMD and BZ. Turicibacter participated in fermentation metabolism, and lactic acid was its main metabolite. The mixture of Turicibacter with Clostridia restored intestinal serotonin in mice to normal levels (Yano et al., 2015). The abundance of Erysipelotrichaceae positively correlated with carbohydrate digestion (Cox et al., 2013; Young et al., 2013; Bermingham et al., 2017). The genus Sphaerochaeta lacked the genes encoding the characteristic flagellar apparatus and, contrary to most other spirochetes, acquired many metabolic and fermentation genes from clostridia (Caro-Quintero et al., 2012). This helped rabbits metabolize carbohydrates after eating high-fiber and low-protein forage. Among the affected genera, the proportion of Alistipes, Sedimentibacter, and Dorea increased only under BMDa treatment. Alistipes, as an intestinal bacterium, could produce short-chain fatty acids, which helped in the fermentation of food without producing alcohol. Generally, it improves the metabolic state of the intestine, which is more conducive to the use of glucose for energy (Yu et al., 2019). Dorea is the main gas-producing bacteria in the human intestine, which uses carbohydrates to produce gas. The proportion of Dorea in IBD significantly reduces, suggesting Lachnospiraceae (Dorea) is important in maintaining intestinal health (Sokol et al., 2017). These results suggested that BMD significantly changed the bacterial microbiome of the cecal content to affect intestinal function.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://doi.org/10.6084/m9.figshare.12609920.

Ethics Statement

The animal study was reviewed and approved by the Institutional Animal Care and Use Committee of Yangzhou University.

Author Contributions

YC was responsible for the collection, analysis of results and wrote the manuscript. SH, JL, and BZ performed experiments. NY, TZ, SL, and SB prepared figures and tables. XW designed the study. All authors read and approved the final manuscript.

Funding

This study was supported by the National Key Research and Development Program of China (Grant No. 2018YFD0502203), the China Agriculture Research System of MOF and MARA (Grant No. CARS-43-A-1), and the National Natural Science Foundation of China (Grant No. 31702081).

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/fmicb.2021.579006/full#supplementary-material

Supplementary Figure 1 | Non-metric multidimensional scaling analysis based on OTUs. Each point in the graph indicates a sample, and the distance between the points indicates the degree of difference. The same group of samples is shown in the same color.

Supplementary Table 1 | Ingredient composition and nutrient levels of the basal diets.

Supplementary Table 2 | Primers used for real-time qRT-PCR.

Supplementary Table 3 | Differential analysis among groups using MRPP. The smaller the observed-delta value, the smaller the difference within the group; the larger the expected-delta value, and the greater the difference among groups. An A value > 0 means that the difference among groups is greater than that within the group; otherwise, it means that the difference within the group is greater than that among groups. A P < 0.05 indicates a significant difference.

Abbreviations

ACE, Abundance-based coverage estimator; BMD, Bacitracin methylene disalicylate; BZ, Bacitracin zinc; CTAB, Hexadecyl trimethyl ammonium bromide; GAPDH, Glyceraldehyde-3-phosphate dehydrogenase; IBD, Inflammatory bowel disease; OD, Optical density; OTUs, Operational taxonomic units; MRPP, Multiresponse permutation procedures; NMDS, Non-metric multidimensional scaling; NCBI, National center for biotechnology information; PGRP, Peptidoglycan recognition protein; QIIME, Quantitative Insights Into Microbial Ecology; RDP, Ribosomal database project; VFA, Volatile fatty acid.

Footnotes

1. ^ http://ccb.jhu.edu/software/FLASH/
2. ^ http://qiime.org/scripts/split_libraries_fastq.html
3. ^ http://drive5.com/uparse/
4. ^ http://sourceforge.net/projects/rdp-classifier/
5. ^ http://greengenes.lbl.gov/cgi-bin/nph-index.cgi

References

Bermingham, E. N., Maclean, P., Thomas, D. G., Cave, N. J., and Young, W. (2017). Key bacterial families (Clostridiaceae, Erysipelotrichaceae and Bacteroidaceae) are related to the digestion of protein and energy in dogs. PeerJ 5:e3019. doi: 10.7717/peerj.3019

PubMed Abstract | CrossRef Full Text | Google Scholar

Brennan, K. M., Graugnard, D. E., Xiao, R., Spry, M. L., Pierce, J. L., Lumpkins, B., et al. (2013). Comparison of gene expression profiles of the jejunum of broilers supplemented with a yeast cell wall-derived mannan oligosaccharide versus bacitractin methylene disalicylate. Br. Poult. Sci. 54, 238–246. doi: 10.1080/00071668.2013.775404

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336.

Google Scholar

Caro-Quintero, A., Ritalahti, K. M., Cusick, K. D., Loffler, F. E., and Konstantinidis, K. T. (2012). The chimeric genome of Sphaerochaeta: nonspiral spirochetes that break with the prevalent dogma in spirochete biology. mBio 3, e00025–12.

Google Scholar

Chen, Y., Zhao, B., Wu, Y., Hu, S., Mu, L., Zhu, C., et al. (2017). Impacts of diarrhea on the immune system, intestinal environment, and expression of PGRPs in New Zealand rabbits. PeerJ 5:e4100. doi: 10.7717/peerj.4100

PubMed Abstract | CrossRef Full Text | Google Scholar

Cox, L. M., Cho, I., Young, S. A., Anderson, W. H. K., Waters, B. J., Hung, S., et al. (2013). The nonfermentable dietary fiber hydroxypropyl methylcellulose modulates intestinal microbiota. FASEB J. 27, 692–702. doi: 10.1096/fj.12-219477

PubMed Abstract | CrossRef Full Text | Google Scholar

Damron, B. L., and Wilson, H. R. (1985). Bacitracin methylene disalicylate in broiler breeder feeds. Poult. Sci. 64, 1451–1454. doi: 10.3382/ps.0641451

PubMed Abstract | CrossRef Full Text | Google Scholar

DeLegge, M. (2008). Nutrition And Gastrointestinal Disease. Totowa: Humana Press.

Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/aem.03006-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Dziarski, R. (2004). Peptidoglycan recognition proteins (PGRPs). Mol. Immunol. 40, 877–886. doi: 10.1016/j.molimm.2003.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Gidenne, T., and Bellier, R. J. L. P. S. (2000). Use of digestible fibre in replacement to available carbohydrates: effect on digestion, rate of passage and caecal fermentation pattern during the growth of the rabbit. Livest. Prod. Sci. 63, 141–152. doi: 10.1016/s0301-6226(99)00123-2

CrossRef Full Text | Google Scholar

Groschwitz, K. R., and Hogan, S. P. (2009). Intestinal barrier function: molecular regulation and disease pathogenesis. J. Allergy Clin. Immunol. 124, 3–20. doi: 10.1016/j.jaci.2009.05.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, Y., Zhao, B., Zhu, J., Chen, Y., Pan, Y., Zhu, C., et al. (2016). Effect of Bacitracin Methylene Disalicylate on production performance and slaughter performance of the New Zealand white rabbit. Chin. J. Rabbit Farming 6, 4–7.

Google Scholar

Jing, X., Zulfiqar, F., Park, S. Y., Nunez, G., Dziarski, R., and Gupta, D. (2014). Peptidoglycan recognition protein 3 and Nod2 synergistically protect mice from dextran sodium sulfate-induced colitis. J. Immunol. 193, 3055–3069. doi: 10.4049/jimmunol.1301548

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, D., Liu, G., Lundstrom, A., Gelius, E., and Steiner, H. (1998). A peptidoglycan recognition protein in innate immunity conserved from insects to humans. Proc. Natl. Acad. Sci. U. S. A. 95, 10078–10082. doi: 10.1073/pnas.95.17.10078

PubMed Abstract | CrossRef Full Text | Google Scholar

Laukoetter, M. G., Nava, P., and Nusrat, A. (2008). Role of the intestinal barrier in inflammatory bowel disease. World J. Gastroenterol. 14, 401–407. doi: 10.3748/wjg.14.401

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Xu, Z., Gupta, D., and Dziarski, R. (2001). Peptidoglycan recognition proteins: a novel family of four human innate immunity pattern recognition molecules. J. Biol. Chem. 276, 34686–34694.

Google Scholar

Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

Noval Rivas, M., Burton, O. T., Wise, P., Zhang, Y. Q., Hobson, S. A., and Garcia Lloret, M. (2013). A microbiota signature associated with experimental food allergy promotes allergic sensitization and anaphylaxis. J. Allergy Clin. Immunol. 131, 201–212. doi: 10.1016/j.jaci.2012.10.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Sims, M. D., Dawson, K. A., Newman, K. E., Spring, P., and Hoogell, D. M. (2004). Effects of dietary mannan oligosaccharide, bacitracin methylene disalicylate, or both on the live performance and intestinal microbiology of turkeys. Poult. Sci. 83, 1148–1154. doi: 10.1093/ps/83.7.1148

PubMed Abstract | CrossRef Full Text | Google Scholar

Sokol, H., Leducq, V., Aschard, H., Pham, H. P., Jegou, S., Landman, C., et al. (2017). Fungal microbiota dysbiosis in IBD. Gut 66, 1039–1048. doi: 10.1136/gutjnl-2015-310746

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, R., Foster, D. N., and Shurson, G. C. (2011). Effects of feeding diets containing bacitracin methylene disalicylate to heat-stressed finishing pigs. J. Anim. Sci. 89, 1830–1843. doi: 10.2527/jas.2010-3218

PubMed Abstract | CrossRef Full Text | Google Scholar

Videhult, F. K., and West, C. E. (2016). Nutrition, gut microbiota and child health outcomes. Curr. Opin. Clin. Nutr. Metab. Care 19, 208–213.

Google Scholar

Videhult, F. K., West, C. E. J. C. O. I. C. N., and Care, M. (2016). Nutrition, gut microbiota and child health outcomes. Curr. Opin. Clin. Nutr. 19, 208–213.

Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. Appl. Environ. Microbiol. 73, 5264–5267.

Google Scholar

Werner, T., Liu, G., Kang, D., Ekengren, S., Steiner, H., and Hultmark, D. (2000). A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc. Natl. Acad. Sci. U. S. A. 97, 13772–13777. doi: 10.1073/pnas.97.25.13772

PubMed Abstract | CrossRef Full Text | Google Scholar

White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 5:e1000352. doi: 10.1371/journal.pcbi.1000352

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, G., Duff, G. C., Hall, L. W., Allen, J. D., Burrows, C. D., Bernal-Rigoli, J. C., et al. (2013). Alteration of digestive tract microbiome in neonatal Holstein bull calves by bacitracin methylene disalicylate treatment and scours. J. Anim. Sci. 91, 4984–4990. doi: 10.2527/jas.2013-6304

PubMed Abstract | CrossRef Full Text | Google Scholar

Yano, J. M., Yu, K., Donaldson, G. P., Shastri, G. G., Ann, P., Ma, L., et al. (2015). Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell 161, 264–276. doi: 10.1016/j.cell.2015.02.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, W., Roy, N. C., Lee, J., Lawley, B., Otter, D., Henderson, G., et al. (2013). Bowel microbiota moderate host physiological responses to dietary konjac in weanling rats. J. Nutr. 143, 1052–1060. doi: 10.3945/jn.113.174854

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J. C., Hale, V. L., Khodadadi, H., and Baban, B. (2019). Whole Body Vibration-Induced Omental Macrophage Polarization and Fecal Microbiome Modification in a Murine Model. Int. J. Mol. Sci. 20:3125. doi: 10.3390/ijms20133125

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaidman-Remy, A., Poidevin, M., Herve, M., Welchman, D. P., Paredes, J. C., Fahlander, C., et al. (2011). Drosophila immunity: analysis of PGRP-SB1 expression, enzymatic activity and function. PLoS One 6:e17231. doi: 10.1371/journal.pone.0017231

PubMed Abstract | CrossRef Full Text | Google Scholar