Multi-omics reveal mechanisms of high enteral starch diet mediated colonic dysbiosis via microbiome-host interactions in young ruminant

Background Although rumen development is crucial, hindgut undertakes a significant role in young ruminants’ physiological development. High-starch diet is usually used to accelerate rumen development for young ruminants, but always leading to the enteral starch overload and hindgut dysbiosis. However, the mechanism behind remains unclear. The combination of colonic transcriptome, colonic luminal metabolome, and metagenome together with histological analysis was conducted using a goat model, with the aim to identify the potential molecular mechanisms behind the disrupted hindgut homeostasis by overload starch in young ruminants. Result Compared with low enteral starch diet (LES), high enteral starch diet (HES)-fed goats had significantly higher colonic pathology scores, and serum diamine oxidase activity, and meanwhile significantly decreased colonic mucosal Mucin-2 (MUC2) protein expression and fecal scores, evidencing the HES-triggered colonic systemic inflammation. The bacterial taxa Prevotella sp. P4-67, Prevotella sp. PINT, and Bacteroides sp. CAG:927, together with fungal taxa Fusarium vanettenii, Neocallimastix californiae, Fusarium sp. AF-8, Hypoxylon sp. EC38, and Fusarium pseudograminearum, and the involved microbial immune pathways including the “T cell receptor signaling pathway” were higher in the colon of HES goats. The integrated metagenome and host transcriptome analysis revealed that these taxa were associated with enhanced pathogenic ability, antigen processing and presentation, and stimulated T helper 2 cell (TH2)-mediated cytokine secretion functions in the colon of HES goats. Further luminal metabolomics analysis showed increased relative content of chenodeoxycholic acid (CDCA) and deoxycholic acid (DCA), and decreased the relative content of hypoxanthine in colonic digesta of HES goats. These altered metabolites contributed to enhancing the expression of TH2-mediated inflammatory-related cytokine secretion including GATA Binding Protein 3 (GATA3), IL-5, and IL-13. Using the linear mixed effect model, the variation of MUC2 biosynthesis explained by the colonic bacteria, bacterial functions, fungi, fungal functions, and metabolites were 21.92, 20.76, 19.43, 12.08, and 44.22%, respectively. The variation of pathology scores explained by the colonic bacterial functions, fungal functions, and metabolites were 15.35, 17.61, and 57.06%. Conclusions Our findings revealed that enteral starch overload can trigger interrupted hindgut host-microbiome homeostasis that led to impaired mucosal, destroyed colonic water absorption, and TH2-mediated inflammatory process. Except for the colonic metabolites mostly contribute to the impaired mucosa, the nonnegligible contribution from fungi deserves more future studies focused on the fungal functions in hindgut dysbiosis of young ruminants. Video Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s40168-024-01760-w.


