Characteristics of Bacterial Microbiota in Different Intestinal Segments of Aohan Fine-Wool Sheep

The microbial community performs vital functions in the intestinal system of animals. Modulation of the gut microbiota structure can indirectly or directly affect gut health and host metabolism. Aohan fine-wool sheep grow in semi-desert grasslands in China and show excellent stress tolerance. In this study, we amplified 16S rRNA gene to investigate the dynamic distribution and adaptability of the gut microbiome in the duodenum, jejunum, ileum, cecum, colon, and rectum of seven Aohan fine-wool sheep at 12 months. The results showed that the microbial composition and diversity of the ileum and the large intestine (collectively termed the hindgut) were close together, and the genetic distance and functional projections between them were similar. Meanwhile, the diversity index results revealed that the bacterial richness and diversity of the hindgut were significantly higher than those of the foregut. We found that from the foregut to the hindgut, the dominant bacteria changed from Proteobacteria to Bacteroidetes. In LEfSe analysis, Succiniclasticum was found to be significantly abundant bacteria in the foregut and was involved in succinic acid metabolism. Ruminococcaceae and Caldicoprobacteraceae were significantly abundant in hindgut, which can degrade cellulose polysaccharides in the large intestine and produce beneficial metabolites. Moreover, Coriobacteriaceae and Eggthellaceae are involved in flavonoid metabolism and polyphenol production. Interestingly, these unique bacteria have not been reported in Mongolian sheep or other sheep breeds. Collectively, the gut microbiota of Aohan fine-wool sheep is one of the keys to adapting to the semi-desert grassland environment. Our results provide new insights into the role of gut microbiota in improving stress tolerance and gut health in sheep.

Introduction

Every part of an individual animal has a microbial community (Neish, 2009). The microbial community exists in a symbiotic relationship with the host (Ley et al., 2006a). This complex collection of microorganisms is called the microbiota, and their genetic material is called the microbiome (Turnbaugh et al., 2007). Animal intestinal microbes are dynamic; during the host’s life cycle, the microflora undergo significant changes due to factors, such as diet, environment, and disease state, among others (Shreiner et al., 2015). Microbes contribute to energy homeostasis, metabolism, intestinal epithelial health, immune activity, and nerve development in animals (Cho and Blaser, 2012). In addition, intestinal microorganisms play an important role in the development of animals. The addition of corn bran to the diet of weaned piglets induced changes in the ratio of Firmicutes to Bacteroidetes changing the microbial diversity and enhancing the anti-inflammatory response of the organism (Liu et al., 2018). Similarly, the addition of caragana to sheep diets changes the composition of intestinal microbes in the body, improving the meat quality of sheep in terms of tenderness and fatty acid content (Zhang et al., 2021).

At present, it is generally considered that the distribution of intestinal bacteria in different intestinal segments is different, and is related to the function of the particular intestinal segment. In a study on broiler chickens, Firmicutes was the dominant genus in the intestinal tract, and Bacteroidetes, as an important type of phylum, only occupies 50% of the abundance in the cecum (Xiao et al., 2017). A study on pigs revealed that Proteobacteria was the dominant phyla in the small intestine, while Firmicutes was the dominant phylum in the large intestine, accounting for 80% of the microbial population, followed by Proteobacteria (Zhao et al., 2015). In Mongolian sheep growing in the Gansu province of China, the proportion of Firmicutes and Bacteroidetes among the intestinal microbes reached 80%, and there was no consistent bacterial ratio among intestinal segments (Zeng et al., 2017). However, the intestinal microbial composition of a camel’s ileum, cecum, and colon was relatively similar, and Firmicutes and Bacteroidetes accounted for 50–60% of the intestinal microbiota (He et al., 2018). This shows that there are obvious differences in intestinal microbes between monogastric animals and ruminants. Meanwhile, different types of ruminants have different compositions of intestinal microbes. These differences are caused by the characteristics of the animal and its natural environment (Guo et al., 2020b). In our analysis, we considered the duodenum and jejunum as the foregut, and the ileum and large intestine as the hindgut based on the microbial composition.

The current analysis of sheep intestinal microbes focuses on the exploration of fecal microbes and microbes in individual intestinal segments. Some studies have reported that the fecal microbiota cannot represent the microbial composition of intestinal segments (Zhao et al., 2015; Donaldson et al., 2016). Moreover, for the study of a single intestinal segment, the correlation between the intestinal segments is usually lacking, and only the characteristics of a single intestinal segment can be displayed (Al-Masaudi et al., 2017; Zhang et al., 2018; Lin et al., 2021). To study the composition of the gut microbiome of Aohan fine-wool sheep, we determined the complete microbiome spectrum of the adult sheep gut. We also systematically studied the changes in the intestinal microbiota between the foregut and hindgut to explore the stability of microorganisms in different intestinal segments.

Aohan fine-wool sheep is produced in Northeast China. It is a wool-meat dual-purpose breed developed by crossing Chinese Mongolian sheep as female parent and Soviet Caucasian sheep and Gustav sheep as male parent. Aohan fine-wool sheep have the advantages of roughage tolerance and high adaptability (Xiao-ping, 2009; Cui et al., 2014). Aohan sheep farms are located in regions with low annual rainfall, low winter temperatures, and strong winds accompanied by sandstorms and dust problems. Therefore, the quantity and quality of pastures in Aohan sheep farms are poor (Masters et al., 1990). Aohan fine-wool sheep have been selected and bred for generations, becoming an excellent breed that can be grown in arid desert areas, for wool and meat. In this study, we analyzed the overall gut microbial composition and differential bacterial genera in Aohan fine-wool sheep and compared the differences in microbial composition between the foregut and hindgut. In addition, we performed the functional analysis of microorganisms and analyzed the importance of intestinal microorganisms in influencing breed characteristics. The results showed the relationship between intestinal microorganisms of Aohan fine-wool sheep and their breed characteristics and environmental adaptability.

