Unique rumen micromorphology and microbiota–metabolite interactions: features and strategies for Tibetan sheep adaptation to the plateau

The rumen microbiota—a symbiont to its host and consists of critical functional substances—plays a vital role in the animal body and represents a new perspective in the study of adaptive evolution in animals. This study used Slide Viewer slicing analysis system, gas chromatography, RT-qPCR and other technologies, as well as 16S and metabolomics determination methods, to measure and analyze the microstructure of rumen epithelium, rumen fermentation parameters, rumen transport genes, rumen microbiota and metabolites in Tibetan sheep and Hu sheep. The results indicate that the rumen nipple height and cuticle thickness of Tibetan sheep are significantly greater than those of Hu sheep (p < 0.01) and that the digestion and absorption of forage are greater. The levels of carbohydrate metabolism, lipid metabolism, and protein turnover were increased in Tibetan sheep, which enabled them to ferment efficiently, utilize forage, and absorb metabolic volatile fatty acids (VFAs). Tibetan sheep rumen metabolites are related to immune function and energy metabolism, which regulate rumen growth and development and gastrointestinal homeostasis. Thus, compared with Hu sheep, Tibetan sheep have more rumen papilla and cuticle corneum, and the synergistic effect of the microbiota and its metabolites is a characteristic and strategy for adapting to high-altitude environments.

1 Introduction

Tibetan sheep constitute one of the three major rough-wooled sheep breeds in China (the others include Kazakh sheep and Mongolian sheep), are distributed mainly in the alpine areas of Tibet and Qinghai on the Tibetan Plateau at altitudes of more than 2,000 m, are important livestock germplasm resources on the Tibetan Plateau, and have low oxygen, rough feeding, and disease resistance, In high-altitude environments, Tibetan sheep have developed particular biological habits and digestive, and metabolic mechanisms adapted to the alpine climate, demonstrating strong environmental adaptability. The rumen, a unique digestive organ in ruminants, is a complex ecosystem that plays a vital role in the metabolism, immune system, and health of the host (Lin et al., 2019). The rumen is inhabited by a wide variety of microorganisms, including bacteria, fungi, and protozoa These microorganisms make up a large microbiota that converts the ingested feed and other nutrients into microbial proteins and volatile fatty acids (VFAs) and synthesizes amino acids and proteins that cannot be synthesized by the animal itself, thus increasing the metabolic reservoirs of the host (Hess et al., 2011; Muegge et al., 2011). Therefore, rumen microorganisms play a vital role in the performance and health of ruminants. VFAs, as significant metabolites of rumen microorganisms, can provide approximately 75% of the energy requirements of ruminants (Bergman, 1990). They are also signaling molecules that regulate rumen epithelial cell proliferation and immune responses (Hosseinkhani et al., 2021; Li et al., 2021). Therefore, understanding the factors that affect the rumen microbiota and its metabolite composition will help reveal the biological mechanisms by which microorganisms and metabolites alter the body’s adaptability and elucidate the mechanisms of host–microbiota–metabolite interactions, thereby providing new coping strategies and methods.

Interactions between the gut microbiota and between the gut microbiota and the host are an essential part of the dynamic equilibrium of the gut microbiota, and changes in metabolites are closely related to changes in the gut microbiota. Factors that influence microbial composition and function include temperature, host species, diet, exercise, age, and geographic location (Wang et al., 2019). Still, these effects account for only 10 ~ 20% of microbial diversity, and the vast majority of microbial diversity between individuals remains unexplained. Therefore, the study of interactions between the integrated host genome, transcriptome, metabolome, and microbiome has received increasing attention. The gut microbiota and host coevolution coregulate host phenotypes. Under the same environmental feeding conditions, the composition of the gut microbiota of different species of Hydra varies widely. The composition of the gut microbiome of each of the two species of Hydra in the wild is very similar to that of the same species under laboratory conditions in that the host selectively determines the gut microbiome (Fraune and Bosch, 2007). Host genetic and environmental factors combine to influence the gut microbiota; In a study of two family lines of chickens with high and low body weights, the host genetic background has an essential influence on the composition and abundance of gut microorganisms. Some gut microbiota showed significant heritability and genetic correlation (Meng et al., 2014). In addition, there is evidence that host factors (such as breed, genetic variation, etc.) influence the composition and function of the bovine rumen microbiome (Zang et al., 2022). For example, the abundance of the rumen microbiota associated with CH₄ emissions in dairy cows is influenced by host genotype to some extent (Difford et al., 2018; Zhang et al., 2020). Host genetic factors also influence rumen microbiota formation in beef cattle (Abbas et al., 2020). In recent years, elucidating the mechanisms related to host biological indicators and phenotypic traits via multiomics has become a research hotspot (Wang et al., 2019). Zhang et al. reported that yaks and Tibetan sheep can adapt to extreme environments through coevolution of the rumen microbiota and host genes (Zhang et al., 2016). The regulatory mechanisms by which the rumen epithelial genome of Tibetan sheep interacts with the microbiota and its metabolites to adapt to the plateau cold season have been investigated previously via RNA sequencing (RNA-Seq) and microbiome analyses (Liu et al., 2022). All these studies indicate that microorganisms metabolically regulate the host phenotype. In addition, Tibetan sheep live at high altitudes throughout the year and have a strong adaptive ability to extreme plateau environments. In contrast, Hu sheep live at low altitudes throughout the year and have the advantages of high fertility, early sexual maturity, roughage tolerance, environmental adaptability (EEr et al., 2020; Li et al., 2022), and good adaptability to semiarid and semidesert areas, and they have been widely introduced in Northern China and high-altitude regions (Guo et al., 2023). However, only a few studies have been performed on rumen microbial–host interactions regulating plateau adaptation in Tibetan sheep and Hu sheep. Therefore, in this study, we compared the morphology of the rumen epithelium and fermentation function with microorganisms and their metabolites by using Tibetan sheep and Hu sheep reared in the same plateau environment to analyze and explore the phenotypic and genetic differences in plateau adaptation between Tibetan sheep and Hu sheep, to reveal the mechanisms by which the Tibetan sheep rumen microbiota and metabolites on rumen fermentation function and host interactions, and providing a basis for the study of plateau adaptation in Tibetan sheep.

2 Materials and methods

2.1 Experimental animals and sample collection

The experimental samples were selected from the same herder’s flock in Hezuo City, Gannan Tibetan Autonomous Prefecture, Gansu Province (103°E, 35°N), which has an altitude of 3,000 m, a humid and high-altitude climate, and a relatively fragile ecological environment. Six Tibetan sheep and 6 Hu sheep weighing 34 kg (± 0.5 kg) and aged 1 year (± 1 month) were selected for the experiment. The experimental sheep exhibited regular feeding and rumination, shiny fur, good physical condition, oval-shaped feces, and adherence to each other after landing, indicating good health. They were all in the local area’s traditional natural grazing management state. The types of forage mainly included Poa poophagorum Bor, Poaceae, Carex coninux, Argentina anserina, and Geranium platyanthum Duthie, without any supplementary feeding. The nutritional components of forage are shown in Supplementary Table S1. The test sheep were slaughtered at Xinling Livestock Product Development Co., Ltd. in Hezuo City, Gannan Tibetan Autonomous Prefecture, Gansu Province, and the rumen contents and tissues were collected. Before grazing in the morning, the rumen fluid was collected via a sheep gastric tube-type rumen sampler. Three tubes were collected from each sheep, and the collected rumen fluid was quickly put into a liquid nitrogen tank for freezing, then returned to the laboratory for preservation at −80°C and subsequent 16S rRNA analysis and metabolite identification. Immediately after the jugular vein was bled to death, the rumen was isolated, and a block of rumen vesicle tissue (1 cm²) was collected, lightly rinsed with saline to remove the contents, and fixed in 4% paraformaldehyde for histomorphological analysis. At the same time, a small piece of the rumen vesicle was clipped, and the contents were rinsed off rapidly with precooled saline at 4°C and PBS buffer, after which the epithelial tissues were separated with blunt scissors and rapidly placed in liquid nitrogen for subsequent total RNA extraction.

2.2 Morphological analysis

