- Key Laboratory of Animal Nutrition and Feed Science at Universities of Inner Mongolia Autonomous Region, College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
Introduction: This study investigated the effects of dietary energy level on the antioxidant capability, immune function, and rectal microbiota in donkey jennets during the last 60 days of gestation.
Methods: Fifteen pregnant DeZhou donkeys with age of 6.0 ± 0.1 years, body weight of 292 ± 33 kg, parity of 2.7 ± 0.1 parities and similar expected date of confinement (74 ± 4 days) were randomly allocated to three groups and feed three diets: high energy (10.92 MJ/kg, H), medium energy (10.49 MJ/kg, M), and low energy (9.94 MJ/kg, L).
Results and Discussion: The serum activity of catalase (CAT), total superoxide dismutase (T-SOD), glutathione peroxidase (GSH-Px), and total antioxidant capacity (T-AOC) in group M was significantly higher, whereas the concentrations of malondialdehyde (MDA), interleukin 1 (IL-1), IL-2, and IL-6 were lower than those recorded for groups H and L (p ≤ 0.05). The dietary energy level significantly affected rectal microbial community structure in the jennet donkeys 35 days and 7 days before the parturition (p ≤ 0.05). The abundances of norank_f_norank_o_Coriobacteriales genus was significantly higher (p ≤ 0.05) in group H, and the abundances of norank_f_norank_o_Mollicutes_RF39 and the Candidatus_Saccharimonas were higher in group L (p ≤ 0.05). The abundance of Fibrobacter in group M was significantly increased (p ≤ 0.05). The abundance of norank_f_norank_o_Coriobacteriales was positively correlated with average daily gain (ADG) and tumor necrosis factor-α (TNF-α) concentrations (p ≤ 0.05). The abundance of norank_f_norank_o_Mollicutes_RF39 was positively correlated with IL-2 and IL-6 concentrations. The abundance of Candidatus_Saccharimonas was positively correlated with CAT, T-SOD and GSH-Px activities (p ≤ 0.05). The abundance of Fibrobacter was positively correlated with CAT and T-SOD activities (p ≤ 0.05), but negatively correlated with IL-2 concentration (p ≤ 0.05). In conclusion, an appropriate dietary with an energy content of 10.49 MJ/kg for jennet donkeys during late gestation increased the prenatal antioxidant capacity, reduced inflammatory cytokines, and promoted fetal growth, and these changes were related to diet-induced changes in rectal microbiota compositions.
1 Introduction
The market for donkey meat and milk continues growing in China as consumers pursue healthier foods (Liang et al., 2022). However, intensive donkey farming systems are not as well developed as those for other farm animal species. The feeding standards and corresponding supporting management systems have not been fully established. Gestation in donkey jennets, especially in the last quarter of pregnancy (3 months), is a critical period in achieving a successful reproduction outcome. The nutritional status and immune function of pregnant jennets at this stage of pregnancy is closely related to dietary energy level, and can determine foal birth weight. In dairy cows overfeeding during the late stage of pregnancy can result in oxidative stress and inflammation (Zhou et al., 2015), and similar results have been shown in ewes 7 days (d) before parturition (Moradi et al., 2021). However, the effects of dietary energy levels on antioxidant capability and immune function in donkeys during late pregnancy has not been reported.
Gut microbes are an “additional genome” of the host and associated with a host’s health through multiple pathways. Diet is a predominant factor that determines the gut microflora in animals. Changes in gut bacterial species are usually described using the Firmicutes/Bacteroidetes ratio as a marker of microbiome dynamics, and an increase in the Firmicutes/Bacteroidetes ratio is associated with high dietary energy levels in animal and human models (Malesza et al., 2021). Previous research found that a high-energy diet can increase lipopolysaccharides derived from the intestinal microbiota, and these derivatives act as inflammation trigger (Creely et al., 2007). The pro-inflammatory response is mediated by activating the nuclear factor, kappa-B pathway, that is involved in an over-expression of pro-inflammatory cytokines, such as tumor necrosis factor-α (TNF-α), interferon-γ, and inducible NO synthase (Dalvi et al., 2017). Due to the increased expression of the latter, the overproduction of nitric oxide (NO) is followed by the accumulation of reactive nitrogen species, and this is addition to reactive oxygen species (ROS) production (Tapias et al., 2017). Therefore, dietary energy level associated changes in the gut microbiota are closely related to the antioxidant and immune status of animals.
The donkey (Equus africanus asinus) has a well-developed hindgut structure. Therefore, microbial fermentation in the cecum and colon have always played an important role in equine species. Yin et al. (2021) study concluded that a dietary digestible energy (DE) level of 10.49 MJ/kg had the highest nutrient digestibility in female donkeys (Yin et al., 2021). It is hypothesized that higher or lower levels of dietary energy than 10.49 MJ/kg may result in changes in intestinal microflora and antioxidant status that are un-favor for donkey digestion systems.
Generally, the preferred sampling site for hindgut fermentation in animals is directly from the cecum or colon, but this requires animal slaughter or surgical cecum fistulation. In contrast, fecal samples contain most of cecal microbiota, and reflect the bacterial profile in the cecum (Kim and Isaacson, 2015; Stanley et al., 2015). Grimm et al. (2017) reported that changes in microbial communities in the cecum, colon and rectum followed similar patterns in response to dietary changes, thus changes in rectal microbiota can be used as an indicator of bacterial changes, energy metabolism and health of the host. Therefore, fecal samples were collected in this study to determine the correlations between rectal microflora diversity and antioxidant and immune indicators in late pregnant jennets in response to dietary energy levels. Changes in antioxidant capability, immune function, and rectal microbiota were also be measured along with foal birth weights and body size. The results will help to develop hypotheses to model the influence of dietary energy levels on gut microbiota and the healthy status of pregnant donkeys.
2 Materials and methods
The experiment was conducted in the Research Station of Inner Mongolia Agricultural University (Hohhot, China). The animal experimental procedures were performed in accordance with the National Standard Guidelines for Ethical Review of Animal Welfare (GB/T 35892-2018).
2.1 Experimental design, diet, and feeding management
Fifteen pregnant DeZhou jennets with an age of 6 ± 0.1 years, body weight of 292 ± 34 kg, parities of 2.7 ± 0.1 and an expected 74 ± 4 days to foaling were selected and randomly allocated to three group of five animals each. Three groups were then randomly allocated to three diets, containing a DE at 10.92 MJ/kg (H), 10.49 MJ/kg (M), or 9.94 MJ/kg (L) for the rest of gestation period. The ingredients and nutrient compositions of diets are shown in Table 1. The experiment lasted for 74 days (referred to as day −74 in reference to the days before foaling): the initial 14 days was for acclimation on the M diet and also to the environmental conditions, and in the following 60 days the groups were feed their respective diets and the required animal parameters for the 60 day experiment period measured. The jennets were kept in individual pens and fed twice a day at 07:30 and 14:00, respectively. The concentrate to forage ratio of the ration was maintained constant at 30 to 70. The jennets were free access to the diets and water.
2.2 Measurements and sample collection
Feed offered and refusal were recorded daily for the experimental period. The amount of feed offered was adjusted daily to allow for 5% of refusals on the next day, so that voluntary feed intake could be measured. The average dry matter feed intake (DMI) for the whole experimental period was calculated.
The jennets were weighed before the morning feeding on days −60 and − 7 to calculate the average daily gain (ADG) over the 54 day period. Foal birth weight and sizes were recorded immediately after birth. Body length was the distance from the anterior margin of the scapula to the end of the sciatic bone; body height was the vertical distance from the top of withers to the ground; heart girth was the distance around the rib cage directly behind the forelimb; tube girth is the circumference of the finest part of the left forelimb cannon bone.
Fecal samples were collected via the rectum from each animal before the morning feeding on days −35 and − 7 and sample placed into sterile cryopreservation tubes. All samples were immediately snap frozen in liquid nitrogen and stored at −80°C until analysed. Blood samples (each 10 mL) were obtained via jugular puncture before the morning feeding on days −35, −21, −14, −7, and − 1. Serum was collected by centrifuging the blood tubes at 2054 × g for 10 min and the serum stored at −20°C.
2.3 Chemical assays
The serum concentrations of glucose (GLU), triglycerides (TG), cholesterol (CHO) and β-hydroxybutyric acid (BHBA) was analyzed by an automatic biochemical analyzer (Hitachi Co., Ltd., Tokyo, Japan, L-8900) using commercial kits (Lepu Diagnostics Co., Ltd., Beijing, China). The concentrations of non-esterified fatty acids (NEFA) and malondialdehyde (MDA) and the activity of catalase (CAT), total antioxidant capacity (T-AOC), total superoxide dismutase (T-SOD), glutathione, and peroxidase (GSH-Px) was analyzed using commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The concentrations of interleukin1 (IL-1), interleukin 2 (IL-2), interleukin 6 (IL-6), and TNF-α were assayed using enzyme-linked immunosorbent assay (ELISA) kits (Shangbao Biotechnology Co., Ltd., Shanghai, China). All analyses were carried out according to the manufacturer’s protocols.
2.4 16S rRNA gene sequencing and operational taxonomic units picking
Total microbial DNA in five fecal samples from each dietary treatment were extracted using the E.Z.N.A.® soil DNA kit (Omega Bio-tek, Norcross, GA, United States) according to the manufacturer’s protocol. The final DNA concentration and OD values at 260 nm and 280 nm were determined using a NanoDrop 2000 UV–vis spectrophotometer (Thermo Scientific, Wilmington, NC, United States), and DNA quality was checked using 1% agarose gel electrophoresis. Negative extraction controls were included in duplicates to control for extraction contamination. The extracted DNA was stored at −20°C until sequenced.
All PCR amplification and sequencing steps were carried out by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Libraries of 16S rRNA gene sequence were generated using a two-step PCR protocol. The V3-V4 region of the 16S rRNA gene was amplified using the universal bacterial and archaeal primers 338F (5’-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5’-GGACTACHVGGGTWTCTAAT-3′) in the first step of PCR. The PCR reaction mix was 20 μL, including 4 μL of 5× FastPfu buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.4 μL of FastPfu polymerase and 10 ng of template DNA. The PCR product was extracted from agarose gel and further purified using AxyPrep DNA Gel Extraction kit (Axygen Biosciences, Union City, CA, United States) and quantified using QuantiFluor™-ST (Promega, Madison, WI, United States) according to the protocol. Sequencing of 16S rRNA gene amplicons was carried out according to the manufacturer’s protocol using an Illumina MiSeq (2 × 300 bp) and the MiSeq Reagent kit v3 (Illumina).
Raw data was filtered and analyzed using QIIME (Quantitative Insights into Microbial Ecology, version 1.9.1) software, quality-filtered by Trimmomatic and merged by FLASH (Fast Length Adjustment of Short Reads). Low-quality reads were removed using the following criteria: (i) The reads were truncated at any site receiving an average quality score < 20 over a 50 bp sliding window, (ii) primers’ matching allowed 2-nucleotide mismatching, and reads containing ambiguous bases were removed, and (iii) sequences with an overlap longer than 10 bp were merged according to their overlap sequence. The assembled sequences were assigned to operational taxonomic units (OTUs) at 97% similarity cutoff using UPARSE (Highly Accurate OTU Sequences from Microbial Amplicon Reads, version 7.1, http://drive5.com/uparse/, 30 September 2013) and chimeric sequences were identified and removed using UCHIME (Chimera Prediction for Amplicon Sequencing). The OTUs were used for α-diversity indexes (Coverage, Chao, Shannon and Simpson) analysis, and OTUs were taxonomically analyzed by the Ribosomal Database Project Classifier algorithm (http://rdp.cme.msu.edu/, 30 September 2016). The rarefaction curves were analysis with Mothur v.1.21.1 to identify sequence depth. Principal coordinate analysis (PCoA) was measured using the Bray-Curtis distance with the R software suite.
2.5 Statistical analysis
The data of DMI, ADG, birth weight and body sizes were analyzed using the analysis of variance (ANOVA) in SAS (SAS Software, Version 9.1; SAS Institute, Cary, NC, United States) and Duncan’s multiple range test set at the 0.05 level of significance (p ≤ 0.05). Blood biochemical indicators, antioxidant and immune indicators were analyzed in a complete randomized design model for repeated measures using the Proc Mixed procedure in SAS. The Mixed model included fixed effects for treatment (DE), prenatal time (PT, days before foaling) and DE × PT interactions and the random effects of pen. The differences in bacterial diversity indexes (Coverage, Chao, Shannon and Simpson) at phylum and genus classification levels were analyzed using ANOVA and Duncan’s multiple range tests in SAS. Spearman correlations was used to correlate nutrient digestibility, blood antioxidant indexes, and blood immune indexes with the bacteria enrichment at the genus classification level using the pheatmap package in the R software suite. Data for diet nutrient digestibility was obtained from a previous study (Yin, et al., 2021, see attached Supplementary Table S1). The threshold |R| ≥ 0.5 and p ≤ 0.05 were considered as a significant for the Spearman correlations. Data was presented as mean and standard error of means (SEM).
3 Results
3.1 ADG, DMI, birth weight and body size
The initial body weight of the donkey jennets was not different among three groups. DMI (p = 0.042) and ADG (p = 0.014) in group H were significantly higher than those in groups M and L, with no significant differences between groups M and L (Table 2). The differences in the nutrient intakes were consistent with DMI.
Table 2. Effects of dietary energy level on DMI, ADG and nutrient intake of donkey jennets during late gestation.
All foals survived the prenatal feeding regimes. The birth weight (p = 0.020), body length (p = 0.008), and body height for M and L group foals were significantly greater in H group (p = 0.004) (Table 3).
Table 3. Effects of dietary energy level for donkey jennets during late gestation on the birth weight and body sizes of donkey foals.
3.2 Serum biochemical indicators
GLU (p < 0.001) and CHO (p < 0.001) concentrations for groups H and M were significantly higher than group L, the CHO concentration for group H was significantly higher than that of group M (p < 0.001) (Table 4). BHBA (p < 0.001) and NEFA (p < 0.001) concentrations were significantly lower in group M than in groups H and L, and group H was significantly lower than group L. The concentration of TG was significantly higher in group H than in groups M and L (p < 0.01). Prior to parturition, there was a significant increase in TG, GLU, BHBA, and NEFA concentrations (p < 0.05).
Table 4. Effects of dietary energy level on the biochemical indexes in serum of donkey jennets during late gestation.
There were significant interactions between DE and PT for GLU, BHBA, NEFA, and TG concentrations (p < 0.001). GLU was significantly higher on day −1 in all three groups, on day −7 in group M, but lower on day −35 in group L, compared with the other combinations of DE and PT (p < 0.001). BHBA concentrations was greatest on day −1 in group L, but lower on day −14 in group H and d − 35 in group M (p < 0.001). NEFA concentrations were greater on d − 1 in groups H and L, but lower on day −35 in all three groups (p < 0.001). The highest concentration of TG was observed on day −1 in group H, followed by day −7 in group H, on day −1 in groups M and L, with significantly lower concentrations on day −35 in group M (p < 0.001).
3.3 Serum blood antioxidants
The serum concentrations of the measured antioxidant indicators are shown in Table 5. CAT, T-SOD, GSH-Px, and T-AOC concentrations in group M were significantly higher than those in groups H and L (p < 0.001), while MDA concentration showed the opposite pattern. Approaching parturition, CAT, T-SOD, GSH-Px, and T-AOC concentrations decreased significantly, but MDA concentrations showed an increase (p < 0.001).
Table 5. Effects of dietary energy level on the antioxidant indexes in serum of donkey jennets during late gestation.
There were significant interactions between DE and PT and CAT, T-SOD, GSH-Px, and T-AOC serum levels (p < 0.001). CAT activity was greater on days −35 and − 21 in group M, but lower on day −1 in groups H and L, and day −7 in group H (p < 0.001). T-AOC concentration was greater on day −35 in group M, but lower on days −14, −7, and − 1 in all three groups (p < 0.001). T-SOD concentration was greater on day −35 in all three groups, and on day −7 in group M, but lower on day −1 in group L (p < 0.001). GSH-Px serum level was greater on days −35, −21, and − 14 in group M, but lower on day −1 in group H and on days −21, −14, −7, and − 1 in group L (p = 0.003).
3.4 Serum immune indicators
The concentration of IL-1 (p < 0.001) and IL-2 (p < 0.001) in group M were significantly lower than those in group L, and group L were significantly lower than group H (Table 6). The concentration of IL-6 (p < 0.001) and TNF-α (p < 0.001) in group M were significantly lower than in groups H and L, and there was no difference between groups H and L. prior to parturition, the concentration of TNF-α increased significantly (p = 0.001).
Table 6. Effects of dietary energy level on the immune indexes in serum of donkey jennets during late gestation.
There were significant interactions between DE and PT and IL-1, IL-2, and TNF-α concentrations (p < 0.001). The IL-1 concentration was greater on day −14 in group H than on the other days (p < 0.001). The concentration of IL-2 was greater on days −7 and − 1 in group H, and lower on days −35, −14, −7, and − 1 in group M (p < 0.001). Elevated TNF-α concentrations were observed on days −21 and − 7 in group H, and lower concentrations were measured on days −35, −21, −7 and − 1 in group M, on days −35 and − 7 in group H, and on days −35 and -21in group L (p < 0.001).
3.5 Rectal microbiota diversity
3.5.1 Sampling depth
On days −35, 2,673, 2,588, and 2,590 OTUs were separately obtained, on the basis of 97% species similarity, from rectal samples in groups H, M, and L, respectively (Table 7). The average of 61,758, 61,528, and 54,499 sequences were recorded for each sample from the groups H, M and L, respectively. On day −7, 2,530, 2,381, and 2,485 OTUs were obtained, and an average of 56,273, 53,261, and 57,794 sequences were measured for each sample from groups H, M, and L, respectively. Both the rarefaction curves (Figure 1) and the high coverage values (Table 8) obtained showed that sampling depths met the requirements for estimating bacterial diversity.
Figure 1. The OTU rarefaction curves of rectal microbiota in pregnant donkey jennets at 35 day (d-35, A) and 7 day (−7, B) before foaling. Curves were drawn using the least sequenced sample as upper limit for the rarefactions. Red curves and H: high dietary energy (10.92 MJ/kg); Blue curves and M: medium dietary energy (10.49 MJ/kg); Green curves and L: low dietary energy (9.94 MJ/kg).
3.5.2 Rectal microbiota diversity indexes
On day −35, Shannon and Chao indexes were significantly greater in group H than those recorded in group L (p ≤ 0.05), and there was no significant differences in Shannon and Chao indexes between groups H and M, as well as, between groups M and L (Table 8). Simpson index was significantly higher in group L than those in groups H and M (p = 0.017), and there was no significant difference between groups H and M. On day −7, Shannon and Chao indexes were significantly greater in group H than those observed in groups M and L (p ≤ 0.05), and there was no significant difference between groups M and L. Simpson index was significantly greater in group L than that in group H (p = 0.025), and there was no significant difference between groups M and H, and between groups M and L.
A total of 3,057 OTUs were obtained for all the samples on day −35 (Figure 2A), of which 2,159 defined as the core OTUs present in all three groups. The core OTUs consisted of 70.6% of the total OTUs. In addition, 152, 132, and 138 OTUs were uniquely identified in groups H, M and L, respectively. On day −7 (Figure 2B), a total of 2,955 OTUs were obtained for all the samples, of which the core OTUs were 1,968, accounting for 66.6% of the total OTUs. In addition, 191, 132, and 159 OTUs were uniquely identified in groups H, M, and L, respectively. The PCoA plots demonstrated dissimilarities between groups H and M, as well as between groups H and L on day −35 (Figures 2C,D), but no specific clustering was identified between groups M and L. However, the results of rectal microbiota sampled on day −7 showed dissimilarities between groups H, M, and L. Samples in group M occupied the top left of PC1, samples in group H occupied the bottom left of PC1, and samples in group L occupied the right side of PC1.
Figure 2. Venn plot of OTUs showing the percent of observations for each OTU (> 0.5%) (A,B) and Principal coordinate analysis (PCoA) of bacterial OTUs (C,D) in the rectal samples of pregnant donkey jennets at 35 days (d-35, A,C) and 7 days (−7, B,D) before foaling. Red color and H: high dietary energy (10.92 MJ/kg); Blue color and M: medium dietary energy (10.49 MJ/kg); Green color and L: low dietary energy (9.94 MJ/kg).
3.5.3 Taxonomic classifications levels of the bacterial communities
At the phylum level (Table 9), on day −35, Firmicutes (> 56% of the total), Bacteroidetes, Spirochaetes, Kiritimatiellaeota, Actinobacteria, Patescibacteria, Fibrobacteres, Proteobacteria, Verrucomicrobia, Tenericutes, were the predominant phyla (> 99% of total), and overall there were not differences between three treatment groups except for Verrucomicrobia. The abundance of Verrucomicrobia was significantly lower in group M than those in groups H and L (p = 0.013). The abundance of these predominant phyla (> 99% of the total) did not change significantly by day −7. The abundance of Patescibacteria and Tenericutes in group M were significantly lower than those of groups H and L, whereas the abundance of Fibrobacteres was higher (p < 0.05). The abundance of Actinobacteria in group L was significantly lower than in group H (p = 0.049), but no significant differences was found between groups M and H, and between M and L.
Table 9. Effects of dietary energy level on the relative abundance (%) of rectal microbiota at Phylum level in donkey jennets during late gestation.
Linear discriminant analysis (LDA, Figure 3A) used for bacterial family discrimination found the abundances of Eubacteriaceae, Peptococcaceae, unclassified_o_Clostridiales, Leptospiraceae and unclassified_k_norank_d_Bacteria were all significantly elevated on day −35 in group H compared with the abundances in groups M and L. The abundance Akermansiacea, Coriobacteriales_Incertae_Sedis, unclassified_c_Clostridia, and unclassified_c_Actinobacteria families were significantly higher on day −35 in group L compared the abundances in groups H and M.
Figure 3. Bacterial families’ differentiation indicated by linear discriminant analysis (LDA) effect size (LEfSe) in rectal samples of pregnant donkey jennets at 5 days (A) and 7 days (B) before foaling. Red color and H: high dietary energy (10.92 MJ/kg); Blue color and M: medium dietary energy (10.49 MJ/kg); Green color and L: low dietary energy (9.94 MJ/kg).
On day −7 (Figure 3B), the abundances of Ruminococcaceae, norank_o_Coriobacteriales, and unclassified_o_Clostridiales were higher in group H than the abundances in groups M and L. However, the abundance of Clostridiaceae_1, Acidaminococcaceae, Muribaculaceae, Fibrobacteraceae, and Atopobiaceae were significantly higher in group M in comparison with those recorded in groups H and L. The abundances of Rikenellaceae F082, Leptospiraceae, Saccharimonadaceae, unclassified_c_Clostridia, norank_o_Mollicutes_RF39, and norank_o_norank_c_MVP-15, were significantly higher in group L than those for groups H and M.
At the genus level (Figure 4A), the genera Ruminococcaceae_NK4A214_group, Ruminococcaceae_UCG-010, Ruminococcaceae_UC-G-014, Howardella, Rs-H88_termite_group, Eubacterium, unclassified_o_Clostridiales, norank_f_Peptococcaceae, Lachnospiraceae_UCG-008, Johnsonella, and unclassified_k_norank_d_Bacteria were greater on day −35 in group H compared with those in groups M and L. The abundances of Lachnospiraceae_NK4B4_group and Family_XIII_UCG-001 were significantly higher in group M. The abundance of Akkermansia, Slackia, unclassified_c_Actinobacteria, and Phoenicibacte were significantly higher in group L than those recorded for the H and M treatments.
Figure 4. Bacterial genus differentiation indicated by linear discriminant analysis (LDA) effect size (LEfSe) in rectal samples of pregnant donkey jennets at 35 days (A) and 7 days (B) before foaling. Red color and H: high dietary energy (10.92 MJ/kg); Blue color and M: medium dietary energy (10.49 MJ/kg); Green color and L: low dietary energy (9.94 MJ/kg).
On day −7 (Figure 4B), the abundance of 13 genera in group H significantly increased compared with those for groups M and L, and the genera were Ruminococcaceae_NK4A214_group, [Eubacterium]_coprostanoligenes_group, Marvinbryantia, Ruminococcaceae_UCG-014, Catenisphaera, Coprococcus_1, norank_f_Family_XIII, unclassified_o_Clostridiales, Caproiciproducens, norank_f_Christensenellaceae, Coprococcus_3, norank_f_norank_o_Coriobacteriales, and [Eubacterium]_oxidoreducens_group. Groups H and L had significantly lower abundances of Clostridium_sensustricto_1, Phascolarctobacterium, Prevotellaceae_UCG-001, norank_f_Muribaculaceae, Enterorhabdus, Sarcina, Blautia, Fibrobacter, [Eubacterium]_ruminantium_group, and norank_f_Atopobiaceae generas than group M. The abundances of Rikenella-ceae_RC9_gut_group, norank_f_F082, Candidatus Saccharimonas, Candidatus_Soleaferrea, norank_f_norank_o_Mollicutes_RF39, Lachnospiraceae_NK4B4_group, Papillibacter, Mogibacterium, norank_f_norank_o_norank_-c_MVP-15, and norank_f_norank_o_Rhodospirillales in group L were significantly evaluated compared with those abundances for groups H and M.
3.5.4 Spearman correlation analysis
Figure 5 shows that crude protein (CP) digestibility was positively correlated with Clostridium_sensu_stricto_1 (Spearman R = 0.58027, Pcorr = 0.02334) and Sarcina (Spearman R = 0.62366, Pcorr = 0.01298), but negatively correlated with Rikenellaceae_RC9_gut_group (Spearman R = −0.60395, Pcorr = 0.01711) and Candidatus_Saccharimonas (Spearman R = −0.6828, Pcorr = 0.00503). Dry matter (DM) digestibility was positively correlated with [Eubacterium]_coprostanoligenes_group (Spearman R = 0.59857, Pcorr = 0.0184), but negatively correlated with norank_f_norank_o_Bacteroidales (Spearman R = −0.5991, Pcorr = 0.01827). Acid detergent fiber (ADF) digestibility was positively correlated with Clostridium_sensustricto_1 (Spearman R = 0.54857, Pcorr = 0.03422), but negatively correlated with Phoenicibacte (Spearman R = −-0.67389, Pcorr = 0.00587). Neutral detergent fiber (NDF) digestibility was positively correlated with Clostridium sensustricto_1 (Spearman R = 0.66728, Pcorr = 0.00657) and norank_f_Muribaculaceae (Spearman R = 0.55167, Pcorr = 0.03301), but negatively correlated with Phoenicibacte (Spearman R = −-0.58537, Pcorr = 0.02187). The ADG was significantly correlated with norank_f_Muribaculaceae (Spearman R = 0.61649, Pcorr = 0.01438) and norank_f_norank_o_Coriobacteriales (Spearman R = 0.70927, Pcorr = 0.00307), but negatively correlated with norank_f_Bacteroidales_UCG-001 (Spearman R = −0.56631, Pcorr = 0.02774), norank_f_Clostridiales_vadinBB60_group (Spearman R = −0.56808, Pcorr = 0.02715), Rikenellaceae_RC9_gut_group (Spearman R = −0.59857, Pcorr = 0.0184), unclassified_c_Clostridia (Spearman R = −0.71709, Pcorr = 0.00262), Papillibacter (Spearman R = −0.73836, Pcorr = 0.00167), and Lachnospiraceae_NK-4B4_group (Spearman R = −0.61721, Pcorr = 0.01423).
Figure 5. The Spearman correlation analysis of the rectal microbial composition at genus levels with nutrient (DM, EE, ADF, NDF) digestibility and ADG in donkey jennets during late gestation. Red color: positive correlation; Blue color: negative correlation. *: p ≤ 0.05; **: p ≤ 0.01.
Figure 6 shows that CAT activity was positively correlated with Clostridium_sensustricto_1 (Spearman R = 0.5716, Pcorr = 0.02601), Sarcina (Spearman R = 0.55571, Pcorr = 0.03149), norank_f_p-251-o5 (Spearman R = 0.57562, Pcorr = 0.02474), Phascolarctobacterium (Spearman R = 0.56838, Pcorr = 0.02706), and Fibrobacter (Spearman R = 0.61363, Pcorr = 0.01497), but negatively correlated with [Eubacte-rium]_coprostanoligenes_group (Spearman R = −0.51407, Pcorr = 0.04995), Ruminococcaceae_NK4A214_group (Spearman R = −0.60277, Pcorr = 0.01739), unclassified_o_Clostridiales (Spearman R = −0.5245, Pcorr = 0.04472), Candidatus_Saccharimonas (Spearman R = −0.63716, Pcorr = 0.01063), Ruminococcaceae_UCG-014 (Spearman R = −0.61363, Pcorr = 0.01497), Candidatus_Soleaferrea (Spearman R = −0.70595, Pcorr = 0.00327), norank_f_F082 (Spearman R = −0.67759, Pcorr = 0.00551), Rikenellaceae_RC9_gut_group (Spearman R = −0.55933, Pcorr = 0.03017), and Lachnospiraceae_NK4-B4_group (Spearman R = −0.54479, Pcorr = 0.03573). T-SOD was positively correlated with Fibrobacter (Spearman R = 0.55158, Pcorr = 0.03304), but negatively correlated with Ruminococcaceae_NK4A214_group (Spearman R = −0.67456, Pcorr = 0.00581), Ruminococcaceae_UCG-014 (Spearman R = −0.63115, Pcorr = 0.01163), Candidatus_Saccharimonas (Spearman R = −0.53892, Pcorr = 0.03818), unclassified_o_Clostridiales (Spearman R = −0.66973, Pcorr = 0.00631), unclassified_c_Actinobacteria (Spearman R = −0.53748, Pcorr = 0.0388), and Phoenicibacter (Spearman R = −0.5382, Pcorr = 0.03849). GSH-Px activity was positively correlated with Clostridium_sensustricto-1 (Spearman R = 0.64933, Pcorr = 0.0088) and Sarcina (Spearman R = 0.58423, Pcorr = 0.02219), but negatively correlated with Candidatus_Saccharimonas (Spearman R = −0.71685, Pcorr = 0.00263) and norank_f_M2PB4-65_termite_group (Spearman R = −0.5178, Pcorr = 0.04803). MDA concentration was positively correlated with [Eubacterium]_coprostanoligenes_group (Spearman R = 0.62705, Pcorr = 0.002621235), Ruminococcaceae_NK4A214_group (Spearman R = 0.56939, Pcorr = 0.02672), and unclassified_o_Clostridiales (Spearman R = 0.65917, Pcorr = 0.00752), but negatively correlated with Clostridium_sensustricto_1 (Spearman R = −0.51399, Pcorr = 0.04999). T-AOC was positively correlated with norank_f_norank_o_Bacteroidales (Spearman R = 0.53134, Pcorr = 0.04152), but negatively correlated with Lachnospiraceae_AC2044_group (Spearman R = −0.58335, Pcorr = 0.02245).
Figure 6. The Spearman correlation analysis of the rectal microbial composition at genus levels with antioxidant (CAT, T-SOD, GSH-Px, MDA, T-AOC) and immune (IL-1, IL-2, IL-6, and TNF-α) indexes in donkey jennets during late gestation. Red color: positive correlation; Blue color: negative correlation. *: p ≤ 0.05; **: p ≤ 0.01.
IL-2 content was positively correlated with Ruminococcaceae_UCG-014 (Spearman R = 0.56269, Pcorr = 0.02898), unclassified_o_Clostridiales (Spearman R = 0.52686, Pcorr = 0.0436), and norank_f_norank_o_Mollicutes_RF39 (Spearman R = 0.68067, Pcorr = 0.00522), but negatively correlated with norank_f_norank_o_Bacteroidales (Spearman R = −0.63599, Pcorr = 0.01082) and Fibrobacter (Spearman R = −0.5645, Pcorr = 0.02836). IL-6 concentration was positively correlated with Lachnospiraceae_NK4B4_group (Spearman R = 0.63397, Pcorr = 0.01115) and norank_f_norank_o_Mollicutes_RF39 (Spearman R = 0.61512, Pcorr = 0.01466), but negatively correlated with Sarcina (Spearman R = −0.60433, Pcorr = 0.01702) and damnorank_f_Muribaculaceae (Spearman R = −0.56836, Pcorr = 0.02706). TNF-α concentration was positively correlated with [Eubacterium]_coprostanoligenes_group (Spearman R = 0.55368, Pcorr = 0.03225) and norank_f_norank_o_Coriobacteriales (Spearman R = 0.73099, Pcorr = 0.00196), but negatively correlated with norank_f_Bacteroidales_UCG-001 (Spearman R = −0.52302, Pcorr = 0.04544), unclassified_c_Clostridia (Spearman R = −0.7605, Pcorr = 0.001), norank_f_M2PB4-65_termite_group (Spearman R = −0.56829, Pcorr = 0.02709), Papillibacter (Spearman R = −0.58795, Pcorr = 0.02116), and norank_f_Clostridiales_vadinBB60_group (Spearman R = −0.56895, Pcorr = 0.02687).
4 Discussion
In late gestation, fetal growth rate is at its highest, demanding more nutrients from the dam, with glucose consumption by the fetus accounting for 46% of maternal glucose production (Weber et al., 2013). The growth and development of the fetus in the uterus occupies most of the space in the abdominal cavity, suppressing the gestating female’s feed intake (Du et al., 2018). Thus animals in late pregnancy are prone to chronic negative energy balance, leading to an increase in fat mobilization and release of large amounts of NEFA into the blood and liver (White, 2015). When NEFA is incompletely oxidized to ketones, forming BHBA, there is decreased production efficiency, but there is also increased blood ketone concentrations, as well as, as an increased risk of fatty liver disease (Weber et al., 2013). In the present study, the donkey jennets during the last two months of gestation showed an increase in serum concentrations of glucose, BHBA, and NEFA with the progression of gestation, and the concentrations reached the highest values prior to parturition. Our results conform to literature reports that donkey jennets during this late stage of pregnancy are subject to energy deficiency problems. However, the status of the negative energy balance in the donkey jennets can be modulated with increased dietary energy levels, with the jennets on the M diet having greater glucose levels but lower NEFA and BHBA concentrations compared to the donkeys on the H and L diets. These results indicate that feeding donkey jennets during the late stage of gestation with a diet containing an appropriate level of dietary energy, i.e., 10.49 MJ/kg, improved energy utilization in the body and reduced the impact of negative energy balance on serum parameters.
DMI and ADG generally increased with increasing dietary energy levels and in this experiment diet H had higher DMI and ADGs than diets M and L. However, birth weight and body sizes of the foals were the same on the M and L diets, and the foals from the jennets on the H diet had lower birth weights. The results indicate that the energy supplied from the M and L diets plus the body fat mobilization in the jennets could meet the energy demand for fetal growth, even at the lower level of dietary energy supplied in the present study. Studies in beef cattle have shown the importance of pre-partum dietary energy levels for neonatal calves on calf birth weight (Chen et al., 2022) with only severe maternal energy deficiency impairing the birth weight of the offspring. This was not the case in this donkey study.
β-oxidation of NEFA produces a large amounts of ROS, and excessive free radicals breaks the balance between oxidation and antioxidant systems, potentially leading to oxidative stress. Oxidative stress is an important cause of increased inflammatory responses and immune suppression (Sordillo and Aitken, 2009). This is likely to be one of the mechanisms for the increased occurrence of oxidative stress and immune suppression in undernourished animals in late stage of gestation. Therefore, it is crucial to increase dietary energy intake levels for dams in late gestation. However, the inevitable accumulation of excessive energy in the body induces cell enlargement and stress, with the release of inflammatory factors (Maury and Brichard, 2010; Meijer et al., 2011). Thus, an appropriate dietary energy intake during gestation is essential. Liang et al. (2020) reported that dairy cows with high body condition scores pre-partum were more likely to be susceptible to oxidative stress. In the present study, serum concentrations of MDA, IL-1, IL-2, and TNF-α were greater in group H, whereas the activities of CAT, T-AOC, and GSH-Px were lower, compared with the jennets on the M and L diets, indicating that the supply of a high level of dietary energy to the jennets during the late gestation can cause oxidative stress and inflammation. The negative energy balance in the jennets on the low energy diet also resulted in a detrimental effect on antioxidant capacity and immune responses. Therefore, the dietary energy level in group M was the more appropriate diet for donkey jennets during late gestation.
The present study found similar microbial diversities in donkey rectums on days 35 and 7 prior to parturition, as shown by the similar values for Shannon and Chao indexes. However, the microbial diversity changed with the dietary energy levels. The Shannon and Chao indexes in group H were higher, but Simpson index was lower than those in the other two groups, reflecting higher microbial richness and more diversified microbial community in pregnant jennets on the high-energy diet. Similarly, intestinal flora diversity was higher in obese horses compared to the flora diversity in lean horses (Biddle et al., 2018). Previous studies have shown that Firmicutes and Bacteroidetes are the dominant microbial communities in the guts of rodents, swine, horses, and cattle (Al Jassim and Andrews, 2009; Callaway et al., 2010; Holman and Chénier, 2014), and are similar to the results in this study. Firmicutes and Bacteroidetes are the main phyla, respectively, for fiber degradation and carbohydrate catabolism in the gastrointestinal tract of herbivores (Spence et al., 2006; Dearing and Kohl, 2017). Compared with group L, there was a greater abundance of Actinobacteria in group H as well as, a greater abundance of Actinobacteria in group M. Detailed analysis identified that norank_o_Coriobacteriales and norank_f_norank_o_Coriobacteriales were the main genera that resulted in the relative abundance of Actinobacteria in group H. The abundance of norank_f_norank_o_Coriobacteriales was positively correlated with ADG and serum TNF-α concentration in this study. Previous studies have shown that Actinobacteria were correlated with lipid content that promotes animal growth and development. These resultant bioactive substances participate in host metabolism, but their abundances increase the prevalence of inflammation (Bredholdt et al., 2007; Dougal et al., 2014). The mechanism of action of Actinobacteria includes the following three aspects: Firstly, acetic acid produced by Actinobacteria, when absorbed by the colon, crosses the blood–brain barrier to activate hypothalamic neurons, activating the parasympathetic nerves, stimulating the secretion of growth hormone that increases feed intake and body weight gain (Koeth et al., 2013; Warrier et al., 2015); Secondly, acetic acid stimulates adipocyte storage by stimulating islet B cells to release insulin (Vrieze et al., 2012); Thirdly, hydrogen sulfide produced by Actinobacteria increases intestinal wall permeability and impairs intestinal integrity, thus some harmful bacteria and their products in the intestine can pass into the intestinal mucosa, leading to increased levels of inflammation and decreased levels of anti-inflammatory factors, triggering a chronic inflammatory state and insulin resistance (Hui et al., 2012; Kilian, 2018). In addition, this experiment also reported the abundance of Firmicutes in group H was numerically increased compared to group M and L, and especially the genus of [Eubacterium]_coprostanoligenes_group. The abundance of [Eubacterium]_coprostanoligenes_group was positively correlated with DM digestibility and serum MDA and TNF-α concentrations, but negatively correlated with CAT activity. Therefore, the greater abundance of [Eubacterium]_coprostanoligenes_group may be one of the factors involved in lowering the concentrations of antioxidant indicators and higher inflammatory factors in group H.
The phyla of Tenericutes cannot produce the anti-oxidant, indole propionic acid and the increase in its abundance reduces antioxidant ability and evokes inflammation in the chorioamniotic membrane (Doyle et al., 2017; Menni et al., 2019). Patescibacteria is a phylum associated with a number of viruses (Flint et al., 2012), and a decrease in its abundance enhances the antioxidant ability in elderly people (Zhang et al., 2021). The low antioxidant indicators and nutrient digestibility, as well as, the great levels of inflammatory factors in group L were strongly associated with the higher abundances of Tenericutes and Paterscibacteria in this study. Further analysis of the differences at the family and genus levels found that the families of norank_f_norank_o_Mollicutes_RF39 and Saccharimonadaceae and genera of norank_f_norank_o_Mollicutes_RF39 and Candidatus_Saccharimonas were significantly enriched in group L compared to group H and M, and were the main reason for the increased abundances of Tenericutes and Paterscibacteria. Moreover, the genus of norank_f_norank_o_Mollicutes_RF39 was positively correlated with the concentrations of IL-2 and IL-6, and the genus of Candidatus_Saccharimonas was negatively correlated with CAT, T-SOD, and GSH-Px activity.
Compared to group H and M, the phylum of Bacteroidetes in group L was numerically increased, and the family of Rikenellaceae and the genus of Rikenellaceae_RC9_gut_group were significantly enriched in group L. Previous studies have found that Rikenellaceae was highly correlated with the lean body mass (BMI < 25) in humans (Flint et al., 2012). Similar results were obtained in this study with the genus Rikenellaceae_RC9_gut_group being inversely correlated with CP digestibility, ADG, and serum CAT activity. The aforementioned finding suggests that the lower antioxidant ability and nutrient digestibility, as well as the greater inflammatory responses of donkey jennets on the L diet may be the result of increases in norank_f_norank_o_Mollicutes_RF39, Candidatus_Saccharimonas, and Rikenellaceae_RC9_gut_group.
Fibrobacteres is the main ligninolytic phylum and is closely correlated with antioxidant abilities in animals. Studies in piglets have shown a positive correlation between Fibrobacteres with both T-AOC and CAT activity in the jejunum (Xu, 2021; Xu et al., 2021). In this study the Fibrobacteres phylum, the family of Fibrobacteraceae and genus of Fibrobacter, were increased in the pregnant jennets on the M diet compared to group H and L. Furthermore, the abundance of Fibrobacter was positively correlated with serum CAT and T-SOD activity and inversely correlated with IL-2 concentration.
Additionally, phylum of Firmicutes in the group M also increased numerically compared to group H and L, and the families of Achiminococcaceae and Clostridiaceae_1 and genera of Phascolarctobacterium and Clostridium_sensustricto_1 were enriched. Phascolarctobacterium produces organic acids that can inhibit Escherichia coli and regulate the gastrointestinal health (Xu et al., 2021). The decrease in the abundance of Clostridium_sensustricto_1 is often a marker of intense inflammation, and a significant decrease in the abundance of Clostridium_sensustric-to_1 has been found in patients with immunoglobulin A nephropathy and membranous nephropathy (Dong et al., 2020). In this study, the abundance of Phascolarctobacterium was positively correlated with CAT activity in serum, and Clostridium_sensustricto_1 was positively correlated with CP, ADF and NDF digestibilities, as well as serum CAT and GSH-Px concentrations, but negatively correlated with the MDA concentration. The above results show that the great antioxidant status of pregnant jennets on the M diet was closely related to the abundances of these bacteria, with Fibrobacter playing a key role.
The results of this experiment demonstrates that dietary energy level can affect pregnant jennets’ health by altering the gut microbial community. Previous studies have found that energy intake indirectly changes the amount of NO produced by altering amino acid metabolism. Nitrites and nitrates oxidized by NO affect the composition and function of the microbiota (Alemany, 2013). It is thus proposed that dietary energy level may further affect the composition of intestinal microflora by influencing amino acid metabolism, and these specific mechanisms need to be further exploration.
In addition, this experiment found that as the parturition approached, mobilization of body fat in the body increased, aggravating the physiological function of the liver with a significant increases in NEFA and BHBA concentrations, especially within the 7 days before parturition. This exacerbated the oxidative stress and inflammatory responses of the pregnant jennets, as indicated by a decrease in antioxidant enzyme activity and an increase in inflammatory factors. Meanwhile, the relative abundance of Lachnospiraceae_NK4B4 increased and norank_f_Bacteroidales_UCG-001 decreased, apparently being associated with changes in the energy metabolism of the host. Therefore, it is critical to ensure that an appropriate level of dietary energy is available for pregnant donkeys during late gestation, and that the feeding and management within 7 days before the parturition deserves special attention.
5 Conclusion
In conclusion, dietary energy level of 10.49 MJ/kg for pregnant donkeys during late gestation increased the antioxidant capacity, reduces inflammatory cytokines, and promoted nutrients digestion, and fetal growth. These changes were associated with the alterations in the rectal microbiota composition. Superfluous or insufficient dietary energy supply during late gestation led to an increase of oxidative stress and inflammation in pre-partum donkeys.
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 studies were approved by all animal procedures were performed under the national standard Guidelines for Ethical Review of Animal Welfare (GB/T 35892-2018). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
YG: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Validation, Visualization, Writing – original draft. GY: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. FH: Data curation, Formal analysis, Investigation, Writing – original draft. XG: Writing – review & editing. BS: Supervision, Writing – review & editing. YZ: Writing – review & editing, SY: Conceptualization, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants of Program for Improving the Scientific Research Ability of Youth Teachers of Inner Mongolia Agricultural University (project no. BR220142), High-level Talent Start-up Research Project of Inner Mongolia Agricultural University (project no. NDYB2020-4), Basic Research Fund for Universities in Inner Mongolia Autonomous Region (project no. BR22-13-13).
Acknowledgments
For their help during laboratory and data analysis, the authors express deep appreciation to Weiqi Wu, Yuanxi Yue, and Yan Zhou from College of Animal Science, Inner Mongolia Agricultural University, China.
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/fmicb.2024.1308171/full#supplementary-material
References
Al Jassim, R. A., and Andrews, F. M. (2009). The bacterial community of the horse gastrointestinal tract and its relation to fermentative acidosis, laminitis, colic, and stomach ulcers. Vet. Clin. North Am. Equine Pract. 25, 199–215. doi: 10.1016/j.cveq.2009.04.005
Alemany, M. (2013). Relationship between energy dense diets and white adipose tissue inflammation in metabolic syndrome. Nutr. Res. 33, 1–11. doi: 10.1016/j.nutres.2012.11.013
Biddle, A. S., Tomb, J. F., and Fan, Z. (2018). Microbiome and blood analyte differences point to community and metabolic signatures in lean and obese horses. Front. Vet. Sci. 5:225. doi: 10.3389/fvets.2018.00225
Bredholdt, H., Galatenko, O. A., Engelhardt, K., Fjaervik, E., Terekhova, L. P., and Zotchev, S. B. (2007). Rare actinomycete bacteria from the shallow water sediments of the Trondheim fjord, Norway: isolation, diversity and biological activity. Environ. Microbiol. 9, 2756–2764. doi: 10.1111/j.1462-2920.2007.01387.x
Callaway, T. R., Dowd, S. E., Edrington, T. S., Anderson, R. C., Krueger, N., Bauer, N., et al. (2010). Evaluation of bacterial diversity in the rumen and feces of cattle fed different levels of dried distillers grains plus solubles using bacterial tag-encoded FLX amplicon pyrosequencing. J. Anim. Sci. 88, 3977–3983. doi: 10.2527/jas.2010-2900
Chen, H., Wang, C., Huasai, S., and Chen, A. (2022). Effect of prepartum dietary energy density on beef cow energy metabolites, and birth weight and antioxidative capabilities of neonatal calves. Sci. Rep. 12:4828. doi: 10.1038/s41598-022-08809-6
Creely, S. J., McTernan, P. G., Kusminski, C. M., Fm, F., Da Silva, N. F., Khanolkar, M., et al. (2007). Lipopolysaccharide activates an innate immune system response in human adipose tissue in obesity and type 2 diabetes. Am. J. Physiol. Endocrinol. Metab. 292, E740–E747. doi: 10.1152/ajpendo.00302.2006
Dalvi, P. S., Chalmers, J. A., Luo, V., Han, D. Y., Wellhauser, L., Liu, Y., et al. (2017). High fat induces acute and chronic inflammation in the hypothalamus: effect of high-fat diet, palmitate and TNF-α on appetite-regulating NPY neurons. Int. J. Obes. 41, 149–158. doi: 10.1038/ijo.2016.183
Dearing, M. D., and Kohl, K. D. (2017). Beyond fermentation: other important services provided to endothermic herbivores by their gut microbiota. Integr. Comp. Biol. 57, 723–731. doi: 10.1093/ICB/ICX020
Dong, R., Bai, M., Zhao, J., Wang, D., Ning, X., and Sun, S. (2020). A comparative study of the gut microbiota associated with immunoglobulin a nephropathy and membranous nephropathy. Front. Cell. Infect. Microbiol. 10:557368. doi: 10.3389/fcimb.2020.557368
Dougal, K., de la Fuente, G., Harris, P. A., Girdwood, S. E., Pinloche, E., Geor, R. J., et al. (2014). Characterisation of the faecal bacterial community in adult and elderly horses fed a high fibre, high oil or high starch diet using 454 pyrosequencing. PLoS One 9:e87424. doi: 10.1371/journal.pone.0087424
Doyle, R. M., Harris, K., Kamiza, S., Harjunmaa, U., Ashorn, U., Nkhoma, M., et al. (2017). Bacterial communities found in placental tissues are associated with severe chorioamnionitis and adverse birth outcomes. PLoS One 12:e0180167. doi: 10.1371/journal.pone.0180167
Du, X., Liu, G., Loor, J. J., Fang, Z., Bucktrout, R., Yang, Y., et al. (2018). Impaired hepatic autophagic activity in dairy cows with severe fatty liver is associated with inflammation and reduced liver function. J. Dairy Sci. 101, 11175–11185. doi: 10.3168/jds.2018-15120
Flint, H. J., Scott, K. P., Duncan, S. H., Louis, P., and Forano, E. (2012). Microbial degradation of complex carbohydrates in the gut. Gut Microbes 3, 289–306. doi: 10.4161/gmic.19897
Grimm, P., Philippeau, C., and Julliand, V. (2017). Faecal parameters as biomarkers of the equine hindgut microbial ecosystem under dietary change. Animal 11, 1136–1145. doi: 10.1017/S1751731116002779
Holman, D. B., and Chénier, M. R. (2014). Temporal changes and the effect of sub-therapeutic concentrations of antibiotics in the gut microbiota of swine. FEMS Microbiol. Ecol. 90, 599–608. doi: 10.1111/1574-6941.12419
Hui, H., Cai, D., Le, K., Li, Y., Xu, X., Liu, T., et al. (2012). Alterations in fecal lactobacillus and bifidobacterium species are associated with differences in cholesterol metabolism in type 2 diabetic patients. Endocr. Abstr. 29:558.
Kilian, M. (2018). The oral microbiome—friend or foe? Eur. J. Oral Sci. 126, 5–12. doi: 10.1111/eos.12527
Kim, H. B., and Isaacson, R. E. (2015). The pig gut microbial diversity: understanding the pig gut microbial ecology through the next generation high throughput sequencing. Vet. Microbiol. 177, 242–251. doi: 10.1016/j.vetmic.2015.03.014
Koeth, R. A., Wang, Z., Levison, B. S., Buffa, J. A., Org, E., Sheehy, B. T., et al. (2013). Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat. Med. 19, 576–585. doi: 10.1038/nm.3145
Liang, Y., Alharthi, A. S., Bucktrout, R., Elolimy, A. A., Lopreiato, V., Martinez-Cortés, I., et al. (2020). Body condition alters glutathione and nuclear factor erythroid 2-like 2 (NFE2L2)-related antioxidant network abundance in subcutaneous adipose tissue of periparturient Holstein cows. J. Dairy Sci. 103, 6439–6453. doi: 10.3168/jds.2019-17813
Liang, X. S., Yue, Y. X., Zhao, Y. L., Guo, Y. M., Guo, X. Y., Shi, B. L., et al. (2022). Effects of dietary concentrate to forage ratio on milk performance, milk amino acid composition and milk protein synthesis of lactating donkeys. Anim. Feed Sci. Tech. 292:115444. doi: 10.1016/j.anifeedsci.2022.115444
Malesza, I. J., Malesza, M., Walkowiak, J., Mussin, N., Walkowiak, D., Aringazina, R., et al. (2021). High-fat, western-style diet, systemic inflammation, and gut microbiota: a narrative review. Cell 10:3164. doi: 10.3390/cells10113164
Maury, E., and Brichard, S. M. (2010). Adipokine dysregulation, adipose tissue inflammation and metabolic syndrome. Mol. Cell. Endocrinol. 314, 1–16. doi: 10.1016/j.mce.2009.07.031
Meijer, K., de Vries, M., Al-Lahham, S., Bruinenberg, M., Weening, D., Dijkstra, M., et al. (2011). Human primary adipocytes exhibit immune cell function: adipocytes prime inflammation independent of macrophages. PLoS One 6:e17154. doi: 10.1371/journal.pone.0017154
Menni, C., Hernandez, M. M., Vital, M., Mohney, R. P., Spector, T. D., and Valdes, A. M. (2019). Circulating levels of the anti-oxidant indoleproprionic acid are associated with higher gut microbiome diversity. Gut Microbes 10, 688–695. doi: 10.1080/19490976.2019.1586038
Moradi, M., Chashnidel, Y., Teimouri Yansari, A., Khazari, B., and Mansouryar, M. (2021). Moderate overfeeding of different sources of metabolizable energy and protein. II: effects on inflammatory status of sheep in late pregnancy and growth trajectory of the offspring. Theriogenology 176, 115–121. doi: 10.1016/j.theriogenology.2021.09.021
Sordillo, L. M., and Aitken, S. L. (2009). Impact of oxidative stress on the health and immune function of dairy cattle. Vet. Immunol. Immunopathol. 128, 104–109. doi: 10.1016/j.vetimm.2008.10.305
Spence, C., Wells, W. G., and Smith, C. J. (2006). Characterization of the primary starch utilization operon in the obligate anaerobe Bacteroides fragilis: regulation by carbon source and oxygen. J. Bacteriol. 188, 4663–4672. doi: 10.1128/JB.00125-06
Stanley, D., Geier, M. S., Chen, H., Hughes, R. J., and Moore, R. J. (2015). Comparison of fecal and cecal microbiotas reveals qualitative similarities but quantitative differences. BMC Microbiol. 15:51. doi: 10.1186/s12866-015-0388-6
Tapias, V., Hu, X., Luk, K. C., Sanders, L. H., Lee, V. M., and Greenamyre, J. T. (2017). Synthetic alpha-synuclein fibrils cause mitochondrial impairment and selective dopamine neurodegeneration in part via iNOS-mediated nitric oxide production. Cell. Mol. Life Sci. 74, 2851–2874. doi: 10.1007/s00018-017-2541-x
Vrieze, A., Van Nood, E., Holleman, F., Salojärvi, J., Kootte, R. S., Bartelsman, J. F., et al. (2012). Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome. Gastroenterol 143, 913–916.e7. doi: 10.1053/j.gastro.2012.06.031
Warrier, M., Shih, D. M., Burrows, A. C., Ferguson, D., Gromovsky, A. D., Brown, A. L., et al. (2015). The TMAO-generating enzyme flavin monooxygenase 3 is a central regulator of cholesterol balance. Cell Rep. 10, 326–338. doi: 10.1016/j.celrep.2014.12.036
Weber, C., Hametner, C., Tuchscherer, A., Losand, B., Kanitz, E., Otten, W., et al. (2013). Hepatic gene expression involved in glucose and lipid metabolism in transition cows: effects of fat mobilization during early lactation in relation to milk performance and metabolic changes. J. Dairy Sci. 96, 5670–5681. doi: 10.3168/jds.2012-6277
White, H. M. (2015). The role of TCA cycle anaplerosis in ketosis and fatty liver in periparturient dairy cows. Animals (Basel) 5, 793–802. doi: 10.3390/ani5030384
Xu, J. (2021) Study on the effect of microbial oligosaccharides on gut microbial metabolism. PhD Thesis. Wuxi: Jiangnan University. (In Chinese).
Xu, B., Qin, W., Xu, Y., Yang, W., Chen, Y., Huang, J., et al. (2021). Dietary quercetin supplementation attenuates diarrhea and intestinal damage by regulating gut microbiota in weanling piglets. Oxidative Med. Cell. Longev. 2021:6221012. doi: 10.1155/2021/6221012
Yin, G., Wu, W., Liang, X., Zhao, Y., Guo, X., Zhang, Y., et al. (2021). Effect of dietary energy level on body weight change and nutrient digestibility of female donkeys in late pregnancy. Feed Res. 44, 79–82. doi: 10.13557/j.cnki.issn1002-2813.2021.01.022
Zhang, J. Y., Lo, H. C., Yang, F. L., Liu, Y. F., Wu, W. M., and Chou, C. C. (2021). Plant-based, antioxidant-rich snacks elevate plasma antioxidant ability and alter gut bacterial composition in older adults. Nutrients 13:3872. doi: 10.3390/nu13113872
Keywords: gestation, dietary energy, donkey jennet, gut microbiota, oxidative stress, immunity
Citation: Guo Y, Yin G, Hui F, Guo X, Shi B, Zhao Y and Yan S (2024) Effects of dietary energy level on antioxidant capability, immune function and rectal microbiota in late gestation donkeys. Front. Microbiol. 15:1308171. doi: 10.3389/fmicb.2024.1308171
Edited by:
Guolong Zhang, Oklahoma State University, United StatesReviewed by:
Tarique Hussain, Nuclear Institute for Agriculture and Biology, PakistanXiaoling Zhou, Tarim University, China
Qiugang Ma, China Agricultural University, China
Copyright © 2024 Guo, Yin, Hui, Guo, Shi, Zhao and Yan. 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: Sumei Yan, eWFuc21pbWF1QDE2My5jb20=