Materials and Methods

Ethics Statement

All experimental designs and operations were approved by the Animal Care and Use Committee of China Agricultural University and performed in accordance with the “Guidelines for Experimental Animals” of the Ministry of Science and Technology, Beijing, China (Permit number: SKLAB-2012- 04-07).

Animals and Sample Collection

During the research, seven Aohan fine wool rams were obtained from JinFeng Animal Husbandry Co., Ltd. (Chifeng City, Inner Mongolia Autonomous Region, China) (42°15′28.1″N, 118°53′12.7″E) (Figure 1A). All sheep were raised according to standard livestock management methods, including the same temperature and humidity (The annual average temperature is 4.9–7.4°C, and the annual precipitation is 218–595 mm), the same feed (Lambs are fed with the milk of ewes, and as adults, the feed is mainly maize straw and alfalfa hay), the same grazing time, and sufficient drinking water (Figure 1B). At the age of 12 months, seven Aohan fine wool sheep were euthanized and slaughtered. Fresh intestinal contents (∼10 g) were collected from duodenum, jejunum, ileum, cecum, colon, and rectum. Samples were placed in liquid nitrogen immediately after collection, and then transferred to –80°C ultra-low temperature freezer for storage.

FIGURE 1

Figure 1. The sampling site of this study was Aohan Banner. (A) The geographic location of raising Aohan fine-wool sheep is indicated by red coordinates. Map from the China Department of Natural Resources. Approval number: GS(2019)1652. (B) The red marks indicate the geographical areas where the Aohan fine-wool sheep live; these areas have a semi-desertified grassland environment. The scale bar is 20 km, source: Google Maps.

DNA Extraction, PCR Amplification, and Illumina MiSeq Sequencing

Microbial genome was extracted according to experimental steps using TIAN amp Genomic DNA Kit (TIANGEN Bio-Tek Co., Ltd., Beijing, China). The quality of DNA was detected by 0.8% agarose gel electrophoresis, and DNA was quantified by UV spectrophotometer. V3-V4 region of 16S rRNA gene was amplified using the following primers: ACTCCTACGGGAGGCAGCA (forward) and GGACTACHVGGGTWTCTAAT (reverse). NEB’s Q5 high-fidelity DNA polymerase (NEB Bio-Tek Inc., Ipswich, MA, United States) was used for amplification, the system is 25 μl: 5× Q5 Reaction Buffer 5 μl, 10 mM dNTPs 0.5 μl, 10 μM Forward and Reverse Primer 1.25 μl for each, Template DNA 1 μl, ddH₂O 16 μl. The thermal cycle includes initial denaturation at 98°C for 2 min, denaturation at 98°C for 15 s, annealing at 55°C for 30 s, extension at 72°C for 30 s, and final extension at 72°C for 5 min, holding at 4°C, and cycling 30 times. AXYGEN gel recovery kit (Corning Bio-Tek Inc., NY, United States) was used for gel cutting and recovery, and fluorescence quantification of PCR amplification and recovery products was performed. The fluorescence reagent is Quant-iT PicoGreen dsDNA Assay Kit (Thermo Fisher Scientific Inc., West Palm Beach, FL, United States), and the quantitative instrument is FLx800 Microplate reader (Agilent BioTek Inc., Santa Clara, CA, United States). After fluorescence quantification, the samples were mixed proportionally according to the sequencing requirements of each sample. The sequencing library was prepared using TruSeq Nano DNA LT Library Prep Kit (Illumina BioTek Inc., San Diego, CA, United States). The constructed libraries, inspected for quality using Agilent High Sensitivity DNA Kit (Agilent BioTek Inc.), presented only a single peak and no linker. MiSeq sequencer (Illumina BioTek Inc.) was used to carry out 2 × 300 bp paired-end sequencing, and the corresponding reagent is MiSeq Reagent Kit V3 (600 cycles) (Illumina BioTek Inc.). In order to ensure quality, the insert range for sequencing was 200–450 bp.

Bioinformatics Analysis and Statistics

The raw data goes through the steps of remove primer, quality filtering, denoise, splicing and de-chimerism to obtain clean reads. This process is realized through the QIIME2 DADA2 platform (Callahan et al., 2016; Bolyen et al., 2018, p. 2). Each deduplicated sequence generated using DADA2 quality control is called an operational taxonomic unit (OTU) representative sequence. By performing statistics on the leveled OTU table, the specific composition table of the microbial communities in each sample at each classification level can be obtained. Through this table, calculate the number of classification units contained in different samples at each classification level. In order to analyze the indicators of microbial richness, diversity and evenness, we conducted an Alpha diversity analysis. Next, we used bray_curtis distance algorithm to performed principal coordinate analysis (PCoA) and cluster dendrogram, and the unweighted pair-group method was used to calculate arithmetic mean. Linear discriminant analysis effect size (LEfSe) analysis realizes simultaneous differential analysis of all classification levels of microorganisms (Segata et al., 2011). Meanwhile, it searches for robust different species between groups.

PICRUSt2 is a software that predicts the functional abundance of samples based on the abundance of marker gene sequences in the samples (Douglas et al., 2020). We used the MetaCyc functional database¹ to predict 16S rRNA gene sequences. And based on the data results, we obtained different metabolic pathways between different subgroups.

GraphPad Prism (version 8.0) was used for statistical analysis. The Repeated Measures ANOVA test method was used to detect the differences of microbial abundance between intestinal segments. F > 1 indicated that difference of the mean squares between groups and within the groups is statistically significant. The criterion of significance was conducted at P < 0.05.

Results

Sequencing Results and Bacterial Diversity of Different Intestinal Segments

