Study on the Plateau Adaptive Synergistic Mechanism of Rumen Microbiome-Metabolome-Resistome in Tibetan Sheep

Tibetan sheep are an important livestock breed adapted to the extreme environment of the Qinghai-Tibet Plateau (QTP). Their energy metabolism and environmental adaptability are highly dependent on the rumen microbiome. However, systematic comparisons of the rumen microbiome, its functions, and the resistome between plateau-adapted breeds and lowland breeds remain insufficient. In this study, 6 Tibetan sheep (TS) and 6 Hu sheep (HS) were selected. All the selected sheep had a body weight of 34 kg (±0.5 kg) and an age of 1 year (±1 month) and were all managed under local traditional natural grazing (without supplementary feeding). Using metagenomics and metabolomics techniques, systematic comparative analysis was conducted on the differences in rumen microbial community structure, functions, resistome, and metabolites between the two breeds. Metagenomic analysis showed that at the phylum level, the abundance of Bacteroidetes in the rumen of TS was significantly higher than that in HS (p < 0.05); at the genus level, the abundance of Bacteroides in TS was also significantly higher (p < 0.05). Carbohydrate-active enzymes (CAZy) analysis indicated that the abundance of Glycosyltransferases (GTs) and Carbohydrate-Binding Modules (CBMs) in the rumen of TS were significantly upregulated (p < 0.05), while HS was rich in various Glycoside Hydrolases (GHs). Comprehensive Antibiotic Resistance Database (CARD) analysis revealed that more than 60% of the Antibiotic Resistance Genes (ARGs) in the rumen of HS were present at higher levels than those in TS. Metabolomics identified a large number of differential metabolites, among which metabolites such as 2E,6Z,8Z,12E-hexadecatetraenoic acid, Leukotriene F4, and Tenurin were significantly upregulated in the rumen of TS. Correlation analysis showed that rumen microbial flora and their metabolites jointly act to regulate rumen ARGs. Specifically, microorganisms including Firmicutes and Succiniclasticum had a significantly positive correlation with ARGs such as rpoB2 (p < 0.05), while differential metabolites like endomorphin-1 and Purothionin AII exhibited a significantly negative correlation with ARGs such as rpoB2 (p < 0.05). Therefore, compared with HS, the synergistic effect of the rumen microbial flora, their metabolites, and the resistome in TS is an important characteristic and strategy for their adaptation to the hypoxic environment of the QTP, and also contributes to the formation of their unique rumen resistome. Despite being reared in the same plateau environment, the rumen microbiome of HS still retains low-altitude characteristics, which are manifested as high GHs activity and high ARGs abundance.