Altitude is the main external environmental pressure affecting the production performance of Tibetan sheep, and the adaptive evolution of many years has formed a certain response mechanism. However, there are few reports on the response of ruminal microbiota and host genomes of Tibetan sheep to high-altitude environments. Here, we conducted an integrated analysis of volatile fatty acids (VFAs), microbial diversity (16S rRNA), epithelial morphology, and epithelial transcriptome in the rumen of Tibetan sheep at different altitudes to understand the changes in ruminal microbiota–host interaction in response to high altitude. The differences in the nutritional quality of forage at different altitudes, especially the differences in fiber content (ADF/NDF), led to changes in rumen VFAs of Tibetan sheep, in which the A/P value (acetic acid/propionic acid) was significantly decreased (p < 0.05). In addition, the concentrations of IgA and IgG in Middle-altitude (MA) and High-altitude Tibetan sheep (HA) were significantly increased (p < 0.05), while the concentrations of IgM were significantly increased in MA (p < 0.05). Morphological results showed that the width of the rumen papilla and the thickness of the basal layer increased significantly in HA Tibetan sheep (p < 0.05). The 16S rRNA analysis found that the rumen microbial diversity of Tibetan sheep gradually decreased with increasing altitude, and there were some differences in phylum- and genus-level microbes at the three altitudes. RDA analysis found that the abundance of the Rikenellaceae RC9 gut group and the Ruminococcaceae NK4A214 group increased with altitudes. Furthermore, a functional analysis of the KEGG microbial database found the “lipid metabolism” function of HA Tibetan sheep to be significantly enriched. WGCNA revealed that five gene modules were enriched in “energy production and conversion”, “lipid transport and metabolism”, and “defense mechanisms”, and cooperated with microbiota to regulate rumen fermentation and epithelial immune barrier function, so as to improve the metabolism and immune level of Tibetan sheep at high altitude.

As an important ruminant on the Qinghai-Tibet Plateau, Tibetan sheep can maintain their population reproduction rate in the harsh high-altitude environment of low temperature and low oxygen, which relies on their special plateau adaptations mechanism that they have formed for a long time. Microbiomes (known as “second genomes”) are closely related to the nutrient absorption, adaptability, and health of the host. In this study, rumen fermentation characteristics, the microbiota, and rumen epithelial gene expression of Tibetan sheep in various months were analyzed. The results show that the rumen fermentation characteristics of Tibetan sheep differed in different months. The total SCFAs (short-chain fatty acids), acetate, propionate, and butyrate concentrations were highest in October and lowest in June. The CL (cellulase) activity was highest in February, while the ACX (acid xylanase) activity was highest in April. In addition, the diversity and abundance of rumen microbes differed in different months. Bacteroidetes (53.4%) and Firmicutes (27.4%) were the dominant phyla. Prevotella_1 and Rikenellaceae_RC9_gut_group were the dominant genera. The abundance of Prevotella_1 was highest in June (27.8%) and lowest in December (17.8%). In addition, the expression of CLAUDIN4 (Claudin-4) and ZO1 (Zonula occludens 1) was significantly higher in April than in August and December, while the expression of SGLT1 (Sodium glucose linked transporter 1) was highest in August. Correlation analysis showed that there were interactions among rumen fermentation characteristics, the microbiota, and host gene expression, mainly by adjusting the amino acid metabolism pathway and energy metabolism pathway to improve energy utilization. At the same time, we adjusted the balance of the rumen “core microbiota” to promote the development of rumen and maintain the homeostasis of rumen environment, which makes Tibetan sheep better able to adapt to the harsh environment in different periods of the Qinghai-Tibet Plateau.