Background
During the early life, due to the dramatic nutritional and physiological changes in this period, the underdeveloped digestive tracts of young ruminants could undergo gastrointestinal dysfunction that significantly affects their growth and health.The dysfunction of gastrointestine and immunosuppression in this period leads to 14.17% of morbidity of young ruminants in China [1] and causes about 32% of deaths in the US [2], contributing to significant economic loss in the livestock industry.Thus, promoting the maturation of the digestive tract has become the major strategy to cope with this problem.At the preweaning phase, due to the immature digestive and immune system in neonatal and young ruminants, colon provides the main ecological niches for microbial colonization and fermentation [3][4][5] and is vulnerable to external stimuli [6][7][8].Colon undertakes a significant role in young ruminants' energy metabolism, immune system, and physiological development during this period [9,10].The implementation of starch to feed young ruminants to stimulate their growth and the rumen development has become a widely accepted practice.A large amount of incompletely degraded starch bypasses the rumen and flows into the hindgut, thereby triggering excessive fermentation in the hindgut [11], which is often accompanied by microbial dysbiosis, hindgut acidosis, colonic mucosal damage, and inflammation [12][13][14].Although hindgut immunity plays an important role in maintaining ruminants' health during this period which can directly influence the ruminants' growth and performance [15,16], the exact mechanism on its responses to high enteral starch-mediated dysbiosis has not been well-defined, preventing the effective intervention strategies.
The hindgut microbiome plays crucial roles in animal health, nutrient absorption and metabolism, gastrointestinal development, and immune function [17,18].Although these have been well studied in human, mice, and other monogastric animals, the understanding of immunoregulatory functions in dysbiosis in hindgut of ruminants during early life is still limited.Several studies have demonstrated the linkage between the changes in the gut microbiome and gastrointestinal dysfunction in ruminants through affecting the microbial metabolites such as short-chain fatty acids (SCFA), bile acids, and other metabolites [19,20].However, the majority of the host-microbiome research during early life focused on rumen development and immune functions in the small intestine.Recent study has shown that colon microbiota differs from rumen and small intestine [21], and colon microbial fermentation plays an important role in providing energy for young ruminants [22].In addition, the colonic microbiome and transcriptome can be significantly affected by the early life feeding strategies such as delayed colostrum feeding [5,23].However, these studies only investigated the colonic microbiota, and/or host transcriptome of calves at 2 days of age, and lack comprehensive knowledge about long-term follow-up assessments on colonic microbial-host interactions.
In this study, we hypothesized that, in young ruminants, the enteral starch can lead to the interrupted colon homeostasis via affecting host-microbe interactions; by altering colonic microbiota and its metabolism, thereby affecting the colonic mucosal barrier and immune functions at gene and protein levels which could result in systemic inflammation of the young ruminants.Therefore, we assessed the host-microbe interactions by profiling colonic microbiome, metabolome, and host transcriptome using integrated omics and quantitative gene and protein expressions in young ruminants using a goat model fed with high enteral starch diet and evaluated their contributions to the host phenotypes differences including fecal scores, MUC2 biosynthesis, pathology scores, serum lipopolysaccharide (LPS) concentration, and serum diamine oxidase (DAO) activity.

Ethics approval statement
This study was conducted in accordance with the recommendations of the Administration of Affairs Concerning Experimental Animals (Ministry of Science and Technology, China, revised in 2004).The animal using protocol was approved by the Institutional Animal Care and Use Committee of the Northwest A&F University.

Animal experiment
The flow chart for the overall experimental design was demonstrated in the supplementary materials (Fig. S1).Forty healthy, weaned 3-month-old male goats (BW = 13.6 ± 0.23 kg) were randomly enrolled to receive either a whole corn grain diet (n = 20) or crushed corn grain diet (n = 20).Considering the processing methods of the corn affects the starch digestion in the gut, the hindgut starch content showed that the crushed corn grain diet was the high hindgut-enteral starch content (HES) group, while the whole corn grain diet was the low hindgut-enteral starch content (LES) group (Fig. S1).
All goats were housed individually in tie-stall barns and fed twice a day (0800 and 1600) with free access to water.The forage-to-concentrate ratio of the two diets was both set as 60:40 and the detailed ingredient compositions and nutrition levels of the feed are presented in Table S1.Feed samples and the dry matter intake were continually collected for 3 days every week throughout the trial.The feeding trial lasted for 90 days, and the fecal scores were evaluated for continuously 7 days (84-90 days) using a 5-grade fecal scores evaluation system that adapted from Woolsoncroft et al. [24] and Ireland-Perry and Stallings [25] (Fig. S2).The standard of this system is as follows: 5 = hard, dry pellets in a small, hard mass; 4 = hard, formable stool like pinecone that retains the shape; 3 = soft, formable stool like stone that retains its shape; 2 = soft, deformable stool that likes human normal feces; 1 = watery, liquid stool.Based on the fecal scores difference, 8 goats from each group (in total 16 goats) were selected for slaughter and sampling at the end of the trial (on day 90).

