Abstract
Introduction:
Feeding low protein (LP) diet to animals impose severe challenge to animals' immune homeostasis. However, limited knowledge about the underlying adaption mechanism of host and ruminal microbiota responding to LP diet were well understood. Herein, this study was performed to examine the changes in relative abundance of ruminal microbiota and host ruminal mucosal transcriptome profiles in response to a LP diet.
Methods:
A total of twenty-four female Xiangdong balck goats with similar weight (20.64 ± 2.40 kg) and age (8 ± 0.3 months) were randomly assigned into two groups, LP (5.52% crude protein containing diet) and CON (10.77% crude protein containing diet) groups. Upon completion of the trial, all goats were slaughtered after a 16-hour fasting period in LiuYang city (N 28°15′, E 113°63′) in China. HE staining, free amino acids measurement, transcriptome analysis and microbiome analysis were applied to detect the morphology alterations, free amino acids profile alterations and the shift in host ruminal mucosal transcriptome and ruminal microbiota communities.
Results:
Firstly, the results showed that feeding LP diet to goats decreased the rumen papilla width (P = 0.043), surface area (P = 0.013) and total ruminal free amino acids concentration (P = 0.016). Secondly, microbiome analysis indicated that 9 microbial genera, including Eubacterium and Prevotella, were enriched in LP group while 11 microbial genera, including Butyrivibrio and Ruminococcus, were enriched in CON group. Finally, in terms of immune-related genes, the expression levels of genes involved in tight junction categories (e.g., MYH11, PPP2R2C, and MYL9) and acquired immunity (e.g., PCP4 and CXCL13) were observed to be upregulated in the LP group when compared to the CON group.
Conclusion:
Under the LP diet, the rumen exhibited increased relative abundance of pathogenic microbiota and VFA-degrading microbiota, leading to disruptions in immune homeostasis within the host's ruminal mucosa. These findings indicate that the ruminal microbiota interacts with host results in the disruption in animals' immune homeostasis under LP diet challenge.
1. Introduction
Feeding low-protein (LP) diets to goats represent a common strategy for reducing the cost of breeding and nitrogen losses to the environment (Gebeyew et al., ). However, as the amino acids play a vital role in maintaining the procedure of animal growth and act as the basal component of the immunity system of the body, low-protein content in the feed imposes negative effects on ruminants' growth performance and homeostasis of the gastrointestinal tract (GIT) (Wang et al., ). Furthermore, inadequate protein supply for animals is closely linked to metabolic and immune disorders, such as inflammatory bowel diseases and cancer development (Lan et al., ; Soares et al., ).
Numerous microorganisms colonize the GIT of mammals and play vital roles in the procedures of extracting energy from diets, regulating the homeostasis of the immune system and producing hormones to regulate a series of body activities (Macpherson et al., ; Shanahan et al., ). Various factors, including diet composition, host genetics, and antibiotics, have been reported to impose significant effects on the diversity of microbiota in the GIT, and finally may impose severe challenges to host health (Gilbert et al., ). In humans and mice, a lower protein diet decreased the abundance of Bifidobacterium in gut digesta, which produced antimicrobial peptides to protect the gut from pathogenic bacteria infection and decreased the risk of gut inflammation (Liu et al., ; Feng et al., ). In addition to the microbiota residing in the digesta, the microbiota colonizing at the mucosal border is also vital to maintain host metabolic and immune homeostasis (Jiao et al., ). Despite the significance, its interactions with the host are still not well interpreted.
As the rumen represents the first site for interplay among ingested nutrients, host GIT, and the microbiota, maintaining metabolic, and immune homeostasis in the rumen is of great importance to ruminants. In lambs and beef steers, limited protein supply caused drastic alterations in bacteria diversity and fermentation parameters in the rumen digesta (Wang et al., ; Lv et al., ). Thus, we hypothesized that the ruminal mucosal microbiota of low-protein diets would be different from that of control diets, and such a variation could manipulate the molecular adaptation of the rumen. To validate this, we used a combination of 16S rRNA sequencing and mucosal transcriptome analysis to dissect bacterial diversity and the host responses in the rumen mucosa.
2. Materials and methods
2.1. Animals, housing, and experimental design
All procedures for animal experimentation were performed according to the protocol ISA-2019-0115 approved by the Animal Care Committee, Institute of Subtropical Agriculture, Chinese Academy of Sciences, Changsha, Hunan, China.
A total of 24 female Xiangdong black goats with similar weight (20.6 ± 2.4 kg) and age (8.0 ± 0.3 months) were used in this experiment (Gebeyew et al., ). These goats were randomly assigned into two groups, a control (CON) group and a low-protein (LP) group. The experimental diets were designed according to the feeding standards of goats in China (NY/T 816-2004) and met 1.3 times the maintenance requirement of metabolic energy (ME) on the basis of previous studies (Tang et al., ), and consisted of 70% of rice straw and 30% of concentrate. The detailed ingredients of experimental diets are shown in Supplementary Table 1. The low-protein diet was formulated to contain half of the CP content of the control diet, and this was achieved by replacing soybean meal with corn meal. The total feeding trial period was 70 days, including 25 days of adaptation and 45 days of experiments. These goats were maintained individually in metabolic cages with access to fresh water offered in equal portions at 07:00 and 17:00 h daily.
2.2. Sample collection
After the trial, all goats were slaughtered following a 16-h fasting period. Immediately after all the goats were killed by bleeding the jugular vein by a registered veterinarian in Liuyang city (N 28°15′, E 113°63′), China, ruminal tissues were collected from the ventral sac of the rumen and washed by cold PBS for three times. Afterward, the mucosa samples were separated from the underlying muscular layer, cut into 1 × 1 cm fragments, frozen in liquid nitrogen, and stored at −80°C for subsequent DNA extraction, transcriptome analysis, and free amino acid profile measurements. Moreover, two 2 × 2 cm mucosa samples from the ventral sac of the rumen were collected and fixated in formaldehyde for further HE staining analysis.
2.3. HE staining
After fixation in formaldehyde for 24 h, ruminal mucosa was embedded in paraffin wax (Sangon, Shanghai, China) and then successively sliced into 10 pieces of 4–5 μm thick sections followed by mounting onto poly-L-lysine-coated glass slides (Hailun, Changsha, China). The slides were stained with hematoxylin for 5 min followed by immersing the slides into hydrochloric acid and ammonia water for 10 s separately. After washing with running water for 1 h, all the slides were stained with eosin for 5 min. All the images were acquired under a 10× objective lens and a 10× ocular lens by a fluorescence microscope (Olympus, Tokyo, Japan) equipped with DP2-BSW software.
2.4. Free amino acid profile
The free amino acid profiles of ruminal mucosa were determined using the procedures described in previous studies (Wu et al., ). In brief, 0.4 g of ruminal mucosa was weighed and ground into powder with liquid nitrogen. The free amino acids were extracted from the mucosa by adding 2 ml of 8% sulfosalicylic acid to the mucosa powder. After quiescence at 4°C for 12 h, the samples were centrifuged at 8,000 × g at 4°C for 10 min, and the supernatant fluid was collected and filtered through a 0.22-μm membrane. The samples were analyzed by an automatic amino acid analyzer (L-8900; Hitachi Global Inc., Hitachi, Japan).
2.5. Transcriptome analysis
Total RNA was extracted from ruminal mucosa following the instructions of the TaKaRa MiniBEST Universal RNA Extraction Kit (TaKaRa, Dalian, China; code no. 9767). The extracted RNA was reverse-transcribed to build the cDNA library using an mRNA-Seq Sample Preparation Kit (Illumina, San Diego, USA). After verifying and quantifying the cDNA, the libraries were sequenced on the Illumina HiSeq 4000 platform.
The procedures for filtering raw sequencing data were detailed in previous studies (Jiao et al., ), and the cleaned data were stored in FASTA format. HISAT2 (v2.04) was applied to map the clean reads to the reference genome (Kim et al., ). The clean reads were aligned to the reference coding gene set by Bowtie (Langmead and Salzberg, ), and the expression level was measured by RSEM (Li and Dewey, ). The analysis of differentially expressed genes (DEGs) was carried out using DESeq2, using thresholds with false discovery rate (FDR) <0.05 and an absolute value of fold change >1.2.
2.6. Microbiome analysis
The total DNA of the ruminal mucosa was extracted using the bead-beating method, as detailed in our previous study (Jiao et al., ). The quantity of total DNA was measured by NanoDrop ND1000 (NanoDrop Technologies, Inc., Wilmington, DE, USA), and the quality of total DNA was observed by gel electrophoresis. Four samples (1/CON, 3/LP) exhibited significant degradation, failing to meet the standards for library sequencing. After amplifying the V3-V4 region of 16S rRNA (341F, ACTCCTACGGGAGGCAGCAG and 806R, GGACTACHVGGGTWTCTAAT) and running an agarose gel, the bands were purified with a QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), which was prior to sequencing by the Illumina MiSeq PE250 platform.
The QIIME (Quantitative Insights into Microbial Ecology) pipeline was utilized to perform quality control on the raw data, with a default quality threshold of 20 (Zhang et al., ). The quality control process involved a sliding window approach to identify low-quality reads alongside a length screening where reads longer than 200 bases were truncated to 200 bases. Afterward, the short reads were assembled into tags using FLASH, and tags were clustered into amplicon sequence variants (ASVs) of 99% similarity using the unoise3 command implemented in USEARCH (Edgar, ). The representative sequences were kept with a total sequencing count >3, and the ASVs were retained that were present in at least two sample replicates to mitigate incidental factors. Taxonomic assignments were performed against the RDP database with a 0.80 confidence threshold (Wang Q. et al., ). Alpha and beta diversities were analyzed using the QIIME pipeline. The Bray–Curtis similarity index was used to calculate the distance matrix for principal coordinate analysis (PCoA). The significance of grouping in the PCoA plots was tested by analysis of dissimilarity (ADONIS) with 999 permutations. To identify the differential microbial communities in the ruminal mucosa between these two groups, a linear discriminate analysis (LDA) effect size (LEfSe) method (Segata et al., ) was applied with an LDA threshold value of 2.0 using the R software (version 4.0.5).
2.7. Statistical analysis
Data were carried out with a one-way analysis of variance (ANOVA) using the R software (version 4.0.5). Student's t-test was applied to determine the effects of a low-protein diet on the ruminal mucosal morphology and free amino acid profiles, with a P-value of < 0.05 considered significant.
3. Results
3.1. Morphological analysis
As shown in Table 1, feeding LP diet to goats significantly decreased ruminal papilla width (P = 0.043) and surface area (P = 0.013), while the length of ruminal papilla showed no significant difference (P = 0.268) between these two groups.
Table 1
| Items | Treatments | SEM | P-value | |
|---|---|---|---|---|
| CON | LP | |||
| Papilla length (cm) | 1.18 | 1.16 | 0.38 | 0.268 |
| Papilla width (cm) | 0.48 | 0.40 | 0.09 | 0.043 |
| Papilla area (cm2) | 0.60 | 0.47 | 0.17 | 0.013 |
Morphological changes of the rumen papilla in goats fed with LP vs. CON diet.
LP, low-protein group; CON, control group.
3.2. Free amino acid profile
Feeding LP diet to goats significantly decreased the total amino acid concentration (P = 0.016). The concentration of glutamic acid (P = 0.011) and aspartic acid (P = 0.037) was significantly decreased in the ruminal mucosa under the LP diet. Similarly, feeding LP diet to goats significantly decreased the concentration of isoleucine (P = 0.033) and valine (P = 0.040). The concentration of leucine tended to decrease in the ruminal mucosa under the LP diet (P = 0.061) (Table 2).
Table 2
| Items (μg/g) | Treatments | SEM | p-value | |
|---|---|---|---|---|
| CON | LP | |||
| Aspartic acid | 46.52 | 36.65 | 12.12 | 0.037 |
| Threonine | 123.04 | 109.20 | 25.36 | 0.121 |
| Serine | 13.37 | 12.63 | 4.21 | 0.356 |
| Glutamic acid | 414.11 | 328.99 | 83.74 | 0.011 |
| Glycine | 243.55 | 223.06 | 39.91 | 0.136 |
| Alanine | 49.74 | 41.79 | 9.361 | 0.030 |
| Cysteine | 2.13 | 2.53 | 1.43 | 0.277 |
| Valine | 22.06 | 18.21 | 5.39 | 0.040 |
| Methionine | 3.93 | 3.78 | 1.35 | 0.408 |
| Isoleucine | 8.71 | 7.04 | 2.01 | 0.033 |
| Leucine | 14.17 | 12.07 | 2.96 | 0.061 |
| Tyrosine | 11.33 | 10.75 | 1.85 | 0.25 |
| Phenylalanine | 10.14 | 8.42 | 1.85 | 0.252 |
| Lysine | 16.18 | 13.44 | 4.21 | 0.079 |
| Histidine | 1.76 | 9.66 | 1.55 | 0.064 |
| Arginine | 8.73 | 6.75 | 3.74 | 0.128 |
| Proline | 23.09 | 21.08 | 5.89 | 0.230 |
| NEAA | 809.98 | 644.87 | 150.64 | 0.005 |
| EAA | 211.63 | 184.81 | 36.54 | 0.044 |
| BCAA | 42.18 | 35.44 | 9.59 | 0.06 |
| TAA | 1021.61 | 866.07 | 161.35 | 0.016 |
Effects of different dietary protein levels on rumen mucosa and free amino acid profiles in goats.
TAA, total amino acids; BCAA, branch chain amino acids; EAA, essential amino acids; NEAA, non-essential amino acids; LP, low-protein group; CON, control group.
3.3. Transcriptome profile and functional analysis
The transcriptome results showed that 41 differentially expressed (DE) genes were observed between the LP and CON groups, of which the expression of 28 genes was upregulated, while 13 genes were downregulated in the LP group in comparison with the CON group. The majority of DEGs were related to immune and muscle categories. For immune-related genes, the expression of genes involved in the categories of tight junction such as MYH11, PPP2R2C, and MYL9, cancer development such as TAGLN and CALD1, and acquired immunity such as PCP4 and CXCL13 was upregulated in the LP group when compared with that in the CON group (Table 3). For muscle-related genes, the expression of genes involved in the categories of muscle contraction such as LMOD1, SYNM, CNN1, TPM2, MYH11, MYL9, MYLK, and TPM1, muscle cell proliferation such as PDLIM3 and RBPMS2, and muscle development such as LDB3 and DMPK was upregulated in the LP group compared with that in the CON group (Table 4). The expression of gene encoding desmin, which acts as a biomarker of muscle, was downregulated in the LP group when compared with that in the CON group.
Table 3
| Gene symbol | Gene name | Category | logFC (LP vs. CON) | Ln (p-value) |
|---|---|---|---|---|
| MYH11 | myosin-11 isoform X2 | Tight junction cancer development | 2.16 | −97.6 |
| PPP2R2C | Serine/threonine-protein phosphatase 2A 55 kDa regulatory subunit B gamma isoform | Tight junction | −1.88 | −57.5 |
| MYL9 | Myosin regulatory light polypeptide 9 isoform X1 | Tight junction | 2.11 | −91.6 |
| MYLK | Myosin light chain kinase, smooth muscle | Organ damage | 2.26 | −106.27 |
| PCP4 | Purkinje cell protein | Anti-virus | 1.46 | −28.2 |
| CALD1 | Caldesmon isoform X1 | Inflammation Stomach cancer | 1.36 | −40.3 |
| CXCL13 | C-X-C motif chemokine 13 | Acquired immunity (T cell, B cell) | 1.46 | −38.3 |
| TAGLN | Transgelin | Cancer development inflammation | 1.58 | −52.4 |
Expression profile of DE immune-related genes.
CON, control group; LP, low-protein group.
Table 4
| Gene symbol | Gene name | category | logFC (LP vs. CON) | Ln (p-value) |
|---|---|---|---|---|
| DES | Desmin | Muscle contraction | −1.57 | −252.1 |
| LMOD1 | Leiomodin-1 | Muscle contraction Muscle cell proliferation | 1.31 | −33.8 |
| SYNM | Synemin isoform X1 | Muscle contraction | 2.47 | −46.5 |
| CNN1 | Calponin-1 | Muscle contraction Actin binding | 1.97 | −37.2 |
| TPM2 | Tropomyosin beta chain | Muscle contraction Actin binding | 1.17 | −132.2 |
| MYH11 | Myosin-11 isoform X2 | Muscle contraction Actin binding | 2.16 | −97.6 |
| MYL9 | Myosin regulatory light polypeptide 9 isoform X1 | Muscle contraction Actin binding | 2.11 | −91.6 |
| MYLK | myosin light chain kinase, smooth muscle | Muscle contraction | 2.26 | −106.27 |
| TPM1 | Tropomyosin alpha-1 chain isoform X9 | Muscle contraction Actin binding | 1.15 | −81.1 |
| PDLIM3 | PDZ and LIM domain protein 3 isoform X1 | Muscle cell proliferation | 1.21 | −28.9 |
| RBPMS2 | RNA-binding protein with multiple splicing 2 | Muscle cell proliferation | 2.31 | −8.97 |
| DMPK | Myotonin-protein kinase isoform X1 | Muscle cell survival | 1.52 | −39.0 |
| LDB3 | LIM domain-binding protein 3 isoform X1 | Muscle development | 1.87 | −34.3 |
Expression profile of DE muscle-related genes.
CON, control group; LP, low-protein group.
3.4. Microbial diversity and composition
The alpha diversity analysis in ruminal mucosa showed lower ACE (P = 0.013), Chao1 (P = 0.015), and Shannon (P = 0.031) diversity indexes in the LP group than in the CON group, while the Simpson index showed no significant differences (P = 0.098, Table 5). Beta diversity analysis showed a clear distance separation between the LP and CON groups, and the major variance observed between these two groups was 14.8% (Figure 1). To explore these findings, we performed a LEfSe analysis to identify the differential microbial taxa between the LP and CON groups. The results showed that 9 microbial genera, namely, Intestinimonas, Atopobium, Sphingomonas, Fretibacterium, Eubacterium, Treponema, Prevotella, Desulfobulbus, and Campylobacter were enriched in the LP group while 11 microbial genera, namely, Desulfovibrio, Fibrobacter, Succiniclastium, Saccharofermentans, Butyrivibrio, Ruminococcus, Ruminobacter, Succinimonas, Clostridium-XIVa, Pseudobutyrivibrio, and Pyramidobacter were enriched in the CON group (Figure 2, LDA score >2).
Table 5
| Items | Treatments | SEM | p-value | |
|---|---|---|---|---|
| CON | LP | |||
| Chao1 | 6,460.43 | 5,584.84 | 791.57 | 0.013 |
| ACE | 6,522.94 | 5,662.03 | 820.71 | 0.015 |
| Shannon | 6.81 | 6.35 | 0.52 | 0.031 |
| Simpson | 1.00 | 0.99 | 0.003 | 0.098 |
| Coverage | 0.96 | 0.97 | 0.005 | 0.006 |
Alpha diversity analysis of ruminal mucosal microbiota in goats fed with LP vs. CON diets.
LP, low-protein group; CON, control group.
Figure 1
Figure 2
3.5. Interactions between bacterial biomarkers and DEGs
We investigated the correlations between differentially expressed genes (DEGs) and the prominent microbes identified by the LEfSe analysis, in order to gain insight into the interaction between the host and microbiota in the rumen (Figure 3). Using Spearman correlations, we identified a strong negative association between the majority of immune and muscle-related DEGs and the prominent microbes, such as Prevotella and Eubacterium.
Figure 3
4. Discussion
It is usually considered that the protein ingested from the diet is converted into NH3-N in the rumen and then utilized by the ruminal microbiota to synthesize microbial protein (MCP). In addition to MCP, peptides and free amino acids digested from the dietary protein and synthesized by the ruminal microbiota (Russell et al., ) could be absorbed in the rumen mucosa (Poole et al., ). Not surprisingly, low-protein content in the diet induced a decline in the ruminal mucosal free amino acid concentrations, especially branched-chain amino acids (BCAAs, such as valine and isoleucine) and non-essential amino acids (NEAAs, such as glutamic acid and aspartic acid). It has been reported that the scanty of glutamate and BCAAs could induce the suppression of cell proliferation and protein synthesis (Sancak et al., ; Li et al., ). Hence, it is reasonable to observe a decline in the morphology of the rumen papilla.
Another notable change in response to the LP diet was genes involved in the muscle, which was likely the result of cells from the muscle layer that remained adhered to the dissected mucosal samples. Previous studies indicated that feeding a low-protein diet to mammals induces a decline in muscle layer thickness and muscle destruction in the GIT of animals (Peng et al., ). Although TMP1, TMP2, DMPK, PDLIM3, RBPMS2, and LDB3, which were reported to engage in muscle repairment (Jauvin et al., ), muscle cell proliferation (Notarnicola et al., ; Yin et al., ), and muscle regeneration (Knight et al., ), were upregulated in LP diet, the upregulation of MYLK, MYH11, and MYL9 may induce muscle relaxation and GIT movement retention (Iwasaki et al., ) and finally impose negative effects on nutrient digestion and immune homeostasis (Kashyap et al., ). Muscle repairment and muscle regeneration are closely related processes. Muscle repairment refers to the physiological mechanisms involved in restoring damaged muscle tissue, while muscle regeneration refers to the process by which new muscle cells are formed to replace damaged or lost muscle tissue (Notarnicola et al., ; Yin et al., ). Movement retention is linked to both muscle repairment and muscle regeneration. When muscles are properly repaired and regenerated, movement retention is improved, allowing individuals to maintain their ability to perform various physical activities (Iwasaki et al., ). Additionally, the decline in muscle layer thickness and muscle destruction in the rumen may account for the downregulation of DES, which was considered the major component of muscle fibers in the LP group (Cizkova et al., ). The LMOD1, which promoted muscle relaxation (Nanda et al., ), was also upregulated in the ruminal mucosa and involved in the procedures of rumen contraction shunt.
In addition to the function of muscle contraction, the MYH11 also plays vital roles in a series of biological processes, such as intracellular signaling transduction and cell adhesion (Derycke et al., ; Bowers et al., ). It has been reported that the promotion of MYH11 expression was related to cell mobility and would cause the disruption of organ integrity and barrier function (Poninska et al., ). As anticipated, we observed a decline in the width and surface area of the ruminal papilla. Synchronously, LP inclusion resulted in the upregulation of PCP4, which played a vital role in the procedure of pathological microbiota infection resistance (Ge, ; Sui et al., ). It is not surprising to find that CALD1 and TAGLN, which promoted inflammation, were upregulated in the LP group, implying an elevated inflammation risk in the GIT (Shen et al., ; Liu et al., ). Notably, the procedures of antigen presentation to immune cells were upregulated via the upregulation of CXCL13 (Ohmatsu et al., ; Wang X. B. et al., ), to enhance the immune system and finally alleviate the negative effect of limited protein supply for goats.
Feeding LP diet to goats decreased the ruminal microbiota diversity, which was in line with previous reports in beef steers (Wang et al., ). It is noteworthy that several pathogenic bacteria genera propagated in ruminal mucosa during LP inclusion, namely, Fretibacterium, Treponema, Intestinimonas, Sphingomonas, and Campylobacter (Brooks et al., ; Gubert et al., ; Fan and Pedersen, ). These pathogens possessed the capacity to produce toxic substances such as endotoxin in the GIT (Zhao et al., ; Fan and Pedersen, ), thereby aggravating the above-mentioned disruption of immune homeostasis in the ruminal mucosa. In addition, the higher proportion of starch as a primary carbohydrate source in the LP group may provide a more suitable growth and metabolic environment for Fretibacterium and Campylobacter. This could potentially lead to a relatively higher abundance of Fretibacterium in the LP group.
Additionally, the results of our analysis revealed a significant negative association between DEGs related to immune response and muscle function and the prominent microbes, such as Prevotella and Eubacterium, in the rumen, suggesting that these microbial taxa may have a suppressive effect on the immune system and muscle development of the host. Meanwhile, the surge of H2S-producing Desulfobulbus in the mucosal microbiota during LP intervention (Cabrera et al., ) might also contribute to the disruption of immune stability in the ruminal mucosa. Intriguingly, LP inclusion declined the abundances of generally accepted carbohydrate-degrading bacteria in the ruminal mucosa, namely, fiber-degrading Fibrobacter and Ruminococcus (Yeoman et al., ), acetate producers Desulfovibrio and Saccharofermentans (Rettenmaier et al., ), succinate producers Succinimonas and Pyramidobacter (Gilbert et al., ), propionate producers Succiniclastium, and butyrate producers Butyrivibrio and Pseudobutyrivibrio (Pidcock et al., ). These highlighted the significance of ruminal microbiota in short-chain fatty acid (SCFA) production for energy supply and indicated a decline of carbohydrate fermentation in ruminal mucosa during LP intervention. Since the SCFAs not only act as the energy source for ruminants but also act as vital signaling molecules in maintaining immune homeostasis (Yao et al., ), the decline in the abundance of VFA-producing microbiota in the LP diet may also contribute to the host mucosal immune disruption.
5. Conclusion
Feeding LP diets to goats significantly decreased NEAA and BCAA concentration in the rumen mucosa, restrained the ruminal papilla growth, altered the gene expression related to muscle contraction, and depressed the gene expression related to immune homeostasis. Moreover, the low-protein diet also induced a decline in the relative abundance of carbohydrate-degrading microbiota while the promotion induced the abundance of pathogenic microbiota in the ruminal mucosa. Taken together, feeding LP diets to goats imposed severe challenges to goats as a result of the combination of pathogenic microbiota enrichment and the disruption of ruminal immune homeostasis.
Statements
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 below: https://ngdc.cncb.ac.cn/gsa—CRA011528 and CRA011517.
Ethics statement
The animal studies were approved by the Animal Care Committee, Institute of Subtropical Agriculture, Chinese Academy of Sciences. 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
JW and ZT designed the research. CT and JJ performed the research and analyzed the samples. CZ and CT contributed intellectually to the analysis and interpretation of the data. CT, JW, and ZT wrote the manuscript. All authors read and approved the final version of the manuscript.
Funding
This work was supported by grants from Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDA26040304) and the National Natural Science Foundation of China (No. 31730092).
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.2023.1237955/full#supplementary-material
Abbreviations
AAs, amino acids; BCAAs, branch-chain amino acids; CALD1, Caldesmon; CXCL13, C-X-C motif chemokine 13; CNN1, calponin-1; DEG, differential expression genes; FDMPK, Myotonin-protein kinase isoform X1; DSG1, Desmoglein-1-alpha; EAA, essential amino acids; GIT, gastrointestinal tract; LDB3, LIM domain-binding protein 3 isoform X1; MCP, Microbial protein; MYH11, Myosin-11 isoform X2; MYL9, Myosin regulatory light polypeptide 9 isoform X1; MYLK, Myosin light chain kinase, SM, smooth muscle; MYH11, Myosin-11; MYL9, Myosin regulatory light polypeptide 9; PCP4, Purkinje cell protein; PDLIM3, PDZ and LIM domain protein 3 isoform X1; PPP1R1A, Protein phosphatase 1 regulatory subunit 1A; PPP1R12B, Protein phosphatase 1 regulatory subunit 12B; PPP2R2C, serine/threonine-protein phosphatase 2A 55 kDa regulatory subunit B gamma isoform; RBPMS2, RNA-binding protein with multiple splicing 2; SYNM, synemin isoform X1; TAGLN, Transgelin; TPM1, Tropomyosin alpha-1 chain; TAA, total amino acids; VFA, volatile fatty acids.
References
1
BowersR. R.ManevichY.TownsendD. M.TewK. D. (2012). Sulfiredoxin redox-sensitive interaction with S100A4 and non-muscle myosin IIA regulates cancer cell motility. Biochemistry51, 7740–7754. 10.1021/bi301006w
2
BrooksP. T.BrakelK. A.BellJ. A.BejcekC. E.GilpinT.BrudvigJ. M.et al. (2017). Transplanted human fecal microbiota enhanced Guillain Barre syndrome autoantibody responses after Campylobacter jejuni infection in C57BL/6 mice. Microbiome5, 1–22. 10.1186/s40168-017-0284-4
3
CabreraG.PerezR.GomezJ. M.AbalosA.CanteroD. (2006). Toxic effects of dissolved heavy metals on Desulfovibrio vulgaris and Desulfovibrio sp strains. J. Hazar. Mater.135, 40–46. 10.1016/j.jhazmat.2005.11.058
4
CizkovaD.SoukupT.MokryJ. (2009). Expression of nestin, desmin and vimentin in intact and regenerating muscle spindles of rat hind limb skeletal muscles. Histochem. Cell Biol.131, 197–206. 10.1007/s00418-008-0523-7
5
DeryckeL.StoveC.Vercoutter-EdouartA. S.De WeverO.DolleL.ColpaertN.et al. (2011). The role of non-muscle myosin IIA in aggregation and invasion of human MCF-7 breast cancer cells. Int. J. Dev. Biol.55, 835–840. 10.1387/ijdb.113336ld
6
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics26, 2460–2461. 10.1093/bioinformatics/btq461
7
FanY.PedersenO. (2021). Gut microbiota in human metabolic health and disease. Nat. Rev. Microbiol.19, 55–71. 10.1038/s41579-020-0433-9
8
FengY. Q.DuanY. F.XuZ. J.LyuN.LiuF.LiangS. H.et al. (2019). An examination of data from the American Gut Project reveals that the dominance of the genus Bifidobacterium is associated with the diversity and robustness of the gut microbiota. Microbiol. Open8, 939. 10.1002/mbo3.939
9
GeL. (2003). Identification of the delta-6 desaturase of human sebaceous glands: expression and enzyme activity. J. Invest. Dermatol.121, 434–434. 10.1046/j.1523-1747.2003.12123.x
10
GebeyewK.ChenW. X.YanQ. X.HeZ. X.TanZ. L. (2021). Growth of pancreas and intestinal enzyme activities in growing goats: influence of a low-protein diet. Agric. Basel11, 1155. 10.3390/agriculture11111155
11
GilbertJ. A.BlaserM. J.CaporasoJ. G.JanssonJ. K.LynchS. V.KnightR. (2018). Current understanding of the human microbiome. Nat. Med.24, 392–400. 10.1038/nm.4517
12
GubertC.KongG.RenoirT.HannanA. J. (2020). Exercise, diet and stress as modulators of gut microbiota: Implications for neurodegenerative diseases. Neurobiol. Dis.134, 104621. 10.1016/j.nbd.2019.104621
13
IwasakiT.Murata-HoriM.IshitobiS.HosoyaH. (2001). Diphosphorylated MRLC is required for organization of stress fibers in interphase cells and the contractile ring in dividing cells. Cell Struct. Funct.26, 677–683. 10.1247/csf.26.677
14
JauvinD.ChretienJ.PandeyS. K.MartineauL.RevillodL.BassezG.et al. (2017). Targeting DMPK with antisense oligonucleotide improves muscle strength in myotonic dystrophy type 1 mice. Molec. Ther. Nucl. Acids7, 465–474. 10.1016/j.omtn.2017.05.007
15
JiaoJ.WuJ.WangM.ZhouC.ZhongR.TanZ. (2018). Rhubarb supplementation promotes intestinal mucosal innate immune homeostasis through modulating intestinal epithelial microbiota in goat kids. J. Agric. Food Chem.66, 1047–1057. 10.1021/acs.jafc.7b05297
16
JiaoJ.ZhangX.WangM.ZhouC.YanQ.TanZ. (2019). Linkages between epithelial microbiota and host transcriptome in the ileum during high-grain challenges: implications for gut homeostasis in goats. J. Agric. Food Chem.67, 551–561. 10.1021/acs.jafc.8b05591
17
JiaoJ. Z.HuangJ. Y.ZhouC. S.TanZ. L. (2015). Taxonomic identification of ruminal epithelial bacterial diversity during rumen development in goats. Appl. Environ. Microbiol.81, 3502–3509. 10.1128/AEM.00203-15
18
KashyapP. C.MarcobalA.UrsellL. K.LaraucheM.DubocH.EarleK. A.et al. (2013). Complex interactions among diet, gastrointestinal transit, and gut microbiota in humanized mice. Gastroenterology144, 967–977. 10.1053/j.gastro.2013.01.047
19
KimD.LandmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods12, 357–U121. 10.1038/nmeth.3317
20
KnightR. F.BaderD. M.BackstromJ. R. (2003). Membrane topology of Bves/Pop1A, a cell adhesion molecule that displays dynamic changes in cellular distribution during development. J. Biol. Chem.278, 32872–32879. 10.1074/jbc.M301961200
21
LanA.AndriamihajaM.BlouinJ. M.LiuX.DescatoireV.Desclee de MaredsousC.et al. (2015). High-protein diet differently modifies intestinal goblet cell characteristics and mucosal cytokine expression in ileum and colon. J. Nutr. Biochem.26, 91–98. 10.1016/j.jnutbio.2014.09.007
22
LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods9, 357–U354. 10.1038/nmeth.1923
23
LiB.DeweyC. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform.12, 323. 10.1186/1471-2105-12-323
24
LiX. G.SuiW. G.GaoC. Q.YanH. C.YinY. L.LiH. C.et al. (2016). L-Glutamate deficiency can trigger proliferation inhibition via down regulation of the mTOR/S6K1 pathway in pig intestinal epithelial cells. J. Animal Sci.94, 1541–1549. 10.2527/jas.2015-9432
25
LiuX. X.BlouinJ. M.SantacruzA.LanA.AndriamihajaM.WilkanowiczS.et al. (2014). High-protein diet modifies colonic microbiota and luminal environment but not colonocyte metabolism in the rat model: the increased luminal bulk connection. Am. J. Physiol. Gastrointest. Liver Physiol.307, G459–G470. 10.1152/ajpgi.00400.2013
26
LiuY. X.XieS. H.ZhuK. Y.GuanX. L.GuoL.LuR. Q. (2021). CALD1 is a prognostic biomarker and correlated with immune infiltrates in gastric cancers. Heliyon7, e07257. 10.1016/j.heliyon.2021.e07257
27
LvCuiK.QiM. L.WangS. Q.DiaoQ. Y.ZhangN. F. (2020). Ruminal microbiota and fermentation in response to dietary protein and energy levels in weaned lambs. Animals10, 109. 10.3390/ani10010109
28
MacphersonA. J.de AgueroM. G.Ganal-VonarburgS. C. (2017). How nutrition and the maternal microbiota shape the neonatal immune system. Nat. Rev. Immunol.17, 508–517. 10.1038/nri.2017.58
29
NandaV.WangT.IyerD.PjanicM.LiuB. X.NguyenT.et al. (2017). Functional regulatory mechanism of smooth muscle cell-restricted LMOD1 coronary artery disease locus. Circulation136, A21021–A21021. 10.1161/circ.136.suppl_1.21021
30
NotarnicolaC.RouleauC.Le GuenL.VirsolvyA.RichardS.FaureS.et al. (2012). The RNA-Binding protein RBPMS2 regulates development of gastrointestinal smooth muscle. Gastroenterology143, 687–U204. 10.1053/j.gastro.2012.05.047
31
OhmatsuH.SugayaM.KadonoT.TamakiK. (2007). CXCL13 and CCL21 are expressed in ectopic lymphoid follicles in cutaneous Lymphoproliferative disorders. J. Investigat. Dermatol.127, 2466–2468. 10.1038/sj.jid.5700873
32
PengY.YuK. F.MuC. L.HangS. Q.CheL. Q.ZhuW. Y. (2017). Progressive response of large intestinal bacterial community and fermentation to the stepwise decrease of dietary crude protein level in growing pigs. Appl. Microbiol. Biotechnol.101, 5415–5426. 10.1007/s00253-017-8285-6
33
PidcockS. E.SkvortsovT.SantosF. G.CourtneyS. J.Sui-TingK.CreeveyC. J.et al. (2021). Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate- active enzyme family isoforms. Microbial. Genom.7, 638. 10.1099/mgen.0.000638
34
PoninskaJ. K.BilinskaZ. T.TruszkowskaG.MichalakE.PodgorskaA.Stepien-WojnoM.et al. (2022). Good performance of the criteria of American College of Medical Genetics and Genomics/Association for Molecular Pathology in prediction of pathogenicity of genetic variants causing thoracic aortic aneurysms and dissections. J. Translat. Med20, 42. 10.1186/s12967-022-03251-8
35
PooleC. A.WongE. A.McElroyA. P.VeitH. P.WebbK. E. (2003). Ontogenesis of peptide transport and morphological changes in the ovine gastrointestinal tract. Small Ruminant Res.50, 163–176. 10.1016/S0921-4488(03)00103-2
36
RettenmaierR.ThiemeN.StreubelJ.Di BelloL.KowollikM. L.HuangL. R.et al. (2021). Variimorphobacter saccharofermentans gen. nov., sp. nov., a new member of the family Lachnospiraceae, isolated from a maize- fed biogas fermenter. Int. J. System. Evolut. Microbiol.71, 5044. 10.1099/ijsem.0.005044
37
RussellJ. B.MuckR. E.WeimerP. J. (2009). Quantitative analysis of cellulose degradation and growth of cellulolytic bacteria in the rumen. Fems Microbiol. Ecol.67, 183–197. 10.1111/j.1574-6941.2008.00633.x
38
SancakY.Bar-PeledL.ZoncuR.MarkhardA. L.NadaS.SabatiniD. M. (2010). Ragulator-rag complex targets mTORC1 to the lysosomal surface and is necessary for its activation by amino acids. Cell141, 290–303. 10.1016/j.cell.2010.02.024
39
SegataN.IzardJ.WaldronL.GeversD.MiropolskyL.GarrettW. S.et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol.12, R60. 10.1186/gb-2011-12-6-r60
40
ShanahanF.van SinderenD.O'TooleP. W.StantonC. (2017). Feeding the microbiota: transducer of nutrient signals for the host. Gut66, 1709–1717. 10.1136/gutjnl-2017-313872
41
ShenJ.YangM.JuD.JiangH.ZhengJ. P.XuZ.et al. (2010). Disruption of SM22 promotes inflammation after artery injury via nuclear factor kappaB activation. Circ. Res.106, 1351–1362. 10.1161/CIRCRESAHA.109.213900
42
SoaresM. H.RodriguesG. A.BarbosaL. M. R.Valente JuniorD. T.SantosF. C.RochaG. C.et al. (2020). Effects of crude protein and lactose levels in diets on growth performance, intestinal morphology, and expression of genes related to intestinal integrity and immune system in weaned piglets. Animal Sci. J.91, e13429. 10.1111/asj.13429
43
SuiB. K.ChenD.LiuW.TianB.LvL.PeiJ.et al. (2021). Comparison of lncRNA and mRNA expression in mouse brains infected by a wild-type and a lab-attenuated Rabies lyssavirus. J. General Virol.102, 1538. 10.1099/jgv.0.001538
44
TangS. X.HeY.ZhangP. H.JiaoJ. Z.HanX. F.YanQ. X.et al. (2019). Nutrient digestion, rumen fermentation and performance as ramie (Boehmeria nivea) is increased in the diets of goats. Animal Feed Sci. Technol.247, 15–22. 10.1016/j.anifeedsci.2018.10.013
45
WangC.LiuQ.GuoG.HuoW. J.LiangY.PeiC. X.et al. (2017). Effects of different dietary protein levels and rumen-protected folic acid on ruminal fermentation, degradability, bacterial populations and urinary excretion of purine derivatives in beef steers. J. Agric. Sci.155, 1477–1486. 10.1017/S0021859617000533
46
WangN.WuC. F.QinY.LiuS. M.ZhangR. (2021). Multi-Angle Investigation of the Fractal Characteristics of Nanoscale Pores in the Lower Cambrian Niutitang Shale and Their Implications for CH4 Adsorption. J. Nanosci. Nanotechnol.21, 156–167. 10.1166/jnn.2021.18463
47
WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol.73, 5261–5267. 10.1128/AEM.00062-07
48
WangX. B.YulingH.YanpingJ.XintiT.YaofangY.FengY.et al. (2007). CCL19 and CXCL13 synergistically regulate interaction between B cell acute lymphocytic leukemia CD23(+) CD5(+) B cells and CD8(+) T cells. J. Immunol.179, 2880–2888. 10.4049/jimmunol.179.5.2880
49
WuJ.ZhangX. L.WangR.WangM.HeZ. X.TanZ. L.et al. (2020). Replacing corn grain with corn gluten feed: Effects on the rumen microbial protein synthesis, functional bacterial groups and epithelial amino acid chemosensing in growing goats. Animal Feed Sci. Technol.270, 114684. 10.1016/j.anifeedsci.2020.114684
50
YaoY.CaiX. Y.FeiW. D.YeY. Q.ZhaoM. D.ZhengC. H. (2022). The role of short-chain fatty acids in immunity, inflammation and metabolism. Crit. Rev. Food Sci. Nutr.62, 1–12. 10.1080/10408398.2020.1854675
51
YeomanC. J.FieldsC. J.LepercqP.RuizP.ForanoE.WhiteB. A.et al. (2021). In Vivo competitions between fibrobacter succinogenes, ruminococcus flavefaciens, and ruminoccus albus in a gnotobiotic sheep model revealed by multi-omic analyses. Mbio12, e03533. 10.1128/mBio.03533-20
52
YinH. D.ZhaoJ.HeH. R.ChenY. Q.WangY.LiD. Y.et al. (2020). Gga-miR-3525 TargetsPDLIM3through the MAPK signaling pathway to regulate the proliferation and differentiation of skeletal muscle satellite cells. Int. J. Molec. Sci.21, 5573. 10.3390/ijms21155573
53
ZhangX. L.WuJ.HanX. F.TanZ. L.JiaoJ. Z. (2019). Effects of rumen-protected glucose on ileal microbiota and genes involved in ileal epithelial metabolism and immune homeostasis in transition dairy cows. Animal Feed Sci. Technol.254, 114199. 10.1016/j.anifeedsci.2019.06.003
54
ZhaoF. J.ZhangX. H.LiuS. Q.ZengT. B.YuJ.GuW. M.et al. (2013). Assessment of the immune responses to Treponema pallidum Gpd DNA vaccine adjuvanted with IL-2 and chitosan nanoparticles before and after Treponema pallidum challenge in rabbits. Sci. China-Life Sci.56, 174–180. 10.1007/s11427-012-4434-4
Summary
Keywords
immune homeostasis, low protein diet, microbiome, transcriptome, goat
Citation
Wu J, Tian C, Jiao J, Yan Q, Zhou C and Tan Z (2023) The epithelial transcriptome and mucosal microbiota are altered for goats fed with a low-protein diet. Front. Microbiol. 14:1237955. doi: 10.3389/fmicb.2023.1237955
Received
10 June 2023
Accepted
14 August 2023
Published
04 September 2023
Volume
14 - 2023
Edited by
Benoit St-Pierre, South Dakota State University, United States
Reviewed by
Jun Zhang, Northwest A&F University, China; Yuying Li, Chinese Academy of Agricultural Sciences, China
Updates
Copyright
© 2023 Wu, Tian, Jiao, Yan, Zhou and Tan.
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: Changxin Tian 505632433@qq.com
Disclaimer
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.