Rumen and cecum microbial dynamics following narasin inclusion in Nellore cattle diets
Johnny M. Souza, Leandro A. F. Silva, Daniel M. Casali, Janaina C. S. M. Souza, Luke A. Wolfe, Joseph H. Skarlupka, Ibrahim Zuniga-Chaves, Andrew J. Steinberger, Courtney L. Deblois, Andrew J. Scheftgen, Ana L. J. Lelis, Tiago Leiva, Murilo Chuba Rodrigues
TL;DR
This study shows that adding narasin to Nellore cattle diets selectively affects gut bacteria in the rumen and cecum without broadly changing microbial diversity.
Contribution
The study reveals compartment-specific microbial responses to narasin in both the rumen and cecum of beef cattle.
Findings
Narasin reduced Chao richness in the rumen but not Shannon diversity in either the rumen or cecum.
Narasin altered phylum-level composition, decreasing Firmicutes and increasing Bacteroidetes in the rumen.
In the cecum, narasin increased Firmicutes, especially Clostridiales, without affecting overall diversity.
Abstract
This study investigated the effects of narasin supplementation on the ruminal and cecal bacterial communities of feedlot Nellore cattle. We hypothesized that narasin would selectively modulate microbial populations in distinct gastrointestinal compartments without causing broad-scale disruption of overall community diversity. Sixty-four Nellore bulls (393 ± 24 kg) were assigned to a completely randomized block design and fed finishing diets containing either 0 or 20 ppm of narasin for 112 days. Rumen and cecal contents were collected at slaughter and analyzed using 16S rRNA gene sequencing to characterize bacterial community structure and composition. Overall, the rumen exhibited greater bacterial diversity and richness than the cecum, regardless of dietary treatment. Narasin supplementation did not affect Shannon diversity in either the rumen (p = 0.182) or the cecum (p = 0.298);…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9| Diets | A1 | A2 | A3 | T |
|---|---|---|---|---|
| Concentrate level (%) | 75.00 | 79.00 | 83.00 | 87.00 |
| Ingredients (% DM) | ||||
| Peanut Hull | 10.00 | 10.00 | 10.00 | 10.00 |
| Hay | 15.00 | 11.00 | 7.00 | 3.00 |
| Ground corn | 51.00 | 56.00 | 61.00 | 66.00 |
| Soybean meal | 10.00 | 8.00 | 6.00 | 4.00 |
| Cottonseed meal | 11.00 | 12.00 | 13.00 | 14.00 |
| Urea | 0.50 | 0.50 | 0.50 | 0.50 |
| Premix | 2.50 | 2.50 | 2.50 | 2.50 |
| Nutritional content | ||||
| Dry matter (DM) (%) | 89.00 | 84.00 | 84.00 | 84.00 |
| Total digestible nutrients (% DM) | 69.00 | 71.00 | 73.00 | 75.00 |
| Crude protein (% DM) | 14.70 | 14.70 | 14.60 | 14.70 |
| Neutral detergent fiber (% DM) | 33.20 | 30.20 | 27.20 | 24.10 |
| Non-fibrous carbohydrates (% DM) | 46.00 | 49.00 | 52.00 | 55.00 |
| Ether extract (% DM) | 2.80 | 2.90 | 3.00 | 3.00 |
| peNDF | 23.00 | 20.00 | 17.00 | 13.00 |
| NEL (Mcal/kg DM) | 1.01 | 1.05 | 1.11 | 1.17 |
| Ca (% DM) | 0.57 | 0.55 | 0.54 | 0.56 |
| P (% DM) | 0.39 | 0.40 | 0.41 | 0.42 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsRuminant Nutrition and Digestive Physiology · Oral microbiology and periodontitis research · Milk Quality and Mastitis in Dairy Cows
Introduction
1
Feedlot diets that are high in energy, particularly those rich in starch, often lead to an increase in short-chain fatty acids (SCFA) in the rumen of cattle. This can result in temporary drops in ruminal pH, which may cause metabolic issues such as ruminal acidosis, subsequently affecting cattle feed intake (Owens et al., 1998). Ruminal acidosis is one of the primary health challenges faced by dairy and feedlot cattle (Silvestre and Millen, 2021; Golder and Lean, 2024). Moreover, numerous studies have focused on preventing or mitigating acidosis to enhance the performance and efficiency of feedlot cattle (Owens et al., 1998; Brown et al., 2006).
Nellore cattle (Bos indicus) represent the predominant beef cattle breed in Brazil and account for a substantial proportion of beef production in tropical and subtropical regions worldwide. Compared with Bos taurus breeds, Nellore cattle exhibit distinct physiological and metabolic characteristics, including differences in feed efficiency, ruminal fermentation patterns, and host–microbiome interactions. Previous studies have demonstrated that both host genetics and breed type can influence ruminal microbial community composition and function, highlighting the importance of considering zebu cattle as a distinct biological model in microbiome research (Henderson et al., 2015). Despite their global relevance, most rumen microbiome studies have focused on taurine cattle, leaving important knowledge gaps regarding the microbial ecology of Nellore cattle under intensive feedlot conditions.
Feed additives, such as ionophores, have been investigated worldwide over the past couple of decades to control ruminal acidification and enhance the performance of feedlot cattle (Silvestre et al., 2023). The primary mechanism of most ionophores is to regulate the growth of gram-positive rumen bacteria, including those that produce lactate in the rumen (Nagaraja and Titgemeyer, 2007). Ionophores alter transmembrane ion gradients in susceptible microorganisms, resulting in shifts in microbial activity and fermentation end-products, such as increased propionate production and improved feed efficiency (Russell and Strobel, 1989; Russell and Houlihan, 2003).
Narasin is one of the ionophores that has been researched for its ability to alter rumen fermentation dynamics (Miszura et al., 2018), resulting in changes in plasma metabolites by increasing glucose levels (Sardinha et al., 2020) and reducing urea concentration (Polizel et al., 2020), which leads to positive results in animal performance (Silva et al., 2015; Polizel et al., 2020). In this context, Limede et al. (2021) reported an increase in ruminal propionate concentration, a reduction in the ruminal acetate-to-propionate ratio, and an increase in total SCFA concentration in response to the inclusion of narasin in the diet of Nellore cattle. Concerning performance, the same authors found that narasin enhanced dry matter intake, average daily gain, and final live weight compared to the control treatment, demonstrating the potential of this ionophore for animal health performance. Although several studies have documented the effects of narasin and other ionophores on ruminal fermentation parameters, information regarding their direct impact on ruminal microbial community structure remains limited and is largely inferred from studies with monensin (Marques and Cooke, 2021). Moreover, the potential effects of narasin on the cecal microbiota—a key site for post-ruminal fermentation—remain largely unexplored, despite evidence that microbial communities differ markedly along the gastrointestinal tract (Wang et al., 2017).
However, previous studies assessing the inclusion of narasin in cattle diets did not examine its effects on bacterial communities in the rumen and cecum. Clarifying narasin’s mode of action is essential to support the reasons for improved performance in cattle consuming this feed additive. Therefore, this study aims to identify bacterial communities through 16S rRNA gene sequencing in the rumen and cecum of Nellore cattle finished in a feedlot on diets containing narasin.
Given the functional differences between the rumen and the cecum and the limited information regarding the effects of narasin on microbial communities across gastrointestinal compartments, a comprehensive evaluation of both sites is warranted. Previous studies have shown that microbial composition and interactions vary substantially along the gastrointestinal tract, reflecting differences in substrate availability and fermentation dynamics (Wang et al., 2017; Malmuthuge and Guan, 2017). We hypothesized that narasin would exert selective pressure on specific bacterial taxa in both compartments, leading to compartment-specific shifts in microbial composition without causing broad-scale disruption of overall community diversity. Therefore, this study aimed to characterize the ruminal and cecal bacterial communities of feedlot Nellore cattle supplemented with narasin using 16S rRNA gene sequencing.
Material and methods
2
Animals and experimental design
2.1
Animal care and handling procedures used in this experiment adhered to the guidelines of the Animal Use Ethics Committee (CEUA) and were approved by the Ethics Committee on Animal Use of São Paulo State University (UNESP), Dracena campus, Brazil (Protocol CEUA 002/2022).
This study was conducted at the experimental feedlot of the Faculty of Agricultural and Technological Sciences at UNESP, Dracena campus, Brazil. Sixty-four Nellore bulls (393 ± 24 kg) were utilized and randomly assigned to 16 pens (n = 4 animals per pen) in a completely randomized block design, with initial live weight serving as the blocking criterion. The study lasted 112 days, after which all animals were slaughtered in a commercial setting slaughterhouse.
The animals were pre-adapted to the facility and were fed Tifton hay for 7 days to standardize the ruminal microbiota population. Subsequently, the animals were randomly allocated within each block to one of the following treatments: 0 ppm (Control) or 20 ppm of Narasin (Elanco Animal Health, São Paulo, SP, Brazil). The additive was provided according to the manufacturer’s instructions. Narasin supplementation was initiated at the beginning of the adaptation period and was maintained continuously throughout the entire experimental period, including the finishing phase, until slaughter.
Management, feeding, and animal care
2.2
Upon arrival, the animals were weighed, treated with anthelmintics, and vaccinated against viral and bacterial diseases, including rotavirus, coronavirus, tetanus, botulism, and seven types of Clostridium spp. (Cattlemaster, Pfizer Animal Health). Animals were housed in partially concreted pens, with a stocking density of four animals per pen (18 m^2^ and 1.5 m of linear feed bunk per animal), with unrestricted access to water via automatic drinkers.
Feeding occurred three times a day: at 08:00 (35% of the total), 11:00 (20% of the total), and 16:00 (45% of the total). The study included a 14-day adaptation period following a step-up protocol with three adaptation diets: adaptation 1 (75% concentrate) for 5 days, adaptation 2 (79% concentrate) for 4 days, and adaptation 3 (83% concentrate) for 5 days. The finishing diet comprised 87% concentrate ingredients (Table 1). Experimental diets were formulated according to the Large Ruminant Nutrition System (LRNS; Fox et al., 2004), level 2. Diets were provided ad libitum and adjusted daily based on feed bunk refusals, which were maintained at 5% of the offered amount.
After 112 days in the feedlot, animals were slaughtered at a commercial slaughterhouse, reaching an average final live weight of approximately 520 kg. The contents of the rumen and cecum were collected post-slaughter from all animals and stored at −80 °C for subsequent bacterial community analysis.
DNA extraction, PCR amplification, and library construction
2.3
Total DNA was extracted through mechanical disruption and a phenolic extraction protocol (Stevenson and Weimer, 2007), with a modification of the PCSA method (phenol:chloroform with bead beating II), as previously described (Henderson et al., 2013), which have proven effective for obtaining high-yield DNA that accurately represents ruminal microbial communities.
Briefly, sample preparation for DNA extraction was initiated using a stomacher apparatus (Seward Stomacher 400; Seward Inc., West Sussex, UK), which employs dual paddles to mechanically agitate samples and detach fiber-associated microbial cells from solid particles while simultaneously fragmenting them into smaller components. Samples (10 mL by volume) was combined with 40 mL of 1 × Phosphate-Buffered Saline (PBS) within stomacher bags with 0.5 mm filter insert (BA6040/CLR/STR; Seward Inc.). After 5 min, the corner of the bag was cut, and the liquid was removed. The resulting buffer suspension underwent centrifugation at 5,000 × g for 1 h to concentrate the cells. Subsequently, the supernatant was discarded, and the resulting pellet was reconstituted in 5 mL of PBS. Then, samples underwent treatment with extraction buffer containing Tris–HCl (0.1 M), EDTA (0.01 M), and NaCl (0.15 M). The extraction buffer served as a negative control throughout the DNA isolation procedure. Duplicate aliquots of 1 mL from the reconstituted sample were transferred to microcentrifuge tubes containing 700 μL of equilibrated phenol (Sigma-Aldrich, St. Louis, MO), 50 μL of 20% sodium dodecyl sulfate solution (Sigma-Aldrich, St. Louis, MO), and 0.5 g of zirconium/silica beads (0.1 mm diameter, BioSpec Products, Bartlesville, OK). Cellular disruption was achieved through bead-beating methodology, conducted in two 2-min intervals with an intermediate 10-min incubation at 60 °C between beating cycles. Following mechanical lysis, tubes were subjected to centrifugation at 14,000 × g for 10 min, and DNA recovery was performed from 850 μL of the resulting supernatant. DNA purification was accomplished through sequential vortexing of the supernatant with equivalent volumes of phenol-chloroform-isoamyl alcohol [PCI, 25:24:1 (vol/vol/vol)], followed by centrifugation at 14,000 × g for 10 min. DNA precipitation was achieved by gentle mixing of 500 μL supernatant with 50 μL of 2 M sodium acetate and 300 μL isopropyl alcohol. The resulting mixture underwent incubation for 2 h or overnight at −20 °C, followed by centrifugation at 14,000 × g for 20 min at 4 °C. The final step involved washing the pellet with 1 mL of ice-cold 70% ethanol and air-drying for 90 min or overnight before reconstitution in elution buffer (10 mM Tris and 1 mM EDTA).
DNA concentrations were determined using Qubit fluorometer reagents (Invitrogen, Waltham, MA). Amplification of the 16S rRNA gene V4 hypervariable region was conducted through polymerase chain reaction (PCR) employing universal bacterial primers (515F: 5’-GTGCCAGCMGCCGCGGTAA-3′; 806R: 5’-GGACTACHVGGGTWTCTAAT-3′), according to Kozich et al. (2013). The primer sequences incorporated adapters designed for compatibility with Illumina sequencing technology (F- AATGATACGGCGACCACCGAGATCTACAC; R- CAAGCAGAAGACGGCATACGAGAT) and contained distinct barcoding sequences to enable sample multiplexing: forward primers carried 16 distinct 8-bp barcodes, while reverse primers contained 24 distinct 8-bp barcodes.
PCR amplification reactions were assembled by combining 25–50 ng of template DNA with 0.2 μmol/L of each primer in a final reaction volume of 25 μL containing 2 × KAPA HiFi HotStart ReadyMix (KAPA Biosystems, Wilmington, MA). Thermal cycling was performed using a Bio-Rad S1000 thermocycler (Bio-Rad Laboratories, Hercules, CA, USA) under the following parameters: initial denaturation at 95 °C for 3 min, followed by 25 amplification cycles consisting of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, with a final extension step at 72 °C for 5 min.
Amplification products were analyzed by electrophoresis on 1% (w/v) low-melting-point agarose gels prepared with AquaPor low-melt agarose (National Diagnostics, Atlanta, GA) and visualized using SYBRSafe DNA gel stain (Invitrogen, Waltham, CA). Successful amplification was confirmed by the presence of bands at approximately 380 bp. Target bands were excised from gels, recovered, and purified using the Zymoclean Gel DNA Recovery Kit (Zymo Research, Irving, CA). No-template negative controls were incorporated into each PCR batch, and any detection of amplification in negative controls resulted in complete re-processing of all samples in that batch beginning from PCR setup. For negative controls showing no amplification, the corresponding gel region (~380 bp) was excised and processed for sequencing to provide additional confirmation of contamination absence.
Purified DNA from gel extraction was quantified using a Qubit fluorometer and 96-well plate spectrophotometer. Library construction involved creating a 4 nmol/L equimolar pool of all PCR products. The resulting library underwent sequencing on an Illumina MiSeq platform following standard Illumina protocols, utilizing a MiSeq v2 2 × 250 sequencing kit at 10 pmol/L concentration with 10% PhiX control. All sequences associated with this study were deposited into the National Center for Biotechnological Information’s Sequence Read Archive and is available under BioProject accession number PRJNA1307616.
Sequence processing and analysis
2.4
Raw fastq files generated by the sequencer underwent quality control and bioinformatics analysis using mothur v1.39.0 (Schloss et al., 2009). Initial sequence screening was performed (maxambig = 0, maxhomop = 8, minlength = 200, maxlength = 500), followed by grouping of identical sequences using the unique.seqs command. Sequences were subsequently aligned (align.seqs) against the SILVA 16S rRNA gene reference alignment database (Release 132, Pruesse et al., 2007) and subjected to additional screening to retain only sequences aligned to the target region (screen.seqs, start = 13,862, end = 23,444). Following alignment, sequences underwent filtering (filter.seqs) and re-grouping of identical sequences (unique.seqs). Highly similar sequences were consolidated using the pre.cluster command (diffs = 2), while chimeric sequences were identified and eliminated using chimera.uchime and remove.seqs commands, respectively. Taxonomic assignment was conducted using the SILVA database, with subsequent removal of non-bacterial sequences (those unclassified at the Kingdom level, Archaea, Eukaryota, cyanobacteria, and mitochondria).
Singleton sequences (appearing only once in the dataset) were eliminated (split.abund) to reduce bias from sequencing errors, and uncorrected pairwise distances between sequences were computed (dist.seqs). Sequence clustering into operational taxonomic units (OTUs) was performed using cluster.split (method = opti, cutoff = 0.03) at 97% sequence similarity threshold. Sample coverage was assessed using Good’s index (Good, 1953), while OTU taxonomic classification was determined using the GreenGenes database (August 2013 release) (DeSantis et al., 2006). Sample normalization was conducted by subsampling to the minimum sequence count observed across all samples, which corresponded to 9,100 sequences. A 0.1% total abundance threshold was implemented on the normalized dataset using filter.shared (minpercent = 0.001). Community structure analysis was subsequently performed using non-metric multidimensional scaling (NMDS) based on Bray–Curtis dissimilarity index (Bray and Curtis, 1957).
Microbial network analysis was performed to identify co-occurrence patterns among bacterial genera. To investigate the structure of microbial interactions within the ruminal and cecal compartments, co-occurrence networks were constructed at the genus level for both the Control and Narasin treatment groups. Spearman’s rank correlation coefficients (ρ) were calculated for all pairwise comparisons of genera. Significant correlations, defined as those with an absolute coefficient |ρ| > 0.3 and a p-value < 0.05, were used to construct the networks. Network visualization and analysis were performed using the igraph package in R. In the resulting graphs, nodes represent genera and edges represent significant correlations. Nodes were sized proportionally to the mean relative abundance of each genus, and colored according to their respective phylum. Edges were colored to represent positive (green) or negative (red) correlations. Finally, network centrality metrics (degree, closeness, and betweenness) were calculated for each node to characterize its topological importance within the microbial community.
Statistical analysis
2.5
All statistical analyses were performed in R version 4.3.3 (R Core Team, 2025). Prior to statistical analysis, all data were tested for normality using the Shapiro–Wilk test and homogeneity of variances using Levene’s test. Bacterial sequences were grouped into operational taxonomic units (OTUs) with 97% sequence similarity. The coverage of Good (1953) was calculated in Mothur (v.1.48.1) for all samples. All samples achieved Good’s coverage ≥93% (mean = 0.94). The OTU counts were normalized to 10,000 sequences per sample, and these normalized OTU counts were used for further analysis.
Alpha diversity was assessed using both Chao’s Richness estimate (Chao, 1984) and Shannon’s Diversity Index (Shannon, 2001), a measure of richness and evenness. Alpha diversity metrics were calculated from normalized OTU count data using the vegan package v.2.6–10 in R v.4.3.3. Shannon’s diversity index was computed using the diversity() function, while species richness was estimated using the Chao1 estimator. Statistical comparisons between treatment groups were performed separately for each compartment using Student’s t-test. Statistical significance was declared at p ≤ 0.05.
Beta diversity was assessed using Bray-Curtis (Bray and Curtis, 1957) and Jaccard (Jaccard, 1912) dissimilarity matrices calculated with vegdist(). Non-metric multidimensional scaling (NMDS) ordination was performed using metaMDS() with k = 2 dimensions and trymax = 100 iterations. Statistical significance of treatment effects was evaluated using PERMANOVA implemented through adonis2(). Ordination plots display 95% confidence ellipses and stress values for interpretation quality.
Taxonomic composition analysis was performed at the phylum level using phyloseq package (McMurdie and Holmes, 2013). The seven most abundant phyla were selected using taxa_sums(). For visualization, samples were grouped by treatment using merge_samples(), followed by relative abundance calculation using transform_sample_counts(). Phyla <1% abundance were grouped as ‘<1% abundance’. Statistical differences in phylum-level abundances between treatment groups were assessed using the Wilcoxon rank-sum test (Mann–Whitney U test) for each individual phylum, comparing control versus narasin-treated animals. p-values were adjusted for multiple comparisons using the false discovery rate (FDR) method (Benjamini and Hochberg, 1995). Phyla with FDR-adjusted p-values < 0.05 were considered significantly different between treatments.
Results
3
The results for the Shannon diversity index and Chao species richness are presented in Figure 1. Overall, the rumen demonstrated greater diversity and richness compared to the cecum, regardless of treatment. No differences were observed in the Shannon index for either the rumen (p = 0.182) or the cecum (p = 0.298). However, narasin reduced Chao richness only in the rumen (p = 0.032), not in the cecum (p = 0.274).
Shannon’s diversity index and Chao richness index of rumen and cecum bacteria from Nellore cattle fed narasin.
Beta-diversity analysis using NMDS ordination revealed no significant differences in overall microbial community structure between control and narasin-treated groups (p > 0.201), either the rumen or the cecum (Figure 2). The ordination provided reliable representation of community patterns (stress values: 0.091–0.11), indicating that while individual taxa may be affected by narasin treatment, the overall community structure remained stable. These results suggest that narasin exerts selective effects on specific bacterial taxa rather than causing broad-scale community disruption.
Difference of beta diversity in the two groups using principal coordinate analysis (PCoA) based on Bray-Curtis distance and Jaccard in rumen and cecum content from Nellore cattle fed narasin.
Regarding relative abundance in the rumen (Figure 3), there was a reduction in the abundance of the phylum Firmicutes, accompanied by an increase in Bacteroidetes in cattle-fed narasin (p = 0.05). Regarding relative abundance in the cecum (Figure 4), the phylum Firmicutes was predominant in the cecum, contrasting sharply with the rumen environment. Furthermore, there was a significant effect (p = 0.05) of narasin on the cecal microbiota abundance, with supplementation increasing the abundance of Firmicutes, particularly within the Clostridiales order (Figure 5). The heat map of microbial abundance for the main genera in the rumen and cecum is shown in Figures 6, 7, respectively. No significant differences due to narasin were observed in the rumen and cecum (p > 0.05).
Relative abundance of rumen bacteria, at phylum level, in Nellore cattle fed narasin. Values bearing different superscripts (a, b) in a row, differ significantly (p = 0.05).
Relative abundance of cecum bacteria, at phylum level, in Nellore cattle fed narasin. Values bearing different superscripts (a, b) in a row, differ significantly (p = 0.05).
Relative abundance of cecum bacteria (Clostridiales) in Nellore cattle fed narasin.
Heatmap of microbial abundance showing the main genera in the rumen of feedlot Nellore cattle fed narasin.
Heatmap of microbial abundance showing the main genera in the cecum of feedlot Nellore cattle fed narasin.
To investigate the impact of narasin on ruminal microbial interactions, we constructed genus-level co-occurrence networks for the Control and Narasin groups (Figure 8). The network from the Control group displayed greater complexity, featuring a total of 40 significant correlations, of which 21 (52.5%) were positive and 19 (47.5%) were negative. This network comprised 27 distinct genera, representing 10 different phyla. In contrast, narasin supplementation resulted in a more simplified interaction network. The Narasin group network exhibited only 30 significant correlations, a 25% reduction in the total number of connections compared to the Control group. In this network, a slight predominance of negative interactions was observed, with 16 correlations (53.3%) being negative and 14 (46.7%) positive. Despite the lower connectivity, the diversity of genera (28) and phyla (10) involved in the network was similar to that of the Control group.
Co-occurrence networks of the ruminal microbiota in control and narasin-treated animals. Genus-level interaction networks based on 16S rRNA gene sequence data, displaying statistically significant correlations (p < 0.05) with an absolute Spearman’s |ρ| > 0.3. Each node represents a bacterial genus, with node size scaled according to its mean relative abundance within each treatment group. Node color indicates the phylum. An edge represents a significant correlation between two genera: green edges indicate a positive correlation, and red edges indicate a negative correlation. Control network – 40 correlations (52.5% positive) 27 genera from 10 phyla; narasin network – 30 correlations (46.7% positive) 28 genera from 10 phyla.
Regarding the cecum, the network of the Control group was characterized by 19 significant correlations among 15 genera from 5 different phyla (Figure 9). This network was dominated by positive interactions, which accounted for 15 of the correlations (78.9%), while 4 were negative (21.1%). In contrast, the cecal network in the Narasin-treated group showed a modest but notable increase in connectivity, with 23 significant correlations among 23 genera from 5 phyla. This network exhibited an even stronger positive bias, with 19 positive correlations (82.6%) and 4 negative correlations (17.4%).
Co-occurrence networks of the cecal microbiota in Control and Narasin-treated animals. Genus-level interaction networks based on 16S rRNA gene sequence data, displaying statistically significant correlations (p < 0.05) with an absolute Spearman’s |ρ| > 0.3. Each node represents a bacterial genus, with node size scaled according to its mean relative abundance within each treatment group. Node color indicates the phylum. An edge represents a significant correlation between two genera: green edges indicate a positive correlation, and red edges indicate a negative correlation. The networks for the Control (left) and Narasin (right) groups are shown, illustrating the dramatic increase in network density and the shift to an exclusively positive interaction structure following narasin supplementation. Control Network – 19 correlations (78.9% positive) 15 genera from 5 phyla; Narasin Network – 23 correlations (82.6% positive) 23 genera from 5 phyla.
Discussion
4
The cecum plays a crucial role in post-ruminal fermentation, yet it remains an underexplored aspect of beef cattle microbiome research. Most studies focus on the rumen because of its primary function in fiber degradation and fermentation (Henderson et al., 2015). However, the cecum, although smaller in volume, can still influence overall digestive efficiency and microbial metabolite production (Li et al., 2019). Alterations in cecal microbiota composition may have physiological implications, potentially affecting nutrient absorption, hindgut fermentation, and even rumen function through enterohepatic circulation (Malmuthuge and Guan, 2017). Thus, understanding the impact of dietary additives such as narasin on the cecum is essential for a comprehensive evaluation of their effects on ruminant digestion and health.
Given the well-established selective mode of action of ionophores, the absence of broad-scale shifts in microbial diversity was expected, with narasin primarily exerting targeted effects on specific bacterial populations. The lack of significant differences in the Shannon diversity index suggests that narasin did not affect the evenness of the microbial communities in either the rumen or the cecum. This aligns with findings from earlier research indicating that ionophores typically modify specific bacterial populations rather than cause widespread diversity shifts (Russell and Strobel, 1989; McCann et al., 2017). However, the observed reduction in Chao richness in the rumen indicates a decrease in species diversity, which aligns with the known inhibitory effects of ionophores on Gram-positive bacteria (Callaway et al., 2003). This reduction is likely driven by the selective suppression of Firmicutes, a major bacterial group in the rumen, which comprises key fiber-digesting species (Russell and Houlihan, 2003). Moreover, narasin did not significantly affect beta-diversity in either the rumen or the cecum. The NMDS ordination provided reliable representation of community patterns (stress values: 0.091–0.11), indicating that while individual taxa may be affected by narasin treatment, the overall community structure remained stable. These results suggest that narasin exerts selective effects on specific bacterial taxa rather than causing broad-scale community disruption, consistent with the known selectivity of ionophores toward Gram-positive bacteria (Russell and Strobel, 1989; Russell and Houlihan, 2003; Marques and Cooke, 2021).
The decrease in Firmicutes and the increase in Bacteroidetes in the rumen of cattle-fed narasin is consistent with the mode of action of ionophores, which preferentially target Gram-positive bacteria by disrupting their ion gradients (Russell and Houlihan, 2003). Given that narasin’s mechanism of action is similar to monensin, a higher relative abundance of Bacteroidetes (Gram-negative) and a lower relative abundance of Firmicutes (Gram-positive) are expected in steers fed narasin. This shift is anticipated, considering that Bacteroidetes are Gram-negative and generally more resistant to ionophores (Mao et al., 2016). However, the absence of significant differences in other phyla suggests that functional redundancy within the microbial community may compensate for compositional shifts, enabling the rumen microbiota to sustain its overall fermentative capacity (Henderson et al., 2015).
Interestingly, a different pattern was observed in the cecal microbiota abundance. First, the phylum Firmicutes dominated in the cecum, which is quite distinct from the rumen environment. This result is consistent with the fact that the rumen serves as the primary fermentation chamber for forage, so little fiber reaches the cecum for fermentation, which may explain the low abundance of Bacteroidetes there (Deusch et al., 2017). On the other hand, due to the rumen liquid passage rate, a significant amount of concentrate can reach the cecum and be fermented, potentially explaining the dominance of Firmicutes in that area (Mao et al., 2016).
Moreover, the cecum exhibited a distinct response to narasin supplementation. Unlike the rumen, where Firmicutes showed a reduction, the cecum demonstrated a significant increase in Firmicutes abundance, particularly within the Clostridiales order. This finding is novel, as few studies have examined the effects of ionophores on the cecal microbiota of ruminants. The predominance of Firmicutes in the cecum is consistent with differences in substrate availability and fermentation dynamics. While the rumen primarily ferments fiber, the cecum ferments residual starch and soluble carbohydrates that escape ruminal digestion (Deusch et al., 2017). The increase in Clostridiales suggests a shift toward microbial populations that may be better adapted to utilizing these substrates under ionophore influence, though functional validation would be required to confirm this hypothesis. Thus, these results demonstrate that narasin supplementation induces specific changes in the cecal microbiota, the functional significance of which warrants further investigation.
Based on existing literature, cecal microbiota modifications have the potential to extend beyond local fermentation processes. There is increasing evidence indicating that changes in hindgut microbial composition can affect systemic metabolism and immune function through microbial metabolites such as SCFA and secondary bile acids (Malmuthuge and Guan, 2017; Huws et al., 2018). Additionally, changes in the cecal microbiome could potentially affect rumen function through enterohepatic recycling of microbial metabolites and microbial cross-talk between gut compartments (Mao et al., 2016). However, this study did not assess metabolite production or functional outcomes, and future research should investigate whether shifts in cecal microbial populations align with functional changes that may have systemic implications.
Our co-occurrence network analysis revealed that narasin supplementation induced a significant restructuring of the ruminal microbial community architecture. The 25% reduction in total interactions within the narasin group (from 40 to 30 edges) suggests that the ionophore acts not only by altering the abundance of specific taxa but also by simplifying the ecosystem by reducing the interdependence among genera. This loss of connectivity is a classic indicator of an ecosystem responding to a chemical stressor, potentially leading to the destabilization of established ecological niches, affecting both synergistic (cooperative) and antagonistic (competitive) relationships (Shade et al., 2012). Such simplification of microbial networks has been previously observed in ruminal environments under dietary pressures, suggesting a shift toward a less resilient but more specialized community (Hernandez-Sanabria et al., 2012). The shift in the balance of correlations, from a slight positive dominance in the control group (52.5%) to a slight negative dominance in the narasin group (53.3%), is particularly noteworthy. This may imply that under the selective pressure of narasin, competitive or exclusionary relationships become more prevalent or detectable compared to mutualistic ones. Ionophores like narasin selectively target Gram-positive bacteria by disrupting their transmembrane ion gradients, an action that can open ecological niches for more resistant microorganisms to establish new competitive relationships for substrates (Russell and Strobel, 1989; Callaway et al., 2003). Analyzing specific interactions highlights these dynamics. In the Control group, the strong positive correlation (r = 0.90) between Clostridium and Coprococcus, both members of the Firmicutes phylum, suggests a cooperative relationship, likely linked to fiber degradation and butyrate production (Pryde et al., 2002). This specific interaction was absent in the narasin network, indicating that narasin may have disrupted this synergy.
In the cecum, the microbial interaction network of the narasin-treated group was different from that observed in the rumen. Rather than simplifying the network, oral narasin supplementation was associated with a moderate increase in microbial connectivity and a shift toward more cooperative interactions. The increase in total interactions (from 19 to 23) and the enhanced positive correlation structure (78.9 to 82.6% positive) indicate the establishment of a more cooperative and potentially more stable community (Coyte et al., 2015). We hypothesize that this is an indirect effect, driven by the primary action of the ionophore in the upper gastrointestinal tract. By altering ruminal fermentation, narasin likely modifies the composition of digesta flowing into the cecum (Duffield et al., 2012). This ruminal outflow, containing different microbial cells, metabolites, and fermentation byproducts, would create a new ecological foundation for the hindgut community. This altered environment appears to have selected for a highly specialized consortium of bacteria, which is consistent with our finding that narasin supplementation increased the abundance of Firmicutes, particularly within the Clostridiales order. Strong positive correlations were maintained among key cecal fermenters, including Coprococcus, Oscillospira, and other members of the Clostridiales order were observed. This pattern suggests that narasin selectively modulates specific microbial interactions rather than uniformly promoting cooperation, resulting in a more structured and potentially more efficient cecal ecosystem (Berry and Widder, 2014). Moreover, the narasin treatment was associated with the emergence of specific competitive interactions, particularly negative correlations involving Clostridium with both Lactobacillus and Bifidobacterium, suggesting altered competitive dynamics within the Firmicutes phylum. This finding is consistent with well-documented antagonistic relationships between these genera in the intestinal environment. Bifidobacterium species have been shown to exert direct inhibitory effects against Clostridium through multiple mechanisms, including acidification of the intestinal environment, competition for nutrients and adhesion sites, and production of antimicrobial compounds (Wei et al., 2018). Similarly, Lactobacillus strains demonstrate competitive exclusion against various Clostridium species, particularly through the production of organic acids and bacteriocins that create unfavorable conditions for clostridial growth (Li et al., 2017). This difference in network complexity between rumen and cecum likely reflects the distinct functional roles of these compartments, with the rumen serving as the primary fermentative chamber processing diverse dietary components, while the cecum represents a more specialized environment focused on secondary fermentation (Wang et al., 2017; Mi et al., 2018).
The present findings contribute to the limited body of literature on cecal microbiota in ruminants and highlight the necessity for further investigation into its role in digestive physiology. To establish causal relationships between taxonomic shifts and functional outcomes, future studies should incorporate metabolomic and transcriptomic approaches to clarify the functional consequences of microbial shifts induced by dietary interventions such as narasin. Understanding the full extent of microbial interactions across gut compartments is crucial for optimizing feed strategies that enhance both rumen efficiency and overall animal health.
Conclusion
5
Narasin appears to exert a selective effect on specific bacterial populations, such as Firmicutes, particularly within the Clostridiales order, without significantly altering overall microbial diversity. This effect may be linked to reduced richness in the rumen and shifts in the cecal microbiota, potentially impacting feed efficiency by favoring certain fermentation pathways. Moreover, our results demonstrate that narasin remodels the ruminal microbial ecology by reducing the complexity of inter-genus interactions. This network reconfiguration, characterized by fewer connections and a shift toward negative correlations, represents a fundamental mechanism explaining the ionophore’s effects on the rumen environment. However, further studies are necessary to understand the functional implications of these microbial shifts and how they translate into production benefits.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Benjamini Y. Hochberg Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B, 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb 02031.x · doi ↗
- 2Berry D. Widder S. (2014). Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. 5:219. doi: 10.3389/fmicb.2014.00219, 24904535 PMC 4033041 · doi ↗ · pubmed ↗
- 3Bray J. R. Curtis J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268 · doi ↗
- 4Brown M. S. Ponce C. H. Pulikanti R. (2006). Adaptation of beef cattle to high-concentrate diets: Performance and ruminal metabolism. Journal of Animal Science, 84, E 25–E 33. doi: 10.2527/2006.8413_suppl E 25x 16582090 · doi ↗ · pubmed ↗
- 5Callaway T. R. Edrington T. S. Rychlik J. L. Nisbet D. J. (2003). Ionophores: their use as ruminant growth promotants and impact on food safety. Curr. Issues Intest. Microbiol. 4, 43–51.14503688 · pubmed ↗
- 6Chao A. (1984). Nonparametric estimation of the number of classes in a population. Scandinavian Journal of Statistics, 11, 265–270.
- 7Coyte K. Z. Schluter J. Foster K. R. (2015). The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666. doi: 10.1126/science.aad 2602, 26542567 · doi ↗ · pubmed ↗
- 8De Santis T. Z. Hugenholtz P. Larsen N. Rojas M. Brodie E. L. Keller K. . (2006). Greengenes, a chimera-checked 16S r RNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05, 16820507 PMC 1489311 · doi ↗ · pubmed ↗