Rumen ventral sac tissues were fixed in 4% paraformaldehyde for 24 h. The fixed tissue samples were immersed in dehydrating agents in the following order: (1) 50% alcohol for 2 h; (2) 75% alcohol for 2 h, or could overnight; (3) 85% alcohol for 2 h; and (4) 95% alcohol for 1 h. The dehydrated tissue samples were washed with alcohol and then quickly transferred to xylene for clearing, and the tissue was transparent after approximately 15 min. After 15 min, the tissue was transparent. The transparent samples were transferred to melted paraffin, allowing the paraffin to contact and immerse the tissue fully, and then embedded in a wax block after 2 h. The wax block was fixed in a sectioning folder. After the wax block in the sectioning folder was fixed, the tissues were sectioned at a thickness of 3 μm, and the sections were spread and mounted in distilled water at 45°C and dried by baking at 60°C. The paraffin sections were stained in the following order: (1) Xylene I dewaxing, 10 ~ 15 min; (2) xylene II dewaxing, 1 ~ 2 min; (3) xylene: Anhydrous ethanol = 1:1 mixture, 1 ~ 2 min; (4) anhydrous ethanol, 1 ~ 2 min; (5) downward-gradient alcohol hydration: 95, 85, 75, and 50% alcohol, 1 ~ 2 min each; (6) gentle washing with distilled water, 1 ~ 2 min; (7) hematoxylin staining, 10 ~ 15 min; (8) gentle washing with distilled water to remove excess dye; (9) alcohol separation in hydrochloric acid, 30 ~ 60 s; (10) rinsing with tap water for bluing, 15 ~ 20 min; (11) gentle washing with distilled water, 1 ~ 2 min; (12) dehydration in upward-graded alcohols: 1 ~ 2 min each in 50%, 75, and 85% alcohols; (13) eosin staining, 2 ~ 3 min; (14) upward gradient alcohol dehydration: 75, 85, 95% alcohol, and anhydrous ethanol dehydration, 1 ~ 2 min each; and (15) xylene I, xylene II transparent, 15 min. Finally, neutral gum sealing, drying, observation under the microscope, and imaging were performed. Rumen muscle layer thickness, nipple height, nipple width, cuticle thickness, granular layer thickness, spinous layer thickness, and basal layer thickness were determined via Slide Viewer.

2.3 16S ribosomal RNA (rRNA) sequencing analysis of rumen microorganisms

Microbiome DNA was extracted from the rumen fluid of Tibetan sheep and Hu sheep using the MN NucleoSpin 96 Soi bacterial DNA extraction kit (Omega, Shanghai, China). Polymerase chain reaction (PCR) amplification of the V3–V4 region of the highly variable region of the 16S rRNA gene was performed using universal primers (forward primer 338F: 5’-ACTCCTACGGGGAGGCAGCA-3’ and reverse primer 806R: 5’-GGACTACHVGGGTWTCTAAT-3’). For library sequencing, small fragment libraries were constructed using a two-step library construction method with double-end sequencing, the amplified products were analyzed by library sequencing on the Illumina MiSeq 2,500 platform (Illumina, San Diego, CA, USA), and bioinformatics analysis was performed using BMKCloud¹ (accessed on 01 June 2024). To evaluate the quality of the raw data obtained from sequencing, paired-end splicing (FLASH version 1.2.7), filtering (Trimmomatic version 0.33), and chimeras (UCHIME version 4.2) were removed to obtain optimized sequences (Tags). The operational taxonomic units (OTUs) were obtained using USEARCH software. OTU classification and annotation were analyzed using the Silva (Bacteria) taxonomy database for further taxonomic analysis, and the community structure at different taxonomic levels (phylum, class, order, family, genus, and species) was obtained. The species diversity was analyzed via α-diversity, the α-diversity indices Ace, Chao1, Shannon, and Simpson were obtained, and a sample rarefaction curve was drawn. Furthermore, β-diversity analysis was used to obtain principal coordinates components (principal coordinate analysis [PCoA]) and boxplots based on multiple distances according to the distance matrix (Anosim). Line discriminant analysis effect size (LEfSe) analysis was used to find biomarkers with significant differences between groups. Metastats was used to conduct t-tests on species abundance data between groups. The species that caused the difference in the composition of the two groups of samples according to the q value were screened out, and 16S gene function analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups of proteins (COG) (Sha et al., 2022).

2.4 Metabolomic analysis of rumen microorganisms

Microbial metabolic profiling of rumen fluid from Tibetan sheep and Hu sheep (n = 12) was performed using a liquid chromatography–mass spectrometry (LC–MS) platform. After the samples were thawed at room temperature, 100 μL of each sample was weighed, and 500 μL of the extraction solution (methanol/acetonitrile volume ratio of 1:1, internal standard concentration of 2 mg/L) containing the internal standard (1,000:2) was added and vortexed for 30 s. The samples were then sonicated in an ice–water bath for 10 min, stood at −20°C for 1 h, and then centrifuged at 4°C at 12,000 rpm for 15 min. Then, 500 μL of the supernatant was placed in an EP tube, and the extract was dried in a vacuum concentrator, 150 μL of extract (acetonitrile−water volume ratio of 1:1) was added to the dried metabolite for redissolution, vortexing was continued for 30 s, sonication was performed for 10 min in an ice–water bath, and centrifugation was performed for 15 min at 12,000 rpm at 4°C. Finally, 120 μL of the supernatant was removed from a 2-ml injection bottle, and 10 μL of each sample was mixed with a QC sample for testing. The LC–MS system for metabolomics analysis was composed of a Waters Acquity I-Class PLUS Ultra High-Performance Liquid Tandem Waters Xevo G2-XS QToF High-Resolution Mass Spectrometer, and the column used here was purchased from a Waters Acquity UPLC HSS T3 column (1.8 μm 2.1 × 100 mm). The samples were eluted using positive (ESI+) and negative (ESI−) mobile phases consisting of water and 5% acetonitrile, 0.1% formic acid as solvent A and acetonitrile and 0.1% formic acid as solvent B at flow rates of 0.35 mL/min and 400 μL/min, respectively. The subsequent mobile phase (A:B) elution gradient was 0–0.25 min 98–2%, 10.0–13.0 min 2–98%, and 13.1–15.0 min 98–2%, followed by an ion source temperature of 150°C and a desolvation temperature of 500°C. The flow rates of the reverse gases and desolvent gas were 50 and 800 L/h, respectively. The raw data collected via MassLynx (version 4.2) were processed via Progenesis QI software. The Progenesis QI software, online databases such as METLIN and the self-built databases of BMG were used for metabolite identification according to sample type. BMKCloud was used to conduct a subsequent bioinformatics analysis of the identified metabolites. The screening criteria for differentially abundant metabolites were FC > 1, p < 0.05, and VIP > 1, and the differentially abundant metabolites were analyzed by KEGG functional annotation and enrichment.

2.5 Measurement of rumen VFAs

Thaw rumen fluid samples at room temperature, centrifuge at 15,000 rpm for 15 min at 4°C, aspirate 1 mL supernatant into a 1.5 mL centrifuge tube, add 0.2 mL of 25% metaphosphoric acid solution containing internal standard 2 EB, mix well and leave on ice for more than 30 min. Centrifuge at 15,000 rpm for 15 min, aspirate the supernatant with a 2 mL syringe and filter the supernatant with a 0.22-μm filter tip. The supernatant was aspirated with a 2-ml syringe and filtered with a 0.22-μm filter tip, and the resulting filtrate was transferred to a new collection tube for testing.

An Agilent 7890B gas chromatograph was used for the determination. The chromatographic column was an AT-FFAP capillary column (30 m × 0.32 mm × 0.50 μm), the temperature of the inlet (SSL): 250°C; the temperature of the detector (flame ionization detector [FID]): 250°C; the carrier gas was high-purity nitrogen (99.999%), the total pressure was 100 kPa, and the splitting ratio was 5:1; the gas flow rate was as follows: air: 400 mL·min⁻¹, H₂: 35 mL·min⁻¹, and N₂: 40 mL·min⁻¹; injection volume was 1 μL; and the heating procedure was as follows: 120°C held for 3 min, 10°C·min⁻¹ heating to 180°C, held for 1 min.

2.6 Measurement of rumen VFA transporter-related gene expression

Total RNA was extracted from the rumen tissues of Tibetan sheep and Hu sheep using the TRIzol reagent method (DP762-T1C; Life Technologies, CA, USA). The RNA concentration and purity were determined using an ultramicro spectrophotometer (Therm Nano Drop-2000; Thermo Scientific, MA, USA), and the OD260:OD280 ranged 1.8–2.1, indicating that the purity of the extracted RNA was good. The integrity of the extracted RNA was examined using an agarose gel electrophoresis instrument (Agient2100, LabChip GX), and three single and well-defined bands were visualized using an agarose gel imaging system (Ready Agarose, Bio-Rad, Hercules, CA, USA), indicating that the integrity of the extracted RNA was good. cDNA synthesis was performed using a reverse transcription kit (HiScript® II Q RT SuperMix for qPCR; Nanjing, China). The primers for anion exchanger 2 (AE2), downregulated in adenoma (DRA), monocarboxylic acid transporter protein1 (MCT1), monocarboxylic acid transporter protein (MCT4), Na⁺/H⁺ exchanger1 (NHE1) and Na⁺/H⁺ exchanger2 (NHE2) genes were designed using Primer5.0 software (β-actin was used as an internal reference gene; see Supplementary Table S2 for detailed primer information). The relative expression of genes involved in the transport of VFAs in the rumen tissues of Tibetan sheep and Hu sheep was investigated using an Applied Biosystems Q6 real-time fluorescence quantitative PCR instrument. Real-time quantitative PCR (RT-qPCR) was performed using a 20-μl reaction system with the following reaction conditions: predenaturation at 95°C for 30 s; cycling at 95°C for 10 s and 60°C for 30 s for 40 cycles; and lysis curve (95°C for 15 s, 60°C for 60 s, and 95°C for 15 s). The relative gene expression was calculated by correcting for β-actin as an internal reference gene, and the data were analyzed by the 2^−∆∆CT method.