In total, we obtained 2,455,588 raw reads from 6 intestinal segments of seven sheep. After filtering out, denoising, chimera checking, and singleton checking, we obtained clean reads of 1,302,080 sequences. And 17,379 operational taxonomic units (OTUs) were obtained by clustering at 97% identity (Supplementary Materials 1, 2). To study the commonality and peculiarity of the microbes, we used a Venn diagram to calculate the number of OTUs shared by different intestinal segments between the groups. The frequency of the screening samples was 50%. First, we analyzed the intestinal segments of the hindgut. The number of core OTUs shared by the four hindgut segments was 380 (Figure 2A). The number of OTUs observed in the foregut was 124, which was less than that observed in the hindgut (Figure 2B). The results showed that the microbial abundance of the foregut was lower than that of the hindgut.

FIGURE 2

Figure 2. Sequencing results and statistical analysis of diversity. (A) Venn diagram showing the OTUs shared among the hindgut segments. (B) Venn diagram showing the OTUs shared between the foregut segments. (C) The Chao1 and Shannon indices of six intestinal segments. Significantly different indices were tested by Kruskal-Wallis test with adjusted *P value of < 0.05, **P > 0.01. (D) Rarefaction curve of Good’s coverage index. Each curve represents the mean within the group.

To prove the accuracy of the analysis results, we conducted a diversity analysis of intestinal microbes and used Chao1 and Shannon indexes to analysis. The results showed that the diversity and richness of intestinal microbes were higher in the hindgut (P < 0.01, F > 1.00), and there was no significant difference in microflora among the hindgut segments (ileum, cecum, colon and rectum; P > 0.05, F = 0.78; Figure 2C). The Good’s coverage index of each sample was more than 90%, and the curve tended to be flat. Among them, the index of the foregut was steep, indicating that the abundance of the detected OTUs was higher in the foregut. Meanwhile, the results show that the sequencing data were sufficient to cover all the bacterial communities (Figure 2D).

Cluster Analysis and Microbial Composition Among Different Intestinal Segments

In previous studies on the cluster analysis of intestinal microbes in sheep, the correlation between the ileum and other intestinal segments was not accurately defined (Zeng et al., 2015; Wang et al., 2017). Therefore, we assumed that the bacterial compositions of the ileum and large intestine were similar. To confirm this hypothesis, we analyzed the microbial community composition of the six intestinal segments. And found that the microbial composition of the foregut and the hindgut was very different, and that of the intestinal segment ileum and the large intestine was similar. In addition, the degree of correlation in the microflora of the duodenum and jejunum in the foregut was not too high, with 29.0 and 7.7% variations explained by principal component 1 (PC1) and PC2, respectively (Figure 3A). Next, we conducted ANOSIM similarity analysis, and the results showed that the intestinal microbial structures of the foregut and hindgut were significantly different (P < 0.01, R = 0.54) (Supplementary Material 3).

FIGURE 3

Figure 3. Cluster analysis of Aohan fine-wool Sheep. (A) Principal coordinate analysis (PCoA) based on all samples. (B) The hierarchical tree shows the UPGMA clustering result. The abscissa indicates the distance between samples, the number after the group abbreviation represents the individual number, and the branch length indicates similarity. On the right is the stacked histogram of the top 10 abundant bacterial families in the sheep intestine. The abscissa indicates the proportion of bacteria. (C) The phylum-level microbial composition of each intestinal segment. (D) The genus-level microbial composition of each intestinal segment.

In the hindgut, for the similar microbial composition clusters, the similarity between individuals in the same intestinal segment is higher than that between different intestinal segment within individuals (Figure 3B). The microorganisms in the hindgut mainly included Ruminococcaceae, Christensenellaceae, Lachnospiraceae, Rikenellaceae, and Akkermansiaceae. These five types of bacteria accounted for 65% of all the gut microbes (Figure 3B). The similarity in microbial composition may be related to functional consistency. In the foregut, the duodenum and jejunum were not clustered among individuals, and the jejunum was quite different among individuals. Burkholderiaceae, Ruminococcaceae, and Saccharimonadaceae were the dominant bacterial families. To judge whether the degree of dispersion of samples within a group is different among different groups, we used the ADONIS permutation test. The degree of dispersion of the foregut was significantly higher than that of the hindgut, and there was no significant difference in the foregut group (Supplementary Material 4).

In the analysis of microbial composition, we focused on exploring the top 10 bacteria at the phylum and genus levels. We analyzed the composition of microorganisms at each taxonomic level. At the phylum level, Firmicutes (average 52.88%) and Bacteroidetes (average 30.23%) were the dominant bacteria in the hindgut. In the foregut, Proteobacteria (average 31.70%) and Firmicutes (average 47.44%) were the dominant bacterial phyla. The proportion of Bacteroidetes (hindgut 30.16%, foregut 1.95%), Verrucomicrobia (hindgut 8.02%, foregut 0.19%), and Spirochaetes (hindgut 2.08%, foregut 0.01%) in the foregut and hindgut was significantly different (P < 0.01, F > 1.00) (Figure 3C). This shows that these bacteria are more active in the hindgut, indicating their role in the fermentation and water absorption processes of the hindgut. At the genus level, Ruminococcaceae_UCG-005, Christensenellaceae_R-7_group, Akkermansia (sum 26.14%) were the dominant genera in the hindgut. The proportion of Christensenellaceae_R-7_group (hindgut 7.43%, foregut 7.06%) remained the same throughout the intestine (P > 0.50, F > 1.00) (Figure 3D). Compared to the foregut, most of the annotated bacteria were more abundant in the hindgut. Additionally, Candidatus-Saccharimonas was a unique bacterium in the foregut.

Microbial Communities of Different Intestinal Segments