Plateau adaptation in animals involves genetic mechanisms as well as coevolutionary mechanisms of the microbiota and metabolome of the animal. Therefore, the characteristics of the rumen microbiome and metabolome, transcriptome, and serum metabolome of Tibetan sheep at different altitudes (4500 m, 3500 m, and 2500 m) were analyzed. The results showed that the rumen differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and carbohydrate metabolism, and there was a significant correlation with microbiota. The differentially expressed genes and metabolites at middle and high altitudes were coenriched in asthma, arachidonic acid metabolism, and butanoate and propanoate metabolism. In addition, the serum differential metabolites at 3500 m and 4500 m were mainly enriched in amino acid metabolism, lipid metabolism, and metabolism of xenobiotics by cytochrome P450, and they were also related to microbiota. Further analysis revealed that rumen metabolites accounted for 7.65% of serum metabolites. These common metabolites were mainly enriched in metabolic pathways and were significantly correlated with host genes (p < 0.05). This study found that microbiota, metabolites, and epithelial genes were coenriched in pathways related to lipid metabolism, energy metabolism, and immune metabolism, which may be involved in the regulation of Tibetan sheep adaptation to plateau environmental changes.

This study explores the effects of the coevolution of the host genome (the first genome) and gut microbiome (the second genome) on nutrition stress in Tibetan sheep during the cold season. The rumen epithelial tissue of six Tibetan sheep (Oula-type) was collected as experimental samples during the cold and warm seasons and the study lasted for half a year. The cDNA library was constructed and subjected to high-throughput sequencing. The circRNAs with significant differential expression were identified through bioinformatics analysis and functional prediction, and verified by real-time quantitative PCR (qRT-PCR). The results showed that a total of 56 differentially expressed (DE) circRNAs of rumen epithelial tissue were identified using RNA-seq technology, among which 29 were significantly upregulated in the cold season. The circRNA-miRNA regulatory network showed that DE circRNAs promoted the adaptation of Tibetan sheep in the cold season by targeting miR-150 and oar-miR-370-3p. The results of correlation analysis among circRNAs, microbiota, and metabolites showed that the circRNA NC_040275.1:28680890|28683112 had a very significant positive correlation with acetate, propionate, butyrate, and total volatile fatty acid (VFA) (p < 0.01), and had a significant positive correlation with Ruminococcus-1 (p < 0.05). In addition, circRNA NC_040256.1:78451819|78454934 and metabolites were enriched in the same KEGG pathway biosynthesis of amino acids (ko01230). In conclusion, the host genome and rumen microbiome of Tibetan sheep co-encoded a certain glycoside hydrolase (β-glucosidase) and coevolved efficient VFA transport functions and amino acid anabolic processes; thus, helping Tibetan sheep adapt to nutrient stress in the cold season in high-altitude areas.

Based on the serious phenomenon of Tibetan sheep “growing strong in warm seasons and losing weight in cold seasons”, this study explores the regulation of lncRNAs, microbiota, and metabolites in the cold season adaptation of Tibetan sheep from the perspective of the co-evolution of the host genome (first genome) and microbiome (second genome). RNA-seq results showed that 172 DE lncRNAs were identified in the rumen epithelium of Tibetan sheep in warm and cold seasons, of which 87 DE lncRNAs were significantly up-regulated in cold seasons. KEGG enrichment showed that target genes of up-regulated lncRNAs were significantly enriched in TNF signaling and oxidative phosphorylation pathways. LncRNA-mRNA regulatory network indicated that DE lncRNAs were involved in nutrient stress in the cold season by targeting ATP1B2, CADPS, TLR5, and UGT1A6. Correlation analysis showed some lncRNAs were significantly correlated with acetic acid, propionic acid, butyric acid, and rumen epithelial histomorphology and had a negative correlation with Butyrivibrio-2 and Succiniclasticum (p < 0.05). In addition, differential metabolites bilirubin and lncRNAs were co-enriched in the bile secretion pathway. lncRNAs played an important role in the adaptation process of Tibetan sheep in the cold season, and mediate the host to participate in nutrient absorption, energy utilization, and immune response, indicating that the host genome and microbial genome promote Tibetan sheep to adapt to nutrient stress in the cold season through co-evolution.