2.7 Statistical analysis of the data

The experimental data were initially organized using Excel 2016, and the analyzed data were expressed as “mean ± standard error” with a statistical significance level of p < 0.05. Then the data were analyzed using the Statistical Package for the Social Sciences (SPSS) version 24.0 software. The content of VFAs was analyzed using an independent samples t-test, and the relative expression of transporter genes was analyzed using the 2^−∆∆CT method. Correlation analyses were performed using Spearman’s correlation test with a screening criterion of p < 0.05.

3 Results

3.1 Morphological structure of the rumen epithelium in Tibetan sheep and Hu sheep

As shown in Table 1, the nipple height of the rumen epithelium of Tibetan sheep is significantly greater than that of Hu sheep (p < 0.01), and the nipple width of Hu sheep is significantly greater than that of Tibetan sheep (p < 0.05). The sheep rumen epithelial layer is divided from outside to inside into the cuticle (a), granular layer (b), sphenoid layer (c), and basal layer (d) (Figures 1B,E). As shown in Figure 1, the differences in the rumen epithelial layer between Tibetan sheep and Hu sheep included a significantly greater thickness of the granular layer and sphenoid layer in Hu sheep than in Tibetan sheep (p < 0.01), a significantly greater thickness of the cuticle layer in Tibetan sheep than in Hu sheep (p < 0.01), and a non-significant difference in the thickness of the basal layer (p > 0.05). In addition, the differences in the morphology of the muscular layer of the rumen wall between Tibetan sheep and Hu sheep, in which the thickness of the muscular layer is significantly greater than that of Tibetan sheep (p < 0.01), and the number of connective tissues in the muscular layer of the Hu sheep is greater than that of the Tibetan sheep (Figures 1A,C,D,F).

Table 1

Table 1. Morphological analysis of rumen epithelium of Tibetan sheep and Hu sheep.

Figure 1

Figure 1. Morphology of rumen epithelial tissue of Tibetan sheep and Hu sheep. (A–C) The morphology of rumen epithelial tissue of Tibetan sheep at 1×, 10×, and 2× magnifications, respectively; (D–F) The morphology of rumen epithelial tissue of Hu sheep at 1×, 10×, and 2× magnifications, respectively. a: Cuticle; b: Granular layer; c: Spinous layer; d: Basal layer.

3.2 Sheep rumen 16S measurement

16S rRNA sequencing of the rumen fluid from Tibetan sheep and Hu sheep (n = 12) and PCoA reveal that there are significant differences in the rumen microbes between Tibetan sheep and Hu sheep, with good within-group reproducibility (Figure 2A). A total of 16,334 OTUs are obtained, including 7,246 OTUs in Tibetan sheep and 8,386 OTUs in Hu sheep (Figure 2B). The number of OTUs specific to Hu sheep is significantly greater than that specific to Tibetan sheep (p < 0.05). Dilution curves are used to verify that the amount of sequencing data is sufficient to reflect the diversity of species in the samples and the abundance of species in the samples. As shown in Figure 2C, the curve flattens at 40,000 Reads, indicating that the amount of sequencing data is sufficient. To compare their differences more systematically, the rumen microbial compositions are analyzed. α-Diversity analysis reveals that the ACE, Chao1, and Shannon indices in the rumen microorganisms of Hu sheep are significantly greater than those in the rumen microorganisms of Tibetan sheep (p < 0.05), and the difference in the Simpson index is not significant (p > 0.05) (Figures 2D–G). The microbial diversity index of Hu sheep is significantly higher than that of Tibetan sheep. The higher the α-diversity is, the more complex and stable the composition of the gut microbiota, indicating that Hu sheep are more resistant and adaptable to external interference, which is beneficial for host health.

Figure 2

Figure 2. Analysis of microbial diversity in Tibetan sheep and Hu sheep. (A) Scatter plot of principal coordinates components (PCoA)-score showing similarity of the microbiological composition of Tibetan sheep and Hu sheep based on UniFrac distance, principal coordinates (PCs) 1 and 2 explained 14.51 and 12.37% of the variance, respectively. (B) OTU Venn diagram. (C) Sample dilution curve. (D–G) Microbial diversity indicators—(D): ACE index; (E): Chao1 index; (F): Shannon index; and (G): Simpson index.

At the phylum level, Firmicutes, Bacteroidetes, Patescibacteria, Spirochaetota, and Verrucomicrobiota are the dominant phyla (Figure 3A). The abundances of Firmicutes and Patescibacteria are significantly greater (p < 0.05) in the Hu sheep than in the Tibetan sheep. In contrast, the abundance of Bacteroidetes is considerably greater (p < 0.05) in the Tibetan sheep than in the Hu sheep (Figure 3B). At the genus level, uncultured_rumen_bacterium, Prevotella, Rikenellaceae_RC9_gut_group, unclassified_F082, and Succiniclasticum are the dominant genera (Figure 3C), and in Tibetan sheep Prevotella and Rikenellaceae_RC9_gut_group abundance is significantly greater than that of Hu sheep (p < 0.05) (Figure 3D).

Figure 3

Figure 3. Analysis of the microbiological composition of Tibetan sheep and Hu sheep. (A,C) Microbial composition at the level of phylum and genus. (B,D) Analysis of species differences at phylum and genus level.

Ruminal differential microbial KEGG functional enrichment reveals that Membrane transport, Translation, Replication, and repair are significantly greater in Hu sheep compared to Tibetan sheep, while carbohydrate metabolism, glycan biosynthesis and metabolism, biosynthesis of other secondary metabolites, lipid metabolism, and signal transduction are significantly higher in Tibetan sheep than in Hu sheep (Figure 4A). COG functional analysis reveals that posttranslational modification, protein turnover, chaperones, and cell wall/membrane/envelope biogenesis functions were all significantly higher in Tibetan sheep than in Hu sheep, whereas in the rumen of Hu sheep, the ribosomal structure and biogenesis and amino acid transport and metabolism functions are greater in Hu sheep than in Tibetan sheep (Figure 4B).

Figure 4

Figure 4. Intestinal microbial function analysis of Tibetan sheep and Hu sheep. (A) Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analysis. (B) Clusters of Orthologous Groups of proteins (COG) functional analysis.

3.3 Rumen microbial metabolomics assay

Analysis of rumen microbiota metabolites in Tibetan sheep and Hu sheep yielded a total of 959,773 pairs of reads. After double-ended reads quality control and splicing, a total of 877,838 clean reads are generated, with at least 72,663 clean reads produced per sample and an average of 73,153 clean reads produced. A total of 3,982 metabolites are identified. Principal component analysis (PCA) can reduce high-dimensional data to two-dimensional or three-dimensional space, making it easier to visualize and analyze the data. It can also identify the most representative information in the data, providing a visual understanding of its distribution characteristics and differences. The PCA results of this study reveal significant differences in the rumen microbial species composition between the Tibetan sheep and Hu sheep groups, with good intragroup reproducibility (Figure 5A). Orthogonal projections to late structures discriminant analysis (OPLS-DA) can be used to screen differentially abundant metabolites effectively in metabolomics analysis. OPLS-DA can be used to obtain intergroup differential information better by establishing a model between metabolite expression levels and grouping relationships and screening differentially abundant metabolites that contribute significantly to distinguishing different groups. OPLS-DA reveals that the R2X, R2Y, and Q2 values are all close to 1, with Q2 > 0.9, further validating the reliability of the OPLS-DA model (Figure 5B). Further cluster analysis is conducted to obtain the number of OTUs in each sample. At the taxonomic level, a total of 26 phyla, 58 classes, 124 orders, 210 families, 368 genera, and 450 species are detected. Differentially abundant metabolites are analyzed using FC > 1, p < 0.05, and VIP > 1 as the screening criteria, and a total of 1,802 differentially abundant metabolites are identified, of which 948 are upregulated and 854 are downregulated (Figure 5C). Further screening of the top 10 up- and downregulated metabolites with multiplicative differences reveal that phalloidin, hexahydro-4-methylphthalic anhydride, 12,15-epoxy-13,14-dimethyleicosa-12,14,16-trienoic acid, 2-methoxynaphthalene, 18-fluoro-octadecanoic acid, 6”-O-carbamoylkanamycin A, and Temurin are significantly upregulated in the rumen of Tibetan sheep and very low in Hu sheep, these metabolites play an important role in lipid metabolism, and can also promote bile secretion, inhibit cancer cell proliferation, and participate in immune regulation. Tabersonine, N-lauroyl arginine, epidermin, glucoconvallasaponin B, miltefosine, N-oleoyl phenylalanine, and ajmalicine are significantly upregulated in the rumen of Hu sheep, which are involved in the regulation of inflammation to ensure normal growth and development of Hu sheep (Figure 5D).