To further explore the differences between the samples, we conducted a Kaplan-Meier analysis. According to the results of the cluster analysis, we performed the Wilcoxon test and LEfSe analysis on the foregut and hindgut, respectively. (Supplementary Material 5). The genus of more abundant bacteria in the cecum was Ruminococcaceae_UCG_010, and Caldicoprobacter in the ileum was a significantly abundant bacteria belonging to the Firmicutes phylum (LDA > 2, P < 0.05; Figure 4A). Coriobacteriales in the colon belongs to the phylum Actinobacteria. The rectum included significantly abundant bacteria Family_XIII_UCG_001 under Clostridium, Eggerthellaceae under Actinomycetes, and p_251_o5 under Bacteroides (LDA > 2, P < 0.05; Figure 4B). In the analysis of the foregut, we found that the significantly abundant bacteria of the duodenum mainly included Succiniclasticum under Acidaminococcaceae. The abundance of Ruminococcaceae, and Defluviitaleaceae under the class of Clostridia was also significantly higher than that in the jejunum (LDA > 2, P < 0.05; Figure 4C). Escherichia_Shigella, Akkermansiaceae, Veillonellaceae, and Butyricicoccus were more abundant in the jejunum compared to duodenum, which may be related to the specific digestive function of the jejunum (LDA > 2, P < 0.05; Figure 4D).

FIGURE 4

Figure 4. Linear discriminant analysis (LDA) effect size (LEfSe) analysis of the intestinal segments of Aohan fine-wool sheep. The LEfSe analysis histogram of hindgut (A,C) foregut. The ordinate is the taxa with significant differences between groups, and the abscissa is a bar graph displaying the LDA logarithmic score value of each taxon. The longer the length, the more significant the difference of the taxon, and the color of the bar graph indicates the sample group with the highest abundance corresponding to the taxon. The LEfSe analysis branch diagram of hindgut (B,D) foregut. The node size corresponds to the average relative abundance of the taxa, and the hollow nodes represent taxa with insignificant differences between groups. The letters identify the names of taxa that differ significantly between the groups.

Microbial Function Prediction and Intestinal Metabolic Pathways

To investigate the functional differences in the intestinal microbiota of Aohan fine-wool sheep, we performed a functional analysis of microbiota using PICRUST2 (Supplementary Material 6). First, we counted the abundance of metabolic pathways in six intestinal segments using the MetaCyc database as a reference. The abundance statistics of functions revealed that more microbial functions were related to biosynthesis and metabolism (Figure 5A). The abundance of the biosynthesis processes was significantly higher than that of the other metabolic pathways. The main biological pathways in the intestine included amino acid biosynthesis, nucleoside and nucleotide biosynthesis, biosynthesis of cofactors, prosthetic groups, electron carriers, and vitamins, and the biological processes of biosynthesis, and fatty acid and lipid biosynthesis.

FIGURE 5

Figure 5. Predictive analysis and statistics of foregut and hindgut microbial function. (A) The abundance of the differential metabolic pathways based on the MetaCyc database. The abscissa is the abundance count of the classification, the ordinate is the functional pathway of MetaCyc’s second classification level, and the rightmost is the first-level classification to which this pathway belongs. (B) The differential analysis of metabolic pathways based on the metagenomeSeq method. Light blue represents the foregut and light yellow represents the hindgut. The right ordinate is the corrected q-value and the left ordinates are different pathway labels.

To investigate the differences between the intestinal metabolic pathways, we used q-value < 0.05 as the standard for the differential enrichment analysis. The results showed that the biological pathways of catechol degradation, salicylic acid degradation, aromatic compound degradation, and the citric acid cycle pathways are more active in the foregut. In contrast, the hindgut mainly includes the biological processes of chondroitin sulfate degradation and protein N-glycosylation synthesis (Figure 5B).

Discussion

This study conducted a comprehensive analysis of the intestinal microbes of Aohan fine-wool sheep. The 16S rRNA gene sequencing method was used to analyze the microbial structure, composition, significantly abundant bacteria, and potential functions of each intestinal segment.

The small intestine is an important organ for digestion and absorption. The large intestine mainly absorbs water and absorbable nutrients. It can also produce acetate, propionate, and butyrate as sources of nutrients for intestinal cells (Miller and Wolin, 1996). The results of the PCoA analysis revealed that the foregut was significantly separated of hindgut. Our results are consistent with the analysis results in Mongolian sheep (Zeng et al., 2017). Our study showed that the intestinal microbial diversity and richness of the foregut of Aohan fine-wool sheep were significantly lower than those of the hindgut. Previous studies have confirmed that the microbial diversity in the small intestine is low (Eckburg et al., 2005). The microbial population in the small intestine is usually affected by a host of adverse factors, such as low pH, faster transit time, and exposure to bile acids and antimicrobial peptides. Studies have shown that Firmicutes and Proteobacteria, which dominate the small intestine, are more tolerant of these factors (Angelakis et al., 2015; Leone et al., 2015). This explains the abundance of Proteobacteria in the foregut of Aohan fine-wool sheep. Proteobacteria are also the dominant phylum on other species, including rumen of calves (Rey et al., 2014) and Yimeng black goats (Li et al., 2021). Interestingly, in previous studies on Mongolian sheep (raised in the Gansu province of China), Proteobacteria was not found to be the dominant phylum (Zeng et al., 2017). However, in camels (bred in Xilin Gol, Inner Mongolia Autonomous Region, China), the results were similar to those from our analysis (He et al., 2018). This indicates that the intestinal microflora is influenced by differences in rearing environments. The Bacteroidetes phylum is the dominant phylum in the large intestine, and the same results have been found in other mammals (Donaldson et al., 2016). The large intestine is characterized by a slow flow rate and a neutral to slightly acidic pH. It is beneficial for the colonization of Bacteroidetes (Flint et al., 2012). In general, our results indicate that in Aohan fine-wool sheep, the intestinal environment behind the ileum is more uniformly characterized by mild pH and low-speed transport. At the phylum level, the dominant phylum of gut microbes in ruminants and monogastric animals were consistent.

