Haibi Zhao
Feng Lv
Xiaojuan Wang- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
Starch and NDF are the main components in the diets of ruminants worldwide and are the main energy source for rumen microorganisms and hosts. The purpose of this study was to investigate the effects of different NDF/starch ratios on rumen fermentation parameters, rumen development and rumen microbes in lambs and to predict the function of rumen microbes by metagenomic techniques. In this study, 30 lambs with birth weights of (3.0 ± 0.5) kg were selected. The lambs of Hu sheep were randomly divided into two groups, fed starter with an NDF/starch ratio of 0.5 (group A) or 1.0 (group C). Samples of the rumen tissue and contents were collected after slaughter. The results showed that the ADG and ADFI of group A were significantly higher than those of group C (P < 0.05), but there was no significant difference in the FCR (P > 0.05). Therefore, from the perspective of feed-related economic benefits, group C showed greater economic value; the A/P of group C was significantly lower than that of group A (0.05 < P < 0.1), and the TVFA showed no significant difference (P > 0.05); The lengths of the rumen papillae of group C was significantly higher than that of group A (0.05 < P < 0.1). There was no significant difference in the abundance of the top 10 species at the phylum level and genus level (P > 0.05). CAZymes gene enrichment was observed in the rumen microbial community of lambs in group C (P < 0.05). In conclusion, group C, fed with starter with a higher NDF/starch ratio, had a higher feeding value. This study provides comprehensive insights into the composition of NDF and starch in lamb starter.
Introduction
The Hu sheep is an extremely inbred species and the most popular species in China due to its excellent prolificacy (3–4 lambs per parturition), rapid growth rate, and adaptation to poor-quality feeds. Young ruminants are born with underdeveloped rumens and are considered functional monogastric animals before weaning, but rumen development in young ruminants plays a crucial role in animal health and nutrition (1). Colonization of the gastrointestinal tract of young ruminants begins at birth and continues through successive colonization and floristic changes until the microbiota reaches a stable state later in life, and changes in early succession of primary colonizers affect the composition and function of the mature animal community (2). The rumen epithelium is a key barrier between the host and the rumen environment. It can provide energy for ruminants by using the metabolites of rumen microorganisms, such as acetic acid, propionic acid and butyric acid, which are produced by fermentation of microorganisms (3). The rumen epithelium influences the use of nutrients throughout the body by absorbing VFAs (3, 4). VFAs can provide energy for the host and stimulate the development of the rumen epithelium. Early dietary intervention is crucial to the formation of the microflora and has a great impact on improving the performance of young ruminants. The early use of starter is very important for rumen functional development and optimal growth of lambs. Determination of the appropriate nutrient level and preparation of the appropriate feed for rumen development are beneficial for rumen development and for complete rumen microbial function in lambs (5). Neutral detergent fiber (NDF), as the main component of structural carbohydrates in forage, is a vital factor impacting the performance of young ruminants (6). Studies have shown that rumen fermentation levels and the promotion of rumen development can be positively regulated by adjusting dietary NDF in the early life stages of calves and lambs, including adjustments in the levels and sources of NDF (7). The main components of starch are non-structural carbohydrates, which can be rapidly degraded in the rumen to provide energy for lambs and rumen microbes and produce volatile fatty acids (VFAs) to stimulate the development of rumen papilla (7). In ruminants, starch is fermented to VFAs in the rumen, which may provide energy for epithelia and the host animal (8). Early use of starter is very important for the development of rumen function and the optimal growth of lambs. Determining the appropriate feed level and modulating the appropriate carbohydrate level for rumen development are beneficial to the development of lambs and promotion of rumen microbial function (5). Diet is one of the main factors leading to changes in the rumen microbial community (9). Considerable research has focused on cataloging the adult rumen microbiome and its relationship to complex diets (9, 10), and studies on the lamb rumen microbiota have been restricted to 16 S rRNA amplicon sequencing and/or different diets and hosts (11, 12). Although emerging studies have been able to use metagenomics to identify the rumen microbiome at the structural and functional levels, few studies have used metagenomics to study the effects of starters on rumen development in lambs, especially the effects of feeding with diets containing different NDF/starch ratios on the rumen microbiome of lambs. Metagenomics is a powerful new tool for understanding the composition and function of complex microbial communities, and it allows the abundance of all genes present in the microbial community to be determined and metabolic pathways to be predicted (13). In this paper, we applied metagenomics to analyze the effects of starters with different NDF/starch ratios on the growth, development and rumen microbes of lambs.
Materials and methods
Animals and feeding management
Experimental animals were provided by Lanzhou Tianxin Breeding Co., Ltd., and the experimental animals were raised by Defu Agricultural Technology Co., Ltd., Minqin County, Gansu Province. Feeding and management were conducted as follows: we adopted a single-cage feeding mode and immunized all lambs in strict accordance with the normal management procedures of the feeding company. Before experiments were conducted, the farm and all available instruments were sterilized with sanitizine and alcohol.
Experimental design
Thirty male lambs of Hu sheep, with a birth weight of 3.0 ± 0.5 kg, were randomly divided into two groups, with 15 lambs in each group. The NDF/starch ratio of the starter was 0.5 (group A) or 1.0 (group C). A total of 100 g of starter was weighed out in the morning of the 1st day, and then the appropriate feeding amount for the day was weighed out according to the amount remaining from the previous day every morning (ensuring that the lamb received enough starter). Starting at 10 days of age, each group was fed the same milk replacer (Beijing Precision Animal Nutrition Research Center; nutrition levels: 96.91% dry matter (DM), 23.22% crude protein (CP), and 13.20% ether extract; patent number: ZL02128844.5) and starters with different NDF/starch ratios every day. The milk replacer was fed until 35 days of age, and the starter was fed to the lambs until 56 days of age. During the feeding period, all lambs were free access to drinking water, and the amount of milk replacer was 1.5% of the weight of the lambs as a reference. The starter used in this experiment was produced by Gansu Runmu Biological Engineering Co., Ltd., and was formulated according to the Feeding Standard for Sheep and Goats for Meat (NY/T816-2004). The composition and nutritional composition of the starter are shown in Table 1.
Table 1. Composition and nutrient levels of starter (air-dry basis) %.
Sample collection
During the feeding process, the feed amount remaining from the previous day and the feed amount for the day were measured every morning. Lambs were weighed before the start of the experiment, every 7 days thereafter, and again at the end of the experiment before slaughter. After lambs were fed to 56 days of age, six lambs with similar body weights in each group were slaughtered, and the rumen pH was measured by a pH meter immediately after slaughter. Rumen ventral sac tissue samples that were 1 cm2 in size were collected after rinsing with normal saline and fixed in 10–15 volumes of 10% neutral formaldehyde solution for tissue sectioning. Moreover, the rumen contents and rumen ventral sac tissues were collected, placed in cryo-storage tubes and immediately placed in liquid nitrogen for temporary storage. After returning to the laboratory, the samples were transferred to a −80°C freezer for storage and subsequent experimental operations. Since there was no significant difference in the macro phenotype production data between the two groups, the rumen contents of the two groups were extracted by metagenomic sequencing for microscopic analysis.
Sample processing
Determination of NH3-N and VFA contents and MCP
The NH3-N concentration was determined by the phenol-sodium hypochlorite colorimetric method according to Wang et al. (14). The VFA content in rumen fluid was determined by gas chromatography (GC-2010 Plus; Shimadzu, Kyoto, Japan). An internal standard method was used, and the internal standard-based was 2-ethylbutyric acid (2EB). The chromatographic column was an AT-FFAP capillary column (30 m × 0.32 mm × 0.25 μm). The temperature program of the chromatographic column was as follows: 60°C for 1 min, followed by increasing at 5°C/min to 115°C, then at 15°C/min to 180°C. The temperature of the detector was 260°C, and the temperature of the injector was 250°C. MCP concentrations were measured using a kit (Nanjing Jiancheng Institute of Biological Engineering) with a microplate reader (Thermo Fisher Scientific).
Calculation of ADFI, ADG, and FCR
ADG (kg/d) = [Weight at the end of the experiment (56 d)-Weight at the beginning of the experiment (10 d)]/experiment period (56 d-10 d);
ADFI (kg/d) = total feed intake during the experiment/experiment period; FCR = ADFI (kg/d)/ADG (kg/d).
Preparation and observation of tissue sections
The rumen ventral sac tissue fixed with formaldehyde solution was dehydrated, trimmed, embedded, sliced, stained and sealed. The specific operation procedure was as follows: the fixed tissue was dehydrated by an automatic dehydrator (dehydration duration: 75% alcohol, 4 h; 85% alcohol, 2 h; 95% alcohol, 1 h; 100% alcohol, 0.5 h; 100% alcohol, 0.5 h; 100% alcohol, 0.5 h; 100% alcohol, 0.5 h; 100% alcohol, 0.5 h; 100% alcohol, 0.5 h; xylene, 10 min; xylene, 10 min; paraffin, 1 h; paraffin, 2 h; paraffin, 3 h), embedded, and sectioned. Then, the following operations were conducted: dewaxing with xylene, hematoxylin staining for 10–20 min, rinsing with running water for 1–3 min, differentiation (1% hydrochloric acid and 99% anhydrous ethanol) for 5–10 s, rinsing with running water for 1–3 min, placing into warm (50°C) water or weakly alkaline aqueous solution to turn the solution blue, and rinsing with running water for 1–3 min. Then 85% alcohol was added for 3–5 min, eosin staining was conducted for 3–5 min, and then the samples were rinsed with running water for 3–5 s before alcohol dehydration. Then, the samples were subjected to xylene dewaxing till transparent and transferred to a neutral gum sealing sheet. The slices were dried and photographed with a light microscope (Nikon, Japan). Nipple length, nipple width and base thickness were then measured using Image-Pro Plus software (Media Cybernetics, Bethesda, MD, USA).
Extraction and detection of total DNA and metagenomic sequencing of rumen contents
The total DNA was extracted by CTAB method after grinding the samples with tissue grinder (Jingxin, Shanghai). Total DNA of the rumen contents was extracted by the CTAB method as follows: (1) a total of 1,000 μl CTAB lysate was added to a 2.0 ml EP tube, 20 μl of lysozyme was added, an appropriate amount of sample was added into the lysate, and the sample was placed in a 65°C water bath (for rumen content samples, the time in the water bath was 2 h). The mixture was homogenized several times during this period to achieve complete lysis. (2) A total of 950 μl of the supernatant was centrifuged, and the same volume of 25:24:1 phenol (pH = 8.0): chloroform:isopentyl alcohol was added. The sample was mixed by inversion and centrifuged at 12,000 rpm for 10 min. (3) The supernatant was collected, isoamyl alcohol (24:1) was added, and the mixture was mixed by inversion and then centrifuged at 12,000 rpm for 10 min. (4) The supernatant was transferred to a 1.5 ml centrifuge tube, 3/4 the volume of isopropyl alcohol was added to the supernatant, and then the mixture was shaken and precipitated at −20°C. (5) The mixture was centrifuged at 12,000 rpm for 10 min, the liquid was removed while taking care to retain the precipitate, the precipitate was washed twice with 1 ml of 75% ethanol, and the remaining small amount of liquid was collected by centrifugation again and removed with a pipette. (6) The sample was blown dry on an ultraclean table or air dried at room temperature (DNA samples were not completely dried to aid in resuspension). (7) A total of 51 μl of ddH2O was added to dissolve the DNA samples, which were incubated at 55–60°C for 10 min for dissolution if necessary. (8) A total of 1 μl of RNase A was added to digest RNA and placed at 37°C for 15 min. The purity and integrity of DNA were analyzed by agarose gel electrophoresis (AGE). A Qubit was used to accurately quantify the DNA concentration. For qualified DNA samples, a Covaris ultrasonic fragmentation instrument was used to randomly interrupt fragments with a growth degree of 350 bp, and the whole library was prepared by terminal repair, adding an A tail and sequencing adaptors, and subjected to purification and PCR amplification. Illumina PE150 sequencing was performed by pooling different libraries according to effective concentration and target data volume.
Bioinformatic analyzes
Based on raw data returned by the Illumina PE150 sequencing platform, data quality control was first carried out: removal of reads with low-quality bases (mass value ≤38) exceeding a certain percentage (default: 40 bp), removal of a certain proportion of reads containing N bases (the default was 10 bp), and removal of reads with overlaps with the adaptor exceeding a certain threshold (the default was 15 bp) were conducted. If the sample was contaminated by the host, it needed to be compared with the host sequence, and the filtered samples could come from the host (15, 16). After pretreatment, clean data were obtained. Second, SOAPdenovo was utilized (17). Assembly software was used for assembly analysis. The assembled scaffolds were interrupted from the N connection to obtain the n-free sequences known as scaftigs (18) (i.e., continuous sequences within Scaffolds). Bowtie 2 software was used to compare clean data of each sample after quality control to scaftigs of each sample after assembly to obtain unused PE reads. Reduced properties were obtained with the following parameters (15, 16): –end-to-end, –sensitive, -I 200, -x 400; unutilized reads of each sample were combined, and k-MER = 55 was selected for mixed assembly (15, 19, 20). Other assembly parameters were the same as those of a single sample. The mixed assembly Scaffolds were interrupted from the N connector to obtain the scaftigs sequences without N. For scaftigs generated from single samples and mixed assemblies, fragments smaller than 500 bp (15, 18, 21, 22) were filtered out, and statistical analysis and subsequent gene prediction were performed.
Starting from a single sample and scaftigs after mixed assembly, MetaGeneMark was used for gene prediction, and genes generated by each sample and mixed assembly prediction were combined to remove redundancy and construct a gene catalog. Starting from the gene catalog (16, 22, 23) and integrating the clean data of each sample, the abundance information of the gene catalog in each sample could be obtained for basic information statistics and core-PAN gene analysis. Species annotation was conducted as follows: DIAMOND software was used to compare genes with various functional databases, and then the LCA algorithm was used to obtain species annotation information (15, 24). Commonly used functional databases were compared, and based on the gene catalog, functional annotation and abundance analysis of metabolic pathways (KEGG), homologous gene clusters (eggNOG) and carbohydrate enzymes (CAZy) were carried out (15, 25–27). Based on the species abundance and functional abundance tables, abundance clustering analysis, NMDS dimensionality reduction analyses, Metastat and linear discriminant analysis effect size (LEfSe) multivariate statistical analyses and metabolic pathway comparative analysis were conducted to mine the differences in species composition and functional composition among samples.
Data analysis
SPSS software (SPSS Version 25.0, SPSS, Inc.) was used for variance analysis of rumen fermentation parameters and rumen papilla development data; was P < 0.05 is considered significant, and 0.05 < P < 0.1 was considered a significant trend of difference.
Results and analysis
Determination of rumen fermentation parameters
As shown in Table 2, the ADFI, ADG, and MCP in group A was significantly higher than those in group C (P < 0.05); the FCR, NH3-N, and pH showed no significant difference between groups A and C (P > 0.05); the papilla width and muscular thickness of the rumen in group A were significantly higher than those in group C (P < 0.05), and the papilla length in group C was higher than that in group A but not significantly (P > 0.05). The acetic acid percentage in group A was significantly higher than that in group C (P < 0.05), impacting the ratio of acetic acid to propionic acid (A/P).
Table 2. Effects of different NDF/starch starter on production performance, rumen fermentation parameters, and rumen papilla development in lambs.
Gene prediction results and abundance analysis
Open reading frame (ORF) prediction of scaftigs (≥500 bp) was performed using MetaGeneMark and the prediction results and basic information of the gene catalog obtained are shown in Table 3. A total of 1,315,345 ORFs were predicted, with an average of 109,612 ORFs per sample. After redundancy removal, a total of 407,675 genes in the gene catalog were obtained for subsequent species annotation and functional analysis. The total length of the non-redundant gene catalog was 331.68 Mbp, the average length was 813.58 bp and the GC content was 47.93%.
Table 3. Gene catalog.
A Venn diagram (Figure 1A) showed 266,703 common genes between the two groups, 24,818 unique genes in group A and 101,856 unique genes in group C. Based on the abundance information of each gene in the gene catalog of each sample, core gene analysis (Figure 1B) was performed. As the number of samples increased, the number of genes in sample combinations gradually stabilized, and additional samples only increased the number of genes in a small number of sample combinations, indicating that the sequencing was reasonable and could meet the requirements of subsequent analysis. The NMDS plot (Figure 1C) shows that the two groups of samples are located on both sides of the horizontal axis, and the distance between the samples in each group is close, indicating that there is a significant difference between the two groups. To investigate the different gene numbers between groups, a box diagram of gene number difference between groups was drawn, and the results are shown in Figure 1D.
Figure 1. Venn diagram of number of genes (A); Core gene dilution curve (B); NMDS plot based on phylum level (C); Box plot of gene number difference between groups (D).
Species annotation and differential analysis
In this study, DIAMOND software was used to compare unigenes with each functional database. For the comparison results of each sequence, the results with the highest score were selected for subsequent analysis, and the relative abundance of different functional levels was calculated. The unigenes for each sample (407,675) were compared with the NR database, and the number of genes that could be annotated in the NR database was 343,192 (84.18%). In the NR database that could be annotations, the proportions of genes annotated at the bacterial kingdom, phylum, class, order, family, genus and species levels were 90.07, 86.29, 82.88, 82.40, 74.59, 69.91, and 47.75%, respectively. For the top 10 species in terms of relative abundance at the phylum and genus levels, a histogram of relative abundance was generated (Figure 2). Figure 2 shows that Bacteroidetes (43.07%), Firmicutes (34.25%), Actinobacteria (2.82%) and Chlamydiae (1.17%) were the main Bacteroides in group A. Bacteroidetes (38.46%), Firmicutes (35.32%), Proteobacteria (2.58%), and Actinobacteria (1.98%) were the main Bacteroidetes in group C. Prevotella (21.34%), Alistipes (6.42%,) Bacteroides (3.98%), Clostridium (3.16%), and Dialister (2.76%) were the main genera in group A. Prevotella (24.99%), Clostridium (9.14%), Bacteroides (2.68%), and Alistipes (6.42%) were the main genera in group C. These were the top 10 bacteria in relative abundance at the genus level, among which the abundance of seven species in group A was higher than that in group C, and the abundance of the other three species in Group C was higher than that in group A. The abundance of Alistipes and Bacteroides in group A was significantly higher than that in group C (P < 0.05); the abundance of Clostridium was significantly higher in group C than in group A (P < 0.05).
Figure 2. Histogram of horizontal relative abundance of phylum (A) and genus (B). The horizontal axis represents the sample name; the vertical axis represents the relative proportion of species annotated to a particular type; the legend on the right shows the species category corresponding to each color block.
LEfSe analysis of differentially abundant species between groups
LEfSe analysis (Figure 3) showed that the biomarker of group A and group C were significantly different (LDA score > 4), there were 14 differentially abundant biomarkers in group A, and there were significant differences between groups A and C.
Figure 3. The left figure shows the LDA value distribution of different species. The histogram of LDA value distribution shows the species whose LDA score is greater than the set value (the default value is 4), that is, the Biomarker with statistical difference between groups. The length of the bar graph represents the impact of different species (LDA score).
Metastat analysis of intergroup differentially abundant species
As shown in Figure 4, there were four phyla with significant differences between group A and group C at the phylum level. The relative abundances of Chloroflexi, Candidatus Peregrinibacteria, Thermotogae, and Candidatus Doudnabacteria in group C were both significantly higher than those in group A (P < 0.05).
Figure 4. Phylum level differential species. The horizontal axis is the sample grouping, and the vertical axis is the relative abundance of corresponding species; the horizontal lines represent two groups with significant differences. *Represents significant difference between the two groups (q < 0.05), **Represents extremely significant difference between the two groups (q < 0.01).
KEGG database functional annotation and abundance analysis
Unigenes (407,675) of each sample were matched with the KEGG database, 261,718 (64.20%) genes were matched with the KEGG database, and 141,546 (34.72%) genes were matched with 3,891 Kos (KEGG ortholog groups). The statistical diagram of the number of genes annotated into KEGG metabolism pathways in rumen microbes of lambs in group A and group C is shown in Figure 5. As shown in the figure, metabolism-related genes accounted for a large proportion of the genes, among which carbohydrate metabolism was the most abundant functional category, and 16,327 genes were annotated in this category. The abundance of genes involved in amino acid metabolism was high, with 13,743 genes being annotated in this category. Nucleotide metabolism, metabolism of cofactors and vitamins and energy metabolism were also abundant. There were 8,982 genes annotated for cellular processes, 11,447 genes for environmental information processing, 19,747 genes for genetic information processing, 6,435 genes for human diseases, and 3,095 genes for organismal systems. In addition, pathway with the highest gene abundance in environmental information processing was membrane transport, which had 8,174 genes. In genetic information processing, translation had the highest gene abundance, with 8,807 genes.
Figure 5. (A) Statistical map of gene number. The number of genes annotated into the KEGG metabolic pathways in rumen microorganisms of the two groups was calculated. (B) Functional annotation abundance graph (level 1).
After analyzing the abundance of each function at level 1 (Figure 5), in group A metabolism accounted for 13.16% of gene function, genetic information processing accounted for 4.99%, environmental information processing accounted for 2.60%, cellular processes accounted for 1.84%, human diseases accounted for 1.40%, and organismal systems accounted for 0.63%. In group C metabolism accounted for 12.56% of gene function, genetic information processing accounted for 5.40%, environmental information processing accounted for 2.41%, cellular processes accounted for 1.83%, human diseases accounted for 1.43%, and organismal systems accounted for 0.64%. These were not statistically significant differences in abundance (P < 0.05).
Functional annotation and abundance analysis based on the eggNOG database
The functional annotation results of the eggNOG database are shown in Figure 6A, and the main functional genes included 12,173 genes for energy production and conversion; 18,668 genes for amino acid transport and metabolism; 17,768 genes for carbohydrate transport and metabolism; 16,301 genes for translation, ribosomal structure and biogenesis; 12,592 genes for transcription; 22,471 genes for replication, recombination and repair; 17,671 genes for cell wall/membrane/envelope biogenesis; and 11,935 genes for inorganic ion transport and metabolism. In addition, unknown functional genes accounted for a large proportion of the total genes.
Figure 6. (A) Statistics of the number of gene annotations based on the eggNOG database; (B) the histogram of functional abundance on level 1 based on the KEGG database.
The functional abundance of lamb intestinal microbes was analyzed at level 1, and the analysis results are shown in Figure 6B. In group A, genes of unknown function accounted for 20.58% of gene function; replication, recombination and repair accounted for 7.25%; carbohydrate transport and metabolism accounted for 4.66%; translation, ribosomal structure and biogenesis accounted for 4.73%; amino acid transport and metabolism accounted for 4.71%; cell wall/membrane/envelope biogenesis accounted for 4.44%; transcription accounted for 3.24%; energy production and conversion accounted for 3.01%; inorganic ion transport and metabolism accounted for 3.03%; and posttranslational modification, protein turnover, chaperones accounted for 2.07%. For group C, genes of unknown function accounted for 19.99% of gene function; replication, recombination and repair accounted for 7.08%, carbohydrate transport and metabolism accounted for 4.00%, translation, ribosomal structure and biogenesis accounted for 5.01%; amino acid transport and metabolism accounted for 4.15%; cell wall/membrane/envelope biogenesis accounted for 4.44%; transcription accounted for 3.15%; energy production and conversion accounted for 3.01%; inorganic ion transport and metabolism accounted for 2.74%; and posttranslational modification, protein turnover, chaperones accounted for 2.15%. There was no significant difference between the two groups (P > 0.05).
CAZy database functional annotation and abundance analysis
We annotated the carbohydrate enzymes contained in groups A and C, and the annotation results are shown in Figure 7A. Seven genes were attributed to auxiliary activities (AAs), 1,355 genes to carbohydrate-binding modules (CBMs), 1,000 genes to carbohydrate esterases (CEs), 8,237 genes to glycoside hydrolases (GHs), 3,980 genes to glycosyl transferases (GTs), and 293 genes to polysaccharide lyases (PLs). The functional abundance of rumen microorganisms in lambs was analyzed at level 1, and the analysis results are shown in Figure 7B.
Figure 7. (A) Statistics of the number of gene annotations based on the CAZy database; (B) the histogram of functional abundance on level 1 based on the CAZy database.
In group A, AA genes accounted for 0.0003% of the abundance, CBM genes accounted for 0.27%, CE genes accounted for 0.22%, GH genes accounted for 2.25%, GT genes accounted for 1.08%, and PL genes accounted for 0.05%. In group C, AA genes accounted for 0.0026% of the abundance, CBM genes accounted for 0.29%, CE genes accounted for 0.25%, GH genes accounted for 1.97%, GT genes accounted for 0.99%, and PL genes accounted for 0.08%. AA genes in group C were significantly enriched in group A (P < 0.05), while the abundances of other genes were not significant between the two groups (P > 0.05).
Differential function analysis
Based on MetaStats analysis, we screened genes with P < 0.05 and |log2(FC)| ≥ 1 as significantly different genes. Table 4 shows the statistics of significantly different genes in the two groups. At KEGG level 2, the cell motility of group C was significantly higher than group A (P < 0.05). At CAZy level 2, GT39, GT1, CBM58, and CBM66 were significantly higher in group A than in group C (P < 0.05), and PL16, CBM77, GT84, CBM61, CBM51, GT39, GH84, PL9, CE8, CBM42, GH119, CBM69, and PL10 were significantly more abundant in group C (P < 0.05).
Table 4. Different functions of rumen microbiota at level 1 or level 2 based on KEGG and CAZy database.
MetaStats analysis found that there was no significant difference in metabolic pathways between groups A and C. In the metabolic pathway diagram (Supplementary Figure 1) of starch and sucrose metabolism (map00500), maltose phosphorylase [EC:2.4.1.8] is a unique enzyme to group C and is produced by Faecalibacterium sp. CAG:74_58_120.
Discussion
Effects of starters with different NDF/starch ratios on the production performance, rumen fermentation characteristics, and rumen tissue morphology of Hu sheep
In this experiment, we found that starter feeding in group A improved the ADG and ADFI of lambs, and there was no difference in FCR between the two groups, but the lambs in group A needed for weight gain, and the lambs in group A needed more concentrate feed, so the starter in group C has more economic value in terms of feed benefits. The pH, NH3-N content, and VFA content are important indicators of rumen fermentation in ruminants (28). pH can be affected by many factors, such as high doses of starch (29, 30). Starch and cellulose are the main dietary components of ruminants, and the degradation of starch and cellulose in ruminants is the key to their high yield (31). In this study, the low starch content in Group C led to a high pH value, which was consistent with the results of Liu et al. (28). The NH3-N concentration comprehensively reflected the degradation of dietary nitrogen in the rumen and the utilization of ammonia by microorganisms, mainly reflecting the utilization of nitrogen by ruminants (32–34). Contrary to a previous study (28), in this study, high concentrate feed levels did not increase the NH3-N concentration, but Prevotella degraded peptides in the rumen and therefore may promote the formation of NH3-N, which may be the reason for the high NH3-N concentration in Group C. Many studies have shown that the acetic acid concentration in diets with high concentrate levels is lower and the propionic acid concentration is significantly increased (28, 35). The results of this study contradict these findings, since animals with high feed efficiency are generally thought to produce more VFAs and less methane (36, 37), while the high acetic acid levels in group A is more likely to cause methane emissions, which is not conducive to feed conversion ratio. At the same time, the decrease in the A/P-value of group C reflects the improvement in the energy utilization efficiency of group C starter. Rumen development affects nutrient digestion and performance in adulthood, and can be measured from the developmental status (length and width) of the rumen papilla, affecting the digestive capacity and yield of adult lambs (38). In this study, the rumen papilla length in group A was lower than that in group C, indicating that group C starter can promote the growth of rumen papillae, which is consistent with our previous research results showing that high NDF feeding can promote rumen papilla development (33).
Effects of starters with different NDF/starch ratios on rumen microbial species abundance in Hu lambs
Metagenomic sequencing analysis revealed the complex and diverse microbial population in the rumen. Recent studies have shown that rumen microbial community diversity has a direct impact on ruminant performance and the rumen environment (39, 40). Ruminants rely on rumen microorganisms to degrade roughage, which provides energy and protein nutrients for the maintenance of growth and lactation in host animals (41). The rumen microflora structure of lambs fed starters with different NDF/starch ratios was analyzed to understand the relationship between host and microbial community, which is conducive to screening the starter for early development of lambs. There was no significant difference in the horizontal abundance ofphyla between groups A and C. Four main phyla existed in the rumen, more than 70% of which was accounted for by Bacteroidetes, which was the most abundant, followed by Firmicutes, Actinobacteria, Proteobacteria, and Chlamydia, which is consistent with previous reports (42–44). In this study, Bacteroidetes and Firmicutes had the highest microbial abundance in the rumen of lambs, which was consistent with previous reports (45, 46).
The species annotation showed that the abundance of Bacteroidetes in group A was higher than that in group C. Bacteroidetes can degrade carbohydrates (47) and are mainly responsible for protein hydrolysis, carbohydrate degradation and amino acid fermentation to acetic acid. It is concluded that the NDF/starch ratio of group A as more effective than that of group C in promoting the growth of Bacteroidetes in the rumen of lambs and in degrading carbohydrates in feed. However, the abundance of Firmicutes in group C was higher than that in group A, which played an important role in energy conversion (47, 48). Therefore, the diet in group C could promote the propagation of Firmicutes and increase the energy conversion efficiency. Consistent with the results of this study, Prevotella has been reported as the most abundant rumen genus (49), and the abundance of Prevotella in group A and group C were 21.34 and 24.99%, respectively. Prevotella was positively correlated with VFA concentration, indicating that Prevotella plays an important role in VFA biosynthesis (50). Although this has not been fully demonstrated in the microbial profile of feed efficiency, Prevotella has been shown to be more abundant in inefficient animals (51). The study also showed that Prevosella abundance was related to propionic acid production and positively correlated with NH3-N (52) content, which was also consistent with our experimental results. Because Prevotella degrades peptides in the rumen, it may promote the formation of NH3-N (53).
Among the four genera with significant differences at the phylum level, Chloroflexi is known to be related to methanogenesis (54), but Chloroflexi has a low abundance in the rumen and produces limited methane compared to methanogens, which can produce methane using acetic acid. Members of Candidatus Peregrinibacteria are involved in carbohydrate metabolism and amino acid synthesis (55), which may be because high NDF starter promoted the mass propagation of microorganisms involved in carbohydrate metabolism. Thermotogae members can grow on simple and complex carbohydrates and have a high H2 production capacity, thereby reducing greenhouse gas emissions (56), which inhibits Chloroflexi methanogenesis. However, the function of Candidatus Doudnabacteria is currently unclear.
Functional research based on different databases
We used comparisons of different databases to predict the function of the lamb rumen microbiome. The results showed that the rumen microbial genes of lambs in the two groups were mainly involved in carbohydrate metabolism and amino acid metabolism, and there was no significant difference in the metabolic process between the two groups, which may be because in addition to the feeding environment, the type and source of starter were also important influencing factors. In addition, compared with previous studies (57, 58), the current study indicate that the overall function of the rumen bacterial community is mainly related to membrane transport, carbohydrate metabolism, amino acid metabolism, replication and repair, translation and energy metabolism. Our data also showed that the difference in rumen bacterial community function between the two groups mainly occurred in metabolism. Specifically, compared with group A, there were more genes related to genetic information processing in group C, and the abundance of genes related to replication and repair was significantly higher than that in group A, while genes related to metabolism and environmental information processing were enriched in group A. Therefore, we speculated that some bacteria in group C regulate their own growth and cell differentiation through the growth and apoptosis of endocrine hormone cells to meet the needs of animals for rich nutrient metabolism and that rumen bacteria in group A can help lambs accumulate more energy.
Carbohydrates are important organic compounds and the main energy source for living cells. Carbohydrates in the gastrointestinal tract play an important role in providing nutrition for hosts and microorganisms or regulating the complex relationship between them (59). Carbohydrate-active enzymes can decompose macromolecular carbohydrates, and various carbohydrate-active enzymes cooperate with each other to degrade oligosaccharides and polysaccharides. The more complex the polysaccharide is, the greater the number of carbohydrate-active enzymes required (60). Studies have shown that the abundance of the gastrointestinal flora increases with increasing dietary fiber intake (61), and a large amount of dietary fiber cannot be digested and decomposed by lambs themselves, requiring a series of enzymes to synergistically aid in degradation. For example, the main chain of polysaccharides is mainly degraded by GHs and partially degraded PLs, and the side chain is degraded by hydrolases such as xylosidase mannosidase and CEs (62). CBMs themselves do not show enzyme activity but help in the binding of GHs to polysaccharides and enhance their activity (63, 64). CE (65) is responsible for the cleavage of glycosidic bonds, and PLs (66) and AAs (67), which target insoluble polymers, are in the CAZymes class with relatively little annotation in the metagenome studied. PL9 and PL10 are members of the pectinase family (68), members of which catalyze the hydrolysis of A-1, 4-glycosidic bonds (69, 70) and play an important role in the degradation and modification of pectin (71, 72). Glycosyltransferases (GTs), as natural biocatalysts, can catalyze the transfer of glycosylates from activated sugar donors to different receptors (72, 73). The GT1 family is a reversal enzyme (74). GT1 family enzymes are common in most organisms, such as bacteria and animals (75), but members from different fields show different functions. For example, “antibiotic glycosylation” is common in microorganisms (76). The GT1 family plays a crucial role in antibiotic biosynthesis and antibiotic resistance (75). In view of these results, the enrichment of CAZyme genes (GH, CE, PL, AA, and CBM) encoding carbohydrate decomposition in the rumen microbial community of group C demonstrated that the group C lambs had a better ability to degrade complex substrates.
In starch and sucrose metabolism, both groups had multiple metabolic pathways and 56 enzymes. The unique enzyme in group C was maltose phosphorylase [EC:2.4.1.8], which is mainly responsible for the degradation of maltose into glucose (77), and this additional metabolic pathway is beneficial for host decomposition of dietary nutrients.
Conclusions
In conclusion, in group C, with an NDF/starch ratio of 1.0, the rumen papillae of lambs were well-developed, the energy utilization efficiency was improved, and CAZyme genes were enriched in the rumen microbial community of lambs. Moreover, group C had greater economic value from the perspective of feed benefits.
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 Faculty Animal Policy and Welfare Committee of Gansu Agricultural University.
Author contributions
HZ was responsible for pilot implementation, sample collection and analysis, manuscript preparation, and manuscript submission and revision. FL was involved in data analysis, statistical analysis, manuscript language revision, journal selection, and manuscript submission and revision. GL participated in the experimental design, experimental implementation, and sample collection. XP performed sample analysis, scanning electron microscopy, and data collection and analysis. XW contributed to supervision, management, assisted students in managing animals, and collecting and analyzing samples. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Natural Science Foundation project of Gansu effects of different NDF levels on intestinal function in lambs (20JR10RA538).
Acknowledgments
We thank the staff at the sheep farm for providing us with the infrastructure. We also thank the staff of Defu Agricultural Technology Co., Ltd. in Minqin County, Gansu province, for providing infrastructure. We thank Gansu Runmu Biological Engineering Co., Ltd. for completing the feed processing work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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/fvets.2023.1064774/full#supplementary-material
References
1. Khan MA, Bach A, Weary DM, Keyserlingk MAGV. Invited review: Transitioning from milk to solid feed in dairy heifers. J Dairy Sci. (2016) 99:885–902. doi: 10.3168/jds.2015-9975
2. Saro C, Hohenester UM, Bernard M, Lagrée M, Martin C, Doreau M, et al. Effectiveness of interventions to modulate the rumen microbiota composition and function in pre-ruminant and ruminant lambs. Front Microbiol. (2018) 9:1273. doi: 10.3389/fmicb.2018.01273
3. Malmuthuge N, Guan LL. Understanding host-microbial interactions in rumen: Searching the best opportunity for microbiota manipulation. J Anim Sci Biotechnol. (2017) 8:8. doi: 10.1186/s40104-016-0135-3
4. Holmes E, Kinross J, Gibson GR, Burcelin R, Jia W, Pettersson S, et al. Therapeutic modulation of microbiota-host metabolic interactions. Sci Transl Med. (2012) 4:137rv6. doi: 10.1126/scitranslmed.3004244
5. Norouzian MA, Valizadeh R, Vahmani P. Rumen development and growth of Balouchi lambs offered alfalfa hay pre- and post-weaning. Trop Anim Health Prod. (2011) 43:1169–74. doi: 10.1007/s11250-011-9819-z
6. Xie B, Huang W, Zhang C, Diao Q, Cui K, Chai J, et al. Influences of starter NDF level on growth performance and rumen development in lambs fed isocaloric and isonitrogenous diets. J Anim Sci. (2020) 98:skaa093. doi: 10.1093/jas/skaa093
7. Nemati M, Amanlou H, Khorvash M, Mirzaei M, Moshiri B, Ghaffari MH. Effect of different alfalfa hay levels on growth performance, rumen fermentation, and structural growth of Holstein dairy calves. J Anim Sci. (2016) 94:1141–8. doi: 10.2527/jas.2015-0111
8. Sun D, Li H, Mao S, Zhu W, Liu J. Effects of different starch source of starter on small intestinal growth and endogenous GLP-2 secretion in preweaned lambs. J Anim Sci. (2018) 96:306–17. doi: 10.1093/jas/skx029
9. Parmar NR, Solanki JV, Patel AB, Shah TM, Patel AK, Parnerkar S, et al. Metagenome of Mehsani buffalo rumen microbiota: an assessment of variation in feed-dependent phylogenetic and functional classification. J Mol Microbiol Biotechnol. (2014) 24:249–61. doi: 10.1159/000365054
10. Noel SJ, Attwood GT, Rakonjac J, Moon CD, Waghorn GC, Janssen PH. Seasonal changes in the digesta-adherent rumen bacterial communities of dairy cattle grazing pasture. PLoS ONE. (2017) 12:e0173819. doi: 10.1371/journal.pone.0173819
11. Lyons T, Boland T, Storey S, Doyle E. Linseed oil supplementation of lambs' diet in early life leads to persistent changes in rumen microbiome structure. Front Microbiol. (2017) 8:1656. doi: 10.3389/fmicb.2017.01656
12. Salami SA, Valenti B, Bella M, O'Grady MN, Luciano G, Kerry JP, et al. Characterisation of the ruminal fermentation and microbiome in lambs supplemented with hydrolysable and condensed tannins. FEMS Microbiol Ecol. (2018) 94:fiy061. doi: 10.1093/femsec/fiy061
13. Wallace RJ, Rooke JA, McKain N, Duthie CA, Hyslop JJ, Ross DW, et al. The rumen microbial metagenome associated with high methane production in cattle. BMC Genom. (2015) 16:839. doi: 10.1186/s12864-015-2032-0
14. Wang Q, Gong J, Huang X, Yu H, Xue F. In vitro evaluation of the activity of microencapsulated carvacrol against Escherichia coli with K88 pili. J Appl Microbiol. (2009) 107:1781–8. doi: 10.1111/j.1365-2672.2009.04374.x
15. Karlsson FH, Tremaroli V, Nookaew I, Bergström G, Behre CJ, Fagerberg B, et al. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature. (2013) 498:99–103. doi: 10.1038/nature12198
16. Karlsson FH, Fåk F, Nookaew I, Tremaroli V, Fagerberg B, Petranovic D, et al. Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat Commun. (2012) 3:1245. doi: 10.1038/ncomms2266
17. Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: An empirically improved memory-efficient short-read de novo assembler. Gigascience. (2012) 1:18. doi: 10.1186/2047-217X-1-18
18. Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat Biotechnol. (2014) 32:822–8. doi: 10.1038/nbt.2939
19. Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, et al. Alterations of the human gut microbiome in liver cirrhosis. Nature. (2014) 513:59–64. doi: 10.1038/nature13568
20. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. (2010) 464:59–65. doi: 10.1038/nature08821
21. Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, et al. Ocean plankton. Structure and function of the global ocean microbiome. Science. (2015) 348:1261359. doi: 10.1126/science.1261359
22. Li J, Jia H, Cai X, Zhong H, Feng Q, Sunagawa S, et al. An integrated catalog of reference genes in the human gut microbiome. Nat Biotechnol. (2014) 32:834–41. doi: 10.1038/nbt.2942
23. Oh J, Byrd AL, Deming C, Conlan S, Kong HH, Segre JA. Biogeography and individuality shape function in the human skin metagenome. Nature. (2014) 514:59–64. doi: 10.1038/nature13786
24. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. (2015) 12:59–60. doi: 10.1038/nmeth.3176
25. Huson DH, Mitra S, Ruscheweyh HJ, Weber N, Schuster SC. Integrative analysis of environmental sequences using MEGAN4. Genome Res. (2011) 21:1552–60. doi: 10.1101/gr.120618.111
26. Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, et al. From genomics to chemical genomics: New developments in KEGG. Nucleic Acids Res. (2006) 34:D354–7. doi: 10.1093/nar/gkj102
27. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. (2017) 45:D353–61. doi: 10.1093/nar/gkw1092
28. Liu H, Xu T, Xu S, Ma L, Han X, Wang X, et al. Effect of dietary concentrate to forage ratio on growth performance, rumen fermentation and bacterial diversity of Tibetan sheep under barn feeding on the Qinghai-Tibetan plateau. PeerJ. (2019) 7:e7462. doi: 10.7717/peerj.7462
29. Lettat A, Nozière P, Silberberg M, Morgavi DP, Berger C, Martin C. Experimental feed induction of ruminal lactic, propionic, or butyric acidosis in sheep. J Anim Sci. (2010) 88:3041–6. doi: 10.2527/jas.2010-2926
30. Minuti A, Ahmed S, Trevisi E, Piccioli-Cappelli F, Bertoni G, Jahan N, et al. Experimental acute rumen acidosis in sheep: Consequences on clinical, rumen, and gastrointestinal permeability conditions and blood chemistry. J Anim Sci. (2014) 92:3966–77. doi: 10.2527/jas.2014-7594
31. Hua D, Hendriks WH, Xiong B, Pellikaan WF. Starch and cellulose degradation in the rumen and applications of metagenomics on ruminal microorganisms. Animals. (2022) 12:213020. doi: 10.3390/ani12213020
32. Laarman AH, Oba M. Short communication: Effect of calf starter on rumen pH of Holstein dairy calves at weaning. J Dairy Sci. (2011) 94:5661–4. doi: 10.3168/jds.2011-4273
33. Lv F, Wang X, Pang X, Liu G. Effects of supplementary feeding on the rumen morphology and bacterial diversity in lambs. PeerJ. (2020) 8:e9353. doi: 10.7717/peerj.9353
34. Stanley CC, Williams CC, Jenny BF, Fernandez JM, Bateman HG 2nd, Nipper WA, et al. Effects of feeding milk replacer once versus twice daily on glucose metabolism in Holstein and Jersey calves. J Dairy Sci. (2002) 85:2335–43. doi: 10.3168/jds.S0022-0302(02)74313-0
35. Polyorach S, Wanapat M, Cherdthong A. Influence of Yeast Fermented Cassava Chip Protein (YEFECAP) and roughage to concentrate ratio on ruminal fermentation and microorganisms using in vitro gas production technique. Asian-australas J Anim Sci. (2014) 27:36–45. doi: 10.5713/ajas.2013.13298
36. Li F, Guan LL. Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl Environ Microbiol. (2017) 83:17. doi: 10.1128/AEM.00061-17
37. Shabat SK, Sasson G, Doron-Faigenboim A, Durman T, Yaacoby S, Berg Miller ME, et al. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. (2016) 10:2958–72. doi: 10.1038/ismej.2016.62
38. Lesmeister KE, Heinrichs AJ, Gabler MT. Effects of supplemental yeast (Saccharomyces cerevisiae) culture on rumen development, growth characteristics, and blood parameters in neonatal dairy calves. J Dairy Sci. (2004) 87:1832–9. doi: 10.3168/jds.S0022-0302(04)73340-8
39. Jami E, White BA, Mizrahi I. Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS ONE. (2014) 9:e85423. doi: 10.1371/journal.pone.0085423
40. Schären M, Frahm J, Kersten S, Meyer U, Hummel J, Breves G, et al. Interrelations between the rumen microbiota and production, behavioral, rumen fermentation, metabolic, and immunological attributes of dairy cows. J Dairy Sci. (2018) 101:4615–37. doi: 10.3168/jds.2017-13736
41. Sadan T, Aravindakshan TV, Radhika G, Anand LF, Ally K. Metagenomic analysis exploring taxonomic diversity of rumen microbial communities in Vechur and crossbred cattle of Kerala state, India. J Appl Genet. (2020) 61:287–97. doi: 10.1007/s13353-020-00547-7
42. Jami E, Mizrahi I. Composition and similarity of bovine rumen microbiota across individual animals. PLoS ONE. (2012) 7:e33306. doi: 10.1371/journal.pone.0033306
43. Lopes LD, de Souza Lima AO, Taketani RG, Darias P, da Silva LR, Romagnoli EM, et al. Exploring the sheep rumen microbiome for carbohydrate-active enzymes. Antonie Van Leeuwenhoek. (2015) 108:15–30. doi: 10.1007/s10482-015-0459-6
44. Mao S, Zhang M, Liu J, Zhu W. Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: Membership and potential function. Sci Rep. (2015) 5:16116. doi: 10.1038/srep16116
45. Dodd D, Mackie RI, Cann IK. Xylan degradation, a metabolic property shared by rumen and human colonic Bacteroidetes. Mol Microbiol. (2011) 79:292–304. doi: 10.1111/j.1365-2958.2010.07473.x
46. Fan Q, Wanapat M, Yan T, Hou F. Altitude influences microbial diversity and herbage fermentation in the rumen of yaks. BMC Microbiol. (2020) 20:370. doi: 10.1186/s12866-020-02054-5
47. Thomas F, Hehemann JH, Rebuffet E, Czjzek M, Michel G. Environmental and gut bacteroidetes: The food connection. Front Microbiol. (2011) 2:93. doi: 10.3389/fmicb.2011.00093
48. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. (2006) 444:1027–31. doi: 10.1038/nature05414
49. Stevenson DM, Weimer PJ. Dominance of Prevotella and low abundance of classical ruminal bacterial species in the bovine rumen revealed by relative quantification real-time PCR. Appl Microbiol Biotechnol. (2007) 75:165–74. doi: 10.1007/s00253-006-0802-y
50. Xue MY, Sun HZ, Wu XH, Liu JX, Guan LL. Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance. Microbiome. (2020) 8:64. doi: 10.1186/s40168-020-00819-8
51. Carberry CA, Kenny DA, Han S, McCabe MS, Waters SM. Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle. Appl Environ Microbiol. (2012) 78:4949–58. doi: 10.1128/AEM.07759-11
52. McLoughlin S, Spillane C, Claffey N, Smith PE, O'Rourke T, Diskin MG, et al. Rumen microbiome composition is altered in sheep divergent in feed efficiency. Front Microbiol. (2020) 11:1981. doi: 10.3389/fmicb.2020.01981
53. Wallace RJ. Ruminal microbial metabolism of peptides and amino acids. J Nutr. (1996) 126:1326s−34s. doi: 10.1093/jn/126.suppl_4.1326S
54. Bovio P, Cabezas A, Etchebehere C. Preliminary analysis of Chloroflexi populations in full-scale UASB methanogenic reactors. J Appl Microbiol. (2019) 126:667–83. doi: 10.1111/jam.14115
55. Anantharaman K, Brown CT, Burstein D, Castelle CJ, Probst AJ, Thomas BC, et al. Analysis of five complete genome sequences for members of the class Peribacteria in the recently recognized Peregrinibacteria bacterial phylum. PeerJ. (2016) 4:e1607. doi: 10.7717/peerj.1607
56. Lanzilli M, Esercizio N, Vastano M, Xu Z, Nuzzo G, Gallo C, et al. Effect of cultivation parameters on fermentation and hydrogen production in the phylum thermotogae. Int J Mol Sci. (2020) 22:10341. doi: 10.3390/ijms22010341
57. Liu J, Li H, Zhu W, Mao S. Dynamic changes in rumen fermentation and bacterial community following rumen fluid transplantation in a sheep model of rumen acidosis: implications for rumen health in ruminants. FASEB J. (2019) 33:8453–67. doi: 10.1096/fj.201802456R
58. Pitta DW, Indugu N, Kumar S, Vecchiarelli B, Sinha R, Baker LD, et al. Metagenomic assessment of the functional potential of the rumen microbiome in Holstein dairy cows. Anaerobe. (2016) 38:50–60. doi: 10.1016/j.anaerobe.2015.12.003
59. Hooper LV, Midtvedt T, Gordon JI. How host-microbial interactions shape the nutrient environment of the mammalian intestine. Annu Rev Nutr. (2002) 22:283–307. doi: 10.1146/annurev.nutr.22.011602.092259
60. Martens EC, Lowe EC, Chiang H, Pudlo NA, Wu M, McNulty NP, et al. Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts. PLoS Biol. (2011) 9:e1001221. doi: 10.1371/journal.pbio.1001221
61. Sonnenburg ED, Smits SA, Tikhonov M, Higginbottom SK, Wingreen NS, Sonnenburg JL. Diet-induced extinctions in the gut microbiota compound over generations. Nature. (2016) 529:212–5. doi: 10.1038/nature16504
62. White BA, Lamed R, Bayer EA, Flint HJ. Biomass utilization by gut microbiomes. Annu Rev Microbiol. (2014) 68:279–96. doi: 10.1146/annurev-micro-092412-155618
63. Du L, Ma L, Ma Q, Guo G, Han X, Xiao D. Hydrolytic boosting of lignocellulosic biomass by a fungal lytic polysaccharide monooxygenase, AnLPMO15g from Aspergillus niger. Ind Crops Prod. (2018) 126:309–15. doi: 10.1016/j.indcrop.2018.10.029
64. Bernardes A, Pellegrini VOA, Curtolo F, Camilo CM, Mello BL, Johns MA, et al. Carbohydrate binding modules enhance cellulose enzymatic hydrolysis by increasing access of cellulases to the substrate. Carbohydr Polym. (2019) 211:57–68. doi: 10.1016/j.carbpol.2019.01.108
65. Bohra V, Dafale NA, Purohit HJ. Understanding the alteration in rumen microbiome and CAZymes profile with diet and host through comparative metagenomic approach. Arch Microbiol. (2019) 201:1385–97. doi: 10.1007/s00203-019-01706-z
66. Bertucci M, Calusinska M, Goux X, Rouland-Lefèvre C, Untereiner B, Ferrer P, et al. Carbohydrate hydrolytic potential and redundancy of an anaerobic digestion microbiome exposed to acidosis, as uncovered by metagenomics. Appl Environ Microbiol. (2019) 85:19. doi: 10.1128/AEM.00895-19
67. Singhania RR, Agarwal RA, Kumar RP, Sukumaran RV. Chapter 19: Energy, environment, and sustainability. In: Waste to Wealth. Cellulosic Biomass-Hydrolyzing Enzymes (2018). p. 441–56. doi: 10.1007/978-981-10-7431-8
68. Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucl Acids Res. (2014) 42:D490–5. doi: 10.1093/nar/gkt1178
69. Marín-Rodríguez MC, Orchard J, Seymour GB. Pectate lyases, cell wall degradation and fruit softening. J Exp Bot. (2002) 53:2115–9. doi: 10.1093/jxb/erf089
70. Hugouvieux-Cotte-Pattat N, Condemine G, Shevchik VE. Bacterial pectate lyases, structural and functional diversity. Environ Microbiol Rep. (2014) 6:427–40. doi: 10.1111/1758-2229.12166
71. Patidar MK, Nighojkar S, Kumar A, Nighojkar A. Pectinolytic enzymes-solid state fermentation, assay methods and applications in fruit juice industries: a review. 3 Biotech. (2018) 8:199. doi: 10.1007/s13205-018-1220-4
72. Herron SR, Benen JA, Scavetta RD, Visser J, Jurnak F. Structure and function of pectic enzymes: Virulence factors of plant pathogens. Proc Natl Acad Sci USA. (2000) 97:8762–9. doi: 10.1073/pnas.97.16.8762
73. Blanchard S, Thorson JS. Enzymatic tools for engineering natural product glycosylation. Curr Opin Chem Biol. (2006) 10:263–71. doi: 10.1016/j.cbpa.2006.04.001
74. Liu J, Mushegian A. Three monophyletic superfamilies account for the majority of the known glycosyltransferases. Protein Sci. (2003) 12:1418–31. doi: 10.1110/ps.0302103
75. Quirós LM, Aguirrezabalaga I, Olano C, Méndez C, Salas JA. Two glycosyltransferases and a glycosidase are involved in oleandomycin modification during its biosynthesis by Streptomyces antibioticus. Mol Microbiol. (1998) 28:1177–85. doi: 10.1046/j.1365-2958.1998.00880.x
76. Zhang P, Zhang Z, Zhang L, Wang J, Wu C. Glycosyltransferase GT1 family: Phylogenetic distribution, substrates coverage, and representative structural features. Comput Struct Biotechnol J. (2020) 18:1383–90. doi: 10.1016/j.csbj.2020.06.003
Keywords: starter, lambs, rumen microbes, metagenomics, NDF/starch
Citation: Zhao H, Lv F, Liu G, Pang X, Han X and Wang X (2023) Effects of starters with different NDF/starch ratio on rumen fermentation parameters and rumen microorganisms in lambs. Front. Vet. Sci. 10:1064774. doi: 10.3389/fvets.2023.1064774
Received: 08 October 2022; Accepted: 10 January 2023;
Published: 26 January 2023.
Edited by:
Mirja Wilkens, Leipzig University, GermanyReviewed by:
Melanie Brede, University of Veterinary Medicine Hannover, GermanySayed Haidar Abbas Raza, South China Agricultural University, China
Copyright © 2023 Zhao, Lv, Liu, Pang, Han and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaojuan Wang, 
wangxj@gsau.edu.cn