Figure 5

Figure 5. Quality control diagram of rumen microbial metabolome data of Tibetan sheep and Hu sheep. (A) Scatter plot of principal component analysis (PCA)-score showing the similarity of the rumen metabolites of Tibetan sheep and Hu sheep based on UniFrac distance, principal components (PCs) 1 and 2 explained 43.59 and 24.61% of the variance, respectively. (B) Orthogonal projections to late structures discriminant analysis (OPLS-DA) based on operational taxonomic units (OTUs) of rumen metabolites in Hu sheep and Tibetan sheep. (C) Volcanic maps of differential metabolites in Tibetan sheep and Hu sheep. (D) Column chart of positive ion difference multiples.

The results of KEGG functional annotation and enrichment of different metabolites in the rumens of Tibetan sheep and Hu sheep are shown in Figure 6, and the different metabolites are enriched mainly in the lipid metabolism, digestive system, nucleotide metabolism, and amino acid metabolism pathways. Among them, the upregulated differentially abundant metabolites are enriched mainly in the caffeine metabolism and lysine degradation pathways, and the downregulated differentially abundant metabolites are enriched in a total of 18 metabolic pathways, with significant enrichment in the purine metabolism, amino sugar, and nucleotide sugar metabolism, folate biosynthesis, fatty acid biosynthesis, and mineral absorption pathways.

Figure 6

Figure 6. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analysis diagram of differential metabolites in rumen microbiomes of Tibetan sheep and Hu sheep. (A) Differential metabolite KEGG pathway classification annotation. (B) Differential metabolite differential abundance score diagram.

3.4 Rumen VFAs concentration and transporter gene expression

As shown in Table 2, there is a significant difference in the concentration of rumen VFAs between Tibetan sheep and Hu sheep, in which the total rumen VFAs concentration of Hu sheep is significantly greater than that of Tibetan sheep (p < 0.05). The concentrations of acetate, propionate, butyrate, and valerate in the Hu sheep rumen are greater, and the differences are significant (p < 0.05) for acetate, propionate, and valerate. The concentrations of isobutyrate and isovalerate acids are significantly greater (p < 0.05) in Tibetan sheep than in Hu sheep, and the Acetate/Propionate (A/P) values are significantly greater (p < 0.05) in Tibetan sheep than in Hu sheep. The results of the relative expression measurements of rumen VFAs transporter genes (AE2, DRA, MCT1, MCT4, NHE1, and NHE2) are shown in Figure 7. The expression levels of AE2, DRA, NHE1, and NHE2 are significantly greater in Tibetan sheep than in Hu sheep (p < 0.05), whereas the expression levels of MCT1 and MCT4 are significantly greater in Hu sheep than in Tibetan sheep (p < 0.05).

Table 2

Table 2. Contents of volatile fatty acids (VFAs) in the rumen of Tibetan sheep and Hu sheep.

Figure 7

Figure 7. Rumen volatile fatty acids (VFAs) transport gene expression results of Tibetan sheep and Hu sheep. **Indicates highly significant difference (p < 0.01).

3.5 Analysis of rumen microbe–metabolite interactions

Analysis of rumen microbiota and its metabolites reveals a certain correlation between the two (Supplementary Figure S3). Further correlation analysis was conducted, and the heatmap of microbiota and metabolites at the phylum level (Figure 8A) revealed that microorganisms are correlated with 12 metabolite modules, of which Bacteroidota and Synergistota are highly significantly positively correlated with three metabolite modules, MEturquoise, MEmagenta, and MEpurple (p < 0.01); Firmicutes is highly significantly positively correlated with MEbrown and MEyellow; two modules are highly significantly positively correlated (p < 0.01); and Patescibacteria is highly significantly positively correlated (p < 0.01) with the MEbrown and MEgreen modules, whereas Bacteroidetes and Synergistota are highly significantly negatively correlated (p < 0.01) with the MEbrown, MEyellow and MEgreen modules. Further correlation and chord plots of microorganisms and metabolites at the genus level revealed that rumen microorganisms and metabolites are also inextricably linked. Fretibacterium is significantly negatively correlated (p < 0.05) with Glutamine glutamate aspartate, Aldoxorubicin, Gyrocyanin, and α-Fluoromethylhistamine, whereas Pediococcus is significantly positively correlated (p < 0.01) with glutamine, glutamate, and aspartate (p < 0.01). Acetitomaculum and Pediococcus are strongly positively correlated with Aldoxorubicin, Gyrocyanin, and α-Fluoromethylhistamine (p < 0.01) (Figure 8B).

Figure 8

Figure 8. Correlation analysis of rumen microbiomes and metabolites in Tibetan sheep and Hu sheep. (A) Microbiomes and metabolites correlation heat map. (B) Microbiomes and metabolites correlation chord diagram. *p < 0.05, **p < 0.01, ***p < 0.001.

3.6 Analysis of rumen microbe-VFAs and transporter gene–epithelial morphology interactions

Interaction analysis of the top 10 rumen microorganisms at the genus level with VFAs and their transporter genes and rumen epithelial morphology reveals that Rikenellaceae_RC9_gut_group is significantly positively correlated (p < 0.05) with nipple height and the AE2 gene and negatively correlated (p < 0.05) with the MCT1 and MCT4 genes and nipple width, Succiniclasticum and Saccharofermentans are strongly positively correlated with propionate, butyrate, valerate, and MCT4 genes, as well as with nipple width, granular layer thickness, and sphenoid layer thickness (p < 0.01) and Succiniclasticum is also highly significantly positively correlated with acetate and MCT1 genes (p < 0.01), and highly significantly negatively correlated with the AE2, DRA, NHE1, and NHE2 genes, as well as with nipple height and cuticle thickness, which are highly significantly negatively correlated (p < 0.01), (Figures 9A,B). Spearman correlation analysis of rumen VFAs, transporter genes, and their rumen epithelial micromorphology in Tibetan sheep and Hu sheep. As shown in Figure 9C, the AE2, DRA, and NHE2 genes are significantly and positively correlated (p < 0.05) with the contents of acetate and total acids, and the DRA and NHE1 genes and their rumen epithelial cuticle thickness are significantly and positively correlated (p < 0.05) with the contents of propionic acid and valerates. The AE2 and NHE1 genes and the nipple height are significantly positively correlated with the butyrate content (p < 0.05), whereas the MCT1 and MCT4 genes and the nipple width are significantly negatively correlated with the butyrate content (p < 0.05). These findings suggest that rumen transporter genes and the rumen epithelial structure are inextricably linked to the production and transporter uptake of VFAs.

Figure 9

Figure 9. Correlation heat map between Tibetan sheep and Hu sheep. (A) Heat map of correlation between horizontal rumen microbiomes and VFAs and transport genes. (B) Heat map of morphological correlation between rumen microbiomes and rumen epithelium. (C) Heat map of correlation between rumen VFAs and transport genes.

3.7 Weighted correlation network analysis (WGCNA) of rumen microbial metabolite-VFAs transporter genes

Based on the results of the rumen microbial metabolomics analysis, further WGCNA of 1,802 differentially abundant metabolites reveal that 10 modules are differentially correlated with rumen epithelial morphology and transporter genes related to VFAs in Tibetan sheep and Hu sheep (Figure 10). There is a strong correlation between the five modules, MEgreen, MEbrown, MEyellow, MEblack, and MEturquoise, and rumen epithelial morphology and VFAs-related transporter genes. The MEgreen, MEbrown, and MEyellow modules are significantly and positively correlated with the MCT1 and MCT4 genes and the thicknesses of the granular and sphenoid layers. The MEblack and MEturquoise modules are positively associated with the AE2, DRA, NHE1, and NHE2 genes, nipple height, and cuticle thickness, whereas they are negatively correlated with the MCT1 and MCT4 genes. Further analysis of the metabolites with the top 10 kME values in the five modules revealed that the metabolites in the MEgreen, MEbrown, and MEyellow modules are enriched in the biosynthesis of unsaturated fatty acids and Fructose and mannose metabolisms such as downregulated metabolites such as alginic acid and linoleic acid; In contrast, in the MEblack and MEturquoise modules, some upregulated metabolites, such as 3α7α,12α-trihydroxy-5β-cholestanoate, 5-hydroxyectoine, and dihydrotestosterone, which are found in Tibetan sheep, were significantly more abundant than those in Hu sheep and are enriched mainly in the pathways of primary bile acid biosynthesis, glycine, serine and threonine metabolism, and steroid hormone biosynthesis.