The results of our study indicated that the same microbial population exists in the foregut and hindgut and performs functions in the gut. Firmicutes and Bacteroidetes were the dominant flora in the hindgut. It participates in energy metabolism and affects obesity (Ley et al., 2006b; Komaroff, 2017). In addition, we found that the abundance of Christensenellaceae and Lachnospiraceae in the foregut and hindgut were relatively stable (P > 0.05), indicating that they play an important role in the intestinal tract. The Christensenellaceae_R-7_group is a member of the Christensenellaceae family. The Christensenellaceae family is a relatively new bacterial family that has previously been related to the host’s health (Waters and Ley, 2019). Moreover, Christensenellaceae is positively correlated with protein catabolism and intestinal metabolites of dietary animal proteins (Beaumont et al., 2017; Manor et al., 2018). In a study on ruminants, Christensenellaceae_R-7_group improved the development of the rumen and increased the absorption and digestion of nutrients (Couch et al., 2020; Huang et al., 2021). These results indicate that Christensenellaceae may be an important part of the gastrointestinal tract of ruminants. Lachnospiraceae is the main component of the intestinal microbiota of ruminants (Kittelmann et al., 2013), and is closely related to the production of butyrate (Haas and Blanchard, 2017). Our results are consistent with those of previous studies (Freetly et al., 2020; Li et al., 2020) and indicate that these stable bacteria are involved in the growth of ruminants.

In this study, we found significantly abundant bacteria in each intestinal segment. Interestingly, previous reports on the gut bacterial of Mongolian sheep were not focused on these bacteria, possibly due to differences in feeding environment and breed adaptation (Zeng et al., 2015, 2017). Butyricicoccus is a type of bacteria that produces butyric acid, providing butyrate as the main nutrient for intestinal epithelial cells (Geirnaert et al., 2014). Akkermansiaceae is related to gastrointestinal homeostasis and metabolic balance (Clarke et al., 2014). These foregut bacteria may be responsible for imparting the characteristics of strong adaptability and tolerance to rough feeding in Aohan fine-wool sheep. The significantly abundant bacteria in the hindgut mainly included Caldicoprobacteraceae and Ruminococcaceae. These bacteria can ferment a variety of nutrients and produce volatile fatty acids. In addition, indigestible polysaccharides can be used to produce metabolites that are beneficial to the intestinal tract (Deng et al., 2017; Guo et al., 2020a). A study showed that Ruminococcaceae-UCG-010 and UCG-005 (at the genus level) are related to the degradation of starch and fiber in ruminants (Kim et al., 2011). These communities may contribute to further fermentation of feed in the cecum. However, these core bacteria were not highly abundant in the cecum of Small Tail Han sheep (Zhang et al., 2018). This difference may be related to the characteristics of the breed and the distribution area. Small Tail Han sheep are distributed in the Shandong Province of China and are known for their very high rates of reproduction and extremely high fecundity. However, the meat quality of Small Tail Han sheep is worse than that of Mongolian sheep (Cannas, 2011). Compared with Mongolian sheep, the bacteria identified in this study has not been reported before (Zeng et al., 2017). However, in other ruminants, such as sika deer, Ruminococcaceae-UCG-010 is reported to be the main bacterial species in the hindgut, where it degrades cellulose and produces short-chain fatty acids (SCFAs) (Li et al., 2019). Similar results have been reported in newborn calves (Elolimy et al., 2020). These observations are consistent with the results of our study. Coriobacteriaceae and Eggerthellaceae belong to the phylum Actinobacteria. Studies have indicated that the abundance of Eggerthellaceae is positively correlated with feed efficiency (Bach et al., 2019). Previously, research on the function of these highly abundant bacteria revealed that strains of Coriobacteriaceae and Eggerthellaceae are particularly involved in the metabolism of daidzein and genistein, and they can convert to food polyphenols (Soukup et al., 2021). However, The Coriobacteriaceae and Eggerthellaceae have not been observed and analyzed in other breeds of fine-wool sheep (Yang et al., 2021). In this study, the ability of Aohan fine-wool sheep to metabolize fiber polysaccharides and flavonoids may be related to the grazing environment. The grassland composition of semi-desert pastures is complex. In the current study, these bacteria found in the hindgut were not of interest in other breeds of sheep. Increased concentrations of these bacteria, such as Caldicoprobacteraceae and Ruminococcaceae in the hindgut, enhance the digestibility of crude fiber while producing beneficial metabolites.

The microbial potential function analysis revealed that the foregut plays an important role in the metabolism of biomass and produces important biosynthetic precursors. The methyl citrate cycling pathway in the foregut can metabolize propionate to pyruvate (Brock et al., 2002). Meanwhile, the metabolism of catechol and the degradation of salicylic acid is more active in the foregut. This pathway produces the energy substrates succinate and acetyl CoA, which are involved in energy metabolism (MacLean et al., 2006; Cámara et al., 2007). Succiniclasticum seems to be involved in the metabolite transformation process, which provides energy to ruminants during the conversion of succinic acid to propionic acid (Van Gylswyk, 1995). Therefore, the combined action of these substances and bacteria provides ruminants with additional energy. Our analysis showed that the metabolism of the hindgut involves the degradation of polysaccharides and the biosynthesis of carbohydrates (Olson, 2002). However, PICRUSt2 still has drawbacks in predicting the potential functions. Since 16S rRNA gene amplification is based on sequencing and analysis of hypervariable regions, functional prediction cannot provide resolution to distinguish strain-specific functionality, this leads to an important limitation of PICRUSt2 and any amplicon-based analysis (Douglas et al., 2020). These results demonstrated the bacteria involvement in the process of metabolic conversion and predicted the potential relationship with the roughing tolerance of Aohan Fine Wool Sheep.

Conclusion