Sample collection
Before the slaughter on day 90, blood samples were collected from jugular vein into 5-mL endotoxin free vacuum blood collection tubes and were then put into 37 °C water bath for 15 min and centrifuged at 3500 × g for 15 min at 4 °C for serum collection.Based on the mean values of fecal scores of two groups, eight represented LES goats and eight represented HES goats, were selected for slaughter and sampling.After slaughter, the rumen fluid samples were filtered and collected to determine the pH and SCFA content following the methods of Shen et al. [26].The colonic digesta were collected following the procedures as described in Ye et al. [27].Colonic digesta of each goat was mixed well and pH was assessed immediately according to the description of Liu et al. [28].The mixed colonic digesta were then divided into three portions; (1) stored in liquid nitrogen for microbial DNA, metabolite extraction, and SCFA analysis [26,29]; (2) processed and stored at − 20 °C for LPS analysis [27]; (3) used to analyze the nutrients' composition [28].For colonic tissues, about 6 cm of colonic tissue of each goat was harvested at the middle of the colon and rinsed 3 times with ice-cold sterile 0.01 M PBS buffer (pH7.2-7.4,Beijing Solarbio Science & Technology Co., Ltd, China) within 5 min of slaughter.The colonic tissue was then divided into three parts: (1) scraped using sterile slides and stored in liquid nitrogen for RNA isolation [30]; (2) fixed in 4% paraformaldehyde to make hematoxylin and eosin (H&E) staining sections for pathologic observation; (3) fixed with methanol-Carnoy's fixative [60% (v/v) dry methanol, 30% (v/v) chloroform, 10% (v/v) glacial acetic acid] for Mucin-2 (MUC2) immunofluorescence and immunohistochemistry (IHC) detection.

Colonic histopathology
Colonic H&E-stained sections were used for the grading of colonic inflammation as reported by Neurath et al. [31].In brief, 0, no signs of inflammation; 1, very low level; 2, low level of leukocytic infiltration; 3, high level of leukocytic infiltration, high vascular density, thickening of the colon wall; 4, transmural infiltrations, loss of goblet cells, high vascular density, thickening of the colon wall.About 9 images per sample were taken with Nikon ECLIPSE Ni-U (Nikon Co., Ltd, Japan) at × 100 magnification, and the average scores of each sample were calculated.

Colonic MUC2 detection using immunofluorescence and IHC staining
Immunofluorescence analysis of colon tissue was performed following the procedures described by Earle et al. [32] and Kong et al. [33].Briefly, after deparaffinization, rehydration, and antigen retrieval, slides were blocked with bovine serum albumin for 30 min at room temperature following incubation with Anti-MUC2 antibody (EPR6145, ab134119, Abcam, UK) at 4 °C overnight.Then, slides were washed three times in PBS and incubated with secondary antibody (ZF-0311, FITC-conjugated goat anti-rabbit, ZSGB-BIO, INC, China) at 37 °C for 30 min.After incubation with the secondary antibody and DAPI, slides were washed three times in PBS, dried, and mounted with anti-fade reagent (AR1109; Boster Bioengineering, INC, China).The images were collected using FluoView FV1000 confocal laser microscope with the FV1000 software (Olympus Optical Technology Co., Ltd, Japan).IHC staining for MUC2 was performed with an initial incubation with MUC2 rabbit polyclonal (Anti-MUC2 Rabbit pAb, GB11344, Servicebio, INC, China) at a dilution of 1:100 followed by the application of the secondary biotinylated goat anti-rabbit IgG (SP-9001, ZSGB-BIO, INC, China).The slides were captured with BA400 Digital microscopic camera system (Motic China Group Co., Ltd, China) and scanned at × 400 magnification.Semi-quantitative image analysis was performed with the open-source software ImageJ, and the IHC profiler plug-in developed by Varghese et al. [34].Software operation and data collection were performed as previously described by Alagaratnam et al. [35].