Figure 10

Figure 10. The results of WGCNA. (A) Cluster dendrogram of rumen metabolites in Tibetan sheep and Hu sheep. Each branch in the figure represents one metabolite, and each color represents one coexpression module. (B) Correlation between the rumen metabolite modules, transport genes and rumen epithelial morphology in Tibetan sheep and Hu sheep.

4 Discussion

This study aimed to analyze the breed differences in the rumen epithelial micromorphology, microbiome, metabolome, and expression levels of VFAs and their transporter genes between Tibetan sheep and Hu sheep under the same feeding conditions and to investigate the mechanism of plateau adaptation in Tibetan sheep on the basis of the rumen microbiome and its metabolome. The rumen is one of the most important sites for digestion, metabolism, and nutrient absorption in ruminants, and previous studies have shown that nutrient absorption across the rumen epithelium is mainly dependent on the degree of keratinization of the cuticle cells (Baldwin and Jesse, 1992). In this study, the thickness of the rumen nipple and cuticle is greater in Tibetan sheep than in Hu sheep because of their greater digestibility of alpine forage of coarse quality and high fiber content, which require strong rumen contraction and diastole to achieve mechanical digestion. The sphenoid layer is the leading site of VFA metabolism in rumen epithelial tissue. Compared with Tibetan sheep, Hu sheep need to metabolize more VFAs for energy to adapt to the plateau environment. The lower thickness of the rumen sphenoid layer in Tibetan sheep indicates that Tibetan sheep inhibit further metabolism of VFAs, and store them for their energy supply under conditions of hypoxia and severe nutritional deficiency. The height of the rumen nipple in Tibetan sheep is greater than that in Hu sheep when stimulated by plateau-wilting grasses. Therefore, a higher rumen nipple is a trait and strategy for the digestion and utilization of alpine pastures in Tibetan sheep and their adaptation to the plateau environment.

Rumen microorganisms play an important role in the digestive metabolism of proteins, sugars, and fats in ruminants. In this study, Firmicutes and Patescibacteria are more abundant in the rumen of Hu sheep at the phylum level, and Firmicutes promote fiber degradation, growth, and development (Lin et al., 2018). In contrast, the abundance of Bacteroidetes is greater in Tibetan sheep and has been observed to be much lower than Firmicutes in the rumen microbiota of low-altitude ruminants (Kim et al., 2011) and goats from semiarid regions (Cunha et al., 2011). Bacteroidetes can efficiently catabolize dietary proteins and carbohydrates into VFAs (Huo et al., 2014) and promote increased rumen growth and volume. The large number of bacteria in Bacteroidetes in plateau ruminants highlights the vital role of these bacteria, which may be beneficial in providing nutrients to the host by degrading the limited resources in the Qinghai-Tibet Plateau (QTP) (Liu et al., 2020). Furthermore, Bacteroidetes play an important role in nutrient uptake and utilization and host disease immunity (Zhang et al., 2018), thus playing a more important role than Firmicutes in the adaptation of animals to high-altitude environments. These findings suggest that Tibetan sheep have a stable energy utilization mechanism when the pasture supply is insufficient at high altitudes, which is important for maintaining homeostasis in the rumen and adapting to harsh environments (Jami et al., 2014). At the genus level, Rikenellaceae_RC9_gut_group, which degrades polysaccharides of plant origin (Peng et al., 2015), is significantly more abundant (p < 0.05) in the rumen of Tibetan sheep than in that of Hu sheep. In addition, Prevotella is also a dominant genus in Tibetan sheep, and Prevotella plays a role in the degradation and utilization of plant non-cellulosic polysaccharides, proteins, starches, and xylans (Liu et al., 2019). A high abundance of Prevotella promotes forage fermentation, resulting in high concentrations of VFAs (Shi et al., 2014), which facilitate the ability of rumen-fermented forage to provide more energy to the organism. Furthermore, KEGG functional enrichment analysis reveals that Carbohydrate metabolism, Glycan biosynthesis and metabolism, and Lipid metabolism are significantly enhanced in Tibetan sheep, which are able to satisfy their energy needs in the alpine hypoxic environment (Fan et al., 2021), whereas functions such as Membrane transport, which can mediate interactions between the intestinal microbiota and host cells (Konishi et al., 2015), are significantly elevated in the rumen of Hu sheep. The COG results also reveal that posttranslational modification, protein turnover, chaperones, and cell wall/membrane/envelope biogenesis are increased in Tibetan sheep. Protein turnover breaks down proteins into amino acids during food digestion for growth and development and is also involved in cellular signaling to adapt to environmental changes (Gonzalez and Cullen, 2022). The above results suggest that, compared with Hu sheep under the same feeding conditions, Tibetan sheep may have more powerful energy metabolism mechanisms to adapt to the Tibetan Plateau environment than do Hu sheep because their genetic factors lead to increased abundances of Bacteroidetes and Prevotella, effective fermentation and utilization of forage, and further enhancement of carbohydrate and lipid metabolism.

Metabolomics plays an important role in understanding the physiological and biochemical state of animals (Xue et al., 2020).In the present study, some differences in rumen are detected metabolites between Tibetan sheep and Hu sheep. Tabersonine, Epidermin, and Glucoconvallasaponin B levels are increased in Hu sheep, whereas metabolites such as phalloidin and 12,15-epoxy-13,14-dimethyleicosa-12,14,16-trienoic acid are increased in the rumens of Tibetan sheep. Phalloidin promotes bile secretion(Rahman and Coleman, 1986), facilitates fat digestion, and inhibits bacterial proliferation, and 12,15-epoxy-13,14-dimethyleicosa-12,14,16-trienoic acid is a furan fatty acid derivative that has been shown to inhibit cancer cell proliferation and block cell recovery. These results suggest that differentially abundant metabolites play important roles in fatty acid catabolism and the intestinal immune barrier in Tibetan sheep, but the specific regulatory mechanisms need to be further verified. Further analysis of the KEGG functional enrichment of the differentially abundant metabolites reveals that the differentially abundant metabolites of Tibetan sheep and Hu sheep are enriched mainly in the Lipid metabolism, Digestive system, nucleotide metabolism, and amino acid metabolism pathways. Among them, lipid metabolism plays an important role in energy metabolism as well as in various aspects of biofilm structure, signaling and other functions, and disorders of lipid metabolism can lead to a variety of diseases (Petrenko et al., 2023), and its enrich metabolites, such as cholesterol, the expression of which is upregulated in the rumen of Tibetan sheep, are precursors for the synthesis of a variety of hormones, such as steroid hormones and glucocorticosteroids, which are related to immunity and energy metabolism (Cortes et al., 2014). In addition, caffeine metabolism and Lysine degradation are upregulated in Tibetan sheep, and its upregulated metabolites, Theophylline and Xanthine, are associated with the regulation of inflammation (Rabe and Dent, 1998; Yagi et al., 2022). Pyruvate and alanine produced by Lysine degradation can further participate in energy metabolism and amino acid metabolism processes. Purine metabolism, amino sugar and nucleotide sugar metabolism, and fatty acid biosynthesis are significantly enriched in Hu sheep and participate in various biological processes to ensure normal growth and development. The above analyses reveal significant differences in rumen microbial metabolites and their enrichment pathways between Tibetan sheep and Hu sheep under the same feeding environment, which is due mainly to breed differences, with Tibetan sheep producing more metabolites related to organismal energy provisioning and immunomodulation than Hu sheep do.