This study revealed the unique gut microbial compositions of Aohan fine-wool sheep growing in a semi-desert environment. The structure and composition of gut microbes in Aohan fine-wool sheep are associated with stronger environmental adaptation and gut health. Our results revealed a robust energy metabolism in the foregut. The hindgut has a strong ability to digest crude fiber and produce SCFAs. Therefore, these findings provide a baseline for understanding the complex intestinal microbiota adapt to the living environment and provides new insights into improving stress resistance and gut health in sheep through microbes.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Animal Care and Use Committee of China Agricultural University. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

XuD: experiment conception and ranch work arrangement. YM: ranch sampling, complete analysis, write article. XiD: complete the preliminary analysis of the article, ranch sampling. JW, XY, and TL: mainly involve in the sampling of pastures. DH participate in pasture sampling and provide assistance in experimental sample processing. LD, RL, and WR: provide laboratory animals and help with sampling. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Programs of National Key R&D Program of China, (grant no. 2021YFF1000700), the Major Project for Cultivation Technology of New Varieties of Genetically Modified Organisms of the Ministry of Agriculture (grant no. 2016ZX08008-001), the National Natural Science Foundation of China (grant no. 32002145), and the Ningxia Agricultural Breeding Project (grant no. NXNYYZ20150103).

Conflict of Interest

RL was employed by Inner Mongolia Grassland Jinfeng Animal Husbandry Co., Ltd.

The remaining 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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.874536/full#supplementary-material

Footnotes

1. ^ https://metacyc.org/

References

Al-Masaudi, S., El Kaoutari, A., Drula, E., Al-Mehdar, H., Redwan, E. M., Lombard, V., et al. (2017). A metagenomics investigation of carbohydrate-active enzymes along the gastrointestinal tract of Saudi sheep. Front. Microbiol. 8:666. doi: 10.3389/fmicb.2017.00666

PubMed Abstract | CrossRef Full Text | Google Scholar

Angelakis, E., Armougom, F., Carrière, F., Bachar, D., Laugier, R., Lagier, J.-C., et al. (2015). A metagenomic investigation of the duodenal microbiota reveals links with obesity. PLoS One 10:e0137784. doi: 10.1371/journal.pone.0137784

PubMed Abstract | CrossRef Full Text | Google Scholar

Bach, A., López-García, A., González-Recio, O., Elcoso, G., Fàbregas, F., Chaucheyras-Durand, F., et al. (2019). Changes in the rumen and colon microbiota and effects of live yeast dietary supplementation during the transition from the dry period to lactation of dairy cows. J. Dairy Sci. 102, 6180–6198. doi: 10.3168/jds.2018-16105

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaumont, M., Portune, K. J., Steuer, N., Lan, A., Cerrudo, V., Audebert, M., et al. (2017). Quantity and source of dietary protein influence metabolite production by gut microbiota and rectal mucosa gene expression: a randomized, parallel, double-blind trial in overweight humans. Am. J. Clin. Nut. 106, 1005–1019. doi: 10.3945/ajcn.117.158816

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C., Al-Ghalith, G. A., et al. (2018). QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. PeerJ Preprints 6:e27295v2. doi: 10.7287/peerj.preprints.27295v2

CrossRef Full Text | Google Scholar

Brock, M., Maerker, C., Schütz, A., Völker, U., and Buckel, W. (2002). Oxidation of propionate to pyruvate in Escherichia coli. Involvement of methylcitrate dehydratase and aconitase. Eur. J. Biochem. 269, 6184–6194. doi: 10.1046/j.1432-1033.2002.03336.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Cámara, B., Bielecki, P., Kaminski, F., dos Santos, V. M., Plumeier, I., Nikodem, P., et al. (2007). A gene cluster involved in degradation of substituted salicylates via ortho cleavage in Pseudomonas sp. strain MT1 encodes enzymes specifically adapted for transformation of 4-methylcatechol and 3-methylmuconate. J. Bacteriol. 189, 1664–1674. doi: 10.1128/JB.01192-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Cannas, A. (2011). The International Sheep and Wool Handbook. Ital. J. Anim. Sci. 10:e45.

Google Scholar

Cho, I., and Blaser, M. J. (2012). The human microbiome: at the interface of health and disease. Nat. Rev. Genet. 13, 260–270. doi: 10.1038/nrg3182

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, S. F., Murphy, E. F., O’Sullivan, O., Lucey, A. J., Humphreys, M., Hogan, A., et al. (2014). Exercise and associated dietary extremes impact on gut microbial diversity. Gut 63, 1913–1920. doi: 10.1136/gutjnl-2013-306541

PubMed Abstract | CrossRef Full Text | Google Scholar

Couch, C. E., Arnold, H. K., Crowhurst, R. S., Jolles, A. E., Sharpton, T. J., Witczak, M. F., et al. (2020). Bighorn sheep gut microbiomes associate with genetic and spatial structure across a metapopulation. Sci. Rep. 10:6582. doi: 10.1038/s41598-020-63401-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, K., Yang, Z., Darwish, H., Zhang, Y., Ge, Y., Zhang, X., et al. (2014). Molecular cloning and characterization of the β-catenin gene from fine-wool sheep. Gene 546, 277–282. doi: 10.1016/j.gene.2014.05.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Y., Huang, Z., Zhao, M., Ruan, W., Miao, H., and Ren, H. (2017). Effects of co-inoculating rice straw with ruminal microbiota and anaerobic sludge: digestion performance and spatial distribution of microbial communities. Appl. Microbiol. Biotechnol. 101, 5937–5948. doi: 10.1007/s00253-017-8332-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Donaldson, G. P., Lee, S. M., and Mazmanian, S. K. (2016). Gut biogeography of the bacterial microbiota. Nat. Rev. Microbiol. 14, 20–32. doi: 10.1038/nrmicro3552

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688. doi: 10.1038/s41587-020-0548-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635–1638. doi: 10.1126/science.1110591

PubMed Abstract | CrossRef Full Text | Google Scholar