Transcriptome analysis of colonic mucosa
Total RNA was extracted from colonic tissue using Trizol reagent (Invitrogen, Carlsbad, INC, USA).The RNA amount and purity were quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, INC, USA) and RNA integrity was assessed using Bioanalyzer 2100 (Agilent, INC, USA) and confirmed with denaturing agarose gel electrophoresis.The total RNAs of twelve goats (with RIN number > 7.0, 6 from each group) with differed fecal scores were subjected for transcriptome analysis using RNA-seq with the 2 × 150 bp paired-end sequencing (PE150) on an Illumina Novaseq ™ 6000 (LC-Bio Technology CO., Ltd, China) according to a previous study [36].The capra hircus V1 (ARS1 (GCF_001704415.1))gene annotation list was used as the reference genome.The expression level of each transcript was calculated as the fragments per kilobase of exon per million mapped reads (FRKM) using StringTie and Ballgown [37].

Metabolomic analysis of colonic digesta
The same 12 colon digesta were thawed at 4 °C, and 1000 μL methanol/ acetonitrile/ H 2 O (2:2:1, v/v/v) were added to homogenized solution for metabolite extraction.Then they were re-dissolved in 100 μL acetonitrile/ water (1:1, v/v) solvent and subjected to LC-MS analysis using an UHPLC (1290 Infinity LC, Agilent Technologies, INC, Germany), coupled to a quadrupole time-of-flight (AB Sciex TripleTOF 6600) in Shanghai Applied Protein Technology Co., Ltd.The raw MS data (wiff.scanfiles) were converted to MzXML files using ProteoWizard MSConvert before importing into freely available XCMS software [48].For peak picking, the following parameters were used: centWave m/z = 25 ppm, peakwidth = c (10, 60), prefilter = c (10, 100).For peak grouping, bw = 5, mzwid = 0.025, minfrac = 0.5 were used.CAMERA (Collection of Algorithms of Metabolite Profile Annotation) [49] was used for annotation of adducts.In the extracted ion features, only the variables having more than 50% of the nonzero measurement values in at least one group were kept.Compound identification of metabolites was performed by comparing of accuracy m/z value (< 25 ppm), and MS/MS spectra with a published database established with available authentic standards [50,51].After normalized to total peak intensity, the processed data were analyzed by R package ropls [53], where it was subjected to multivariate data analysis, including Pareto-scaled principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA).According to Li et al. [53], the seven-fold cross-validation and response permutation testing were used to evaluate the robustness of the model.The variable importance in the projection (VIP) value of each variable in the OPLS-DA model was calculated to indicate its contribution to the classification.Origin assessment of colonic metabolites was based on MetOrigin database (2022-08 version, Capra hircus, https:// metor igin.metbioin forma tics.cn) [53].

Analysis of LPS concentration and DAO activity
LPS of colonic digesta and serum was detected using Bioendo EC Endotoxin Test Kit (End-point Chromogenic Assay, Xiamen Bioendo Technology Co., Ltd, China) following the manufacturer's procedure and as reported by Chen et al. [53].The concentration of digesta and serum LPS was calculated using the standard curve and expressed as EU/g and EU/mL, respectively.Serum DAO activity was measured with DAO assay kit (A088-1-1, Njjcbio, INC, China) with ultraviolet chromatometry.The serum DAO activity was expressed as U/L.All measurements were conducted with three replicates.

Omics-based explainability of phenotype variations
The variation in host phenotypes that could be explained by the colonic microbial variance, metabolic variance, and KOs variance were further evaluated using the methods described in Xue et al. [53].Colonic microbiome, metabolites, and KOs were normalized to construct the relationship matrix, respectively [53].The liner mixed model (LMM) was conducted as y ij = μ + starch i + x j + e ij with lme4 package in R [53].y ijk is the fecal scores/ MUC2 expression/ pathology scores/ serum LPS; μ is the model intercept; starch i is the fixed effect; x j is the colonic bacteria/bacterial functions/ fungi/ fungal functions/ metabolites random effect for the jth animal∼NID (0, Xσ 2 X ), σ 2 X is the colonic bacteria/ bacterial functions/ fungi/ fungal functions/ metabolites variance.The X is the colonic bacteria/ bacterial functions/ fungi/ fungal functions/ metabolites relationship matrix; and e ij is the residual effects.

Statistical analysis
Except for the omics data, dry matter intake (DMI), digesta nutrients' compositions, pH, SCFAs, fecal scores, histopathology scores, MUC2 IHC assessment, LPS, DAO, and the mRNA relative abundance (qPCR) were analyzed by using Grubbs' test and assessed by the independent sample T-test in SPSS 20.0 (SPSS INC, USA).For the omics data, P value was corrected by using the false discovery rate method (FDR) [53].For host transcriptome, the differentially expressed genes (DEGs) were determined with fold change ≥ 2 or ≤ 0.5 and FDR < 0.05 using DESeq2 in R package [33,53], and then DEGs were subjected to KEGG enrichment analysis.0.05 < FDR < 0.1 was referred to as approaching significance.
Based on metagenomics data, the microbial taxa that were present in more than 60% of the samples for either group were retained for downstream analysis.Colonic bacterial and fungal species were compared by LEfSe with LDA > 2 and P < 0.05 [53].Abundances of microbial functions, including, but not limited to, KEGG Orthology (KO) terms, pathways, and CAZymes were normalized into reads per kilobase million (RPKM) in at least 60% of the samples for either group.The differentially functional genes were selected by using LEfSe methods (LDA > 2, P < 0.05).Metabolomics data was log2 transformed [53].Metabolites with the VIP > 1 were further applied to T-test at univariate level to measure the significance of each metabolite.The FDR < 0.05 was considered as statistically significant.Considering the high diversity, high dimensionality, and discrete distribution of omics data, we used the Spearman correlation test in the present study to explore the relationships among the phenotypes, nutrients, microbiome, metabolites, and host transcripts.|R|> 0.5, P < 0.05 was used to identify significant correlations.Cytoscape 3.6.1 was used to construct the network.

Enteral starch altered the phenotypes of hindgut in goats
Primarily, through the fecal assessment of 40 goats, we found that HES significantly decreased fecal scores of goats (P = 0.001; Fig. 1A).Based on the mean values of fecal scores of two groups, eight LES goats (BW = 19.77± 2.73 kg, fecal scores = 4.42 ± 0.13) and eight HES goats (BW = 17.69 ± 5.55 kg, fecal scores = 3.42 ± 0.30) were selected for slaughter and sampling.The starch content of postrumen digesta (0.05 < P < 0.1, Fig. S3A), colonic digesta (0.05 < P < 0.1, Fig. S3B and C), and colonic pathology scores (P < 0.05; Fig. 1B) were higher in the HES goats.Immunofluorescence and IHC results showed the lower MUC2 protein expression in the HES goats (Fig. 1C).Further, the serum level of DAO (P < 0.05) and LPS (0.05 < P < 0.1) in the HES goats were higher than those in LES goats (Fig. 1D).Moreover, no significant differences in the ruminal pH, total SCFA concentration, and proportion of acetate, propionate, and butyrate were identified between the two groups (P > 0.05; Table S3).
Due to the upregulation of c-Maf in the colon of HES goats, we further detected the expression levels of key genes in colonic mucosa that are associated with T H 1 and T H 2 cells (Fig. 4F) and found that the expression of IL-2 of T H 1 cells was significantly downregulated in the HES goats (P < 0.05).The expression of the master transcription factor GATA3, c-Maf inducing protein (CMIP), IL-5, and IL-13 of T H 2 cells were significantly upregulated in the HES goats (P < 0.05).Similarly, the expression of mucosal barrier and cytokine genes was evaluated using qPCR (Fig. 4A).The results showed upregulated DUOX2 (P < 0.05) and downregulated Claudin4 (0.05 < P < 0.1) and Occludin (P < 0.05) in the colon of HES goats.HES increased the expression of IL-6 (P < 0.05), IL-22 (P < 0.05), IFN-γ (0.05 < P < 0.1), and C3 (P < 0.01) in colonic mucosa.Besides, we also found that HES decreased the expression of AQP1 (P < 0.05) and AQP8 (P < 0.05) (Fig. 1A).
Using linear mixed effect model (see "Methods"), the variation of phenotypes explained by the colonic microbiota, microbial functions, and metabolites are shown in Fig. 6.The variation of fecal scores explained by the colonic bacteria and fungi were 40.51 and 28.12%, respectively.The variation of MUC2 biosynthesis explained by the colonic bacteria, bacterial functions, fungi, fungal functions, and metabolites were 21.92, 20.76, 19.43, 12.08, and 44.22%, respectively.The variation of pathology scores explained by the colonic bacterial functions, fungal functions, and metabolites were 15.35, 17.61, and 57.06, respectively.The variation of serum LPS explained by the colonic bacteria was 33.56%.

Discussion
Colon plays an important role in energy metabolism and intestinal health of ruminants during early life and other physiological phases [13,14,20,53,53].By using a high-grain or high-starch diet, colonic inflammation has been linked to an altered microbiota [13,53,53], but the possibility that colonic dysbiosis may be derived from the potential effects of a high-starch diet on ruminal development and fermentation disorders has not been excluded.By maintaining the same dietary starch content and changing only the corn processing methods, which simultaneously changed the rate of starch digestion in the rumen and passed through the rumen [53], the HES diet with crushed corn increased the hindgut starch content but did not significantly change the ruminal fermentation indices, development, or health conditions compared with LES diets with whole corn.Hence, the identified colonic inflammation was mainly derived from colonic microbiota disorders and may not be affected by ruminal dysfunction.Notably, our study is one of the few studies that demonstrated altered colonic microbiome, metabolome, and host transcriptome in response to HES feeding and its associated pathologies in young ruminants.Metagenomics sequencing showed the dynamic changes of colonic microbiome in the goats and found the impacts of colonic microbial dysbiosis in response to different enteral starch diets.The observed significant shifts in host-microbiome interactions led to impaired mucosa and systemic inflammation.
First, our study revealed that changes in the microbiome at both taxonomic and functional levels were associated with decreased colonic MUC2 expression in response to HES feeding.Lower MUC2 expression indicated the depletion of major structural mucus components or reduction of the response of goblet cells to microbial challenge, abnormal penetrability, and breakdown of intestinal barrier function induces colonic inflammation [53,53].Generally, Muribaculaceae, Bifidobacterium, and Bacteroidales can degrade all or most O-glycan mucosal sugars [53].But in the present study, no microbial association with colonic MUC2 expression was found.The analysis of omics basing explainability of phenotype variations showed that the variation of MUC2 biosynthesis explained by the colonic bacteria, bacterial functions, fungi, fungal functions, and metabolites were 21.92, 20.76, 19.43, 12.08, and 44.22%, respectively.Thus exhibits a comprehensive effect of multiple factors in contributing or positively interacting with other microorganisms responsible for colonic MUC2 depletion.Previous studies showed that the DCA and CDCA breakdown the epithelial tight junctions and contribute to inflammatory bowel disease (IBD) and irritable bowel syndrome (IBS) [53,53,53,53].The DCA and CDCA were enriched in HES group, implying that the presence of higher level of DCA and CDCA increases the risk of developing mucosal barrier dysfunction and inflammation.Together with the higher serum LPS in HES goats, it suggests that enteric starch increased mucosal permeability by selectively increasing the abundance of bacteria with enhanced mucin degradation capacity.
Our study further focused on the roles of altered microbiome in affecting colonic mucosal immunologic state between LES and HES goats.First, overall analyzed the key differential genes, TRAV and c-Maf, that separately severing as an antigen receptor of CD4 + /8 + T cells or the promoter of the T H 2 cell maturation [53], were linked to upregulated antigen processing and presentation and inflammatory bowel disease (IBD) pathways in colonic mucosa of HES goats.Further, since the IL-6 activates both T H 1 and T H 2 cell continually during inflammation, the role of T H 1/T H 2 polarization in mucosal immunity was then tested.IL-12 induces T H 1 cell differentiation by activation and phosphorylation of the signal transducer  [53].As an opportunistic pathogen, Prevotella has been proven to infect the host with other microorganisms when the ecological conditions change [53].Previous study showed that the colonization of the Prevotella and Muribaculaceae could exacerbate intestinal inflammation [53].Our study demonstrated that Bacteroides sp.CAG:927, 2 species of Prevotella (Prevotella sp.PINT and Prevotella sp.P4-67), and 6 differential bacterial functions (ko04215, ko04660, ko04664, ko05170, ko05210, and ko05323) were closely related to the expression of immune-related DEGs that includes TRAV, KIR3DL, and c-MAF for antigen presentation; c-MAF, GATA3, IL-5, and IL-13 for T H 2-mediated inflammatory process; A1AT1 and C3 for complement and coagulation cascades, indicating that the T H 2-mediated colonic inflammation of the HES goats could be attributed to these bacterial species and functions [53].
Besides, the present study also revealed the roles of colonic fungal community and fungal functions in affecting the colonic mucosal immunologic state differences between LES and HES goats.Our study showed that the Hypoxylon sp.EC38, Fusarium vanettenii, Fusarium sp.AF-8, Neocallimastix californiae, and fungal functions (ko05203, ko05034, and ko05322) had positive correlations with immune-related DEGs that includes TRAV, KIR3DL, and c-MAF for antigen presentation; c-MAF, GATA3, and IL-5 for T H 2-mediated inflammatory process; as well as A1AT1 and C3 for complement and coagulation cascades.Of these fungi, Fusarium vanettenii and Fusarium sp.AF-8 both belong to the genus of Fusarium, which was widely suggested to produce the Fusarium mycotoxins that includes the Zearalenone [53].Hence, the increased potential pathogenic roles of these fungi under the HES diet may be attributed to their produced Fusarium mycotoxins, which were also widely suggested to induce the gut inflammation through activating T H 2-mediated inflammatory process [53,53].
Last, in the present study, we evaluated the proportions of variation in host phenotypes that explained by colonic microbiota, microbial functions, and luminal metabolome [53].The contributions of the colonic bacteria and bacterial functions were higher than the fungi and fungal functions on the fecal scores and MUC2, indicating that the colonic bacteria occupied a dominant position in the connection between environmental factors and host responses in the colon.The omics-explainability of the colonic luminal metabolome in our study suggests that the colonic metabolism makes greater contributions to MUC2 and pathologic scores compared with the colonic microbiota and functions [53].Therefore, fungi and metabolites should be incorporated into the future systemic studies aimed at hindgut dysbiosis in young ruminants.Our data could help in advancing our understanding of the interactions between the microbiome and host in the hindgut of young ruminants.However, we acknowledge that some limitations should be noted.The present study exhibits that the HES diet had a comprehensive effect of multiple factors in mediating colonic inflammation of goats.Investigating the causality of HES diet and hindgut dysbiosis should be continually conducted in the future even though our findings cannot be evidence of causation.Further systematic research, e.g., metagenome assembled genomes, proteome, or other multi-omics research, is needed to complement our findings due to the few studies that have investigated the influence of hindgut homeostasis on intestinal functions and health in young ruminants.

Conclusions
This is one of the most comprehensive studies with the aim of elucidating the molecular mechanisms behind the colonic dysbiosis in response to HES diet.In summary, our study demonstrated that HES caused BA accumulation, and weakened host mucosal MUC2 biosynthesis and epithelial tight junction, resulting in luminal macromolecule breaching the physical barrier.Colonic microbiome and its metabolites stimulated colonic inflammation and tissue injury via promoting the antigen presentation and promoting T H 2-mediated inflammatory process.Furthermore, the capacity of colonic water absorption was also repressed (Fig. 7).The present study provides systemic insight into the disrupted colonic homeostasis in a goat model.The nonnegligible contribution of fungi is worth more future studies in hindgut microbehost interactions.It will be instrumental in developing new nutritional strategies to relieve the overload starchmediated gut dysfunction and hindgut dysbiosis in young ruminants.

Fig. 1
Fig.1Host phenotypes and related gene expression of two groups of goats fed with low enteral starch (LES) and high enteral starch (HES) diets.A Fecal scores (n = 20) and gene expression of water absorption and ion transport (n = 8) of two groups.B Representative H&E-stained colonic sections and colonic histopathology scores of two groups (n = 8).C Representative colonic MUC2 immunofluorescence and IHC sections with semi-quantitative analysis of MUC2 biosynthesis in colonic mucosa (%) of two groups (n = 8).For immunofluorescence sections, samples were stained with DAPI (blue) and MUC2 antibody (green).For IHC sections, brown staining shows MUC2 expression, hematoxylin (blue) was used as a counterstain for nuclei.D Serum level of LPS and DAO between two groups of growing goats.The difference between two groups was identified by independent sample T-test.Symbols indicate significance ( * * , P < 0.01; * , P < 0.05; #, 0.05 < P < 0.1)

Fig. 4
Fig. 4 Effects of enteral starch content on colonic mucosal barrier and transcriptional profile (n = 8 for mucosal defense parameters, and n = 6 for host transcriptome).A Gene expression of colonic mucosal defense of two groups of growing goats.The difference between two groups was identified by independent sample T-test.B Heatmap of colonic mucosal DEGs.Compared with LES, 64 DEGs were upregulated and 81 DEGs were downregulated significantly in HES goats, respectively.C KEGG pathway enrichment analysis of DEGs (top10).D DEGs related to differential pathways in the colon of HES goats (ko04610, ko04612, and ko05321).ko04610, Complement and coagulation cascades; ko04612, Antigen processing and presentation; ko05321, Inflammatory bowel disease (IBD).The upregulated DEGs are shown in red, whereas the downregulated DEGs are presented in blue.E DEGs (RNA-Seq) related to the mucosal immunologic pathways of the HES goats compared with the LES goats (the difference was identified by Mann-Whitney U test.).qPCR validation of DEGs in the colonic mucosa of HES goats compared with the LES goats (the difference was identified by independent sample T-test).The values are presented as log 2 (fold change).F Effects of enteral starch on the gene expression which involved in colonic mucosal T helper (T H ) cells' function.The difference between two groups was identified by independent sample T-test.Symbols indicate significance ( * * , P < 0.01; * , P < 0.05; #, 0.05 < P < 0.1) The proportion of variance in host phenotypes explained by the colonic microbiota, microbial functions, and metabolome.Only omics-based explainability ≥ 10% were displayed and IL-13 in the colonic mucosa of HES goats confirmed that the colonic inflammation in growing goats caused by HES was attributed to the T H 2 cytokine pattern, which might be related to the T H 2-mediated immune disorder.On this basis, colon microbiome and their metabolome in affecting colonic mucosal immunologic state were studied in the present study.The Bacteroides sp.
CAG:927, Prevotella sp.P4-67, and Prevotella sp.PINT were enriched in HES group.Bacteroides sp.CAG:927 can cause endogenous infection during the disruption of intestinal microbiota.Prevotella is saccharolytic or moderately saccharolytic that promotes the production of acetate, succinate, and a few isobutyrate, isovalerate, and lactate