Under the regulation of rumen microorganisms and their metabolites, VFAs also change accordingly. Energy supply is the most important role of VFAs in ruminants, and the rumen acetate content accounts for 70 ~ 75% of the total VFA content. Increasing the level of acetate in the rumen not only alters the fermentation pattern in the rumen, improving the efficiency of energy metabolism and reducing energy losses in ruminants, but also improves the utilization of nitrogen. Propionate is largely converted to glucose, and in sheep, the glucose produced from propionate provides 30 ~ 50% of the requirement. The concentrations of both acetate and propionate in the rumen of Hu sheep are significantly higher than those in the rumen of Tibetan sheep, indicating that the rumen of Hu sheep has a greater fermentation capacity for forage than that of Tibetan sheep. Most of the butyrate in the rumen is converted by rumen epithelial cells into ketones or CO₂, which are byproducts of fatty acid oxidation (McGarry and Foster, 1980). The A/p value is related to the efficiency of forage energy utilization; the lower the ratio is, the greater the energy utilization efficiency of the organism (Liu et al., 2019). The A/p value of Hu sheep is significantly lower than that of Tibetan sheep, which indicates that Hu sheep have a higher utilization rate of forage and can produce more VFAs; however, in plateau environments, the metabolism of energy substances in Tibetan sheep is greater than that in Hu sheep. VFAs, as important products of rumen fermentation, are absorbed and transported in the rumen epithelium by passive diffusion and anion exchange (del Bianco Benedeti et al., 2018), and the substances that play these roles are mainly transporter genes, such as the AE2, DRA, NHE, and MCT genes. Both AE2 and DRA encode anion-exchange proteins that exchange VFA⁻ for HCO₃⁻ and play important roles in regulating the stability of the internal environment. Zhang et al. reported that AE2 may be involved in maintaining the intracellular pH balance during the transport of VFAs (Zhang et al., 2023), which indirectly affects the transport and metabolism of VFAs, and high expression of AE2 may indicate that the cellular capacity for bicarbonate transport is increased, which contributes to maintaining the stability of the intracellular environment and supports the further metabolism of VFAs. Bilk et al. (2005) reported that the rumen epithelium exchanges HCO₃⁻ and Cl⁻ via AE2, DRA, and PAT1 and is associated with the uptake and transport of VFAs in this process. NHE is a membrane transporter of Na⁺ and H⁺ that transports VFAs mainly by passive diffusion (Müller et al., 2002), and it has been demonstrated that NHE1 and NHE2, among others, are expressed in the rumen epithelium and that their activities can influence the uptake of VFAs by the rumen epithelium. Yang et al. (2012) reported that the expression and activity of rumen epithelial NHEs were increased when goats were fed highly fermented diets, and in the present study, the expression of the above four genes in the rumen of Tibetan sheep was significantly greater than that in the rumen of Hu sheep, indicating that the absorption and metabolism of VFAs by Tibetan sheep in the plateau environment is stronger than that of Hu sheep and that the content of VFAs is lower than that of Hu sheep. Monocarboxylic acid transporter (MCT) proteins, which are distributed at the plasma membrane surface of cells in the basal layer of the rumen epithelium, are involved in regulating the efflux of H⁺ lactate and ketone bodies and maintaining intercellular pH homeostasis (Kirat et al., 2006). The results of the present study revealed that the expression of MCT1 and MCT4 in the rumen epithelium of Hu sheep was significantly greater than that in the rumen epithelium of Tibetan sheep. Pérez de Heredia et al. (2010) reported that hypoxia can upregulate the expression of MCT1 and MCT4 in human adipocytes and that the impression of MCT4, similar to other enzymes involved in glucose catabolism, is increased upregulated under hypoxic stimulation through the regulatory mechanism of HIF-1α, and its upregulated expression enables the rapid efflux of lactic acid from the cell due to sugar catabolism during hypoxia, which is in line with the results of the present study. Therefore, efficient uptake and metabolism of VFAs constitute another characteristic and effective strategy for Tibetan sheep to cope with the harsh alpine environment.

On the basis of the above analyses, we identified a certain synergistic relationship between the rumen microbiome and metabolome and their fermentation functions. In this study, the modules that are significantly positively correlated with Bacteroidetes are mainly related to immune-related metabolites such as Nylidrin, Physalin E acetate, and syringic acid. Among them, syringic acid, which has a structural unit similar to that of lignin, acts as an antioxidant to reduce oxidative stress and exerts antioxidant and anti-inflammatory effects (Ferah Okkay et al., 2022). Bacteroidota, a dominant genus of bacteria in Tibetan sheep, has been reported to be associated with microbial dysbiosis and chronic intestinal inflammation (Parker et al., 2020), suggesting that Bacteroidota and these metabolites are collectively involved in the regulation of intestinal immunity in Tibetan sheep, thus ensuring their health. Pediococcus spp. strains have the potential to inhibit mycotoxin-producing molds as producers of potent antifungal metabolites (Fugaban et al., 2022), which are significantly and positively correlated with metabolites such as Glutamine glutamate aspartate, which is a major metabolic fuel for the small intestine and can improve hepatic energy metabolism (Qi et al., 2020). In addition, there was an association between the rumen microbiota flora, rumen epithelial micromorphology, and VFAs in Tibetan sheep. Rikenellaceae_RC9_gut_group is significantly and positively correlated with nipple height and the AE2 gene (p < 0.05), while Rikenellaceae_RC9_gut_group may promote the digestion and absorption of carbohydrates in the intestinal tract (Berry et al., 2015) and inhibit the production of proinflammatory cytokines and other damaging factors (Bian et al., 2019), and Rikenellaceae_RC9_gut_group and Succiniclasticum may be involved in methane formation and VFAs production. In this study, Rikenellaceae_RC9_gut_group and Succiniclasticum as the dominant genera, increased the content of VFAs, indicating that Tibetan sheep regulate the digestion and absorption of VFAs through the upregulation of genes, such as Rikenellaceae_RC9_gut_group, to promote the growth of the rumen nipple and AE2. This, in turn, enhances energy metabolism and immune regulation in Tibetan sheep to ensure their growth and development. Therefore, rumen microbes and their metabolites synergistically regulate VFAs and rumen epithelial morphology to improve energy metabolism and reduce oxidative stress in Tibetan sheep in the same feeding environment.

Furthermore, WGCNA of rumen metabolites reveals that metabolite modules associated with highly expressed VFAs transporter genes and rumen epithelial micromorphology act mainly through saturated and unsaturated fatty acid biosynthesis and Fructose and mannose metabolism pathways to maintain normal physiological activities in Hu sheep. Several metabolites related to energy metabolism and immune barriers have been identified in Tibetan sheep, among which 3α,7α,12α-trihydroxy-5α-cholestanoate is an important precursor for bile acid synthesis(Mihalik et al., 2002). Bile acid promotes both the digestion and absorption of dietary lipids and acts as a hormone that activates specific receptors. The activation of these receptors can alter gene expression in a wide range of tissues, leading to changes not only in bile acid metabolism but also in glucose homeostasis, lipid and lipoprotein metabolism, energy expenditure, intestinal motility, and bacterial growth, and inflammation (de Aguiar et al., 2013). 5-Hydroxyectoine can bind to some 5-hydroxyectoine receptors on the surface of immune cells to regulate the immune system and participate in immune metabolism. In addition, many metabolites related to amino acids, such as Gly, Ser, and Thr, are involved in immune function, anti-inflammatory processes, and antioxidant responses, in addition to meeting the normal developmental and amino acid requirements of Tibetan sheep (Imenshahidi and Hossenzadeh, 2022). The interactions between the above metabolites and VFAs and epithelial morphology further provide insight into plateau adaptation in Tibetan sheep (Figure 11).

Figure 11

Figure 11. Modeling of rumen microorganisms and their fermentative metabolism in Tibetan sheep and Hu sheep.

5 Conclusion

In the case of rough grass quality on the plateau, the high rumen nipple and thick cuticle of Tibetan sheep are the main differences in their adaptation to high-altitude environments compared to Hu sheep, Second, the rumen microbiota of Tibetan sheep has increased abundances of Bacteroidetes, Prevotella and Rikenellaceae_RC9_gut_group, which efficiently ferment and utilize forage and absorb the metabolize VFAs, promote the expression of VFAs transporter genes, and improve carbohydrate, amino acid and lipid metabolism, thus giving Tibetan sheep a more powerful energy metabolism mechanism to adapt to the environment. In addition, Tibetan sheep can produce more metabolites related to the organismal energy supply and immune regulation, and these differential metabolites are associated with purine metabolism, amino sugar, and nucleotide sugar metabolism, Fatty acid biosynthesis, and mineral absorption pathways, which are involved in the regulation of rumen growth and development and gastrointestinal homeostasis. Correlation analysis revealed that Tibetan sheep rumen microbes and their metabolites synergistically regulated VFAs and epithelial morphology, enhanced their energy metabolism, and reduced oxidative stress.

Data availability statement

The data sets presented in this study can be found in the NCBI Sequence Read Archive (SRA) under accession number PRJNA1135557.

Ethics statement

The animal study was approved by the Livestock Care Committee of Gansu Agricultural University (Approval No. GAU-LC-2020-27). The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

QC: Writing – original draft. YS: Writing – review & editing, Conceptualization. XL: Writing – review & editing. YH: Writing – review & editing, Conceptualization. XC: Writing – review & editing, Data curation. WY: Writing – review & editing, Data curation. MG: Writing – original draft, Data curation. WH: Writing – review & editing, Data curation. JW: Writing – review & editing, Conceptualization. JH: Writing – review & editing, Data curation. LW: Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was funded by the National Natural Science Foundation of China (32260820), the Discipline Team Project of Gansu Agricultural University (GAU-XKTD-2022-21), the Gansu Agricultural University Youth Mentor Support Fund project (GAU-QDFC-2022-06), the Gansu Province Science and Technology Plan Project (Key R&D Program) (23YFFA0027), and the Gansu Province Higher Education Innovation Fund Project (2021A-441).

Acknowledgments

The author thanked the Gansu Provincial Key Laboratory of Herbivore Biotechnology for providing a research platform.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1471732/full#supplementary-material

Footnotes

1. ^www.biocloud.net

References

