- 1College of Animal Science, Guizhou University, Guiyang, China
- 2Key Laboratory of Animal Diseases and Veterinary Public Health of Guizhou Province, College of Animal Science, Guizhou University, Guiyang, China
Fermented distiller's grains (FDG)-based diets are nutritious and can improve the growth and intestinal immunity in livestock. However, there is limited research examining the effect of feeding FDG-based diets on changes in intestinal metabolites and related pathways in livestock. In this study, nine Guanling crossbred cattle (Guizhou Guanling Yellow cattle × Simmental cattle) were selected and randomly divided into a basal diet (BD) group and two experimental groups fed with FDG replacing 15% and 30% of the daily ration concentrates (FDG-Case A and FDG-Case B), respectively, with three cattle in each group. Fresh jejunum (J) and cecum (C) tissues were collected for metabolomic analysis. Differential metabolites and metabolic pathways were explored by means of univariate and multivariate statistical analysis. Compared with the J-BD group, 30 and 100 differential metabolites (VIP > 1, p < 0.05) were obtained in the J-FDG-Case A group and J-FDG-Case B group, respectively, and the J-FDG-Case B vs. J-FDG-Case A comparison revealed 63 significantly differential metabolites, which were mainly divided into superclasses including lipids and lipid-like molecules, organoheterocyclic compounds, and organic acids and derivatives. Compared with the C-BD, 3 and 26 differential metabolites (VIP > 1, p < 0.05) were found in the C-FDG-Case A group and C-FDG-Case B group, respectively, and the C-FDG-Case B vs. C-FDG-Case A comparison revealed 21 significantly different metabolites, which were also mainly divided into superclasses including lipids and lipid-like molecules, organoheterocyclic compounds, and organic acids and derivatives. A total of 40 metabolic pathways were identified, with a significance threshold set at p < 0.05. Among them, 2, 14, and 18 metabolic pathways were significantly enriched in the J-FDG-Case A vs. J-BD, J-FDG-Case B vs. J-BD, and J-FDG-Case B vs. J-FDG-Case A comparisons, respectively. Meanwhile, 1, 2, and 3 metabolic pathways were obtained in the C-FDG-Case A vs. C-BD, C-FDG-Case B vs. C-BD, and C-FDG-Case B vs. C-FDG-Case A comparisons, respectively. Furthermore, four significant metabolic pathways, namely insulin resistance, biosynthesis of unsaturated fatty acids, linoleic acid metabolism, and primary bile acid biosynthesis, were significantly enriched in Guanling crossbred cattle fed FDG diets. These results suggest that feeding FDG diets may promote the growth and intestinal immunity of Guanling crossbred cattle by regulating metabolic patterns of lipid compounds and related metabolic pathways. This study sheds light on the potential metabolic regulatory mechanisms of FDG diets and offers some references for their use in livestock feed.
1. Introduction
With the rapid advancement of China's agricultural industry in recent years, there has been a surge in costs for essential raw materials such as corn and bran, which in turn has resulted in elevated animal breeding expenses (1). Distillers' grains (DGs) are a plentiful, readily available, and cost-effective alternative resource in China, with a production of 20–30 million tons in 2020. DGs are worth reusing due to their high content of crude protein, fiber, fat, and vitamins, making them an excellent alternative resource for animal feed (2). The growing use of fermented distillers' grains (FDGs) has been evidenced in the livestock industry, including cattle, pigs, and broilers, which can improve the production performance and meat quality of the animals (3–5). Evidence has shown that the supplement of fermented distillers' grains (FDGs) to cattle diet can increase weight gain, improve slaughter performance, including crude protein content, and reduce the feed-to-weight ratio in cattle (6, 7). Additionally, supplementing beef cattle diets with FDGs can increase the protein content of the resulting beef (8). These studies suggest a promising potential for the utilization of FDGs in animal feeds, which promotes animal growth and development and improves beef quality. A recent study demonstrated that feeding Moutai FDG enhanced meat quality by increasing back muscle thickness, intramuscular protein content, and fat content, as well as upregulating the expression of various lipid metabolism molecules in crossbred cattle, suggesting a positive regulation of the Moutai FDG diet on slaughter performance and muscle lipid metabolism (9).
The intestine is an organ with a special tubular structure that is primarily responsible for digesting food and absorbing water, inorganic salts, and nutrients (10). Physiological and pathological events that trigger disturbances of intestinal homeostasis can negatively affect intestinal health and cause intestinal diseases that affect organism growth and development (11). The intestine also serves as the largest immune organ and the frontline for immune defense in the body, whose health is related to the occurrence of intestinal diseases and the multi-system regulation of the entire body (12, 13). The maintenance of gut balance is associated with a variety of factors, including the immune system, the microbial ecosystem, and the diet. Additionally, the gut plays an important role in the regulation of systemic physiology, metabolism, and immunity (14). Evidence has shown that fermented diets can produce short-chain fatty acids (SCFAs), phenolics, and other small metabolic molecules that not only contribute to maintaining the balance between the animal's intestinal barrier function and microbial community but also provide energy to the intestinal epithelial cells, and such substances also have antioxidant properties, which reduce oxidative stress and inflammatory responses in the intestine (15–20). To date, most of the studies mentioned above have focused on the regulation of intestinal flora structure and metabolism in intestinal contents, with little attention to intestinal (wall) metabolism. The change of indices in the intestinal contents reflects to some extent the combined action of the body and the intestinal microorganisms, while the metabolic profile of the intestinal wall is more biased toward indicating changes in the body's metabolism. In this study, we investigated the effect of feeding FDG diets on the intestinal metabolic profile in Guanling crossbred cattle via LC-MS-based non-targeted metabolomics. The results will provide new insights into the regulatory effects of dietary FDG supplementation on intestinal metabolites and their metabolic pathways in ruminants, which will be beneficial for further development and utilization of FDGs as a potentially valuable feed resource.
2. Materials and methods
2.1. Origin and treatment of FDGs
The DGs used in this study were produced by Kweichow Moutai Group in Moutai Town, Renhuai City, Guizhou Province, China. The main ingredients of the fermentation substrate include 92% fresh DG, 3% corn flour, 3% rapeseed meal, and 2% wheat bran. After the fermentation substrate was evenly mixed with a feed mixer, 50 kg of microbial starter and water mixture were added, which was prepared by mixing 100 g of biological fermentation agent with 5 kg of water containing 0.1% brown sugar. The mixture is added while stirring until it is completely added. After activation for 3–5 days, the DGs were bagged and sealed and then fermented for 9 days to obtain the FDGs (the actual cost of FDGs in either the 15% or 30% FDG group is approximately 0.4 RMB/kg) to be used for subsequent experiments. The components of the microbial starter included lactic acid bacteria, yeast, Bacillus, Bifidobacterium, Clostridium butyricum, amylase, protease, cellulase, and lipase.
2.2. Animal experiment
This study was conducted at the Yueyuan cattle farm located in Guanling County, Anshun City, Guizhou Province (105° 58′ E, 25° 98′ N, altitude, approximately 1863 m). The experimental period was from September 2020 to November 2020, during which the average minimum outdoor temperature was 16.9°C and the maximum temperature was 25.6°C. Nine Guanling crossbred cattle (Guizhou Guanling Yellow cattle × Simmental cattle) in good health and of similar age (18 months old), weighing 300 ± 25 kg, were provided by the beef cattle fattening farm of Guanling County Yellow Cattle Group. Laboratory tests showed that these cattle were negative for infection with Brucella, Mycobacterium tuberculosis, foot-and-mouth disease virus, and lumpy skin disease virus. These cattle were ear-marked and randomly divided into a basal diet (BD) group, a Case A group, and a Case B group (containing 0%, 15%, and 30% FDG, respectively, to replace part of the concentrates), with three cattle in each group. The composition and nutrient content of the diets are shown in Supplementary Table S1 (21). The diets were thoroughly mixed to form a complete mixed diet and fed twice daily (9:00 a.m. and 4:30 p.m.). Cattle in each group were housed in separate pens and fed free water under identical housing conditions throughout the duration of the experiment.
2.3. Collection and processing of samples
On the 60th day of the formal feeding period, the animals underwent overnight fasting, and then tissues were collected from the jejunum, namely jejunum-BD (J-BD), jejunum-FDG-Case A (J-FDG-Case A), and jejunum-FDG-Case B (J-FDG-Case B), and from the cecum, namely cecum-BD (C-BD), cecum-FDG-Case A (C-FDG-Case A), and cecum-FDG-Case B (C-FDG-Case B). The collected tissues were then rinsed with RNase-free saline. QC samples were obtained by mixing equal volumes of the extracts from all collected samples. Prior to the metabolomics analysis, 30 mg of intestinal tissue samples were added to a 1.5-mL tube with 20 μL of L-2-chlorophenylalanine (0.3 mg/mL) dissolved in methanol as internal standard, and then 400 μL mixture of methanol/water (1/4, v/v) was added and pre-chilled at −20°C for 2 min and then ground in a grinder for 2 min. All the samples were subjected to ultrasonic extraction for 10 min, stored at −20°C for 30 min, and then centrifuged at 13,000 × g for 10 min at 4°C. A total of 300 μL of supernatants were dried in a freeze-concentration centrifugal dryer and resolubilized by 200 μL of methanol/water (1/4), vortexed for 30 s, and subjected to ultrasonic extraction for 3 min. After placement for 2 h at −20°C, the samples were centrifuged at 4°C (13,000 × g) for 10 min, and 150 μL of supernatants were filtered through 0.22-μm microfilters and then transferred to liquid chromatography (LC) vials (22).
2.4. LC-MS analysis
The LC-MS analysis of intestinal samples in both ESI positive and ESI negative ion modes was conducted using an ACQUITY UPLC I Class system (Waters Corporation Milford, USA), coupled with a VION IMS Q-TOF mass spectrometer (Waters Corporation Milford, USA). A 1-μL aliquot of intestinal tissue samples was injected into an ACQUITY UPLC™ BEH C18 column (50 × 2.1 mm i.d., 1.7 μm; Waters Corporation, Milford, MA, USA). Mobile phase A consisted of 0.1% formic acid in water, and mobile phase B comprised a mixture of acetonitrile/methanol (2/3) containing 0.1% formic acid (v/v). The flow rate was 0.4 mL/min. The conditions for UPLC separation and ESI-VIONIMS Q-TOF detection are shown in Supplementary Tables S2, S3. The QC samples were utilized to evaluate the reproducibility and reliability of the LC-MS system, which contains metabolic information that enables the assessment of the stability of the mass spectrometry system. The ion peaks with a relative standard deviation of >0.4 for the QC group samples were removed. The stability of the system was then assessed using principal component analysis (PCA). The PCA model plot was obtained after seven cycles of cross-validation to verify whether the QC samples were closely clustered together, thus judging the stability of the instrument's detection.
2.5. Data analysis
For metabolomics data processing, the obtained LC-MS raw data were initially analyzed using the Progenesis QI v2.3 software (Non-linear Dynamics, Newcastle, UK). The metabolites were qualitatively analyzed by utilizing the Human Metabolome Database (HMDB), Lipid Maps (v2.3), the METLIN database, and a self-developed database. For the extracted data, ion peaks with missing values (0 value)>50% were eliminated, the 0 value with half of the minimum value was replaced, and the qualitatively obtained compounds were screened according to the qualitative result evaluation (score). The qualitative results were screened between 36 and 60 points, below which they were considered inaccurate and deleted. Finally, the positive and negative ion data were combined into a data matrix containing all the information extracted from the original data necessary for analysis. Subsequently, multivariate statistical analyses were carried out using SIMCA 14.0 (Umetrics AB, Umeå, Sweden).
The PCA and orthogonal partial least squares discriminant analysis (OPLS-DA) were performed to reveal the global metabolic changes between the BD and FDG groups using SIMCA 14.0. An unsupervised PCA was first used to examine both the metabolic differences between the groups and the individual metabolic differences within samples. To maximize the differences within the model, a supervised OPLS-DA was performed on metabolites between groups. In addition, to prevent overfitting of the model, the quality of the model was examined by response permutation testing with 200 responses. Then, univariate analyses were performed on the sample data, which consisted of a t-test and fold change analysis. The fold change (FC) value was calculated through a fold analysis of variance to assess differences in metabolite expression between the two groups. The p-values were obtained using a t-test, and then volcano plotting was performed using Origin 2020 (Origin Lab, Northampton, MA, USA) to visualize the p-values and FC to facilitate the rapid identification of metabolites with significant differences. Multivariate and univariate statistical analyses were combined to screen the differential metabolites between groups.
2.6. Key biomarkers and metabolic pathway analysis
To reveal the mechanism of metabolic pathway variation in intestinal samples, an analysis of differential metabolites was carried out with metabolic pathway enrichment analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.kegg.jp/kegg/pathway.html). Their KEGG IDs and pathways were obtained, and then the number of metabolites enriched in each corresponding pathway was calculated. The pathway with a p < 0.05 was selected as an enriched pathway. The formula for calculating p-values is provided as follows:
where N: the total number of metabolites, n: the number of differential metabolites, M: the number of metabolites annotated as a specific pathway, and m: the number of differential metabolites annotated as a specific pathway. A lower p-value indicates a more significant difference within the metabolic pathway.
2.7. Statistical analyses
All continuous variables were expressed as the mean ± standard deviation (M ± SD). The differences between the groups were analyzed using Student's t-test, where the significant difference was set at a p < 0.05. All statistical analyses were performed using SPSS Statistics 26.0 (Chicago, IL, USA).
3. Results
3.1. Quality control of the metabolomics of jejunum and cecum tissues
An LC-MS full-scan detection method was used to illustrate the metabolic alterations in the jejunum and cecum tissues of Guanling crossbred cattle fed the FDG diets. In total, 8,464 features in ESI+ mode and 2,564 ion features in ESI− mode were detectable, and 1,399 metabolites were identified. PCA scoring plots were initially generated using the processed data obtained from the metabolomics analysis of the jejunum and cecum to evaluate systematic errors in the samples and the trends within each comparison. As shown in Supplementary Figure S1, the quality control (QC) samples clustered in the center of the PCA score plots suggest that the analyses were robust and reproducible. Clear separation trends were observed in the FDG-Case A vs. BD, FDG-Case B vs. BD, and FDG-Case B vs. FDG-Case A comparisons, respectively (Figures 1A–F). To maximize the identification of differential metabolites between BD and FDG groups, OPLS-DA was used to observe the tendency for separation on the metabolic spectrum (Figures 2A–F). The validity of the model was evaluated through permutation testing (n = 200), with R2 values close to 1 and intercepts between the regression line and the vertical axis of Q2 <0. Moreover, Q2 gradually decreased with reduced displacement retention, affirming the absence of overfitting in the original model (Figures 2G, H). These results provide assurance for subsequent tests and analyses.
Figure 1. PCA model score of the untargeted metabolomics approach. (A) J-FDG-Case A vs. J-BD comparison (R2X = 0.786), (B) J-FDG-Case B vs. J-BD comparison (R2X = 0.781), (C) J-FDG-Case B vs. J-FDG-Case A comparison (R2X = 0.810), (D) C-FDG-Case A vs. C-BD comparison (R2X = 0.828), (E) C-FDG-Case B vs. C-BD comparison (R2X = 0.835), and (F) C-FDG-B vs. C-FDG-Case A comparison (R2X = 0.818; p < 0.05). Each dot, triangle, or diamond on the plot represents a sample in the corresponding group. BD, basal diet group; FDG-Case A and FDG-Case B: 15% FDG group and 30% FDG group.
Figure 2. OPLS-DA for data distribution. 2D score scatter plot of the PCA model for the (A) J-FDG-Case A vs. BD (R2X = 0.786), (B) J-FDG-Case B vs. BD (R2X = 0.781), (C) J-FDG-Case B vs. J-FDG-Case A (R2X = 0.810), (D) C-FDG-Case A vs. BD (R2X = 0.828), (E) C-FDG-Case B vs. BD (R2X = 0.835), and (F) C-FDG-B vs. C-FDG-Case A comparisons, respectively (R2X = 0.818; p < 0.05). Permutation plots for the OPLS-DA model show R2 (blue) and Q2 (green) values. The results of the permutation test strongly indicate that the original model was valid [(G, H), R2 intercept = 0.888/0.644, Q2 intercept = −0.805/−0.385]. BD, basal diet group; FDG-Case A and FDG-Case B: 15% FDG group and 30% FDG group.
3.2. Differential metabolites in the intestinal tissues of FDG and BD groups
To visualize the fold change (FC) analysis of 1,399 metabolites for quick identification of significantly differential metabolites, a volcano plot was created by transferring the FC values of each metabolite to log2(FC) and by transferring the p-value (p = 0.05) of Student's t-test to –log (p) while simultaneously meeting the variable importance in projection (VIP) values of > 1 for the first principal component. As shown in Figures 3A–F, a total of 178 metabolites [78 metabolites upregulated with log2(FC) > 0 and 100 metabolites downregulated with log2(FC) <0], 227 metabolites (121 upregulated and 106 downregulated), and 179 metabolites (174 upregulated and 5 downregulated) were identified in the J-FDG-Case A vs. J-BD, J-FDG-Case B vs. J-BD, and J-FDG-Case B vs. J-FDG-Case A comparisons, respectively, meeting the criteria of p < 0.05 and VIP value of >1. In addition, 12 metabolites (10 upregulated and 2 downregulated), 58 metabolites (57 upregulated and 1 downregulated), and 55 metabolites (43 upregulated and 12 downregulated) were identified in the C-FDG-Case A vs. C-BD, C-FDG-Case B vs. C-BD, and C-FDG-Case B vs. C-FDG-Case A comparisons, respectively.
Figure 3. Volcano plots based on LC-MS of the jejunal and cecal metabolites. (A–F) Volcano plot based on LC-MS analyses of J-FDG-Case A vs. J-BD (A), J-FDG-Case B vs. J-BD (B), J-FDG-Case B vs. J-FDG-Case A (C), C-FDG-Case A vs. C-BD (D), C-FDG-Case B vs. C-BD (E), and C-FDG-Case B vs. C-FDG-Case A (F) comparisons, respectively. The 1,399 significantly altered metabolites in the model group. Red and blue represent upregulated and downregulated metabolites, respectively, with VIP > 1, fold change> 1.3, and a p < 0.05. The gray area indicates unchanged metabolites with a fold change of < 1.3 and a p > 0.05. BD, basal diet group; FDG-Case A and FDG-Case B: 15% FDG group and 30% FDG group.
The levels of the top 50 differential metabolites were visualized by a heatmap depicted in Figures 4A–F, in which colors represent increased (red) or decreased (blue) abundance, with the intensity reflecting the corresponding concentration. Of these, the relative richness (referring to the ratio of FDG/BD) of 23, 49, and 50 metabolites were decreased in the J-FDG-Case A vs. J-BD, J-FDG-Case B vs. J-BD, and J-FDG-Case B vs. J-FDG-Case A comparisons, respectively. Conversely, the richness of 7, 1, and 0 metabolites was significantly increased in the aforementioned jejunum comparisons, respectively. As for cecum samples, the relative richness of 3, 26, and 21 metabolites was decreased in the C-FDG-Case A vs. C-BD, C-FDG-Case B vs. C-BD, and C-FDG-Case B vs. C-FDG-Case A comparisons, respectively. However, no upregulated differential metabolites were found in the aforementioned cecum comparisons, respectively.
Figure 4. Heatmap of the top 50 differential metabolites in the jejunum and cecum. Heatmap plot based on LC-MS analyses of J-FDG-Case A vs. J-BD (A), J-FDG-Case B vs. J-BD (B), J-FDG-Case B vs. J-FDG-Case A (C), C-FDG-Case A vs. C-BD (D), C-FDG-Case B vs. C-BD (E), and C-FDG-Case B vs. C-FDG-Case A (F) comparisons, respectively. The graph displays the differential metabolites on the Y-axis and the sample names on the X-axis. The color gradient from blue to red represents the expression abundance of the metabolites, with red indicating higher expression abundance and blue indicating lower expression abundance. BD, basal diet group; FDG-Case A and FDG-Case B: 15% FDG group and 30% FDG group.
The classification of all significantly differential metabolites is shown in Figures 5A, B, and Tables 1, 2 offer detailed information on the top 20 differential metabolites. Of these, in the J-FDG-Case A vs. J-BD comparison, they were mainly classified into lipids and lipid-like molecules (fatty acyls, steroids and steroid derivatives, and prenol lipids) and organoheterocyclic compounds (imidazopyrimidines). In the J-FDG-Case B vs. J-BD comparison, the differential metabolites were mainly classified into lipids and lipid-like molecules (fatty acyls, steroids and steroid derivatives, and glycerophospholipids), organic acids and derivatives (carboxylic acids and derivatives), and organic oxygen compounds (organooxygen compounds). In the J-FDG-Case B vs. J-FDG-Case A comparison, the differential metabolites were mainly classified into lipids and lipid-like molecules (glycerophospholipids, fatty acyls, and prenol lipids), organic acids and derivatives (carboxylic acids and derivatives), and organoheterocyclic compounds (tetrahydroisoquinolines). As for cecum, in the C-FDG-Case A vs. C-BD comparison, the differential metabolites were mainly classified into lipids and lipid-like molecules (steroids and steroid derivatives, and prenol lipids) and organoheterocyclic compounds (diazines). In the C-FDG-Case B vs. C-BD, the differential metabolites were mainly classified into lipids and lipid-like molecules (fatty acyls, steroids and steroid derivatives, and glycerophospholipids), organoheterocyclic compounds (imidazopyrimidines and diazines), and benzenoids (phenols, benzene, and substituted derivatives). In C-FDG-Case B vs. C-FDG-Case A, the differential metabolites were mainly classified into lipids and lipid-like molecules (fatty acyls, glycerophospholipids, and steroids and steroid derivatives), organoheterocyclic compounds (imidazopyrimidines), and organic acids and derivatives (carboxylic acids and derivatives). Collectively, these results indicate that the metabolite profiles in the jejunum and cecum are predominantly affected by the 15% FDG and 30% FDG diets.
Figure 5. Pie chart shows the classification of total significant metabolites. (A) Classify total metabolites into superclasses. (B) The total metabolites are further divided into classes.
Table 1. Top 20 differential metabolites in the jejunum of experimental Guanling crossbred cattle determined by ultra-performance liquid chromatography-mass spectrometry (UPLC-MS)-based non-targeted metabonomics.
Table 2. Top 20 differential metabolites in the cecum of experimental Guanling crossbred cattle determined by ultra-performance liquid chromatography-mass spectrometry (UPLC-MS)-based non-targeted metabonomics.
3.3. Metabolic pathway analysis
Next, we investigated potential metabolic pathways responsible for the observed changes in the intestinal metabolic profile associated with FDG diets. The identified metabolites were annotated with KEGG and HMDB. Results showed that a total of 40 metabolic pathways were identified (p < 0.05). Among them, there were 2 (insulin resistance and bile secretion), 14 (including biosynthesis of unsaturated fatty acids and terpenoid backbone biosynthesis), and 18 (including biosynthesis of unsaturated fatty acids, linoleic acid metabolism, and GnRH signaling pathway) significant metabolic pathways (p < 0.05) enriched in the J-FDG-Case A vs. J-BD, J-FDG-Case B vs. J-BD, and J-FDG-Case B vs. J-FDG-Case A comparisons, respectively, and 1 (primary bile acid biosynthesis), 2 (primary bile acid biosynthesis and adrenergic signaling in cardiomyocytes), and 3 (linoleic acid metabolism, biosynthesis of unsaturated fatty acids, and choline metabolism in cancer) significantly enriched metabolic pathways (p < 0.05) were obtained in the C-FDG-Case A vs. C-BD, C-FDG-Case B vs. C-BD, and C-FDG-Case B vs. C-FDG-Case A comparisons, respectively (Figures 6A–F). Furthermore, four significant pathways, including biosynthesis of unsaturated fatty acids, linoleic acid metabolism, insulin resistance, and primary bile acid biosynthesis [p < 0.05, –log(p)>1.3] were significantly enriched in Guanling crossbred cattle fed with FDG diets (Table 3). Among them, the insulin resistance pathway was found to be specific to the jejunum, while the primary bile acid biosynthesis pathway was specific to the cecum. The biosynthesis of unsaturated fatty acids and linoleic acid metabolism pathways were shared between the jejunum and cecum.
Figure 6. KEGG pathway analysis. Differentially altered metabolic pathways enriched in the J-FDG-Case A vs. J-BD (A), J-FDG-Case B vs. J-BD (B), J-FDG-Case B vs. J-FDG-Case A (C), C-FDG-Case A vs. C-BD (D), C-FDG-Case B vs. C-BD (E), and C-FDG-Case B vs. C-FDG-Case A comparisons (F), respectively. The graph displays the metabolic pathways on the Y-axis and the enrichment factor (Rich factor = The number of significantly different metabolites/The total number of metabolites in the pathway) on the X-axis. A higher Rich factor indicates a higher degree of enrichment, calculated by the number of significantly different metabolites divided by the total number of metabolites in the pathway. The color gradient from green to red represents decreasing p-values. Larger dots indicate more metabolites enriched in the pathway. BD, basal diet group; FDG-Case A and FDG-Case B: 15% FDG group and 30% FDG group.
Table 3. List of metabolic pathways with significant difference in the FDG group compared to the BD group.
4. Discussion
Guizhou province of China is well-known for the production of Jiang-flavor Moutai liquor, and a large number of distillers' grain (DG) resources are generated as byproducts with an annual output exceeding 19 million tons (23). DGs are nutrition-rich byproducts and are abundant in various crude substances (24, 25). Evidence has revealed that incorporating fermented distillers' grains (FDGs) as a feed additive is beneficial for the growth and development of ruminants, such as cattle and sheep. This is because fermentation treatment of DG increases feed palatability and intramuscular protein and fat content and stimulates rumen fermentation and intestinal development in animals (9, 26, 27). Additionally, the fermentation of DG improves its nutrients, thus promoting their utilization in feed production. Alterations in animal intestinal flora structure have been identified in several FDG diet-related studies (28–31). However, few studies have focused on the metabolic alterations of animal intestinal tissues to date. The present study showed that the intestinal metabolomic patterns of Guanling crossbred cattle fed the FDG diets were significantly different from those of BD diet-feeding cattle. FDG diets exert growth-promoting and anti-inflammatory potential through the regulation of fatty acyls, steroids, and steroid derivatives. The levels of several metabolites, including l-acetylcarnitine, 8,11,14-eicosatrienoic acid, adrenic acid, arachidonic acid, clupanodonic acid, gamma-Linolenic acid, docosahexaenoic acid, 9-OxoODE, chenodeoxycholic acid, and chenodeoxycholic acid glycine conjugate, were notably increased in the FDG group. Based on the identified differential metabolites, insulin resistance, biosynthesis of unsaturated fatty acids, linoleic acid metabolism, and primary bile acid biosynthesis metabolic pathways were significantly enriched in the FDG group, where biosynthesis of unsaturated fatty acids and linoleic acid metabolism are shared by the jejunum and cecum samples, suggesting potentially specific functions in the improvement of intestinal inflammation in cattle.
Insulin resistance is a specific metabolic pathway significantly enriched in the jejunum of the FDG-Case A group. The J-FDG-Case A group exhibits substantial upregulation of L-acetylcarnitine, which is an amino acid derivative capable of regulating insulin resistance and impacting gut health (Figure 6A). As a key site for food digestion and nutrient absorption, the intestinal tract also plays a vital role in influencing insulin sensitivity and metabolic function and can simultaneously be regulated by insulin (32). L-acetylcarnitine is an intermediate metabolite in insulin resistance, which increases fatty acid oxidation during lipid metabolism and promotes fatty acid utilization when the body's response to insulin is reduced. This results in elevated blood glucose and insulin overproduction utilization and reduced fat deposition and adipose tissue inflammation, thus improving insulin resistance and reducing oxidative stress and inflammatory responses caused by insulin resistance, suggesting the antioxidant and anti-inflammatory properties of L-acetylcarnitine (33). Our study suggests that FDG diets may help maintain intestinal epithelial cell integrity and function and sustain intestinal health by improving insulin resistance through L-acetylcarnitine, thereby reducing oxidative stress and promoting the growth and development of Guanling crossbred cattle. The complex interplay between insulin resistance and intestinal metabolism entails the regulation of multiple factors, and the underlying mechanisms remain to be further studied to gain new insights into the onset and progression of insulin resistance and devise new intervention strategies to improve insulin sensitivity and reduce the risk of metabolic disorders.
Primary bile acid biosynthesis is a metabolic pathway specific to the cecum, wherein chenodeoxycholic acid (CDCA) and chenodeoxycholic acid glycine conjugate (CDCA-Gly) are the main metabolites. The levels of CDCA and CDCA-Gly were significantly higher in both the C-FDG-Case A vs. C-BD and the C-FDG-Case B vs. C-BD comparisons (Figures 6D, E). CDCA and CDCA-Gly are produced through primary bile acid biosynthesis and play a critical role in the absorption and digestion of fat in the gallbladder and small intestine (34). Higher concentrations of CDCA and CDCA-Gly have been evidenced to promote gallbladder contraction, enhance the formation of bile and fatty acids to form an emulsion in the intestinal flow, which in turn improves the intestinal absorption of fatty acids, glycerol, and cholesterol, and influence the growth and composition of intestinal microorganisms, which ultimately promotes intestinal digestion and absorption (35, 36). The present study suggests that FDG diets may increase the concentrations of CDCA and CDCA-Gly, thereby stimulating the intestinal bile circulation system to reduce bile acid synthesis. Such promotion of intestinal digestion and absorption is beneficial for the intestinal metabolism of Guanling crossbred cattle. However, it is important to note that over-regulation may adversely affect the synthesis and circulation of bile acids. Therefore, it is necessary to control the metabolic concentrations of CDCA and CDCA-Gly.
The biosynthesis of unsaturated fatty acids and linoleic acid metabolism pathways are shared by the jejunum and cecum, and both of them belong to fatty acid metabolism. 8,11,14-Eicosatrienoic acid, adrenic acid, arachidonic acid, clupanodonic acid, gamma-linolenic acid, docosahexaenoic acid, gamma-linolenic acid, and 9-OxoODE are the initial substrates for the biosynthesis of unsaturated fatty acids and linoleic acid metabolism. These metabolites were significantly upregulated in J-FDG-Case B vs. J-BD, J-FDG-Case B vs. J-FDG-Case A, and C-FDG-Case B vs. C-FDG-Case A comparisons, respectively. In the biosynthesis of unsaturated fatty acids, 8,11,14-eicosatrienoic acid, clupanodonic acid, and docosahexaenoic acid (DHA) are precursors of ω-3 series fatty acids, which are involved in the synthesis of ω-3 polyunsaturated fatty acids and the regulation of nerve and brain development, immune function, and inflammatory response (37). Arachidonic acid (AA), gamma-linolenic acid, and adrenic acid are the precursors of ω-6 fatty acids and are involved in the synthesis of ω-6 polyunsaturated fatty acids, which are involved in the synthesis of cell membranes, cell signaling, and important physiologically active substances such as prostaglandins and leukotrienes possessing anti-inflammatory and immunomodulatory effects (38–43). Overall, these metabolites work together in the biosynthesis process to create a variety of unsaturated fatty acids that play important physiological functions and regulatory roles. During linoleic acid metabolism, 8,11,14-eicosatrienoic acid acts as a precursor of AA and forms interrupted aromatic terminal ketone (EETA) catalyzed by cytochrome P450 (44). 9-OxoODE is formed by cytochrome P450 enzymes and cyclooxygenases, which can inhibit EET hydrolases, leading to an increased concentration of EETs. This has an impact on the angiotensin pathway, inflammatory processes, and apoptosis. AA and γ-linolenic acid are important metabolites of linoleic acid, both of which are converted into prostaglandins (PGs) and leukotrienes through the activity of cyclo-oxidases and lipoxygenases, respectively (45–47). The aforementioned studies demonstrate that γ-linolenic acid, AA, and DHA possess anti-inflammatory properties. γ-linolenic acid and DHA are the principal fatty acids present in brain and retinal tissues, which can promote the differentiation function of immunomodulatory cells, enhance immune response, maintain the integrity of the intestinal mucosal barrier, and reduce mucosal oxidative damage to promote intestinal health. This study suggests that feeding FDG diets may stimulate the biosynthesis of unsaturated fatty acids and linoleic acid metabolism in Guanling crossbred cattle, resulting in elevated levels of these metabolites. These may have anti-inflammatory, immunomodulatory, and intestinal health-promoting effects, thereby benefiting the growth and development of Guanling crossbred cattle. However, it is important to note that the specific effects of the FDG diets may vary according to individual differences, background health status, and other environmental factors. Therefore, further studies are still needed to investigate the effects and potential mechanisms of action of these metabolites on the gut health of animals.
5. Conclusion
In conclusion, our results suggest that the FDG diet may promote the gut health and immunity of Guanling crossbred cattle by modulating the metabolic patterns of lipid compounds and lipid-related metabolic pathways, such as biosynthesis of unsaturated fatty acids and linoleic acid metabolism, and primary bile acid biosynthesis. However, this study is merely a pilot study, and further in-depth analysis with a large number of animals is needed to confirm the specific regulatory functions and their potential mechanisms of action. Our findings contribute to further understanding of the nutritive and regulatory functions of the FDG diet for Guanling crossbred cattle and their underlying mechanisms as well as provide some references for FDG applications as a feed resource.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The animal studies were approved by EAE-GZu-2020-E018, 2 September 2020. 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
XL, EZ, ZCheng, and CC conceived the study. XL, TZ, DX, and MZ performed the experiments. XL, JZ, RZ, GH, EZ, ZCheng, and CC analyzed experimental results and data. GH, BZ, KW, CC, and SM assisted with animal experiments. XL wrote the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by grants from the Guizhou Provincial Department of Agriculture Project (No. [2022] 163), the Guizhou Provincial Support Fund of Science and Technology (No. [2021]5646), the Guizhou Provincial Science and Technology Projects (QKHJC-ZK [2022] YB158), the Cultivating Project of Guizhou University (GDPY[2020]85), and the Scientific Research Fund for talents recruiting of Guizhou University (GDRJHZ [2020]63).
Acknowledgments
We kindly acknowledge the Kweichow Moutai Group for providing distillers' grains and Guizhou Cattle Industry Group Co., Ltd. for providing the Guanling crossbred cattle.
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.1238064/full#supplementary-material
References
1. Xu D, Chen C, Zhu M, Wang Y. Effect of fermented distiller's grains feeding level on growth performance, blood biochemical index and muscle amino acid composition of Guangling crossbred cattle. Chin J Anim Sci. (2023) 59:212–6.
2. Chen Y, Wei Y, Li S, Yin M, Zuo F, Huang M, et al. Effect of different proportions of peanut seedlings and corn flour on nutritional value in anaerobic fermentation of highly type baijiu distiller's grains. Chin J Anim Nutr. (2023) 35:2690–701.
3. Cao G, Wang Z, Hu R, Wang X, Ma J, Zhou X, et al. Effect of fermented distiller's grains on growth performance, nutrient apparent digestibility and rumen fermentation of finishing Yaks. Chin J Anim Nutr. (2020) 32:3255–62.
4. Li Z, Zhang T, Su J, Wang D, Kong X. Effect of fermented distiller's grains on growth performance, digestibiliy metabolism and nitrogen excretion of weaned piglets. Chin J Anim Nutr. (2018) 30:2807–16.
5. Ma Y, Zhou J, Zhang H, Wang J, Wu S, Qi G, et al. Effects of dietary yeast culture of white distiller's grains in different energy level diets on growth performance, carcass composition and serum biochemical indices of broilers. Chin J Anim Nutr. (2018) 30:274–83.
6. Wang Y, Li Z. Effect of different levels of probiotic fermented distiller's grains on growth performance, slaughter performance and meat quality of fattening cattle. China Feed. (2020) 11:117–20. doi: 10.15906/j.cnki.cn11-2975/s.20201127
7. Zhou X, Wang Z, Hui R, Wang X, Zhu Y, Cao G, et al. Effects of different types of white distiller's grains on growth performance, nutrient apparent digestibility, serum biochemical indices and rumen fermentation parameters of dairy bulls. Chin J Anim Nutr.. (2020) 32:3708–16.
8. Liu W. Effects of Vinasse Fattened on the Nutritional Component Meat Quality of the Simmtental Beef. Hohhot: Inner Mongolia Agricultural University. (2015). Available online at: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=1015427840.nh&DbName=CMFD2015 (accessed December 15, 2015).
9. Yang Y, Wu X, Hu R, Wang C, Xu E. Effect of different proportions of Moutai distillers' grains bio-material on slaughter performance and muscle lipid metabolomics of crossbred cattle. China Feed. (2023) 2023:1–7.
10. Franchimont D, Louis E. Genetics of inflammatory bowel disease. Bull Mem Acad R Med Belg. (2010) 165:435–449.
11. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. (2008) 320:1647–51. doi: 10.1126/science.1155725
12. Sato T, Vries RG, Snippert HJ, van de Wetering M, Barker N, Stange DE, et al. Single Lgr5 stem cells build crypt-villus structures in vitro without a mesenchymal niche. Nature. (2009) 459:262–5. doi: 10.1038/nature07935
13. Spence JR, Mayhew CN, Rankin SA, Kuhar MF, Vallance JE, Tolle K, et al. Directed differentiation of human pluripotent stem cells into intestinal tissue in vitro. Nature. (2011) 470:105–9. doi: 10.1038/nature09691
14. Zhang R, Zhang J, Dang W, Irwin DM, Wang Z, Zhang, S. Unveiling the biogeography and potential functions of the intestinal digesta and mucosa-associated microbiome of donkeys. Front Microbiol. (2020) 11:596882. doi: 10.3389/fmicb.2020.596882
15. Qing J, Li C, Hu X, Song W, Tirichen H, Yaigoub H, et al. Differentiation of T helper 17 cells may mediate the abnormal humoral immunity in IgA nephropathy and inflammatory bowel disease based on shared genetic effects. Front Immunol. (2022) 13:916934. doi: 10.3389/fimmu.2022.916934
16. Zhuang Y, Huang H, Liu S, Liu F, Tu Q, Yin Y, et al. Resveratrol improves growth performance, intestinal morphology, and microbiota composition and metabolism in mice. Front Microbiol. (2021) 12:726878. doi: 10.3389/fmicb.2021.726878
17. Deng F, Lin ZB, Sun QS, Min Y, Zhang Y, Chen Y, et al. The role of intestinal microbiota and its metabolites in intestinal and extraintestinal organ injury induced by intestinal ischemia reperfusion injury. Int J Biol Sci. (2022) 18:3981–92. doi: 10.7150/ijbs.71491
18. Xing PY, Pettersson S, Kundu P. Microbial metabolites and intestinal stem cells tune intestinal homeostasis. Proteomics. (2020) 20:e1800419. doi: 10.1002/pmic.201800419
19. Koh A, Vadder De, Kovatcheva-Datchary F, Backhed P. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. (2016) 165:1332–45. doi: 10.1016/j.cell.2016.05.041
20. Shi H, Yang J, Li J. Pomegranate peel polyphenols interaction with intestinal flora and its metabolic transformation. Xenobiotica. (2022) 52:442–52. doi: 10.1080/00498254.2022.2073291
21. Cheng Q, Xu D, Chen Y, Zhu M, Fan X, Li M, et al. Influence of fermented-moutai distillers' grain on growth performance, meat quality, and blood metabolites of finishing cattle. Front Vet Sci. (2022) 9:874453. doi: 10.3389/fvets.2022.874453
22. Zhang J, Cheng Z, Zhang T, Xu D, Zhu M, Luo X, et al. Plasma metabolomic profiling after feeding dried distiller's grains with solubles in different cattle breeds. Int J Mol Sci. (2023) 24:13. doi: 10.3390/ijms241310677
23. Zhang Y, Guo L, Luo Y, Zhang Z, Li Y, Zhu R, et al. Effects of high-doses of distillers' grains on rumen fermentation, microflora structure and inflammatory genes. Chin J Vet Sci. (2023) 43:748–56.
24. Khusro A, Aarti C. Metabolic heterogeneity and techno-functional attributes of fermented foods-associated coagulase-negative staphylococci. Food Microbiol. (2022) 105:104028. doi: 10.1016/j.fm.2022.104028
25. Singh AK, Rehal J, Kaur A, Jyot G. Enhancement of attributes of cereals by germination and fermentation: a review. Crit Rev Food Sci Nutr. (2015) 55:1575–89. doi: 10.1080/10408398.2012.706661
26. Sun R. Effect of different treatment methods of wine lees on the weight gain effect of beef cattle. China Herbi Sci. (2017) 37:66–8.
27. Zheng H, Yin X, Liu Z, Zhang Y, Liu D, Wang Y, et al. Effect of fermented distillers grains bran mixture on growth performance, slaughter performance and immune organ index of broilers. China Poultry. (2023) 19:1–7.
28. Ito KR, Sato T, Goto H, Sato K, Watanabe J, Yokoo M, et al. Utilization of sake lees as broiler feedstuff and its effects on growth performance and intestinal immunity. J Poult Sci. (2022) 59:247–59. doi: 10.2141/jpsa.0210087
29. Liu Y, Feng J, Wang Y, Lv J, Li J, Guo L, et al. Fermented corn-soybean meal mixed feed modulates intestinal morphology, barrier functions and cecal microbiota in laying hens. Animals. (2021) 11:11. doi: 10.3390/ani11113059
30. Rho Y, Wey D, Zhu C, Kiarie E, Moran K, van Heugten E, et al. Growth performance, gastrointestinal and digestibility responses in growing pigs when fed corn-soybean meal-based diets with corn DDGS treated with fiber degrading enzymes with or without liquid fermentation. J Anim Sci. (2018) 96:5188–97. doi: 10.1093/jas/sky369
31. Zhong-Hua L, Xiao-Ge H, Jin-Hui Z, Le HE. Liquid fermentation of ganoderma applanatum and antioxidant activity of exopolysaccharides. Open Biomed Eng J. (2015) 9:224–7. doi: 10.2174/1874120701509010224
32. Angelini G, Salinari S, Castagneto-Gissey L, Bertuzzi A, Casella-Mariolo J, Ahlin S, et al. Small intestinal metabolism is central to whole-body insulin resistance. Gut. (2021) 70:1098–109. doi: 10.1136/gutjnl-2020-322073
33. Li X, Cai Z, Liu J, Wang N, Zhu X, Lu Z, et al. Antiobesity effect of L-arabinose via ameliorating insulin resistance and modulating gut microbiota in obese mice. Nutrition. (2023) 111:112041. doi: 10.1016/j.nut.2023.112041
34. Xu L, Li Y, Wei Z, Bai R, Gao G, Sun W, et al. Chenodeoxycholic acid (CDCA) promoted intestinal epithelial cell proliferation by regulating cell cycle progression and mitochondrial biogenesis in ipec-j2 cells. Antioxidants. (2022) 11:11. doi: 10.3390/antiox11112285
35. Du Xiang J X, Li Y, Ji R, Xu H, Mai K, et al. Molecular cloning and characterization of farnesoid X receptor from large yellow croaker (Larimichthys crocea) and the effect of dietary cdca on the expression of inflammatory genes in intestine and spleen. Comp Biochem Physiol B-Biochem Mol Biol. (2018) 216:10–7. doi: 10.1016/j.cbpb.2017.09.007
36. Li X, Yuan Z, Liu R, Hassan HM, Yang H, Sun R, et al. UDCA and CDCA alleviate 17alpha-ethinylestradiol-induced cholestasis through PKA-AMPK pathways in rats. Toxicol Appl Pharmacol. (2016) 311:12–25. doi: 10.1016/j.taap.2016.10.011
37. Klievik BJ, Tyrrell AD, Chen CT, Bazinet RP. Measuring brain docosahexaenoic acid turnover as a marker of metabolic consumption. Pharmacol Ther. (2023) 248:108437. doi: 10.1016/j.pharmthera.2023.108437
38. Ahmad F, Ali SS, Usmanghani K, Ali M. Lipid contents of marine fish: carcharhinus melanopterus (black shark) and lutjanus johnii (hira). Pak J Pharm Sci. (1991) 4:91–101.
39. Howton DR, Mead JF. Metabolism of the essential fatty acids. X. Conversion of 8,11,14-eicosatrienoic acid to arachidonic acid in the rat. Annu Prog Rep Univ Calif Lab Nucl Radiat Biol. (1960) 456:1–11. doi: 10.1016/S0021-9258(18)64476-X
40. Xu C, Gu L, Hu L, Jiang C, Li Q, Sun L, et al. FADS1-arachidonic acid axis enhances arachidonic acid metabolism by altering intestinal microecology in colorectal cancer. Nat Commun. (2023) 14:2042. doi: 10.1038/s41467-023-37590-x
41. Zhao J, Nishiumi S, Tagawa R, Yano Y, Inoue J, Hoshi N, et al. Adrenic acid induces oxidative stress in hepatocytes. Biochem Biophys Res Commun. (2020) 532: 620–5. doi: 10.1016/j.bbrc.2020.08.102
42. Zou Q, Wei M, Zhang N, Niu X, Weng C, Deng ZY, et al. Different influences of trans fatty acids on the phospholipase A2 and arachidonic acid metabolic pathway in hepatocytes. J Agric Food Chem. (2021) 69:4120–33. doi: 10.1021/acs.jafc.1c01097
43. Sergeant S, Rahbar E, Chilton FH. Gamma-linolenic acid, dihommo-gamma linolenic, eicosanoids and inflammatory processes. Eur J Pharmacol. (2016) 785:77–86. doi: 10.1016/j.ejphar.2016.04.020
44. Kikut J, Komorniak N, Zietek M, Palma J, Szczuko, M. Inflammation with the participation of arachidonic (AA) and linoleic acid (LA) derivatives (HETEs and HODEs) is necessary in the course of a normal reproductive cycle and pregnancy. J Reprod Immunol. (2020) 141:103177. doi: 10.1016/j.jri.2020.103177
45. Wang Q, Lin Y, Sheng X, Xu J, Hou X, Li Y, et al. Arachidonic acid promotes intestinal regeneration by activating wnt signaling. Stem Cell Rep. (2020) 15:374–8. doi: 10.1016/j.stemcr.2020.06.009
46. Andres CG, Koster J, Souza J, Laguna J, Mavangira V, Nelli RK, et al. Lipolysis modulates the biosynthesis of inflammatory lipid mediators derived from linoleic acid in adipose tissue of periparturient dairy cows. J Dairy Sci. (2020) 103:1944–955. doi: 10.3168/jds.2019-17256
Keywords: fermented distiller's grains (FDG), Guanling crossbred cattle, intestinal metabolic profile, metabolites, metabolic pathways
Citation: Luo X, Zhang T, Xu D, Zhu M, Zhang J, Zhang R, Hu Q, Wang Y, He G, Chen Z, Mei S, Zhou B, Wang K, Chen C, Zhu E and Cheng Z (2023) The effect of feeding fermented distillers' grains diet on the intestinal metabolic profile of Guanling crossbred cattle. Front. Vet. Sci. 10:1238064. doi: 10.3389/fvets.2023.1238064
Received: 10 June 2023; Accepted: 25 September 2023;
Published: 20 October 2023.
Edited by:
Christian Visscher, University of Veterinary Medicine Hannover, GermanyReviewed by:
Kang Xu, Chinese Academy of Sciences (CAS), ChinaAlok Wankar, Maharashtra Animal and Fishery Sciences University, India
Copyright © 2023 Luo, Zhang, Xu, Zhu, Zhang, Zhang, Hu, Wang, He, Chen, Mei, Zhou, Wang, Chen, Zhu and Cheng. 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: Erpeng Zhu, zhu13782701756@126.com; Zhentao Cheng, ztcheng@gzu.edu.cn