Elolimy, A., Alharthi, A., Zeineldin, M., Parys, C., and Loor, J. J. (2020). Residual feed intake divergence during the preweaning period is associated with unique hindgut microbiome and metabolome profiles in neonatal Holstein heifer calves. J. Anim. Sci. Biotechnol. 11:13. doi: 10.1186/s40104-019-0406-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Flint, H. J., Scott, K. P., Louis, P., and Duncan, S. H. (2012). The role of the gut microbiota in nutrition and health. Nat. Rev. Gastroenterol. Hepatol. 9, 577–589.

Google Scholar

Freetly, H. C., Dickey, A., Lindholm-Perry, A. K., Thallman, R. M., Keele, J. W., Foote, A. P., et al. (2020). Digestive tract microbiota of beef cattle that differed in feed efficiency. J. Anim. Sci. 98:skaa008. doi: 10.1093/jas/skaa008

PubMed Abstract | CrossRef Full Text | Google Scholar

Geirnaert, A., Steyaert, A., Eeckhaut, V., Debruyne, B., Arends, J. B., Van Immerseel, F., et al. (2014). Butyricicoccus pullicaecorum, a butyrate producer with probiotic potential, is intrinsically tolerant to stomach and small intestine conditions. Anaerobe 30, 70–74. doi: 10.1016/j.anaerobe.2014.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, P., Zhang, K., Ma, X., and He, P. (2020a). Clostridium species as probiotics: potentials and challenges. J. Anim. Sci. Biotechnol. 11:24. doi: 10.1186/s40104-019-0402-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Lei, H., Zhang, K., Ke, F., and Song, C. (2020b). Diversification of animal gut microbes and NRPS gene clusters in some carnivores, herbivores and omnivores. Biotechnol. Biotechnol. Equip. 34, 1280–1287. doi: 10.1080/13102818.2020.1835536

CrossRef Full Text | Google Scholar

Haas, K., and Blanchard, J. L. (2017). Kineothrix alysoides, gen. Nov., sp. nov., a saccharolytic butyrate-producer within the family Lachnospiraceae. Int. J. Syst. Bacteriol. 67, 402–410. doi: 10.1099/ijsem.0.001643

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J., Yi, L., Hai, L., Ming, L., Gao, W., and Ji, R. (2018). Characterizing the bacterial microbiota in different gastrointestinal tract segments of the Bactrian camel. Sci. Rep. 8:654. doi: 10.1038/s41598-017-18298-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Ge, F., Yao, X., Guo, X., Bao, P., Ma, X., et al. (2021). Microbiome and Metabolomics Reveal the Effects of Different Feeding Systems on the Growth and Ruminal Development of Yaks. Front. Microbiol. 12:682989. doi: 10.3389/fmicb.2021.682989

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Morrison, M., and Yu, Z. (2011). Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol. Ecol. 76, 49–63. doi: 10.1111/j.1574-6941.2010.01029.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kittelmann, S., Seedorf, H., Walters, W. A., Clemente, J. C., Knight, R., Gordon, J. I., et al. (2013). Simultaneous Amplicon Sequencing to Explore Co-Occurrence Patterns of Bacterial, Archaeal and Eukaryotic Microorganisms in Rumen Microbial Communities. PLoS One 8:e47879. doi: 10.1371/journal.pone.0047879

PubMed Abstract | CrossRef Full Text | Google Scholar

Komaroff, A. L. (2017). The microbiome and risk for obesity and diabetes. Jama 317, 355–356. doi: 10.1001/jama.2016.20099

PubMed Abstract | CrossRef Full Text | Google Scholar

Leone, V., Gibbons, S. M., Martinez, K., Hutchison, A. L., Huang, E. Y., Cham, C. M., et al. (2015). Effects of diurnal variation of gut microbes and high-fat feeding on host circadian clock function and metabolism. Cell Host Microbe 17, 681–689. doi: 10.1016/j.chom.2015.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ley, R. E., Peterson, D. A., and Gordon, J. I. (2006a). Ecological and Evolutionary Forces Shaping Microbial Diversity in the Human Intestine. Cell 124, 837–848. doi: 10.1016/j.cell.2006.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Ley, R. E., Turnbaugh, P. J., Klein, S., and Gordon, J. I. (2006b). Human gut microbes associated with obesity. Nature 444, 1022–1023. doi: 10.1038/4441022a

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, A., Yang, Y., Zhang, Y., Lv, S., Jin, T., Li, K., et al. (2021). Microbiome analysis reveals the alterations in gut microbiota in different intestinal segments of Yimeng black goats. Microb. Pathog. 155:104900. doi: 10.1016/j.micpath.2021.104900

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Li, R., Chen, H., Gao, J., Wang, Y., Zhang, Y., et al. (2020). Effect of different seasons (spring vs summer) on the microbiota diversity in the feces of dairy cows. Int. J. Biometeorol. 64, 345–354. doi: 10.1007/s00484-019-01812-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Si, H., Nan, W., Wang, X., Zhang, T., and Li, G. (2019). Bacterial community and metabolome shifts in the cecum and colon of captive sika deer (Cervus nippon) from birth to post weaning. FEMS Microbiol. Lett. 366:fnz010. doi: 10.1093/femsle/fnz010

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, L., Trabi, E. B., Xie, F., and Mao, S. (2021). Comparison of the fermentation and bacterial community in the colon of Hu sheep fed a low-grain, non-pelleted, or pelleted high-grain diet. Appl. Microbiol. Biotechnol. 105, 2071–2080. doi: 10.1007/s00253-021-11158-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Zhao, J., Wang, W., Guo, P., Lu, W., Wang, C., et al. (2018). Dietary corn bran altered the diversity of microbial communities and cytokine production in weaned pigs. Front. Microbiol. 9, 2090. doi: 10.3389/fmicb.2018.02090

PubMed Abstract | CrossRef Full Text | Google Scholar