Abbas, W., Howard, J. T., Paz, H. A., Hales, K. E., Wells, J. E., Kuehn, L. A., et al. (2020). Influence of host genetics in shaping the rumen bacterial community in beef cattle. Sci. Rep. 10:15101. doi: 10.1038/s41598-020-72011-9

PubMed Abstract | Crossref Full Text | Google Scholar

Baldwin, R. T., and Jesse, B. W. (1992). Developmental changes in glucose and butyrate metabolism by isolated sheep ruminal cells. J. Nutr. 122, 1149–1153. doi: 10.1093/jn/122.5.1149

PubMed Abstract | Crossref Full Text | Google Scholar

Bergman, E. N. (1990). Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol. Rev. 70, 567–590. doi: 10.1152/physrev.1990.70.2.567

PubMed Abstract | Crossref Full Text | Google Scholar

Berry, D., Mader, E., Lee, T. K., Woebken, D., Wang, Y., Zhu, D., et al. (2015). Tracking heavy water (D2O) incorporation for identifying and sorting active microbial cells. Proc. Natl. Acad. Sci. USA 112, E194–E203. doi: 10.1073/pnas.1420406112

PubMed Abstract | Crossref Full Text | Google Scholar

Bian, X., Wu, W., Yang, L., Lv, L., Wang, Q., Li, Y., et al. (2019). Administration of Akkermansia muciniphila ameliorates dextran sulfate sodium-induced ulcerative colitis in mice. Front. Microbiol. 10:2259. doi: 10.3389/fmicb.2019.02259

PubMed Abstract | Crossref Full Text | Google Scholar

Bilk, S., Huhn, K., Honscha, K. U., Pfannkuche, H., and Gäbel, G. (2005). Bicarbonate exporting transporters in the ovine ruminal epithelium. J. Comp. Physiol. B 175, 365–374. doi: 10.1007/s00360-005-0493-1

PubMed Abstract | Crossref Full Text | Google Scholar

Cortes, V. A., Busso, D., Maiz, A., Arteaga, A., Nervi, F., and Rigotti, A. (2014). Physiological and pathological implications of cholesterol. Front. Biosci. 19, 416–428. doi: 10.2741/4216

Crossref Full Text | Google Scholar

Cunha, I. S., Barreto, C. C., Costa, O. Y. A., Bomfim, M. A., Castro, A. P., Kruger, R. H., et al. (2011). Bacteria and Archaea community structure in the rumen microbiome of goats (Capra hircus) from the semiarid region of Brazil. Anaerobe 17, 118–124. doi: 10.1016/j.anaerobe.2011.04.018

Crossref Full Text | Google Scholar

de Aguiar Vallim, T. Q., Tarling, E. J., and Edwards, P. A.. (2013). Pleiotropic roles of bile acids in metabolism. Cell Metab. 17, 657–669. doi: 10.1016/j.cmet.2013.03.013

PubMed Abstract | Crossref Full Text | Google Scholar

del Bianco Benedeti, P., Silva, B. C., Pacheco, M. V. C., Serão, N. V. L., Carvalho Filho, I., Lopes, M. M., et al. (2018). Effects of grain processing methods on the expression of genes involved in volatile fatty acid transport and pH regulation, and keratinization in rumen epithelium of beef cattle. PLoS One 13:e0198963. doi: 10.1371/journal.pone.0198963

PubMed Abstract | Crossref Full Text | Google Scholar

Difford, G. F., Plichta, D. R., Løvendahl, P., Lassen, J., Noel, S. J., Højberg, O., et al. (2018). Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cows. PLoS Genet. 14:e1007580. doi: 10.1371/journal.pgen.1007580

PubMed Abstract | Crossref Full Text | Google Scholar

EEr, H., Ma, L., Xie, X., Ma, J., Ma, X., Yue, C., et al. (2020). Genetic polymorphism association analysis of SNPs on the species conservation genes of Tan sheep and Hu sheep. Trop. Anim. Health Prod. 52, 915–926. doi: 10.1007/s11250-019-02063-1

PubMed Abstract | Crossref Full Text | Google Scholar

Fan, Q., Cui, X., Wang, Z., Chang, S., Wanapat, M., Yan, T., et al. (2021). Rumen microbiota of Tibetan sheep (Ovis aries) adaptation to extremely cold season on the Qinghai-Tibetan plateau. Front Vet Sci 8:673822. doi: 10.3389/fvets.2021.673822

PubMed Abstract | Crossref Full Text | Google Scholar

Ferah Okkay, I., Okkay, U., Gundogdu, O. L., Bayram, C., Mendil, A. S., Ertugrul, M. S., et al. (2022). Syringic acid protects against thioacetamide-induced hepatic encephalopathy: behavioral, biochemical, and molecular evidence. Neurosci. Lett. 769:136385. doi: 10.1016/j.neulet.2021.136385

PubMed Abstract | Crossref Full Text | Google Scholar

Fraune, S., and Bosch, T. C. (2007). Long-term maintenance of species-specific bacterial microbiota in the basal metazoan Hydra. Proc. Natl. Acad. Sci. USA 104, 13146–13151. doi: 10.1073/pnas.0703375104

PubMed Abstract | Crossref Full Text | Google Scholar

Fugaban, J., Vazquez Bucheli, J. E., Park, Y. J., Suh, D. H., Jung, E. S., Franco, B. D. G. M., et al. (2022). Antimicrobial properties of Pediococcus acidilactici and Pediococcus pentosaceus isolated from silage. J. Appl. Microbiol. 132, 311–330. doi: 10.1111/jam.15205

PubMed Abstract | Crossref Full Text | Google Scholar

Gonzalez, B., and Cullen, P. J. (2022). Regulation of Cdc42 protein turnover modulates the filamentous growth MAPK pathway. J. Cell Biol. 221, 1307–1316. doi: 10.1083/jcb.202112100

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, R., Wang, H., Meng, C., Gui, H., Li, Y., Chen, F., et al. (2023). Efficient and specific generation of MSTN-edited Hu sheep using C-CRISPR. Genes (Basel) 14, 1216–1231. doi: 10.3390/genes14061216

PubMed Abstract | Crossref Full Text | Google Scholar

Hess, M., Sczyrba, A., Egan, R., Kim, T. W., Chokhawala, H., Schroth, G., et al. (2011). Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science 331, 463–467. doi: 10.1126/science.1200387

PubMed Abstract | Crossref Full Text | Google Scholar

Hosseinkhani, F., Heinken, A., Thiele, I., Lindenburg, P. W., Harms, A. C., and Hankemeier, T. (2021). The contribution of gut bacterial metabolites in the human immune signaling pathway of non-communicable diseases. Gut Microbes 13, 1–22. doi: 10.1080/19490976.2021.1882927

PubMed Abstract | Crossref Full Text | Google Scholar

Huo, W., Zhu, W., and Mao, S. (2014). Impact of subacute ruminal acidosis on the diversity of liquid and solid-associated bacteria in the rumen of goats. World J. Microbiol. Biotechnol. 30, 669–680. doi: 10.1007/s11274-013-1489-8

PubMed Abstract | Crossref Full Text | Google Scholar

Imenshahidi, M., and Hossenzadeh, H. (2022). Effects of glycine on metabolic syndrome components: a review. J. Endocrinol. Investig. 45, 927–939. doi: 10.1007/s40618-021-01720-3

PubMed Abstract | Crossref Full Text | Google Scholar

Jami, E., White, B. A., and Mizrahi, I. (2014). Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS One 9:e85423. doi: 10.1371/journal.pone.0085423

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, M., Morrison, M., and Yu, Z. (2011). Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol. Ecol. 76, 49–63. doi: 10.1111/j.1574-6941.2010.01029.x

PubMed Abstract | Crossref Full Text | Google Scholar

Kirat, D., Masuoka, J., Hayashi, H., Iwano, H., Yokota, H., Taniyama, H., et al. (2006). Monocarboxylate transporter 1 (MCT1) plays a direct role in short-chain fatty acids absorption in caprine rumen. J. Physiol. 576, 635–647. doi: 10.1113/jphysiol.2006.115931

PubMed Abstract | Crossref Full Text | Google Scholar

Konishi, H., Fujiya, M., and Kohgo, Y. (2015). Host-microbe interactions via membrane transport systems. Environ. Microbiol. 17, 931–937. doi: 10.1111/1462-2920.12632

Crossref Full Text | Google Scholar

Li, Y., Chen, Z., Fang, Y., Cao, C., Zhang, Z., Pan, Y., et al. (2022). Runs of homozygosity revealed reproductive traits of Hu sheep. Genes (Basel) 13, 1848–1864. doi: 10.3390/genes13101848

PubMed Abstract | Crossref Full Text | Google Scholar

Li, G., Lin, J., Zhang, C., Gao, H., Lu, H., Gao, X., et al. (2021). Microbiota metabolite butyrate constrains neutrophil functions and ameliorates mucosal inflammation in inflammatory bowel disease. Gut Microbes 13:1968257. doi: 10.1080/19490976.2021.1968257