MacLean, A. M., MacPherson, G., Aneja, P., and Finan, T. M. (2006). Characterization of the β-ketoadipate pathway in Sinorhizobium meliloti. Appl. Environ. Microbiol. 72, 5403–5413. doi: 10.1128/aem.00580-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Manor, O., Zubair, N., Conomos, M. P., Xu, X., Rohwer, J. E., Krafft, C. E., et al. (2018). A multi-omic association study of trimethylamine N-oxide. Cell Rep. 24, 935–946. doi: 10.1016/j.celrep.2018.06.096

PubMed Abstract | CrossRef Full Text | Google Scholar

Masters, D. G., Purser, D. B., Yu, S. X., Wang, Z. S., Yang, R. Z., Liu, N., et al. (1990). Production from fine wool sheep in three areas in northern China. Asian Australasian J. Anim. Sci. 3, 305–312. doi: 10.5713/ajas.1990.305

CrossRef Full Text | Google Scholar

Miller, T. L., and Wolin, M. J. (1996). Pathways of acetate, propionate, and butyrate formation by the human fecal microbial flora. Appl. Environ. Microbiol. 62, 1589–1592. doi: 10.1128/aem.62.5.1589-1592.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Neish, A. S. (2009). Microbes in Gastrointestinal Health and Disease. Gastroenterology 136, 65–80. doi: 10.1053/j.gastro.2008.10.080

PubMed Abstract | CrossRef Full Text | Google Scholar

Olson, L. (2002). Clearing a path for nerve growth. Nature 416, 589–590. doi: 10.1038/416589a

PubMed Abstract | CrossRef Full Text | Google Scholar

Rey, M., Enjalbert, F., Combes, S., Cauquil, L., Bouchez, O., and Monteils, V. (2014). Establishment of ruminal bacterial community in dairy calves from birth to weaning is sequential. J. Appl. Microbiol. 116, 245–257. doi: 10.1111/jam.12405

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome. Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text | Google Scholar

Shreiner, A. B., Kao, J. Y., and Young, V. B. (2015). The gut microbiome in health and in disease. Curr. Opin. Gastroenterol. 31, 69–75. doi: 10.1097/MOG.0000000000000139

PubMed Abstract | CrossRef Full Text | Google Scholar

Soukup, S. T., Stoll, D. A., Danylec, N., Schoepf, A., Kulling, S. E., and Huch, M. (2021). Metabolism of Daidzein and Genistein by Gut Bacteria of the Class Coriobacteriia. Foods 10:2741. doi: 10.3390/foods10112741

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbaugh, P. J., Ley, R. E., Hamady, M., Fraser-Liggett, C. M., Knight, R., and Gordon, J. I. (2007). The Human Microbiome Project. Nature 449, 804–810. doi: 10.1038/nature06244

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Gylswyk, N. O. (1995). Succiniclasticum ruminis gen. nov., sp. nov., a ruminal bacterium converting succinate to propionate as the sole energy-yielding mechanism. Int. J. Syst. Evol. Microbiol. 45, 297–300. doi: 10.1099/00207713-45-2-297

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Fan, H., Han, Y., Zhao, J., and Zhou, Z. (2017). Characterization of the microbial communities along the gastrointestinal tract of sheep by 454 pyrosequencing analysis. Asian Australasian J. Anim. Sci. 30, 100–110. doi: 10.5713/ajas.16.0166

PubMed Abstract | CrossRef Full Text | Google Scholar

Waters, J. L., and Ley, R. E. (2019). The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 17:83. doi: 10.1186/s12915-019-0699-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Y., Xiang, Y., Zhou, W., Chen, J., Li, K., and Yang, H. (2017). Microbial community mapping in intestinal tract of broiler chicken. Poult. Sci. 96, 1387–1393. doi: 10.3382/ps/pew372

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao-ping, Z. (2009). Determination of Breeding Objective Traits of Aohan Merino Sheep and Calculation of Their Marginal Benefit. J. Anhui Agric. Sci. 2009:13.

Google Scholar

Yang, J., Yu, Z., Wang, B., and Ndayisenga, F. (2021). Gut region induces gastrointestinal microbiota community shift in Ujimqin sheep (Ovis aries): from a multi-domain perspective. Environ. Microbiol. 23, 7603–7616. doi: 10.1111/1462-2920.15782

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Y., Zeng, D., Ni, X., Zhu, H., Jian, P., Zhou, Y., et al. (2017). Microbial community compositions in the gastrointestinal tract of Chinese Mongolian sheep using Illumina MiSeq sequencing revealed high microbial diversity. Amb. Exp. 7:75. doi: 10.1186/s13568-017-0378-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Y., Zeng, D., Zhang, Y., Ni, X., Tang, Y., Zhu, H., et al. (2015). Characterization of the cellulolytic bacteria communities along the gastrointestinal tract of Chinese Mongolian sheep by using PCR-DGGE and real-time PCR analysis. World J. Microbiol. Biotechnol. 31, 1103–1113. doi: 10.1007/s11274-015-1860-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Shao, M., Huang, H., Wang, S., Ma, L., Wang, H., et al. (2018). The Dynamic Distribution of Small-Tail Han Sheep Microbiota across Different Intestinal Segments. Front. Microbiol. 9:32. doi: 10.3389/fmicb.2018.00032

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Qian, Q., Mao, Y., Xu, Y., Yang, Y., Chen, Y., et al. (2021). Characterization of growth phenotypes and gastrointestinal tract microbiota in sheep fed with caragana. J. Appl. Microbiol. 131, 2763–2779. doi: 10.1111/jam.15138

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, W., Wang, Y., Liu, S., Huang, J., Zhai, Z., He, C., et al. (2015). The dynamic distribution of porcine microbiota across different ages and gastrointestinal tract segments. PLoS One 10:e0117441. doi: 10.1371/journal.pone.0117441

PubMed Abstract | CrossRef Full Text | Google Scholar