Crossref Full Text | Google Scholar

Lin, X., Wang, J., Hou, Q., Wang, Y., Hu, Z., Shi, K., et al. (2018). Effect of hay supplementation timing on rumen microbiota in suckling calves. Microbiology 7:e00430. doi: 10.1002/mbo3.430

Crossref Full Text | Google Scholar

Lin, L., Xie, F., Sun, D., Liu, J., Zhu, W., and Mao, S. (2019). Ruminal microbiome-host crosstalk stimulates the development of the ruminal epithelium in a lamb model. Microbiome 7:83. doi: 10.1186/s40168-019-0701-y

Crossref Full Text | Google Scholar

Liu, H., Hu, L., Han, X., Zhao, N., Xu, T., Ma, L., et al. (2020). Tibetan sheep adapt to plant phenology in alpine meadows by changing rumen microbial community structure and function. Front. Microbiol. 11:587558. doi: 10.3389/fmicb.2020.587558

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Sha, Y., Lv, W., Cao, G., Guo, X., Pu, X., et al. (2022). Multi-omics reveals that the rumen transcriptome, microbiome, and its metabolome co-regulate cold season adaptability of Tibetan sheep. Front. Microbiol. 13:859601. doi: 10.3389/fmicb.2022.859601

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Xu, T., Xu, S., Ma, L., Han, X., Wang, X., et al. (2019). Effect of dietary concentrate to forage ratio on growth performance, rumen fermentation and bacterial diversity of Tibetan sheep under barn feeding on the Qinghai-Tibetan plateau. PeerJ 7:e7462. doi: 10.7717/peerj.7462

PubMed Abstract | Crossref Full Text | Google Scholar

McGarry, J. D., and Foster, D. W. (1980). Regulation of hepatic fatty acid oxidation and ketone body production. Annu. Rev. Biochem. 49, 395–420. doi: 10.1146/annurev.bi.49.070180.002143

Crossref Full Text | Google Scholar

Meng, H., Zhang, Y., Zhao, L., Zhao, W., He, C., Honaker, C. F., et al. (2014). Body weight selection affects quantitative genetic correlated responses in gut microbiota. PLoS One 9:e89862. doi: 10.1371/journal.pone.0089862

PubMed Abstract | Crossref Full Text | Google Scholar

Mihalik, S. J., Steinberg, S. J., Pei, Z., Park, J., Kim, D. G., Heinzer, A. K., et al. (2002). Participation of two members of the very long-chain acyl-CoA synthetase family in bile acid synthesis and recycling. J. Biol. Chem. 277, 24771–24779. doi: 10.1074/jbc.M203295200

Crossref Full Text | Google Scholar

Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., González, A., Fontana, L., et al. (2011). Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science 332, 970–974. doi: 10.1126/science.1198719

PubMed Abstract | Crossref Full Text | Google Scholar

Müller, F., Huber, K., Pfannkuche, H., Aschenbach, J. R., Breves, G., and Gäbel, G. (2002). Transport of ketone bodies and lactate in the sheep ruminal epithelium by monocarboxylate transporter 1. Am. J. Physiol. Gastrointest. Liver Physiol. 283, G1139–G1146. doi: 10.1152/ajpgi.00268.2001

PubMed Abstract | Crossref Full Text | Google Scholar

Parker, B. J., Wearsch, P. A., Veloo, A. C. M., and Rodriguez-Palacios, A. (2020). The genus Alistipes: gut Bacteria with emerging implications to inflammation, Cancer, and mental health. Front. Immunol. 11:906. doi: 10.3389/fimmu.2020.00906

PubMed Abstract | Crossref Full Text | Google Scholar

Peng, B., Huang, S., Liu, T., and Geng, A. (2015). Bacterial xylose isomerases from the mammal gut Bacteroidetes cluster function in Saccharomyces cerevisiae for effective xylose fermentation. Microb. Cell Factories 14:70. doi: 10.1186/s12934-015-0253-1

PubMed Abstract | Crossref Full Text | Google Scholar

Pérez de Heredia, F., Wood, I. S., and Trayhurn, P. (2010). Hypoxia stimulates lactate release and modulates monocarboxylate transporter (MCT1, MCT2, and MCT4) expression in human adipocytes. Pflugers Arch. 459, 509–518. doi: 10.1007/s00424-009-0750-3

PubMed Abstract | Crossref Full Text | Google Scholar

Petrenko, V., Sinturel, F., Riezman, H., and Dibner, C. (2023). Lipid metabolism around the body clocks. Prog. Lipid Res. 91:101235. doi: 10.1016/j.plipres.2023.101235

Crossref Full Text | Google Scholar

Qi, M., Wang, J., Tan, B.'., Li, J., Liao, S., Liu, Y., et al. (2020). Dietary glutamine, glutamate, and aspartate supplementation improves hepatic lipid metabolism in post-weaning piglets. Anim. Nutr. 6, 124–129. doi: 10.1016/j.aninu.2019.12.002

PubMed Abstract | Crossref Full Text | Google Scholar

Rabe, K. F., and Dent, G. (1998). Theophylline and airway inflammation. Clin. Exp. Allergy 28 Suppl 3, 35–41

PubMed Abstract | Google Scholar

Rahman, K., and Coleman, R. (1986). Selective biliary lipid secretion at low bile-salt-output rates in the isolated perfused rat liver. Effects of phalloidin. Biochem J 237, 301–304. doi: 10.1042/bj2370301

PubMed Abstract | Crossref Full Text | Google Scholar

Sha, Y., Ren, Y., Zhao, S., He, Y., Guo, X., Pu, X., et al. (2022). Response of ruminal microbiota-host gene interaction to high-altitude environments in Tibetan sheep. Int. J. Mol. Sci. 23:12430. doi: 10.3390/ijms232012430

PubMed Abstract | Crossref Full Text | Google Scholar

Shi, W., Moon, C. D., Leahy, S. C., Kang, D., Froula, J., Kittelmann, S., et al. (2014). Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome. Genome Res. 24, 1517–1525. doi: 10.1101/gr.168245.113

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q., Wang, K., Wu, W., Giannoulatou, E., Ho, J. W. K., and Li, L. (2019). Host and microbiome multi-omics integration: applications and methodologies. Biophys. Rev. 11, 55–65. doi: 10.1007/s12551-018-0491-7

PubMed Abstract | Crossref Full Text | Google Scholar

Xue, M. Y., Sun, H. Z., Wu, X. H., Liu, J. X., and Guan, L. L. (2020). Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance. Microbiome 8:64. doi: 10.1186/s40168-020-00819-8

PubMed Abstract | Crossref Full Text | Google Scholar

Yagi, C., Kusunoki, Y., Tsunoda, T., Murase, T., Nakamura, T., Osugi, K., et al. (2022). Xanthine oxidoreductase activity is correlated with hepatic steatosis. Sci. Rep. 12:12282. doi: 10.1038/s41598-022-16688-0

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, W., Shen, Z., and Martens, H. (2012). An energy-rich diet enhances expression of Na(+)/H(+) exchanger isoform 1 and 3 messenger RNA in rumen epithelium of goat. J. Anim. Sci. 90, 307–317. doi: 10.2527/jas.2011-3854

PubMed Abstract | Crossref Full Text | Google Scholar

Zang, X. W., Sun, H. Z., Xue, M. Y., Zhang, Z., Plastow, G., Yang, T., et al. (2022). Heritable and nonheritable rumen Bacteria are associated with different characters of lactation performance of dairy cows. Microbial Systems 7:e0042222. doi: 10.1128/msystems.00422-22

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Q., Difford, G., Sahana, G., Løvendahl, P., Lassen, J., Lund, M. S., et al. (2020). Bayesian modeling reveals host genetics associated with rumen microbiota jointly influence methane emission in dairy cows. ISME J. 14, 2019–2033. doi: 10.1038/s41396-020-0663-x

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Q., Jian, L., Yao, D., Rao, B., Xia, Y., Hu, K., et al. (2023). The structural basis of the pH-homeostasis mediated by the cl(−)/HCO(3)(−) exchanger, AE2. Nat. Commun. 14:1812. doi: 10.1038/s41467-023-37557-y

Crossref Full Text | Google Scholar

Zhang, H., Shao, M., Huang, H., Wang, S., Ma, L., Wang, H., et al. (2018). The dynamic distribution of small-tail Han sheep microbiota across different intestinal segments. Front. Microbiol. 9:32. doi: 10.3389/fmicb.2018.00032

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Z., Xu, D., Wang, L., Hao, J., Wang, J., Zhou, X., et al. (2016). Convergent evolution of rumen microbiomes in high-altitude mammals. Curr. Biol. 26, 1873–1879. doi: 10.1016/j.cub.2016.05.012

PubMed Abstract | Crossref Full Text | Google Scholar